Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical mechanics”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Towards satellite tests combining general relativity and quantum mechanics through quantum optical interferometry: progress on the deep space quantum link

The Deep Space Quantum Link (DSQL) is a space-mission concept that aims to explore the interplay between general relativity and quantum mechanics using quantum optical interferometry. This mission concept was formally presented to the United States National Academy of Science Decadal Survey as a research campaign for Fundamental Physics in 2022. Since then, advances have been made in the space-based quantum optical technologies required to conduct a DSQL-type mission. In addition, other research efforts have defined alternative measurement concepts to explore the same scientific questions motivating the DSQL mission. This paper serves as an update to the community on the status of the DSQL mission concept and related research and technology development efforts.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

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

Impact of Laboratory-Accelerated Aging Methods to Study Alkali–Silica Reaction and Reinforcement Corrosion on the Properties of Concrete

This study focuses on two separate investigations of the main aging mechanisms: alkali–silica reactivity (ASR) and the corrosion of reinforcing steel (rebar) concrete, both of which may result in a premature failure to meet the serviceability or strength requirements of a concrete structure. However, these processes occur very slowly, spanning decades. The impact of direct chemical additives to fresh concrete to accelerate ASR and the corrosion of reinforcing steel on the fresh and hardened properties of the ensuing material are investigated to inform the potential use of chemicals in large-scale studies. The deterioration of reinforced concrete (RC) is determined by means of expansion, cracking, bulk diffusivity and surface resistivity measurements, and compressive, split tensile and flexural strength tests. The results indicate that the addition of sodium hydroxide and calcium chloride can effectively accelerate the crack formation and propagation in concrete due to ASR and the corrosion of rebar, respectively. The ASR-induced cracks maintained a constant crack width from 0.05 mm to 0.1 mm over the measurement period regardless of the intensity of aging acceleration. Adding 4% chloride by weight of cement for accelerating rebar corrosion resulted in an average crack that was 82% larger than in the case of ASR accelerated with the addition of sodium hydroxide. The addition of alkali resulted in an increase in early-age (7-day) strength. At a total alkali loading of 2.98 kg/m3, 3.84 kg/m3 and 5.57 kg/m3, the 28-day compressive strength of concrete decreased by 3%, 10% and 24%, respectively. Similarly, a higher early-age strength and a lower later-age strength was observed for the concrete in the presence of corrosive calcium chloride. The results from this research are expected to inform future studies on the long-term performance of RC structures under accelerated ASR and corrosion.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Single-Step Conversion of Ethanol to n-Butene over Ag-ZrO2/SiO2 Catalysts

Ethanol is a promising platform molecule for production of a variety of fuels and chemicals. Of particular interest is producing middle distillate fuels (i.e., jet and diesel blendstock) from renewable ethanol feedstock. State-of-the-art alcohol-to-jet technology requires multiple process steps based on catalytic dehydration of ethanol to form ethylene, followed by sometimes a multi-step oligomerization, and then hydrotreatment and distillation. Here we report on a new catalytic route in which ethanol is directly converted to n-butene (1- and 2-butene mixtures) over Ag-ZrO2/SBA-16, thus offering the potential for a reduction in the number of required processing steps versus conventional alcohol-to-jet technology. This catalyst system provides the balanced metal and Lewis acid sites required to selectively facilitate a cascading sequence of reactions that includes dehydrogenation, aldol condensation, Meerwein–Ponndorf–Verley reduction, dehydration, and hydrogenation. High conversion and selectivity toward either n-butene or 1,3-butadiene is achieved by tuning the hydrogen feed partial pressure and other process/catalyst parameters. With sufficient hydrogen partial pressure 1,3-butadiene is completely and selectively hydrogenated to form n-butene. The reaction mechanism was elucidated through operando-nuclear magnetic resonance investigations coupled with reactivity measurements. Combined experimental-computational investigation reveals how changes in silver and zirconium composition and the silver oxidation state affects reactivity under controlled hydrogen partial pressures and after prolonged run times. Finally, catalyst effectiveness also was demonstrated when using wet ethanol feed, thus highlighting process flexibility in terms of feedstock purity requirements. This work was financially supported by the U.S. Department of Energy (DOE), Office of Energy Efficiency and Renewable Energy, Bioenergy Technologies Office, and was performed at the Pacific Northwest National Laboratory (PNNL) under Contract No. DE-AC05-76RL01830 and the National Renewable Energy Laboratory under Contract No. DE-AC36- 08GO28308. Part of the work conducted by S. A. Akhade was performed under the auspices of the U.S. DOE at Lawrence Livermore National Laboratory under Contract No. DE-AC52-07NA27344. This work was partly supported through the PNNL-WSU Distinguished Graduate Research Program for ADW. NMR and XPS experiments were performed using EMSL (grid.436923.9), a DOE Office of Science User Facility sponsored by the Office of Biological and Environmental Research

