Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational chemistry”

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

MBX V1.2: Accelerating Data-Driven Many-Body Molecular Dynamics Simulations

The MBX software provides an advanced platform for molecular dynamics simulations, leveraging state-of-the-art MB-pol and MB-nrg data-driven many-body potential energy functions. Developed over the past decade, these potential energy functions integrate physics-based and machine-learned many-body terms trained on electronic structure data calculated at the "gold standard" coupled-cluster level of theory. Recent advancements in MBX have focused on optimizing its performance, resulting in the release of MBX v1.2. While the inherently many-body nature of MB-pol and MB-nrg ensures high accuracy, it poses computational challenges. MBX v1.2 addresses these challenges with significant performance improvements, including enhanced parallelism that fully harnesses the power of modern multicore CPUs. In conclusion, these advancements enable simulations on nanosecond time scales for condensed-phase systems, significantly expanding the scope of high-accuracy, predictive simulations of complex molecular systems powered by data-driven many-body potential energy functions.

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

Ultrafast Dynamics and Rearrangement of the EUV Photoacid Generator Phenyl Triflate

Extreme UV photolithography is a high-energy process crucial for modern computer chip fabrication, where photoacid generators (PAGs) enhance radiation exposure by producing reactive chemical species. This study investigates the ultrafast dissociative dynamics of phenyl triflate, a widely used PAG, under high-energy ionizing conditions. The dissociative ionization of the phenyl triflate cation, which includes a molecular rearrangement that releases SO 2 and forms the phenyl trifluoromethyl ether cation, is measured over the first 5 ps. These dynamics reveal vibrational coherence corresponding to a torsional mode of the cation, involving the twisting of the phenyl group into the O–S–C plane. Electronic structure calculations for the radical cation are in good agreement with the experimentally observed vibrational coherence. These findings provide valuable insights into the behavior of phenyl triflate under ionizing conditions similar to its industrial usage.

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

Is the Matrix Completion of Reduced Density Matrices Unique?

Reduced density matrices are central to describing observables in many-body quantum systems. In electronic structure theory, the two-particle reduced density matrix (2-RDM) suffices to determine the energy and other key properties. Recent work has used matrix completion, leveraging the low-rank structure of RDMs and approximate theoretical models, to reconstruct the 2-RDM from partial data and thus reduce the computational cost. However, matrix completion is, in general, an under-determined problem. Revisiting Rosina’s theorem (Rosina, M. Queen’s Papers on Pure and Applied Mathematics , 1968, No. 11, 369), we here show that the matrix completion is unique under certain conditions, identifying the subset of 2-RDM elements that enables its exact reconstruction from incomplete information. Building on this, we introduce a hybrid quantum–stochastic algorithm that achieves exact matrix completion, demonstrated through applications to the Fermi–Hubbard model.

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

Renormalization of states and quasiparticles in many-body downfolding

We explore the principles of many-body Hamiltonian complexity reduction via downfolding on an effective low-dimensional representation. We show that the renormalization factor provides a unique measure of the quality of the compression as it directly represents the projection between the approximate stationary state of the many-body Hamiltonian and the full many-body wavefunction. Hence, the renormalization factor is a measure of fidelity between the effective (reduced-rank) description and the full many-body treatment for arbitrary (i.e., ground and excited) states. When the entire problem is mapped on a system of interacting quasiparticles [Romanova et al., npj Comput. Mater. 9, 126 (2023)], the effective Hamiltonians can faithfully reproduce the physics only when a clear energy scale separation exists between the subsystems and their environment. We also demonstrate that it is necessary to include quasiparticle renormalization at distinct energy scales, capturing the distinct interaction between subsystems and their surrounding environments. Numerical results from simple, exactly solvable models highlight the limitations and strengths of this approach, particularly for ground and low-lying excited states. This work lays the groundwork for applying dynamical downfolding techniques to problems concerned with (quantum) interfaces.

Green-functions technique↗
DOE OSTI2025

