Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel”

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 OSTI2025

A two-fluid solar-wind model with intermittent Alfvénic turbulence

In one of the leading theories for the origin of the solar wind, photospheric motions launch Alfvén waves (AWs) that propagate along open magnetic-field lines through the solar atmosphere and into the solar wind. The radial variation in the Alfvén speed causes some of the AWs to reflect, and counter-propagating AWs subsequently interact to produce Alfveńic turbulence, in which AW energy cascades from long wavelengths to short wavelengths and dissipates, heating the plasma. In this paper we develop a one-dimensional two-fluid solar-wind model that includes Alfvénic turbulence, proton temperature anisotropy and a novel method for apportioning the turbulent heating rate between parallel proton heating, perpendicular proton heating and electron heating. We employ a turbulence model that accounts for recent observations from NASA’s Parker Solar Probe, which find that AW fluctuations in the near-Sun solar wind are intermittent and less anisotropic than in previous models of anisotropic magnetohydrodynamic turbulence. Our solar-wind model reproduces a wide range of remote observations of the corona and in-situ measurements of the solar wind, and our turbulent heating model consists of analytic equations that could be usefully incorporated into other solar-wind models and numerical models of more distant astrophysical plasmas.

astrophysical plasmas↗
DOE OSTI2025

Thermodynamics and collisionality in firehose-susceptible high- β plasmas

We study the evolution of collisionless plasmas that, due to their macroscopic evolution, are susceptible to the firehose instability, using both analytic theory and hybrid-kinetic particle-in-cell simulations. We establish that, depending on the relative magnitude of the plasma β, the characteristic time scale of macroscopic evolution and the ion-Larmor frequency, the saturation of the firehose instability in high-β plasmas can result in three qualitatively distinct thermodynamic (and electromagnetic) states. By contrast with the previously identified ‘ultra-high-beta’ and ‘Alfvén-inhibiting’ states, the newly identified ‘Alfvén-enabling’ state, which is realised when the macroscopic evolution time τ exceeds the ion-Larmor frequency by a β-dependent critical parameter, can support linear Alfvén waves and Alfvénic turbulence because the magnetic tension associated with the plasma’s macroscopic magnetic field is never completely negated by anisotropic pressure forces. We characterise these states in detail, including their saturated magnetic-energy spectra. The effective collision operator associated with the firehose fluctuations is also described; we find it to be well approximated in the Alfvén-enabling state by a simple quasi-linear pitch-angle scattering operator. The box-averaged collision frequency is ν eff ∼ β/τ, in agreement with previous results, but certain subpopulations of particles scatter at a much larger (or smaller) rate depending on their velocity in the direction parallel to the magnetic field. Our findings are essential for understanding low-collisionality astrophysical plasmas including the solar wind, the intracluster medium of galaxy clusters and black hole accretion flows. We show that all three of these plasmas are in the Alfvén-enabling regime of firehose saturation and discuss the implications of this result.

astrophysical plasmas↗
DOE OSTI2024

The development and applications of multidimensional biomolecular spectroscopy illustrated by photosynthetic light harvesting

The parallel and synergistic developments of atomic resolution structural information, new spectroscopic methods, their underpinning formalism, and the application of sophisticated theoretical methods have led to a step function change in our understanding of photosynthetic light harvesting, the process by which photosynthetic organisms collect solar energy and supply it to their reaction centers to initiate the chemistry of photosynthesis. The new spectroscopic methods, in particular multidimensional spectroscopies, have enabled a transition from recording rates of processes to focusing on mechanism. We discuss two ultrafast spectroscopies – two-dimensional electronic spectroscopy and two-dimensional electronic-vibrational spectroscopy – and illustrate their development through the lens of photosynthetic light harvesting. Both spectroscopies provide enhanced spectral resolution and, in different ways, reveal pathways of energy flow and coherent oscillations which relate to the quantum mechanical mixing of, for example, electronic excitations (excitons) and nuclear motions. The new types of information present in these spectra provoked the application of sophisticated quantum dynamical theories to describe the temporal evolution of the spectra and provide new questions for experimental investigation. While multidimensional spectroscopies have applications in many other areas of science, we feel that the investigation of photosynthetic light harvesting has had the largest influence on the development of spectroscopic and theoretical methods for the study of quantum dynamics in biology, hence the focus of this review. We conclude with key questions for the next decade of this review.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Crystal structure of perfluorononanoic acid, C 9 HF 17 O 2