Dagle, Vanessa↗
DOE OSTI2021

Broadband 2DES detection of vibrational coherence in the S x state of canthaxanthin

The nonadiabatic mechanism that mediates nonradiative decay of the bright S 2 state to the dark S 1 state of carotenoids involves population of a bridging intermediate state, S x , in several examples. The nature of S x remains to be determined definitively, but it has been recently suggested that S x corresponds to conformationally distorted molecules evolving along out-of-plane coordinates of the isoprenoid backbone near a low barrier between planar and distorted conformations on the S 2 potential surface. Here, the electronic and vibrational dynamics accompanying the formation of S x in toluene solutions of the ketocarotenoid canthaxanthin (CAN) are characterized with broadband two-dimensional electronic spectroscopy (2DES) with 7.8 fs excitation pulses and detection of the linear polarization components of the third-order nonlinear optical signal. A stimulated-emission cross peak in the 2DES spectrum accompanies the formation of S x in <20 fs following excitation of the main absorption band. S x is prepared instantaneously, however, with excitation of hot-band transitions associated with distorted conformations of CAN’s isoprenoid backbone in the low frequency onset of the main absorption band. Vibrational coherence oscillation maps and modulated anisotropy transients show that S x undergoes displacements from the Franck–Condon S 2 state along out-of-plane coordinates as it passes to the S 1 state. The results are consistent with the conclusion that CAN’s carbonyl-substituted β-ionone rings impart an intramolecular charge-transfer character that frictionally slows the passage from S x to S 1 compared to carotenoids lacking carbonyl substitution. Despite the longer lifetime, the S 1 state of CAN is formed with retention of vibrational coherence after passing through a conical intersection seam with the S x state.

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

Projector-Free Capped-Fragment Scheme within Density Functional Embedding Theory for Covalent and Ionic Compounds

Quantum-mechanics-(QM)-based simulations now routinely aid in understanding and even discovering new chemistries involving molecules and materials exhibiting desired functionalities. Ab initio correlated wavefunction (CW) theories systematically improve QM methods, with many exhibiting high accuracy. However, execution of CW methods requires expensive computations that typically scale poorly with system size. Divide-and-conquer approaches partition large systems into smaller fragments; a lower level of theory treats fragment interactions while a preferred higher level of theory describes the important fragment. These methods offer ways to incorporate CWs into chemical simulations of large systems, e.g., biomolecules, surfaces, large inorganic clusters, bulk crystals, etc. Here we propose a partitioning protocol that utilizes capping atoms to saturate severed covalent bonds at fragment interfaces and density functional embedding theory (DFET) to describe fragment interactions. The capping groups in each fragment provide an ad-hoc potential that approximates the effects of the environment. An embedding potential optimized via DFET then serves as an augmentation of the capping group to simulate the effects of the environment. We concurrently use an auxiliary fragment (a separate system comprised of only the combined capping groups) to account for, and thereby correct, the electron density contributions of all the capping groups added to all of the fragments. This method depends only on the capped-subsystem and auxiliary-fragment electron densities, forgoing, as with the original DFET developed for metallic systems, orbital-based projector approaches that determine a nonlocal action of the embedding potential onto the fragment electrons. By using an auxiliary fragment, the method maintains a purely electron-density-dependent embedding potential, substantially lessening the cost and leading to simpler implementation. In this work, we demonstrate the utility of our capped-DFET and ensuing capped embedded CW method in two contrasting systems, namely, an organic molecule and an ionic metal oxide cluster.

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

