Dynamic measurement of impurity ion transport in runaway electron plateaus in DIII-D
Explore the source record for details and available documents.
Hollmann, E. M.↗
SEARCH · Engineering Papers
Results for “TIME DEPENDENCE”
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.
Ultrafast infrared transient absorption spectroscopy of gas-phase Ni(CO) 4 photodissociation at 261 nm
Here, we employ ultrafast mid-infrared transient absorption spectroscopy to probe the rapid loss of carbonyl ligands from gas-phase nickel tetracarbonyl following ultraviolet photoexcitation at 261 nm. Here, nickel tetracarbonyl undergoes prompt dissociation to produce nickel tricarbonyl in a singlet excited state; this electronically excited tricarbonyl loses another CO group over tens of picoseconds. Our results also suggest the presence of a parallel, concerted dissociation mechanism to produce nickel dicarbonyl in a triplet excited state, which likely dissociates to nickel monocarbonyl. Mechanisms for the formation of these photoproducts in multiple electronic excited states are theoretically predicted with one-dimensional cuts through the potential energy surfaces and computation of spin–orbit coupling constants using equation of motion coupled cluster methods (EOM-CC) and coupled cluster theory with single and double excitations (CCSD). Bond dissociation energies are calculated with CCSD, and anharmonic frequencies of ground and excited state species are computed using density functional theory (DFT) and time-dependent density functional theory (TD-DFT).
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Laser-driven, ion-scale magnetospheres in laboratory plasmas. II. Particle-in-cell simulations
Ion-scale magnetospheres have been observed around comets, weakly magnetized asteroids, and localized regions on the Moon and provide a unique environment to study kinetic-scale plasma physics, in particular in the collision-less regime. In this work, we present the results of particle-in-cell simulations that replicate recent experiments on the large plasma device at the University of California, Los Angeles. Using high-repetition rate lasers, ion-scale magnetospheres were created to drive a plasma flow into a dipolar magnetic field embedded in a uniform background magnetic field. The simulations are employed to evolve idealized 2D configurations of the experiments, study highly resolved, volumetric datasets, and determine the magnetospheric structure, magnetopause location, and kinetic-scale structures of the plasma current distribution. We show the formation of a magnetic cavity and a magnetic compression in the magnetospheric region, and two main current structures in the dayside of the magnetic obstacle: the diamagnetic current, supported by the driver plasma flow, and the current associated with the magnetopause, supported by both the background and driver plasmas with some time-dependence. From multiple parameter scans, we show a reflection of the magnetic compression, bounded by the length of the driver plasma, and a higher separation of the main current structures for lower dipolar magnetic moments.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Momentum conservation in current drive and alpha-channeling-mediated rotation drive
Alpha channeling uses waves to extract hot ash from a fusion plasma, transferring energy from the ash to the wave. It has been proposed that this process could create a radial electric field, efficiently driving E × B rotation. However, existing theories ignore the nonresonant particles, which play a critical role in enforcing momentum conservation in quasilinear theory. Because cross field charge transport and momentum conservation are fundamentally linked, this non-consistency throws the rotation drive into question. This paper has two main goals. First, we provide a pedantic and cohesive introduction to the recently developed simple, general, self-consistent quasilinear theory for electrostatic waves that explains the torques which allow for current drive parallel to the magnetic field, and charge extraction across it; a theory that has largely resolved the question of rotation drive by alpha channeling. We show how the theory reveals a fundamental difference between the reaction of nonresonant particles to plane waves that grow in time vs steady-state waves that have a nonuniform spatial structure, allowing rotation drive in the latter case while precluding it in the former, and we review the local and global conservation laws that lead to this result. Second, we provide two new results in support of the theory. First, we provide a novel two-particle Hamiltonian model that rigorously establishes the relationship between charge transport and momentum conservation. Second, we compare the new quasilinear theory to the oscillation-center theories of ponderomotive forces, showing how the latter often obscure the time-dependent nonresonant recoil, but ultimately lead to similar results.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Deep modeling of plasma and neutral fluctuations from gas puff turbulence imaging
The role of turbulence in setting boundary plasma conditions is presently a key uncertainty in projecting to fusion energy reactors. To robustly diagnose edge turbulence, we develop and demonstrate a technique to translate brightness measurements of HeI line radiation into local plasma fluctuations via a novel integrated deep learning framework that combines neutral transport physics and collisional radiative theory for the 3 3 D - 2 3 P transition in atomic helium with unbounded correlation constraints between the electron density and temperature. The tenets for experimental validity are reviewed, illustrating that this turbulence analysis for ionized gases is transferable to both magnetized and unmagnetized environments with arbitrary geometries. Based on fast camera data on the Alcator C-Mod tokamak, we present the first two-dimensional time-dependent experimental measurements of the turbulent electron density, electron temperature, and neutral density, revealing shadowing effects in a fusion plasma using a single spectral line.
47 OTHER INSTRUMENTATION↗
Role of self-generated magnetic fields in the inertial fusion ignition threshold
Magnetic fields spontaneously grow at unstable interfaces around hot-spot asymmetries during inertial confinement fusion implosions. Although difficult to measure, theoretical considerations and numerical simulations predict field strengths exceeding 5 kT in current National Ignition Facility experiments. Magnetic confinement of electrons then reduces the rate of hot-spot heat loss by >5%. We demonstrate this via magnetic post-processing of two-dimensional xRAGE hydrodynamic simulation data at bang time. We then derive a model for the self-magnetization, finding that it varies with the square of the hot-spot temperature and inversely with the areal density. The self-magnetized Lawson analysis then gives a slightly reduced ignition threshold. Time-dependent hot-spot energy balance models corroborate this finding, with the magnetic field quadrupling the fusion yield for near-threshold parameters. The inclusion of magnetized multi-dimensional fluid instabilities could further alter the ignition threshold and will be the subject of future work.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
An optimally tuned range-separated hybrid starting point for ab initio GW plus Bethe–Salpeter equation calculations of molecules
The ab initio GW plus Bethe–Salpeter equation (GW-BSE, where G is the one particle Green's function and W is the screened Coulomb interaction) approach has emerged as a leading method for predicting excitations in both solids and molecules with a predictive power contingent upon several factors. Among these factors are the (1) generalized Kohn–Sham eigensystem used to construct the GW self-energy and to solve the BSE and (2) the efficacy and suitability of the Tamm–Dancoff approximation. Here, we present a detailed benchmark study of low-lying singlet excitations from a generalized Kohn–Sham (gKS) starting point based on an optimally tuned range-separated hybrid (OTRSH) functional. We show that the use of this gKS starting point with one-shot G0W0 and G0W0-BSE leads to the lowest mean absolute errors (MAEs) and mean signed errors (MSEs), with respect to high-accuracy reference values, demonstrated in the literature thus far for the ionization potentials of the GW100 benchmark set and for low-lying neutral excitations of Thiel’s set molecules in the gas phase, without the need for self-consistency. The MSEs and MAEs of one-shot G0W0-BSE@OTRSH excitation energies are comparable to or lower than those obtained with other functional starting points after self-consistency. Additionally, we compare these results with linear-response time-dependent density functional theory (TDDFT) calculations and find GW-BSE to be superior to TDDFT when calculations are based on the same exchange-correlation functional. This work demonstrates tuned range-separated hybrids used in combination with GW and GW-BSE can greatly suppress starting point dependence for molecules, leading to accuracy similar to that for higher-order wavefunction-based theories for molecules without the need for costlier iterations to self-consistency.
McKeon, Caroline A. (ORCID:0000000217373503)↗
Mason–Weaver theory: Revised and extended for a semi-infinite domain
Mason and Weaver developed equations to describe small particles settling under the influence of gravity and Brownian motion, including the limiting case for an infinitely deep suspension. We encountered this common convection–diffusion equation and no-flux boundary conditions in a model for dynamics of adsorbed polymers in dead end pores of a depolymerization catalyst. Close examination reveals that the Mason–Weaver solution is not correct for the infinite domain with a non-uniform initial condition. In this paper, we obtain the time dependent Green’s function for the no flux boundary condition and also for a more general reactive boundary condition. We demonstrate how the results provide solutions, via superposition, which provide solutions for several boundary conditions and all initial conditions.
Chen, Ziqiu↗
Linear stability analysis via simulated annealing and accelerated relaxation
Simulated annealing (SA) is a kind of relaxation method for finding equilibria of Hamiltonian systems. A set of evolution equations is solved with SA, which is derived from the original Hamiltonian system so that the energy of the system changes monotonically while preserving Casimir invariants inherent to noncanonical Hamiltonian systems. The energy extremum reached by SA is an equilibrium. Since SA searches for an energy extremum, it can also be used for stability analysis when initiated from a state where a perturbation is added to an equilibrium. The procedure of the stability analysis is explained, and some examples are shown. Because the time evolution is computationally time consuming, efficient relaxation is necessary for SA to be practically useful. An acceleration method is developed by introducing time dependence in the symmetric kernel used in the double bracket, which is part of the SA formulation described here. An explicit formulation for low-beta reduced magnetohydrodynamics (MHD) in cylindrical geometry is presented. In conclusion, since SA for low-beta reduced MHD has two advection fields that relax, it is important to balance the orders of magnitude of these advection fields.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Measurements of multiple heat flux components at the divertor target by using surface eroding thermocouples (invited)
The Surface Eroding Thermocouple (SETC) is a robust diagnostic utilized in DIII-D to provide fast, edge-localized modes (ELMs) resolved heat flux measurements, in particular in geometric regions that are too shadowed for traditional infrared thermography. In order to further investigate the power dissipation in the divertor region, a combination of flush-mounted and recessed SETCs was developed to assess the effect on surface heating from non-charged particles at the divertor target. Utilizing the Divertor Materials Evaluation System sample exposure platform, the first demonstration of the feasibility of using this new method to distinguish between the heat flux from charged particles and that from neutrals and radiative heating was achieved. This paper details the process of using the combination of flush SETCs and recessed SETCs to measure the multiple heat flux components at the divertor target and further discusses how to determine two important ratios, α (ratio of heat flux from charged particles deposit on recessed SETC to that deposit on flush SETC) and β (ratio of heat flux from non-charged particles deposit on recessed SETC to that deposit on flush SETC), in the estimation of the heat flux from non-charged particle sources. Using a time dependent ratio α, it was found that ~50% of the total incident heat flux is attributable to the non-charged particles in the fully detached open divertor in DIII-D. Finally, the new application of similar SETC diagnostics in the Small Angle Slot divertor with a V-like configuration and partial tungsten coated surface (SAS-VW) is also introduced.
47 OTHER INSTRUMENTATION↗
Electron temperature gradient driven transport model for tokamak plasmas
A new model for electron temperature gradient (ETG) modes is developed as a component of the multi-mode anomalous transport module [Rafiq et al., Phys Plasmas 20, 032506 (2013)] to predict a time-dependent electron temperature profile in conventional and low aspect ratio tokamaks. This model is based on two-fluid equations that govern the dynamics of low-frequency short- and long-wavelength electromagnetic toroidal ETG driven drift modes. A low collisionality NSTX discharge is used to scan the plasma parameter dependence on the ETG real frequency, growth rate, and electron thermal diffusivity. Electron thermal transport is discovered in the deep core region where modes are more electromagnetic in nature. Several previously reported gyrokinetic trends are reproduced, including the dependencies of density gradients, magnetic shear, β and gradient of β (β'), collisionality, safety factor, and toroidicity, where β is the ratio of the plasma pressure to the magnetic pressure. Furthermore, the electron heat diffusivity associated with the ETG mode is discovered to be on a scale consistent with the experimental diffusivity determined by power balance analysis.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Lagrangian particle simulation of hydrogen pellets and SPI into runaway electron beam in ITER
Numerical studies of the ablation of pellets and shattered pellet injection (SPI) fragments into a runaway electron beam in ITER have been performed using a time-dependent pellet ablation code [Samulyak et al., Nucl. Fusion, 61(4), 046007 (2021)]. The code resolves detailed ablation physics near pellet fragments and large-scale expansion of ablated clouds. The study of a single-fragment ablation quantifies the influence of various factors, in particular, the impact ionization by runaway electrons and cross-field transport models, on the dynamics of ablated plasma and its penetration into the runaway beam. Simulations of SPI performed using different numbers of pellet fragments study the formation and evolution of the ablation clouds and their large-scale dynamics in ITER. In conclusion, the penetration depth of the ablation clouds is found to be of the order of 50 cm.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Core spectroscopy of oxazole
We have measured, analyzed and simulated the ground state valence photoelectron (PES), X-ray absorption (XAS), X-ray photoelectron (XPS) as well as normal and resonant Auger-Meitner electron (AES) spectra of oxazole at the carbon, oxygen and nitrogen K-edge in order to understand its electronic structure. Experimental data is compared to theoretical calculations performed at the coupled cluster, restricted active space perturbation theory to second order and time-dependent density functional levels of theory. We demonstrate that both N and O K-edge XAS are sensitive to the amount of dynamical electron correlation included in the theoretical description, and that for a complete description of XPS, additional orbital correlation and orbital relaxation effects need to be considered. The normal AES are dominated by a singlet excitation channel and well described by theory. Furthermore, the resonant AES, however, are more complicated. While the participator decay channels, dominating at higher kinetic energies, are well described by coupled cluster theory, spectator channels can only be described satisfactorily using a method that combines restricted active space perturbation theory to second order for the bound part and a one-center approximation for the continuum.
74 ATOMIC AND MOLECULAR PHYSICS↗
Validation and parameterization of a novel physics-constrained neural dynamics model applied to turbulent fluid flow
We report, in fluid physics, data-driven models to enhance or accelerate time to solution are becoming increasingly popular for many application domains, such as alternatives to turbulence closures, system surrogates, or for new physics discovery. In the context of reduced order models of high-dimensional time-dependent fluid systems, machine learning methods grant the benefit of automated learning from data, but the burden of a model lies on its reduced-order representation of both the fluid state and physical dynamics. In this work, we build a physics-constrained, data-driven reduced order model for Navier–Stokes equations to approximate spatiotemporal fluid dynamics in the canonical case of isotropic turbulence in a triply periodic box. The model design choices mimic numerical and physical constraints by, for example, implicitly enforcing the incompressibility constraint and utilizing continuous neural ordinary differential equations for tracking the evolution of the governing differential equation. We demonstrate this technique on a three-dimensional, moderate Reynolds number turbulent fluid flow. In assessing the statistical quality and characteristics of the machine-learned model through rigorous diagnostic tests, we find that our model is capable of reconstructing the dynamics of the flow over large integral timescales, favoring accuracy at the larger length scales. More significantly, comprehensive diagnostics suggest that physically interpretable model parameters, corresponding to the representations of the fluid state and dynamics, have attributable and quantifiable impact on the quality of the model predictions and computational complexity.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Impact of photoexcitation on secondary electron emission: A Monte Carlo study
Understanding the transport of photogenerated charge carriers in semiconductors is crucial for applications in photovoltaics, optoelectronics, and photo-detectors. While recent experimental studies using scanning ultrafast electron microscopy (SUEM) have demonstrated that the local change in the secondary electron emission induced by photoexcitation enables direct visualization of the photocarrier dynamics in space and time, the origin of the corresponding image contrast still remains unclear. Here, we investigate the impact of photoexcitation on secondary electron emissions from semiconductors using a Monte Carlo simulation aided by time-dependent density functional theory. Particularly, we examine two photoinduced effects: the generation of photocarriers in the sample bulk and the surface photovoltage (SPV) effect. Using doped silicon as a model system and focusing on primary electron energies below 1 keV, we found that both the hot photocarrier effect immediately after photoexcitation and the SPV effect play dominant roles in changing the secondary electron yield (SEY), while the distribution of photocarriers in the bulk leads to a negligible change in SEY. Our work provides insights into electron–matter interaction under photo-illumination and paves the way toward a quantitative interpretation of the SUEM contrasts.
Physics↗
An efficient protocol for excited states of large biochromophores
Efficient energy transport in photosynthetic antenna is a long-standing source of inspiration for artificial light harvesting materials. However, characterizing the excited states of the constituent chromophores poses a considerable challenge to mainstream quantum chemical and semiempirical excited state methods due to their size and complexity and the accuracy required to describe small but functionally important changes in their properties. In this paper, we explore an alternative approach to calculating the excited states of large biochromophores, exemplified by a specific method for calculating the Q y transition of bacteriochlorophyll a, which we name Chl-xTB. Using a diagonally dominant approximation to the Casida equation and a bespoke parameterization scheme, Chl-xTB can match time-dependent density functional theory’s accuracy and semiempirical speed for calculating the potential energy surfaces and absorption spectra of chlorophylls. We demonstrate that Chl-xTB (and other prospective realizations of our protocol) can be integrated into multiscale models, including concurrent excitonic and point-charge embedding frameworks, enabling the analysis of biochromophore networks in a native environment. We exploit this capability to probe the low-frequency spectral densities of excitonic energies and interchromophore interactions in the light harvesting antenna protein LH2 (light harvesting complex 2). The impact of low-frequency protein motion on interchromophore coupling and exciton transport has routinely been ignored due to the prohibitive costs of including it in simulations. Our results provide a more rigorous basis for continued use of this approximation by demonstrating that exciton transition energies are unaffected by low-frequency vibrational coupling to exciton interaction energies.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Tuning nonequilibrium phase transitions with inertia
In striking contrast to equilibrium systems, inertia can profoundly alter the structure of active systems. Here, we demonstrate that driven systems can exhibit effective equilibrium-like states with increasing particle inertia, despite rigorously violating the fluctuation–dissipation theorem. Increasing inertia progressively eliminates motility-induced phase separation and restores equilibrium crystallization for active Brownian spheres. This effect appears to be general for a wide class of active systems, including those driven by deterministic time-dependent external fields, whose nonequilibrium patterns ultimately disappear with increasing inertia. The path to this effective equilibrium limit can be complex, with finite inertia sometimes acting to accentuate nonequilibrium transitions. The restoration of near equilibrium statistics can be understood through the conversion of active momentum sources to passive-like stresses. Unlike truly equilibrium systems, the effective temperature is now density dependent, the only remnant of the nonequilibrium dynamics. Further, this density-dependent temperature can in principle introduce departures from equilibrium expectations, particularly in response to strong gradients. Our results provide additional insight into the effective temperature ansatz while revealing a mechanism to tune nonequilibrium phase transitions.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Studying excited-state-specific perturbation theory on the Thiel set
Here, we explore the performance of a recently introduced N 5 -scaling excited-state-specific second order perturbation theory (ESMP2) on the singlet excitations of the Thiel benchmarking set. We find that, without regularization, ESMP2 is quite sensitive to π system size, performing well in molecules with small π systems but poorly in those with larger π systems. With regularization, ESMP2 is far less sensitive to π system size and shows a higher overall accuracy on the Thiel set than CC2, equation of motion-coupled cluster with singles and doubles, CC3, and a wide variety of time-dependent density functional approaches. Unsurprisingly, even regularized ESMP2 is less accurate than multi-reference perturbation theory on this test set, which can, in part, be explained by the set’s inclusion of some doubly excited states but none of the strong charge transfer states that often pose challenges for state-averaging. Beyond energetics, we find that the ESMP2 doubles norm offers a relatively low-cost way to test for doubly excited character without the need to define an active space.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