Bronin et al. [Phys. Rev. E 108, 045209 (2023)2470-004510.1103/PhysRevE.108.045209] recently reported molecular-dynamics simulations of ultracold neutral plasmas expanding in a quadrupole magnetic field. While the main results are in agreement with prior experimental measurements, we present data showing oscillations not captured in the simulations of Bronin et al. Plasmas formed using pulsed or continuous-wave ionization processes have similar confinement times.
Bronin et al. [Phys. Rev. E 108, 045209 (2023)2470-004510.1103/PhysRevE.108.045209] recently reported molecular-dynamics simulations of ultracold neutral plasmas expanding in a quadrupole magnetic field. While the main results are in agreement with prior experimental measurements, we present data showing oscillations not captured in the simulations of Bronin et al. Plasmas formed using pulsed or continuous-wave ionization processes have similar confinement times.
In this Reply, we respond to the Comment by Schlitters et al. on our recent work [Phys. Rev. E 108, 045209 (2023)10.1103/PhysRevE.108.045209], where we present simulation results of ultracold Sr plasma expansion in a quadrupole magnetic field using a molecular-dynamics method. In the Comment, Schlitters et al. present their experimental results, some of which, from their point of view, contradict our simulation results and others that confirm the experiments of Gorman et al. [Phys. Rev. Lett. 126, 085002 (2021)10.1103/PhysRevLett.126.085002] and our results. In addition, Schlitters et al. also provide results that were not described in our work. Here we show that there is not a contradiction but a misunderstanding of our results. This is due to the fact that simulations allow the use of processing methods that are difficult to obtain in experiment. We also present different simulation results related to expansion velocity variations that reflect the experimental results of Schlitters et al. and provide an explanation for the origin of these variations.
We present simulation results of ultracold Sr plasma expansion in a quadrupole magnetic field by means of molecular dynamics. An analysis of plasma evolution influenced by a magnetic field is given. Plasma confinement time behavior under variation of magnetic field strength is estimated. Similarity of the time dependence of the concentration and distribution of ion velocities against the parameters of the plasma and magnetic field is established. Simulation results are in agreement with the experimental ones.
We demonstrate magnetic confinement of an ultracold neutral plasma (UCNP) created at the null of a biconic cusp, or quadrupole magnetic field. Initially, the UCNP expands due to electron thermal pressure. As the plasma encounters stronger fields, expansion slows and the density distribution molds to the field. UCNP electrons are strongly magnetized over most of the plasma, while ion magnetization is only significant at the boundaries. Observations suggest that electrons and ions are predominantly trapped by magnetic mirroring and ambipolar electric fields, respectively. Confinement times approach 0.5 ms, while unmagnetized plasmas dissipate on a timescale of a few tens of microseconds.
Abstract The results of calculations using the molecular dynamics method of low-frequency ion oscillations in ultracold plasma, previously observed in the experiment, are presented. Oscillations of the radial velocity of ions and changes in the direction of their flows during the expansion of ultracold plasma in a quadrupole magnetic field are detected. This is due to the influence of the magnetic field on the charge imbalance, which leads to the simultaneous existence of potential wells for both electrons and ions. Ion temperature oscillations with a frequency close to the ion plasma frequency are also detected. Satisfactory agreement with the experiment is obtained.
In the recent work Fletcher et al. [Phys. Rev. Lett., v.99, p.145001 (2007)] reported on the novel experimental technique, enabling to measure the temperature of the expanding ultracold plasmas over a considerable time interval (up to 60-70 μs). It was unexpectedly found that the electron temperature dropped with time as T_e(t) ~ t^{-α} with α= 1.2 +/- 0.1 \approx 1 instead of α= 2, which would be expected for the adiabatic cooling of electrons in the cloud expanding linearly in time. The above-cited authors supposed that 'the difference is likely due to the significant heating effects from 3-body recombination' (i.e. the inelastic processes), but they did not provide sufficient quantitative estimates supporting such a conclusion. The aim of the present comment is to mention that the experimentally revealed t^{-1}-dependence can be explained under quite general assumptions by the purely elastic processes in the ultracold plasma, as it was done a few years ago in our work [Yu.V. Dumin, J. Low Temp. Phys., v.119, p.377 (2000)]. Therefore, t^{-1}-dependence revealed in the above-cited experiment should be primarily a manifestation of "virialization" of the electron velocities in the regime of strong electron-ion correlations, while the heat release by inelastic processes should be of secondary importance (it might be responsible, particularly, for changing the exponent from -1 to -1.2). This is confirmed also by our ab initio molecular-dynamic simulations taking into account the strong electron-ion correlations (i.e. not based on the PIC method, Vlasov approximation for electrons, etc., which are applicable only to small-angular scattering).
Plasmas, often called the fourth state of matter and the most common one in the universe, have parameters varying by many orders of magnitude, from temperature of a few hundred kelvin in the Earth's ionosphere to 10 16 K in the magnetosphere of a pulsar. Ultracold plasmas, produced by photoionizing a sample of laser-cooled and trapped atoms near the ionization limit, have extended traditional neutral plasma parameters by many orders of magnitude, to electron temperatures below 1 K and ion temperatures in the tens of &mu K to a few Kelvin, and densities of 10 5 cm -3 to 10 10 cm -3 . These plasmas thus provide a testing ground to study basic plasma theory in a clean and simple system with or without a magnetic field. Previous studies of ultracold plasmas have primarily concentrated on temperature measurements, collective modes and expansion dynamics in the absence of magnetic fields. This thesis presents the first study of ultracold plasma dynamics in a magnetic field. The presence of a magnetic field during the expansion can initiate various phenomena, such as plasma confinement and plasma instabilities. While the electron temperatures are very low in ultracold plasmas, we need only tens of Gauss of magnetic field to observe significant effects on the expansion dynamics. To probe the ultraocold plasma dynamics in a magnetic field, we developed a new diagnostic - projection imaging, which images the ion distribution by extracting the ions with a high voltage pulse onto a position-sensitive detector. Early in the lifetime of the plasma (< 20 &mu s), the size of the image is dominated by the time-of-flight Coulomb explosion of the dense ion cloud. For later times, we measure the 2-D Gaussian width of the ion image, obtaining the transverse expansion velocity as a function of magnetic field (up to 70 G),and observe that the transverse expansion velocity scales as B &minus1/2 , explained by a nonlinear ambipolar diffusion model that involes anisotropic diffusion in two different directions. We also present the first observation of a plasma instability in an expanding ultracold plasma. We observe periodic emission of electrons from an ultracold plasma in weak, crossed magnetic and electric fields, and a strong perturbed electron density distribution in electron time-of-flight projection images. We identify this instability as a high-frequency electron drift instability due to the coupling between the electron drift wave and electron cyclotron harmonic, which has large wavenumbers corresponding to wavelengths close to the electron gyroradius.