This paper presents the calculation of elastic properties of three-dimensional Yukawa or dust crystals using molecular dynamics simulations. The elastic properties are computed by deforming (compressing/expanding) the dust crystals along different directions. The stress and strain of the deformed crystal are used to calculate elastic properties. The bulk modulus, shear modulus, and Poisson's ratio are determined as a function of shielding parameter κ and strong coupling parameter Γ. The bulk and shear modulus values at 0 K temperature are consistent with the previous literature results, while the finite-temperature results are new. The finite-temperature bulk modulus of Yukawa crystals is found to be higher than that of 0 K crystals. The shear modulus of the Yukawa solids decreases nonlinearly near the solid–liquid boundary in the premelting region. The Poisson's ratio of Yukawa crystals changes sharply at the solid–liquid boundary, emphasizing its potential for identifying phase transitions and assessing incompressibility in Yukawa systems. The bulk and shear moduli calculated in this paper are useful for determining accurate values of sound and shear velocity in Yukawa systems across a wide range of the (κ, Γ) parameter space.

We present a study on the transport and materials properties of aluminum spanning from ambient to warm dense matter conditions using a machine-learned interatomic potential (ML-IAP). Prior research has utilized ML-IAPs to simulate phenomena in warm dense matter, but these potentials have often been calibrated for a narrow range of temperature and pressures. In contrast, we train a single ML-IAP over a wide range of temperatures, using density functional theory molecular dynamics (DFT-MD) data. Our approach overcomes computational limitations of DFT-MD simulations, enabling us to study transport and materials properties of matter at higher temperatures and longer time scales. We demonstrate the ML-IAP transferability across a wide range of temperatures using molecular-dynamics (MD) by examining the thermal conductivity, diffusion coefficient, viscosity, sound velocity, and ion-ion structure factor of aluminum up to about 60,000 K, where we find good agreement with previous theoretical data.

The small semi-major axes of Hot Jupiters lead to high atmospheric temperatures of up to several thousand Kelvin. Under these conditions, thermally ionised metals provide a rich source of charged particles and thus build up a sizeable electrical conductivity. Subsequent electromagnetic effects, such as the induction of electric currents, Ohmic heating, magnetic drag, or the weakening of zonal winds have thus far been considered mainly in the framework of a linear, steady-state model of induction. For Hot Jupiters with an equilibrium temperature T_eq > 1500 K, the induction of atmospheric magnetic fields is a runaway process that can only be stopped by non-linear feedback. For example, the back-reaction of the magnetic field onto the flow via the Lorentz force or the occurrence of magnetic instabilities. Moreover, we discuss the possibility of self-excited atmospheric dynamos. Our results suggest that the induced atmospheric magnetic fields and electric currents become independent of the electrical conductivity and the internal field, but instead are limited by the planetary rotation rate and wind speed. As an explicit example, we characterise the induction process for the hottest exoplanet, KELT-9b by calculating the electrical conductivity along atmospheric P-T-profiles for the day- and nightside. Despite the temperature varying between 3000 K and 4500 K, the resulting electrical conductivity attains an elevated value of roughly 1 S/m throughout the atmosphere. The induced magnetic fields are predominately horizontal and might reach up to a saturation field strength of 400 mT, exceeding the internal field by two orders of magnitude.

We study ionization and transport processes in partially ionized multicomponent plasmas. The plasma composition is calculated via a system of coupled mass-action laws. The electronic transport properties are determined by the electron-ion and electron-neutral transport cross sections. The influence of electron-electron scattering is considered via a correction factor to the electron-ion contribution. Based on these data, the electrical and thermal conductivities as well as the Lorenz number are calculated. For the thermal conductivity, we consider also the contributions of the translational motion of neutral particles and of the dissociation, ionization, and recombination reactions. We apply our approach to a partially ionized plasma composed of hydrogen, helium, and a small fraction of metals (Li, Na, Ca, Fe, K, Rb, and Cs) as typical for atmospheres of hot Jupiters. We present results for the plasma composition and the transport properties as a function of density and temperature and then along typical P-T profiles for the outer part of the hot Jupiter HD 209458b. The electrical conductivity profile allows revising the Ohmic heating power related to the fierce winds in the planet's atmosphere. We show that the higher temperatures suggested by recent interior models could boost the conductivity and thus the Ohmic heating power to values large enough to explain the observed inflation of HD 209458b.