The crystal structure of perfluorononanoic acid (PFNA) was solved via parallel tempering using synchrotron powder diffraction data obtained from the Brockhouse X-ray Diffraction and Scattering (BXDS) Wiggler Lower Energy (WLE) beamline at the Canadian Light Source. PFNA crystallizes in monoclinic space group P2 1 /c (#14) with lattice parameters a = 26.172(1) Å, b = 5.6345(2) Å, c = 10.9501(4) Å, and β = 98.752(2)°. The crystal structure is composed of dimers, with pairs of PFNA molecules connected by hydrogen bonds via the carboxylic acid functional groups. The Rietveld-refined structure was compared to a density functional theory-optimized structure, and the root-mean-square Cartesian difference was larger than normally observed for correct powder structures. The powder data likely exhibited evidence of disorder which was not successfully modeled.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Steric and Electronic Influence of Excited-State Decay in Cu(I) MLCT Chromophores

For the past eleven years, a dedicated effort in our research group focused on fundamentally advancing the photophysical properties of cuprous bis-phenanthroline-based metal-to-ligand charge transfer (MLCT) excited states. We rationalized that by gaining control over the numerous factors limiting the more widespread use of Cu I MLCT photosensitizers, they would be readily adopted in numerous light-activated applications given the earth-abundance of copper and the extensive library of 1,10-phenanthrolines developed over the last century. Significant progress has been achieved by recognizing valuable structure-property concepts developed by other researchers in tandem with detailed ultrafast and conventional time-scale investigations, in-silico-inspired molecular designs to predict spectroscopic properties, and applying novel synthetic methodologies. Ultimately, we achieved a plateau in exerting cooperative steric influence to control Cu I MLCT excited state decay. This led to combining sterics with π-conjugation and/or inductive electronic effects to further exert control over molecular photophysical properties. The lessons gleaned from our studies of homoleptic complexes were recently extended to heteroleptic bis(phenanthrolines) featuring enhanced visible light absorption properties and long-lived room-temperature photoluminescence. Furthermore, this Account navigates the reader through our intellectual journey of decision-making, molecular and experimental design, and data interpretation in parallel with appropriate background information related to the quantitative characterization of molecular photophysics using Cu I MLCT chromophores as prototypical examples.

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

Improved Characterization of Soil Organic Matter by Integrating FT-ICR MS, Liquid Chromatography Tandem Mass Spectrometry, and Molecular Networking: A Case Study of Root Litter Decay under Drought Conditions

Understanding of how soil organic matter (SOM) chemistry is altered in a changing climate has advanced considerably; however, most SOM components remain unidentified, impeding the ability to characterize a major fraction of organic matter and predict what types of molecules, and from which sources, will persist in soil. Here we present a novel approach to better characterize SOM extracts by integrating information from three types of analyses, and we deploy this method to characterize decaying root-detritus soil microcosms subjected to either drought or normal conditions. To observe broad differences in composition, we employed direct infusion Fourier-transform ion cyclotron resonance mass spectrometry (DI-FT-ICR MS). We complemented this with liquid chromatography tandem mass spectrometry (LC-MS/MS) to identify components by library matching. Since libraries contain only a small fraction of SOM components, we also used fragment spectral cosine similarity scores to relate unknowns and library matches through molecular networks. This integrated approach allowed us to corroborate DI-FT-ICR MS molecular formulas using library matches, which included fungal metabolites and related polyphenolic compounds. We also inferred structures of unknowns from molecular networks and improved LC-MS/MS annotation rates from ~5 to 35% by considering DI-FT-ICR MS molecular formula assignments. Under drought conditions, we found greater relative amounts of lignin-like vs condensed aromatic polyphenol formulas and lower average nominal oxidation state of carbon, suggesting reduced decomposition of SOM and/or microbes under stress. Our integrated approach provides a framework for enhanced annotation of SOM components that is more comprehensive than performing individual data analyses in parallel.

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

An Antibody Cocktail-Based Immunoaffinity-LC-MS Method Enabled Ultra-Sensitive and Robust Quantification of Circulating Proinsulin Proteoforms and C-Peptide

Accurately measuring circulating proinsulin proteoforms is crucial for clinical investigation of diabetes, but was previously not feasible owing to limited assay specificity/sensitivity. Here, in this study, we devised a highly sensitive LC-MS-based strategy to quantify intact proinsulin, des-31,32 and des-64,65 proinsulin, and C-peptide in circulation. The method involves: (i) quantitative, robust affinity capture using an optimized antibody cocktail, eliminating the severe quantitative bias across multiple proteoforms typically introduced when using a single antibody; (ii) Lys-C digestion producing unique signature peptides for each proteoform, and (iii) trapping-nano-LC coupled with FAIMS/dCV-MS for an ultrasensitive analysis. The selective trapping/delivery ensured sensitive/selective analysis of the targets while achieving excellent analytical robustness that is critical for clinical assays, and the FAIMS/dCV substantially reduces baseline noise/interferences, further enhancing S/N. The assay achieved exceptional sensitivity, with serum LOQs of 1.7, 2.3, and 3.6 pg/mL respectively for intact-proinsulin, des-31,32 and des-64,65, representing the first assay capable of sensitively quantifying these major circulating proinsulin proteoforms. We applied this assay to 78 subjects, including autoantibody positive (n = 20) and new-onset type 1 diabetes (T1D, n = 19) with respective age/sex/BMI-matched controls, enabling the first accurate profiling of proinsulin proteoforms in clinical groups. The assay results demonstrated a clear separation of control and new-onset T1D groups that a parallel total-proinsulin ELISA assay fails to capture. Furthermore, distinct expression patterns in relative abundance ratios among proteoforms were observed across clinical groups. This assay may provide valuable insights into the β-cell functions and the onset/progression of diabetes and other associated conditions. Moreover, the strategy is broadly applicable to targeted measurement of other biomarker proteoforms.

Shen, Qingqing [University at Buffalo, NY (United ↗
DOE OSTI2026

Integrating Flow Imaging Analysis and Single-Particle ICP-TOFMS for Comprehensive Micro- and Nanoplastic Characterization

Flow imaging analysis (FIA), provides composition-agnostic morphological characterization. These measurements of particle size and shape are valuable to mass-based analysis, such as single particle inductively coupled plasma time-of-flight mass spectrometry (sp-ICP-TOFMS), which provides quantitative data on elements within particles. Using these two methods together enables informed use of geometric assumptions required by sp-ICP-TOFMS, as particle mass is typically converted to a particle diameter using assumed-spherical geometry. To validate this concept, parallel measurements to determine particle diameters were performed by FIA and sp-ICP-TOFMS on four particle suspensions: 300 nm polystyrene Eu-doped nanoparticles, 1 μm Fe-rich beads, 3 μm four element calibration polystyrene beads and 5 μm polystyrene beads. The Fe-particles obtained the highest percent difference from the manufacturer’s nominal diameter, as the mean diameter obtained by FIA was overestimated by 21% and sp-ICP-TOFMS underestimated the mean diameter by 20.7%. Two types of particles were selected to test the effect of varying the particle number concentrations (PNC) on sizing accuracy, and both methods accurately sized each particle population at the PNC expected. Single particle analysis of carbon has continued to be a popular research topic, with direct applications to environmental pollutants in terms of nano- and micro- plastics. Real-world plastic particles were studied, and FIA’s measured circularity values demonstrated that the particles deviated from spherical geometries, therefore sp-ICP-TOFMS data should be interpreted as mass-based rather than size-based. Combining these techniques enables improved interpretation of particle populations and evaluation of particle sizes.

Szakas, Sarah [ORNL] (ORCID:0000000241332197)↗
DOE OSTI2024

Ectopic Production of 3,4-Dihydroxybenzoate in Planta Affects Cellulose Structure and Organization

Lignocellulosic biomass is a highly sustainable and largely carbon dioxide neutral feedstock for the production of biofuels and advanced biomaterials. Although thermochemical pretreatment is typically used to increase the efficiency of cell wall deconstruction, genetic engineering of the major plant cell wall polymers, especially lignin, has shown promise as an alternative approach to reduce biomass recalcitrance. Poplar trees with reduced lignin content and altered composition were previously developed by overexpressing bacterial 3-dehydroshikimate dehydratase (QsuB) enzyme to divert carbon flux from the shikimate pathway. In this work, three transgenic poplar lines with increasing QsuB expression levels and different lignin contents were studied using small-angle neutron scattering (SANS) and wide-angle X-ray scattering (WAXS). SANS showed that although the cellulose microfibril cross-sectional dimension remained unchanged, the ordered organization of the microfibrils progressively decreased with increased QsuB expression. This was correlated with decreasing total lignin content in the QsuB lines. WAXS showed that the crystallite dimensions of cellulose microfibrils transverse to the growth direction were not affected by the QsuB expression, but the crystallite dimensions parallel to the growth direction were decreased by ~20%. Cellulose crystallinity was also decreased with increased QsuB expression, which could be related to high levels of 3,4-dihydroxybenzoate, the product of QsuB expression, disrupting microfibril crystallization. In addition, the cellulose microfibril orientation angle showed a bimodal distribution at higher QsuB expression levels. Altogether, this study provides new structural insights into the impact of ectopic synthesis of small-molecule metabolites on cellulose organization and structure that can be used for future efforts aimed at reducing biomass recalcitrance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Introduction: Neuromorphic Materials

The explosive growth in data collection and the need to process it efficiently, as well as the desire to automate increasingly complex tasks in transportation, medical care, manufacturing, security and many other fields have motivated a growing interest in neuromorphic computing. Unlike the binary, transistorbased ON/OFF logic gates and separate logic and memory functionalities employed in digital computing, neuromorphic computing is inspired by animal brains that use interconnected synapses and neurons to perform processing, storage and transmission of information at the same location, while only consuming ~20 W or less of power. Motivated by the brain’s efficiency, adaptability, self-learning and resiliency qualities, neuromorphic computing can be broadly defined as an approach to processing and storing information using hardware and algorithms inspired by models of biological neural systems. Present research in neuromorphic computing encompasses approaches that vary significantly in their degree of neuro-inspiration, from systems that only incorporate features such as asynchronous, event-driven operation or use crossbar arrays of non-volatile memory (NVM) elements to accelerate deep neural networks (DNNs), to designs that embrace the extreme parallelism, sparsity, reconfigurability, adaptability, complexity and stochasticity observed in nervous systems. The term ‘neuromorphic’ computing is often credited to Carver Mead, who in the 1980s investigated Si-based analog electronics to replicate functions of the animal retina. Earlier important advances in this field include the work of Frank Rosenblatt, who proposed the concept of the perceptron, Bernard Widrow, who used this concept to build one of the first analog neural networks, the Adaline and many other researchers (see ref. 6 for an historical perspective on neuromorphic computing). With the recent increase in the use of artificial intelligence and large language models, and rising concerns over the associated energy costs, interest in neuromorphic hardware has expanded rapidly. According to some estimates, driven largely by the drastic growth in the training use of artificial intelligence (AI) models using the current computing architectures, the energy cost of computing is projected to reach the energy supply worldwide by 2045. Furthermore, while this is not a realistic outcome, it means that, if more efficient computing technologies are not developed -- soon -- the world will soon become one where demand for energy and market constraints limit the continued increase of societal access to AI and cloud services from data centers. Data centers used for training and use of these models consume hundreds of terawatt hours of electricity, already past 4% of the US electricity demand.

Circuits↗
DOE OSTI2026

Function, Structure, and Regulation of Nitrogen Fixation-like Metalloproteins for Nitrogen, Energy, Carbon, and Sulfur Metabolism

Nitrogenases (N 2 ases) and nitrogen fixation-like (NFL) systems play distinct roles in nitrogen, carbon, sulfur, and energy metabolism based on their fundamental differences in structure and metallocofactor identity. As new NFL systems have recently been identified and characterized, striking parallels and differences compared to N 2 ase structure, catalysis, and regulation have emerged. NFL systems use metallocofactors that span from simple [4Fe-4S] clusters to complex clusters akin to FeMo-co, previously only thought to occur in N 2 ase. This review describes the present state of knowledge on the function, structure, catalytic mechanisms, and regulation of NFL systems that perform distinct biological roles across all three domains of life. Recent advancements in N 2 ase spectroscopic techniques for probing metallocofactor structure and electronic states guide current and future work on how each NFL system catalyzes its specific biological reaction(s). Key knowledge gaps and needed areas of research for uncovering the specific metallocofactors and structural motifs that are at the heart of NFL system reaction specificity, along with how these systems are regulated, are discussed.

Bacteria↗
DOE OSTI2025

Effect of Impurities on the Redox Properties of Goethite

Iron oxide minerals regulate the flux of electrons in the environment and are important hosts for trace and minor, yet critical, elements. Here, we present the first evidence of a direct link between the local coordination environments of Ni and Zn and the redox properties of their host phase goethite (α-FeOOH), the most abundant Fe(III) (oxyhydr)oxide at Earth’s surface. Here, we used aqueous redox measurements to show that the redox potential E H 0 , and hence the mineral’s stability, follows the order: pure goethite ≥ Zn-goethite > Ni-goethite. Parallel X-ray absorption and scattering measurements demonstrate, using quantum-informed analysis, that the local coordination environment of the smaller impurity, Ni, causes more bulk strain energy than Zn, which nearly accounts for the difference in E H 0 between Ni- and Zn-goethite. Our theory-informed, experimental study reveals how two common impurities affect the stability of goethite with implications for the biogeochemical reactivity of Fe(III) (oxyhydr)oxide in mediating elemental and electron fluxes in the environment.

42 ENGINEERING↗
DOE OSTI2025

Effects of Particle Mixing and Gravitational Settling on Charge Transport in Carbon Flow-Electrode Cells

Flowable carbon slurries are actively studied and under development for charge transport in various electrochemical systems including flow capacitors, capacitive deionization cells, semi-solid flow batteries, and lithium extraction. However, much less is known about in operando slurry flow dynamics and their corresponding effect on charge transport. We performed an experimental study of mixing and settling dynamics of slurry electrodes within an electrochemical flow cell during continuous operations. The electrochemical cell consisted of two horizontal co-flowing channels, separated by a cation-exchange membrane (CEM). Here we used high-speed optical imaging of planes parallel to gravity and simultaneous electrochemical measurements. At low flow rates, dense yet dynamic particle beds formed on the bottom electrode in each channel, which unexpectedly yielded the highest currents. This approach enables the operation of the flow cell at low system-average particle concentrations while leveraging gravity-driven particle settling to locally enhance carbon concentrations precisely at the current collector sites. Conversely, high flow rates were characterized by thin particle beds and well-mixed particle flows. In the latter case, the electrodes in closest proximity (located on either side of the CEM) achieved a current higher than the other electrode pairs. The observations have implications for slurry control and electrode designs in electrochemical systems.

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

Global Impacts of Marine Methanethiol Emissions and Chemistry in the Atmosphere

Oceanic emissions of dimethyl sulfide (DMS) have long been known to influence aerosol particle composition, cloud condensation nuclei (CCN) concentration, and Earth’s radiative budget. However, the impact of oceanic emissions of methanethiol (MeSH), a sulfur compound produced by the same oceanic precursor as DMS, has been relatively less explored. The gas-phase oxidation of MeSH has a higher effective yield of SO 2 and a shorter oxidative lifetime compared to DMS, highlighting the relevance of this pathway for the modeled representation of particle formation, growth, and CCN abundance in the marine atmosphere. Here, we use the global chemical transport model GEOS-Chem to explore possible scenarios representative of specific environmental conditions and MeSH emission schemes based on previous experimental studies. We further implement and test previously reported chemical mechanisms for MeSH oxidation, along with additional improvements, highlighting key uncertainties and sensitivities for regional and global sulfur budgets. We place our results in the context of recent modeling updates to DMS chemistry and cloud processing, which further impact SO 2 production in the marine atmosphere in parallel with MeSH oxidation. Within the overall marine sulfur budget, our findings highlight that MeSH plays a significant role in SO 2 production in the marine atmosphere, contributing to regional surface layer concentration increases of up to 40–60%. These results point to the importance of MeSH for efforts aimed at improving the modeled representation of sulfur spatiotemporal patterns relevant to air quality predictions and climate impact assessments.

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

Development of Steady-State and Dynamic Mass and Energy Constrained Neural Networks for Distributed Chemical Systems Using Noisy Transient Data

The paper presents the development of algorithms for mass and energy constrained neural network models that can exactly conserve the overall mass and energy of distributed chemical process systems, even though the noisy transient data used for optimal model training violate the same. In contrast to approximately satisfying mass and energy balance constraints of a system by soft penalization of objective function, algorithms have been developed for solving equality-constrained nonlinear optimization problems, thus providing the guarantee of exactly satisfying the system mass and energy conservation laws. For developing dynamic mass-energy constrained network models for distributed systems, hybrid series and parallel dynamic-static neural networks have been leveraged. The developed algorithms for solving both the training and forward problems are validated using both steady-state and dynamic data in the presence of various noise characteristics. The developed data-driven algorithms are flexible to exactly satisfy mass and energy balance constraints for dynamic chemical processes if the system holdup information is available. The proposed network structures and algorithms are applied to the development of data-driven lumped and distributed models of an adiabatic superheater/reheater system, a nonisothermal continuous stirred tank reactor, as well as an electrically heated plug-flow reactor system where one form of energy gets transformed to another. It has been observed that the mass-energy constrained neural networks yield a root mean squared error of <1% with respect to the system truth for the case studies evaluated in this work.

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

Magnetic and EPR Spectroscopic Studies of Thiolate Bridged Divalent Ni, Pd, and Pt Ions Capped with VO(N 2 S 2 ) Metalloligands

Reactions of the metallodithiolate complex VO(bme-dach) (hereafter abbreviated as V, where bme-dach = N,N ′-bis(2-mercaptoethyl)-1,4-diazacycloheptane) with [Pd II (CH 3 CN) 4 ](BF 4 ) 2 and [Pt II (CH 3 CN) 4 ](BF 4 ) 2 yield the V–M–V trimetallic compounds [VPdV](BF 4 ) 2 (2) and [VPtV](BF 4 ) 2 (3). Reaction of a similar metalloligand, VO(bme-daco) (hereafter abbreviated as V′ where bme-daco = N,N′-bis(2-mercaptoethyl)-1,5-diazacyclooctane) with [Ni II (CH 3 CN) 6 ](BF 4 ) 2 afforded the related salt [V′NiV′](BF 4 ) 2 (1). X-ray structural analyses revealed that cations in 1, 2, and 3 adopt a stairstep C 2h structure consisting of two terminal VO(N 2 S 2 ) moieties bridged via thiolate sulfur to the group 10 metal ions. Weak ferromagnetic superexchange coupling (J = 0.282 cm –1 for 1, 0.954 cm –1 for 2, and 1.372 cm –1 for 3) was observed between the two S = 1/2 V IV centers separated by distances in the range of 5.9–6.3 Å, with J values varying following the order Ni < Pd < Pt. Frozen-solution EPR spectra measured on the more soluble [BArF 24 ] − (BArF 24 – = tetrakis((3,5-trifluoromethyl)phenyl)borate) analogues revealed that the [VPtV] 2+ cation exhibits a 15-line hyperfine splitting of 225 MHz at g = 4 in parallel mode, confirming exchange coupling between the two 51 V, I = 7/2 nuclei. Density-functional theory (DFT) calculations indicate an S = 1 ground state for 1–3. These results demonstrate that the choice of paramagnetic metallodithiolate ligand and diamagnetic bridge in such trimetallic species influences the sign and magnitude of magnetic interactions.

anions↗
DOE OSTI2025

Simulation of 24,000 Electron Dynamics: Real-Time Time-Dependent Density Functional Theory (TDDFT) with the Real-Space Multigrids (RMG)

Here, we present the theory, implementation, and benchmarking of a real-time time-dependent density functional theory (RT-TDDFT) module within the RMG code, designed to simulate the electronic response of molecular systems to external perturbations. Our method offers insights into nonequilibrium dynamics and excited states across a diverse range of systems, from small organic molecules to large metallic nanoparticles. Benchmarking results demonstrate excellent agreement with established TDDFT implementations and showcase the superior stability of our time integration algorithm, enabling long-term simulations with minimal energy drift. The scalability and efficiency of RMG on massively parallel architectures allow for simulations of complex systems, such as plasmonic nanoparticles with thousands of atoms. Future extensions, including nuclear and spin dynamics, will broaden the applicability of this RT-TDDFT implementation, providing a powerful toolset for studies of photoactive materials, nanoscale devices, and other systems where real-time electronic dynamics is essential.

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

GPU-Accelerated Solution of the Bethe–Salpeter Equation for Large and Heterogeneous Systems

We present a massively parallel GPU-accelerated implementation of the Bethe–Salpeter equation (BSE) for the calculation of the vertical excitation energies (VEEs) and optical absorption spectra of condensed and molecular systems, starting from single-particle eigenvalues and eigenvectors obtained with density functional theory. The algorithms adopted here circumvent the slowly converging sums over empty and occupied states and the inversion of large dielectric matrices through a density matrix perturbation theory approach and a low-rank decomposition of the screened Coulomb interaction, respectively. Further computational savings are achieved by exploiting the nearsightedness of the density matrix of semiconductors and insulators to reduce the number of screened Coulomb integrals. We scale our calculations to thousands of GPUs with a hierarchical loop and data distribution strategy. The efficacy of our method is demonstrated by computing the VEEs of several spin defects in wide-band-gap materials, showing that supercells with up to 1000 atoms are necessary to obtain converged results. We discuss the validity of the common approximation that solves the BSE with truncated sums over empty and occupied states. In conclusion, we then apply our GW-BSE implementation to a diamond lattice with 1727 atoms to study the symmetry breaking of triplet states caused by the interaction of a point defect with an extended line defect.

Absorption spectra↗