Search NASA⌕ Search

SEARCH · Search NASA

Results for “diatomic molecule”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Revealing the Bacterial Quorum-Sensing Effect on the Biofilm Formation of Diatom Cylindrotheca sp. Using Multimodal Imaging

Diatoms contribute to carbon fixation in the oceans by photosynthesis and always form biofouling organized by extracellular polymeric substances (EPS) in the marine environment. Bacteria-produced quorum-sensing signal molecules N-acyl homoserine lactones (AHLs) were found to play an important role in the development of Cylindrotheca sp. in previous studies, but the EPS composition change was unclear. This study used the technology of alcian blue staining and scanning electron microscopy (SEM), confocal laser scanning microscopy (CLSM), and time-of-flight secondary ion mass spectrometry (ToF-SIMS) to directly observe the biofilm formation process. The results showed that AHLs promote the growth rates of diatoms and the EPS secretion of biofilm components. AHLs facilitated the diatom-biofilm formation by a forming process dependent on the length of carbon chains. AHLs increased the biofilm thickness and the fluorescence intensity and then altered the three-dimensional (3D) structures of the diatom-biofilm. In addition, the enhanced EPS content in the diatom-biofilm testified that AHLs aided biofilm formation. This study provides a collection of new experimental evidence of the interaction between bacteria and microalgae in fouling biofilms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Breaking covalent bonds in the context of the many-body expansion (MBE). I. The purported “first row anomaly” in XH n (X = C, Si, Ge, Sn; n = 1–4)

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.

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

Competing ionization and dissociation: Extension of the energy-dependent frame transformation to the gerade symmetry of H 2

This article solves two major tasks that frequently arise in the theory of electron collisions with a target molecular cation. First, it extends the energy-dependent frame transformation (EDFT) treatment, which is needed to map fixed-nuclei electron-molecule scattering matrices into an energy-dependent laboratory-frame scattering matrix with vibrational channel indices. The EDFT mapping can now be carried out even when the target molecule possesses multiple low-energy potential curves, significantly transcending previous applications. Second, it implements a method to extract the rest of the full laboratory-frame scattering matrix, i.e., the columns and rows describing input and/or output dissociation channels. The treatment is benchmarked in this article against the essentially exact solution of a refined two-dimensional model of the singlet gerade Σ symmetry of H 2 . Our tests demonstrate that the theory accurately maps fixed-nuclei scattering information, of the type provided by existing electron-molecule computer codes, into a laboratory-frame scattering matrix that includes both ionization and dissociation. Furthermore, this treatment can provide a general framework applicable to a broad class of electron collision processes involving diatomic target ions, suitable for an accurate description of challenging processes such as dissociative recombination.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Annotation of DOM metabolomes with an ultrahigh resolution mass spectrometry molecular formula library

Current approaches to analyzing metabolomic data often rely on matching MS/MS fragmentation data to sparse libraries or databases. This approach results in limited identification of features, often with less than 10% of the dataset being annotated. A complementary approach is to assign molecular formula to features based on accurate mass measurements, but the platforms commonly used for metabolomics do not have the needed accuracy or resolving power to do this robustly, particularly for larger molecules. Using our newly modified analysis tool, CoreMS, we generated a library of molecular formula from pooled samples analyzed with LC-21T FT-ICR MS. This library successfully annotated approximately 53.2% of features identified from the exometabolome of marine diatom Phaeodactylum tricornutum – a nearly ten-fold increase over the 5.9% annotation rate achieved using a conventional MS/MS library matching approach. Using this FT-ICR MS library approach, we were able to differentiate differences in the exometabolome of P. tricornutum in iron replete and iron limited conditions, with 668 metabolites being differentially expressed (p < 0.05, 2 x intensity difference) under these conditions. The traditional MS/MS fragmentation-based annotation approach only annotated 61 of these metabolites, while our novel pipeline annotated 450 metabolites and revealed 12 metabolites that were significantly more abundant under low iron conditions. Our results demonstrate the utility of ultrahigh resolution mass spectrometry for generating more comprehensive and confident molecular annotations.

21T-FTICR-MS, CoreMS↗
DOE OSTI2024

Bottom-up design of actinide materials from molecular clusters: Demonstration of a general-purpose simulation capability leveraging machine-learned atomic potentials

