Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

IEEE International Conference on Plasma Science (ICOPS)

Simulating The Expansion Of Ultracold Neutral Plasmas Using A Kinetic Code In Spherical Coordinates

Lucas J. Stanek, J. Haack, Michael S. Murillo

Published May 22, 2022• 0 citations
Publisher

Abstract

Ultracold neutral plasmas (UNPs) are created in a vacuum by photoionizing laser-cooled atoms. The expansion of UNPs begins as a non-equilibrium system that quickly relaxes to a local equilibrium state [1] . We aim to quantify and compare the transport processes that govern the relaxation to equilibrium in the expanding UNP. Using a 1D-1V kinetic code in spherical coordinates, we model the radial expansion of an UNP. We compare our simulation results against analytic expressions to determine where the expansion of the UNP is adiabatic or isothermal. We carry out simulations for various initial density distributions such as an exponentially decaying density distribution [2] . We ultimately aim to extend this approach to a 1D-3V multi-component kinetic code that incudes collisions [3] .

DOI
10.1109/icops45751.2022.9813239

Related papers

Entropy Generation in Ultracold Neutral Plasmas

Lucas J. Stanek, Michael S. Murillo, J. Haack

Ultracold neutral plasmas 1 (UNPs) are strongly coupled, non- equilibrium systems created by photoionizing laser cooled atoms in a vacuum. As a result, experimental platforms are developed to elucidate the properties of UNPs such as their spatial expansion and temperature relaxation 2 . Data from these experiments are used to validate theories and computational models that aim to model UNPs. One such computational model derived from a multi-species BGK kinetic theory 3 .

IEEE International Conference on Plasma Science (ICOPS)•Sep 12, 2021•0 citations

Calculation of Ion Equilibrium Temperature in Ultracold Neutral Plasmas

Jinxing Li, Mingtao Cao, Liang Han +5 more

We provide a fast iteration method to calculate the ion equilibrium temperature in an ultracold neutral plasma (UNP). The temperature as functions of electron initial temperature and ion density is obtained and compared with the recent UNP experimental data. The theoretical predictions agree with the experimental results very well. The calculated ion equilibrium temperature by this method can be applied to study the UNP expansion process more effectively.

Chinese Physics Letters•Dec 1, 2011•4 citations

Expansion of an ultracold Rydberg plasma

Gabriel Forest, Y. Li, Edwin Ward +2 more

We report a systematic experimental and numerical study of the expansion of ultracold Rydberg plasmas. Specifically, we have measured the asymptotic expansion velocities, ${v}_{0}$, of ultracold neutral plasmas (UNPs) which evolve from cold, dense samples of Rydberg rubidium atoms using ion time-of-flight spectroscopy. From this, we have obtained values for the effective initial plasma electron temperature, ${T}_{e,0}={m}_{\mathrm{ion}}{v}_{0}^{2}/{k}_{B}$ (where ${m}_{\mathrm{ion}}$ is the ${\mathrm{Rb}}^{+}$ ion mass), as a function of the original Rydberg atom density and binding energy, ${E}_{b,i}$. We have also simulated numerically the interaction of UNPs with a large reservoir of Rydberg atoms to obtain data to compare with our experimental results. We find that for Rydberg atom densities in the range ${10}^{7}--{10}^{9} {\mathrm{cm}}^{\ensuremath{-}3}$, for states with principal quantum number $n>40, {T}_{e,0}$ is insensitive to the initial ionization mechanism which seeds the plasma. In addition, the quantity ${k}_{B}\phantom{\rule{0.16em}{0ex}}{T}_{e,0}$ is strongly correlated with the fraction of atoms which ionize, and is in the range $0.6\ifmmode\times\else\texttimes\fi{}|{E}_{b,i}|\ensuremath{\lesssim}{k}_{B}{T}_{e,0}\ensuremath{\lesssim}2.5\ifmmode\times\else\texttimes\fi{}|{E}_{b,i}|$. On the other hand, plasmas from Rydberg samples with $n\ensuremath{\lesssim}40$ evolve with no significant additional ionization of the remaining atoms once a threshold number of ions has been established. The dominant interaction between the plasma electrons and the Rydberg atoms is one in which the atoms are deexcited, a heating process for electrons that competes with adiabatic cooling to establish an equilibrium where ${T}_{e,0}$ is determined by their Coulomb coupling parameter, ${\mathrm{\ensuremath{\Gamma}}}_{e}\ensuremath{\sim}0.01$.

