Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Free electron model”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A deep potential model with long-range electrostatic interactions

Machine learning models for the potential energy of multi-atomic systems, such as the deep potential (DP) model, make molecular simulations with the accuracy of quantum mechanical density functional theory possible at a cost only moderately higher than that of empirical force fields. However, the majority of these models lack explicit long-range interactions and fail to describe properties that derive from the Coulombic tail of the forces. To overcome this limitation, we extend the DP model by approximating the long-range electrostatic interaction between ions (nuclei + core electrons) and valence electrons with that of distributions of spherical Gaussian charges located at ionic and electronic sites. The latter are rigorously defined in terms of the centers of the maximally localized Wannier distributions, whose dependence on the local atomic environment is modeled accurately by a deep neural network. In the DP long-range (DPLR) model, the electrostatic energy of the Gaussian charge system is added to short-range interactions that are represented as in the standard DP model. The resulting potential energy surface is smooth and possesses analytical forces and virial. Missing effects in the standard DP scheme are recovered, improving on accuracy and predictive power. By including long-range electrostatics, DPLR correctly extrapolates to large systems the potential energy surface learned from quantum mechanical calculations on smaller systems. We illustrate the approach with three examples: the potential energy profile of the water dimer, the free energy of interaction of a water molecule with a liquid water slab, and the phonon dispersion curves of the NaCl crystal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Integrated Strategies for Overcoming Resolution Limits in Electron Beam Lithography of Chemically Amplified Resists

Electron beam lithography (EBL) of chemically amplified resists (CARs) faces fundamental challenges, including stochastic electron scattering and acid diffusion, that limit resolution and reproducibility. Using SU-8 as a model CAR, this study systematically investigated complementary strategies to address these challenges, combining multipass exposure, proximity effect correction (PEC) with midrange correction factors, base quencher incorporation, and post-exposure bake (PEB) suppression. Monte Carlo simulations and calibrated PEC modeling revealed that extending the point spread function to include a midrange scattering component significantly improved critical dimension (CD) control across varying pattern densities, correcting deviations that conventional two-term PEC failed to capture. Multipass exposure, particularly 4-pass writing with a 25% offset, redistributed the dose to average stochastic beam and scattering fluctuations, reducing line-width roughness by more than 50% and yielding more uniform nanoscale features. Photoacid confinement was investigated by adding urea as a base quencher, which successfully reduced acid diffusion but introduced substantial sensitivity penalties without improving ultimate resolution or Z-factor performance, underscoring the trade-offs of chemical versus physical confinement. Suppressing PEB most directly minimized acid diffusion, resulting in improved Z-factors and reproducible 30 nm half-pitch dense line/space patterns. Overall, these results demonstrated that PEC with midrange correction, multipass strategies, quencher additives, and PEB-free processing addresses different aspects of the EBL process window and that their integration provides a comprehensive framework for managing stochastic scattering, diffusion, and chemical amplification effects. This framework advances dense nanoscale patterning in CARs and establishes guiding principles for optimizing resist design and process strategies in high-resolution EBL and potentially other advanced lithographies, such as extreme ultraviolet (EUV) lithography.

36 MATERIALS SCIENCE↗
DOE OSTI2024

2D kinetic-ion simulations of inverted corona fusion targets

Laser-driven “inverted corona” fusion targets have attracted interest as a low-convergence neutron source and platform for studying kinetic physics. The scheme consists of a hollow or gas-filled spherical shell made of deuterated plastic. The shell has one or more laser entrance holes (LEH), resembling a spherical hohlraum. The laser passes through the LEH’s and illuminates the interior surface of the shell, ablating a plasma that travels inward towards the target center. Long ion mean free paths in the converging plasma can lead to significant interpenetration, atomic mix, and other kinetic effects. Here, in this work we report on numerical simulations of inverted corona targets using the kinetic-ion, fluid–electron hybrid particle-in-cell (PIC) approach in 2D RZ geometry. 2D simulations suggest that shape effects do not have a significant impact on plasma evolution and observed yield trends are primarily the result of 1D kinetic mix mechanisms. Simulations are also compared against available experimental data recorded at the OMEGA laser facility. In particular, synthetic x-ray emission images show good qualitative agreement with experimental results, albeit with an apparent timing discrepancy for the two-sided vacuum target. More generally, we demonstrate the potential of hybrid-PIC simulations for full-system modeling and experimental design, including collisional absorption of laser energy, plasma evolution, mix, and fusion burn.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Reanalysis of the 𝛽− $\overline{v}$ 𝑒 Angular Correlation Measurement from the aSPECT Experiment with New Constraints on Fierz Interference