Actinide thin-film coatings such as uranium dioxide (UO 2 ) play an important role in nuclear reactors and other mission-relevant applications, but realization of their potential requires a deep fundamental understanding of the chemical vapor deposition (CVD) processes used for their growth. The slow experimental progress can be attributed, in part, to the standard safety guidelines associated with handling uranium byproducts, which are often corrosive, toxic, and radioactive. Accurate simulation techniques, when used in concert with experiment, can improve laboratory safety, material durability, and deliverable timeframes. However, state-of-the-art computational methods are either insufficiently accurate or intractably expensive. To remedy this situation, in this project we suggested a machine-learning (ML) accelerated workflow for simulating molecular clustering toward deposition. As a benchmark test case, we considered molecular clustering in steam and assessed independent components of our workflow by comparing with measured thermodynamic properties of water. After analyzing each component individually and finding no fundamental barrier to realization of the workflow, we attempted to integrate the ML component, a Sandia-developed tool called FitSNAP. As this was the first application of FitSNAP to atoms and molecules in the gas phase at Sandia, the method required more fitting data than was originally anticipated. Systematic improvements were made by including in the fit data diatomic potentials, molecular single-bond-breaking curves, and symmetry-constrained intermolecular potentials. We concluded that our strategy provides a feasible pathway toward modeling CVD and related processes, but that extensive training data must be generated before it can be of practical use.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A bloom of a single bacterium shapes the microbiome during outdoor diatom cultivation collapse

Algae-dominated ecosystems are fundamentally influenced by their microbiome. We lack information on the identity and function of bacteria that specialize in consuming algal-derived dissolved organic matter in high algal density ecosystems such as outdoor algal ponds used for biofuel production. Here, we describe the metagenomic and metaproteomic signatures of a single bacterial strain that bloomed during a population-wide crash of the diatom, Phaeodactylum tricornutum, grown in outdoor ponds. 16S rRNA gene data indicated that a single Kordia sp. strain (family Flavobacteriaceae) contributed up to 93% of the bacterial community during P. tricornutum demise. Kordia sp. expressed proteins linked to microbial antagonism and biopolymer breakdown, which likely contributed to its dominance over other microbial taxa during diatom demise. Analysis of accompanying downstream microbiota (primarily of the Rhodobacteraceae family) provided evidence that cross-feeding may be a pathway supporting microbial diversity during diatom demise. In situ and laboratory data with a different strain suggested that Kordia was a primary degrader of biopolymers during algal demise, and co-occurring Rhodobacteraceae exploited degradation molecules for carbon. An analysis of 30 Rhodobacteraceae metagenome assembled genomes suggested that algal pond Rhodobacteraceae commonly harbored pathways to use diverse carbon and energy sources, including carbon monoxide, which may have contributed to the prevalence of this taxonomic group within the ponds. These observations further constrain the roles of functionally distinct heterotrophic bacteria in algal microbiomes, demonstrating how a single dominant bacterium, specialized in processing senescing or dead algal biomass, shapes the microbial community of outdoor algal biofuel ponds.

Kordia↗
DOE OSTI2023

ORNL_AISD-Ex: Quantum chemical prediction of UV/Vis absorption spectra for over 10 million organic molecules

