Search NASA⌕ Search

SEARCH · Search NASA

Results for “Analytic linearization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Axial H-Bonding Solvent Controls Inhomogeneous Spectral Broadening, While Peripheral H-Bonding Solvent Controls Vibronic Broadening: Cresyl Violet in Methanol

The dynamics of the nuclei of both a chromophore and its condensed-phase environment control many spectral features, including the vibronic and inhomogeneous broadening present in spectral line shapes. For the cresyl violet chromophore in methanol, we here analyze and isolate the effect of specific chromophore–solvent interactions on simulated spectral densities, reorganization energies, and linear absorption spectra. Employing both chromophore and its condensed-phase environment control many spectral features, including the vibronic and inhomogeneous broadening present in spectral line shapes. For the cresyl violet chromophore in methanol, we here analyze and isolate the effect of specific chromophore–solvent interactions on simulated spectral densities, reorganization energies, and linear absorption spectra. Employing both force field and ab initio molecular dynamics trajectories along with the inclusion of only certain solvent molecules in the excited-state calculations, we determine that the methanol molecules axial to the chromophore are responsible for the majority of inhomogeneous broadening, with a single methanol molecule that forms an axial hydrogen bond dominating the response. Furthermore, the strong peripheral hydrogen bonds do not contribute to spectral broadening, as they are very stable throughout the dynamics and do not lead to increased energy-gap fluctuations. We also find that treating the strong peripheral hydrogen bonds as molecular mechanical point charges during the molecular dynamics simulation underestimates the vibronic coupling. Including these peripheral hydrogen bonding methanol molecules in the quantum-mechanical region in a geometry optimization increases the vibronic coupling, suggesting that a more advanced treatment of these strongly interacting solvent molecules during the molecular dynamics trajectory may be necessary to capture the full vibronic spectral broadening.

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

RT-EZ: A Golden Gate Assembly Toolkit for Streamlined Genetic Engineering of Rhodotorula toruloides

For economic and sustainable biomanufacturing, the oleaginous yeast Rhodotorula toruloides has emerged as a promising platform for producing biofuels, pharmaceuticals, and other valuable chemicals. However, genetic manipulation of R. toruloides has been limited by its high GC content and the lack of a replicating plasmid, necessitating gene integration into the genome of the yeast. To address these challenges, we developed the RT-EZ (R. toruloides Efficient Zipper) toolkit, a versatile tool based on Golden Gate assembly, designed to streamline R. toruloides engineering with improved efficiency and flexibility. The RT-EZ toolkit simplifies vector construction by incorporating new features such as bidirectional promoters and 2A peptides, color-based screening using RFP, and sequences optimized for both Agrobacterium tumefaciens-mediated transformation (ATMT) and easy linearization, enabling straightforward selection and transformation. Notably, the RT-EZ kit can be used to construct an expression cassette with four different genes in one assembly reaction, significantly improving vector construction speed and efficiency. The utility of the RT-EZ toolkit was demonstrated through the successful synthesis of arachidonic acid in R. toruloides by coexpressing fatty acid elongases and desaturases. Furthermore, this result underscores the potential of the RT-EZ toolkit to advance synthetic biology in R. toruloides, providing a streamlined method for addressing genetic engineering challenges in the yeast.

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

Twisted Tin‐Chloride Perovskite Single‐Crystal Heterostructures

Self-assembly affords simpler synthetic routes to heterostructures compared with manual layer-by-layer stacking, yet controlling interlayer twist angles in a bulk solid remains an outstanding challenge. We report two new single-crystal heterostructures: (Sn 2 Cl 2 )(CYS) 2 SnCl 4 (CYS = + NH 3 (CH 2 ) 2 S – ; Sn_CYS) and (Sn 2 Cl 2 )(SeCYS) 2 SnCl 4 (SeCYS = + NH 3 (CH 2 ) 2 Se – ; Sn_SeCYS) synthesized in solution, with alternating perovskite and intergrowth layers. Notably, compared to the recently reported lead analog, (Pb 2 Cl 2 )(CYS) 2 PbCl 4 (Pb_CYS), the tin heterostructures feature a twist between the perovskite and intergrowth layers. We trace this twist to local distortions at the Sn centers, which change the interfacial lattice-matching requirements compared to those of the Pb analog. Electronic band structure calculations show that the striking differences in the relative energies of perovskite- and intergrowth-derived bands in Sn_CYS and Pb_CYS arise from structural and not compositional differences. The structural anisotropy of Sn_CYS is also reflected in a large in-plane photoluminescence linear anisotropy ratio. Interfacial strain further affords differential incorporation of Pb into the perovskite and intergrowth layers of the Sn heterostructures, resulting in redshifted optical absorption onsets. Thus, we posit that local structural distortions may be exploited to manipulate the twist angle and interfacial strain in bulk heterostructures, providing a new handle for tuning the band alignments of bulk quantum-well electronic structures.

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

