Ultrafast dynamics of strongly coupled plasmas

Temperature relaxation in strongly-coupled binary ionic mixtures

New facilities such as the National Ignition Facility and the Linac Coherent Light Source have pushed the frontiers of high energy-density matter. These facilities offer unprecedented opportunities for exploring extreme states of matter, ranging from cryogenic solid-state systems to hot, dense plasmas, with applications to inertial-confinement fusion and astrophysics. However, significant gaps in our understanding of material properties in these rapidly evolving systems still persist. In particular, non-equilibrium transport properties of strongly-coupled Coulomb systems remain an open question. Here, we study ion-ion temperature relaxation in a binary mixture, exploiting a recently-developed dual-species ultracold neutral plasma. We compare measured relaxation rates with atomistic simulations and a range of popular theories. Our work validates the assumptions and capabilities of the simulations and invalidates theoretical models in this regime. This work illustrates an approach for precision determinations of detailed material properties in Coulomb mixtures across a wide range of conditions.

Temperature Equilibration due to Charge State Fluctuations in Dense Plasmas

The charge states of ions in dense plasmas fluctuate due to collisional ionization and recombination. Here, we show how, by modifying the ion interaction potential, these fluctuations can mediate energy exchange between the plasma electrons and ions. Moreover, we develop a theory for this novel electron-ion energy transfer mechanism. Calculations using a random walk approach for the fluctuations suggest that the energy exchange rate from charge state fluctuations could be comparable to direct electron-ion collisions. This mechanism is, however, predicted to exhibit a complex dependence on the temperature and ionization state of the plasma, which could contribute to our understanding of significant variation in experimental measurements of equilibration times.

Sculpted ultracold neutral plasmas

Ultracold neutral plasma (UNP) experiments allow for careful control of plasma properties across Coulomb coupling regimes. Here, we examine how UNPs can be used to study heterogeneous, nonequilibrium phenomena, including nonlinear waves, transport, hydrodynamics, kinetics, stopping power, and instabilities. Through a series of molecular dynamics simulations, we have explored UNPs formed with spatially modulated ionizing radiation. We have developed a computational model for such sculpted UNPs that includes an ionic screened Coulomb interaction with a spatiotemporal screening length, and Langevin-based spatial ion-electron and ion-neutral collisions. We have also developed a hydrodynamics model and have extracted its field quantities (density, flow velocity, and temperature) from the molecular dynamics simulation data, allowing us to investigate kinetics by examining moment ratios and phase-space dynamics; we find that it is possible to create UNPs that vary from nearly perfect fluids (Euler limit) to highly kinetic plasmas. We have examined plasmas in three geometries: a solid rod, a hollow rod, and a gapped slab; we have studied basic properties of these plasmas, including the spatial Coulomb coupling parameter. By varying the initial conditions, we find that we can design experimental plasmas that would allow the exploration of a wide range of phenomena, including shock and blast waves, stopping power, two-stream instabilities, and much more. Using an evaporative cooling geometry, our results suggest that much larger Coulomb couplings can be achieved, possibly in excess of 10.

Ion friction at small values of the Coulomb logarithm

Transport properties of high-energy-density plasmas are influenced by the ion collision rate. Traditionally, this rate involves the Coulomb logarithm, lnΛ. Typical values of lnΛ are ≈10-20 in kinetic theories where transport properties are dominated by weak-scattering events caused by long-range forces. The validity of these theories breaks down for strongly coupled plasmas, when lnΛ is of order one. We present measurements and simulations of collision data in strongly coupled plasmas when lnΛ is small. Experiments are carried out in the first dual-species ultracold neutral plasma (UNP), using Ca^{+} and Yb^{+} ions. We find strong collisional coupling between the different ion species in the bulk of the plasma. We simulate the plasma using a two-species fluid code that includes Coulomb logarithms derived from either a screened Coulomb potential or a the potential of mean force. We find generally good agreement between the experimental measurements and the simulations. With some improvements, the mixed Ca^{+} and Yb^{+} dual-species UNP will be a promising platform for testing theoretical expressions for lnΛ and collision cross-sections from kinetic theories through measurements of energy relaxation, stopping power, two-stream instabilities, and the evolution of sculpted distribution functions in an idealized environment in which the initial temperatures, densities, and charge states are accurately known.

Laser cooling of ions in a neutral plasma