Binding Free Energy Analysis of Colicin D, E3 and E8 to Their Respective Cognate Immunity Proteins Using Computational Simulations

Colicins are antimicrobial proteins produced by bacteria for the purpose of destroying neighboring bacteria. Colicin activity is neutralized by a specific cognate immunity protein in order to protect the host. This study investigates the structural and binding mechanisms underlying the interaction of colicin-D, -E3 and -E8 to their respective immunity proteins (ImD, Im3 and Im8) using structure prediction, molecular dynamics (MD) simulations and MM-PBSA approach of free energy calculations. High-confidence colicin-immunity (Col-Im) complex structures predicted using AlphaFold2 were subjected to MD simulations of 150 ns with GROMACS and were analyzed for the binding free energy calculation using gmx_MMPBSA. Results showed that the complex of Col_E3-Im3 exhibited the most favorable binding free energy, driven by strong van der Waals and electrostatic interactions. Col_D-ImD and Col_E8-Im8 also showed the favorable binding. Electrostatics and hydrogen bonding emerged as a key factor driving binding and stability, while polar solvation acted as a destabilizing factor across all systems. These outcomes provide an understanding of the molecular mechanisms of Col-Im systems, with potential applications for developing natural antimicrobials for food safety.

Biochemistry & Molecular Biology↗
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

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry

Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.

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

Size-Transferable Prediction of Excited State Properties for Molecular Assemblies with a Machine Learning Exciton Model

Computational modeling of the excited states of molecular aggregates faces significant computational challenges and size heterogeneity. Current machine learning (ML) models, typically trained on specific-sized aggregates, struggle with scalability. We found that the exciton model Hamiltonian of large aggregates can be decomposed into dimer pairs, allowing an ML model trained on dimers to reconstruct Hamiltonians for aggregates of any size. We also proposed a new method to address the phase-correction problem by introducing coupling terms’ approximations. Our model accurately predicted the excitation energies of the trimer and tetramer of perylene and tetracene and estimated S1 oscillator strengths of perylene aggregates. Leveraging our ML model, the optical gaps of nanosized perylene aggregates with up to 50 monomers are analyzed, qualitatively revealing the role of different couplings on their size dependency. Future work will explore transferability across different monomers to predict optical properties in heterogeneous assemblies.

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

Structurally Driven Selective Adsorption of Hydrocarbons by Metal Substitution in Isostructural Rare-Earth Metal–Organic Frameworks

The design and realization of highly selective nanoporous materials are necessary to target critical separations across industries. By leveraging pore size, pore shape, and linker functionalization, the design of nanoporous solid adsorbents will enable the rapid production of energy efficient separation materials for high-value gas mixtures. This study uses a combination of modeling, synthesis, and gas adsorption testing to investigate a new class of small-pore isostructural rare-earth (RE) 2,5-dihydroxyterephthalic acid (DOBDC) metal–organic frameworks (MOFs) (RE: Pr-, Gd-, Er-, Yb; DOBDC = 2,5-dihydroxyterephthalic acid) and their adsorption selectivity for acetylene/ethylene mixtures. Density functional theory simulations identified that selective binding of acetylene over ethylene in the Gd-, Er-, and Yb-DOBDC MOFs was due to hydrogen-bonding between acetylene and the linker hydroxyl. Adsorption experiments validated the computational results by identifying mechanisms that control the acetylene/ethylene adsorption selectivity and high acetylene adsorption. Furthermore, dynamic column breakthrough experiments with the Gd-DOBDC MOF validated the simulations and indicated that ethylene can be separated from acetylene in a mixture containing 1 vol % acetylene and 39 vol % ethylene (balance argon). In conclusion, the results highlight the complexity of gas binding in functional porous materials and how combining modeling and experiment enables a fundamental understanding of gas–framework interactions that can be leveraged for the design of future separation materials.

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

Anion–Cation–Anion Ion Triplet Characterization by Computation and Photoelectron Spectroscopy