Differential absorption of circularly polarized light by a centrosymmetric crystal

Crystalline solids are governed by universal structure-property relationships derived from their crystal symmetry, leading to paradigmatic rules on what properties they can and cannot exhibit. A long-held structure-property relationship is that centrosymmetric crystals cannot differentially absorb circularly polarized light. In this study, we demonstrate the design, synthesis, and characterization of the centrosymmetric material Li 2 Co 3 (SeO 3 ) 4 , which violates this relationship not by defying symmetry-imposed selection rules but by invoking a photophysical process not previously characterized for crystalline solids. This process originates from an interference between linear dichroism and linear birefringence, referred to as LD-LB, and involves strong chiroptical signals that invert upon sample flipping. Here, in addition to enabling a chiroptical response under centrosymmetry, this process opens up photonic engineering opportunities based on crystalline solids.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Machine-Learned Force Field for Molecular Dynamics Simulations of Nonequilibrium Ammonia Synthesis on Iron Catalysts

Ammonia (NH 3 ) is one of the most important industrial chemicals. The conventional NH 3 synthesis method-the Haber–Bosch process-converts atmospheric nitrogen (N 2 ) into NH 3 using H 2 with an iron (Fe) catalyst. However, this process requires high pressures (100–200 atm) and temperatures (700–800 K) near thermal equilibrium. Recently, Fe-based nanocatalysts have been reported to produce promising NH 3 yields under atmospheric pressures and temperature-modulated nonequilibrium conditions. Understanding the mechanism of nonequilibrium catalysis with programmed temperature variation could help to optimize this fully electrified and less energy-intensive process. Although reactive molecular dynamics (RMD) simulations can be a useful tool to model nonequilibrium catalytic processes, they require the development of accurate force fields (i.e., interatomic potentials). Here, we present a machine-learned (ML) force field within the Deep Potential MD (DPMD) framework, trained using periodic density functional theory (DFT) calculations, to model NH 3 synthesis on Fe catalysts with various surface adsorbates such as *N, *H, *N 2 , *H 2 , *NH, *NH 2 , and *NH 3 . Here, we generated the DFT data from static models of elementary reactions on the most stable (110) surface of body-centered cubic Fe, which then were augmented by data from constant number of particles–volume–temperature (NVT) DFT-MD trajectories at various temperatures. Finally, we utilized the fully optimized ML force field to investigate reaction dynamics at an Fe(110) surface at linearly increasing temperatures using NVT-DPMD simulations. Our simulations indicate that pulsed temperature ramping could prove favorable for NH3 synthesis. For example, we conducted ramping under multiple sets of conditions: (i) from 900 to 1200 K over periods of 0.1–0.3 ns for Fe surfaces precovered with N or NH along with H; and (ii) from 300 to 600 K over 0.1–0.3 ns for Fe surfaces precovered with NH 3 . While our simulations so far are limited to short time scales (very rapid heating), these observations shed light on the mechanism of the high NH 3 synthesis rate achieved in a novel temperature-modulated nonequilibrium catalytic reactor using pulsed heating and cooling.

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

Efficient Parameterization of Density Functional Tight-Binding for 5 f -Elements: A Th–O Case Study

