SEARCH · Engineering Papers

Results for “orbit design”

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.

SourceResearch connectionFrom yearThrough year

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

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Advances and Challenges of SCAN and r 2 SCAN Density Functionals in Transition‐Metal Compounds

ABSTRACT Transition‐metal compounds (TMCs) with open‐shell d ‐electrons are characterized by a complex interplay of lattice, charge, orbital, and spin degrees of freedom, giving rise to various fascinating applications. Often exhibiting exotic properties, these compounds are commonly classified as correlated systems due to strong inter‐electronic interactions called Hubbard U . This inherent complexity presents significant challenges to Kohn‐Sham density functional theory (KS‐DFT), the most widely used electronic structure method in condensed matter physics and materials science. While KS‐DFT is, in principle, exact for the ground‐state total energy, its exchange‐correlation energy must be approximated in practice. The mean‐field nature of KS implementations, combined with the limitations of current exchange‐correlation density functional approximations, has led to the perception that DFT is inadequate for correlated systems, particularly TMCs. Consequently, a common workaround involves augmenting DFT with an on‐site Hubbard‐like U correction. In recent years, the strongly constrained and appropriately normed (SCAN) density functional, along with its refined variant r 2 SCAN, has achieved remarkable progress in accurately describing the structural, energetic, electronic, magnetic, and vibrational properties of TMCs, challenging the traditional perception of DFT's limitations. This review explores the design principles of SCAN and r 2 SCAN, highlights their key advancements in studying TMCs, explains the mechanisms driving these improvements, and addresses the remaining challenges in this evolving field.

