Ultracold Plasma Research Articles

Theoretical Investigation of Rate Coefficients and Dynamical Mechanisms for N + N + N Three-Body Recombination Based on Full-Dimensional Potential Energy Surfaces

Three-body recombination reactions, in which two particles form a bound state while a third one bounces off after the collision, play significant roles in many fields, such as cold and ultracold chemistry, astrochemistry, atmospheric physics, and plasma physics. In this work, the dynamics of the recombination reaction for the N3 system over a wide temperature range (5000–20,000 K) are investigated in detail using the quasi-classical trajectory (QCT) method based on recently developed full-dimensional potential energy surfaces. The recombination products are N2(X) + N(4S) in the 14A″ state, N2(A) + N(4S) in the 24A″ state, and N2(X) + N(2D) in both the 12A″ and 22A″ states. A three-body collision recombination model involving two sets of relative translational energies and collision parameters and a time-delay parameter is adopted in the QCT calculations. The recombination process occurs after forming an intermediate with a certain lifetime, which has a great influence on the recombination probability. Recombination processes occurring through a one-step three-body collision mechanism and two distinct two-step binary collision mechanisms are found in each state. And the two-step exchange mechanism is more dominant than the two-step transfer mechanism at higher temperatures. N2(X) formed in all three related states is always the major recombination product in the temperature range from 5000 K to 20,000 K, with the relative abundance of N2(A) increasing as temperature decreases. After hyperthermal collisions, the formed N2(X/A) molecules are distributed in highly excited rotational and vibrational states, with internal energies mainly distributed near the dissociation threshold. Additionally, the rate coefficients for this three-body recombination reaction in each state are determined and exhibit a negative correlation with temperature. The dynamic insights presented in this work might be very useful to further simulate non-equilibrium dynamic processes in plasma physics involving N3 systems.

N-changing population of Rydberg states by low-energy electron-Rydberg collisions.

Inelastic n-changing collisions play an important role in the evolution of Rydberg atoms into ultracold plasmas. However, for the initially intermediate n (n ∼ 40) Rydberg states, these collisions can hardly be observed due to the low electron temperature in ultracold plasmas. In this work, we designed an experimental scheme to facilitate collisions between free electrons at 1.5eV and intermediate n Rydberg atoms. Using the field ionization technique, we measured the state distributions resulting from the evolution of initially cold rubidium atoms in the 45P3/2 Rydberg state. The experimentally obtained probability of inelastic collisions excitation agrees well with the Monte Carlo simulation results. In addition, our experimental results indicate that the n-changing population induced by hot electrons is significant for lower nP Rydberg states. Our work plays a significant role in calculating the rates of electron-ion three-body recombination in ultracold plasmas.

Evolution of shielding cloud under oscillatory external forcing in strongly coupled ultracold neutral plasma

This paper investigates the dynamics of crystalline clusters observed in Molecular Dynamics (MD) studies conducted earlier (Yadav et al., 2023) for ultra-cold neutral plasmas. An external oscillatory forcing is applied for this purpose, and the evolution is tracked with the help of MD simulations using the open-source LAMMPS software. Interesting observations relating to cluster dynamics are presented. The formation of a pentagonal arrangement of particles is also reported.

Reply to "Comment on 'Ultracold plasma expansion in quadrupole magnetic field'".

In this Reply, we respond to the Comment by Schlitters etal. 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 etal. present their experimental results, some of which, from their point of view, contradict our simulation results and others that confirm the experiments of Gorman etal. [Phys. Rev. Lett. 126, 085002 (2021)10.1103/PhysRevLett.126.085002] and our results. In addition, Schlitters etal. 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 etal. and provide an explanation for the origin of these variations.

Comment on "Ultracold plasma expansion in quadrupole magnetic field".