arXiv:1702.01463•Physical Review A•Apr 2, 2018•4 citations

Expansion and electron temperature evolution in an ultracold neutral plasma

Priya Gupta

This work describes the evolution of an ultracold neutral plasma as it expands freely in vacuum. It presents a comprehensive study of the electron temperature evolution under different initial conditions. Ultracold neutral plasmas are created by photoionizing laser-cooled neutral atoms in ultrahigh vacuum. The ions are typically at a temperature of ∼ 1K while the electron temperature can be set from 1--1000 K. After photoionization, some of the highly energetic electrons escape from the cloud, leaving a net positive charge in the cloud. This creates a Coulomb well which traps the rest of the electrons, and a plasma is formed. Since the electrons have a lot of kinetic energy, they tend to leave the cloud, however, the Coulomb force from the ion pulls the electrons back into the cloud. This exerts a recoil force on the ions, and the whole plasma starts expanding radially outwards. Since the expansion is caused by the thermal pressure of the electrons, a study of the plasma expansion unravels the complicated electron temperature evolution, under different initial conditions. Many collisional processes become significant as a plasma expands. These physical processes tend to heat or cool the ions and electrons, leading to very different kinds of evolution depending on the initial conditions of the plasma. This work demonstrates three different regions of parameter space where the degree of significance of these physical processes is different during the ultracold neutral plasma evolution. The experimental results are verified by theoretical simulations, performed by Thomas Pohl, which untangle the complicated electron temperature evolution.

•Jan 1, 2007•2 citations

Expansion and Shockwave Development in Ultracold Neutral Plasmas with an Initial Exponentially Decaying Density Distribution

Mackenzie Kate Warrens

Ultracold neutral plasmas (UNPs), which are formed by photoionizing a cloud of cold atoms, have been a powerful platform for studying a wide variety of plasma phenomena. Two areas of interest include expansion into vacuum and the development of shockwaves. Previous UNP experiments have used plasmas with an initial Gaussian density distribution. UNPs with a radially exponentially decaying density distribution, or exponential UNPs, can be formed by photoionizing atoms out of a purely magnetic trap. This dissertation studies the expansion in free space of exponential UNPs as well as the development and characterization of shockwaves in exponential UNPs. The free space expansion of an exponential UNP is compared to the well-known self-similar hydrodynamic expansion of a Gaussian plasma. Similar scaling laws describe the size evolution well for both cases. The evolution of the plasma size and velocity show a characteristic ion acoustic time scale, indicating that the plasma is generally well-described by a hydrodynamic description. While the exponential UNP is predominantly in the hydrodynamic regime, heating at the central peak for low-density, high-electron temperature suggests significant local non-neutrality at early times. This thesis presents the first observation of a shockwave in a UNP. The evolution of the density and velocity for exponential UNPs show signs of wave steepening. A significant density and velocity jump over a narrow region develops. These occur at the same location in the plasma. Additionally, a large spike in the ion temperature occurs at this front. The relative ion velocity across the front modestly surpasses the local sound speed. The initial conditions were varied to characterize the development of shockwaves. Varying the density and electron temperature did not impact shock formation for the ranges used. Changing the shape caused the shocks to disappear as the density gradient decreased. This suggests that the electron thermal pressure gradient plays a significant role in the shock formation. The establishment of shock formation in UNPs opens a new avenue of research for UNPs.

•Apr 24, 2025•0 citations

High Resolution Measurement and Modeling of Ion Dynamics in an Ultracold Neutral Plasma

Patrick McQuillen

