Search NASA⌕ Search

SEARCH · Search NASA

Results for “ambipolar diffusion”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Gravito-turbulence and dynamo in poorly ionized protostellar discs – I. Zero-net-flux case

In their early stages, protoplanetary discs are sufficiently massive to undergo gravitational instability (GI). This instability is thought to be involved in mass accretion, planet formation via gas fragmentation, the generation of spiral density waves, and outbursts. A key and very recent area of research is the interaction between the GI and magnetic fields in young protoplanetary discs, in particular whether this instability is able to sustain a magnetic field via a dynamo. We conduct 3D, stratified shearing-box simulations using two independent codes, PLUTO and Athena++, to characterize the GI dynamo in poorly ionized protostellar discs subject to ambipolar diffusion. We find that the dynamo operates across a large range of ambipolar Elssaser number Am (which characterizes the strength of ambipolar diffusion) and is particularly strong in the regime Am = 10–100, with typical magnetic to thermal energy ratios of order unity. The dynamo is only weakly dependent on resolution (at least for Am ≲ 100), box size, and cooling law. The magnetic field is produced by the combination of differential rotation and large-scale vertical roll motions associated with spiral density waves. Our results have direct implications for the dynamo process in young protoplanetary discs and possibly some regions of active galactic nucleus discs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Formation and evolution of protostellar accretion discs – I. Angular-momentum budget, gravitational self-regulation, and numerical convergence

Here, we investigate the formation and early evolution of a protostellar disc from a magnetized prestellar core using non-ideal magnetohydrodynamic (MHD) simulations including ambipolar diffusion and Ohmic dissipation. The dynamical contraction of the prestellar core ultimately leads to the formation of a first hydrostatic core, after ambipolar diffusion decouples the magnetic field from the predominantly neutral gas. The hydrostatic core accumulates angular momentum from the infalling material, evolving into a rotationally supported torus; this ‘first hydrostatic torus’ then forms an accreting protostar and a rotationally supported disc. The disc spreads out by gravitational instability, reaching ~30 au in diameter at ~3 kyr after protostar formation. The total mass and angular momentum of the protostar–disc system are determined mainly by accretion of gas from an infalling pseudo-disc, which has low specific angular momentum because of magnetic braking; their removal from the protostar–disc system by outflow and disc magnetic braking are negligible, in part because the magnetic field is poorly coupled there. The redistribution of angular momentum within the protostar–disc system is facilitated mainly by gravitational instability; this allows formation of relatively large discs even when the specific angular momentum of infalling material is low. We argue that such discs should remain marginally unstable as they grow (with Toomre Q ~ 1–2), an idea that is broadly consistent with recent observational estimates for Class 0/I discs. We discuss the numerical convergence of our results, and show that properly treating the inner boundary condition is crucial for achieving convergence at an acceptable computational cost.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Can we observe the ion-neutral drift velocity in prestellar cores?

Given the low-ionization fraction of molecular clouds, ambipolar diffusion is thought to be an integral process in star formation. However, chemical and radiative-transfer effects, observational challenges, and the fact that the ion-neutral drift velocity is inherently very small render a definite detection of ambipolar diffusion extremely non-trivial. Here, we study the ion-neutral drift velocity in a suite of chemodynamical, non-ideal magnetohydrodynamic (MHD), two-dimensional axisymmetric simulations of prestellar cores where we alter the temperature, cosmic-ray ionization rate, visual extinction, mass-to-flux ratio, and chemical evolution. Subsequently, we perform a number of non-local thermodynamic equilibrium (non-LTE) radiative-transfer calculations considering various idealized and non-idealized scenarios in order to assess which factor (chemistry, radiative transfer, and/or observational difficulties) is the most challenging to overcome in our efforts to detect the ion-neutral drift velocity. We find that temperature has a significant effect in the amplitude of the drift velocity with the coldest modelled cores (T = 6 K) exhibiting drift velocities comparable to the sound speed. Against expectations, we find that in idealized scenarios (where two species are perfectly chemically co-evolving) the drift velocity ‘survives’ radiative-transfer effects and can in principle be observed. However, we find that observational challenges and chemical effects can significantly hinder our view of the ion-neutral drift velocity. Finally, we propose that HCN and HCNH + ⁠, being chemically co-evolving, could be used in future observational studies aiming to measure the ion-neutral drift velocity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Persistent hot carrier diffusion in boron arsenide single crystals imaged by ultrafast electron microscopy