Losses of CO and CO 2 upon collision-activated dissociation of substituted 2-methoxyphenoxides after methyl radical loss

Some deprotonated, substituted 2-methoxyphenols fragment by methyl radical loss followed by competing losses of CO and CO 2 upon collisional activation in linear quadrupole ion trap mass spectrometers. These reactions were examined experimentally and computationally in order to determine their mechanisms. For deprotonated vanillin, the CO loss was found to involve ring contraction, with the highest free energy barrier of 58.0 kcal mol -1 (for the syn rotamer; here, syn refers to the relationship between the aldehyde oxygen atom and the oxygen atom of the methoxy group) for the entire process. The atoms lost in this fragmentation are the oxygen atom that was bound to the first eliminated methyl group and the aromatic carbon atom bound to this oxygen. Examination of carbon-13 labeled vanillin supports these assignments. Examination of several model compounds revealed that this reaction requires the presence of an electron-withdrawing substituent in the para-position relative to the phenol moiety. In contrast, the CO 2 loss from deprotonated vanillin occurs via ring opening followed by re-cyclization and then ring contraction, leading to the loss of CO 2 in a process wherein the highest free energy barrier is 84.5 kcal mol -1 (for the syn isomer). The atoms lost in this fragmentation are the carbon and oxygen atoms from the phenoxide group and the oxygen atom that was bound to the first eliminated methyl group, which is supported by examination of carbon-13 labeled vanillin. Despite the higher total free energy requirement, the CO 2 loss is competitive with CO loss, possibly due to favorable entropy and low activation energy (free energy barrier of 48.9 kcal mol -1 ) for the first reaction step (the analogous value for CO loss is 58.0 kcal mol -1 ). The extent of CO 2 loss is strongly affected by substituents – it either is not observed or is very slow for the other compounds studied here. For example, it is substantially less favorable than CO loss (highest free energy barrier 60.9 kcal mol -1 ) for deprotonated acetovanillone for which the first reaction step for CO 2 loss has a substantially greater free energy barrier (60.6 kcal mol-1) than for vanillin.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Illuminating the Photophysics of Oxygen Atom Transfer in Rare Earth Metal-Organic Complexes

Industrial separations of rare earth metals are reliant on selectivity based on ionic radii, which are time, energy, and resource consuming.1 The unique and discontinuous manifold of energy states associated with the f-electron configurations of the rare earth ions2 offer an opportunity to overcome some of these limitations by incorporating light into the separation process. Here, we use transient absorption spectroscopy to understand the differences in mechanism of light-initiated oxygen atom transfer (OAT) in yttrium (Y) and dysprosium (Dy) rare earth complexes. These metal-organic complexes are functionalized with hexfluoroacetyl-acetonate (hfac) and 4-methylmorpholine-N-oxide (NMMO). After light absorption in the presence of triphenyl phosphine (TPP), the NMMO in the Y-NMMO complex is replaced by the TPPO OAT reaction product. Because the Dy f-electron manifold is accessible to the ligand excited states, while there are no accessible states associated with f-electrons for Y, there are significant differences in reactivity. Transient absorption spectroscopy (TAS) was used to monitor changes after 340 nm photoexcitation of Dy/Y-NMMO or the product complex Dy/Y-TPPO. TAS is a pump-probe technique where a white light probe pulse is delayed relative to a visible pump pulse, which photoexcites a sample. Differences between the ground and excited states are monitored over time. With TAS, we access time scales from 100 fs-400 us to determine the photophysics of the complexes of interest. Photoexcitation of Y-NMMO (Fig. 1A) creates a broad photoinduced absorption (PIA) that evolves into a sharp positive feature over 4.2 ps centered at 450 nm. A secondary feature centered at 550 nm grows in with a 200 ps time constant and remains through the 5 ns time window of our experiment. Control experiments for Y-TPPO reveal that the feature at 450 nm is due to internal dynamics of the hfac ligand, while the 550 nm shoulder is due to the presence of the NMMO. Comparing kinetics at 550 nm for Y-NMMO and Dy-NMMO (Fig. 1B) shows that while Y exhibits a distinct growth in the NMMO feature, the Dy decays almost completely within 5 ns. Control experiments with Dy-TPPO prove that the ligand excited states centered on the hfac ligand are deactivated via energy transfer to the Dy manifold within 200 ps. This deactivation leads to significant differences in OAT reactivity between Y and Dy.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗
DOE OSTI2020