We performed calculations of electronic excitation energies and associated oscillator strengths based on the time-dependent density-functional tight-binding (TD-DFTB) method [1]. The SMILES (Simplified molecular-input line-entry system) strings of the molecules from the AISD HOMO-LUMO database [2] were converted to a 3D atomistic structure and stored in a PDB file after preliminary geometry optimization using the Merck Molecular Force Field (MMFF94) in RDKit [3,4]. The primary information stored in the PDB file archive consists of Cartesian coordinates for each atom of the molecule in their 3D location in space, along with summary information about the structure, sequence, and experiment. We then performed molecular geometry optimization using the density-functional tight-binding (DFTB) method [5] in the electronic ground state, followed by single-point excited states calculations, as described below. We note that, since RDKit employs a random choice for the generation of molecular conformers, the molecular geometries obtained in this dataset could be different from the ones that were generated when the AISD HOMO-LUMO dataset was generated. The computed excitation energies and associated oscillator strengths can be converted to predict UV/Vis absorption spectra, where excitation energies correspond to absorption peak positions, and oscillator strengths are a good measure of the probability of absorption of visible or UV light in transitions between electronic ground and excited states. The conversion of SMILES strings to 3D Cartesian coordinates of fully DFTB-optimized molecules was successful for 10,502,904 out of 10,502,917 molecules. For these molecules, both geometry optimizations and excited states calculations were successful. The DFTB calculations did not complete for 13 molecules of the original AISD HOMO-LUMO dataset. We still provide information about the geometry of these molecules. The molecules are diverse for chemical compositions (which span 5 non-hydrogen elements: oxygen, carbon, nitrogen, fluorine, sulfur) and molecular size (the smallest molecule contains 5 non-hydrogen atoms, and the largest molecule contains 71 non-hydrogen atoms). The DFTB method [5] is an approximation to density functional theory (DFT), utilizing a minimal basis set in conjunction with a two-center approximation to the electronic Hamiltonian and overlap matrix elements. The DFTB total energy is the sum of an electronic and a repulsive energy contribution, and their calculation requires optimized electronic parameters and diatomic repulsive potential energy functions. All DFTB calculations were performed using the DFTB+ code [6] (version 21.2) and the wrapper for DFTB+ in the Atomic Simulation Environment (ASE) (version 3.22.1) [7], which performed an internal conversion of Cartesian coordinates from PDB to the .gen file format. For the geometry optimizations on the electronic ground state potential energy surface of the molecules, we have chosen the third-order DFTB (DFTB3) method [5c] and employed the matching 3ob set of electronic parameters and repulsive potentials [8]. The empirical γ-damping for hydrogen bond correction, and Grimme's D3 empirical dispersion correction with Becke-Johnson damping (D3(BJ)) [9] dispersion correction was included to improve the description of non-covalent interactions. For excited states single-point energy calculations, we employed the TD-DFTB method in conjunction with the DFTB2 method [5b] and the matching mio [5b,10] and halorg [11] parameter sets. We opted to request the simultaneous calculation of 50 excited states for singlet transition to investigate sufficient number of excited states, based on linear response theory using the Casida equation [Ref: T. A. Niehaus, S. Suhai, F. Della Sala, P Lugli, M. Elstner, G. Seifert, and Th. Frauenheim. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 63:085108, 2001] and the ARPACK diagonalizer [R. B. Lehoucq, D. C. Sorensen, and C. Yang. Arpack users guide: Solution of large-scale eigenvalue problems by implicitly restarted arnoldi methods, 1997. 46, 51]. The dataset contains 1001 tar.gz files. Tar files are named as “ornl_aisd_ex_1.tar.gz†through “ornl_aisd_ex_1000.tar.gzâ€. Additionally, the 13 failed molecules are in “ornl_aisd_ex_unprocessed.tar.gzâ€. Except for the tar files listed below, each tar file contains 10,500 molecules. Tar files numbered 34, 121, 128, 352, 360, 429, 495, 509, 518, 627, 676, 668, and 862 contain 10,499 molecules each. The last tar file numbered 1000 contains 13,417 molecules. The total size of the uncompressed dataset is over 283 Gigabytes. The code for calculating the electronic excitation energies and statistical analysis of the dataset is provided at the following GitLab repository: https://github.com/ORNL/Analysis-of-Large-Scale-Molecular-Datasets-with-Python Calculating the UV spectrum of a molecule requires performing 3 main operations: 1. Converting the smiles string representation of a molecule into a geometric structure where each atom is assigned XYZ coordinates. The geometric structure is written to the file smiles.pdb. 2. Using smiles.pdb to compute the relaxed geometry of the molecule, which corresponds with the position of the atoms at the position of equilibrium at the ground state. This generates the files band.out, detailed.out, and geo_end.gen. 3. Using geo_end.gen to calculate the UV spectrum of the molecule which is written into the file EXC.DAT. Every molecule in the dataset has its own directory. The files contained in each molecule directory are as follows: 1. geo_end.gen 2. detailed.out 3. band.out 4. EXC.DAT 5. smiles.pdb REFERENCES [1] Niehaus, T. A.; Suhai, S.; Della Salla, F.; Lugli, P.; Elstner, M.; Seifert, G.; Frauenheim, Th. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 2001, 63, 085108/1-9. [2] Blanchard, A.; Gounley, J.; Metha, K.; Yoo, P.; Irle, S. AISD HOMO-LUMO. DOI: 10.13139/ORNLNCCS/1869409 [3] RDKit: Cheminformatics and Machine Learning Software. 2013, [http://www.rdkit.org] [4] Tosco, P.; Stiefl, N. and Landrum, G. Bringing the MMFF force field to the RDKit: implementation and validation. J Cheminform. 2014, 6, 1–4. [5] a) Porezag, D.; Frauenheim, T.; Kohler, T.; Seifert, G.; Kaschner, Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon, R. Phys. Rev. B 1995, 51, 12947-12957; b) Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, Th.; Suhai, S.; Seifert, G.; Phys. Rev. B 1998, 58, 7260-7268; c) Gaus, M.; Cui, Q.; Elstner, M. DFTB3: Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method (SCC-DFTB), J. Chem. Theory Comput. 2011, 7, 931-948; d) Cui, Q.; Elstner, M. Density functional tight binding: values of semi-empirical methods in an ab initio era, Phys. Chem. Chem. Phys. 2014, 16, 14368-14377. [6] Hourahine, B. et al. DFTB+, a software package for efficient approximate density functional theory based atomistic simulations, J. Chem. Phys. 2020, 152, 124101/1-19. [7] Larsen, A. H. et al. The atomic simulation environment—a Python library for working with atoms. J. Phys.: Cond. Matter 2017, 29, 273002. [8] Kubillus, M.; Kubar, T.; Gaus, M.; Rezac, J.; Elstner, M. Parameterization of the DFTB3 Method for Br, Ca, Cl, F, I, K, and Na in Organic and Biological Systems, J. Chem. Theory Comput. 2015, 11, 332-342. [9] Brandenburg, J. G.; Grimme, S. Accurate Modeling of Organic Molecular Crystals by Dispersion-Corrected Density Functional Tight Binding (DFTB), J. Phys. Chem. Lett. 2014, 5, 1785−1789. [10] a) Niehaus, T. A.; Elstner, M.; Frauenheim, Th.; Suhai, S. Application of an approximate density-functional method to sulfur containing compounds. J. Mol. Struct.: THEOCHEM 2001, 541, 185-94; b) Elstner, M.; Hobza, P.; Frauenheim, Th.; Suhai, S.; Kaxiras, E. Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment. J. Chem. Phys. 2001, 114, 5149-55. [11] Kubar, T.; Bodrog, Z.; Gaus, M.; Köhler, C.; Aradi, B.; Frauenheim, Th.; Elstner, M. Parametrization of the SCC-DFTB Method for Halogens. J. Chem. Theory Comput. 2013, 9, 2939-49.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Evolutionary diversity of proton and water channels on the oxidizing side of photosystem II and their relevance to function