Cubic boron arsenide (BAs) is promising for microelectronics thermal management because of its high thermal conductivity. Recently, its potential as an optoelectronic material is also being explored. However, it remains challenging to measure its photocarrier transport properties because of small sizes of available high-quality crystals. Here, we use scanning ultrafast electron microscopy (SUEM) to directly visualize the diffusion of photocarriers in BAs single crystals. Surprisingly, we observed ambipolar diffusion at low optical fluence with persistent hot carrier dynamics for above 200 ps, which can likely be attributed to the large frequency gap between acoustic and optical phonons, the same feature that is responsible for the high thermal conductivity. At higher optical fluence, we observed spontaneous electron-hole separation. Our results show BAs is an attractive optoelectronic material combining high thermal conductivity and excellent photocarrier transport properties. Finally, our study also demonstrates the capability of SUEM to probe photocarrier transport in emerging materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Production of warm ions in electron beam generated E × B plasma

Several recent experiments have demonstrated low-damage processing of 2D materials, such as graphene and single crystal diamond, using electron beam (e-beam) generated plasmas with applied crossed electric and magnetic (E × B) fields. Here, the low damage of these sensitive materials is commonly attributed to the low energy of ions incident to the substrate surface and the ion confinement in E x B fields. In this work, measurements of atom and ion velocity distribution functions in an e-beam E × B plasma at sub-mTorr argon pressures using a laser-induced fluorescence diagnostic revealed the presence of a warm population of ions with temperatures of ~ 1 eV that are sufficient to destroy the ion confinement in E × B fields and drive the ion flux by cross field diffusion in the direction opposite to the applied electric field, toward the plasma-bounded walls or substrate. Thus, it is this nearly ambipolar diffusion process that is responsible for the flux of charged particles impinging on the wall/substrate surface.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Electron kinetics in a positive column of AC discharges in a dynamic regime

Abstract We have performed hybrid kinetic-fluid simulations of a positive column in alternating current (AC) argon discharges over a range of driving frequencies f and gas pressure p for the conditions when the spatial nonlocality of the electron energy distribution function (EEDF) is substantial. Our simulations confirmed that the most efficient conditions of plasma maintenance are observed in the dynamic regime when time modulations of mean electron energy (temperature) are substantial. The minimal values of the root mean square electric field and the electron temperature have been observed at f/p values of about 3 kHz Torr −1 in a tube of radius R = 1 cm. The ionization rate and plasma density reached maximal values under these conditions. The numerical solution of a kinetic equation allowed accounting for the kinetic effects associated with spatial and temporal nonlocality of the EEDF. Using the kinetic energy of electrons as an independent variable, we solved an anisotropic tensor diffusion equation in phase space. We clarified the role of different flux components during electron diffusion in phase space over surfaces of constant total energy. We have shown that the kinetic theory uncovers a more exciting and rich physics than the classical ambipolar diffusion (Schottky) model. Non-monotonic radial distributions of excitation rates, metastable densities, and plasma density have been observed in our simulations at pR > 6 Torr cm. The predicted off-axis plasma density peak in the dynamic regime has never been observed in experiments so far. We hope our results stimulate further experimental studies of the AC positive column. The kinetic analysis could help uncover new physics even for such a well-known plasma object as a positive column in noble gases.

Physics↗
DOE OSTI2026

Charge Carrier and Spin Diffusion in a Polycrystalline and Single Crystal Lead Halide Perovskite Semiconductor

The directed movement of spin-bearing excitations in semiconductors can enable many potential schemes for spintronic applications. Understanding the mechanism of spin motion, as opposed to charge carrier or exciton motion, may require specialized techniques and sample conditions. Here we employ the noncontact and time-resolved technique of light-induced transient grating spectroscopy (LITG) to measure both carrier and spin motion in a perovskite semiconductor with varying crystallinity. The carrier motion aligns with expectations from past studies undertaken at low fluence, revealing an ambipolar diffusion coefficient on the order of 1 cm 2 /sec and diffusion length of roughly 1.5 μm for single crystals that is strongly attenuated as crystallite size decreases. The spin diffusion is measured with cross-polarized excitation beams and uncovers an intensity-dependent coefficient rising above 100 cm 2 /sec. The fast room-temperature spin relaxation limits the spin diffusion length, but the remarkable speedup of spin over carrier diffusion suggests a mechanism involving a combination of exchange-mediated and doping effects that enhance spin transport.

14 SOLAR ENERGY↗
DOE OSTI2020

Learning compact physics‐aware delayed photocurrent models using dynamic mode decomposition