Zhang, Yubo [Minjiang Collaborative Center for The↗

DOE OSTI2020

Electron Dropout Events and Flux Enhancements Associated With Geomagnetic Storms Observed by PROBA-V/Energetic Particle Telescope From 2013 to 2019

Electron flux variations for E > 500 keV during geomagnetic storms are investigated using the Energetic Particle Telescope (EPT). This detector launched in May 2013 on board the satellite PROBA-V at an altitude of 820 km was designed to provide uncontaminated spectra of electrons, protons, and alpha particles. Electron flux dropout events are observed during the main phase of each storm and even during substorms: a rapid reduction of the electron flux is noted throughout the outer electron radiation belt at all energies above about 0.5 MeV on timescales of a few hours. The electron spectrograms measured by the EPT between 2013 and 2019 show that after each geomagnetic storm, dropout events are followed by a flux enhancement starting first at low L values, and reaching the slot or even the inner belt for the strongest storms. We determine the link between Disturbed Storm Time (Dst) and the minimum value of the L-shell where the dropouts deplete the outer belt, as well as the nonlinear relation between Dst and the minimum L-shell where the flux penetrates in the slot region or even the inner belt during the storms. Dropouts appear at all energies measured by EPT and penetrate down to $\textit{L}$~3.5 for the strongest events. Dropouts are observed at Low Earth Orbit each time Dst has an inverted peak < –40 nT. Flux enhancements appear at lower L only for big storm events with Dst < –50 nT. They penetrate down to an impenetrable barrier with a minimum L-shell related to Dst and to the energy. For E > 1 MeV, this limit is also linked to the plasmapause position.

58 GEOSCIENCES↗

DOE OSTI2021

Strong Ferromagnetic Exchange Coupling and Single-Molecule Magnetism in MoS 4 3– -Bridged Dilanthanide Complexes

We report the synthesis and characterization of the trinuclear 4d-4f compounds [Co(C 5 Me 5 ) 2 ][(C 5 Me 5 ) 2 Ln(μ-S) 2 Mo(μ-S) 2 Ln(C 5 Me 5 ) 2 ], 1-Ln (Ln = Y, Gd, Tb, Dy), containing the highly polarizable MoS 4 3 - bridging unit. UV-Vis-NIR diffuse reflectance spectra and DFT calculations of 1-Ln reveal a low-energy metal-to-metal charge transfer transition assigned to charge transfer from the singly occupied 4d z 2 orbital of Mo V to the empty 5d orbitals of the lanthanides (4d in the case of 1-Y ), mediated by sulfur-based 3p orbitals. Electron paramagnetic resonance spectra collected for 1-Y in a tetrahydrofuran solution show large 89 Y hyperfine coupling constants of A ⊥ = 23 MHz and A || = 26 MHz, indicating the presence of significant yttrium-localized unpaired electron density. Magnetic susceptibility data support similar electron delocalization and ferromagnetic Ln-Mo exchange for 1-Gd , 1-Tb , and 1-Dy . This ferromagnetic exchange gives rise to an S = 15/2 ground state for 1-Gd and one of the largest magnetic exchange constants involving Gd III observed to date, with J Gd-Mo = +16.1(2) cm -1 . Additional characterization of 1-Tb and 1-Dy by ac magnetic susceptibility measurements reveals that both compounds exhibit slow magnetic relaxation. Although a Raman magnetic relaxation process is dominant for both 1-Tb and 1-Dy , an extracted thermal relaxation barrier of U eff = 68 cm -1 for 1-Dy is the largest yet reported for a complex containing a paramagnetic 4d metal center. Furthermore, these results provide a potentially generalizable route to enhanced n d-4f magnetic exchange, revealing opportunities for the design of new n d-4f single-molecule magnets and bulk magnetic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Cyclometallated Co(III) Complexes with Lowest-Energy Charge Transfer Excited States Accessible with Visible Light

The Co(III) complexes, cis-[Co(ppy) 2 (L)]PF 6 , where ppy = 2-phenylpyridine and L = bpy (2,2′-bipyridine; 1), phen (1,10-phenanthroline; 2), and DAP (1,12-diazaperylene; 3), are reported and their photophysical properties were investigated to evaluate their potential as sensitizers for applications that include solar energy conversion schemes and photoredox catalysis. Calculations show that cyclometallation in the cis-[Co(ppy) 2 (L)]PF 6 series affords strong Co(dπ)/ppy(π) orbital interactions that result in a Co/ppy(π*) highest occupied molecular orbital (HOMO) and a lowest unoccupied molecular orbital (LUMO) localized on the diimine ligand, L(π*). Complexes 1–3 exhibit relatively invariant oxidation potentials, whereas the reduction event is dependent on the identity of the diimine ligand, L, consistent with the theoretical predictions. For 3 a broad Co/ppy(π*) → L(π*) metal/ligand-to-ligand charge transfer (ML-LCT) absorption band is observed in CH 3 CN with a maxima at 507 nm, extending beyond 600 nm. Upon excitation of the 1ML-LCT transition, transient absorption features consistent with the population of a 3 ML-LCT excited state with lifetimes, τ, of 3.0 ps, 4.6 and 42 ps for 1, 2 and 3 in CH 3 CN respectively are observed. Upon irradiation with 505 nm, 3 is able to reduce methyl viologen (MV 2+ ), an electron acceptor commonly in photocatalytic schemes. Here, to our knowledge, 3 represents the first heteroleptic molecular Co(III) complex that combines cyclometallation with a diimine ligand with lowest-lying metal-to-ligand charge transfer excited states able to undergo photoinduced charge transfer with low-energy green light. As such, the structural design of 3 represents an important step toward d6 photosensitizers based on earth abundant metals.

Energy↗

DOE OSTI2022

Record-quality GaAs two-dimensional hole systems

The complex band structure, large spin-orbit induced band splitting, and heavy effective mass of two-dimensional (2D) hole systems hosted in GaAs quantum wells render them rich platforms to study many-body physics and ballistic transport phenomena. Here we report ultra-high-quality (001) GaAs 2D hole systems, fabricated using molecular beam epitaxy and modulation doping, with mobility values as high as 5.8 × 10 6 cm 2 /(V s) at a hole density of p = 1.3 × 10 11 /cm 2 , implying a mean free path of ≃ 27 μm. In the low-temperature magnetoresistance trace of this sample, we observe high-order fractional quantum Hall states up to the Landau level filling ν = 12/25 near ν = 1/2. Furthermore, we see a deep minimum develop at ν = 1/5 in the magnetoresistance of a sample with a much lower hole density of p = 4.0 × 10 10 /cm 2 where we measure a mobility of 3.6 × 10 6 cm2/(V s). Finally, these improvements in sample quality were achieved by the reduction of residual impurities both in the GaAs channel and in the AlGaAs barrier material, as well as optimization in the design of the sample structure.

2-dimensional systems↗

DOE OSTI2023

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Localization of Nonblinking Point Sources Using Higher-Order-Mode Detection and Optical Heterodyning: Developing a Strategy for Extending the Scope of Molecular, Super-resolution Imaging

While the stochastic, “blinking” nature of fluorescent systems has enabled the super-resolution of their localization by the fitting of their point-spread functions (PSFs), this strategy cannot be exploited for similar resolution of “nonblinking” systems, such as those that might be encountered in a coherent Raman experiment. An alternative method for subdiffraction-limited imaging lies in the exploitation of optical heterodyning. For example, if a Gaussian PSF (a TEM00 mode) of a point emitter is displaced with respect to the origin of the optical system, photons in the higher-order TEM modes carry information about that displacement. Information concerning the displacement can be extracted from photons in these higher-order modes. These photons can be collected by optical heterodyning, which exploits the large gain in a detector’s response to an optical signal from an emitter coupled to a local oscillator, which is prepared in the TEM of interest, e.g., TEM10. We have generalized and developed the heterodyning technique to localize point emitters via the detection of higher-order spatial modes. We have developed a theoretical approach to find a practical estimation limit of the localization parameters using a realistic model that accounts for shot noise, background noise, and Gaussian noise. To demonstrate the applicability of the method, we designed experiments in which a laser is a surrogate for one and two point emitters. Using the Fisher information and its accompanying Cramér-Rao lower bound, we demonstrate super-resolution localization in these cases: we show that objects can be localized to roughly 2–3 orders of magnitude of their point-spread function’s size for a given optical system. Finally and most importantly, it is suggested that the results will ultimately be generalizable to multiple emitters and, most importantly, to “nonblinking” molecular systems, which will be essential for broadening the scope of super-resolution measurements beyond the limits of fluorescence-based techniques.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Chlorine gas and anion radical reactivity in molten salts and the link to chlorobasicity

Next-generation nuclear power plants may include exciting novel designs in which molten salts are the coolant or a combination of the coolant and fuel. Whereas it is straightforward to see why having a low volatility coolant can be advantageous for safety, much is not understood about the production of volatile halogen gases as a result of radiation and even less is known about the distribution of these species at and away from interfaces. Using first principles molecular dynamics simulations, we investigate the product of the disproportionation reaction between chlorine anion radicals (nominally Cl 2 ˙ − ) in the bulk and slab configurations. We find that the product depends on the chlorobasicity of the medium. For example, in ZnCl 2 , Cl 2 forms, but in a eutectic mixture of LiCl and KCl, Cl 3 − is formed as a product. We also find that Cl 3 − prefers to form at the vapor interface and this may have implications for corrosion and reactivity. Furthermore, the mechanisms of the mobility of Cl 2 and Cl 3 − are radically different, the first one being vehicular and the second Grotthus-like. Chlorobasicity is linked to the electronic structure of the host melt; ZnCl 2 forms extended networks along which metal ions and anionic counterions have significant electronic orbital overlap forming long, linear, molecular-like constructs; the opposite is true for the alkali metal eutectic salt.

Nguyen, Hung H. [Univ. of Iowa, Iowa City, IA (Uni↗

DOE OSTI2022

Effects of neoclassical tearing modes and toroidal field ripple on lost alpha power in the SPARC tokamak

Here, using the SPIRAL Monte Carlo, full particle-orbit simulation code, we investigate the effects of neoclassical tearing modes (NTMs) and toroidal field (TF) ripple on alpha power losses during steady-state operation of the SPARC primary reference discharge. Model perturbations for TF ripple and the $m/n = 2/1$ and $3/2$ NTMs with exaggerated widths selected based on an H-mode plasma approaching thermal quench are added to a simulated SPARC magnetic equilibrium through which marker particles are tracked. The 3/2 and 2/1 NTMs are located at $\rho_{\mathrm{pol}} \sim 0.76$ and $\rho_\mathrm{pol} \sim 0.86$ respectively, well positioned to increase alpha particle transport into and within an outer lossy region of the plasma beyond $\rho_\mathrm{pol} \sim 0.8$ where over 95% of lost alpha particles are born. Total alpha power losses are shown to increase modestly from 1.73% lost at a minimum to 2.34% lost at a maximum, and alpha particle surface power densities form localized hotspots on the first-wall near the low-field side midplane due to NTMs and TF ripple. We establish a conservative upper limit for first-wall alpha surface power densities on a toroidally symmetric wall for typical, flattop operation and motivate the consideration of NTMs in the design of three dimensional limiter surfaces for SPARC.

alpha power losses↗

DOE OSTI2020

Quantal diffusion approach for multinucleon transfer processes in the 58,64 Ni + 208 Pb reactions: Toward the production of unknown neutron-rich nuclei

In recent years, substantial efforts have been made for the study of multinucleon transfer reactions at energies around the Coulomb barrier both experimentally and theoretically, aiming at the production of unknown neutron-rich heavy nuclei. It is crucial to provide reliable theoretical predictions based on microscopic theories with sufficient predictive power. Purpose: This article aims to clarify the applicability of the quantal diffusion approach based on the stochastic mean-field (SMF) theory for multinucleon transfer processes. Isotope production cross sections are evaluated for the reactions of 64 Ni + 208 Pb at E c.m. = 268 MeV and 58 Ni + 208 Pb at E c.m. = 270 MeV and are compared with available experimental data. Methods: Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are carried out for a range of initial orbital angular momenta with Skyrme SLy4d functional. Quantal diffusion equations, derived based on the SMF theory, for variances and covariance of neutron and proton numbers of reaction products are solved, with microscopic drift and diffusion coefficients obtained from time evolution of occupied single-particle orbitals in TDHF. Secondary de-excitation processes, both particle evaporation and fission, are simulated by a statistical compound-nucleus de-excitation model, gemini++. Dynamics of a fast isospin equilibration process followed by a slow drift toward the mass symmetry are commonly observed, as expected. Various reaction outcomes are evaluated, including average mass and charge numbers of reaction products, total kinetic energy loss (TKEL), scattering angle, contact time, and production cross sections for primary and secondary products. By comparing with the experimental data, we find that SMF and TDHF quantitatively reproduce experimental data for few-nucleon-transfer channels around the average values. In contrast, for many-nucleon-transfer channels, we find that the SMF approach provides much better description of the experimentally measured isotopic distributions. The results underline the importance of beyond-mean-field effects, especially one-body (mean-field) fluctuations and correlations, in describing multinucleon transfer processes. Moreover, through a combined analysis of SMF with a statistical model, gemini++, we find a significant contribution of transfer-induced fission, which is consistent with the experimental observation. In some cases, the SMF approach overestimates the isotopic width, requiring further improvements of the theoretical description. Possible ways to improve the description are discussed. The SMF approach is designed to describe the quantum many-body problem according to an ensemble of mean-field trajectories, taking into account part of many-body correlations in the description. As it requires feasible computational costs comparable to the ordinary TDHF approach, together with further model improvements, it will be a promising tool in the search for optimal reaction conditions to produce yet-unknown neutron-rich heavy nuclei through the multinucleon transfer reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Electronic structure and magnetism of pristine, defected, and strained $\mathrm{Ti_2N}$ $\mathrm{MX}$ene

Here, from first principles electronic structure calculations, we unravel the evolution of structural, electronic, and magnetic properties of pristine, defected, and strained titanium nitride MXene with different functional groups (-F, -O, -H, and -OH). The formation and cohesive energies reveal their chemical stability. The dynamical stability of Ti 2 N mono-layer is also confirmed by phonon calculations. The MAX phase and defect free functionalized MXenes are metallic except for oxygen terminated (Ti 2 NO 2 ) one which is 100% spin polarized half-metallic ferromagnet. The spin–orbit coupling significantly influences the bare MXene (Ti 2 N) to exhibit Dirac topology and band inversion near the high symmetry directions. The strain effect sways the Fermi level thereby shifting it towards lower energy state under compression and towards higher energy state under tensile strain in Ti 2 NH 2 . The Ti 2 NO 2 exhibits exotic electronic structure not only in pristine but also in strained and defected structures. Its half-metallic nature changes to semi-metallic under 1% compression and it is completely destroyed under 2% compression. In single vacancy defect, its band structure remarkably transforms from half-metallic to semi-conducting with large band gap in 12.5% Ti, weakly semi-conducting in 5.5% Ti, and semi-metallic in 12.5% O. The 25% N defect changes its half-metallic characteristic to metallic. Further, the 12.5% Co substitution preserves its half-metallic character, whereas Mn substitution allows it to convert half-metallic characteristic into weak semi-metallic characteristic preserving ferromagnetism. However, Cr substitution converts half-metallic ferromagnetic state to half-metallic anti-ferromagnetic state. The understanding made here on collective structural stability, and electronic band structure, and magnetic phenomena in novel 2D Ti 2 N derived MXenes open up their possibility in designing them for synthesis.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Structure- and Interaction-Based Design of Anti-SARS-CoV-2 Aptamers

Aptamer selection against novel infections is a complicated and time-consuming approach. Synergy can be achieved by using computational methods together with experimental procedures. In this study, we aim to develop a reliable methodology for a rational aptamer in silico et vitro design. The new approach combines multiple steps: (1) Molecular design, based on screening in a DNA aptamer library and directed mutagenesis to fit the protein tertiary structure; (2) 3D molecular modeling of the target; (3) Molecular docking of an aptamer with the protein; (4) Molecular dynamics (MD) simulations of the complexes; (5) Quantum-mechanical (QM) evaluation of the interactions between aptamer and target with further analysis; (6) Experimental verification at each cycle for structure and binding affinity by using small-angle X-ray scattering, cytometry, and fluorescence polarization. By using a new iterative design procedure, structure- and interaction-based drug design (SIBDD), a highly specific aptamer to the receptor-binding domain of the SARS-CoV-2 spike protein, was developed and validated. The SIBDD approach enhances speed of the high-affinity aptamers development from scratch, using a target protein structure. The method could be used to improve existing aptamers for stronger binding. This approach brings to an advanced level the development of novel affinity probes, functional nucleic acids. It offers a blueprint for the straightforward design of targeting molecules for new pathogen agents and emerging variants.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

Wall heating by subcritical energetic electrons generated by the runaway electron avalanche source *

Abstract Subcritical energetic electrons (SEEs) produced by the runaway electron (RE) avalanche source at energies below the runaway threshold are found to be the primary contributor to surface heating of plasma-facing components (PFCs) during final loss events. This finding is supported by theoretical analysis, computational modeling with the Kinetic Orbit Runaway electrons Code (KORC), and qualitative agreement with DIII-D experimental observations. The avalanche source generates significantly more secondary electrons below the runaway threshold, which thermalize rapidly when well-confined. However, during a final loss event, the RE beam impacts the first wall, and SEEs are deconfined before they can thermalize. Additionally, because the energy deposition length decreases faster than energy, the deposited energy density, and thus the maximum PFC surface temperature change, is larger for SEEs than REs. KORC simulations employ an analytic first wall to model particle deconfinement onto a non-axisymmetric wall composed of individual tiles. PFC surface heating is calculated using a 1D model extended to include an energy-dependent deposition length scale. Simulations of DIII-D qualitatively agree with infrared (IR) imaging only when SEEs from the avalanche source are included. These results demonstrate that SEEs are the dominant contributor to PFC surface heating and indicate that the avalanche source plays a critical role in the PFC damage caused during final loss events. The prominence of SEEs also has important implications for interpreting IR imaging, one of the primary diagnostics for RE-wall interaction diagnosis, despite REs dominating the energy and current density. This result improves predictions of wall damage due to post-disruption REs to estimate material lifetime and design RE mitigation systems for ITER and future reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Assessment of BPM options for the EIC Beam Accumulator Ring

The electron injection system for the Electron-Ion Collider (EIC) at BNL is designed to provide a beam of polarized electrons, which is crucial for studying the structure of protons and atomic nuclei. The Beam Accumulator Ring (BAR) is a part of the injection chain between the 750 MeV linear accelerator and the Rapid Cycling Synchrotron (RCS), which accelerates the beam up to full energy (5–18 GeV). The functional role of the BAR is to accumulate the charge injected from the linear accelerator in order to achieve the high intensity of the polarized electron beam required by the specifications for injection into the RCS. This objective will be realized through the sequential injection of bunches with a charge of 1.1 nC at a repetition rate of 30 Hz. Once the charge of a single bunch reaches 28 nC, the beam will be extracted from the BAR and injected into the RCS. The beam instrumentation needs to provide reliable measurements with the required accuracy over the dynamic range from 0.1 nC (one tenth of a typical injected bunch charge) to 32 nC – the maximum expected accumulated current. There are 10 Beam Position Monitors (BPMs) in the ring, 1 extra button-electrode assembly for the RF system, and 8 BPMs in the beam transport lines. The BPM locations are shown in Fig. 1, marked by blue rectangles. The ring BPMs will be capable of both average orbit and turn-by-turn measurement modes. Accurate measurement of the beam position with a large horizontal offset requires polynomial correction of the BPM nonlinearity.

43 PARTICLE ACCELERATORS↗

DOE OSTI2023

How fluctuation intensity flux drives SOL expansion

Abstract Predictions of heat load widths λ q based on particle orbits alone are very pessimistic. This paper shows that pedestal peeling-ballooning (P-B) magnetohydrodynamic (MHD) turbulence broadens the stable scrape-off layer (SOL) by the transport, or spreading, of fluctuation energy from the pedestal. λ q is seen to increase with Γ ε , the fluctuation energy density flux. We elucidate the fundamental physics of the spreading process. Γ ε increases with pressure fluctuation correlation length. P-B turbulence is seen to be especially effective at spreading, on account of its large effective mixing length. Spreading is shown to be a multiscale process, which is enhanced by the synergy of large and small-scale modes. Pressure fluctuation skewness correlates well with the spreading flux–with the zero crossing of skewness and Γ ε spatially coincident–suggesting the role of coherent fluctuation structures and the presence of intermittency in λ q broadening. λ q ∼ B p − 1 scaling persists for the broadened SOL. We show that the spreading flux increases for increasing pedestal pressure gradient ∇ P 0 and for decreasing pedestal collisionality υ ped ∗ . This trend is due to the dominance of peeling modes for large ∇ P 0 and low υ ped ∗ . Ultimately, we see that a state of weak MHD turbulence, as for small ELMs, is very attractive for heat load management. Our findings have transformative implications for future fusion reactor designs and call for experimental investigations to validate the observed trends.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Energy Frontier Research Centers: Center for the Computational Design of Functional Layered Materials (CCDM) August 1, 2014 - July 31, 2018; Center for Complex Materials from First Principles (CCM) August 1, 2018 - July 31, 2021 (Final Report)

The mission of the DOE Energy Frontier Research Centers CCDM (2014-2018) and CCM (2018-2021) was to theoretically develop, computationally apply, and experimentally validate electronic structure methods for all materials, with a focus on the complex materials, especially layered and two-dimensional materials, strongly-correlated materials, and liquid water. This was achieved by over 200 published journal articles authored by about 17 senior investigators from physics and chemistry and from theory, computation, and experiment, plus their collaborators. In particular, the Centers confirmed the predictive power of the SCAN (strongly constrained and appropriately normed) density functional, which was constructed to satisfy 17 known exact constraints and several appropriate norms. Without being fitted to real bonded systems, and at a modest computational cost, SCAN correctly predicted covalent, ionic, metallic, hydrogen, and van der Waals bonds in many challenging materials. SCAN gave an improved description of defects in semiconductors, surface properties of metals, seven phases of ice, liquid water, liquid and supercooled silicon, subtle structural distortions in ferroelectrics, formation energies and structural predictions for solids, and critical pressures for structural phase transitions. Perhaps most remarkably, SCAN correctly described some strongly-correlated materials that were previously believed to be beyond the reach of density-functional approximations. SCAN is the only density functional that correctly predicts the band gap closing under chemical doping of the cuprate high-temperature superconducting materials. SCAN also predicts a landscape of competing stripe and magnetic phases in the cuprates. For some materials with some codes, SCAN has convergence problems that are greatly reduced by the CCM-developed r 2 SCAN, without loss of accuracy or rigor. SCAN and r 2 SCAN still make some self-interaction error, which is greatly reduced by the CCDM/ CCM-developed local orbital scaling correction (LOSC). These Centers further proved that the fundamental energy gaps of a solid from an orbital energy difference and from total energy differences are the same for a large class of generalized Kohn-Sham (GKS) functionals, including SCAN and standard hybrid functionals, and that symmetry breaking arises when a dynamic density fluctuation drops to zero frequency. The Centers identified new mechanisms for catalysis in layered materials with ions intercalated between the layers, investigated charge density waves both in model systems and in real layered materials, studied changes of band gap with the number of layers, and explored topological ultrathin films, bent nanoribbons, and defects.

08 HYDROGEN↗

DOE OSTI2020

Au 3 -to-Ag 3 coordinate-covalent bonding and other supramolecular interactions with covalent bonding strength

An efficient strategy for designing charge-transfer complexes using coinage metal cyclic trinuclear complexes (CTCs) is described herein. Due to opposite quadrupolar electrostatic contributions from metal ions and ligand substituents, [Au(m-Pz-(i-C 3 H7) 2 )] 3 ∙[Ag(m-Tz-(n-C 3 F 7 ) 2 )] 3 (Pz = pyrazolate, Tz = triazolate) has been obtained and its structure verified by single crystal X-ray diffraction – representing the 1st crystallographically-verified M 3 @M' 3 stacked adduct of monovalent coinage metal CTCs. Abundant supramolecular interactions with aggregate covalent bonding strength arise from a combination of M–M' (Au / Ag), metal–π, π–π interactions and hydrogen bonding in this charge-transfer complex, according to density functional theory analyses, yielding a computed binding energy of 66 kcal mol -1 between the two trimer moieties – a large value for intermolecular interactions between adjacent d 10 centres (nearly doubling the value for a recently claimed Au(I) / Cu(I) polar-covalent bond: Proc. Natl. Acad. Sci. U.S.A., 2017, 114, E5042) – which becomes 87 kcal mol -1 with benzene stacking. Surprisingly, DFT analysis suggests that: (a) some other literature precedents should have attained a stacked M 3 @M' 3 product akin to the one herein, with similar or even higher binding energy; and (b) a high overall intertrimer bonding energy by inferior electrostatic assistance, underscoring genuine orbital overlap between M and M' frontier molecular orbitals in such polar-covalent M–M' bonds in this family of molecules. The Au / Ag bonding is reminiscent of classical Werner-type coordinate-covalent bonds such as H 3 N: → Ag in [Ag(NH 3 ) 2 ] + , as demonstrated herein quantitatively. Solidstate and molecular modeling illustrate electron flow from the p-basic gold trimer to the p-acidic silver trimer with augmented contributions from ligand-to-ligand’ (LL'CT) and metal-to-ligand (MLCT) charge transfer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

CRADA Final Report: CRADA Number NFE-22-09330 with General Fusion

General Fusion is developing a magnetized target fusion (MTF) approach that involves compressing an initial magnetically confined plasma inside a cavity formed in liquid metal. This approach builds from concepts initially developed under the Linus program at the U.S. Naval Research Laboratory and combines it with advances from compact toroid experiment (CTX) and sustained spheromak plasma experiment (SSPX) in compact toroid plasmas and coaxial Marshall gun systems. Modeling the tokamak during compression is central to designing a successful MTF device. The plasma is formed by coaxial helicity injection in the General Fusion device. Immediately after formation, the plasma has a diverted tokamak configuration with a single null. As the wall moves inwards, the plasma is repelled from the conducting surface and driven inwards by currents induced by its magnetic field in the liquid metal wall. As the liquid metal closes (or bridges) the opening of the coaxial plasma injector, the magnetic field topology alters to remove the null. Due to this, the plasma moves from a diverted to a wall-limited configuration. The liquid metal liner continues to close in and change shape, reducing in radius by a factor of ten at the peak of plasma compression. A model of the MTF plasma must be able to handle this continually varying geometry, and to be predictive, it must faithfully include the real imperfections arising in the process. In this project, we pursued a Monte Carlo approach to closures for MHD by computing kinetic electron trajectories in an MHD plasma background from simulations of GF devices. This requires enhancing the capabilities of the KORC-T code for running large ensembles of kinetic trajectories by porting it to GPU architectures and enabling workflows for large ensembles on OLCF machines. With these capabilities, it is possible to produce a large library of kinetic calculations of electron orbits evolving in plasma configurations spanning the magnetic configurations and plasma density profiles, including non-axisymmetry, arising in the General Fusion’s existing PI3 spherical tokamak device. Using ensembles will capture particles passing a single point in space in a given magnetic configuration, and the entire dataset will cover a range of global magnetic field geometries. By sampling around many starting points, this dataset will capture the spatial dependence of the plasma parameters. From this large dataset, it is possible to produce a reduced model for the kinetic effects not captured in MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