Ion triplets of the chloride salts of two commonly used weakly coordinating cations are reported (i.e., Cl – ·NMe 4 + Cl – (1 – ) and Cl – ·PPh 4 + Cl – (2 – )). Negative ion photoelectron spectra at 20 K afford vertical and adiabatic detachment energies of 5.18 and 5.00 eV (1 – ) and 5.03 and 4.70 eV (2 – ), respectively. These results are well reproduced by coupled cluster calculations with single, double, and perturbative triple excitations (CCSD(T)) whereas M06-2X is systematically too small by ~0.3 eV (i.e., 7 kcal mol –1 ). The structures of both 1 – and 2 – have five or six C–H···Cl – interactions that stabilize these cluster anions by 32 (1 – ) and 27 (2 – ) kcal mol –1 as given by their chloride dissociation enthalpies. In conclusion, these values drop to 7.4 and 3.8 kcal mol –1 in dichloromethane based up conductor-like polarizable continuum model calculations and suggest that X – ·M + X – ion triplets with a weakly coordinating cation maybe the reactive form of salts under some conditions.

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

Photon counting with intensified charge coupled device (ICCD) – I. In-depth detector characterization

While the adsorption properties of transition metal catalysts have been widely studied, leading to the discovery of various scaling relations, descriptors of catalytic activity, and well-established computational models, a similar understanding of semiconductor catalysts has not yet been achieved. In this work, we present a high-throughput density functional theory investigation into the adsorption properties of 5 oxides of interest to the photocatalytic CO 2 reduction reaction: TiO 2 (rutile and anatase), SrTiO 3 , NaTaO 3 , and CeO 2 . Using a systematic approach, we exhaustively identify unique surfaces and construct adsorption structures to undergo geometry optimizations. We then perform a data-driven analysis, which reveals the presence of weak adsorption energy scaling relations, the propensity of adsorbates of interest to interact with oxygen surface sites, and the importance of slab deformation upon adsorption. Our findings are presented in the context of experimental observations and in comparison to previously studied classes of catalysts, such as pure metals and tellurium-containing semiconductors, and reinforce the need for a comprehensive approach to the study of site-specific surface phenomena on semiconductors.

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

Non-stoichiometry Governs the Pathway from Amorphous to Crystalline Calcium Carbonate

The controlled crystallization of calcium carbonate underlies the elaborate architectures of marine corals, the design of biomimetic materials, and the global carbon cycle. Despite its ubiquity, the chemical mechanisms that govern the crystallization of calcium carbonate from its amorphous precursor, amorphous calcium carbonate (ACC), remain elusive, largely due to the difficulty in resolving the structure and chemistry of this transient amorphous state. Here, we use time-lapse photography, image analysis, spatially resolved pH determination, in situ synchrotron pair distribution function (PDF) analysis, and dynamic nuclear polarization (DNP) solid-state nuclear magnetic resonance (NMR) spectroscopy to reveal pervasive compositional variability in ACC that underpins its metastability and crystallization. By evaluating the kinetics of the amorphous-to-crystalline transition across different solution chemistries, we demonstrate that ACC is non-stoichiometric and CO 3 -deficient. Counterions from the precursor (e.g., NO 3 - from Ca(NO 3 ) 2 ) substitute into the ACC network, displacing CO 3 ions and mediating both ACC stability and transformation. Crystallization proceeds through refinement of the stoichiometry toward CaCO 3 and uptake of free CO 3 2- anions. Furthermore, these findings are consistent across a wide range of concentrations and different carbonate sources, explaining the diverse behaviors observed for ACC and providing a chemical framework for controlling calcium carbonate crystallization.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Surface-Controlled TiO 2 Nanocrystals with Catalytically Active Single-Site Co Incorporation for the Oxygen Evolution Reaction

