SEARCH · Search NASA
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.
A new Co–Ni–N–C electrocatalyst for CO 2 reduction, featuring diatomic metal-nitrogen sites on N-doped carbon, has been developed. It shows high CO yield and faradaic efficiency, promising for various electrochemical reactions.
Intra-molecular interactions or chemical bonds represent one of the main distinguishing characteristic of molecular fluids. Development of accurate and practical methods to treat these effects is one of the long standing problems in classical site density functional theory (SDFT). One particular instance when these issues become particularly severe is the case of classical interactions potentials with auxiliary sites or dummy atoms. In this situation current SDFT implementations, such as three-dimensional reference interaction site model (3D-RISM), lead to nonphysical results. We re-examine this issue in this work using our recent reformulation of SDFT. We put forward a simple practical solution to this problem, and illustrate its utility for the case of spherical solutes in diatomic liquids.
Over many years, we have published vibrational spectra of molecular materials at extreme pressure and temperature obtained using shock compression and coherent Raman spectroscopy methods. For diatomic molecules, we were able to extract vibrational temperatures from the intensities of spectrally resolved vibrational hot bands. Larger molecules in the condensed phase suffer band broadening effects that obscure the vibrational hot bands as the anharmonicities are typically smaller than the widths of the bands. This inability to resolve the hot bands inhibits the extraction of vibrational temperature and the measurement of the vibrational frequency of the fundamental. Here, we use a hot band model based on gas-phase anharmonic coupling coefficients to fit coherent anti-Stokes Raman spectra of the ν 1 vibrational mode of shock compressed condensed phase N 2 O with shock pressures and temperatures estimated from literature equations of state and compare to fits from a model using a single Gaussian peak. We report the resulting vibrational frequency shifts with shock pressure.
Potential energy surfaces for high-energy collisions between an oxygen molecule and a nitrogen atom are useful for modeling chemical dynamics in shock waves. In the present work, we present doublet, quartet, and sextet potential energy surfaces that are suitable for studying collisions of O 2 ( 3 Σ$^{–}_{g}$) with N( 4 S) in the electronically adiabatic approximation. Two sets of surfaces are developed, one using neural networks (NNs) with permutationally invariant polynomials (PIPs) and one with the least-squares many-body (MB) method, where a two-body part is an accurate diatomic potential and the three-body part is expressed with connected PIPs in mixed-exponential-Gaussian bond order variables (MEGs). We find, using the same dataset for both fits, that the fitting performance of the PIP-NN method is significantly better than that of the MB-PIP-MEG method, even though the MB-PIP-MEG fit uses a higher-order PIP than those used in previous MB-PIP-MEG fits of related systems (such as N 4 and N 2 O 2 ). However, the evaluation of the PIP-NN fit in trajectory calculations requires about 5 times more computer time than is required for the MB-PIP-MEG fit.
Here, in this work, we provide a nuanced view of electron correlation in the context of transition metal complexes, reconciling computational characterization via spin and spatial symmetry breaking in single-reference methods with qualitative concepts from ligand-field and molecular orbital theories. These insights provide the tools to reliably diagnose the multi-reference character, and our analysis reveals that while strong (i.e., static) correlation can be found in linear molecules (e.g., diatomics) and weakly bound and antiferromagnetically coupled (monometal-noninnocent ligand or multi-metal) complexes, it is rarely found in the ground-states of mono-transition-metal complexes. This leads to a picture of static correlation that is no more complex for transition metals than it is, e.g., for organic biradicaloids. In contrast, the ability of organometallic species to form more complex interactions, involving both ligand-to-metal σ-donation and metal-to-ligand π-backdonation, places a larger burden on a theory's treatment of dynamic correlation. We hypothesize that chemical bonds in which inter-electron pair correlation is non-negligible cannot be adequately described by theories using MP2 correlation energies and indeed find large errors vs experiment for carbonyl-dissociation energies from double-hybrid density functionals. A theory's description of dynamic correlation (and to a less important extent, delocalization error), which affects relative spin-state energetics and thus spin symmetry breaking, is found to govern the efficacy of its use to diagnose static correlation.
We present a new, novel implementation of the Many-Body Expansion (MBE) to account for the breaking of covalent bonds, thus extending the range of applications from its previous popular usage in the breaking of hydrogen bonds in clusters to molecules. A central concept of the new implementation is the in situ atomic electronic state of an atom in a molecule that casts the one-body term as the energy required to promote it to that state from its ground state. The rest of the terms correspond to the individual diatomic, triatomic, etc., fragments. Its application to the atomization energies of the XH n series, X = C, Si, Ge, Sn and n = 1–4, suggests that the (negative, stabilizing) 2-B is by far the largest term in the MBE with the higher order terms oscillating between positive and negative values and decreasing dramatically in size with increasing rank of the expansion. The analysis offers an alternative explanation for the purported “first row anomaly” in the incremental H n–1 X–H bond energies seen when these energies are evaluated with respect to the lowest energy among the states of the XH n molecules. Due to the “flipping” of the ground/first excited state between CH 2 ( 3 B 1 ground state, 1 A 1 first excited state) and XH 2 , X = Si, Ge, Sn ( 1 A 1 ground state, 3 B 1 first excited state), the overall picture does not exhibit a “first row anomaly” when the incremental bond energies are evaluated with respect to the molecular states having the same in situ atomic states.
The fully correlated frequency-independent Dirac–Coulomb–Breit Hamiltonian provides the most accurate description of electron–electron interaction before going to a genuine relativistic quantum electrodynamics theory of many-electron systems. In this work, we introduce a correlated Dirac–Coulomb–Breit multiconfigurational self-consistent-field method within the frameworks of complete active space and density matrix renormalization group. In this approach, the Dirac–Coulomb–Breit Hamiltonian is included variationally in both the mean-field and correlated electron treatment. Here, we also analyze the importance of the Breit operator in electron correlation and the rotation between the positive- and negative-orbital space in the no-virtual-pair approximation. Atomic fine-structure splittings and lanthanide contraction in diatomic fluorides are used as benchmark studies to understand the contribution from the Breit correlation.
We consider an ensemble of diatomic molecules resonantly coupled to an optical cavity under strong coupling conditions at normal incidence. Photodissociation dynamics is examined via direct numerical integration of the coupled Maxwell–Schrödinger equations with molecular rovibrational degrees of freedom explicitly taken into account. It is shown that the dissociation is significantly affected (slowed down) when the system is driven at its polaritonic frequencies. The observed effect is demonstrated to be of transient nature and has no classical analog. An intuitive explanation of the dissociation slowdown at polaritonic frequencies is proposed.
Time-resolved optical emission spectroscopy is used to explore the early development of anode-initiated flashover in vacuum. An experimental apparatus for collecting flashover self-luminosity is introduced, which enables spatially resolved light collection from regions adjacent to the anode and cathode triple junctions. The emitted light is spectrally and temporally resolved utilizing an imaging spectrograph and a fast electronic shutter camera. The first light detectable from flashovers across both Rexolite and polytetrafluoroethylene (PTFE) insulators begins less than 10 ns prior to the flashover gap's impedance collapse and consists of faint broadband feature in primarily visible wavelengths. This suggests luminosity due to excitation in the solid, possibly occurring as a result of field emitted electrons. In the few nanoseconds leading up to and including the impedance collapse, the broadband spectrum grows in intensity and extends into the UV, eventually accompanied by a few emission lines of the insulator material. This is strong evidence that the early stages of anode-initiated flashover include surface layer breakdown of the insulator, as a contrast from cathode-initiated flashover, which is predominately an above-surface process. Spectra accumulated over longer exposures of PTFE flashovers indicate that, in the first few hundred nanoseconds after the impedance collapse, ions from the insulator material constitute the majority of the emission lines. Later, neutral and ionic metal species from the electrode contribute with similar prominence as well as molecular bands from diatomic carbon. Finally, a comparison is provided of the results of these studies to other spectroscopic investigations of vacuum flashover from the literature.
We present simulations of an oxygen gas puff Z-pinch on a University scale generator at the Weizmann Institute of Science. The work accounts for the detailed geometry of the nozzle, the initial neutral gas density distribution, and the subsequent implosion. The modeling results show significant improvement with data for the current at the time of stagnation in comparison with a previous effort [Rosenzweig et al., Phys. Plasmas 27, 022705 (2020)]. As a first step, we performed simulations of the flow of neutral diatomic oxygen from a plenum through a nozzle within a recessed cathode, across a gap, and into the anode with a recessed grounded honeycomb. These simulations show an agreement with the measured initial gas density profiles within the region not blocked by the recesses and accessible to visible measurements. The computed neutral gas flow profile serves as the initial condition for a radiation magnetohydrodynamic simulation of the implosion using the MACH2-TCRE code. By considering the specific details of the nozzle and chamber geometry, we find agreement with the measured current profile, including the inductive notch. The simulations predict that the plasma undergoes a strong pinch within the hidden anode recess. The simulations also predict the strongest radiation pulse occurs within the anode recess and at the time of the observed inductive notch.
A relativistic magnetic hyperfine interaction Hamiltonian based on the Douglas–Kroll–Hess (DKH) theory up to the second order is implemented within the ab initio multireference methods, including spin–orbit coupling in the Molcas/OpenMolcas package. This implementation is applied to calculate relativistic hyperfine coupling (HFC) parameters for atomic systems and diatomic radicals with valence s or d orbitals by systematically varying active space size in the restricted active space self-consistent field formalism with restricted active space state interaction for spin–orbit coupling. The DKH relativistic treatment of the hyperfine interaction reduces the Fermi contact contribution to the HFC due to the presence of kinetic factors that regularize the singularity of the Dirac delta function in the nonrelativistic Fermi contact operator. This effect is more prominent for heavier nuclei. As the active space size increases, the relativistic correction of the Fermi contact contribution converges well to the experimental data for light and moderately heavy nuclei. The relativistic correction, however, does not significantly affect the spin-dipole contribution to the hyperfine interaction. In addition to the atomic and molecular systems, the implementation is applied to calculate the relativistic HFC parameters for large trivalent and divalent Tb-based single-molecule magnets (SMMs), such as Tb(III)Pc2 and Tb(II)(CpiPr5)2 without ligand truncation using well-converged basis sets. In particular, for the divalent SMM, which has an unpaired valence 6s/5d hybrid orbital, the relativistic treatment of HFC is crucial for a proper description of the Fermi contact contribution. Even with the relativistic hyperfine Hamiltonian, the divalent SMM is shown to exhibit strong tunability of HFC via an external electric field (i.e., strong hyperfine Stark effect).
We report the detection of atomic oxygen and quantitative measurements of its electronic Raman coherence decays in flames and low-temperature plasmas using time-resolved hybrid femtosecond/picosecond coherent anti-Stokes Raman scattering (CARS). Atomic oxygen was detected using the Raman transitions between the spin–orbit coupled triplet ground states. Atomic oxygen was generated in an H 2 /O 2 /Ar diffusion flame and an O 2 /Ar pulsed plasma discharge. Single exponential decays were observed for the O( 3 P 2 )–O( 3 P 1 ) Raman transition at 158.3 cm −1 and the O( 3 P 2 )–O( 3 P 0 ) Raman transition at 227 cm −1 . From the decay measurements, the atomic O Raman linewidths were obtained from 25 to 150 Torr in non-equilibrium plasma and at 760 Torr in a flame. Enhanced signal-to-noise ratios (SNRs) of atomic oxygen and atomic to molecular oxygen signal contrasts were obtained by taking advantage of electronic triplet coherence beating. Enhancement of up to seven times in the atomic O SNR was observed. Furthermore, we also found that the dephasing rates of O 2 (v = 0–3, N = 37) were similar, which provides evidence for the assumption that vibrational excitation does not influence the dephasing of diatomic molecular rotational CARS transitions.
We present a variant of the approximate second order coupled-cluster method (CC2) with a two-parameter size-consistent Brillouin–Wigner (BW-s) partitioning instead of a Møller–Plesset (MP) partitioning for the unperturbed Hamiltonian, which we refer to as BWs-CC2. The computational complexity of this model scales identically to CC2 with molecular size. Conventional CC2 and its regularized BWs-CC2 variants, as well as conventional MP2 and two of its regularized BW-s2 variants, were assessed on a 535 element database spanning thermochemistry, non-covalent interactions, barrier heights, and isomerization energies. To ensure a well-defined model chemistry, the assessment was performed using internally stable spin-polarized Hartree–Fock (HF) orbitals in the finite aug-cc-pVQZ basis without counterpoise corrections. As a result of using stable orbitals, contrary to conventional wisdom, we find that CC2 substantially outperforms MP2 on molecules with significantly spin contaminated reference orbitals without a significant increase in error on systems with a spin-pure reference, showing the value of its single substitutions. While no single choice of regularization parameters can be optimal for all datasets, we find that BWs-CC2 generally outperforms both CC2 and BW-s2 with a single judicious parameter choice. Additional tests on dipole moments and bond lengths of diatomics provide further support for the utility of this choice. Furthermore, the main outliers and poorest performing cases are associated with large amounts of spin-contamination in the HF reference, which is indicative of systems with either strong correlation or extensive artificial symmetry breaking. Overall, these findings argue that the perception of the quality of the CC2 ground state should be reevaluated and that it can be further improved upon by the soundly based BWs-CC2 variant with the recommended parameter choice.
Accurate and efficient gradients of molecular energy with respect to nuclear degrees of freedom are essential for geometry optimization and molecular dynamics, including simulations that go beyond the Born–Oppenheimer regime. A common approach involves deriving analytical formulas for new electronic structure methods, which is often conceptually difficult and requires tedious coding. Here, we implement analytical, semi-numerical, and automatic differentiation (AD)-based gradient pathways for semiempirical Hamiltonian models in the PYSEQM software package, leveraging both graphics processing unit (GPU) and central processing unit (CPU) architectures. We further extend these capabilities to excited states calculated using the configuration interaction singles and time-dependent Hartree–Fock ansätze. We benchmark wall time, peak memory usage, and accuracy across three molecular families of varying chemical complexity, including systems of up to a thousand atoms. For ground-state simulations, analytical and AD gradients achieve near-identical GPU runtimes, while semi-numerical gradients are slower on GPU but remain competitive on CPU. For excited states, both analytical and custom AD approaches using implicit differentiation show similar performance and low memory requirements, whereas gradients with full AD are memory-limited. AD gradients match analytical ones in accuracy across all tested systems, aided by a quaternion-based diatomic frame rotation for two-center quantities that ensures smooth energy surfaces. Overall, automatic differentiation emerges as a practical alternative to analytical gradients in semiempirical quantum chemistry, offering high accuracy while allowing seamless integration in AI-driven workflows and popular packages, such as PyTorch and JAX. Our results provide actionable guidance for selecting optimal gradient strategies in large-scale ground- and excited-state molecular dynamics simulations.
Marine dissolved organic matter, which originates from phytoplankton, holds as much carbon as Earth’s atmosphere; yet, the biological processes governing its fate are primarily studied under idealized laboratory conditions or through indirect measures such as genome sequencing. In this work, we used isotope labeling to directly quantify uptake of complex carbon pools from the two primary sources of marine organic carbon (diatoms and cyanobacteria) by a natural microbial community. Furthermore, our data show that carbon pools are partitioned into distinct microbial lineages whose physiological properties and resource acquisition strategies match the chemical nature of their preferred substrates. Our results provide ecological and functional insights into the patterns of microbial community structure changes that occur during marine phytoplankton blooms.
The cyclopentadiene (C 5 H 6 ) molecule has emerged as a molecular building block of nonplanar polycyclic aromatic hydrocarbons (PAHs) and carbonaceous nanostructures such as corannulene (C 20 H 10 ), nanobowls (C 40 H 10 ), and fullerenes (C 60 ) in deep space. However, the underlying elementary gas-phase processes synthesizing cyclopentadiene from acyclic hydrocarbon precursors have remained elusive. Here, by merging crossed molecular beam experiments with rate coefficient calculations and comprehensive astrochemical modeling, we afford persuasive testimony on an unconventional low-temperature cyclization pathway to cyclopentadiene from acyclic precursors through the reaction of the simplest diatomic organic radical—methylidyne (CH)—with 1,3-butadiene (C 4 H 6 ) representing main route to cyclopentadiene observed in TaurusMolecular Cloud. This facile route provides potential solution for the incorporation of the cyclopentadiene moiety in complex aromatic systems via bottom–up molecular mass growth processes and offers an entry point to the low-temperature chemistry in deep space leading eventually to nonplanar PAHs in our carbonaceous Universe.
Microscale symbioses can be critical to ecosystem functions, but the mechanisms of these interactions in nature are often cryptic. Here, we use a combination of stable isotope imaging and tracing to reveal carbon (C) and nitrogen (N) exchanges among three symbiotic primary producers that fuel a salmon-bearing river food web. Bulk isotope analysis, nanoSIMS (secondary ion mass spectrometry) isotope imaging, and density centrifugation for quantitative stable isotope probing enabled quantification of organism-specific C- and N-fixation rates from the subcellular scale to the ecosystem. After winters with riverbed-scouring floods, the macroalga Cladophora glomerata uses nutrients in spring runoff to grow streamers up to 10 m long. During summer flow recession, riverine N concentrations wane and Cladophora becomes densely epiphytized by three species of Epithemia, diatoms with N-fixing endosymbionts (proto-organelles) descended from a free-living Crocosphaera cyanobacterium. Over summertime epiphyte succession on Cladophora, N-fixation rates increased as Epithemia spp. became dominant, Cladophora C-fixation declined to near zero, and Epithemia C-fixation increased. Carbon transfer to caddisflies grazing on Cladophora with high densities of Epithemia was 10-fold higher than C transfer to caddisflies grazing Cladophora with low Epithemia loads. In response to demand for N, Epithemia allocates high levels of newly fixed C to its endosymbiont. Consequently, these endosymbionts have the highest rates of C and N accumulation of any taxon in this tripartite symbiosis during the biologically productive season and can produce one of the highest areal rates of N-fixation reported in any river ecosystem.
The paper reports development of an atomic mean-field (AMF) spin–orbit approach within the spin-free exact two-component theory in its one-electron variant (SFX2C-1e), a new relativistic quantum-chemical method for treating spin-orbit coupling in molecules. The effective one-electron spin–orbit integral matrix in the four-component representation is assembled as a direct sum of one-centre spin–orbit integral matrices with the mean-field two-electron contributions evaluated using atomic SFX2C-1e Hartree–Fock density matrices. It is then transformed into two-component representation using analytic SFX2C-1e energy derivative formulation. The resulting two-component spin–orbit integral matrix is by design suitable for use in perturbative calculations of spin–orbit coupling, treating SFX2C-1e wavefunctions as unperturbed states. The accuracy of the present AMF approach has been demonstrated using benchmark calculations of spin–orbit splittings for representative diatomic radicals at the equation-of-motion coupled-cluster singles and doubles level. In conclusion, to demonstrate the applicability and accuracy of the present perturbative spin–orbit scheme in calculations of challenging heavy-element containing systems, a thorough computational investigation of six low-lying electronic states of ThO+ is reported.