Benchmarking an Embedded Adaptive Sampling Configuration Interaction Method for Surface Reactions: H 2 Desorption from and CH 4 Dissociation on Cu(111)

Embedded (emb-) correlated wavefunction (CW) theory enables accurate assessments of both ground- and excited-state reaction mechanisms involved in heterogeneous catalysis. Embedded multireference second-order perturbation theory (emb-MRPT2) based on reference wavefunctions generated via embedded complete active space self-consistent field (emb-CASSCF) theory is currently state-of-the-art. However, the factorial scaling of CASSCF limits the size of active space and the complexity of systems that can be studied. In this work, we assess the efficacy of an alternative CW method, adaptive sampling configuration interaction (ASCI)–which enables large active spaces to be used–for studying surface reactions. We couple ASCI with density functional embedding theory (DFET) and benchmark its performance for two reactions: H 2 desorption from and CH 4 dissociation on the Cu(111) surface. Unlike embedded complete active space second-order perturbation theory (emb-CASPT2) that accurately reproduces a measured H 2 desorption barrier, embedded ASCI, using a very large active space (though one that still comprises a small portion of the full set of orbitals) fails to do so. Adding an extra correlation term from embedded Møller–Plesset second-order perturbation theory (emb-MP2) improves the desorption barrier and endothermicity predictions. Thus, the inaccuracy of embedded ASCI comes from the missing dynamic correlation from the many other electrons and orbitals not included in the active space. For CH 4 dissociation, again embedded ASCI overestimates the dissociation barrier compared to emb-CASPT2 predictions. Adding dynamic correlation from emb-MP2 helps correct the barrier. However, this composite approach suffers from double counting of correlation within embedded ASCI followed by emb-MP2 calculations. We therefore conclude that the state-of-the-art emb-MRPT2 based on reference wavefunctions generated via emb-CASSCF remains the method of choice for studying surface reactions. emb-ASCI is useful when large active spaces beyond the limit of emb-CASSCF are essential, such as to study complex surface reactions with significant multiconfigurational character (static correlation) but weak dynamic correlation.

08 HYDROGEN↗
DOE OSTI2020

ChIMES: A Machine-Learned Interatomic Model Targeting Improved Description of Condensed Phase Chemistry in Energetic Materials

