Researcher profile

Scott D. Baalrud

Scott D. Baalrud contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
14works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2020arXiv

Quantum Mean-Force Kinetic Theory: General Formulation and Application to Electron-Ion Transport in Warm Dense Matter

We present an approach to extend plasma transport theory into the Warm Dense Matter (WDM) regime characterized by moderate Coulomb coupling and electron degeneracy. It is based on a recently proposed closure of the BBGKY hierarchy that expands in terms of the departure of correlations from their equilibrium value, rather than in terms of the strength of correlations. This kinetic equation contains modifications to the collision term in addition to a second term that models the non-ideal contributions to the equation of state. An explicit collision operator is derived in the semiclassical limit that is similar to that of the Uehling-Uhlenbeck equation, but where scattering is mediated by the potential of mean force (PMF). As a demonstration, we use this collision integral to evaluate temperature and momentum relaxation rates in dense plasmas. We obtain degeneracy- and coupling-dependent 'Coulomb integrals' that take the place of $\lnΛ$ in the scattering rates. We additionally find a novel difference in the way in which degeneracy influences momentum relaxation in comparison to temperature relaxation. Finally, we evaluate electron-ion relaxation rates for the case of warm dense deuterium over a range of density and temperature spanning the classical to quantum and weak to strong coupling transitions. Results are compared with the Landau-Spitzer rate and rates obtained from the quantum Landau-Fokker-Planck equation and Lee-More model. We find that the models diverge significantly in the degenerate and moderately coupled regime and attribute this difference to how the various models treat the physics of Pauli blocking, correlations, large-angle scattering, and diffraction.

preprint2019arXiv

Interaction of Biased Electrodes and Plasmas: Sheaths, Double Layers and Fireballs

Biased electrodes are common components of plasma sources and diagnostics. The plasma-electrode interaction is mediated by an intervening sheath structure that influences properties of the electrons and ions contacting the electrode surface, as well as how the electrode influences properties of the bulk plasma. A rich variety of sheath structures have been observed, including ion sheaths, electron sheaths, double sheaths, double layers, anode glow, and fireballs. These represent complex self-organized responses of the plasma that depend not only on the local influence of the electrode, but also on the global properties of the plasma and the other boundaries that it is in contact with. This review summarizes recent advances in understanding the conditions under which each type of sheath forms, what the basic stability criteria and steady-state properties of each are, and the ways in which each can influence plasma-boundary interactions and bulk plasma properties. These results may be of interest to a number of application areas where biased electrodes are used, including diagnostics, plasma modification of materials, plasma sources, electric propulsion, and the interaction of plasmas with objects in space.

preprint2016arXiv

Electron Preheaths: The Outsized Influence of Positive Boundaries on Plasmas

Electron sheaths form near the surface of objects biased more positive than the plasma potential, such as in the electron saturation region of a Langmuir probe trace. Generally, the formation of electron sheaths requires that the electron-collecting area be sufficiently smaller ($\sqrt{2.3m_e/M}$ times) than the ion-collecting area. They are commonly thought to be local phenomena that collect the random thermal electron current, but do not otherwise perturb a plasma. Here, using experiments on an electrode embedded in a wall, particle-in-cell simulations and theory, it is shown that under low temperature plasma conditions ($T_e \gg T_i$) electron sheaths are far from local. Instead, a long presheath region (in excess of 70 $λ_D$) extends into the plasma where electrons are accelerated via a pressure gradient to a flow speed exceeding the electron thermal speed at the sheath edge. This fast flow is found to excite instabilities, causing strong fluctuations near the sheath edge.

preprint2016arXiv

Particle-in-cell study of the ion-to-electron sheath transition

The form of a sheath near a small electrode, with bias changing from below to above the plasma potential is studied using 2D particle-in-cell (PIC) simulations. Five cases are studied: (A) an electrode biased more than the electron temperature ($T_e/e$) below the plasma potential, (B) an electrode biased less than $T_e/2e$ below the plasma potential, (C) an electrode biased nearly at the plasma potential, (D) an electrode biased more than $T_i/2e$ but less than $T_e/2e$ above the plasma potential, and (E) an electrode biased much greater than $T_e/2e$ above the plasma potential. In case (A), the electron velocity distribution function (EVDF) is observed to be Maxwellian with a Boltzmann-type exponential density decay through the ion sheath and presheath. In cases (B) and (C), the EVDFs exhibit a loss-cone type truncation due to fast electrons overcoming the small potential difference between the electrode and plasma. No sheath is present in this regime, and the plasma remains quasineutral up to the electrode. The EVDF truncation leads to a presheath-like density and flow velocity gradient. In case (D) an electron sheath is present, and essentially all ions are repelled. Here the truncation driven behavior persists, but is accompanied by a shift in the maximum value of the EVDF that is not present in the negative bias cases. In case (E), the flow shift becomes greater and the loss-cone moves further into the tail of the EVDF. In this case the flow moment has significant contributions from both the flow shift of the EVDF maximum, and the loss-cone truncation.