The design of advanced electrocatalysts is often hindered by uncertainties in identifying and controlling the active surfaces and catalytic centers within heterogeneous materials. Here we present the synthesis of single-site Co catalysts, substitutionally doped into surface-controlled TiO 2 anatase nanocrystals, aimed at enhancing the oxygen evolution reaction (OER). Grand canonical quantum mechanics calculations reveal that the kinetics of the OER, following an adsorbate evolution mechanism, is markedly influenced by the coordination environment of Co. The simulations suggest significantly higher turnover frequencies when Co is doped into the (001) surface of TiO 2 compared to the (101) surface. Consistent with the computational findings, experimental results show that Co-doped TiO 2 (Co-TiO 2 ) nanoplates with selectively exposed {001} surfaces exhibit enhanced current densities and turnover frequencies compared to Co-TiO 2 nanobipyramids with {101} surfaces. This study highlights the synergy between theoretical calculations and precision synthesis in the development of more effective catalysts.

25 ENERGY STORAGE↗
DOE OSTI2025

Predicted Ferromagnetism in Discovered Co–Bi Binary Phases

Binary solid-state materials offer unique insight into how the interplay of factors such as stoichiometry and bonding interactions affects magnetism and electronic properties. In this work, we considered systems where a transition metal provides the spin moment and a heavy main group element bolsters strong spin–orbit coupling. Within this context, cobalt, a known component of permanent magnets, and bismuth, functionally the heaviest element stable to radioactive decay, form a compelling combination. The Co–Bi system has been previously shown to exhibit superconductivity in a phase recovered from high pressure. We expected the Co−Bi system could also be ferromagnetic, resulting in two sets of compounds within one chemical system, one superconducting and one ferromagnetic. Subsequently, we investigated the Co–Bi system through both experimental and theoretical approaches to discover new candidates for permanent magnets. Ab initio random structure searching calculations identified five new compounds with diverse structural motifs that may form at higher pressures than previously reported. Experimental high-pressure synthesis yielded four compounds: α-CoBi, α-CoBi 2 , β-CoBi, and β-CoBi 2 . Three of these phases, α-CoBi 2 , β-CoBi, and β-CoBi 2 , were consistent with the calculated structures, corresponding to a 60% success rate for our structure search and underscoring the strength of combining computation with experiment. Theory predicts β-CoBi and β-CoBi2 are ferromagnetic, with β-CoBi possessing larger magnetocrystalline anisotropy energy than familiar permanent magnets such as CoPt and Nd–Fe–B. These results suggest the Co–Bi system could be a platform for understanding the factors that underpin magnetism and, to an extent, superconductivity in a chemically simple binary system.

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

Charge Transfer and Recombination Pathways through Fullerene Guests in Porphyrin-Based MOFs

Porphyrin-based metal–organic frameworks (MOFs) offer a unique platform for building porous donor–acceptor networks that exhibit long-lived charge separation and transport upon incorporation of electron acceptor guest species. Here, in this work, porphyrin-based MOFs, PCN-222(H 2 ) and PCN-222(Zn), synthesized as nanoparticle suspensions, are successfully infiltrated with fullerene acceptor molecules, C 60 or PC 61 BM, in both polar and nonpolar solvent environments. The location and relative binding strength of these guest species are evaluated through a combination of N 2 physisorption measurements, photoluminescence quenching, and UV–vis absorption titration experiments. Semiempirical tight binding calculations are used to screen potential locations of the fullerene guest within the MOF pores, and hybrid density functional theory (DFT)-computed interaction energies confirm the energetically favorable positions. The fundamental photophysics of these donor–acceptor host–guest combinations are probed using ultrafast transient absorption spectroscopy. Sub-picosecond electron transfer involving initial exciplex population is observed, with slow charge recombination lifetimes on the order of τ ~1 ns for all systems in both dimethylformamide and 1,4-dioxane. Charge recombination occurs through population of fullerene and/or framework porphyrin triplet states depending on the porphyrin metalation status. The photophysics of the fullerene-loaded MOFs are discussed in the context of relevant porphyrin–fullerene donor–acceptor molecules to highlight the unique role of the framework environment in dictating photoinduced electron transfer and decay pathways.

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