In this report we detail completion of a Physics and Engineering Model Level Two Milestone targeting improved reactive interatomic potentials (IAPs) for energetic materials (EM) through machine learning. The specific goals of this milestone were to develop, validate, and document a new reactive molecular dynamics method for EM, based on machine learning by (1) generating databases of first-principles-derived forces, stresses, and energies for HN3 and 3,4-bis(3-nitrofurazan- 4-yl)furoxan (DNTF) (2) generate atomistic force fields from these databases via ML, and (3) benchmark model performance against first principles calculations. These goals were achieved by (1) further developing a machine learned reactive IAP and generation approach (i.e. the Chebyshev Interaction Model for Efficient Simulation or “ChIMES”), for which resulting IAPs can approach the predictive power of quantum-mechanical approaches at a fraction of the computational expense, and (2) applying the ChIMES framework to develop models for HN3 and DNTF. We find that for simple energetic materials like HN3, high accuracy ChIMES models can be obtained through application of a fitting approach that does not use active machine learning. We demonstrate the suitability of ChIMES models for simulations involving EM by using the HN3 model in multiscale shock technique simulations to predict the HN3 Chapman-Jouguet detonation state and investigate chemical evolution out to 1 ns following shock compression. This model is then used in larger direct shock (DS) simulations for a preliminary investigation of how bubbles (i.e. voids) influence material response under shock compression. We find that more complex EM (i.e. DNTF) necessitate a more sophisticated fitting approach, and develop a new active learning method and python tool to meet this challenge. We demonstrate that this fitting approach yields ChIMES models that out-perform commonly used standard reactive IAPs as well as semi-empirical quantum methods, and discuss the systematic improvability of these actively learned ChIMES models. We also describe challenges related to model development for EM such as DNTF, for which few experimental or previous simulation data are available (e.g. which could otherwise inform generation of training data). To overcome this issue, we establish a semi-empirical quantum ChIMES capability which can be used to efficiently map out relevant thermodynamic and configurational space, and generate ChIMES-IAP training data in a multiscale manner. We also show that these semi-empirical quantum ChIMES models can be used to generate predictions for the shock Hugoniot (the Hugoniot is the locus of thermodynamic states found in a shocked material) equation of state, investigate related thermochemistry, and explore carbon condensation following shock compression. This work represents a substantial advance in our atomistic modeling capability for EM that will provide much needed information on the chemistry of detonation for continued development of continuum models based on the Cheetah thermochemical code.

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

Metal-to-Ligand Charge-Transfer Spectrum of a Ru-Bipyridine-Sensitized TiO 2 Cluster from Embedded Multiconfigurational Excited-State Theory

Understanding optical properties of the dye molecule in dye-sensitized solar cells (DSSCs) from first-principles quantum mechanics can contribute to improving the efficiency of such devices. While density functional theory (DFT) and time-dependent DFT have been pivotal in simulating optoelectronic properties of photoanodes used in DSSCs at the atomic scale, questions remain regarding DFT’s adequacy and accuracy to furnish critical information needed to understand the various excited-state processes involved. In this work, we simulate the absorption spectra of a dye-sensitized solar cell analogue, comprised of a Ru-bipyridine (Ru-bpy) dye molecule and a small TiO2 cluster via DFT and via an accurate embedded correlated wavefunction (CW) theory. We generated CW spectra for the adsorbed Ru-bpy dye via a recently introduced capped density functional embedding theory or capped-DFET (to generate the embedding potential that accounts for the interaction of the molecule and the TiO 2 cluster). We then combined capped-DFET with the accurate but expensive multiconfigurational complete active space second-order perturbation theory (CASPT2)—embedded CASPT2. Because the CW theory is conducted on only a portion of the total system in the presence of an embedding potential that describes that portion’s interaction with its environment, we efficiently obtain CW-quality predictions that reflect local properties of the entire system. Specifically, for example, with capped-DFET and embedded CW theory, we can simulate accurately a plethora of metal-to-ligand charge-transfer excited properties at a manageable computational cost. Here, we predict detailed electronic spectra within the visible region, featuring the lowest three singlet and triplet excited states, along with predictions of the singlets’ lifetimes. We illustrated these results using a Jablonski diagram that show the relative energy position of the singlet and longer-lived triplet excited states and analyzed and proposed relaxation paths for the excited state corresponding to the most intense but short-lived absorption (interconversion, intersystem crossing, fluorescence, and phosphorescence) that may lead to longer-lived excited states necessary for efficient charge separation required to generate current in solar cells.

14 SOLAR ENERGY↗
DOE OSTI2024

