1
|
Djienbekov NE, Bastykova NK, Ramazanov TS, Kodanova SK. Non-equilibrium molecular dynamics study of heat transfer parameters in two-dimensional Yukawa systems under uniform magnetic field. Sci Rep 2024; 14:15042. [PMID: 38951572 PMCID: PMC11217476 DOI: 10.1038/s41598-024-64866-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2024] [Accepted: 06/13/2024] [Indexed: 07/03/2024] Open
Abstract
The present study explores the effect of a magnetic field on the thermal conductivity of two-dimensional (2D) Yukawa systems in a wide range of system parameters using the non-equilibrium molecular dynamic method (NEMD). We consider an external magnetic field with Ω = ω c / ω p ≤ 1 (with Ω being the ratio of the cyclotron frequency to plasma frequency) and the coupling parameter values in the range 1 ≤ Γ ≤ 100 (with Γ being the ratio of the Coulomb interaction energy at mean inter-particle distance to the thermal energy of particles). The results show that an external uniform magnetic field results in the reduction of the thermal conductivity at the considered values of the coupling parameter Γ . Additionally, we found that the effect of the magnetic field on thermal conduction is weaker at larger values of the system coupling parameter. To ensure that calculated results for the thermal conductivity are accurate and reliable, we performed a detailed investigation of the convergence of the results with respect to simulation parameters in NEMD with a strong external magnetic field. We believe that the presented results will serve as useful benchmark data for the theoretical models of (2D) Yukawa systems.
Collapse
Affiliation(s)
- N E Djienbekov
- Institute for Experimental and Theoretical Physics, Al-Farabi Kazakh National University, 71 Al-Farabi Ave., 050040, Almaty, Kazakhstan
- Institute of Applied Sciences and IT, 40-48 Shashkin Str., 050038, Almaty, Kazakhstan
| | - N Kh Bastykova
- Institute for Experimental and Theoretical Physics, Al-Farabi Kazakh National University, 71 Al-Farabi Ave., 050040, Almaty, Kazakhstan.
- Institute of Applied Sciences and IT, 40-48 Shashkin Str., 050038, Almaty, Kazakhstan.
| | - T S Ramazanov
- Institute for Experimental and Theoretical Physics, Al-Farabi Kazakh National University, 71 Al-Farabi Ave., 050040, Almaty, Kazakhstan
- Institute of Applied Sciences and IT, 40-48 Shashkin Str., 050038, Almaty, Kazakhstan
| | - S K Kodanova
- Institute for Experimental and Theoretical Physics, Al-Farabi Kazakh National University, 71 Al-Farabi Ave., 050040, Almaty, Kazakhstan
- Institute of Applied Sciences and IT, 40-48 Shashkin Str., 050038, Almaty, Kazakhstan
| |
Collapse
|
2
|
Pak C, Billings V, Schlitters M, Bergeson SD, Murillo MS. Preliminary study of plasma modes and electron-ion collisions in partially magnetized strongly coupled plasmas. Phys Rev E 2024; 109:015201. [PMID: 38366520 DOI: 10.1103/physreve.109.015201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2023] [Accepted: 12/07/2023] [Indexed: 02/18/2024]
Abstract
Magnetic fields influence ion transport in plasmas. Straightforward comparisons of experimental measurements with plasma theories are complicated when the plasma is inhomogeneous, far from equilibrium, or characterized by strong gradients. To better understand ion transport in a partially magnetized system, we study the hydrodynamic velocity and temperature evolution in an ultracold neutral plasma at intermediate values of the magnetic field. We observe a transverse, radial breathing mode that does not couple to the longitudinal velocity. The inhomogeneous density distribution gives rise to a shear velocity gradient that appears to be only weakly damped. This mode is excited by ion oscillations originating in the wings of the distribution where the plasma becomes non-neutral. The ion temperature shows evidence of an enhanced electron-ion collision rate in the presence of the magnetic field. Ultracold neutral plasmas provide a rich system for studying mode excitation and decay.
Collapse
Affiliation(s)
- Chanhyun Pak
- Department of Physics and Astronomy, Brigham Young University, Provo, Utah 84602, USA
| | - Virginia Billings
- Department of Physics and Astronomy, Brigham Young University, Provo, Utah 84602, USA
| | - Matthew Schlitters
- Department of Physics and Astronomy, Brigham Young University, Provo, Utah 84602, USA
| | - Scott D Bergeson
- Department of Physics and Astronomy, Brigham Young University, Provo, Utah 84602, USA
| | - Michael S Murillo
- Department of Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
3
|
Qiu P, Feng Y. Fast particles overtaking shock front in two-dimensional Yukawa solids. Phys Rev E 2022; 106:015203. [PMID: 35974640 DOI: 10.1103/physreve.106.015203] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2022] [Accepted: 07/04/2022] [Indexed: 06/15/2023]
Abstract
High-speed particles overtaking the shock front during the propagation of compressional shocks in two-dimensional (2D) Yukawa solids are investigated using molecular dynamical simulations. When the compressional speed is lower, all particles around the shock front are almost accelerated synchronously. However, when the compressional speed is much higher, some particles penetrate the shock front to enter the preshock region. Around the shock front, it is found that the particle velocity profile at the first peak of the dispersive shock wave (DSW) is able to be described using the Gaussian distribution, so that the amplitudes of the DSW can be well characterized. As the compressional speed increases, the particle velocity corresponding to these DSW's amplitudes increase more substantially than the shock front speed. These amplitudes of the DSW are found to be able to predict the occurrence of the fast particles. Combined with the previous study of the DSW's period, it is demonstrated that the properties of the DSW are nearly not affected by the conditions of the 2D Yukawa systems, but only related to the compressional speed.
Collapse
Affiliation(s)
- Pengwei Qiu
- Institute of Plasma Physics and Technology, School of Physical Science and Technology, Soochow University, Suzhou 215006, China
| | - Yan Feng
- Institute of Plasma Physics and Technology, School of Physical Science and Technology, Soochow University, Suzhou 215006, China
| |
Collapse
|
4
|
Sprenkle RT, Bergeson SD, Silvestri LG, Murillo MS. Ultracold neutral plasma expansion in a strong uniform magnetic field. Phys Rev E 2022; 105:045201. [PMID: 35590663 DOI: 10.1103/physreve.105.045201] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2021] [Accepted: 02/16/2022] [Indexed: 06/15/2023]
Abstract
In strongly magnetized neutral plasmas, electron motion is reduced perpendicular to the magnetic field direction. This changes dynamical plasma properties such as temperature equilibration, spatial density evolution, electron pressure, and thermal and electrical conductivity. In this paper we report measurements of free plasma expansion in the presence of a strong magnetic field. We image laser-induced fluorescence from an ultracold neutral Ca^{+} plasma to map the plasma size as a function of time for a range of magnetic field strengths. The asymptotic expansion velocity perpendicular to the magnetic field direction falls rapidly with increasing magnetic field strength. We observe that the initially Gaussian spatial distribution remains Gaussian throughout the expansion in both the parallel and perpendicular directions. We compare these observations with a diffusion model and with a self-similar expansion model and show that neither of these models reproduces the observed behavior over the entire range of magnetic fields used in this study. Modeling the expansion of a magnetized ultracold plasma poses a nontrivial theoretical challenge.
Collapse
Affiliation(s)
- R Tucker Sprenkle
- Department of Physics and Astronomy, Brigham Young University, Provo, Utah 84602, USA
| | - S D Bergeson
- Department of Physics and Astronomy, Brigham Young University, Provo, Utah 84602, USA
| | - Luciano G Silvestri
- Department of Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, Michigan 48824, USA
| | - Michael S Murillo
- Department of Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
5
|
Scheiner B, Baalrud SD. Viscosity of the magnetized strongly coupled one-component plasma. Phys Rev E 2021; 102:063202. [PMID: 33466065 DOI: 10.1103/physreve.102.063202] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2020] [Accepted: 11/10/2020] [Indexed: 11/07/2022]
Abstract
The viscosity tensor of the magnetized one-component plasma, consisting of five independent shear viscosity coefficients, a bulk viscosity coefficient, and a cross coefficient, is computed using equilibrium molecular dynamics simulations and the Green-Kubo relations. A broad range of Coulomb coupling and magnetization strength conditions are studied. Magnetization is found to strongly influence the shear viscosity coefficients when the gyrofrequency exceeds the Coulomb collision frequency. Three regimes are identified as the Coulomb coupling strength and magnetization strength are varied. The Green-Kubo relations are used to separate kinetic and potential energy contributions to each viscosity coefficient, showing how each contribution depends upon the magnetization strength. The shear viscosity coefficient associated with the component of the pressure tensor parallel to the magnetic field, and the two coefficients associated with the component perpendicular to the magnetic field, are all found to merge to a common value at strong Coulomb coupling.
Collapse
Affiliation(s)
- Brett Scheiner
- Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Scott D Baalrud
- Department of Physics and Astronomy, University of Iowa, Iowa City, Iowa 52242, USA
| |
Collapse
|
6
|
Scheiner B, Baalrud SD. Testing thermal conductivity models with equilibrium molecular dynamics simulations of the one-component plasma. Phys Rev E 2019; 100:043206. [PMID: 31770988 DOI: 10.1103/physreve.100.043206] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2019] [Indexed: 11/07/2022]
Abstract
Equilibrium molecular dynamics simulations are used to calculate the thermal conductivity of the one-component plasma via the Green-Kubo formalism over a broad range of Coulomb coupling strength, 0.1≤Γ≤180. These simulations address previous discrepancies between computations using equilibrium versus nonequilibrium methods. Analysis of heat flux autocorrelation functions show that very long (6×10^{5}ω_{p}^{-1}) time series are needed to reduce the noise level to allow ≲2% accuracy. The simulations provide accurate data for Γ≲1. This enables a test of the traditional Landau-Spitzer theory, which is found to agree with the simulations for Γ≲0.3. It also enables tests of theories to address moderate and strong Coulomb coupling. Two are found to provide accurate extensions to the moderate coupling regime of Γ≲10, but none are accurate in the Γ≳10 regime where potential energy transport and coupling between mass flow and stress dominate thermal conduction.
Collapse
Affiliation(s)
- Brett Scheiner
- Department of Physics and Astronomy, University of Iowa, Iowa City, Iowa 52240, USA
| | - Scott D Baalrud
- Department of Physics and Astronomy, University of Iowa, Iowa City, Iowa 52240, USA
| |
Collapse
|
7
|
Feng Y, Lin W, Murillo MS. Viscosity of two-dimensional strongly coupled dusty plasma modified by a perpendicular magnetic field. Phys Rev E 2018; 96:053208. [PMID: 29347770 DOI: 10.1103/physreve.96.053208] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2017] [Indexed: 11/07/2022]
Abstract
Transport properties of two-dimensional (2D) strongly coupled dusty plasmas have been investigated in detail, but never for viscosity with a strong perpendicular magnetic field; here, we examine this scenario using Langevin dynamics simulations of 2D liquids with a binary Yukawa interparticle interaction. The shear viscosity η of 2D liquid dusty plasma is estimated from the simulation data using the Green-Kubo relation, which is the integration of the shear stress autocorrelation function. It is found that, when a perpendicular magnetic field is applied, the shear viscosity of 2D liquid dusty plasma is modified substantially. When the magnetic field is increased, its viscosity increases at low temperatures, while at high temperatures its viscosity diminishes. It is determined that these different variational trends of η arise from the different behaviors of the kinetic and potential parts of the shear stress under external magnetic fields.
Collapse
Affiliation(s)
- Yan Feng
- Center for Soft Condensed Matter Physics and Interdisciplinary Research, College of Physics, Optoelectronics and Energy, Soochow University, Suzhou 215006, China
| | - Wei Lin
- Center for Soft Condensed Matter Physics and Interdisciplinary Research, College of Physics, Optoelectronics and Energy, Soochow University, Suzhou 215006, China
| | - M S Murillo
- Department of Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
8
|
Baalrud SD, Daligault J. Transport regimes spanning magnetization-coupling phase space. Phys Rev E 2017; 96:043202. [PMID: 29347622 DOI: 10.1103/physreve.96.043202] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2017] [Indexed: 06/07/2023]
Abstract
The manner in which transport properties vary over the entire parameter-space of coupling and magnetization strength is explored. Four regimes are identified based on the relative size of the gyroradius compared to other fundamental length scales: the collision mean free path, Debye length, distance of closest approach, and interparticle spacing. Molecular dynamics simulations of self-diffusion and temperature anisotropy relaxation spanning the parameter space are found to agree well with the predicted boundaries. Comparison with existing theories reveals regimes where they succeed, where they fail, and where no theory has yet been developed.
Collapse
Affiliation(s)
- Scott D Baalrud
- Department of Physics and Astronomy, University of Iowa, Iowa City, Iowa 52242, USA
| | - Jérôme Daligault
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| |
Collapse
|
9
|
Di Cintio P, Livi R, Lepri S, Ciraolo G. Multiparticle collision simulations of two-dimensional one-component plasmas: Anomalous transport and dimensional crossovers. Phys Rev E 2017; 95:043203. [PMID: 28505790 DOI: 10.1103/physreve.95.043203] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2016] [Indexed: 06/07/2023]
Abstract
By means of hybrid multiparticle collsion-particle-in-cell (MPC-PIC) simulations we study the dynamical scaling of energy and density correlations at equilibrium in moderately coupled two-dimensional (2D) and quasi-one-dimensional (1D) plasmas. We find that the predictions of nonlinear fluctuating hydrodynamics for the structure factors of density and energy fluctuations in 1D systems with three global conservation laws hold true also for 2D systems that are more extended along one of the two spatial dimensions. Moreover, from the analysis of the equilibrium energy correlators and density structure factors of both 1D and 2D neutral plasmas, we find that neglecting the contribution of the fluctuations of the vanishing self-consistent electrostatic fields overestimates the interval of frequencies over which the anomalous transport is observed. Such violations of the expected scaling in the currents correlation are found in different regimes, hindering the observation of the asymptotic scaling predicted by the theory.
Collapse
Affiliation(s)
- Pierfrancesco Di Cintio
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata "Nello Carrara" via Madonna del piano 10, I-50019 Sesto Fiorentino, Italy
- Dipartimento di Fisica e Astronomia and CSDC, Universitá di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
- Istituto Nazionale di Fisica Nucleare, Sezione di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
| | - Roberto Livi
- Dipartimento di Fisica e Astronomia and CSDC, Universitá di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
- Istituto Nazionale di Fisica Nucleare, Sezione di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
- Consiglio Nazionale delle Ricerche, Istituto dei Sistemi Complessi via Madonna del piano 10, I-50019 Sesto Fiorentino, Italy
| | - Stefano Lepri
- Istituto Nazionale di Fisica Nucleare, Sezione di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
- Consiglio Nazionale delle Ricerche, Istituto dei Sistemi Complessi via Madonna del piano 10, I-50019 Sesto Fiorentino, Italy
| | | |
Collapse
|
10
|
Ott T, Bonitz M, Hartmann P, Donkó Z. Spontaneous generation of temperature anisotropy in a strongly coupled magnetized plasma. Phys Rev E 2017; 95:013209. [PMID: 28208314 DOI: 10.1103/physreve.95.013209] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2016] [Indexed: 06/06/2023]
Abstract
A magnetic field was recently shown to enhance field-parallel heat conduction in a strongly correlated plasma whereas cross-field conduction is reduced. Here we show that in such plasmas, the magnetic field has the additional effect of inhibiting the isotropization process between field-parallel and cross-field temperature components, thus leading to the emergence of strong and long-lived temperature anisotropies when the plasma is locally perturbed. An extended heat equation is shown to describe this process accurately.
Collapse
Affiliation(s)
- T Ott
- Christian-Albrechts-University Kiel, Institute for Theoretical Physics and Astrophysics, Leibnizstraße 15, 24098 Kiel, Germany
| | - M Bonitz
- Christian-Albrechts-University Kiel, Institute for Theoretical Physics and Astrophysics, Leibnizstraße 15, 24098 Kiel, Germany
| | - P Hartmann
- Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, P.O. Box 49, H-1525 Budapest, Hungary
| | - Z Donkó
- Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, P.O. Box 49, H-1525 Budapest, Hungary
| |
Collapse
|