Density functional tight binding (DFTB) models for f-element species are challenging to parametrize owing to the large number of adjustable parameters. The explicit optimization of the terms entering the semiempirical DFTB Hamiltonian related to f orbitals is crucial to generating a reliable parametrization for f-block elements, because they play import roles in bonding interactions. However, since the number of parameters grows quadratically with the number of orbitals, the computational cost for parameter optimization is much more expensive for the f-elements than for the main group elements. In this work we present a set of efficient approaches for mitigating the hurdle imposed by the large size of the parameter space. A novel group-by-orbital correction functions for two-center bond integrals was developed. With this approach the number of parameters is reduced, and it grows linearly with the number of elements, maintaining the accuracy and the number of parameters, in the case of f elements, by more than 40%. The parameter optimization step was accelerated by means of the mini-batch BFGS method. This method allows parameter optimizations with much larger training sets than other single batch methods. A stochastic optimizer was employed that helped overcome shallow local minima in the objective function. The proposed algorithm was used to parametrize the DFTB Hamiltonian for the Th–O system, which was subsequently applied to the study of ThO 2 nanoparticles. The training set consisted of 6322 unique structures, which is barely feasible with conventional optimization methods. The optimized parameter set, LANL-ThO, displays good agreement with DFT-calculated properties such as energies, forces, and structures for both clusters and bulk ThO 2 . Benefiting from the fewer number of parameters and lower computational costs for objective function evaluations, this new approach shows its potential applications in DFTB parametrization for elements with high angular momentum, which present a challenge to conventional methods.

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

Using Separation-Enhanced Isotope Ratio Mass Spectrometry to Enable Increased Renewable Carbon Content in Transportation Fuels (CRADA 525)

Stable isotope ratio measurements of carbon atoms using isotope ratio mass spectrometry (IRMS) can be an effective tool for quantifying biogenic carbon in co-processed fuels, with results approaching the precision and accuracy of accelerator mass spectrometry (AMS). The lower cost of an IRMS may enable deployment to refineries, improving access and analysis turnaround times (≤2 hours), and, by extension, provide data that can allow process optimization to maximize renewable carbon in desired refinery products. This project explored the integration of chemical separation with IRMS analyses to enable highly detailed tracking of biogenic carbon into fuel product streams separated by boiling point range, chemical class, or specific compound. Forty-nine fuels and fuel components of fossil and biogenic origin, spanning gasoline and diesel boiling point ranges, were received from three refiners and were analyzed for their δ 13 C values via IRMS. Results spanned a 13 C range from ca. 10‰ to 44‰ and reflect materials derived from sustainable sources (e.g., C4 or C3 plants, animal-based pathways, syngas) or from fossil-derived fuels. Common ranges are approximately 18‰ to 9‰ and approximately 30‰ to 20‰ for C4 and C3 plants, respectively, and approximately 34‰ to 24‰ and approximately 70‰ to 33‰ for petroleum-derived fuels and methane, respectively. Fuel-like standards were developed and tested using direct-injection elemental analyzer (EA) IRMS for liquid fuels. This method was compared with the published methods, yielding statistically similar results. Four blend curve sets were produced ranging from 0% to 100% of a fuel containing biogenic carbon, focusing on 0% to 10% biogenic carbon. Linear fits were the most applicable for two of the four blend curve sets; however, two sets were found to exhibit slightly quadratic behavior, which was more pronounced in low biogenic blend samples, necessitating second-order fits. The origin of the slight quadratic behavior remains unclear; however, the discussion points to possible interpretations. CanmetENERGY thoroughly characterized a majority of the samples using one- and two-dimensional gas chromatography (GC and GC×GC, respectively) and other analyses. Selected samples were subjected to solid phase extraction (SPE) for saturate, olefin, aromatic, and polar (SOAP) analysis, and the resulting solvent-diluted fractions containing saturates and aromatics were returned to Pacific Northwest National Laboratory (PNNL), where the solvent was removed via evaporation or physical separation using GC techniques. Characterization and separations provided an understanding of saturate and aromatic content, as well as boiling point ranges for each sample and sample fraction. Samples resulting from SPE were examined using EA-IRMS and gas chromatography combustion IRMS (GC-C-IRMS) analyses. Both approaches suggest that the range in values between end-members can be increased by selecting the paraffinic or aromatic fraction of the end-member or by selecting among individual compounds resulting from GC separation of the paraffinic fractions. Considerable work remains to put these approaches into practice and statistically validate the benefit for using a fraction or individual compound over bulk analysis of a sample. However, initial results suggest that separations provide advantages for samples having blend ratios of less than 10% biogenic blendstocks. 13 C results showed statistically similar biofuel blend results to those obtained at PNNL, although additional work is needed to obtain better reproducibility. Select samples were sent to Los Alamos National Laboratory (LANL) for IRMS measurements and Beta Analytics for AMS measurements. This work suggests that IRMS and AMS yield closely comparable results and in some circumstances, IRMS could serve as a surrogate for AMS. While additional work is needed to better resolve statistical advantages for separations and better show the comparable nature of IRMS and AMS in both the biogenic carbon analysis of bulk chemical classes, initial results from this study suggest that these should be pursued in order to proliferate this approach for quantifying biogenic carbon in transportation fuels to the refinery level, thereby potentially enabling process optimization in co-processing scenarios.