One of the reasons for the high efficiency and selectivity of biological catalysts arise from their ability to control the pathways of substrates and products using protein channels, and by modulating the transport in the channels using the interaction with the protein residues and the water/hydrogen-bonding network. Here, this process is clearly demonstrated in Photosystem II (PS II), where its light-driven water oxidation reaction catalyzed by the Mn 4 CaO 5 cluster occurs deep inside the protein complex and thus requires the transport of two water molecules to and four protons from the metal center to the bulk water. Based on the recent advances in structural studies of PS II from X-ray crystallography and cryo-electron microscopy, in this review we compare the channels that have been proposed to facilitate this mass transport in cyanobacteria, red and green algae, diatoms, and higher plants. The three major channels (O1, O4, and Cl1 channels) are present in all species investigated; however, some differences exist in the reported structures that arise from the different composition and arrangement of membrane extrinsic subunits between the species. Among the three channels, the Cl1 channel, including the proton gate, is the most conserved among all photosynthetic species. We also found at least one branch for the O1 channel in all organisms, extending all the way from Ca/O1 via the ‘water wheel’ to the lumen. However, the extending path after the water wheel varies between most species. The O4 channel is, like the Cl1 channel, highly conserved among all species while having different orientations at the end of the path near the bulk. The comparison suggests that the previously proposed functionality of the channels in T. vestitus (Ibrahim et al., Proc Natl Acad Sci USA 117:12624–12635, 2020; Hussein et al., Nat Commun 12:6531, 2021) is conserved through the species, i.e. the O1-like channel is used for substrate water intake, and the tighter Cl1 and O4 channels for proton release. The comparison does not eliminate the potential role of O4 channel as a water intake channel. However, the highly ordered hydrogen-bonded water wire connected to the Mn 4 CaO 5 cluster via the O4 may strongly suggest that it functions in proton release, especially during the S 0 → S 1 transition (Saito et al., Nat Commun 6:8488, 2015; Kern et al., Nature 563:421–425, 2018; Ibrahim et al., Proc Natl Acad Sci USA 117:12624–12635, 2020; Sakashita et al., Phys Chem Chem Phys 22:15831–15841, 2020; Hussein et al., Nat Commun 12:6531, 2021).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Direction-specific enhanced diffusion of CO 2 in chiral hexagonal boron nitride nanotubes