On the basis of revisions of some of the systematic errors, we reanalyzed the electron-antineutrino angular correlation (𝑎 coefficient) in free neutron decay inferred from the recoil energy spectrum of the protons which are detected in 4⁢𝜋 by the aSPECT spectrometer. With 𝑎=−0.104 02⁢(82), the new value differs only marginally from the one published in 2020. Here, the experiment also has sensitivity to 𝑏, the Fierz interference term. From a correlated (𝑏,𝑎) fit to the proton recoil spectrum, we derive a limit of 𝑏 =−0.0098⁢(193) which translates into a somewhat improved 90% confidence interval region of −0.041 ≤ 𝑏 ≤ 0.022 on this hypothetical term. Tighter constraints on 𝑏 can be set from a combined analysis of the PERKEO III (𝛽 asymmetry) and aSPECT measurement which suggests a finite value of 𝑏 with 𝑏 (𝑐) =−0.0181 ± 0.0065 deviating by 2.82⁢𝜎 from the standard model.

Beta decay↗
DOE OSTI2024

Sex Differences in Odds of Brain Metastasis and Outcomes by Brain Metastasis Status after Advanced Melanoma Diagnosis

Sex differences in cancer are well-established. However, less is known about sex differences in diagnosis of brain metastasis and outcomes among patients with advanced melanoma. Using a United States nationwide electronic health record-derived de-identified database, we evaluated patients diagnosed with advanced melanoma from 1 January 2011–30 July 2022 who received an oncologist-defined rule-based first line of therapy (n = 7969, 33% female according to EHR, 35% w/documentation of brain metastases). The odds of documented brain metastasis diagnosis were calculated using multivariable logistic regression adjusted for age, practice type, diagnosis period (pre/post-2017), ECOG performance status, anatomic site of melanoma, group stage, documentation of non-brain metastases prior to first-line of treatment, and BRAF positive status. Real-world overall survival (rwOS) and progression-free survival (rwPFS) starting from first-line initiation were assessed by sex, accounting for brain metastasis diagnosis as a time-varying covariate using the Cox proportional hazards model, with the same adjustments as the logistic model, excluding group stage, while also adjusting for race, socioeconomic status, and insurance status. Adjusted analysis revealed males with advanced melanoma were 22% more likely to receive a brain metastasis diagnosis compared to females (adjusted odds ratio [aOR]: 1.22, 95% confidence interval [CI]: 1.09, 1.36). Males with brain metastases had worse rwOS (aHR: 1.15, 95% CI: 1.04, 1.28) but not worse rwPFS (adjusted hazard ratio [aHR]: 1.04, 95% CI: 0.95, 1.14) following first-line treatment initiation. Among patients with advanced melanoma who were not diagnosed with brain metastases, survival was not different by sex (rwOS aHR: 1.06 [95% CI: 0.97, 1.16], rwPFS aHR: 1.02 [95% CI: 0.94, 1.1]). This study showed that males had greater odds of brain metastasis and, among those with brain metastasis, poorer rwOS compared to females, while there were no sex differences in clinical outcomes for those with advanced melanoma without brain metastasis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Atomic-Scale Structural Mapping of Active Sites in Monolayer PGM-Free Catalysts by Low-Voltage 4D-STEM