09 BIOMASS FUELS↗
DOE OSTI2025

Detecting isolated resonance curves using fixed frequency voltage control tests

Isolated resonance curves, or isolas, are resonance branches of the harmonically forced system that exist separately from the main nonlinear forced response curve, leading to excessive vibrations. Traditional stepped or swept sine simulations and tests rely on continuation along the frequency parameter, typically resulting in a jump phenomenon along the primary resonance branch, prior to the disconnected isola. The main objective of this research is to propose an approach to identify isolated resonance curves by performing continuation along the input amplitude that initializes the response from a low-amplitude solution in the linear regime. Furthermore, this is achieved with the open-loop fixed frequency voltage control method that continues along the shaker voltage parameter and measures the so-called S-curves, which are theoretically a continuous solution branch that connect to the isola. The methodology is demonstrated on a fixture-wing-pylon assembly with a vibro-impact nonlinearity localized in a pylon subcomponent attachment. Multi-harmonic balance simulations are deployed to compute both the nonlinear forced response curves and S-curves to demonstrate the isola detection strategy on a reduced-order finite element model of the nonlinear system. Swept sine and fixed frequency voltage control tests are then conducted on the physical structure to demonstrate the isola detection experimentally, revealing the existence of the large amplitude vibrations that are undetected in the forces levels and frequencies measured with traditional frequency sweeping.

Characterization and Analytical Technique↗
DOE OSTI2024

In-situ polarization modulation IRRAS investigation of ammonia electrooxidation on Pt-Ir and Pt-Ru nanoparticles prepared on engineered catalyst supports

The catalytic activity and surface reactivity of monometallic Pt and bimetallic Pt-Ir and Pt-Ru nanoparticles, supported on two distinct Engineered Catalyst Supports (ECSs), were investigated for the Ammonia Electrooxidation Reaction (AmER) in alkaline media. XRD measurements confirmed alloy formation between Pt-Ir and Pt-Ru nanoparticles, as indicated by the shift of the (111) reflection to higher 2θ values. Cyclic voltammetry, linear sweep voltammetry, and chronoamperometry experiments were conducted to assess the catalytic activity of the Pt, Pt-Ir, and Pt-Ru electrocatalysts. All bimetallic catalysts exhibited lower onset potentials compared to Pt. The differing Tafel slopes between Pt (74 mV dec⁻¹), Pt-Ir (152 mV dec⁻¹), and Pt-Ru (118–197 mV dec⁻¹) suggest that alloying Pt with Ir or Ru alters the reaction mechanisms. Furthermore, the bimetallic Pt-Ir and Pt-Ru catalysts demonstrated greater tolerance for concentrated ammonia solutions relative to Pt. In-situ Polarization Modulation Infrared Reflection Absorption Spectroscopy (PM-IRRAS) provided insights into the formation of N-H species, azide anions (N₃⁻), and N-O compounds. For the Pt-Ru catalyst, an additional peak around ~3600 cm⁻¹ was observed, corresponding to OH⁻ species. The PM-IRRAS results align with the Gerischer–Mauerer mechanism, indicating that partially dehydrogenated ammonia adsorbates act as active intermediates in the oxidation of ammonia over Pt-Ir and Pt-Ru catalysts.

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

Quantum stabilization of unexpected ordered phases on the honeycomb lattice