Abstract Radiation‐induced photocurrent in semiconductor devices can be simulated using complex physics‐based models, which are accurate, but computationally expensive. This presents a challenge for implementing device characteristics in high‐level circuit simulations where it is computationally infeasible to evaluate detailed models for multiple individual circuit elements. In this work we demonstrate a procedure for learning compact delayed photocurrent models that are efficient enough to implement in large‐scale circuit simulations, but remain faithful to the underlying physics. Our approach utilizes dynamic mode decomposition (DMD), a system identification technique for learning reduced‐order discrete‐time dynamical systems from time series data based on singular value decomposition. To obtain physics‐aware device models, we simulate the excess carrier density induced by radiation pulses by solving numerically the ambipolar diffusion equation, then use the simulated internal state as training data for the DMD algorithm. Our results show that the significantly reduced‐order delayed photocurrent models obtained via this method accurately approximate the dynamics of the internal excess carrier density—which can be used to calculate the induced current at the device boundaries—while remaining compact enough to incorporate into larger circuit simulations.

Hanson, Joshua↗
DOE OSTI2024

Dynamics of carrier injection through V-defects in long wavelength GaN LEDs

The efficiency of high-power operation of multiple quantum well (QW) light emitting diodes (LEDs) to a large degree depends on the realization of uniform hole distribution between the QWs. In long wavelength InGaN/GaN QW LEDs, the thermionic interwell hole transport is hindered by high GaN barriers. However, in polar LED structures, these barriers may be circumvented by the lateral hole injection via semipolar {$10\bar{1}1$} QWs that form on the facets of V-defects. The efficiency of such carrier transfer depends on the transport time since transport in the semipolar QWs is competed by recombination. In this work, we study the carrier transfer from the semipolar to polar QWs by time-resolved photoluminescence in long wavelength (green to red) LEDs. We find that the carrier transfer through the semipolar QWs is fast, a few tens of picoseconds with the estimated room temperature ambipolar diffusion coefficient of ∼5.5 cm 2 /s. With diffusion much faster than recombination, the hole transport from the p-side of the structure to the polar QWs should proceed without a substantial loss, contributing to the high efficiency of long wavelength GaN LEDs.

Crystallographic defects↗
DOE OSTI2021

Striations in moderate pressure dc driven nitrogen glow discharge

Plasma stratification has been studied for more than a century. Despite the many experimental studies reported on this topic, theoretical analyses and numerical modeling of this phenomenon have been mostly limited to rare gases. In this work, a one-dimensional fluid model with detailed kinetics of electrons and vibrationally excited molecules is employed to simulate moderate-pressure (i.e. a few Torrs) dc discharge in nitrogen in a 15.5 cm long tube of radius 0.55 cm. The model also considers ambipolar diffusion to account for the radial loss of ions and electrons to the wall. The proposed model predicts self-excited standing striations in nitrogen for a range of discharge currents. The impact of electron transport parameters and reaction rates obtained from a solution of local two-term and a multi-term Boltzmann equation on the predictions are assessed. In-depth kinetic analysis indicates that the striations result from the undulations in electron temperature caused due to the interaction between ionization and vibrational reactions. Furthermore, the vibrationally excited molecules associated with the lower energy levels are found to influence nitrogen plasma stratification and the striation pattern strongly. A balance between ionization processes and electron energy transport allows the formation of the observed standing striations. Simulations were conducted for a range of discharge current densities from ~0.018 to 0.080 mA cm –2 , for an operating pressure of 0.7 Torr. Parametric studies show that the striation length decreases with increasing discharge current. Furthermore, the predictions from the model are compared against experimental measurements and are found to agree favorably.

42 ENGINEERING↗
DOE OSTI2023

Particle-in-cell modeling of electron beam generated plasma

Plasmas generated using energetic electron beams are well known for their low electron temperature (T e ) and plasma potential, which makes them attractive for atomic-precision plasma processing applications such as atomic layer etch and deposition. A 2-dimensional particle-in-cell model for an electron beam-generated plasma in argon confined by a constant applied magnetic field is described here in this article. Plasma production primarily occurs in the path of the beam electrons in the center of the chamber. The resulting plasma spreads out in the chamber through non-ambipolar diffusion with a short-circuiting effect allowing unequal electron and ion fluxes to different regions of the bounding conductive chamber walls. The cross-field transport of the electrons (and thus the steady-state characteristics of the plasma) are strongly impacted by the magnetic field. T e is anisotropic in the electron beam region, but low and isotropic away from the plasma production zone. The plasma density increases and the plasma becomes more confined near the region of production when the magnetic field strengthens. The magnetic field reduces both electron physical and energy transport perpendicular to the magnetic field. T e is uniform along the magnetic field lines and slowly decreases perpendicular to it. Electrons are less energetic in the sheath regions where the sheath electric field repels and confines the low-energy electrons from the bulk plasma. Even though electron and ion densities are similar in the bulk plasma due to quasi-neutrality, electron and ion fluxes on the grounded chamber walls are unequal at most locations. Electron confinement by the magnetic field weakens with increasing pressure, and the plasma spread out farther from the electron beam region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Magnetorotational instability in a swirling partially ionized gas

