Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ANL”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Ambient Synthesis for Fe(II) Polypyridyl Complexes with an Order of Magnitude Increase in Charge-Transfer Excited-State Lifetimes over [Fe(bpy) 3 ] 2+

Replacing precious metals with abundant metals is an important research focus in photochemical energy conversion and storage to meet global energy demands. However, transition metal complexes (TMCs) based on abundant 3d metals typically possess photochemical disadvantages─such as short charge-transfer excited-state lifetimes─and molecular modifications have focused on optimizing the interplay of structure, dynamics, and energetics to overcome their limitations. One strategy to do so is the use of bespoke ligands that can extend the lifetimes of chemically useful excited states. Here, in this study, we report the synthesis and characterization of novel Fe(II) complexes featuring lengthy polypyridyl ligands that can be readily synthesized. Steady-state and transient absorption spectroscopies indicate that these complexes have desirable properties and their excited metal-to-ligand charge-transfer states live an order of magnitude longer than in the benchmark [Fe(bpy) 3 ] 2+ . This lifetime is largely preserved in the heteroleptic complexes, thereby enabling the preparation of asymmetric complexes. Additionally, we apply nonradiative transition theory to explain the long-time decay kinetics. In light of their ease of preparation and reasonable excited-state lifetimes, we suggest the use of these complexes in Fe(II) dye-sensitized solar cells, where the rate of charge injection would be competitive with increased lifetime.

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

Synthesis and Characterization of Layered Actinide (U, Np, Pu) Oxide and Hydroxide Phases

Systematic structural comparisons across the early actinides remain limited by the scarcity of well-defined transuranic layered oxide and oxyhydroxide phases. Here, we report the synthesis and single-crystal characterization of new layered actinide compounds spanning U, Np, and Pu obtained under mild hydrothermal conditions in concentrated alkali hydroxide media. These include hydrated oxides α-Cs 2 U 2 O 7 •0.5H 2 O and Rb 2 An 2 O 7 •0.5H 2 O (An = Np, Pu), oxy-hydroxides Rb 6 [(AnO 2 ) 6 O 8 (OH) 2 ]•xH 2 O x = 0, 0.5 (An = U, Np) and Rb 4 [(UO 2 ) 5 O 6 (OH) 2 ]•2H 2 O, as well as nitrate-intercalated compounds Cs 6 [(AnO 2 ) 3 O 4 (OH) 2 ](NO 3 ) 2 (An = Np, Pu). Single-crystal X-ray diffraction studies reveal extended two-dimensional architectures constructed from edge- and vertex-sharing actinyl polyhedra, with systematic evolution in equatorial coordination, hydration, and anionic sheet topology across the U–Np–Pu series. Incorporation of nitrate anions within the interlayer region of Cs 6 [(AnO 2 ) 3 O 4 (OH) 2 ](NO 3 ) 2 establishes a previously unobserved structural motif in layered transuranic oxyhydroxides, demonstrating an additional pathway for anion-mediated framework stabilization. Correlation of crystallographic metrics with single-crystal Raman spectroscopy provides new vibrational benchmarks linking differences in An═O yl bond lengths to the equatorial coordination and the interstitial cations. These findings expand the structural hierarchy of layered actinide materials and address clarifying periodic trends governing topology, bonding, and vibrational signatures in high-valent 5ƒ oxide systems.

actinides↗
DOE OSTI2026

KBKit: A Python Toolkit for Kirkwood–Buff Theory from Molecular Dynamics

Thermodynamic properties of liquid mixtures govern processes that range from drug delivery to energy storage, yet extracting these properties from molecular simulations remains challenging. Kirkwood–Buff (KB) theory offers a rigorous route by linking microscopic pair distribution functions to macroscopic free energies, but practical use of the theory has been hindered by two obstacles: (i) the long simulations needed to obtain well-converged Kirkwood-Buff integrals (KBIs) and (ii) the specialized corrections required to translate finite-size data to the thermodynamic limit. $\texttt{KBKit}$ is an open-source Python package that removes these barriers. It automatically computes KBIs and derived thermodynamic quantities from GROMACS input files, applies state-of-the-art finite-size corrections, and provides built-in diagnostic tools to quantify statistical uncertainty. Written with modern software-engineering practices—continuous integration, extensive unit testing, and thorough documentation—$\texttt{KBKit}$ is both reliable and easy to extend. By condensing complex KBI analysis into a few intuitive commands, $\texttt{KBKit}$ enables researchers to incorporate KB theory into routine simulation workflows and accelerate the discovery of solution-phase thermodynamics.

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