Two-dimensional (2D) materials have attracted a large amount of attention in both basic and applied fields, and scanning transmission electron microscopy (STEM) is often uniquely well-suited for characterizing the atomic-scale structure of these materials [1-4]. As a result, STEM is poised to significantly impact progress on platinum group metal (PGM)-free catalysts, which are currently under intense development to enable low-cost, commercially viable hydrogen fuel cells [5]. While recent advancements have resulted in fuel cell performance comparable to Pt catalysts by some measures [6], cell durability remains a significant challenge, limiting practical applications [7]. Catalytically active sites in PGM-free materials are proposed to be FeN4 complexes embedded in a graphene lattice (Fig. 1b) within layered or other larger materials, but this is still under debate largely due to the range of potential actives sites predicted by computational methods and lack of methods for directly validating these models [5]. Fundamental insights into the atomic structure and resulting degradation pathways of proposed active sites are therefore needed to fully understand and control cell performance and durability [6].2D materials typically make ideal samples for STEM, but those within PGM-free catalysts present additional challenges since these materials are often defect-rich, with a high density of edges, dopant atoms, etc., which significantly increase susceptibility to beam damage at standard operating voltages. This makes analysis of potential FeN4 active sites particularly challenging, since a large proportion of Fe exists at edge sites where beam-induced atomic displacements can prohibit high-resolution structural characterization [6]. Conventional dark-field imaging compounds this problem by producing less signal for a given dose and being less sensitive to light elements than dose-efficient phase contrast imaging techniques such as those enabled by four-dimensional (4D)-STEM [8-10] (Fig. 1a). Consequently, active site structural analysis is often left to methods such as low-resolution imaging combined with quantum chemical calculations [6], which hinders accurate determination of reaction and degradation mechanisms.Here, we demonstrate direct atomic-scale structural mapping of FeN4 sites by performing low-voltage 4D-STEM on a model PGM-free catalyst system with many exposed monolayer regions. To accomplish this, we pair a 30 keV aberration-corrected probe with a fast pixelated detector that has optimal performance at low beam voltages [11]. This enables us to simultaneously image light and heavy elements with high signal-to-noise by center-of-mass analysis (Fig. 1c) while minimizing beam-induced atomic displacements at sensitive sites. The monolayer nature of these materials additionally allows for experimental validation by direct comparison with multislice simulations [12] of model structures (Fig. 1d-e). This work demonstrates how low-voltage 4D-STEM will provide new insights into the atomic-scale structure and degradation mechanisms of active sites in PGM-free catalysts, facilitating the development of low-cost hydrogen fuel cells and other energy conversion technologies in the future [13].

Zachman, Michael↗
DOE OSTI2024

Revealing the three-dimensional arrangement of polar topology in nanoparticles

In the early 2000s, low dimensional ferroelectric systems were predicted to have topologically nontrivial polar structures, such as vortices or skyrmions, depending on mechanical or electrical boundary conditions. A few variants of these structures have been experimentally observed in thin film model systems, where they are engineered by balancing electrostatic charge and elastic distortion energies. However, the measurement and classification of topological textures for general ferroelectric nanostructures have remained elusive, as it requires mapping the local polarization at the atomic scale in three dimensions. Here we unveil topological polar structures in ferroelectric BaTiO3 nanoparticles via atomic electron tomography, which enables us to reconstruct the full three-dimensional arrangement of cation atoms at an individual atom level. Our three-dimensional polarization maps reveal clear topological orderings, along with evidence of size-dependent topological transitions from a single vortex structure to multiple vortices, consistent with theoretical predictions. The discovery of the predicted topological polar ordering in nanoscale ferroelectrics, independent of epitaxial strain, widens the research perspective and offers potential for practical applications utilizing contact-free switchable toroidal moments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Poisson_CCD: A dedicated simulator for modeling CCDs

A dedicated simulator, Poisson_CCD, has been constructed which models astronomical CCDs by solving Poisson’s equation numerically and simulating charge transport within the CCD. The potentials and free carrier densities within the CCD are self-consistently solved for, giving realistic results for the charge distribution within the CCD storage wells. The simulator has been used to model the CCDs that are being used to construct the LSST digital camera. The simulator output has been validated by comparing its predictions with several different types of CCD measurements, including astrometric shifts, brighter-fatter induced pixel–pixel covariances, saturation effects, and diffusion spreading. The code is open source and freely available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The 2022 magneto-optics roadmap