ABSTRACT The magnetorotational instability (MRI) has been proposed as the method of angular momentum transport that enables accretion in astrophysical discs. However, for weakly ionized discs, such as protoplanetary discs, it remains unclear whether the combined non-ideal magnetohydrodynamic (MHD) effects of Ohmic resistivity, ambipolar diffusion, and the Hall effect make these discs MRI stable. While much effort has been made to simulate non-ideal MHD MRI, these simulations make simplifying assumptions and are not always in agreement with each other. Furthermore, it is difficult to directly observe the MRI astrophysically because it occurs on small scales. Here, we propose the concept of a swirling gas experiment of weakly ionized argon gas between two concentric cylinders threaded with an axial magnetic field that can be used to study non-ideal MHD MRI. For our proposed experiment, we derive the hydrodynamic equilibrium flow and a dispersion relation for MRI that includes the three non-ideal effects. We solve this dispersion relation numerically for the parameters of our proposed experiment. We find it should be possible to produce a non-ideal MRI in such an experiment because of the Hall effect, which increases the MRI growth rate when the vertical magnetic field is anti-aligned with the rotation axis. As a proof of concept, we also present experimental results for a hydrodynamic flow in an unmagnetized prototype. We find that our prototype has a small, but non-negligible, α-parameter that could serve as a baseline for comparison to our proposed magnetized experiment, which could be subject to additional turbulence from the MRI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Revealing the Drag Instability in One-fluid Nonideal Magnetohydrodynamic Simulations of a 1D Isothermal C-shock

C-type shocks are believed to be ubiquitous in turbulent molecular clouds thanks to ambipolar diffusion. We investigate whether the drag instability in 1D isothermal C-shocks, inferred from the local linear theory of Gu & Chen, can appear in nonideal magnetohydrodynamic simulations. Two C-shock models (with narrow and broad steady-state shock widths) are considered to represent the typical environment of star-forming clouds. The ionization-recombination equilibrium is adopted for the one-fluid approach. In the 1D simulation, the inflow gas is continuously perturbed by a sinusoidal density fluctuation with a constant frequency. The perturbations clearly grow after entering the C-shock region until they start being damped at the transition to the post-shock region. We show that the profiles of a predominant Fourier mode extracted locally from the simulated growing perturbation match those of the growing mode derived from the linear analysis. Moreover, the local growth rate and wave frequency derived from the predominant mode generally agree with those from the linear theory. Therefore, we confirm the presence of the drag instability in simulated 1D isothermal C-shocks. We also explore the nonlinear behavior of the instability by imposing larger-amplitude perturbations to the simulation. We find that the drag instability is subject to wave steepening, leading to saturated perturbation growth. Issues concerning local analysis, nonlinear effects, one-fluid approach, and astrophysical applications are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Charge decay in the spatial afterglow of plasmas and its impact on diffusion regimes

Abstract The spatial afterglow is a region at the boundary of a non-equilibrium plasma where charged species relax into ambient equilibrium. In many applications, the spatial afterglow is the part of the plasma that interacts with surfaces, such as suspended particles or a material substrate. However, compared to the bulk plasma, there has been little effort devoted to studying the properties of the spatial afterglow, and a fundamental analysis has not yet been developed. Here, we apply double Langmuir probe measurements and develop an advection-diffusion-recombination model to provide a detailed description of charged species in the spatial afterglow over a wide range of pressures, temperatures, plasma dimensions, and flow rates. We find that the density of charged species in the spatial afterglow decays by orders of magnitude, which leads to a transition from ambipolar to free diffusion. These insights can be used to explain or predict experimental observations of phenomena, such as the charging of dust grains and the dose of charged species to a biomaterial.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Particle charge distributions in the effluent of a flow-through atmospheric pressure low temperature plasma