Using high-resolution laser-induced fluorescence spectroscopy, ion dynamics in\nan expanding ultracold neutral plasma (UNP) have been studied in unprecedented\ndetail. The evolution of the ion temperature is of great interest because the ions in\nan ultracold neutral plasma are strongly coupled, meaning their Coulomb interaction\nenergy exceeds the thermal energy. This leads to novel plasma properties that are\nhard to describe theoretically. Understanding all the factors that contribute to ion\ntemperature evolution is critical for designing schemes to cool the ions and achieve\neven stronger Coulomb coupling.\nThis work has also observed several phenomena that have not previously been\nstudied in ultracold plasmas. Ion adiabatic cooling has been observed in a UNP for\nthe first time, resulting in ion temperatures as low as 100mK and up to a tripling of\nthe Coulomb coupling parameter. The importance of electron-ion energy transfer is\ndemonstrated with the first observations of electron-ion collisional heating. Inclusion\nof this effect into existing numerical models provides better agreement with experimental\nresults. Results of molecular dynamics simulations of equilibrating plasmas\nhave been compared to early time experimental data, which provides a new method\nof density determination that is much less sensitive to experimental systematics and\nis accurate to within 10%. This improved certainty in the density allows accurate\ncomparison between data and a theoretical model of the plasma dynamics. The\nagreement is excellent, but on a hydrodynamic time scale a small amount of extra\nheat (on order of 100mK), not ascribable to electron-ion collisions or other sources in\nthe model, is detected in the ion component. Potential sources are discussed within\nas well as proposals for further studies and improvements.