Enhanced Production of Λ b 0 Baryons in High-Multiplicity p p Collisions at s = 13 TeV

The production rate of Λ b 0 baryons relative to B 0 mesons in p p collisions at a center-of-mass energy s = 13 TeV is measured by the LHCb experiment. The ratio of Λ b 0 to B 0 production cross sections shows a significant dependence on both the transverse momentum and the measured charged-particle multiplicity. At low multiplicity, the ratio measured at LHCb is consistent with the value measured in e + e − collisions, and increases by a factor of ∼ 2 with increasing multiplicity. At relatively low transverse momentum, the ratio of Λ b 0 to B 0 cross sections is higher than what is measured in e + e − collisions, but converges with the e + e − ratio as the momentum increases. These results imply that the evolution of heavy b quarks into final-state hadrons is influenced by the density of the hadronic environment produced in the collision. Comparisons with several models and implications for the mechanisms enforcing quark confinement are discussed. © 2024 CERN, for the LHCb Collaboration 2024 CERN

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Recent Developments in DFTB+, a Software Package for Efficient Atomistic Quantum Mechanical Simulations

DFTB+ is a flexible, open-source software package developed by its community, designed for fast and efficient atomistic quantum mechanical simulations. It employs various methods that approximate density functional theory (DFT), such as density functional-based tight binding (DFTB) and the extended tight binding (xTB) approach allowing simulations of large systems over extended time scales with reasonable accuracy, while being significantly faster than traditional ab initio methods. In recent years, several new extensions of the DFTB method have been developed and implemented in the DFTB+ program package in order to improve the accuracy and generality of the available simulation results. In this paper, we review those enhancements, show several use case examples and discuss the strengths and limitations of its features.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dynamic, Adaptive, Systems and Materials: Complex, Simple and Emergent Behaviors