In this work, I discuss and showcase the utility of the method called minimally-augmented spin-wave theory (MAGSWT) as a relatively simple semi-analytical approach to the phase diagrams of quantum magnets. It complements numerical methods by providing physical insight into which states are competitive and by yielding approximate phase boundaries that agree well with much more numerically intensive calculations. This approach enabled me to construct the faithful phase diagrams of two paradigmatic honeycomb-lattice models in the quantum S=1/2 limit, the J 1 –J 3 FM-AF and J 1 –J 2 AF models, for the collinear quantum phases that replace or extend the classical ones. The results are in good qualitative and semi-quantitative agreement with state-of-the-art numerical studies for these models, correctly capturing the emergence of unexpected quantum phases and the suppression of classically favored spiral orders by quantum fluctuations. This study provides a much-needed important guidance to the ongoing theoretical and experimental searches of the unconventional quantum states. This work will be of significant and timely interest to both theorists and experimentalists in the field of quantum magnetism, broadly defined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Investigating Particle Phase State Dependencies during a Historic Pacific Northwest Wildfire Event

The function and lifetime of an atmospheric particle are greatly dependent upon its phase state. Studying particle phase states from specific ambient source types is often challenging due to uncontrolled environmental conditions and multiple source contributions. However, a historic Pacific Northwest wildfire event occurred in September 2020, providing a rare opportunity to study a pseudo-single particle type ambient sample (i.e., wildfire) over time. Sampling was performed in Richland, Washington, capturing aged particles from shifting wildfire sources. In addition to source variability, temporal variation in wildfire particle phase state was observed from single particle analyses. Despite this, single particle elemental compositions were homogenous over time (=87% carbonaceous) with no significant temporal variation in organic functional group abundances within individual particles either (e.g., 4% relative standard deviation in alkene contributions). While particle composition is indeed known to affect phase state, water uptake appeared to be the most significant contributor to phase state variability amongst single particles here, as particle aspect ratios (related to phase state) linearly correlated with ambient relative humidity during sampling (R2 = 0.65). Therefore, while limited to a case study, meteorological considerations are found to convey greater comparative importance than particle composition for aged particle phase state predictions.

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

Catalytic closed-loop recycling of polyethylene-like materials produced by acceptorless dehydrogenative polymerization of bio-derived diols

Petroleum-derived polyolefins exhibit diverse properties and are the most important and largest volume class of plastics. However, polyolefins are difficult to efficiently recycle or break down and are now a persistent global contaminant. Broadly replacing polyolefins with bio-derived and degradable polyethylene-like materials is an important yet challenging endeavour towards sustainable plastics. Here, in this study, we report a solution for circular bio-based polyethylene-like materials synthesized by acceptorless dehydrogenative polymerization from linear and branched diols and their catalytic closed-loop recycling. The polymerization and depolymerization processes utilize earth-abundant manganese complexes as catalysts. These materials exhibit a wide range of mechanical properties, encompassing thermoplastics to plastomers to elastomers. The branched diols, produced through a thiol-ene click reaction, can be polymerized to plastics with significantly enhanced tensile properties, toughness and adhesive properties. These materials could be depolymerized back to monomers through hydrogenation and were separatable with a monomer recovery of up to 99%, unaffected by the presence of dyes and additives. Overall, this system establishes a route to more sustainable plastics.

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

Formation of Linear Plasmonic Heterotrimers Using Nanoparticle Docking to DNA Origami Cages