Laser cooling of a neutral plasma is a challenging task because of the high temperatures typically associated with the plasma state. By using an ultracold neutral plasma created by photoionization of an ultracold atomic gas, we avoid this obstacle and demonstrate laser cooling of ions in a neutral plasma. After 135 microseconds of cooling, we observed a reduction in ion temperature by up to a factor of four, with the temperature reaching as low as 50(4) millikelvin. This pushes laboratory studies of neutral plasmas deeper into the strongly coupled regime, beyond the limits of validity of current kinetic theories for calculating transport properties. The same optical forces also retard the plasma expansion, opening avenues for neutral-plasma confinement and manipulation.

Molecular dynamics studies of electron-ion temperature equilibration in hydrogen plasmas within the coupled-mode regime

We use classical molecular dynamics (MD) to study electron-ion temperature equilibration in two-component plasmas in regimes for which the presence of coupled collective modes has been predicted to substantively reduce the equilibration rate. Guided by previous kinetic theory work, we examine hydrogen plasmas at a density of n=10^{26}cm^{-3}, T_{i}=10^{5}K, and 10^{7}K<T_{e}<10^{9}K. The nonequilibrium classical MD simulations are performed with interparticle interactions modeled by quantum statistical potentials (QSPs). Our MD results indicate (i) a large effect from time-varying potential energy, which we quantify by appealing to an adiabatic two-temperature equation of state, and (ii) a notable deviation in the energy equilibration rate when compared to calculations from classical Lenard-Balescu theory including the QSPs. In particular, it is shown that the energy equilibration rates from MD are more similar to those of the theory when coupled modes are neglected. We suggest possible reasons for this surprising result and propose directions of further research along these lines.

Ultracold neutral plasmas

By photoionizing samples of laser-cooled atoms with laser light tuned just above the ionization limit, plasmas can be created with electron and ion temperatures below 10 K. These ultracold neutral plasmas have extended the temperature bounds of plasma physics by two orders of magnitude. Table-top experiments, using many of the tools from atomic physics, allow for the study of plasma phenomena in this new regime with independent control over the density and temperature of the plasma through the excitation process. Characteristic of these systems is an inhomogeneous density profile, inherited from the density distribution of the laser-cooled neutral atom sample. Most work has dealt with unconfined plasmas in vacuum, which expand outward at velocities of order 100 m/s, governed by electron pressure, and with lifetimes of order 100 μs, limited by stray electric fields. Using detection of charged particles and optical detection techniques, a wide variety of properties and phenomena have been observed, including expansion dynamics, collective excitations in both the electrons and ions, and collisional properties. Through three-body recombination collisions, the plasmas rapidly form Rydberg atoms, and clouds of cold Rydberg atoms have been observed to spontaneously avalanche ionize to form plasmas. Of particular interest is the possibility of the formation of strongly coupled plasmas, where Coulomb forces dominate thermal motion and correlations become important. The strongest impediment to strong coupling is disorder-induced heating, a process in which Coulomb energy from an initially disordered sample is converted into thermal energy. This restricts electrons to a weakly coupled regime and leaves the ions barely within the strongly coupled regime. This review will give an overview of the field of ultracold neutral plasmas, from its inception in 1999 to current work, including efforts to increase strong coupling and effects on plasma properties due to strong coupling.

Demonstrating universal scaling for dynamics of Yukawa one-component plasmas after an interaction quench

The Yukawa one-component plasma (OCP) model is a paradigm for describing plasmas that contain one component of interest and one or more other components that can be treated as a neutralizing, screening background. In appropriately scaled units, interactions are characterized entirely by a screening parameter, κ. As a result, systems of similar κ show the same dynamics, regardless of the underlying parameters (e.g., density and temperature). We demonstrate this behavior using ultracold neutral plasmas (UNPs) created by photoionizing a cold (T≤10 mK) gas. The ions in UNP systems are well described by the Yukawa model, with the electrons providing the screening. Creation of the plasma through photoionization can be thought of as a rapid quench of the interaction potential from κ=∞ to a final κ value set by the electron density and temperature. We demonstrate experimentally that the postquench dynamics are universal in κ over a factor of 30 in density and an order of magnitude in temperature. Results are compared with molecular-dynamics simulations. We also demonstrate that features of the postquench kinetic energy evolution, such as disorder-induced heating and kinetic-energy oscillations, can be used to determine the plasma density and the electron temperature.

Strongly-coupled plasmas formed from laser-heated solids

We present an analysis of ion temperatures in laser-produced plasmas formed from solids with different initial lattice structures. We show that the equilibrium ion temperature is limited by a mismatch between the initial crystallographic configuration and the close-packed configuration of a strongly-coupled plasma, similar to experiments in ultracold neutral plasmas. We propose experiments to demonstrate and exploit this crystallographic heating in order to produce a strongly coupled plasma with a coupling parameter of several hundred.

