Search

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.

In this Letter we discuss our approach that makes possible creation of the steady-state ultracold plasma having various densities and temperatures by means of continuous two-step optical excitation of calcium atoms in the magneto-optical trap. A strongly coupled ultracold plasma can be used as an excellent test platform for studying many-body interactions associated with various plasma phenomena. The parameters of the plasma are studied using laser-induced fluorescence of calcium ions. The experimental results are well described by a simple theoretical model involving equilibration of the continuous source of charged particles by the hydrodynamical ion outflux and three-body recombination. The ultracold plasma with the peak ion density of 2.7×10^{6}  cm^{-3} and the minimum electron temperature near 2 K has been prepared. Our steady-state approach in combination with the strong magnetic confinement of the plasma will make it possible to reach extremely strong coupling in such system.

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 present the results of a direct simulation of the expansion of a two-component ultracold plasma for various numbers of particles, densities, and electron temperatures. A description of the expansion process common to all plasma parameters is given. After the escape of fast electrons from the plasma cloud, the excess positive charge is localized at the outer boundary, in a narrow layer. This layer has a characteristic front shape with a sharp drop in the charge concentration. The charged layer retains the remaining electrons during the entire expansion process. As the plasma expands, the speed of movement of the charged layer becomes constant and significantly exceeds the sonic speed of ions. In addition, the dependence of the radial velocity of ions on the radius acquires a self-similar character long before the final stage of expansion. Based on the calculation results, equations and self-similar solutions are obtained. General dependencies on plasma parameters are determined, which are compared with experimental data.

A two-component plasma model, which we called a "shelf Coulomb" model has been developed in this work. A Monte Carlo study has been undertaken to calculate equations of state, pair distribution functions, internal energies, and other thermodynamics properties. A canonical NVT ensemble with periodic boundary conditions was used. The motivation behind the model is also discussed in this work. The "shelf Coulomb" model can be compared to classical two-component (electron-proton) model where charges with zero size interact via a classical Coulomb law. With important difference for interaction of opposite charges: electrons and protons interact via the Coulomb law for large distances between particles, while interaction potential is cut off on small distances. The cut off distance is defined by an arbitrary ɛ parameter, which depends on system temperature. All the thermodynamics properties of the model depend on dimensionless parameters ɛ and γ = βe(2)n(1/3) (where β = 1/kBT, n is the particle's density, kB is the Boltzmann constant, and T is the temperature) only. In addition, it has been shown that the virial theorem works in this model. All the calculations were carried over a wide range of dimensionless ɛ and γ parameters in order to find the phase transition region, critical point, spinodal, and binodal lines of a model system. The system is observed to undergo a first order gas-liquid type phase transition with the critical point being in the vicinity of ɛ(crit) ≈ 13(T(*)(crit) ≈ 0.076), γ(crit) ≈ 1.8(v(*)(crit) ≈ 0.17), P(*)(crit) ≈ 0.39, where specific volume v* = 1/γ(3) and reduced temperature T(*) = ɛ(-1).