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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1972

Impulsive model for reactive collisions

A simple classical mechanical model of the reactive scattering of a structureless atom A and a quasi-diatomic BC is developed which takes full advantage of energy, linear and angular momentum conservation relations but introduces a minimum of further assumptions. These are as follows: (1) the vibrational degree of freedom of the reactant (BC) and product (AB) molecules is suppressed, so the change in vibrational energy is simply a parameter; (2) straight-line trajectories are assumed outside of a reaction shell; (3) within this zone, momentum transfer occurs impulsively (essentially instantaneously) following mass transfer; (4) the impulse, which may be either positive or negative, is directed along the BC axis, which may, however, assume all orientations with respect to the incident relative velocity. The model yields differential and total cross sections and product rotational energy distributions for a given collision exoergicity Q, or for any known distribution over Q. Numerical results are presented for several prototype reactions whose dynamics have been well-studied.

Marron, M. T.↗
NASA NTRS1986

Theoretical dissociation energies for ionic molecules

Ab initio calculations at the self-consistent-field and singles plus doubles configuration-interaction level are used to determine accurate spectroscopic parameters for most of the alkali and alkaline-earth fluorides, chlorides, oxides, sulfides, hydroxides, and isocyanides. Numerical Hartree-Fock (NHF) calculations are performed on selected systems to ensure that the extended Slater basis sets employed for the diatomic systems are near the Hartree-Fock limit. Extended Gaussian basis sets of at least triple-zeta plus double polarization equality are employed for the triatomic system. With this model, correlation effects are relatively small, but invariably increase the theoretical dissociation energies. The importance of correlating the electrons on both the anion and the metal is discussed. The theoretical dissociation energies are critically compared with the literature to rule out disparate experimental values. Theoretical (sup 2)Pi - (sup 2)Sigma (sup +) energy separations are presented for the alkali oxides and sulfides.

Langhoff, S. R.↗
NASA NTRS2001

Neutral Gas Temperature Estimates in an Inductively Coupled CF4 Plasma by Fitting Diatomic Emission Spectra

This work examines the accuracy of plasma neutral temperature estimates by fitting the rotational band envelope of different diatomic species in emission. Experiments are performed in an inductively coupled CF4 plasma generated in a Gaseous Electronics Conference reference cell. Visible and ultraviolet emission spectra are collected at a power of 300 W (approximately 0.7 W/cc) and pressure of 30 mtorr. The emission bands of several molecules (CF, CN, C2, CO, and SiF) are fit simultaneously for rotational and vibrational temperatures and compared. Four different rotational temperatures are obtained: 1250 K for CF and CN, 1600 K for CO, 1800 K for C2, and 2300 K for SiF. The vibrational temperatures obtained vary from 1750-5950 K, with the higher vibrational temperatures generally corresponding to the lower rotational temperatures. These results suggest that the different species have achieved different degrees of equilibration between the rotational and vibrational modes and may not be equilibrated with the translational temperatures. The different temperatures are also related to the likelihood that the species are produced by ion bombardment of the surface, with etch products like SiF, CO, and C2 having higher temperatures than species expected to have formed in the gas phase.

Cruden, Brett A.↗
NASA NTRS1986

Theoretical studies of diatomic and triatomic systems containing the group IB atoms Cu, Ag, and Au

Selected portions of the ground state potential energy surfaces of the Cu3, Ag3, AgCu2, and AuCu2 trimers are studied at the single-reference singles plus doubles configuration interaction and couple pair functional levels correlating 33 electrons. The calculations use the effective core potentials of Hay and Wadt (1985) to replace the deep core levels. The Cu3 and Ag3 molecules are found to have 2B2 obtuse-angled ground states, with low-lying 2A1 acute-angled excited states. The AgCu2 and AuCu2 molecules have 2A1 acute-angled ground states. The Cu3 molecule has a smaller 3d population than Ag3, and Cu3 has a smaller s electron density in the open-shell orbital than does Ag3, in agreement with recent ESR experiments.

Walch, Stephen P.↗
NASA NTRS1991

Theoretical studies of molecular collisions