preprint2015arXiv

Influence of Ion Streaming Instabilities on Transport Near Plasma Boundaries

Plasma boundary layers are susceptible to electrostatic instabilities driven by ion flows in presheaths and, when present, these instabilities can influence transport. In plasmas with a single species of positive ion, ion-acoustic instabilities are expected under conditions of low pressure and large electron-to-ion temperature ratio ($T_e/T_i \gg 1$). In plasmas with two species of positive ions, ion-ion two-stream instabilities can also be excited. The stability phase-space is characterized using the Penrose criterion and approximate linear dispersion relations. Predictions for how these instabilities affect ion and electron transport in presheaths, including rapid thermalization due to instability-enhanced collisions and an instability-enhanced ion-ion friction force, are also briefly reviewed. Recent experimental tests of these predictions are discussed along with research needs required for further validation. The calculated stability boundaries provide a guide to determine the experimental conditions at which these effects can be expected.

preprint2015arXiv

Modified Enskog Kinetic Theory for Strongly Coupled Plasmas

Concepts underlying the Enskog kinetic theory of hard-spheres are applied to include short-range correlation effects in a model for transport coefficients of strongly coupled plasmas. The approach is based on an extension of the effective potential transport theory [S.~D.~Baalrud and J.~Daligault, Phys.~Rev.~Lett.~{\bf 110}, 235001 (2013)] to include an exclusion radius surrounding individual charged particles that is associated with Coulomb repulsion. This is obtained by analogy with the finite size of hard spheres in Enskog's theory. Predictions for the self-diffusion and shear viscosity coefficients of the one-component plasma are tested against molecular dynamics simulations. The theory is found to accurately capture the kinetic contributions to the transport coefficients, but not the potential contributions that arise at very strong coupling ($Γ\gtrsim 30$). Considerations related to a first-principles generalization of Enskog's kinetic equation to continuous potentials are also discussed.

preprint2015arXiv

Theory of the Electron Sheath and Presheath

Electron sheaths are commonly found near Langmuir probes collecting the electron saturation current. The common assumption is that the probe collects the random flux of electrons incident on the sheath, which tacitly implies that there is no electron presheath and that the flux collected is due to a velocity space truncation of the electron velocity distribution function (EVDF). This work provides a dedicated theory of electron sheaths, which suggests that they are not so simple. Motivated by EVDFs observed in Particle-In-Cell (PIC) simulations, a 1D model for the electron sheath and presheath is developed. In the model, under low temperature plasma conditions ($T_e\gg T_i$), an electron pressure gradient accelerates electrons in the presheath to a flow velocity that exceeds the electron thermal speed at the sheath edge. This pressure gradient generates large flow velocities compared to what would be generated by ballistic motion in response to the electric field. It is found that in many situations, under common plasma conditions, the electron presheath extends much further into the plasma than an analogous ion presheath. PIC simulations reveal that the ion density in the electron presheath is determined by a flow around the electron sheath and that this flow is due to 2D aspects of the sheath geometry. Simulations also indicate the presence of ion acoustic waves excited by the differential flow between electrons and ions in the presheath which result in sheath edge fluctuations. The 1D model and time averaged PIC simulations are compared and it is shown that the model provides a good description of the electron sheath and presheath.

preprint2014arXiv

Accurate Determination of the Shear Viscosity of the One-Component Plasma

The shear viscosity coefficient of the one-component plasma is calculated with unprecedented accuracy using equilibrium molecular dynamics simulations and the Green-Kubo relation. Numerical and statistical uncertainties and their mitigation for improving accuracy are analyzed. In the weakly coupled regime, our the results agree with the Landau-Spitzer prediction. In the moderately and strongly coupled regimes, our results are found in good agreement with recent results obtained for the Yukawa one-component plasma using non-equilibrium molecular dynamics. A practical formula is provided for evaluating the viscosity coefficient across coupling regimes, from the weakly-coupled regime up to solidification threshold. The results are used to test theoretical predictions of the viscosity coefficients found in the literature.

preprint2014arXiv

Effective Potential Theory: A Practical Way to Extend Plasma Transport Theory to Strong Coupling