Experimental and Computational Evaluation of Nicotinamide Cofactor Biomimetics

Oxidoreductase enzymes are widely used biocatalysts due to their high enantioselectivity and broad substrate compatibility in useful transformations. Many oxidoreductases require nicotinamide cofactors (i.e., NAD(P)H). To replace this costly natural cofactor, synthetic nicotinamide cofactor biomimetics (NCBs) offer different shapes, binding affinities, and reducing potentials that exceed the capabilities of wild-type NAD(P)H. However, the ill-defined structure–activity relationships (SARs) of various NCBs slow rationally guided innovation, such as customized reducing potentials. Here, we dissect two essential elements of NCB design, holding the nicotinamide invariant. First, the linker length between the nicotinamide and an unconjugated aromatic ring uncovered unexpected benefits to redox activity for two or three carbon linkers. Second, substitution on this unconjugated aryl group (Ring 2) might not be expected to affect activity. However, SAR trends demonstrate substantial benefits to reductive potential conferred by electron-donating functionalities on Ring 2. Furthermore, catalysis by two enzymes demonstrates enzyme-dependent tolerance or sensitivity to the NCB structures. Density functional theory (DFT) and computational modeling provide a theoretical framework to understand and build upon these observations. Ring 2 reaches up to the nicotinamide to stabilize its positive charge after oxidation through π–π stacking and charge transfer. Thus, the systematic examination of NCB’s stability, electrochemical redox potentials, and kinetics uncovers trends for the improved design of NCBs.

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

Development, Application, and Mechanistic Interrogation of a Dual Ni Catalysis Approach to Photoredox-Based C(sp 3 )–C(sp 3 ) Cross-Coupling

The installation of alkyl substituents such as methyl groups is a crucial tactic for the synthesis and diversification of medicinal compounds with improved pharmacological profiles. However, strategies that leverage methyl radical for C(sp 3 ) methylation remain underdeveloped due to the challenge of obtaining cross-selectivity between fleeting aliphatic radicals and alkyl electrophiles. Here, we report the development, application, and interrogation of a conceptually novel mechanistic framework for C(sp 3 )–Me bond formation using two distinct Ni catalysts capable of cross-coupling sterically and electronically diverse alkyl halides (chlorides and bromides) with methyl radical generated photocatalytically from benzaldehyde dimethyl acetal. Furthermore, by modifying the alkyl substituents on the acetal coupling partner, we demonstrate cross-couplings beyond methylation to access an array of 1°–1° and 1°–2° alkyl–alkyl bonds. Experimental and computational mechanistic studies provide support for cooperativity between an in situ-generated (bpy)Ni I (X) catalyst that facilitates XAT and inner-sphere C–C bond formation and a (Tp*)Ni II (acac) cocatalyst that captures methyl radical and engages in concurrent outer-sphere S H 2 coupling.

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

Computational methods based on density functional theory for reactions and processes involving electronic spin (Final Technical Report)

This award supports one post-doctoral researcher for 1.5 years. Publications that acknowledge this grant: Refs. 1–14. Refs. 3,10,12 assess current methodology for the evaluation of magnetic exchange couplings in transition metal complexes. In particular, Ref. 10 validates the use of an approximate (non-iterative) Green’s function approach for the calculation of magnetic exchange couplings and will be the foundation for Thrust 2 in this proposal. Refs. 3 and 12 focus on widely used density functional approaches based on the standard energy differences methodology for the particular case of oxo-bridged Fe(III) complexes. Refs. 2,4–7,11 apply current methodologies to problems of practical interest in molecular magnetism. Ref. 13 presents a methodology to explicitly simulate the dynamics of open quantum systems within density functional theory (DFT) calculations based on the Liouville-von Neumann equation of motion for quantum systems driven out-of-equilibrium. Ref. 8 uses non-collinear spin DFT to explain the mechanical behavior of magnetic mono-atomic Pt wires produced in break-junction experiments in the presence of a magnetic field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