The following subject areas are covered: (1) total integral reactive cross sections and vibrationally resolved reaction probabilities for F + H2 = HF + H; (2) a theoretical study of inelastic O + N2 collisions; (3) body frame close coupling wave packet approach to gas phase atom-rigit rotor inelastic collisions; (4) wave packet study of gas phase atom-rigit motor scattering; (5) the application of optical potentials for reactive scattering; (6) time dependent, three dimensional body frame quantal wave packet treatment of the H + H2 exchange reaction; (7) a time dependent wave packet approach to atom-diatom reactive collision probabilities; (8) time dependent wave packet for the complete determination of s-matrix elements for reactive molecular collisions in three dimensions; (9) a comparison of three time dependent wave packet methods for calculating electron-atom elastic scattering cross sections; and (10) a numerically exact full wave packet approach to molecule-surface scattering.

Kouri, Donald J.↗
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↗
NASA NTRS2010

Comparative Spectra of Oxygen-Rich Versus Carbon-Rich Circumstellar Shells: VY Canis Majoris and IRC(plus)10216 at 215-285 GHz

A sensitive (1sigma rms at 1 MHz resolution approx.3 mK) 1 mm spectral line survey (214.5-285.5 GHz) of VY Canis Majoris (VY CMa) and IRC +10216 has been conducted to compare the chemistries of oxygen- and carbon-rich circumstellar envelopes. This study was carried out using the Submillimeter Telescope of the Arizona Radio Observatory with a new Atacama Large Millimeter Array type receiver. This survey is the first to chemically characterize an O-rich circumstellar shell at millimeter wavelengths. In VY CMa, 128 emission features were detected arising from 18 different molecules; and in IRC +10216, 720 lines were observed, assigned to 32 different species. The 1 mm spectrum of VY CMa is dominated by SO, and SiS; in IRC +10216, C4H and SiC2 are the most recurrent species. Ten molecules were common to both sources: CO, SiS, SiO, CS, CN, HCN, HNC, NaCl, PN, and HCO(+). Sulfur plays an important role in VY CMa, but saturated/ unsaturated carbon dominates the molecular content of IRC +102.16, producing CH2NH, for example. Although the molecular complexity of IRC +10216 is greater, VY CMa supports a unique "inorganic" chemistry leading to the oxides PO, AlO, and AlOH. Only diatomic and triatomic compounds were observed in VY CMa, while species with four or more atoms are common in IRC +10216, reflecting carbon's ability to form multiple strong bonds, unlike oxygen. In VY CMa, a new water maser (v2 = 2) has been found, as well as vibrationally excited NaCl. Toward IRC +10216, vibrationally excited CCH was detected for the first time.

Tenebaum, E. D.↗
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↗
NASA NTRS1988

Quantum mechanical algebraic variational methods for inelastic and reactive molecular collisions

The quantum mechanical problem of reactive or nonreactive scattering of atoms and molecules is formulated in terms of square-integrable basis sets with variational expressions for the reactance matrix. Several formulations involving expansions of the wave function (the Schwinger variational principle) or amplitude density (a generalization of the Newton variational principle), single-channel or multichannel distortion potentials, and primitive or contracted basis functions are presented and tested. The test results, for inelastic and reactive atom-diatom collisions, suggest that the methods may be useful for a variety of collision calculations and may allow the accurate quantal treatment of systems for which other available methods would be prohibitively expensive.

Schwenke, David W.↗
NASA NTRS1997

Extraterrestrial Helium (He@C60) Trapped in Fullerenes in the Sudbury Impact Structure

