Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic structure generation”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Tailoring Growth Interfaces of Virtual Substrates for Power Electronics

Power electronics materials are poised to play a critical role in fulfilling next generation energy needs, with up to 90% of future energy demand predicted to flow through power electronics at some point. AlxGa1-xN ranks high among candidate materials, having bipolar dopability, thermal and chemical stability and an ultra-wide bandgap. However, AlGaN growth is limited by a lack of lattice-matched substrates, ultimately stunting material quality at higher thicknesses needed for power electronics applications. Further, high power applications increasingly call for fully vertical device structures, necessitating a conductive substrate. Recently our group identified the (111) plane of TaC as a conductive surface lattice-matched to Al0.55Ga0.45N, taking inspiration from prior work of AlN and GaN binaries on carbide and boride substrates. In this talk we demonstrate the growth of (111)-oriented TaC by RF sputtering. We investigate the interface of TaC with sapphire and SiC substrates and identify means to suppress competing Ta2C nucleation in order to stabilize (111)-oriented TaC. Potential stacking sequences are identified with respect to crystal structure and observed twinning in the TaC films. We next assess structural changes and film recrystallization that results from face-to-face annealing of TaC thin films at high temperatures above 1500 degrees C. Changes to grain structure and domain size are assessed by x-ray diffraction and surface morphology is explored using atomic force microscopy. Figure 1 shows significant improvements to in- and out-of-plane strain following annealing along with the formation of terraced step edges at the film surface. Strain as a function of material composition and thickness is considered, as this may play a major role in future nucleation of AlGaN layers. (1) R. J. in a face-to-face configuration, as illustrated in the schematic at left. Kaplar et al 2017, ECS J. Solid State Sci. Technol. 6 Q3061; (2) D. M. Roberts et al 2022, https://arxiv.org/abs/2208.11769; (3) T. Aizawa et al 2008, J Crys Growth 310, 1 22; (4) R. Liu et al 2002, Appl. Phys. Lett. 81, 3182-3184.

ENGINEERING↗
DOE OSTI2026

HydraGNN_Predictive_GFM_2026 - Ensemble of predictive graph foundation models for atomistic materials modeling