gRASPA

GPU Monte Carlo Simulation Code with a taste of RASPA We present enhancements in Monte Carlo simulation speed and functionality within an open-source code, gRASPA, which uses graphical processing units (GPUs) to achieve significant performance improvements compared to serial, CPU implementations of Monte Carlo. The code supports a wide range of Monte Carlo simulations, including canonical ensemble (NVT), grand canonical, NVT Gibbs, Widom test particle insertions, and continuous-fractional component Monte Carlo. Implementation of grand canonical transition matrix Monte Carlo (GC-TMMC) and a novel feature to allow different moves for the different components of metal-organic framework (MOF) structures exemplify the capabilities of gRASPA for precise free energy calculations and enhanced adsorption studies, respectively. The introduction of a High-Throughput Computing (HTC) mode permits many Monte Carlo simulations on a single GPU device for accelerated materials discovery. The code can incorporate machine learning (ML) potentials. The open-source nature of gRASPA promotes reproducibility and openness in science, and users may add features to the code and optimize it for their own purposes. The code is written in CUDA/C++ and SYCL/C++ to support different GPU vendors. The gRASPA code is publicly available at https://github.com/snurr-group/gRASPA.

Li, Zhao [Purdue/Northwestern/Notre Dame Universit↗
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↗
DOE OSTI2025

Quantum Time Dynamics Mediated by the Yang–Baxter Equation and Artificial Neural Networks

Quantum computing shows great potential, but errors pose a significant challenge. This study explores new strategies for mitigating quantum errors using artificial neural networks (ANNs) and the Yang–Baxter equation (YBE). Unlike traditional error mitigation methods, which are computationally intensive, we investigate artificial error mitigation. We developed a novel method that combines ANNs for noise mitigation combined with the YBE to generate noisy data. This approach effectively reduces noise in quantum simulations, enhancing the accuracy of the results. The YBE rigorously preserves quantum correlations and symmetries in spin chain simulations in certain classes of integrable lattice models, enabling effective compression of quantum circuits while retaining linear scalability with the number of qubits. This compression facilitates both full and partial implementations, allowing the generation of noisy quantum data on hardware alongside noiseless simulations using classical platforms. By introducing controlled noise through the YBE, we enhance the data set for error mitigation. We train an ANN model on partial data from quantum simulations, demonstrating its effectiveness in mitigating errors in time-evolving quantum states, providing a scalable framework to enhance quantum computation fidelity, particularly in noisy intermediate-scale quantum (NISQ) systems. We demonstrate the efficacy of this approach by performing quantum time dynamics simulations using the Heisenberg XY Hamiltonian on real quantum devices.

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

Enabling Multireference Calculations on Multimetallic Systems with Graphic Processing Units

Modeling multimetallic systems efficiently enables faster prediction of desirable chemical properties and the design of new materials. This work describes an initial implementation for performing multireference wave function method localized active-space self-consistent field (LASSCF) calculations through the use of multiple graphics processing units (GPUs) to accelerate time-to-solution. Density fitting is leveraged to reduce memory requirements, and we demonstrate the ability to fully utilize multi-GPU compute nodes. Performance improvements of 5–10x in total application runtime were observed in LASSCF calculations for multimetallic catalyst systems up to 1200 AOs and an active space of (22e,40o) using up to four NVIDIA A100 GPUs. Furthermore, written with performance portability in mind, a comparable performance is also observed in early runs on the Aurora exascale system using Intel Max Series GPUs.

Algorithms↗
DOE OSTI2026

Reducing the Cost of CCSD Basis Set Extrapolation in Ab Initio Computational Thermochemistry

Here, a series of approximations to CCSD contributions in computational model chemistries is presented in the context of kcal mol –1 , kJ mol –1 , and 20 cm –1 theoretical predictions of total atomization energies, benchmarked within the HEAT+CH 4 test suite. A specific set of circumstances where MP2, without empirical scaling, may be used as an effective intermediate in the first two of these accuracy ranges was determined. However, SDQ-MP4, a method long used in pursuit of kcal mol –1 accuracy but relatively unstudied in the subchemical accuracy community, offers significant improvement over the quality of MP2 as a basis-set intermediate at significantly reduced cost compared to CCSD. Given this, we argue for SDQ-MP4 as the de facto CCSD basis-set intermediate in sub-chemical accuracy calculations when CCSD in a desired basis set becomes unaffordable. We additionally report on a “CBS-like” scheme, where MP2 and SDQ-MP4 are used in conjunction to create a “cheap” three-part approximation of large CCSD basis set limits. The data for the CCSD approximation schemes are organized in such a way that model chemistry developers can locate an analog of their current approach for the CCSD basis set limit and explore alternative intermediates that either decrease computational cost or increase computational accuracy. We also show, for a handful of molecules, that SDQ-MP4 shows promise as an effective basis-set intermediate for harmonic and fundamental frequency computations, allowing for zero-point corrections of nearly CCSD(T)/ANO1 quality using simple composite methods that only require CCSD(T)/ANO0.

Thorpe, James H. [Argonne National Laboratory (ANL↗
DOE OSTI2026

Hierarchical Reinforcement Learning of a Short-Range Bond-Order Potential for Silica: Analytic Embedding of Coordination with Classical Efficiency

Reinforcement learning (RL) has recently emerged as a data-efficient strategy to parametrize short-range interatomic potentials. Building on our past RL optimization of pairwise silica models, we extend the framework to a bond-order (Tersoff-type) potential that provides an analytic embedding of local coordination through a three-body term. A hierarchical RL workflow combining continuous-action Monte Carlo Tree Search and property-based rewards efficiently explores the 26-dimensional parameter space, sequentially optimizing lattice parameters, densities, angles, and cohesive energies of 21 silica polymorphs. The resulting models, Q-Tersoff and ML-Tersoff, reproduce the energetic ordering of low-energy phases and capture the angular correlations and amorphous structure factors of silica with improved fidelity over pairwise force fields, while remaining orders of magnitude faster than high-dimensional machine-learned potentials. Both models underperform for elastic constants and high-energy frameworks, delineating the limits of the current analytic form. The approach establishes a general and interpretable route to angle-aware, short-range potentials that bridge physics-based and machine-learned descriptions of silicate materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Comparison of Electronic Structure Methods for Predicting the Hydrogenation Energies of Candidate Molecules for Hydrogen Storage

The development of novel energy materials and fuels is required to expand current available energy sources. Aiming to reach this goal, there is growing interest in using molecular hydrogen as an energy carrier due to its abundance and high energy density. Liquid organic hydrogen carriers (LOHCs) are a promising route to the large-scale storage and transport of hydrogen for use in the energy economy. The search for thermodynamically viable LOHC molecules for real world use has led to a set of constraints on the dehydrogenation enthalpy and the minimum gravimetric hydrogen capacity. These constraints allow one to formulate the search for an ideal LOHC candidate molecule as an optimization problem well suited to the strengths of machine learning and artificial intelligence computational approaches. A critical barrier to a large-scale, high-throughput screening of LOHC candidate molecules is the lack of reliable training data. Computational electronic structure methods including density functional theory, coupled cluster approximations, and diffusion Monte Carlo can be used to provide training data where experimental data are either unreliable or do not exist. In this work, we use these methods to calculate the dehydrogenation energies and enthalpies of candidate LOHC molecules.

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

Nanosecond Structure of Radical Pair Intermediates from High-Frequency Quantum Oscillations: Insight into the Q A •– to Q B Electron Transfer Step in Purple Bacterial Photosynthesis

We demonstrate the validity of our approach to deduce, from the anisotropy of quantum oscillations, the geometry of short-lived radical pair intermediates in photosynthesis. A global fit of a two-dimensional W-band (94 GHz) electron paramagnetic resonance (EPR) experiment provides the same global minimum values for the geometry of the A-side radical pair P 700 •+ A 1A •− in photosystem I (PSI) as observed in a previous Q-band (34 GHz) EPR study, yet with a significantly increased convergence rate of 62%. This demonstrates that the global fit yields the correct radical pair geometry even at Q-band frequencies. With this information, we revisit our previous Q-band study of the cofactor arrangement of P 865 •+ Q A •− , the stabilized charge-separated state in purple bacterial reaction centers (RCs). Analysis of calculated two-dimensional data sets of P 865 •+ Q A •− reveals that the quantum oscillation technique is unaffected by a mirror ambiguity in disordered solids and thus can provide unambiguous solutions for all five Euler angles of the radical pair geometry. This enables us to elucidate the Q A •− to Q B electron transfer step in purple bacterial photosynthesis, the subject of controversial discussions for more than 25 years. Our results show that this electron transfer step involves a gating mechanism requiring a 60° rotation of the headgroup of Q A •− in its binding pocket.

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

Computational Study of Indium Oxide Photoelectrodes

Using a combination of first-principles molecular dynamics simulations and electronic structure calculations, we characterize the atomistic structure and vibrational properties of a photocatalytic surface of In 2 O 3 , a promising photoelectrode for the production of hydrogen peroxide. We then investigate the surface in contact with water and show that the electronic states of In 2 O 3 are appropriately positioned in energy to facilitate the two-electron water oxidation reaction (WOR) over the competing four-electron oxygen evolution reaction. In conclusion, we further propose that the use of strained thin films interfaced with water is beneficial in decreasing the optical gap of In 2 O 3 and thus in utilizing a wider portion of the solar spectrum for the WOR.

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

Spin Dynamics of SCRPs in ZnO Quantum Dot–Organic Molecule Conjugates Studied with Pulsed EPR

Photogenerated spin-correlated radical pairs (SCRPs) are emerging as promising new candidates for Quantum Information Science applications, where they act as electron spin qubit pairs (SQPs). Historically, these pairs have been mostly investigated in natural photosynthesis and in molecular organic donor–(linker)–acceptor systems. Recently, we have shown that these spin pairs can also be observed by time-resolved electron-paramagnetic resonance (EPR) spectroscopy in hybrid inorganic–organic conjugates. The current study builds on recent work (ACS Nano, 2025, 19, 12194–12207), in which we systematically prepared hybrid inorganic–organic molecule systems with tunable geometries that can host SCRPs/SQPs. Here, we demonstrate using pulsed EPR spectroscopy that we can generate, study, and manipulate the spins in these photogenerated SCRPs hosted on highly tailorable inorganic–organic hybrid systems. Microwave pulse-based spin manipulation is the first step toward developing quantum computing and quantum sensing applications with these promising spin-based qubit materials, which can be reversibly photogenerated in highly spin-polarized states at moderate temperatures. An attractive feature of this new class of hybrid materials is that the g-factor of the unpaired electron spin in the ZnO QD (as part of the SCRP/SQP) can be adjusted using the quantum size effect in these QDs, which allows for the selective addressability of each spin in the SCRP using microwave pulses. Additionally, we also observed a difference between the g-values of the unpaired electron in the ZnO QD in the transient SCRP state and stable photoreduced state. We hypothesize that the unpaired electron in the SCRP is delocalized as a band-like conduction electron, while in the stable photoreduced state, the electron is delocalized in a shallow electron trap state.

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

Investigation of Electrocatalytic CO 2 Reduction on MXene Materials via First-Principles Simulations

Computational studies of CO 2 reduction to yield various products were carried out on the basal plane and edges of different MXene materials. The impact of vacancies upon Mo 2 TiC 2 T x , W 2 TiC 2 T x , and Ti 3 C 2 T x (T x = O and OH) was also examined for both edge and basal sites. Initial calibrations were carried out to generate surfaces with optimal oxide/hydroxide ratios and proper termination sites upon which various vacancy sites were explored to ensure an accurate model of the surfaces’ resting states. From this work, Mo 2 TiC 2 O was determined to exhibit the lowest theoretical overpotential for methane as determined by volcano plot analyses. At a large enough vacancy concentration, the CO 2 reduction reaction (CO 2 RR) is predicted to outcompete the hydrogen evolution reaction (HER) as the predominant reaction on the surface. In conclusion, when examining the edge of Mo 2 TiC 2 O, stronger CO 2 binding was exhibited to such an extent that the reaction was predicted to terminate after the generation of formate on the edge.

Lund, Colton [Argonne National Laboratory (ANL), A↗
DOE OSTI2025

Pressure-Induced Stacking of Rhodamine 6G Molecules

Extended conjugation of molecular π systems are of great significance for solar cells and organic light-emitting diodes as these can serve as new light absorbing and emitting chromophores. This research demonstrates the potential of using pressure and solvent environment to create larger assemblies of molecular π-conjugated systems. Here, we selected rhodamine 6G (Rh6G) to study the pressure-induced stacking of molecules in three different solvents: methanol (MeOH), isopropanol (IPA), and 16:3:1 mixture of methanol: ethanol: water (MEW). Absorption, emission, and fluorescence lifetime were studied at pressures up to 11.38 GPa. We found that pressure-induced stacking mechanisms can present solvent dependence: high pressures tend to create emissive Haggregates in MeOH and MEW, but excimers in IPA. This work guides and demonstrates routes to stacking of organic molecules using high pressures and provides fundamental insights into the formation of extended conjugation π systems.

Johnson, Kade [Univ. of Texas at San Antonio, TX (↗
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 OSTI2025

Tunable Stabilization of Cuprous Ions via Kinetic and Thermodynamic Control of Cu Electrodeposition in Non-Aqueous Media

In this work, significant stabilization of Cu + ions is achieved in non-aqueous electrolytes containing TFSI - and Cl - anions with diglyme as solvent. When Cl - is absent, stabilization of Cu + occurs due to the very low stability constant of Cu 2+ in diglyme. However, in the presence of relatively low Cl - concentrations (i.e., equimolar relative to Cu 2+ ), both Cu + and Cu 2+ form extremely stable Cu 2+ / + [Cl - ] complexes resulting in a similar to 1200 mV negative shift of the Cu + /0 redox couple relative to that of more easily reduced Cu + [TFSI - ][Cl - ] complexes. Furthermore, these results suggest that the lower solvation energies and wider electrochemical stability windows of non-aqueous solvents relative to water enable anions to play a much more significant role in guiding complexation behavior-providing new possibilities for (electro)chemical stabilization of reactive intermediates and highlighting the wealth of unexplored opportunities for electrolyte design in non-aqueous systems.

Anion Association↗
DOE OSTI2025

Characterizing Ultrafast Intersystem Crossing Pathways in Molecular Pt Dimers Using Time-Resolved Wide-Angle X-ray Scattering

Vibronic coupling between transition metal charge transfer states is a potential mechanism for enhancing the intersystem crossing (ISC) rate. Vibronic coupling-driven ISC has been observed in Pt­(II) dimer complexes, where the trajectory across excited-state pathways is tuned by atomic displacements via Pt–Pt stretching vibrations. Time-resolved wide-angle X-ray scattering (TR-WAXS) was utilized to quantify the Pt–Pt contraction following metal–metal-to-ligand charge transfer (MMLCT) excitation in Pt dimers with different bridging ligands. Both complexes exhibit Pt–Pt bond formation with a decrease in Pt–Pt distance of ∼ 0.25 Å and coherent vibrational wavepackets (CVWPs) encoded in the Pt–Pt contraction of both dimers. However, the complexes exhibit different time-dependent evolution of their CVWPs. Analysis of interference patterns between different CVWPs is used to track the trajectory across the excited-state surfaces. Furthermore, this work demonstrates that the interference between CVWPs in ultrafast TR-WAXS encodes indirect information regarding electronic excited-states to reveal the Pt dimer bridge-dependent ISC mechanism.

Computational chemistry↗