To meet performance requirements, the next generation of gas separation membranes will need both high gas permeability and selectivity, attainable if we could coax adsorbates to minimize random Brownian motion and produce direction-specific diffusion along a desired axis. In this atomistic modeling study, we detail how direction-specific diffusion of CO 2 can be achieved in chiral hexagonal boron nitride nanotubes (hBNNTs) by means of a non-Knudsen diffusion mechanism. Our findings detail how this mechanism of diffusion is driven by interactions with the tube walls and enables the CO 2 molecules to diffuse along the nanotube’s z-axis with minimized collisions and directional changes. hBNNTs with chiral indices exhibit CO 2 diffusion rates faster than non-chiral tubes of comparable and larger diameters. Of the hBNNTs studied, a (7,3) tube appears to be ideally sized (3.7 Å radius) exhibiting CO 2 diffusion that is 3.4 times faster than diatomic N 2 . Applying this mechanism of diffusion to hypothetical sheet membranes prepared with aligned chiral (7,3) hBNNTs results in membranes with a calculated CO 2 /N 2 permselectivity of 170 and a CO 2 permeability limit of nearly 1.35 ×10 7 Barrer, readily surpassing the Robeson upper bound for CO 2 /N 2 separations.

CO2↗
DOE OSTI2021

Potential energy surfaces for high-energy N + O 2 collisions

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.

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

Single-cell isotope tracing reveals functional guilds of bacteria associated with the diatom Phaeodactylum tricornutum

Abstract Bacterial remineralization of algal organic matter fuels algal growth but is rarely quantified. Consequently, we cannot currently predict whether some bacterial taxa may provide more remineralized nutrients to algae than others. Here, we quantified bacterial incorporation of algal-derived complex dissolved organic carbon and nitrogen and algal incorporation of remineralized carbon and nitrogen in fifteen bacterial co-cultures growing with the diatom Phaeodactylum tricornutum at the single-cell level using isotope tracing and nanoSIMS. We found unexpected strain-to-strain and cell-to-cell variability in net carbon and nitrogen incorporation, including non-ubiquitous complex organic nitrogen utilization and remineralization. We used these data to identify three distinct functional guilds of metabolic interactions, which we termed macromolecule remineralizers, macromolecule users, and small-molecule users, the latter exhibiting efficient growth under low carbon availability. The functional guilds were not linked to phylogeny and could not be elucidated strictly from metabolic capacity as predicted by comparative genomics, highlighting the need for direct activity-based measurements in ecological studies of microbial metabolic interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Dynamic Phaeodactylum tricornutum exometabolites shape surrounding bacterial communities

Summary Roles of different ecological classes of algal exometabolites in regulating microbial community composition are not well understood. Here, we identify exometabolites from the model diatom Phaeodactylum tricornutum and demonstrate their potential to influence bacterial abundances. We profiled exometabolites across a time course of axenic algal growth using liquid chromatography–tandem mass spectrometry. We then investigated growth of 12 bacterial isolates on individual‐identified exometabolites. Lastly, we compared responses of a P. tricornutum ‐adapted enrichment community to additions of two contrasting metabolites: selective growth substrate 4‐hydroxybenzoic acid and putative signaling/facilitator molecule lumichrome. We identified 50 P. tricornutum metabolites and found distinct temporal accumulation patterns. Two exometabolites (of 12 tested) supported growth of distinct subsets of bacterial isolates. While algal exudates and algal presence drove similar changes in community composition compared with controls, exogenous 4‐hydroxybenzoic acid addition promoted increased abundances of taxa that utilized it in isolation, and also revealed the importance of factors relating to algal presence in regulating community composition. This work demonstrates that secretion of selective bacterial growth substrates represents one mechanism by which algal exometabolites can influence bacterial community composition and illustrates how the algal exometabolome has the potential to modulate bacterial communities as a function of algal growth.

4-hydroxybenzoic acid↗
DOE OSTI2025

The approximate second order coupled-cluster method based on a size-consistent Brillouin–Wigner partitioning

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.

Correlation energy↗
DOE OSTI2020

Automated Construction of Potential Energy Surfaces Suitable to Describe van der Waals Complexes with Highly Excited Nascent Molecules: The Rotational Spectra of Ar–CS( v ) and Ar–SiS( v )

