One of the properties of ultracold plasmas that make them interesting objects of study is that they are cold enough that strong coupling effects can be made manifest at their typical densities. In order to study strong coupling effects, sufficiently low temperatures need to be obtained. In turn, this means that the limitations to the lowest achievable temperatures for the electrons and ions in ultracold plasmas are worth investigating as they determine the degree to which strong coupling can be achieved. In addition, understanding these limitations also illuminates the basic physics of ultracold plasmas. A DC electric field applied during ultracold plasma formation can result in significant heating of the electron component. In the work presented here, we use molecular dynamics simulations to study this heating process and determine its impact as a function of ultracold plasma parameters such as electron temperature and density. We find that this heating can have a significant impact on the lowest achievable temperatures for lower-density ultracold plasmas in particular.
One of the properties of ultracold plasmas that make them interesting objects of study is that they are cold enough that strong coupling effects can be made manifest at their typical densities. In order to study strong coupling effects, sufficiently low temperatures need to be obtained. In turn, this means that the limitations to the lowest achievable temperatures for the electrons and ions in ultracold plasmas are worth investigating as they determine the degree to which strong coupling can be achieved. In addition, understanding these limitations also illuminates the basic physics of ultracold plasmas. A DC electric field applied during ultracold plasma formation can result in significant heating of the electron component. In the work presented here, we use molecular dynamics simulations to study this heating process and determine its impact as a function of ultracold plasma parameters such as electron temperature and density. We find that this heating can have a significant impact on the lowest achievable temperatures for lower-density ultracold plasmas in particular.
Ultracold neutral plasmas provide a useful system for studying extreme parameter regime plasma physics in an accessible laboratory setting. The parameter space of plasma physics can be characterized in part by coupling strength and degree of magnetization. The range of achievable strong coupling is determined in part by the lowest possible temperatures that can be achieved. This work examines the early-lifetime electron heating of moderately coupled, strongly magnetized plasmas. This heating is dominated by disorder-induced heating (DIH) and heating due to Rydberg atom formation. By using experimentally informed simulations, it is found that DIH has a large influence in electron temperature well into the plasma lifetime. Additionally, the dependence of the minimum achievable electron temperature on magnetization and initial electron energy is examined. In this work, we find electron temperatures as low as 0.53 ± 0.05 K (for an electron density, ne, of 6.1×1012 m−3), which determines the maximum coupling strength for the measured experimental conditions.
Ultracold neutral plasmas are created by photoionizing laser-cooled atoms in a magneto-optical trap (MOT). Due to their large electrical potential energies and comparatively small kinetic energies, ultracold plasmas fall into a regime of plasma systems which are called “strongly coupled.” A priority in the field of ultracold plasmas is to generate plasmas with higher values of the strong coupling parameter Γ, which is given as the ratio of the nearest-neighbor Coulomb potential energy to the average kinetic energy. The equilibrium strong coupling in ultracold plasmas is limited by the ultrafast relaxation of the ions due to spatial disorder in the initial system. This heating mechanism is called “disorder-induced heating” (DIH) and it limits the ion strong coupling in ultracold plasmas to order unity. This thesis describes experiments that explore ways to generate higher values of the strong coupling parameter in an ultracold neutral calcium plasma.One way to increase Γ is to mitigate the effects of DIH using electron screening. This thesis describes an experiment in which the initial electron temperature was systematically changed to determine the effect that electron screening has on the ion thermalization. At lower initial electron temperatures, corresponding to a higher degree of electron shielding, it was found that the screening slows the ion thermalization and reduces the equilibrium ion temperature by as much as a factor of two. However, electron screening also reduces the ion interaction strength by the same amount, which has the net effect of leaving the effective Γ unchanged.Another method for increasing the strong coupling of an ultracold plasma is to excite the plasma ions to a higher ionization state. Simulations predict that doubly ionizing the plasma ions can increase the strong coupling in an ultracold plasma by as much as a factor of 4, with the maximum value of Γ depending on the timing of the second ionization relative to the DIH process. This thesis describes an experiment designed to test these predictions in a Ca2+ plasma. Measurements of the change in the Ca+ ion temperature as a function of the timing of the second ionization pulses were made using laser-induced fluorescence. Results of these measurements show that the heating of the Ca+ ions due to the second ionization depends on the timing of the second ionization pulses, as predicted by MD simulations.
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.
Electron heating in a magnetized ultracold neutral plasma is modeled using classical molecular dynamics simulations. A sufficiently strong magnetic field is found to reduce disorder induced heating and three body recombination heating of electrons by constraining electron motion, and therefore heating, to the single dimension aligned with the magnetic field. A strong and long-lasting temperature anisotropy develops, and the overall kinetic electron temperature is effectively reduced by a factor of three. These results suggest that experiments may increase the effective electron coupling strength using an applied magnetic field.
The expansion of ultracold neutral plasmas (UCP) is driven primarily by the thermal pressure of the electron component and is therefore sensitive to the electron temperature. For typical UCP spatial extents, evaporative cooling has a significant influence on the UCP expansion rate at lower densities (less than 108/cm3). We studied the effect of electron evaporation in this density range. Owing to the low density, the effects of three-body recombination were negligible. We modeled the expansion by taking into account the change in electron temperature owing to evaporation as well as adiabatic expansion and found good agreement with our data. We also developed a simple model for initial evaporation over a range of ultracold plasma densities, sizes, and electron temperatures to determine over what parameter range electron evaporation is expected to have a significant effect. We also report on a signal calibration technique, which relates the signal at our detector to the total number of ions and electrons in the ultracold plasma.
In this dissertation I describe a number of different computational investigations which I have undertaken during my time at Colorado State University. Perhaps the most significant of my accomplishments was the development of a general molecular dynamic model that simulates a wide variety of physical phenomena in ultracold plasmas (UCPs). This model formed the basis of most of the numerical investigations discussed in this thesis. The model utilized the massively parallel architecture of GPUs to achieve significant computing speed increases (up to 2 orders of magnitude) above traditional single core computing. This increased computing power allowed for each particle in an actual UCP experimental system to be explicitly modeled in simulations. By using this model, I was able to undertake a number of theoretical investigations into ultracold plasma systems. Chief among these was our lab's investigation of electron center-of-mass damping, in which the molecular dynamics model was an essential tool in interpreting the results of the experiment. Originally, it was assumed that this damping would solely be a function of electron-ion collisions. However, the model was able to identify an additional collisionless damping mechanism that was determined to be significant in the first iteration of our experiment. To mitigate this collisionless damping, the model was used to find a new parameter range where this mechanism was negligible. In this new parameter range, the model was an integral part in verifying the achievement of a record low measured UCP electron temperature of 1.57 ± 0.28K and a record high electron strong coupling parameter, Γ, of 0.35 ± 0.08. Additionally, the model, along with experimental measurements, was used to verify the breakdown of the standard weak coupling approximation for Coulomb collisions. The general molecular dynamics model was also used in other contexts. These included the modeling of both the formation process of ultracold plasmas and the thermalization of the electron component of an ultracold plasma. Our modeling of UCP formation is still in its infancy, and there is still much outstanding work. However, we have already discovered a previously unreported electron heating mechanism that arises from an external electric field being applied during UCP formation. Thermalization modeling showed that the ion density distribution plays a role in the thermalization of electrons in ultracold plasma, a consideration not typically included in plasma modeling. A Gaussian ion density distribution was shown to lead to a slightly faster electron thermalization rate than an equivalent uniform ion density distribution as a result of collisionless effects. Three distinct phases of UCP electron thermalization during formation were identified. Finally, the dissertation will describe additional computational investigations that preceded the general molecular dynamics model. These include simulations of ultracold plasma ion expansion driven by non-neutrality, as well as an investigation into electron evaporation. To test the effects of non-neutrality on ion expansion, a numerical model was developed that used the King model of the electron to describe the electron distribution for an arbitrary charge imbalance. The model found that increased non-neutrality of the plasma led to the rapid expansion of ions on the plasma exterior, which in turn led to a sharp ion cliff-like spatial structure. Additionally, this rapid expansion led to additional cooling of the electron component of the plasma. The evaporation modeling was used to test the underlying assumptions of previously developed analytical expression for charged particle evaporation. The model used Monte Carlo techniques to simulate the collisions and the evaporation process. The model found that neither of the underlying assumption of the charged particle evaporation expressions held true for typical ultracold plasma parameters and provides a route for computations in spite of the breakdown of these two typical assumptions.