The effective potential theory is a physically motivated method for extending traditional plasma transport theories to stronger coupling. It is practical in the sense that it is easily incorporated within the framework of the Chapman-Enskog or Grad methods that are commonly applied in plasma physics and it is computationally efficient to evaluate. The extension is to treat binary scatterers as interacting through the potential of mean force, rather than the bare Coulomb or Debye-screened Coulomb potential. This allows for aspects of many-body correlations to be included in the transport coefficients. Recent work has shown that this method accurately extends plasma theory to orders of magnitude stronger coupling when applied to the classical one-component plasma model. The present work shows that similar accuracy is realized for the Yukawa one-component plasma model and it provides a comparison with other approaches.

preprint2014arXiv

Extensions and Applications of the Bohm Criterion

The generalized Bohm criterion is revisited in the context of incorporating kinetic effects of the electron and ion distribution functions into the theory. The underlying assumptions and results of two different approaches are compared: The conventional `kinetic Bohm criterion' and a fluid-moment hierarchy approach. The former is based on the asymptotic limit of an infinitely thin sheath ($λ_D/l =0$), whereas the latter is based on a perturbative expansion of a sheath that is thin compared to the plasma ($λ_D/l \ll 1$). Here $λ_D$ is the Debye length, which characterizes the sheath length scale, and $l$ is a measure of the plasma or presheath length scale. The consequences of these assumptions are discussed in terms of how they restrict the class of distribution functions to which the resulting criteria can be applied. Two examples are considered to provide concrete comparisons between the two approaches. The first is a Tonks-Langmuir model including a warm ion source [Robertson 2009 \textit{Phys.\ Plasmas} {\bf 16} 103503]. This highlights a substantial difference between the conventional kinetic theory, which predicts slow ions dominate at the sheath edge, and the fluid moment approach, which predicts slow ions have little influence. The second example considers planar electrostatic probes biased near the plasma potential using model equations and particle-in-cell simulations. This demonstrates a situation where electron kinetic effects alter the Bohm criterion, leading to a subsonic ion flow at the sheath edge.

preprint2014arXiv

Instability-Enhanced Friction in the Presheath of Two-Ion-Species Plasmas

The speed at which ions enter a sheath is a fundamental property of a plasma that also provides a useful boundary condition in modeling. A recent theory proposed that this can be significantly influenced by an instability-enhanced friction force arising from two-stream instabilities in the presheath when multiple ion species are present. Although experiments appeared to confirm this theory, recent particle simulations have brought it into question. We reconcile this controversy using direct numerical solutions of the dispersion relation, which show that there is a dependence on the electron-ion temperature ratio that was not considered previously. In addition, particle-in-cell simulations are used to show that ion-ion two-stream instabilities can arise near the sheath edge and generate an enhanced ion-ion friction force. Only by accounting for the instability-enhanced friction force can theory predict the simulated ion speeds at the sheath edge.

preprint2013arXiv

An Effective Potential Theory for Transport Coefficients Across Coupling Regimes

A plasma transport theory that spans weak to strong coupling is developed from a binary collision picture, but where the interaction potential is taken to be an effective potential that includes correlation effects and screening self-consistently. This physically motivated approach provides a practical model for evaluating transport coefficients across coupling regimes. The theory is shown to compare well with classical molecular dynamics simulations of temperature relaxation in electron-ion plasmas, as well as simulations and experiments of self-diffusion in one component plasmas. The approach is versatile and can be applied to other transport coefficients as well.

preprint2013arXiv

The incomplete plasma dispersion function: properties and application to waves in bounded plasmas

The incomplete plasma dispersion function is a generalization of the plasma dispersion function in which the defining integral spans a semi-infinite, rather than infinite, domain. It is useful for describing the linear dielectric response and wave dispersion in non-Maxwellian plasmas when the distribution functions can be approximated as Maxwellian over finite, or semi-infinite, intervals in velocity phase-space. A ubiquitous example is the depleted Maxwellian electron distribution found near boundary sheaths or double layers, where the passing interval can be modeled as Maxwellian with a lower temperature than the trapped interval. The depleted Maxwellian is used as an example to demonstrate the utility of using the incomplete plasma dispersion function for calculating modifications to wave dispersion relations.

preprint2011arXiv

A generalized Coulomb logarithm for strongly coupled plasmas

Transport coefficients for Coulomb collision processes, e.g., friction, energy exchange, and resistivity, are calculated for Debye-screened (Yukawa) plasmas including correlation effects within the binary collision approximation. A generalized Coulomb logarithm, $Ξ(Λ)$, is required, which asymptotes to $\ln Λ$ in the weakly coupled limit, but is significantly modified when $Λ\lesssim 10^2$ due to large angle collisions. In the regime $10^{-2} \lesssim Λ\lesssim 10^2$, $Ξ$ can be substantially larger for attractive collisions than repulsive.