Some reactions produce extremely hot nascent products which nevertheless can form sufficiently long-lived van der Waals (vdW) complexes—with atoms or molecules from a bath gas—as to be observed via microwave spectroscopy. Theoretical calculations of such unbound resonance states can be much more challenging than ordinary bound-state calculations depending on the approach employed. One encounters not just the floppy, and perhaps multiwelled potential energy surface (PES) characteristic of vdWs complexes, but in addition, one must contend with excitation of the intramolecular modes and its corresponding influence on the PES. Straightforward computation of the (resonance) rovibrational levels of interest, involves the added complication of the unbound nature of the wave function, often treated with techniques such as introducing a complex absorbing potential. Here, we have demonstrated that a simplified approach of making a series of vibrationally effective PESs for the intermolecular coordinates—one for each reaction product vibrational quantum number of interest—can produce vdW levels for the complex with spectroscopic accuracy. This requires constructing a series of appropriately weighted lower-dimensional PESs for which we use our freely available PES-fitting code AUTOSURF. The applications of this study are the Ar–CS and Ar–SiS complexes, which are isovalent to Ar–CO and Ar–SiO, the latter of which we considered in a previously reported study. Using a series of vibrationally effective PESs, rovibrational levels and predicted microwave transition frequencies for both complexes were computed variationally. A series of shifting rotational transition frequencies were also computed as a function of the diatom vibrational quantum number. The predicted transitions were used to guide and inform an experimental effort to make complementary observations. Comparisons are given for the transitions that are within the range of the spectrometer and were successfully recorded. Calculations of the rovibrational level pattern agree to within 0.2% with experimental measurements.

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

Low-temperature gas-phase formation of cyclopentadiene and its role in the formation of aromatics in the interstellar medium

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.

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

Performance of an atomic mean-field spin–orbit approach within exact two-component theory for perturbative treatment of spin–orbit coupling

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.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Nuclear cusps and singularities in the nonadditive kinetic potential bifunctional from analytical inversion

The non-additive kinetic potential $v^{\text{NAD}}$ is a key quantity in density-functional theory (DFT) embedding methods, such as frozen density embedding theory and partition DFT. $v^{\text{NAD}}$ is a bi-functional of electron densities $\rho_{\rm B}$ and $\rho_{\rm tot} = \rho_{\rm A} + \rho_{\rm B}$. It can be evaluated using approximate kinetic-energy functionals, but accurate approximations are challenging. The behavior of $v^{\text{NAD}}$ in the vicinity of the nuclei has long been questioned, and singularities were seen in some approximate calculations. In this article, the existence of singularities in $v^{\text{NAD}}$ is analyzed analytically for various choices of $\rho_{\rm B}$ and $\rho_{\rm tot}$, using the nuclear cusp conditions for the density and Kohn-Sham potential. It is shown that no singularities arise from smoothly partitioned ground-state Kohn-Sham densities. We confirm this result by numerical calculations on diatomic test systems HeHe, HeLi$^+$, and H$_2$, using analytical inversion to obtain a numerically exact $v^{\rm NAD}$ for the local density approximation. We examine features of $v^{\rm NAD}$ which can be used for development and testing of approximations to $v^{\rm NAD}[\rho_{\rm B},\rho_{\rm tot}]$ and kinetic-energy functionals.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Nonequilibrium Solvent Effects during Photodissociation in Liquids: Dynamical Energy Surfaces, Caging and Chemical Identity

In the gas phase, potential energy surfaces can be used to provide insight into the details of photochemical reaction dynamics. In solution, however, it is unclear what potential energy surfaces, if any, can be used to describe even simple chemical reactions such as the photodissociation of a diatomic solute. In this paper, we use mixed quantum/classical (MQC) molecular dynamics (MD) to study the photodissociation of $Na$ $^+_2$ in both liquid Ar and liquid tetrahydrofuran (THF). We examine both the gas-phase potential surfaces and potentials of mean force (PMF), which assume that the solvent remains at equilibrium with the solute throughout the photodissociation process and show that neither resemble a nonequilibrium dynamical energy surface that is generated by taking the time integral of work. For the photodissociation of $Na$ $^+_2$ in liquid Ar, the dynamical energy surface shows clear signatures of solvent caging, and the degree of caging is directly related to the mass of the solvent atoms. For $Na$ $^+_2$ in liquid THF, local specific interactions between the solute and solvent lead to changes in chemical identity that create a kinetic trap that effectively prevents the molecule from dissociating. Here, the results show that nonequilibrium effects play an important role even in simple solution-phase reactions, requiring the use of dynamical energy surface to understand such chemical events.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