A pair of counter-rotating symmetric vortices has been observed in the wake behind a stationary obstacle (dust void) in a flowing dusty plasma. A strongly coupled dusty plasma flow with controllable velocity is generated and directed toward the void in a novel experiment. In the unsteady laminar flow regime, the curl of the fluid flow velocity along the boundary layer of the void generates the vortex pair behind the void. Particle image velocimetry analysis of high speed image data clearly depicts the flow pattern and the vorticities. The shear viscosity of the dusty plasma fluid along with the experimental parameters is considered to obtain the Reynolds number range for the evolution of the vortices. Molecular dynamics simulation is also performed to support the experimental observation.

Spiral wave formations are ubiquitous in nature. In the present paper, the excitation of spiral waves in the context of driven two-dimensional dusty plasma (Yukawa system) has been demonstrated at particle level using molecular-dynamics simulations. The interaction amidst dust particles is modeled by the Yukawa potential to take account of the shielding of dust charges by the lighter electron and ion species. The spatiotemporal evolution of these spiral waves has been characterized as a function of the frequency and amplitude of the driving force and dust neutral collisions. The effect of strong coupling has been studied, which shows that the excited spiral wave structures get clearer as the medium gets more strongly coupled. The radial propagation speed of the spiral wave is observed to remain unaltered with the coupling parameter. However, it is found to depend on the screening parameter of the dust medium and decreases when it is increased. In the crystalline phase (with screening parameter κ>0.58), the spiral wavefronts are shown to be hexagonal in shape. This shows that the radial propagation speed depends on the interparticle spacing.

The collective response of the plasma medium is well known and has been explored extensively in the context of dusty plasma medium. On the other hand, the individual particle response associated with the collisional character giving rise to the dissipative phenomena has not been explored adequately. In this paper, two-dimensional molecular dynamics simulation of dust particles interacting via Yukawa potential has been considered. It has been shown that disturbances induced in a dust crystal elicit both collective and single particle responses. Generation of a few particles moving at speeds considerably higher than acoustic and/or shock speed (excited by the external disturbance) is observed. This is an indication of a single particle response. Furthermore, as these individual energetic particles propagate, the dust crystal is observed to crack along their path. Initially when the energy is high, these particles generate secondary energetic particles by the collisional scattering process. However, ultimately as these particles slow down they excite a collective response in the dust medium at secondary locations in a region which is undisturbed by the primary external disturbance. The condition when the cracking of the crystal stops and collective excitations get initiated has been identified quantitatively. The trailing collective primary disturbances would thus often encounter a disturbed medium with secondary and tertiary collective perturbations, thereby suffering significant modification in its propagation. It is thus clear that there is an interesting interplay (other than mere dissipation) between the single particle and collective response which governs the dynamics of any disturbance introduced in the medium.

Spiral waves are observed in many natural phenomena. They have been extensively represented by the mathematical FitzHugh-Nagumo (FHN) model [Barkley et al., Phys. Rev. A 42, 2489 (1990)] of excitable media. In incompressible fluid simulations also excitation of thermal spiral waves have been reported by Li et al. [Phys. of Fluids 22, 011701 (2010)]. In the present manuscript the spatiotemporal development of spiral waves in the context of weak and strong coupling limits have been shown. While the weakly coupled medium has been represented by a simple fluid description, for the strong coupling a generalized visco - elastic fluid description has been employed. The medium has been driven by an external force in the form of a rotating electric field. It is shown that when the amplitude of force is small the density perturbations in the medium are also small. In this case the excitations do not develop as a spiral wave. Only when the amplitude of force is high so as to drive the density perturbations to nonlinear amplitudes does the spiral density wave formation occur. The role of forcing frequency, the effect of strong coupling and sound velocity of medium on the formation and evolution of spiral waves have been investigated in detail.

The excitation and evolution of Korteweg-de Vries (KdV) solitons in a dusty plasma medium are studied using Molecular Dynamics (MD) simulations. The dusty plasma medium is modelled as a collection of dust particles interacting through Yukawa potential, which takes into account dust charge screening due to the lighter electron and ion species. The collective response of such screened dust particles to an applied electric field impulse is studied here. An excitation of a perturbed positive density pulse propagating in one direction along with a train of negative perturbed rarefactive density oscillations in the opposite direction is observed. These observations are in accordance with evolution governed by the KdV equation. Detailed studies of (a) amplitude vs. width variation of the observed pulse, (b) the emergence of intact separate pulses with an associated phase shift after collisional interaction amidst them, etc., conclusively qualify the positive pulses observed in the simulations as KdV solitons. It is also observed that by increasing the strength of the electric field impulse, multiple solitonic structures get excited. The excitations of the multiple solitons are similar to the experimental observations reported recently by Boruah et al. [Phys. Plasmas 23, 093704 (2016)] for dusty plasmas. The role of coupling parameter has also been investigated here, which shows that with increasing coupling parameter, the amplitude of the solitonic pulse increases whereas its width decreases.