This program has been funded by DoE/BES for twenty years. It has moved into and out of various subjects as it has developed, but it has retained its focus on complexity and complex systems. The project has evolved in the following way: Self-Assembly and Biomimetic Self-Assembly: All self-assembling systems depend upon a minimum of two types of interaction: a repulsion and an attraction. For the familiar molecular systems, attractive interactions are typically hydrogen bonds and electrostatic interactions. Repulsive interactions include steric effects, hydrophobic effects (in biological systems), and charge-charge repulsion. We have expanded this repertoire to include surface interactions, magnetic interactions, and others. I list these systems in the order in which we have explored them: i) A key emphasis in current work is in understanding how the movement of ions in a magnetic field (the Lorentz effect) interacts with catalytic systems. We have demonstrated that an acceleration in rate of reduction of CO 2 to CO can be accomplished by applying an external magnetic field. This acceleration is largely due to the application of the Lorentz effect on mass transport at the catalyst’s surface. ii) We have also extensively explored the influence of electrostatics, as exhibited in self-assembling systems, by tribocharging. iii) Another key system involves surface tension effects; examples include interactions between heavy particles floating at a liquid-air interface, and interacting by changes in surface area; interactions of bubbles and bubble rafts, behaviors of bubble trains in microfluidic networks, and behaviors of microorganisms in constraining environments. iv) This work has intentionally de-emphasized biological systems; but it does include some work on protein-ligand interactions and interactions among microorganisms. v) We have also explored applications of some of these effects, these explorations include bubble rafts as diffraction gratings, exploration of the structures that can be obtained by tribocharging and uses of these structures in exploring nucleation and melting of crystals. vi) Although not a major focus of this work, several other topics have emerged and offer opportunities for future work. These include the behavior of bubble trains and bubble rafts in microfluidic systems. A particularly interesting example is the formation of bubble trains that repeat in the alteration of large and small bubbles according to rules we do not presently understand, but are uniquely large-period oscillating systems. These systems offer a new route into understanding the instabilities of the type represented by oscillations. vii) We have also begun exploratory projects on magnetic levitation (especially to determine molecular density), and information storage (in molecules). Magnetic Levitation: Self-assembling and biomimetic systems require both attraction and repulsion. We have used electrostatics (tribocharging), interfacial free-energies (surface tension and related forces) and others. Potential uses include reconfigurable diffraction gratings and liquid lenses; exploration of mechanisms in tribocharging; tunneling in EGaIn junctions; and bubble trains (especially in micro-fluidic systems). Examples of systems representing these topics is included in the following papers: Complexity: Disks rotating at a water-air interface; Benard-Marangoni effects; Vortex-Crystals from spinning magnetic disks (Marangoni effects); EGaIn Electrode to study quantum tunneling; Self-Assembly of electrostatically-charged metallic spheres (electrets); Dynamically reconfigurable lens; Using computational designs of ligands for enzymes; Electrostatic self-assembly by tribocharging; Monodisperse bubble trains in microchannel systems; Inverted dripping faucet; Flames; Printing of micro-organisms to regenerate the “ink” of printing device; Using micro-organisms to move loads (“microoxen”); Motion of bacterial swarms near surfaces; Making monodisperse particles in microfluidic systems; Coding/decoding of information stored in droplet trains in microfluidic networks; Magnetic levitation; and Information storage. i) Tribocharging. The change in focus of this work on electrets from the fundamentals of charging to applications of these materials in studying self-assembly using electrostatic interactions. ii) Bubbles in Microchannels. The realization that systems of bubbles in microchannels represented a major opportunity to study complexity in a very tractable system, and the development of a semi-quantitative theory of this subject. iii) Flames. The growth of “flames” remains an exploratory subject for the research, although their currently relatively little active work involving it ongoing. iv) Systems with Microorganisms. The removal of work in biological systems from this project. Based on work supported in this program, we now have a significant project on the development of microfluidic tools for studying C. elegans (a nematode), but this work was not appropriate for a program focused on complexity, and we developed separate support for it. (It is, however, an example of successful seeding of a new area by BES.) The work on electrets has gone through a period in which a part of the program was the subject of a MURI; the focus of this work was to develop materials that did not charge electrostatically on friction or contact. The MURI is now over, and the work on dynamic self-assembly (supported by BES) is the major focus. “Flames” has also enjoyed synergistic support, in terms of a project supported by DARPA on flame suppression (in the absence of extinguishing agents, using acoustic and electrostatic interactions). This work was helpful in understanding some of the basics of flames, but is entirely distinct from the BES focus in complexity. A growing interest is in the Lorentz effect. The Lorentz effect is the force exerted on charged particles (electrons, ions, charged molecules) when they move through a perpendicular magnetic field. The Lorentz effect is almost ubiquitous in modern technology: examples of applications include electric motors, dynamos, cathode ray tubes, many batteries, and most systems that control electrical currents with magnetic forces. We have begun to explore the Lorentz effect in electrochemical systems and heterogeneous catalytic systems involving charged organic species and inorganic ions. This work is still at an early stage, but initial studies that Lorentz effects can be large when ions move through magnetic fields, or magnetic fields move in the presence of ions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The 2022 magneto-optics roadmap

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

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Investigating the Effect of Water on the Mechanical Properties of Cellulose from Multiscale Molecular Dynamics Simulations

Classical molecular dynamics (MD) simulations provide insight into the structure and physicochemical properties of materials with atomic resolution. However, the length and time scales accessible to atomistic MD are orders of magnitude smaller than many relevant processes such as the response of a bulk material to experimentally accessible strain rates, which presents challenges when comparing models to experimental measurements. Bottom-up coarse-graining provides a means for systematically mapping atomistic information to lower resolution models to increase the length and time scales achievable by simulation. Cellulose is an abundant carbohydrate biopolymer with applications to many fields of research, such as materials science and renewable energy, due to its desirable mechanical properties and viability for conversion into biofuel. The effect of moisture content on the Young's modulus of cellulose is of special interest due to its native environment often being in the hydrated secondary plant cell wall and the grinding energy requirements for biomass feedstock preprocessing. The current work investigates the effects of water solvent on the Young's modulus of cellulose calculated from coarse-grained MD mechanical stress simulations. The coarse-grained model was parametrized from atomistic MD calculations of cellulose-cellulose potentials of mean force using umbrella sampling techniques under vacuum and solvated conditions. The Young's moduli of the coarse-grained cellulose assemblies parametrized from cellulose in vacuum or solvated in water were computed via mechanical stress simulations to highlight the importance of capturing solvent interactions for modeling the mechanical behavior of cellulose.