•Apr 22, 2015•0 citations
Ultracold neutral plasmas (UNPs) are created in a vacuum by photoionizing laser-cooled atoms. The expansion of UNPs begins as a non-equilibrium system that quickly relaxes to a local equilibrium state [1] . We aim to quantify and compare the transport processes that govern the relaxation to equilibrium in the expanding UNP. Using a 1D-1V kinetic code in spherical coordinates, we model the radial expansion of an UNP. We compare our simulation results against analytic expressions to determine where the expansion of the UNP is adiabatic or isothermal. We carry out simulations for various initial density distributions such as an exponentially decaying density distribution [2] . We ultimately aim to extend this approach to a 1D-3V multi-component kinetic code that incudes collisions [3] .
We report a systematic experimental and numerical study of the expansion of ultracold Rydberg plasmas. Specifically, we have measured the asymptotic expansion velocities, v0{v}_{0}v0​, of ultracold neutral plasmas (UNPs) which evolve from cold, dense samples of Rydberg rubidium atoms using ion time-of-flight spectroscopy. From this, we have obtained values for the effective initial plasma electron temperature, Te,0=mionv02/kB{T}_{e,0}={m}_{\mathrm{ion}}{v}_{0}^{2}/{k}_{B}Te,0​=mion​v02​/kB​ (where mion{m}_{\mathrm{ion}}mion​ is the Rb+{\mathrm{Rb}}^{+}Rb+ ion mass), as a function of the original Rydberg atom density and binding energy, Eb,i{E}_{b,i}Eb,i​. We have also simulated numerically the interaction of UNPs with a large reservoir of Rydberg atoms to obtain data to compare with our experimental results. We find that for Rydberg atom densities in the range ${10}^{7}--{10}^{9} {\mathrm{cm}}^{\ensuremath{-}3}$, for states with principal quantum number n>40,Te,0n>40, {T}_{e,0}n>40,Te,0​ is insensitive to the initial ionization mechanism which seeds the plasma. In addition, the quantity kBTe,0{k}_{B}\phantom{\rule{0.16em}{0ex}}{T}_{e,0}kB​Te,0​ is strongly correlated with the fraction of atoms which ionize, and is in the range $0.6\ifmmode\times\else\texttimes\fi{}|{E}_{b,i}|\ensuremath{\lesssim}{k}_{B}{T}_{e,0}\ensuremath{\lesssim}2.5\ifmmode\times\else\texttimes\fi{}|{E}_{b,i}|$. On the other hand, plasmas from Rydberg samples with $n\ensuremath{\lesssim}40$ evolve with no significant additional ionization of the remaining atoms once a threshold number of ions has been established. The dominant interaction between the plasma electrons and the Rydberg atoms is one in which the atoms are deexcited, a heating process for electrons that competes with adiabatic cooling to establish an equilibrium where Te,0{T}_{e,0}Te,0​ is determined by their Coulomb coupling parameter, ${\mathrm{\ensuremath{\Gamma}}}_{e}\ensuremath{\sim}0.01$.
This work describes the evolution of an ultracold neutral plasma as it expands freely in vacuum. It presents a comprehensive study of the electron temperature evolution under different initial conditions. Ultracold neutral plasmas are created by photoionizing laser-cooled neutral atoms in ultrahigh vacuum. The ions are typically at a temperature of ∼ 1K while the electron temperature can be set from 1--1000 K. After photoionization, some of the highly energetic electrons escape from the cloud, leaving a net positive charge in the cloud. This creates a Coulomb well which traps the rest of the electrons, and a plasma is formed. Since the electrons have a lot of kinetic energy, they tend to leave the cloud, however, the Coulomb force from the ion pulls the electrons back into the cloud. This exerts a recoil force on the ions, and the whole plasma starts expanding radially outwards. Since the expansion is caused by the thermal pressure of the electrons, a study of the plasma expansion unravels the complicated electron temperature evolution, under different initial conditions. Many collisional processes become significant as a plasma expands. These physical processes tend to heat or cool the ions and electrons, leading to very different kinds of evolution depending on the initial conditions of the plasma. This work demonstrates three different regions of parameter space where the degree of significance of these physical processes is different during the ultracold neutral plasma evolution. The experimental results are verified by theoretical simulations, performed by Thomas Pohl, which untangle the complicated electron temperature evolution.
Ultracold neutral plasmas (UNPs), which are formed by photoionizing a cloud of cold atoms, have been a powerful platform for studying a wide variety of plasma phenomena. Two areas of interest include expansion into vacuum and the development of shockwaves. Previous UNP experiments have used plasmas with an initial Gaussian density distribution. UNPs with a radially exponentially decaying density distribution, or exponential UNPs, can be formed by photoionizing atoms out of a purely magnetic trap. This dissertation studies the expansion in free space of exponential UNPs as well as the development and characterization of shockwaves in exponential UNPs. The free space expansion of an exponential UNP is compared to the well-known self-similar hydrodynamic expansion of a Gaussian plasma. Similar scaling laws describe the size evolution well for both cases. The evolution of the plasma size and velocity show a characteristic ion acoustic time scale, indicating that the plasma is generally well-described by a hydrodynamic description. While the exponential UNP is predominantly in the hydrodynamic regime, heating at the central peak for low-density, high-electron temperature suggests significant local non-neutrality at early times. This thesis presents the first observation of a shockwave in a UNP. The evolution of the density and velocity for exponential UNPs show signs of wave steepening. A significant density and velocity jump over a narrow region develops. These occur at the same location in the plasma. Additionally, a large spike in the ion temperature occurs at this front. The relative ion velocity across the front modestly surpasses the local sound speed. The initial conditions were varied to characterize the development of shockwaves. Varying the density and electron temperature did not impact shock formation for the ranges used. Changing the shape caused the shocks to disappear as the density gradient decreased. This suggests that the electron thermal pressure gradient plays a significant role in the shock formation. The establishment of shock formation in UNPs opens a new avenue of research for UNPs.
Using high-resolution laser-induced fluorescence spectroscopy, ion dynamics in\nan expanding ultracold neutral plasma (UNP) have been studied in unprecedented\ndetail. The evolution of the ion temperature is of great interest because the ions in\nan ultracold neutral plasma are strongly coupled, meaning their Coulomb interaction\nenergy exceeds the thermal energy. This leads to novel plasma properties that are\nhard to describe theoretically. Understanding all the factors that contribute to ion\ntemperature evolution is critical for designing schemes to cool the ions and achieve\neven stronger Coulomb coupling.\nThis work has also observed several phenomena that have not previously been\nstudied in ultracold plasmas. Ion adiabatic cooling has been observed in a UNP for\nthe first time, resulting in ion temperatures as low as 100mK and up to a tripling of\nthe Coulomb coupling parameter. The importance of electron-ion energy transfer is\ndemonstrated with the first observations of electron-ion collisional heating. Inclusion\nof this effect into existing numerical models provides better agreement with experimental\nresults. Results of molecular dynamics simulations of equilibrating plasmas\nhave been compared to early time experimental data, which provides a new method\nof density determination that is much less sensitive to experimental systematics and\nis accurate to within 10%. This improved certainty in the density allows accurate\ncomparison between data and a theoretical model of the plasma dynamics. The\nagreement is excellent, but on a hydrodynamic time scale a small amount of extra\nheat (on order of 100mK), not ascribable to electron-ion collisions or other sources in\nthe model, is detected in the ion component. Potential sources are discussed within\nas well as proposals for further studies and improvements.