No Arabic abstract
The friction force on a test particle traveling through a plasma that is both strongly coupled and strongly magnetized is studied using molecular dynamics simulations. In addition to the usual stopping power component aligned antiparallel to the velocity, a transverse component that is perpendicular to both the velocity and Lorentz force is observed. This component, which was recently discovered in weakly coupled plasmas, is found to increase in both absolute and relative magnitude in the strongly coupled regime. Strong coupling is also observed to induce a third component of the friction force in the direction of the Lorentz force. These first-principles simulations reveal novel physics associated with collisions in strongly coupled, strongly magnetized, plasmas that are not predicted by existing kinetic theories. The effect is expected to influence macroscopic transport in a number of laboratory experiments and astrophysical plasmas.
Plasmas that are strongly magnetized in the sense that the gyrofrequency exceeds the plasma frequency exhibit novel transport properties that are not well understood. As a representative example, we compute the friction force acting on a massive test charge moving through a strongly coupled and strongly magnetized one-component plasma using a generalized Boltzmann kinetic theory. Recent works studying the weakly coupled regime have shown that strong magnetization leads to a transverse component of the friction force that is perpendicular to both the Lorentz force and velocity of the test charge; in addition to the stopping power component aligned antiparallel to the velocity. Recent molecular dynamics simulations have also shown that strong Coulomb coupling in addition to strong magnetization gives rise to a third component of the friction force in the direction of the Lorentz force. Here, we show that the generalized Boltzmann kinetic theory captures these effects, and generally agrees well with the molecular dynamics simulations over a broad range of Coulomb coupling and magnetization strength regimes. The theory is also used to show that a gyro component of the friction in the direction of the Lorentz force arises due to asymmetries associated with gyromotion during short-range collisions. Computing the average motion of the test charge through the background plasma, the transverse force is found to strongly influence the trajectory by changing the gyroradius and the gyro friction force is found to slightly change the gyrofrequency of the test charge resulting in a phase shift.
Coulomb collisions in plasmas are typically modeled using the Boltzmann collision operator, or its variants, which apply to weakly magnetized plasmas in which the typical gyroradius of particles significantly exceeds the Debye length. Conversely, ONeil has developed a kinetic theory to treat plasmas that are so strongly magnetized that the typical gyroradius of particles is much smaller than the distance of closest approach in a binary collision. Here, we develop a generalized collision operator that applies across the full range of magnetization strength. To demonstrate novel physics associated with strong magnetization, it is used to compute the friction force on a massive test charge. In addition to the traditional stopping power component, this is found to exhibit a transverse component that is perpendicular to both the velocity and Lorentz force vectors in the strongly magnetized regime, as was predicted recently using linear response theory. Good agreement is found between the collision theory and linear response theory in the regime in which both apply, but the new collision theory also applies to stronger magnetization strength regimes than the linear response theory is expected to apply in.
A generalized Ohms law is derived to treat strongly magnetized plasmas in which the electron gyrofrequency significantly exceeds the electron plasma frequency. The frictional drag due to Coulomb collisions between electrons and ions is found to shift, producing an additional transverse resistivity term in the generalized Ohms law that is perpendicular to both the current ($vc{J}$) and the Hall ($vc{J} times vc{B}$) direction. In the limit of very strong magnetization, the parallel resistivity is found to increase by a factor of 3/2, and the perpendicular resistivity to scale as $ln (omega_{ce} tau_e)$, where $omega_{ce} tau_e$ is the Hall parameter. Correspondingly, the parallel conductivity coefficient is reduced by a factor of 2/3, and the perpendicular conductivity scales as $ln(omega_{ce} tau_e)/(omega_{ce} tau_e)^2$. These results suggest that strong magnetization significantly changes the magnetohydrodynamic evolution of a plasma.
We analyze how the turbulent transport of $mathbf{E}times mathbf{B}$ type in magnetically confined plasmas is affected by intermittent features of turbulence. The latter are captured by the non-Gaussian distribution $P(phi)$ of the turbulent electric potential $phi$. Our analysis is performed at an analytical level and confirmed numerically using two statistical approaches. We have found that the diffusion is inhibited linearly by intermittency, mainly via the kurtosis of the distribution $P(phi)$. The associated susceptibility for this linear process is shown to be dependent on the poloidal velocity $V_p$ and on the correlation time $tau_c$ with a maxima at the time-of-flight $tau_{fl}$. Intermittency does not affect the Kubo number scaling in the strong regime.
Molecular dynamics simulations are used to show that strong magnetization significantly increases the space and time scales associated with interparticle correlations. The physical mechanism responsible is a channeling effect whereby particles are confined to move along narrow cylinders with a width characterized by the gyroradius and a length characterized by the collision mean free path. The predominant interaction is $180^circ$ collisions at the ends of the collision cylinders, resulting in a long-range correlation parallel to the magnetic field. Its influence is demonstrated via the dependence of the velocity autocorrelation functions and self-diffusion coefficients on the domain size and run time in simulations of the one-component plasma. A very large number of particles, and therefore domain size, must be used to resolve the long-range correlations, suggesting that the number of charged particles in the collection must increase in order to constitute a plasma. Correspondingly, this effect significantly delays the time it takes to reach a diffusive regime, in which the mean square displacement of particles increases linearly in time. This result presents challenges for connecting measurements in non-neutral and ultracold neutral plasma experiments, as well as molecular dynamics simulations, with fluid transport properties due to their finite size.