Abstract Magneto-optical (MO) effects, viz. magnetically induced changes in light intensity or polarization upon reflection from or transmission through a magnetic sample, were discovered over a century and a half ago. Initially they played a crucially relevant role in unveiling the fundamentals of electromagnetism and quantum mechanics. A more broad-based relevance and wide-spread use of MO methods, however, remained quite limited until the 1960s due to a lack of suitable, reliable and easy-to-operate light sources. The advent of Laser technology and the availability of other novel light sources led to an enormous expansion of MO measurement techniques and applications that continues to this day (see section 1). The here-assembled roadmap article is intended to provide a meaningful survey over many of the most relevant recent developments, advances, and emerging research directions in a rather condensed form, so that readers can easily access a significant overview about this very dynamic research field. While light source technology and other experimental developments were crucial in the establishment of today’s magneto-optics, progress also relies on an ever-increasing theoretical understanding of MO effects from a quantum mechanical perspective (see section 2), as well as using electromagnetic theory and modelling approaches (see section 3) to enable quantitatively reliable predictions for ever more complex materials, metamaterials, and device geometries. The latest advances in established MO methodologies and especially the utilization of the MO Kerr effect (MOKE) are presented in sections 4 (MOKE spectroscopy), 5 (higher order MOKE effects), 6 (MOKE microscopy), 8 (high sensitivity MOKE), 9 (generalized MO ellipsometry), and 20 (Cotton–Mouton effect in two-dimensional materials). In addition, MO effects are now being investigated and utilized in spectral ranges, to which they originally seemed completely foreign, as those of synchrotron radiation x-rays (see section 14 on three-dimensional magnetic characterization and section 16 on light beams carrying orbital angular momentum) and, very recently, the terahertz (THz) regime (see section 18 on THz MOKE and section 19 on THz ellipsometry for electron paramagnetic resonance detection). Magneto-optics also demonstrates its strength in a unique way when combined with femtosecond laser pulses (see section 10 on ultrafast MOKE and section 15 on magneto-optics using x-ray free electron lasers), facilitating the very active field of time-resolved MO spectroscopy that enables investigations of phenomena like spin relaxation of non-equilibrium photoexcited carriers, transient modifications of ferromagnetic order, and photo-induced dynamic phase transitions, to name a few. Recent progress in nanoscience and nanotechnology, which is intimately linked to the achieved impressive ability to reliably fabricate materials and functional structures at the nanoscale, now enables the exploitation of strongly enhanced MO effects induced by light–matter interaction at the nanoscale (see section 12 on magnetoplasmonics and section 13 on MO metasurfaces). MO effects are also at the very heart of powerful magnetic characterization techniques like Brillouin light scattering and time-resolved pump-probe measurements for the study of spin waves (see section 7), their interactions with acoustic waves (see section 11), and ultra-sensitive magnetic field sensing applications based on nitrogen-vacancy centres in diamond (see section 17). Despite our best attempt to represent the field of magneto-optics accurately and do justice to all its novel developments and its diversity, the research area is so extensive and active that there remains great latitude in deciding what to include in an article of this sort, which in turn means that some areas might not be adequately represented here. However, we feel that the 20 sections that form this 2022 magneto-optics roadmap article, each written by experts in the field and addressing a specific subject on only two pages, provide an accurate snapshot of where this research field stands today. Correspondingly, it should act as a valuable reference point and guideline for emerging research directions in modern magneto-optics, as well as illustrate the directions this research field might take in the foreseeable future.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Modeling of injected current stream-induced 3D perturbations in local helicity injection plasmas

Solenoid-free tokamak startup techniques are essential for spherical tokamaks and offer a pathway to cost reduction and design simplification in fusion energy systems. Local helicity injection (LHI) is one such approach, employing compact edge current sources to drive open field line current that initiates and sustains tokamak plasmas. The recently commissioned Pegasus-III Experiment provides a platform for advancing this and other solenoid-free startup methods. This study investigates the effect of LHI on magnetic topology in Pegasus-III plasmas. A helical filament model represents the injected current, and the linear plasma response to its three-dimensional field is calculated with M3D-C1. Poincaré mapping reveals substantial flux surface degradation in all modeled cases. The onset of overlapping magnetic structures and large-scale surface deformation begins near Ψ N ≈ 0.37, indicating a broad region of perturbed topology extending toward the edge. In rotating plasmas, both single-fluid and two-fluid models exhibit partial screening of the n = 1 perturbation, with two-fluid calculations showing stronger suppression near the edge. In contrast, the absence of rotation leads to strong resonant field amplification in the single-fluid case, while the two-fluid case with zero electron rotation mitigates this amplification and preserves edge screening. Magnetic probe measurements indicate that modeling the stream with spatial spreading—representing distributed current and/or oscillatory motion—better reproduces measured magnetic power profiles than a rigid filament model. The results underscore the role of rotation and two-fluid physics in screening stream perturbations and point to plasma flow measurements and refined stream models as key steps toward improving predictive fidelity.