This release contains data and parameters of HydraGNN-based graph foundation models trained as a result of the work published in the pre-print "Exascale Multi-Task Graph Foundation Models for Imbalanced, Multi-Fidelity Atomistic Data" by M. Lupo Pasini et al. (https://arxiv.org/abs/2604.15380). We jointly train on 16 open first-principles datasets (544+ million structures covering 85+ elements) using a multi-task architecture with per-dataset heads and a scalable ADIOS2/DDStore data pipeline. On Frontier, we execute six large-scale DeepHyper hyperparameter optimization campaigns in FP64 and promote the top-performing message-passing models to sustained 2,048-node training, yielding a PaiNN-based lead model. The version of HydraGNN used to generate the outputs provided in this release is HydraGNN v5.0 (https://github.com/ORNL/HydraGNN/releases/tag/v5.0) The list of datasets used for the training of the graph foundation model is the following: 1) Alexandria [1] 2) ANI1x [2] 3) MPTrj [3] 4) Open Catalyst 2020 (OC20) [4] 5) Open Catalyst 2022 (OC22) [5] 6) Open Catalyst 2025 (OC25) [6] 7) Open Direct ir Capture 2023 (ODAC23) [7] 8) Open Materials 2024 (OMat24) [8] 9) Open Molecules 2025 (OMol25) [9] 10) OMol25-neutral (subset of OMol25 that contains only molecules with zero total charge) 11) OMol25-non-neutral (subset of OMol25 that contains only molecules with non-zero total charge) 12) Open Polymers 2026 (OPoly2026) [10] 13) Nabla2DFT [11] 14) QCML [12] 15) QM7X [reference 13] 16) transition1x [14] Dataset references: [1] J. Schmidt et al., “A dataset of 175k stable and metastable materials calculated with the PBEsol and SCAN functionals,” Scientific Data, vol. 9, p. 64, 2022. [2] J. S. Smith et al., “The ANI-1ccx and ANI-1x data sets, coupled-cluster and density functional theory properties for molecules,” Scientific Data, vol. 7, p. 134, 2020. [Online]. Available: https: //www.nature.com/articles/s41597-020-0473-z [3] A. Jain et al., “Commentary: The Materials Project: A materials genome approach to accelerating materials innovation,” APL Materials, vol. 1, no. 1, p. 011002, 07 2013. [Online]. Available: https://doi.org/10.1063/1.4812323 [4] L. Chanussot et al., “Open catalyst 2020 (oc20) dataset and community challenges,” ACS Catalysis, vol. 11, no. 10, pp. 6059–6072, 2021. [Online]. Available: https://doi.org/10.1021/acscatal.0c04525 [5] K. Tran et al., “Open catalyst 2022 (oc22) dataset and challenges for oxidation electrocatalysts,” ACS Catalysis, vol. 13, no. 5, pp. 3066–3084, 2023. [Online]. Available: https://doi.org/10.1021/acscatal.2c05426 [6] S. J. Sahoo et al., “The open catalyst 2025 (oc25) dataset and models for solid-liquid interfaces,” arXiv preprint arXiv:2509.17862, 2025. [Online]. Available: https://arxiv.org/abs/2509.17862 [7] A. Sriram et al., “The open DAC 2023 dataset and challenges for sorbent discovery in direct air capture,” ACS Central Science, vol. 10, no. 5, pp. 923–941, 2024. [8] L. Barroso-Luque et al., “Open materials 2024 (omat24) inorganic materials dataset and models,” 2024. [Online]. Available: https://arxiv.org/abs/2410.12771 [9] D. S. Levine et al., “The open molecules 2025 (OMol25) dataset, evaluations, and models,” 2025. [Online]. Available: https://arxiv.org/abs/2505.08762 [10] D. S. Levine et al., The open polymers 2026 (OPoly26) dataset and evaluations,” arXiv preprint arXiv:2512.23117, 2025. [Online]. Available: https://arxiv.org/abs/2512.23117 [11] K. Khrabrov et al., “Nabla2dft: A universal quantum chemistry dataset of drug-like molecules and a benchmark for neural network potentials,” in NeurIPS 2024 Datasets and Benchmarks Track, 2024. [Online]. Available: https://openreview.net/forum?id=ElUrNM9U8c [12] S. Ganscha et al., “The QCML dataset, quantum chemistry reference data from 33.5M DFT and 14.7B semi-empirical calculations,” Scientific Data, vol. 12, p. 406, 2025. [13] J. Hoja et al., “QM7-X, a comprehensive dataset of quantum-mechanical properties spanning the chemical space of small organic molecules,” Scientific Data, vol. 8, p. 43, 2021. [Online]. Available: https://www.nature.com/articles/s41597-021-00812-2 [14] M. Schreiner et al., “Transition1x - a dataset for building generalizable reactive machine learning potentials,” Scientific Data, vol. 9, p. 779, 2022. The folder "datasets_ADIOS2_format" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "datasets_ADIOS2_format" directory contains 2 sub-directories, one for the version "v1" of the datasets and one for the version "v2" of the datasets. The version "v1" of the datasets provides values of the total energy as they are extracted from the original data as it was released by the respective institutions. The version "v2" of the datasets provides values of the energy that have been realigned. The realignment was performed by training a linear regression model that predicts the total energy as a function of the chemical composition of the atomistic structure, and then subtract such prediction from the original value of the total energy. Both folders "v1" and "v2" contain 16 sub-directories, each corresponding to an ADIOS2-formatted dataset The folder "DeepHyper-results" contains the configurational files and model's parameters for all the 186 HPO trials that were successfully completed by the scalable hyperparameter optimization (HPO) runs on Frontier. The content of the folder "DeepHyper-results" I structured as follows: 1) task-list.txt: list of mpnn name, jobid, and deephyper task id 2) gfm_${MPNN}_${JOBID}_0.${TASKID}: run directory with checkpoint files 3) gfm_${MPNN}: deephyper summary directory (*.csv) for each specific MPNN type 4) deephyper-experiment-${JOBID}: output and error logs for each job The file "deephyper-sorted.csv" contains the details of each HydraGNN model built and tested by HPO, obtained by merging the (*.csv) filed from each HPO run executed. Out of all the HPO trials, we selected 10 to continue the training of the respective HydraGNN models. Due to limited computational budget available in the LRN070 allocation we could not complete the training till convergence for all these 10 selected models. The folder "models" contains multiple sub-folders, one per each HydraGNN model trained. Each model sub-folder contains the parameters of each HydraGNN model, with multiple checkpoint-restarts. The list of sub-folders are as follows: 1) multidataset_hpo-BEST1-fp64 2) multidataset_hpo-BEST2-fp64 3) multidataset_hpo-BEST3-fp64 4) multidataset_hpo-BEST4-fp64 5) multidataset_hpo-BEST5-fp64 6) multidataset_hpo-BEST6-fp64 7) multidataset_hpo-BEST7-fp64 8) multidataset_hpo-BEST8-fp64 9) multidataset_hpo-BEST9-fp64 10) multidataset_hpo-BEST10-fp64 Within each one of these folders, additional auxiliary log files are provided with descriptions about how the training proceeded. The lead PaiNN-model is contained inside "multidataset_hpo-BEST6-fp64". The file "mlp_branch_weights" contains the parameters of the multi-layer perceptron (MLP) used to reconcile the predictions of the 16 output decoding heads of the HydragNN architectures. The MLP takes in input the chemical composition of the atomistic structure and predicts averaging weights to linearly mix the predictions of each output decoding head toward consolidating them into a single one. The folder "1.1billion-structure-inference" contains 1.1 billion atomistic structures randomly generated. Each structures is associated with energy and forces predicted with the lead-PaiNN model combined with the MLP model for reconciliation of the multi-branch predictions generated by the 16 output decoding heads. The folder "1.1billion-structure-inference" contains 9,300 (*.tar.gz) subdirectories, one per Frontier compute node used to execute the inference at exascale. Once uncompressed, each (*.tar.gz) subdirectory contains an ADIOS2 (*.bp) file container, where each atomistic structure is stored as a PyTorch-Geometric Data object. The file "export_dataset_environment_variables.sh" contains the environment variables that need to be set before running the HydraGNN code to reproduce the results provided in this dataset release. The code that can be used to load the ADIOS2 files, load HydraGNN models, and run inference is available at: https://github.com/ORNL/HydraGNN/releases/tag/v5.0

36 MATERIALS SCIENCE↗
DOE OSTI2024

MolSym : A Python package for handling symmetry in molecular quantum chemistry

A consideration of the point group symmetry of molecules is often advantageous from a computational efficiency standpoint and sometimes necessary for the correct treatment of chemical physics problems. Many modern electronic structure software packages include a treatment of symmetry, but these are sometimes incomplete or unusable outside of that program’s environment. Therefore, we have developed the MolSym package for handling molecular symmetry and its associated functionalities to provide a platform for including symmetry in the implementation and development of other methods. Features include point group detection, molecule symmetrization, arbitrary generation of symmetry element sets and character tables, and symmetry adapted linear combinations of real spherical harmonic basis functions, Cartesian displacement coordinates, and internal coordinates. We present some of the advantages of using molecular symmetry as achieved by MolSym, particularly with respect to Hartree–Fock theory, and the reduction of finite difference displacements in gradient/Hessian computations. Furthermore, this package is designed to be easily integrated into other software development efforts and may be extended to further symmetry applications.

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

Heavy boron doping effects on biaxially tensile strained germanium (>1.5%) investigated via structural characterization, effective lifetime assessment and atomistic modeling

Highly tensile strained germanium (ε-Ge) represents an essential material system for emerging electronic and photonics applications. Moreover, adjusting the doping levels to moderate or high concentrations can effectively tailor the properties of ε-Ge for specific applications. This article combines experimental characterization with a theoretical framework to examine the effects of heavy elemental boron (B) doping on pseudomorphic sub-50 nm ε-Ge. High resolution X-ray diffractometry is used to validate tensile strain levels of 1.53% and 1.68% in Ge epilayers, surpassing the indirect-to-direct band gap crossover point at ∼1.5% biaxial tensile strain. Cross-sectional transmission electron microscopy revealed visual evidence of stacking faults and surface roughening in 1.68% ε-Ge, although a coherent and abrupt Ge/III–V heterointerface is observed, devoid of interfacial misfit dislocations. Effective lifetime measurements demonstrated approximately twofold enhancement in 1.53% B-doped ε-Ge (N B ∼7 × 10 19 cm −3 ) compared to its unstrained B-doped counterpart, while no such improvement was observed in 1.68% B-doped ε-Ge. This lack of enhancement is attributed to the presence of stacking faults and surface roughness within the ε-Ge epilayer. Through density functional theory calculations, we independently demonstrate that substitutional B atoms induce local deformation of Ge–Ge bonds in both unstrained Ge and ε-Ge epilayers, resulting in an additive tensile strain. This phenomenon could potentially lead to dynamic reduction and overcoming of the critical layer thickness for the system, facilitating the nucleation and subsequent glide of 90° leading Shockley partial dislocations, thereby generating stacking faults. In essence, these findings establish an upper limit on the B-doping concentration that can be achieved in highly ε-Ge epilayers, and collectively, offer valuable insights into the significance of heavy doping in Ge-based heterostructures. As such, this study delineates a fundamental constraint for integrating heavily doped ε-Ge in high-performance optoelectronic systems, necessitating precise strain-doping co-optimization to avoid performance degradation.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Local distortion driven magnetic phase switching in pyrochlore Yb 2 ( Ti 1 - x Sn x ) 2 O 7

While it is commonly accepted that the disorder induced by magnetic ion doping in quantum magnets usually generates a rugged free-energy landscape resulting in slow or glassy spin dynamics, the disorder/distortion effects associated with nonmagnetic ion sites doping are still illusive. Here, using AC susceptibility measurements, we show that the mixture of Sn/Ti on the nonmagnetic ion sites of pyrochlore Yb 2 ⁢(Ti 1-x⁢ Sn x ) 2⁢ O 7 induces an antiferromagnetic ground state despite both parent compounds, Yb 2 ⁢Ti 2 ⁢O 7 and Yb 2 ⁢Sn 2⁢ O 7 , order ferromagnetically. Local structure studies through neutron total scattering reveals the local distortion in the nonmagnetic ion sites and its strong correlation with the magnetic phase switching. Further, our study demonstrates the local distortion as induced by the nonmagnetic ion site mixture could be a new path to achieve magnetic phase switching, which has been traditionally obtained by external stimuli such as temperature, magnetic field, pressure, strain, light, etc.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

An Investigation of Oxides of Tantalum Produced by Pulsed Laser Ablation and Continuous Wave Laser Heating

Recent progress has seen multiple Ta 2 O 5 polymorphs generated by different synthesis techniques. However, discrepancies arise when these polymorphs are produced in widely varying thermodynamic conditions and characterized using different techniques. This work aimed to characterize and compare Ta 2 O 5 particles formed at high and low temperatures using nanosecond pulsed laser ablation (PLA) and continuous wave (CW) laser heating of a local area of tantalum in either air or an 18 O 2 atmosphere. Scanning electron microscopy (SEM) and Raman spectroscopy of the micrometer-sized particles generated by PLA were consistent with either a localized amorphous Ta 2 O 5 phase or a similar, but not identical, crystalline β-Ta 2 O 5 phase. The Raman spectrum of the material formed at the point of CW laser impingement was in good agreement with the previously established ceramic “H-Ta 2 O 5 ” phase. TEM and electron diffraction analysis of these particles indicated the phase structure matched an oxygen-vacated superstructure of monoclinic H-Ta 2 O 5 . Further from the point of laser impingement, CW heating produced particles with a Raman spectrum that matched β-Ta 2 O 5 . We confirmed that the high-temperature ceramic phase characterized in previous work by Raman spectroscopy was the same monoclinic phase characterized in different work by TEM and could be produced by direct laser heating of metal in air.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Entanglement Generation in Weakly Driven Arrays of Multilevel Atoms via Dipolar Interactions

We investigate the driven-dissipative dynamics of multilevel atomic arrays interacting via dipolar interactions at subwavelength spacings. Unlike two-level atoms in the weakly excited regime, multilevel atoms can become strongly entangled. Here, the entanglement manifests as the growth of spin waves in the ground-state manifold and survives after turning off the drive. We propose the 2.9 μ⁢m transition between 3 P 2 ↔ 3 D 3 in 88 Sr with 389 nm trapping light as a platform to test our predictions and explore many-body physics with light-matter interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Evaluating the Chemical Reactivity of DFT-Simulated Liquid Water with Hydrated Electrons via the Dual Descriptor

Modeling the various properties of liquid water, particularly its reactivity, has been a longstanding challenge for simulation methods. Recently, ab initio simulations based on density functional theory (DFT) have come to the fore as tenable methods for calculating the properties and reactivity of water, with varying degrees of success for different exchange-correlation functionals. In particular, hybrid-GGA and meta-GGA functionals have been shown to reproduce many of the structural, dynamical, and energetic properties of water to a high degree of accuracy relative to their computational cost. Here, we show that the dual descriptor (DD) measure of nucleophilicity and electrophilicity, which is sometimes used to elucidate organic chemistry reaction mechanisms, can also be used to characterize the reactivity of DFTsimulated liquid water. The DD is especially apt for understanding the reactivity of excess electrons with water as its calculation explicitly involves adding and removing an excess electron from a reference system. We use the DD to explore the reactivity of water simulated using three different DFT functionals: the LDA functional (LDA), a hybrid-GGA functional (PBE0), and a hybrid meta- GGA functional (SCAN0). Using the DD, we show that the SCAN0 functional with the standard 25% Hartree–Fock exchange produces simulated liquid water with many regions that are far more reactive than either PBE0 or LDA. To understand the implications of these highly reactive regions, we then add a strong nucleophile in the form of an excess electron and find that although PBE0 and LDA predict stable hydrated electrons, the excess electron reacts nearly instantaneously with SCAN0 water via proton abstraction to form a hydrogen atom and hydroxide ion. We show that the DD provides the ability to not only predict whether or not liquid water will react with a hydrated electron but also which particular waters will be involved solely from analyzing pure water configurations generated with each functional. We rationalize this result in terms of the known trap-seeking behavior of injected hydrated electrons, which are able to find the most electronegative region in bulk water. These results highlight the utility of the dual descriptor as a fast and interpretable method for investigating condensed-phase reactivity with excess electrons.

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

Illuminating the Material World: Autonomous Microscopy to Understand Order, Disorder, and Everything In Between

Artificial intelligence (AI) holds immense promise for revolutionizing microscopy, yet its widespread adoption has been hindered by challenges ranging from user inexperience to limited model transferability and difficulties in operationalizing machine learning. This presentation showcases our approach to developing practical autonomy for materials discovery, aiming to accelerate the integration of AI into everyday microscopy workflows. As shown in Fig. 1, I will focus on three key areas: understanding order-disorder transitions, quantifying point defects, and achieving truly device-scale microscopy. First, I will demonstrate the power of multi-modal knowledge graphs for integrating diverse microscopy data. By combining imaging, spectroscopy, and diffraction data, these graphs provide a holistic view of material behavior, capturing the intricate relationships between different modalities [1,2]. I will present a case study on how these models illuminate the structural and chemical changes associated with irradiation in oxide thin films, revealing critical insights for designing materials for extreme environments like spaceflight and nuclear energy. Specifically, I will show how multi-modal analysis clarifies the evolution of order-disorder transitions under irradiation, a key factor influencing material performance in these applications. Next, I will address the challenge of quantifying point defects in 2D materials. We demonstrate the application of computer vision and transfer learning to accurately identify and classify various defect types, such as vacancies and substitutional atoms, and to quantify their concentrations. This information is crucial for understanding and tailoring the properties of 2D materials for applications in electronics, optoelectronics, and catalysis. For example, I will show how our models can characterize the topological distribution of point defects in MXene transition metal carbides, providing valuable insights for optimizing their performance in energy storage and separation science. Finally, I will discuss our progress toward autonomous device-scale microscopy [3,4]. We are fundamentally redesigning electron microscopes around the principles of machine reasoning, enabling automation beyond basic tasks like sample navigation and data acquisition to include sophisticated experimental design. This approach paves the way for truly reproducible and massively scaled analysis campaigns. I will emphasize the importance of autonomous microscopy platforms for high-throughput materials discovery and characterization, facilitating the rapid screening of materials for a broad range of applications and accelerating the development of next-generation technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Pedestal origin and extrapolation of high-density small edge-localised-modes peak parallel energy fluence in ITER and SPARC

Experimental analysis and simulations with the BOUT++ code show that small edge-localised modes (ELMs) in reactor-relevant high-density regimes originate in a region close to the separatrix and only marginally perturb the pedestal structure. The measured divertor peak parallel energy fluence (ε ∥,peak ) for a database of small ELM scenarios in DIII-D and ASDEX Upgrade can be reproduced, within 40 % accuracy on average, if an ad hoc modification of the Eich peak parallel ELM energy fluence model is applied to account for the small ELM pedestal birth location. This allows for first-order extrapolation of small-ELM divertor ε ∥,peak to ITER and SPARC, resulting in values that satisfy the nominal melting threshold of tungsten monoblocks of 12 MJ m −2 . The findings reported in this study, both via modelling and direct measurements, constitute a step forward in assessing small ELMs in high edge-collisionality scenarios as a viable plasma regime for the operation of next-generation fusion machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A tutorial on high-order harmonic generation in atoms, molecules, and condensed matter

This tutorial introduces strong-field-driven high-order harmonics, their experimental generation and characterization techniques, and their main applications including attosecond pulse generation and ultrafast spectroscopy of the target material. We begin from the use of atomic targets, where the first high-order harmonic generation (HHG) experiments were realized in the late 1980s. Then, we briefly discuss the basics of the microscopic generation mechanism and how various steps of the mechanism were exploited in applications such as generating isolated attosecond pulses and probing molecular orbitals. We introduce and describe the standard experimental approaches for condensed phase HHG, where we discuss unique technical challenges of the use of solid-state materials, such as the mitigation of plasma formation and laser damage. We cover the fundamentals of high-harmonic spectroscopy in condensed matter systems, such as wide bandgap dielectrics, semiconductors, liquid media, and 2D-crystals. We provide some examples of rapidly emerging spectroscopic capabilities, such as for probing crystal symmetries, Berry phases, and associated non-trivial topological properties of the source material. Finally, we provide an overview of the research field, including some of the challenges, opportunities, and open questions.

Attosecond pulses↗
DOE OSTI2026

Disentangling High Harmonic Generation from Surface and Bulk States of a Topological Insulator

The discovery of topological phases has introduced a new dimension to materials science. Three-dimensional (3D) topological insulators (TIs) are a remarkable class of matter that is insulating in the bulk while hosting conductive topological surface states (TSSs) with unique charge and spin properties. High-order harmonic generation (HHG) has emerged as a powerful tool to probe condensed matter systems by providing insights into their electronic structure and dynamic behavior. Here, we investigate HHG in the prototype 3D-TI Bi$_2$Se$_3$. We demonstrate that the contributions of bulk and surface states to the harmonic emission can be controlled by tuning the thickness of thin film samples. An ultrathin (6 nm) film substantially enhances HHG from the surface states, while the bulk states dominate HHG in a thicker (50 nm) film. By applying a quasi-static terahertz perturbing field, we disentangle the bulk and surface responses and reveal the significant impact of the surface states' shift vector and Berry curvature on HHG. Our study provides effective methods for isolating the optical responses of TSSs from those of the bulk, which opens the door to resolving an ongoing debate regarding whether it is possible to reliably extract topological signatures in HHG.

Atomic Physics (physics.atom-ph)↗
DOE OSTI2024

Generation of shear flows induced by AE / EPM in LHD plasma

The generation of shear flows (SFs) by Alfven Eigenmodes (AEs) and energetic particle modes (EPMs) have important effects on the operation of future nuclear fusion reactors, because SFs regulate the saturation of the AEs/EPMs, the transport of EPs and thermal plasma, as well as the formation of transport barriers among other consequences. The aim of this study is the analysis of SFs generation during the saturation phase of AEs and EPMs in LHD plasma. Experiments performed in the 23rd and 24th LHD experimental campaigns are dedicated to explore the destabilization of AEs/EPMs in discharges with different heating patterns, thermal plasma and magnetic field configurations. In particular, the shots 176490 and 179697 show the destabilization of MHD bursts and energetic-ion-driven resistive interchange modes (EIC), respectively. Charge exchange spectroscopy measurements in both discharges indicate that the generation of SFs by AE/EPM is uncorrelated with the perturbation induced by the neutral beam injector (NBI). Nonlinear simulations performed using the gyro-fluid code FAR3d show the generation of zonal structures, especially SFs, induced during the saturation phase of Toroidal Alfven Eigenmodes (TAEs) triggered in the MHD burst as well as by the 1/1 EIC in the bursting phase. The simulations indicate that SFs are caused by the radial electric fields powered by energy transfers from the unstable AE/EPM towards the thermal plasma. The strongest SFs are measured during the EIC bursting phase once the 1/1 EPM overlaps with nearby resonances at the plasma periphery. Likewise, the largest SFs during the MHD burst are observed once TAEs radially overlap in the inner-middle plasma region.

AE↗
DOE OSTI2025

Gyrokinetic simulations of the effects of magnetic islands on microturbulence in KSTAR

Gyrokinetic simulations are utilized to study effects of magnetic islands (MIs) on the ion temperature gradient (ITG) turbulence in the KSTAR tokamak with resonant magnetic perturbations. Simulations show that the transport is controlled by the nonlinear interactions between the ITG turbulence and self-generated vortex flows and zonal flows, leading to an anisotropic structure of fluctuation and transport on the poloidal plane and in the toroidal direction. MIs greatly enhance turbulent transport of both particle and heat. The turbulent transport exhibits variations in the toroidal direction, with transport through the resonant layer near the island X-point being enhanced when the X-point is located at the outer mid-plane. A quantitative agreement is shown between simulations and KSTAR experiments in terms of time frequency and perpendicular wavevector spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Linear and Nonlinear X-ray Spectra of Chiral Molecules: X-ray Circular Dichroism, Sum- and Difference-Frequency Generation of Fenchone and Cysteine

Recent advancements in X-ray light sources enable element-sensitive nonlinear spectroscopies for probing molecular chirality. We simulate X-ray absorption spectroscopy, X-ray circular dichroism (XCD), and nonlinear optical/X-ray SFG and DFG (OX SFG/DFG) signals for two prototypical chiral molecules, fenchone and cysteine. Furthermore, our multireference simulations reproduce experimental data and reveal how novel X-ray spectroscopies exploit the site- and element-sensitivity of X-rays to uncover molecular asymmetry. The XCD spectra show strong asymmetries at chiral centers, while distant atoms contribute weaker signals. The OX SFG/DFG signals, under fixed resonant optical excitation, strongly depend on the core and valence excitations. This allows us to introduce two-dimensional (2D) chirality-sensitive valence-core spectroscopy, which provides insight into the overlap between valence orbitals and local molecular asymmetry. Finally, our estimates using realistic laser and X-ray pulse parameters demonstrate that such nonlinear experiments are feasible at XFELs, offering a promising tool for investigating chiral molecules.

Circular dichroism spectroscopy↗
DOE OSTI2024

Electronic excitations and spin interactions in chromium trihalides from embedded many-body wavefunctions

Although chromium trihalides are widely regarded as a promising class of two-dimensional magnets for next-generation devices, an accurate description of their electronic structure and magnetic interactions has proven challenging to achieve. Here, we quantify electronic excitations and spin interactions in CrX 3 (X = Cl, Br, I) using embedded many-body wavefunction calculations and fully generalized spin Hamiltonians. We find that the three trihalides feature comparable d-shell excitations, consisting of a high-spin 4 A 2 ($t^3_{2g}$$e^0_g$) ground state lying 1.5–1.7 eV below the first excited state 4 T 2 ($t^2_{2g}$$e^1_g$). CrCl 3 exhibits a single-ion anisotropy A sia = - 0.02 meV, while the Cr spin-3/2 moments are ferromagnetically coupled through bilinear and biquadratic exchange interactions of J 1 = - 0.97 meV and J 2 = - 0.05 meV, respectively. The corresponding values for CrBr 3 and CrI 3 increase to A sia = -0.08 meV and A sia = - 0.12 meV for the single-ion anisotropy, J 1 = -1.21 meV, J 2 = -0.05 meV and J 1 = -1.38 meV, J 2 = -0.06 meV for the exchange couplings, respectively. We find that the overall magnetic anisotropy is defined by the interplay between A sia and A dip due to magnetic dipole–dipole interaction that favors in-plane orientation of magnetic moments in ferromagnetic monolayers and bulk layered magnets. The competition between the two contributions sets CrCl 3 and CrI 3 as the easy-plane (A sia + A dip >0) and easy-axis (A sia + A dip <0) ferromagnets, respectively. The differences between the magnets trace back to the atomic radii of the halogen ligands and the magnitude of spin–orbit coupling. Our findings are in excellent agreement with recent experiments, thus providing reference values for the fundamental interactions in chromium trihalides.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Effect of controlled magnetic island bifurcation on electron diffusion

Magnetic islands strongly influence cross-field electron transport in magnetized plasmas. In particular, bifurcations of the island topology modify the number and location of O-points, X-points, and separatrix boundaries, thereby altering diffusion pathways. In recent DIII-D experiments, external magnetic perturbations were used to rotate and periodically bifurcate the island on the q = 2 surface, causing a switchback between a q = 2/1-dominated structure and a narrower q = 4/2-dominated structure. To investigate how this topological change affects electron transport, we employ the field line tracing code TRIP3D with an implemented collisional operator. Thermal, tracer electrons launched from O-points, X-points, and outside separatrix boundaries reveal distinct diffusion regimes, including classical, subdiffusive, and superdiffusive behavior, depending on both the dominant island mode and launch location. These results suggest that island bifurcation can alter electron diffusion across rational surfaces, with direct implications for particle confinement. While the present work emphasizes diffusion as a general framework, the findings provide insight into the conditions under which electron trapping into an island or stochastization of the island's separatrix can enable additional mechanisms, such as the generation of energetic electrons.

DIII-D↗
DOE OSTI2025

A code-to-code benchmark for magneto-convection in a horizontal duct

Liquid metals and magnetic fields are used in many technical applications such as metallurgy, crystal growth and nuclear fusion reactors. When an electrically conducting fluid moves in a magnetic environment, electric currents and electromagnetic forces are generated that affect velocity and pressure losses in the flow. These magnetohydrodynamic (MHD) interactions have to be investigated to optimize the engineering processes. The characteristics of MHD flows depend on the geometrical configuration, the strength of the applied magnetic field, the electrical properties of fluid and structural materials and the thermal conditions. In the so-called blankets for fusion reactors, where liquid metals are used to breed the plasma fuel component tritium and to extract the generated heat, magneto-convective flows play a crucial role in determining heat and mass transfer. Therefore, the availability of numerical codes to simulate this type of flow is mandatory and their validation is a necessary step to guarantee the reliability of the results. For that reason, a benchmark problem has been defined to simulate liquid metal flows in a horizontal rectangular duct heated from below and exposed to a non-uniform magnetic field. Results obtained by five research groups using different codes are compared.

benchmark↗