The fabrication of complex assemblies with interesting collective properties from plasmonic nanoparticles (NPs) is often challenging. While DNA-directed self-assembly has emerged as one of the most promising approaches to forming such complex assemblies, the resulting structures tend to have large variability in gap sizes and shapes, as the DNA strands used to organize these particles are flexible, and the polydispersity of the NPs leads to variability in these critical structural features. Here, we use a new strategy termed docking to DNA origami cages (D-DOC) to organize spherical NPs into a linear heterotrimer with a precisely defined geometrical arrangement. Instead of binding NPs to the exterior of the DNA templates, D-DOC binds the NPs to either the interior or the opening of a 3D cage, which significantly reduces the variability of critical structural features by incorporating multiple diametrically arranged capture strands to tether NPs. Additionally, such a spatial arrangement of the capture strand can work synergistically with shape complementarity to achieve tighter confinement. To assemble NPs via D-DOC, we developed a multistep assembly process that first encapsulates an NP inside a cage and then binds two other NPs to the openings. Microscopic characterization shows low variability in the bond angles and gap sizes. Both UV–vis absorption and surface-enhanced Raman scattering (SERS) measurements showed strong plasmonic coupling that aligned with predictions by electrodynamic simulations, further confirming the precision of the assembly. These results suggest D-DOC could open new opportunities in biomolecular sensing, SERS and fluorescence spectroscopies, and energy harvesting through the self-assembly of NPs into more complex 3D assemblies.

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

Photo-Reactive Amine-Based Direct Air Capture and Conversion of CO2

This study presents a novel photochemical approach to direct air capture (DAC) and CO2 conversion, utilizing a ruthenium-modified mesoporous TiO2 composite infused with linear polyethyleneimine (L-PEI) and enhanced by light absorbing and earth abundant titanium nitride (TiN). This light-activated system operates at ambient pressure and addresses long-standing challenges in conventional amine-based thermal reactive carbon capture (RCC), including amine degradation and CO2 slip. Our findings demonstrate that L-PEI effectively stabilizes CO2 adsorption and facilitates high-yield methanation over 50 cycles with a non-flammable forming gas mixture (approximately 5% H2). Comparatively, composites incorporating small-molecule amines exhibited poor stability under illumination, while highly mobilie branched PEI (B-PEI) formulations suffered from significant catalyst deactivation. A technoeconomic analysis suggests that methane synthesis via this platform could be achieved at costs below $5/kg under current electrocatalytic hydrogen pricing. These findings highlight the potential of this approach to enhance energy security by enabling decentralized, scalable fuel production from atmospheric CO2.

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

Importance of Electron Correlation on the Geometry and Electronic Structure of [2Fe–2S] Systems: A Benchmark Study of the [Fe 2 S 2 (SCH 3 ) 4 ] 2–,3–,4– , [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S- p -tol) 4 ] 2– , and [Fe 2 S 2 (S- o -xyl) 4 ] 2– Complexes

Iron–sulfur clusters are crucial for biological electron transport and catalysis. Obtaining accurate geometries, energetics, manifolds of their excited electronic states, and reduction energies is important to understand their role in these processes. Here, using a [2Fe–2S] model complex with Fe II and Fe III oxidation states, which leads to different charges, i.e., [Fe 2 S 2 (SMe) 4 ] 2–,3–,4– , we benchmarked a variety of computational methodologies ranging from density functional theory (DFT) to post-Hartree–Fock methods, including complete active space self-consistent field (CASSCF), multireference configuration interaction, the second-order N-electron valence state perturbation theory (NEVPT2), and the linearized integrand approximation of adiabatic connection (AC0) approaches. Additionally, we studied three experimentally well-characterized complexes, [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S-o-tol) 4 ] 2– , and [Fe 2 S 2 (S-o-xyl) 4 ] 2– , via DFT methods. We conclude that the dynamic electron correlation is important for accurately predicting the geometry of these complexes. Broken symmetry (BS) DFT correctly predicts experimental geometries of low-spin multiplicity, while CASSCF does not. However, BS-DFT significantly overestimates the difference between the low- and high-spin electronic states for a given oxidation state. At the same time, CASSCF underestimates it but provides relative energies closer to the reference NEVPT2 results. Finally, AC0 provides energetics of NEVPT2 quality with the additional advantage of being able to use large CASSCF sizes. NEVPT2 gives the best estimates of the Fe III /Fe III → Fe II /Fe III (4.27 eV) and Fe II /F III → Fe II /F II (7.72 eV) reduction energies. The results provide insight into the electronic structure of these complexes and assist in the understanding of their physical properties.

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

Ab Initio Polariton Spectra of ZnTPP Molecules Collectively Coupled Inside an Optical Cavity