Fullerenes (C60 and C70) have recently been identified in a shock-produced breccia (Onaping Formation) associated with the 1.85-Ga Sudbury Impact Crater. The presence of parts-per-million levels of fullerenes in this impact structure raises interesting questions about the processes that led to the formation of fullerenes and the potential for delivery of intact organic material to the Earth by a large bolide (e.g., asteroid or comet). Two possible scenarios for the presence of fullerenes in the Sudbury impact deposits are that (1) fullerenes are synthesized within the impact plume from the C contained in the bolide; or (2) fullerenes are already present in the bolide and survived the impact event. The correlation of C and trapped noble gas atoms in meteorites is well established. Primitive meteorites contain several trapped noble gas components that have anomalous isotopic compositions, some of which may have a presolar origin. Several C-bearing phases, including SiC, graphite, and diamond, have been recognized as carriers of trapped noble gases. It has also been suggested that fullerenes (C60 and C70) might be a carrier of noble gas components in carbonaceous chondrites. Recently, fullerenes have been detected in separate samples in the Allende meteorite. Carbon-60 is large enough to enclose the noble gases He, Ne, Ar, Kr, and Xe, but it is too small to contain diatomic gases such as N2 or triatomic gases such as CO2. Recent experimental work has demonstrated that noble gases of a specific isotopic composition can be introduced into synthetic fullerenes at high temperatures and pressures; these encapsulated gases can then be released by the breaking of one or more C bonds during step-heating under vacuum. These thermal-release patterns for He encapsulated within the C60 molecule (He@C60) are similar to the patterns for acid residues of carbonaceous chondrites, suggesting that fullerenes could be an additional carrier of trapped noble gases in acid residues of meteorites. Analysis and Results: In order to characterize the noble gas compositions of the Sudbury fullerenes, we undertook a systematic study of acid-resistant residues throughout the C-rich layer (Black member) of the Onaping Formation. Samples were demineralized and extracted using standard techniques. The Onaping extracts were analyzed using several techniques, including UV-Vis adsorption, electro spray mass spectrometry, and laser desorption (linear and reflectron) time-of-flight (TOF) mass spectrometry (LDMS). The Sudbury fullerenes were then separated and purified using HPLC coupled with a photo diode array detector. The HPLC extracts containing the purified fullerenes were loaded into a metal tube furnace within a glove box under a N atmosphere in preparation for noble gas analyses. The 3-He and 4-He content of the fullerene extracts was measured using previously reported standard techniques . Discussion: Fullerenes (C60 and C70) in the Sudbury Impact Structure have been found to contain trapped He with a 3-He/4-He ratio greater than 5 x 10(exp -4). The 3-He/4-He ratio exceeds the accepted solar value by more than 30% and is more than 10x higher than the maximum reported mantle value. Terrestrial nuclear reactions or cosmic-my bombardment are not sufficient to generate such a high ratio. The 3-He/4-He ratios in the Sudbury fullerenes are similar to those determined for interplanetary dust particles. The greater-than-solar ratios of 3-He/4-He in the Sudbury fullerenes may indicate a presolar origin, although alternative mechanisms occurring in the ISM to explain these high ratios (e.g., spallation reactions, selective He implantation, etc.) cannot be entirely ruled out. We are currently attempting to isolate enough fullerene material to measure anomalous Ne (or Kr or Xe) contained within the C60 (e.g., the "pure" 22-Ne component) and thus determine whether the Sudbury fullerenes are indeed presolar in origin.

Becker, L.↗
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↗
NASA NTRS1973

Equations for the kinetic modeling of supersonically flowing electrically excited lasers

The equations for the kinetic modeling of a supersonically flowing electrically excited laser system are presented. The work focuses on the use of diatomic gases, in particular carbon monoxide mixtures. The equations presented include the vibrational rate equation which describes the vibrational population distribution, the electron, ion and electronic level rate equations, the gasdynamic equations for an ionized gas in the presence of an applied electric field, and the free electron Boltzmann equation including flow and gradient coupling terms. The model developed accounts for vibration-vibration collisions, vibration-translation collisions, electron-molecule inelastic excitation and superelastic de-excitation collisions, charge particle collisions, ionization and three body recombination collisions, elastic collisions, and radiative decay, all of which take place in such a system. A simplified form of the free electron Boltzmann equation is developed and discussed with emphasis placed on its coupling with the supersonic flow. A brief description of a possible solution procedure for the set of coupled equations is then discussed.

Lind, R. C.↗
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↗
NASA NTRS1988

Electron-temperature dependence of the recombination of electrons with O4(+) ions

The coefficient for recombination of O4(+) ions with electrons has been measured as a function of Te by analyzing electron- and ion-density decays in microwave-generated afterglow plasmas in neon-oxygen mixtures. Te values in the afterglow were elevated above those of the neutral species and ions by continuous application of microwave power. The data analysis included effects on the Te arising from inelastic collisions of electrons with the minority oxygen molecules and spatial nonuniformities of Te due to the spatially nonuniform microwave heating field. Over the temperature range Te = 143-5500 K the recombination coefficients were found to have a Te exp -z dependence, with z = 0.48 + 0.1 or - 0.05. While the values of the recombination coefficients for this dimer ion are an order of magnitude larger than those of the monomer ions, the temperature dependence is close to the Te exp -0.5 dependence predicted and often observed for diatomic ions.

Dulaney, J. L.↗