Atmospheric pressure low-temperature plasmas are often utilized to perform particle synthesis, treatment, and removal. It is well-known that dust particles are highly negatively charged in these plasmas; however, little is known about dust particle charging behavior as particles leave the plasma volume and pass through the spatial afterglow region. In this work, monodisperse particles of various sizes and work functions were introduced into an atmospheric pressure radiofrequency capacitively coupled flow-through plasma. Dust particle electrical mobility distributions downstream of the flow-through plasma were measured utilizing a differential mobility analyzer in conjunction with a condensation particle counter at various gas flow velocities. Charge distributions were determined from the measured electrical mobility distributions. Experiments confirm that particles become less negatively charged, and even net-positively charged after leaving the plasma volume, with a distribution that follows a shifted Boltzmann charge distribution. Additionally, particle charge in the effluent of the flow-through plasma is negligibly dependent on work function but highly size and flow velocity dependent. Larger particles were shown to have a higher magnitude of charge under all studied conditions; however, particle polarity was switchable by varying gas flow velocity. The charging dynamics were simulated utilizing a constant number Monte Carlo model that accounts for electron temperature decay and the transition from ambipolar to free diffusion of electrons and ions in the spatial afterglow. Simulation results also suggest that, at the same flow velocity, larger particles obtain a greater magnitude of charge, negative or positive. Furthermore, the decrease in electron mobility and the difference between ion and electron convective loss rates create an ion-rich region in the plasma effluent that promotes ion–particle collisions and drives particle charge removal and even reversal of polarity. Larger particles more favorably collide with energetic species in these environments, which results in higher charge states.

atmospheric pressure low temperature plasma↗
DOE OSTI2020

Long-range exciton diffusion in molecular non-fullerene acceptors

Abstract The short exciton diffusion length associated with most classical organic semiconductors used in organic photovoltaics (5-20 nm) imposes severe limits on the maximum size of the donor and acceptor domains within the photoactive layer of the cell. Identifying materials that are able to transport excitons over longer distances can help advancing our understanding and lead to solar cells with higher efficiency. Here, we measure the exciton diffusion length in a wide range of nonfullerene acceptor molecules using two different experimental techniques based on photocurrent and ultrafast spectroscopy measurements. The acceptors exhibit balanced ambipolar charge transport and surprisingly long exciton diffusion lengths in the range of 20 to 47 nm. With the aid of quantum-chemical calculations, we are able to rationalize the exciton dynamics and draw basic chemical design rules, particularly on the importance of the end-group substituent on the crystal packing of nonfullerene acceptors.

14 SOLAR ENERGY↗
DOE OSTI2023

Toroidal modeling of plasma flow damping and density pump-out by RMP during ELM mitigation in HL-2A

Reduction of both the plasma density and toroidal flow speed, due to application of the predominantly n = 1 (n is the toroidal mode number) resonant magnetic perturbation (RMP) for controlling the edge localized mode in the HL-2A tokamak, is numerically investigated utilizing the quasi-linear initial-value code MARS-Q (Liu et al 2013 Phys. Plasmas 20 042503). Simulation results reveal that the neoclassical toroidal viscosity (NTV) due to three dimensional fields plays the key role in modifying the plasma momentum and particle transport in the HL-2A discharge. By comparing the modeling results with the measured density pump-out in the experiment, the electron NTV particle flux model, in combination with the free-boundary condition for the axisymmetric change of the density at the plasma edge, is found to yield the best agreement in terms of both the pump-out level and the overall time scale. Further sensitivity studies show that the simulated density pump-out level is reasonably robust against variations in the model assumptions, including the particle diffusion model and the non-ambipolar versus ambipolar NTV particle flux. The latter however affects the time scale for reaching the steady state solution. Finally, it is found that the plasma edge-peeling response, the NTV torque, as well as the plasma momentum and particle transport, all are sensitive to the toroidal phase difference between the upper and lower rows of the RMP coil currents in HL-2A, with the 30 ∘ coil phasing producing the minimal side effects on the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Electron transport in a tokamak scrape-off layer: impact of toroidal nonuniformities of divertor targets

An analysis of electron dynamics is carried out in the tokamak scrape-off layer (SOL) region. Small non-axisymmetric perturbations of the divertor target sheath potential affect electron drift orbits in the SOL and may lead to non-ambipolar electron radial transport. The resulting rates of electron convection, diffusion, and mobility are calculated analytically, and the analytic results are compared with direct numerical simulations of electron drift orbits. The proposed mechanism of non-ambipolar electron transport may be relevant to the sustainment of the SOL plasma quasi-neutrality in the ‘heuristic’ model of the SOL width (Eich et al 2011 Phys. Rev. Lett. 107 215001).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