Using higher ionization states to increase Coulomb coupling in an ultracold neutral plasma

We report measurements and simulations of the time-evolving rms velocity distribution in an ultracold neutral plasma. A strongly coupled ultracold neutral Ca+ plasma is generated by photoionizing laser-cooled atoms close to threshold. A fraction of these ions is then promoted to the second ionization state to form a mixed Ca+-Ca2+ plasma. By varying the time delay between the first and the second ionization events, a minimum in ion heating is achieved. We show that the Coulomb strong-coupling parameter Γ increases by a factor of 1.4 to a maximum value of 3.6. A pure Ca2+ plasma would have Γ=6.8, moving these strongly coupled systems closer to the regime of liquid-like correlations.

"Ultracold" neutral plasmas at room temperature

We report a measurement of the electron temperature in a plasma generated by a high-intensity laser focused into a jet of neon. The 15 eV electron temperature is determined using an analytic solution of the plasma equations assuming local thermodynamic equilibrium, initially developed for ultracold neutral plasmas. We show that this analysis method accurately reproduces more sophisticated plasma simulations in our temperature and density range. While our plasma temperatures are far outside the typical "ultracold" regime, the ion temperature is determined by the plasma density through disorder-induced heating just as in ultracold neutral plasma experiments. Based on our results, we outline a pathway for achieving a strongly coupled neutral laser-produced plasma that even more closely resembles ultracold neutral plasma conditions.

Ion acoustic waves in ultracold neutral plasmas

We photoionize laser-cooled atoms with a laser beam possessing spatially periodic intensity modulations to create ultracold neutral plasmas with controlled density perturbations. Laser-induced fluorescence imaging reveals that the density perturbations oscillate in space and time, and the dispersion relation of the oscillations matches that of ion acoustic waves, which are long-wavelength, electrostatic, density waves.

Molecular dynamics simulation of disorder-induced heating in ultracold neutral plasma

Disorder-induced heating (DIH) is one of the main reasons reducing the coupling strength in ultracold plasma. We propose applying an optical lattice as periodic confinement before the ultracold atomic cloud is ionized to eliminate its effect. We demonstrate a numerical simulation for the dynamics of the ultracold plasmas using classical molecular dynamics method with open boundary. DIH is reproduced in the simulation with the random Gaussian initial distribution and is absent in the results with the ordered lattice initial distribution. We further find that the collisional heating from electrons is important for ultracold plasmas with chosen spatial correlations in the optical lattice. Carefully preparing the initial condition (e.g., the ion density, initial electron temperature, and so on), collisional heating for the ions would be significantly reduced and eventually negligible. This allows a much stronger coupling in ultracold plasma to be realized.

Recombination fluorescence in ultracold neutral plasmas

We present the first measurements and simulations of recombination fluorescence from ultracold neutral calcium plasmas. This method probes three-body recombination at times less than 1 micros, shorter than previously published time scales. For the lowest initial electron temperatures, the recombination rate scales with the density as n0(2.2), significantly slower than the predicted n0(3). Recombination fluorescence opens a new diagnostic window in ultracold plasmas. In most cases it probes deeply bound level populations that depend critically on electron energetics. However, a perturbation in the calcium 4snd Rydberg series allows our fluorescence measurements to probe the population in weakly bound levels that result just after recombination.

Trapping and evolution dynamics of ultracold two-component plasmas

We demonstrate the trapping of a strongly magnetized, quasineutral ultracold plasma in a nested Penning trap with a background field of 2.9 T. Electrons remain trapped in this system for several milliseconds. Early in the evolution, the dynamics are driven by a breathing-mode oscillation in the ionic charge distribution, which modulates the electron trap depth. Over longer times scales, the electronic component undergoes cooling. Trap loss resulting from E x B drift is characterized.

Ultracold Neutral Plasmas

Ultracold neutral plasmas occupy an exotic regime of plasma physics in which electrons form a swarming, neutralizing background for ions that sluggishly move in a correlated manner. Strong interactions between the charged particles give rise to surprising dynamics such as oscillations of the average kinetic energy during equilibration and extremely fast recombination. Such phenomena offer stimulating and challenging problems for computational scientists, and the physics can be applied to other environments, such as the interior of gas giant planets and plasmas created by short-pulse laser irradiation of solid, liquid, and cluster targets.