Schaefer, Carolyn E. [Univ. of Wisconsin, Madison,↗
DOE OSTI2023

Colossal transverse magnetoresistance due to nematic superconducting phase fluctuations in a copper oxide

Electronic anisotropy (“nematicity”) has been detected in cuprate superconductors by various experimental techniques. Using angle-resolved transverse resistance (ARTR) measurements, a very sensitive and background-free technique that can detect 0.5% anisotropy in transport, we have observed it also in La 2-x Sr x CuO 4 (LSCO) for 0.02 ≤ x ≤ 0.25. A central enigma in LSCO is the rotation of the nematic director (orientation of the largest longitudinal resistance) with temperature; this has not been seen before in any material. Here, we address this puzzle by measuring the angle-resolved transverse magnetoresistance (ARTMR) in LSCO. We report the discovery of colossal transverse magnetoresistance (CTMR)—an order-of-magnitude drop in the transverse resistivity in the magnetic field of 6 T. We show that the apparent rotation of the nematic director is caused by anisotropic superconducting fluctuations, which are not aligned with the normal electron fluid, consistent with coexisting bond-aligned and diagonal nematic orders. We quantify this by modeling the (magneto-)conductivity as a sum of normal (Drude) and paraconducting (Aslamazov–Larkin) channels but extended to contain anisotropic Drude and Cooper-pair effective mass tensors. Strikingly, the anisotropy of Cooper-pair stiffness is much larger than that of the normal electrons. It grows dramatically on the underdoped side, where the fluctuations become quasi-one-dimensional. Our analysis is general rather than model dependent. Still, we discuss some candidate microscopic models, including coupled strongly-correlated ladders where the transverse (interladder) phase stiffness is low compared with the longitudinal intraladder stiffness, as well as the anisotropic superconducting fluctuations expected close to the transition to a pair-density wave state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

High resolution electronic spectroscopy of uranium mononitride, UN

The isoelectronic molecules UN and UO + are known to have Ω = 3.5 and Ω = 4.5 ground states, respectively (where Ω is the unsigned projection of the electronic angular momentum along the internuclear axis). A ligand field theory model has been proposed to account for the difference [Matthew and Morse, J. Chem. Phys. 138, 184303 (2013)]. The ground state of UO + arises from the U 3+ (5f 3 ( 4 I 4.5 ))O 2− configuration. Owing to the higher nominal charge of the N 3− ligand, the U 3+ ion in UN is stabilized by promoting one of the 5f electrons to the more polarizable 7s orbital, reducing the repulsive interaction with the ligand and rendering U 3+ (5f 2 7s( 4 H 3.5 ))N 3− the lowest energy configuration. In the present work, we have advanced the characterization of the UN ground state through studies of two electronic transitions, [18.35]4.5-X(1)3.5 and [18.63]4.5-X(1)3.5, using sub-Doppler laser excitation techniques with fluorescence detection. Further, spectra were recorded under field-free conditions and in the presence of static electric or magnetic fields. The ground state electric dipole moment [μ = 4.30(2) D] and magnetic g e -factor [2.160(9)] were determined from these data. These values were both consistent with the 5f 2 7s configurational assignment. Dispersed fluorescence measurements were used to determine vibrational constants for the ground and first electronically excited states. Electric dipole moments and magnetic g e -factors are also reported for the higher-energy electronically excited states.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Effects of strain, disorder, and Coulomb screening on free-carrier mobility in doped cadmium oxide

Here, the interplay of stress, disorder, and Coulomb screening dictating the mobility of doped cadmium oxide (CdO) is examined using Raman spectroscopy to identify the mechanisms driving dopant incorporation and scattering within this emerging infrared optical material. Specifically, multi-wavelength Raman and UV-vis spectroscopies are combined with electrical Hall measurements on a series of yttrium (X = Y) and indium (X = In) doped X:CdO thin-films. Hall measurements confirm n-type doping and establish carrier concentrations and mobilities. Spectral fitting along the low-frequency Raman combination bands, especially the TA+TO(X) mode, reveals that the evolution of strain and disorder within the lattice as a function of dopant concentration is strongly correlated with mobility. Coupling between the electronic and lattice environments was examined through analysis of first- and second-order longitudinal–optical phonon–plasmon coupled modes that monotonically decrease in energy and asymmetrically broaden with increasing dopant concentration. By fitting these trends to an impurity-induced Fröhlich model for the Raman scattering intensity, exciton–phonon and exciton–impurity coupling factors are quantified. These coupling factors indicate a continual decrease in the amount of ionized impurity scattering with increasing dopant concentration and are not as well correlated with mobility. This shows that lattice strain and disorder are the primary determining factors for mobility in donor-doped CdO. In aggregate, the study confirms previously postulated defect equilibrium arguments for dopant incorporation in CdO while at the same time identifying paths for its further refinement.

42 ENGINEERING↗
DOE OSTI2026

Sulfite Is Not Required for N 2 Reduction Catalyzed by Mo-Nitrogenase

Mo-nitrogenase catalyzes the reduction of dinitrogen (N 2 ) to two ammonia (NH 3 ) at the active-site FeMo-cofactor. Substrate activation requires the accumulation of three or four electrons and protons as two Fe-bound hydrides and is coupled to obligatory H 2 release through reductive hydride elimination. Subsequent delivery of four or five additional electrons and protons to the bound N 2 yields two NH 3 molecules. Increasing evidence suggests that at least one belt sulfide within FeMo-cofactor is dynamically involved in the catalytic cycle. A recent report further proposed that sulfite (SO 3 2− ) is required for N 2 reduction, with sulfite binding required for NH 3 release and a subsequent six-electron reduction of the bound sulfite to regenerate the resting cofactor. To test this proposal, we conducted turnover studies of Mo-nitrogenase under sulfite-free conditions using a reduced viologen as reductant and protein preparations devoid of dithionite or sulfite. Under these conditions, nitrogenase effectively catalyzed both N 2 reduction and proton reduction, exhibiting steady-state turnover under N 2 for 6 min, with a turnover number exceeding 150, approaching that observed with dithionite as reductant. The same H 2 -formed/N 2 -reduced ratio was observed whether dithionite or the viologen species was used as reductant. Further, EPR spectroscopic analyses showed that the FeMo-cofactor returned to its resting state after multiple catalytic cycles in the absence of sulfite. Finally, physiological bypass of sulfite formation does not affect the capacity for diazotrophic growth of the model nitrogen-fixing organism Azotobacter vinelandii . Finally, these results demonstrate that sulfite is not required for Mo-nitrogenase-catalyzed N 2 reduction either in vitro or in vivo .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Optical properties of a diamond NV color center from capped embedded multiconfigurational correlated wavefunction theory

Diamond defects are among the most promising qubits. Modeling their properties through accurate quantum mechanical simulations can further their development into robust units of information. We use the recently developed capped density functional embedding theory (capped-DFET) with the multiconfigurational n-electron valence second-order perturbation theory to characterize the electronic excitation energies for different spin manifolds of the well-characterized negatively charged substitutional N defect adjacent to a vacancy (V C ) in diamond (N C V C − ). We successfully reproduce vertical excitation energies for both triplet and singlet states of N C V C − with errors < 0.1 eV. Unlike other embedding methods, capped-DFET exhibits robust predictions that are approximately independent of the embedded cluster size: it only requires a cluster to contain the defect atoms and their nearest neighbors (as small as a 40-atom capped cluster). Furthermore, our method is free from slowly converging Coulomb interactions between charged defects, and thus also only weakly dependent on supercell size.

Chemistry↗
DOE OSTI2024

Overview of recent results from the ST40 compact high-field spherical tokamak

ST40 is a compact, high-field ($B_{\mathrm{T0}}\unicode{x2A7D} 2.1\,\,\,\textrm{T}$) spherical tokamak (ST) with a mission to expand the physics and technology basis for the ST route to commercial fusion. The ST40 research programme covers confinement and stability; solenoid-free start-up; high-performance operating scenarios; and plasma exhaust. In 2022, ST40 obtained central deuterium ion temperatures of $9.6 \pm 0.4\ \textrm{keV}$, demonstrating for the first time that pilot plant relevant ion temperatures can be reached in a compact, high-field ST. Analysis of these high-ion temperature plasmas is presented, including a summary of confinement, transport and microstability characteristics, and energetic particle instabilities. Recent scenario development activities have focused on establishing diverted H-mode plasmas across a range of toroidal fields and plasma currents, along with scenarios with high non-inductive current fractions. In future operations, beginning in 2025, a 1 MW dual frequency (104/137 GHz) electron cyclotron (EC) system will be installed to enable the study of EC and electron Bernstein wave plasma start-up and current drive. Predictive modelling of the potential performance of these systems is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

In Situ Characterization of Surface Recombination in p-Si/SiO x Based Photoelectrochemical Cells

Time-resolved infrared (TRIR) and electrochemical impedance spectroscopy (EIS) were utilized to quantify surface states present at silicon oxide (SiO x ) protected crystalline p-Si electrolyte interfaces. The primary goal was to identify p-Si/SiO x photoelectrodes with both low surface recombination rates and efficient multi-electron transfer to an acceptor present in the external electrolyte. Three SiO x layers were investigated: native oxide (nOx), chemical oxide (cOx), and rapid thermal annealed (RTA) thermal oxide (tOx). Comparative study with [Ru(bpy) 3 ](PF 6 ) 2 as the electron acceptor indicated that tOx was most optimal with a small effective recombination rate, multi-electron transfer capability, and photovoltage of 500 ± 50 mV. A secondary goal was to analyze the surface recombination rates with the Shockley–Read–Hall (SRH) kinetic model. Two surface states were identified from this analysis, one closer to the CB edge (V t,1 ) and the other near the midgap (V t,2 ). EIS and SRH analyses revealed that a forming gas (5% H 2 /N 2 ) anneal (FGA) decreased surface recombination for tOx and nOx through a lower density of surface states. In the case of tOx, the infrared data indicated that V t,2 was completely removed. Here, the energetic positions of the band edges were correlated with the surface state density; low densities corresponded to more favorable potentials for inversion layer formation, which is expected to be most optimal for photocatalysis. Collectively this study indicates that the free carrier dynamics provided by TRIR represent a powerful in situ probe of the band edge and the surface state energetics in silicon based photoelectrochemical cells.

Electrochemical Impedance Spectroscopy↗
DOE OSTI2022

Kinetic Landau-fluid closures of non-Maxwellian distributions

New kinetic Landau-fluid closures, based on the cutoff Maxwellian distribution, are derived. A special static case is considered (the frequency ω=0). In the strongly collisional regime, our model reduces to Braginskii's heat flux model, and the transport is local. In the weak collisional regime, our model indicates that the heat flux is non-local and recovers the Hammett–Perkins model while the value of the cutoff velocity approaches to infinity. We compare the thermal transport coefficient χ of Maxwellian, cutoff Maxwellian and super-Gaussian distribution. The results show that the reduction of the high-speed tail particles leads to the corresponding reduction of the thermal transport coefficient χ across the entire range of collisionality, more reduction of the free streaming transport toward the weak collisional regime. In the collisionless limit, χ approaches to zero for the cutoff Maxwellian and the super-Gaussian distribution but remains finite for Maxwellian distribution. χ is complex if the cutoff Maxwellian distribution is asymmetric. The Im(χ) approaches to different convergent values in both collisionless and strongly collisional limit, respectively. It yields an additional streaming heat flux in comparison with the symmetric cutoff Maxwellian distribution. Furthermore, due to the asymmetric distribution, there is a background heat flux q0 though there is no perturbation. Furthermore, the derived Landau-fluid closures are general for fluid moment models, and applicable for the cutoff Maxwellian distribution in an open magnetic field line region, such as the scape-off-layer of Tokamak plasmas, in the thermal quench plasmas during a tokamak disruption, and the super-Gaussian electron distribution function due to inverse bremsstrahlung heating in laser-plasma studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