Bronin etal. [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 etal. Plasmas formed using pulsed or continuous-wave ionization processes have similar confinement times.

Full spectroscopic model and trihybrid experimental-perturbative-variational line list for NH

ABSTRACT Imidogen (NH) is a reactive molecule whose presence in astrochemical environments is of interest due to its role in the formation of nitrogen-containing molecules and as a potential probe of nitrogen abundance. Spectroscopic NH monitoring is useful for Earth-based combustion and photolysis processes of ammonia and other nitrogen-containing species. NH is also relevant to ultracold molecular physics and plasma studies. To enable these diverse applications, high-quality molecular spectroscopic data are required. Here, a new line list with significant advantages over existing data is presented. Most notably, this line list models isotopologue spectroscopy and forbidden transitions (important for NH visible absorption), alongside some overall improvements to accuracy and completeness. This approach takes advantage of existing experimental data (from a previous MARVEL compilation) and perturbative line lists together with new MRCI ab initio electronic data. These are used to produce a novel variational spectroscopic model and trihybrid line list for the main 14N1H isotopologue, as well as isotopologue-extrapolated hybrid line lists for the 14N2H, 15N1H, and 15N2H isotopologues. The new 14N1H ExoMol-style trihybrid line list, kNigHt, comprises 4076 energy levels (1078 experimental) and 327 014 transitions up to 47 500 cm−1 (211 nm) between five low-lying electronic states (X 3Σ−, a 1Δ, b 1Σ+, A 3Π, and c 1Π). For most anticipated applications aside from far-infrared studies, this line list will be of sufficient quality; any improvements should focus on the b 1Σ+ energies, and the a 1Δ – A 3Π and b 1Σ+ – A 3Π spin–orbit couplings.

Suppressing electron disorder-induced heating of ultracold neutral plasma via optical lattice

Disorder-induced heating (DIH) prevents ultracold neutral plasma into the electron strong coupling regime. Here, we propose a scheme to suppress electronic DIH via optical lattice. We simulate the evolution dynamics of ultracold neutral plasma constrained by three-dimensional optical lattice using the classical molecular dynamics method. The results show that for experimentally achievable conditions, electronic DIH is suppressed by a factor of 1.3, and the Coulomb coupling strength can reach 0.8, which is approaching the strong coupling regime. Suppressing electronic DIH via optical lattice may pave a way for the research of electronic strongly coupled plasma.

Steady-State Ultracold Plasma Created by Continuous Photoionization of Laser Cooled Atoms.

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 2K 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.

Electronic resonances in expanding non-neutral ultracold plasma

Calculations of the spectrum of eigen oscillations of an inhomogeneous ultracold plasma are presented both in the absence and in the presence of charge imbalance. The collective modes of these plasma oscillations are recorded in experiments as absorption resonances of a radio frequency electric field. It is shown that in the presence of friction of the electronic component, a discrete spectrum of plasma eigen oscillations is formed. The dependence of the frequency of these resonances on the plasma expansion time in the presence of a charge imbalance is obtained. There is good agreement with the experiments of different authors.

Calibrated heating rate measurements using electric-field-induced electron extraction in ultracold neutral plasmas

The heating rate of plasma electrons induced by external fields or other processes can be used as an experimental tool to measure fundamental plasma properties such as electrical conductivity or electron–ion collision rates. We have developed a technique that can measure electron heating rates in ultracold neutral plasmas (UNPs) with $\sim 10\,\%$ precision while simultaneously referencing the measurement to a calibrated amount of heating. This technique uses a sequence of applied electric fields in four sections: to control the ratio of electrons to ions in the UNP; to provide a time for the application of fields that cause electron heating and subsequent thermalization of the electrons after the application of those fields; to extract electrons from the UNP using a method sensitive to electron temperature that allows the measurement of electron heating; and to extract the remaining electrons to measure the total electron (and therefore ion) number. The primary signal used to measure the heating rate is the measurement of the number of electrons that escape in the third section of the experiment as a larger number of escaping electrons indicates a larger amount of heating. We illustrate the use of this technique by measuring electron heating caused by high-frequency radiofrequency (RF) fields. In addition to the main technique, several subtechniques to calibrate the electron temperature, electron density, amount of heating and applied RF field amplitude were developed as well.

Self-organization in the avalanche, quench and dissipation of a molecular ultracold plasma

Spontaneous avalanche to plasma begins in the core of an ellipsoidal Rydberg gas of nitric oxide. Ambipolar expansion of NO $^+$ draws energy from avalanche-heated electrons. Then, cycles of long-range resonant electron transfer from Rydberg molecules to ions equalize their relative velocities. This sequence of steps gives rise to a remarkable mechanics of self-assembly, in which the kinetic energy of initially formed hot electrons and ions drives an observed separation of plasma volumes. These dynamics adiabatically sequester energy in a reservoir of mass transport, starting a process that anneals separating volumes to form an apparent glass of strongly coupled ions and electrons. Short-time electron spectroscopy provides experimental evidence for complete ionization. The long lifetime of this system, particularly its stability with respect to recombination and neutral dissociation, suggests that this transformation affords a robust state of arrested relaxation, far from thermal equilibrium. We see this most directly in the excitation spectrum of transitions to states in the initially selected Rydberg series, detected as the long-lived signal that survives a flight time of $500\ \mathrm {\mu }$ s to reach an imaging detector. The initial density of electrons produced by prompt Penning ionization, which varies with the selected initial principal quantum number and density of the Rydberg gas, determines a balance between the rising density of ions and the falling density of Rydberg molecules. This Penning-regulated ion-Rydberg molecule balance appears necessary as a critical factor in achieving the long ultracold plasma lifetime to produce spectral features detected after very long delays.

Effects of disorder on electron heating in ultracold plasmas

Starting from the beginning of their research in the early 2000s, the ultracold plasmas were considered as a promising tool to achieve considerable values of the Coulomb coupling parameter for electrons. Unfortunately, this was found to be precluded by a sharp spontaneous increase in temperature, which was often attributed to the so-called disorder-induced heating. It is the aim of the present paper to quantify the effect of spontaneous heating as a function of the initial ionic disorder and, thereby, to estimate the efficiency of its mitigation, e.g., by the Rydberg blockade. As a result of the performed simulations, we found that the dynamics of electrons exhibited a well-expressed transition from the case of the quasi-regular arrangement of ions to the disordered one; the magnitude of the effect being about 30%. Thereby, we can conclude that the two-step formation of ultracold plasmas—involving the intermediate stage of the blockaded Rydberg gas—can really serve as a tool to increase the degree of Coulomb coupling, but the efficiency of this method is moderate.

Ultracold plasma expansion in quadrupole magnetic field.

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.

Velocity-mapped imaging of electron dynamics in an ultracold laser-induced plasma

Velocity-mapped imaging of electron dynamics in an ultracold laser-induced plasma

Ion dynamics in strongly coupled ultracold neutral plasmas at the early stage

We measured the velocity distribution of ions in the ultracold neutral plasma (UNP) during its early evolution by the technique of velocity map imaging. The ion temperature during the phase of ion equilibrium was obtained experimentally through this method. The Coulomb coupling parameter of ions in the UNP after disorder-induced heating was also determined to be 2.1, which agreed well with the prediction from a charged particle tracing simulation. In addition, the ion expansion during the ion equilibration phase was observed. Notably, the experimentally observed expansion speed is larger than the value obtained from the self-similar expansion model, indicating the involvement of additional mechanisms, besides the electron thermal pressure, in driving the ion expansion. We have also discussed the contributions of ion–ion correlations and charge imbalance to the plasma expansion.

Charged particle diagnostic system for lithium ultracold plasma experiments

Charged particle detection systems are widely used in modern atomic physics experiments. In this paper we describe an experimental setup for ultracold plasma (UCP) diagnostics based on the basis of a microchannel plate-based detector (MCP). Two different readout techniques for the MCP output are described in detail. Expansion of the low-density UCP and the frequency-dependent response of the UCP to an external electric field are measured. Preliminary results are obtained.

Magnetic Field Effect on the Formation of Ultracold Plasma

The molecular dynamics method is used to simulate the formation of ultracold plasma under continuous ionizing irradiation in a quadrupole magnetic field with the gradient of the magnetic field along the axis of symmetry equal to 0, 30, 150, and 500 G/cm. An increase in the magnetic field promotes an increase in the plasma density owing to the trapping of some part of fast electrons by the quadrupole magnetic field.

Analysis and simulation of time-of-flight spectrum in Rb<sup>+</sup>-Rb hybrid trap

The time-of-flight mass spectrum of charged particles, which are created through two-step cw-laser photoionization of laser-cooled <sup>87</sup>Rb atoms in an ion-neutral hybrid trap, is quantitatively investigated to further facilitate the study of Rb<sup>+</sup>-Rb reactive collisions. A microchannel plate (MCP) is used to detect charged particles, and two spectral peaks corresponding to the <sup>87</sup>Rb<sup>+</sup> ions and the product <inline-formula><tex-math id="Z-20230423060254-1">\begin{document}$ \rm {}^{87}Rb_2^+ $\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="9-20222273_Z-20230423060254-1.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="9-20222273_Z-20230423060254-1.png"/></alternatives></inline-formula> of the Rb<sup>+</sup>-Rb reaction were observed in the time-of-flight spectrum, respectively. The two peaks overlapped with each other and both showed an asymmetric profile. The information about the intensity, position, and half-width of the peak for a specific ion species was derived by fitting the time-of-flight spectrum with the probability density function of the Gumbel distribution. Then the relative ion intensity was converted into absolute ion number through the following steps. The rate equation of the total number of ions, which includes the number of atoms, the calibration factor of MCP, and the effective decay rate of ions in the ion trap, was established by modeling the photoionization of atoms. Combined with the absolute number of atoms measured by absorption imaging, the calibration factor in converting the ion intensity into the ion number was derived and the relative ion intensity was converted into the absolute number of ions. This provides a method of calibrating the MCP. The reliability of our calibration method was proved by the fact that the calibration factor in converting the intensity measured by MCP into particle number is independent of the duration of photoionization, the intensity and wavelength of the ionizing laser. Moreover, in order to explain the relationship between the peak width and temperature of the corresponding ion species, the time-of-flight spectra of the ions trapped in the ion trap were simulated by using COMSOL Multiphysics. The simulation results demonstrated that the large ion kinetic energy results in a narrow spectral peak. In sum, we quantitatively analyze and simulate the time-of-flight spectrum of the photoionization of cold atoms in the Rb<sup>+</sup>-Rb hybrid trap. The absolute number of ions is obtained by the intensity of the spectral peak, and the width of the spectral peak is related to the temperature of the ions. These results lay a foundation for the in-depth analysis of the ion-atom reaction collision and charged particle temperature relaxation in the photoionization of cold atoms, and thus further elucidating the subsequent collisional dynamics of ultracold plasmas.

General expression for the probability distribution function of electric field in a spatially inhomogeneous non-neutral plasma

We discuss the influence of micro- and macro-fields on spectral lines of ions as it takes place for spatially inhomogeneous plasma. A distribution function of an electric field is obtained. The function accounts for inhomogeneity and non-neutrality of plasma. The results of calculations of this function for various regimes are presented. Experimental results for ultracold plasma are used to compare theory with experiment. Dependence of the absorption coefficient on the function is shown. These results may be useful for diagnostics of various types of plasmas. One of the methods of plasma diagnostics is the analysis of the influence of its electric field on the shape of the spectral lines of atoms and ions (the Stark effect).

Photoionization-cross-section measurement of the Rb85 5 P3/22 excited state using microwave cavity resonance spectroscopy

We present microwave cavity resonance spectroscopy as a technique to determine excited-state photoionization cross sections of laser-cooled atoms. We demonstrate this technique by measuring the photoionization cross section of the $5^{2}P_{3/2}({F}^{\ensuremath{'}}=4)$ excited state of $^{85}\mathrm{Rb}$ to continuum by creating an ultracold plasma inside a 5-GHz resonant microwave cavity. We find a photoionization cross section of $\ensuremath{\sim}5\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}22}\phantom{\rule{4pt}{0ex}}{\mathrm{m}}^{2}$ for photon excess energies of 50, 100, 200, and 500 K above the ionization threshold. The measurements yield a cross section which is approximately 2 to 3 times smaller than the value provided by existing theory for the $5p$ excited state of $^{85}\mathrm{Rb}$ and two earlier measurements performed with other techniques for similar excess energies.