Exciton-polaritons are quasi-particles formed by the quantum mechanical hybridization of electronic and photonic excitations. Despite extensive investigations, a fundamental understanding of molecular polariton spectra and the polariton delocalization from an ab initio theoretical perspective remains elusive. We simulate experimentally measured linear transmission spectroscopy of many Zinc(II) tetraphenylporphyrin (ZnTPP) molecules collectively coupled to a cavity from first principles. Our theoretical approach incorporates many low-lying electronic excitations in ZnTPP molecules, as well as collective light-matter couplings between ZnTPP and the quantized radiation modes, both of which are shown to be the key to accurately recovering the experimental spectra. We further analyzed to what extent the polariton and dark states are delocalized over many molecules, for the first time, using fully ab initio descriptions of the molecules. We finally investigate the line width as a function of detuning, providing new theoretical insights into the experimentally observed motional narrowing behavior. Our work presents first-ofits- kind theoretical studies on molecular polariton spectra, offering a new perspective on molecular polariton formation in realistic ab initio molecular systems whose rich, many-state nature provides spectral features enabled by the high density of electronic states beyond simple quantum optics models.

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

Position-specific kinetic isotope effects for nitrous oxide: a new expansion of the Rayleigh model

Nitrous oxide (N 2 O) is a potent greenhouse gas and the most significant anthropogenic ozone-depleting substance currently being emitted. A major source of anthropogenic N 2 O emissions is the microbial conversion of fixed nitrogen species from fertilizers in agricultural soils. Thus, understanding the enzymatic mechanisms by which microbes produce N 2 O has environmental significance. Measurement of the 15 N/ 14 N isotope ratios of N 2 O produced by purified enzymes or axenic microbial cultures is a promising technique for studying N 2 O biosynthesis. Typically, N 2 O-producing enzymes combine nitrogen atoms from two identical substrate molecules (NO or NH 2 OH). Position-specific isotope analysis of the central (N α ) and outer (N β ) nitrogen atoms in N 2 O enables the determination of the individual kinetic isotope effects (KIEs) for N α and N β , providing mechanistic insight into the incorporation of each nitrogen atom. Previously, position-specific KIEs (and fractionation factors) were quantified using the Rayleigh distillation equation, i.e., via linear regression of δ 15 N α or δ 15 N β against [–f In f / (1 – f)], where f is the fraction of substrate remaining in a closed system. This approach, however, is inaccurate for N α and N β because it does not account for fractionation at N α affecting the isotopic composition of substrate available for incorporation into the β position (and vice versa). Therefore, we developed a new expansion of the Rayleigh model that includes specific terms for fractionation at the individual N 2 O nitrogen atoms. By applying this Expanded Rayleigh model to a variety of simulated N 2 O synthesis reactions with different combinations of normal, inverse, and/or no KIEs at N α and N β , we demonstrate that our new model is both accurate and robust. We also applied this new model to two previously published datasets describing N 2 O production from NH 2 OH oxidation in a methanotroph culture (Methylosinus trichosporium) and N 2 O production from NO by a purified Histoplasma capsulatum (fungal) P450 NOR, demonstrating that the Expanded Rayleigh model is a useful tool in calculating position-specific fractionation for N 2 O synthesis.

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

Learning Latent Representations to Bridge Coarse-Grained and Atomistic Resolutions in Polymer Simulations

We present a machine-learning-based framework for learning reduced-order representations of polymer chain conformations across coarse-grained (CG) and united-atom (UA) fidelities. By employing linear singular value decomposition and nonlinear autoencoders, we compress high-dimensional polymer configurations into latent spaces with minimal loss of structural accuracy. Crucially, we demonstrate a near-perfect linear mapping between CG and UA latent spaces, enabling an efficient super-resolution back-mapping procedure that reconstructs high-fidelity UA configurations from CG simulations. While minor structural inaccuracies occur, they are effectively corrected through a brief molecular dynamics relaxation, forming a practical hybrid machine learning−physics scheme. This approach establishes the key structural prerequisites for accelerated polymer dynamics simulations: a compact and accurate latent encoding of polymer chain conformations and a validated multi-fidelity mapping that permits reconstruction of UA structures from CG configurations. The extension of this framework to explicit time evolution within the latent space, enabling dynamics to be propagated at CG fidelity and decoded to UA resolution only when required, represents a natural and well-motivated direction for future work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