BASIC BIOLOGICAL SCIENCES,RADIATION PROTECTION AND↗
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 OSTI2023

Materials Design Directions for Solar Thermochemical Water Splitting

The sustainable, economical production of molecular hydrogen is a crucial component of a net zero-greenhouse-gas-emissions future. Solar thermochemical water splitting (STWS) offers a renewable route to hydrogen with the potential to help decarbonize several industries, including transportation, manufacturing, mining, metals processing, and electricity generation, as well as provide sustainable hydrogen as a chemical feedstock. STWS uses high temperatures generated from concentrated sunlight or other sustainable means for high-temperature heat to produce hydrogen and oxygen from steam. For example, in its simplest form of a two-step thermochemical cycle, a redox-active metal oxide is heated to ≈1700-2000 K, driving off molecular oxygen while producing oxygen vacancies in the material. The reduced metal oxide then cools (ideally with the extracted heat recuperated for re-use) and, in a separate step, comes into contact with steam, which reacts with oxygen vacancies to produce molecular hydrogen while recovering the original state of the metal oxide. Despite its promising use of the entire solar spectrum to split water thermochemically, the current estimated cost of hydrogen produced via STWS is ≈4-6× the U.S. Department of Energy (DOE) Hydrogen Shot target value of $1/kg. One contributing approach to bridging this cost gap is the design of new materials with improved thermodynamic properties to enable higher efficiencies. The state-of-the-art (SOA) redox-active metal oxide for STWS is ceria (CeO 2 ), due to its close to optimal, although too high, oxygen vacancy formation enthalpy and large configurational and electronic entropy of reduction. However, ceria requires high operating temperatures and its efficiency is insufficient. Therefore, efforts to increase the efficiency of STWS cycles have focused on further optimizing oxygen vacancy formation enthalpies and augmenting the reduction entropy via substitution or doping and materials discovery schemes. Examples of the latter include the perovskites BaCe 0.25 Mn 0.75 O 3 and (Ca,Ce)(Ti,Mn)O 3 . These efforts and others have revealed intuitive chemical principles for the efficient and systematic design of more effective materials, such as the strong correlation between the enthalpies of crystal bond dissociation and solid-state cation reduction with the enthalpy of oxygen vacancy formation, as well as configurational entropy augmentation via the coexistence of two or more redox-active cation sublattices. The purpose of this chapter is to prepare the reader with an up-to-date account of STWS redox-active materials, both the SOA and promising newcomers, as well as to provide chemically intuitive strategies for improving their cycle efficiencies through materials design – in conjunction with ongoing efforts in reactor engineering and gas separations – to reach the cost points for commercial viability. First, we will introduce the thermodynamics of STWS using a two-step, metal-oxide, thermochemical cycle with economics in mind. We also will compare the pros and cons of processes that do or do not involve phase changes. Second, we will describe the qualities that make ceria the SOA STWS redox-active material, as well as its limitations. Third, we will survey some of the most promising candidates to date in the search for materials to supplant ceria, emphasizing the post-ternary, metal-oxide-perovskite alloys. Lastly, we will enumerate and discuss the following materials design directions for STWS redox-active materials: crystal reduction potentials as a proxy for oxygen vacancy formation enthalpies, engineering the electronic and configurational entropy of reduction via f-shells and simultaneous redox, and vetting materials stability via temperature-dependent phase diagrams and melting-point prediction.

08 HYDROGEN↗