Search NASA⌕ Search

SEARCH · Search NASA

Results for “calculating molecular abundances”

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 37 records · Page 2

Multi-scale characterization and simulation of impact welding between immiscible Mg/steel alloys

Vaporizing foil actuator spot welding method is used in this paper to join magnesium alloy AZ31 and uncoated high-strength steel DP590, which are typically considered as un-weldable due to their high physical property disparities, low mutual solubility, and the lack of any intermetallic phases. Characterization results from scanning electron microscopy (SEM) and high-resolution transmission electron microscopy (HRTEM) of the weld interface indicate that the impact creates an Mg nanocrystalline interlayer with abundant Fe particles. The interlayer exhibits intact bonding with both DP590 and AZ31 substrates. To investigate the fundamental bond formation mechanisms at the interface, a finite element (FE)-based process simulation is first performed to calculate the local temperature and deformation at the interface under the given macroscopic experimental condition. Finally, taking the FE results at the interface as inputs, molecular dynamics (MD) simulations are conducted to study the interlayer formation at the Mg/Fe interface during the impact and cooling. The results found a high velocity shearing-induced mechanical mixing mechanism that mixes Mg/Fe atoms at the interface and creates the interlayer, leading to the metallurgical bond between Mg/steel alloys.

36 MATERIALS SCIENCE↗

r-Process Radioisotopes from Near-Earth Supernovae and Kilonovae

The astrophysical sites where r-process elements are synthesized remain mysterious: it is clear that neutron star mergers (kilonovae (KNe)) contribute, and some classes of core-collapse supernovae (SNe) are also possible sources of at least the lighter r-process species. The discovery of 60Fe on the Earth and Moon implies that one or more astrophysical explosions have occurred near the Earth within the last few million years, probably SNe. Intriguingly, 244 Pu has now been detected, mostly overlapping with 60 Fe pulses. However, the 244 Pu flux may extend to before 12 Myr ago, pointing to a different origin. Here, motivated by these observations and difficulties for r-process nucleosynthesis in SN models, we propose that ejecta from a KN enriched the giant molecular cloud that gave rise to the Local Bubble, where the Sun resides. Accelerator mass spectrometry (AMS) measurements of 244 Pu and searches for other live isotopes could probe the origins of the r-process and the history of the solar neighborhood, including triggers for mass extinctions, e.g., that at the end of the Devonian epoch, motivating the calculations of the abundances of live r-process radioisotopes produced in SNe and KNe that we present here. Given the presence of 244 Pu, other r-process species such as 93 Zr, 107 Pd, 129 I, 135 Cs, 182 Hf, 236 U, 237 Np, and 247 Cm should be present. Their abundances and well-resolved time histories could distinguish between the SN and KN scenarios, and we discuss prospects for their detection in deep-ocean deposits and the lunar regolith. We show that AMS 129 I measurements in Fe–Mn crusts already constrain a possible nearby KN scenario.

79 ASTRONOMY AND ASTROPHYSICS↗

SOMAS: a platform for data-driven material discovery in redox flow battery development

Abstract Aqueous organic redox flow batteries offer an environmentally benign, tunable, and safe route to large-scale energy storage. The energy density is one of the key performance parameters of organic redox flow batteries, which critically depends on the solubility of the redox-active molecule in water. Prediction of aqueous solubility remains a challenge in chemistry. Recently, machine learning models have been developed for molecular properties prediction in chemistry and material science. The fidelity of a machine learning model critically depends on the diversity, accuracy, and abundancy of the training datasets. We build a comprehensive open access organic molecular database “Solubility of Organic Molecules in Aqueous Solution” (SOMAS) containing about 12,000 molecules that covers wider chemical and solubility regimes suitable for aqueous organic redox flow battery development efforts. In addition to experimental solubility, we also provide eight distinctive quantum descriptors including optimized geometry derived from high-throughput density functional theory calculations along with six molecular descriptors for each molecule. SOMAS builds a critical foundation for future efforts in artificial intelligence-based solubility prediction models.

25 ENERGY STORAGE↗

Molecular structure determination: Equilibrium structure of pyrimidine (m-C 4 H 4 N 2 ) from rotational spectroscopy (reSE) and high-level ab initio calculation (re) agree within the uncertainty of experimental measurement

The pure rotational spectrum of pyrimidine (m-C 4 H 4 N 2 ), the meta-substituted dinitrogen analog of benzene, has been studied in the millimeter-wave region from 235 GHz to 360 GHz. The rotational spectrum of the ground vibrational state has been assigned and fit to yield accurate rotational and distortion constants. Over 1700 distinct transitions were identified for the normal isotopologue in its ground vibrational state and least-squares fit to a partial sextic S-reduced Hamiltonian. Transitions for all four singly substituted 13 C and 15 N isotopologues were observed at natural abundance and were likewise fit. Deuterium-enriched samples of pyrimidine were synthesized, giving access to all eleven possible deuterium-substituted isotopologues, ten of which were previously unreported. Experimental values of rotational constants and computed values of vibration–rotation interaction constants and electron-mass corrections were used to determine semi-experimental equilibrium structures (r e SE ) of pyrimidine. The r e SE structure obtained using coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] corrections shows exceptional agreement with the re structure computed at the CCSD(T)/cc-pCV5Z level (≤0.0002 Å in bond distance and ≤0.03° in bond angle). Of the various computational methods examined, CCSD(T)/cc-pCV5Z is the only method for which the computed value of each geometric parameter lies within the statistical experimental uncertainty (2σ) of the corresponding semi-experimental coordinate. The exceptionally high accuracy and precision of the structure determination is a consequence of the large number of isotopologues measured, the precision and extent of the experimental frequency measurements, and the sophisticated theoretical treatment of the effects of vibration–rotation coupling and electron mass. Taken together, these demanding experimental and computational studies establish the capabilities of modern structural analysis for a prototypical monocyclic aromatic compound.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Principal Hugoniot of Iron-Bearing Olivine to 1465 GPa

Shock compression experiments on natural compositions are imperative to accurately model planetary accretion and the interior dynamics of planets. Combining shock compression experiments from the Sandia Z Machine and the OMEGA EP laser facility with density functional theory-based molecular dynamics calculations, we report the first pressure-density-temperature (P-ρ-T) relationship of natural iron (Fe)-bearing olivine ((Mg 0.91 Fe 0.09 ) 2 SiO 4 ) on the principal Hugoniot between 166 and 1,465 GPa. Additionally, we report the first reflectivities of natural olivine liquid in this pressure range. Compared to the magnesium-endmember forsterite (Mg 2 SiO 4 ), the presence of Fe in typical mantle abundance (~9 wt% FeO) alters the U S -u P relation of olivine. Furthermore, the shock temperature and reflectivity of olivine are indistinguishable from forsterite where experimental conditions overlap. Both forsterite and olivine increase in reflectivity (and hence optical conductivity) with increasing temperature, with a maximum reflectivity of ~31% at shock velocities greater than 22 km/s (~800 GPa).

58 GEOSCIENCES↗

High-resolution spectroscopy of barium monofluoride: Odd isotopologues, hyperfine structure, and isotope shifts

Barium monofluoride (BaF) is a promising molecular species for precision tests of fundamental symmetries and interactions. We present a combined theoretical and experimental study of BaF spectra and isotope shifts, focusing in particular on the poorly understood odd isotopologues 137 BaF and 135 BaF. By comparing state-of-the-art ab initio calculations with high-resolution fluorescence and absorption spectroscopy data, we provide a benchmark for electronic structure theory and disentangle the hyperfine and rovibrational spectra of the five most abundant isotopologues, from 138 BaF to 134 BaF. The comprehensive knowledge gained enables a King plot analysis of the isotope shifts that reveals the odd-even staggering of the barium nuclear charge radii. It also paths the way for improved laser cooling of rare BaF isotopologues and crucially supports future measurements of nuclear anapole and Schiff moments.

Kogel, Felix [Univ. of Stuttgart (Germany)] (ORCID↗

Atomic structure and grain shape evolution of nanodiamond during annealing in oxidizing atmosphere from neutron diffraction and MD simulations

In this work, real space pair-distribution function (PDF) analysis in conjunction with molecular dynamics (MD) simulations was applied to characterize crystal structure, shape and surface structure of diamond nanoparticles annealed in air up to 880 °C. MD simulated models of individual nanodiamond particles of various shapes and sizes were used to calculate theoretical radial distribution functions G(r) for powders with different grain size distributions. The calculated curves, showing distinct differences for the explored parameter space, were compared to the experimental G(r) functions obtained from neutron scattering for annealed detonated nanodiamond samples in a search for a model best describing the actual powders. It was found that as-synthesized nanodiamond grains are terminated by equally abundant (100), (110) and (111) crystallographic faces and their structure does not change upon annealing in air up to 280 °C. The shape and size of grains annealed above 480 °C changes due to surface etching. The etching rate differs between crystals faces and depends on temperature, resulting in the observed differences in grain shape. Below 680 °C the share of the (110) facets gradually increases and at 680 °C the grains become rhombic dodecahedra terminated solely by (110) surfaces. Further increase of the annealing temperature promotes formation of (111) surfaces and at 880 °C the preferred grain shape becomes an octahedron terminated by a variant of (111) surfaces possessing three dangling sp3 bonds per atom. At lower temperatures all grains become smaller and the size distribution moves towards the lower values. Furthermore, as the annealing temperature increases, the smallest grains disappear entirely and the average grain size grows back.

36 MATERIALS SCIENCE↗

Collisional Excitation of HCN by CO to Refine the Modeling of Cometary Comae

Here, we present the first dataset of collisional (de)-excitation rate coefficients of HCN induced by CO, one of the main perturbing gases in cometary atmospheres. The dataset spans the temperature range of 5–50 K. It includes both state-to-state rate coefficients involving the lowest ten and nine rotational levels of HCN and CO, respectively, and the so-called “thermalized” rate coefficients over the rotational population of CO at each kinetic temperature. The derivation of these coefficients exploited the good performance of the statistical adiabatic channel model (SACM) on top of an accurate interaction potential computed at the CCSD(T)-F12b/CBS level of theory. The reliability of the SACM approach was validated by comparison with full quantum calculations restricted at the lowest total angular momentum of the system. These results provide essential input to accurately model the distribution among the rotational energy levels and the abundance of HCN in cometary atmospheres, accounting for deviations from local thermodynamic equilibrium that typically occurs in such environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ultraviolet H 2 luminescence in molecular clouds induced by cosmic rays

Context. Galactic cosmic rays (CRs) play a crucial role in ionisation, dissociation, and excitation processes within dense cloud regions where UV radiation is absorbed by dust grains and gas species. CRs regulate the abundance of ions and radicals, leading to the formation of more and more complex molecular species, and determine the charge distribution on dust grains. A quantitative analysis of these effects is essential for understanding the dynamical and chemical evolution of star-forming regions. Aims. The CR-induced photon flux has a significant impact on the evolution of the dense molecular medium in its gas and dust components. This study evaluates the flux of UV photons generated by CRs to calculate the photon-induced dissociation and ionisation rates of a vast number of atomic and molecular species, as well as the integrated UV photon flux. Methods. To achieve these goals, we took advantage of recent developments in the determination of the spectra of secondary electrons, in the calculation of state-resolved excitation cross sections of H 2 by electron impact, and of photodissociation and photoionisation cross sections. Results. We calculated the H 2 level population of each rovibrational level of the X, B, C, B′, D, B″, D′, and a states. We then computed the UV photon spectrum of H 2 in its line and continuum components between 72 and 700 nm, with unprecedented accuracy, as a function of the CR spectrum incident on a molecular cloud, the H 2 column density, the isomeric H 2 composition, and the dust properties. The resulting photodissociation and photoionisation rates are, on average, lower than previous determinations by a factor of about 2, with deviations of up to a factor of 5 for the photodissociation of species such as AlH, C 2 H 2 , C 2 H 3 , C 3 H 3 , LiH, N 2 , NaCl, NaH, O 2 + , S2, SiH, l-C 4 , and l-C 5 H. A special focus is given to the photoionisation rates of H 2 , HF, and N 2 , as well as to the photodissociation of H 2 , which we find to be orders of magnitude higher than previous estimates. We give parameterisations for both the photorates and the integrated UV photon flux as a function of the CR ionisation rate, which implicitly depends on the H 2 column density, as well as the dust properties.

79 ASTRONOMY AND ASTROPHYSICS↗

Ultraviolet H 2 luminescence in molecular clouds induced by cosmic rays

Context. Galactic cosmic rays (CRs) play a crucial role in ionisation, dissociation, and excitation processes within dense cloud regions where UV radiation is absorbed by dust grains and gas species. CRs regulate the abundance of ions and radicals, leading to the formation of more and more complex molecular species, and determine the charge distribution on dust grains. A quantitative analysis of these effects is essential for understanding the dynamical and chemical evolution of star-forming regions. Aims. The CR-induced photon flux has a significant impact on the evolution of the dense molecular medium in its gas and dust components. This study evaluates the flux of UV photons generated by CRs to calculate the photon-induced dissociation and ionisation rates of a vast number of atomic and molecular species, as well as the integrated UV photon flux. Methods. To achieve these goals, we took advantage of recent developments in the determination of the spectra of secondary electrons, in the calculation of state-resolved excitation cross sections of H 2 by electron impact, and of photodissociation and photoionisation cross sections. Results. We calculated the H 2 level population of each rovibrational level of the X, B, C, B', D, B", D', and a states. We then computed the UV photon spectrum of H 2 in its line and continuum components between 72 and 700 nm, with unprecedented accuracy, as a function of the CR spectrum incident on a molecular cloud, the H 2 column density, the isomeric H 2 composition, and the dust properties. The resulting photodissociation and photoionisation rates are, on average, lower than previous determinations by a factor of about 2, with deviations of up to a factor of 5 for the photodissociation of species such as AlH, C 2 H 2 , C 2 H 3 , C 3 H 3 , LiH, N 2 , NaCl, NaH, O$^+_2$, S 2 , SiH, l-C 4 , and l-C 5 H. A special focus is given to the photoionisation rates of H 2 , HF, and N 2 , as well as to the photodissociation of H 2 , which we find to be orders of magnitude higher than previous estimates. We give parameterisations for both the photorates and the integrated UV photon flux as a function of the CR ionisation rate, which implicitly depends on the H 2 column density, as well as the dust properties.

79 ASTRONOMY AND ASTROPHYSICS↗

nmrrr : A Reproducible Workflow for Binning and Visualizing NMR Spectra From Environmental Samples

Nuclear magnetic resonance (NMR) spectroscopy is a useful tool for detection and identification of molecular structural information, with increasing applications in environmental sciences. NMR instrument outputs are however heterogeneous and require extensive post-processing, creating barriers to their use and application by non-specialists. Here, we report on a new open-source R package, nmrrr, that processes and visualizes spectral data obtained from one-dimensional solution-state and solid-state NMR experiments; the package also performs relevant calculations commonly applied in natural organic matter communities, such as computing the relative abundance of various functional groups. We document the package's installation, dependencies, and functions; and provide a standard workflow for processing NMR data. This package is currently available on CRAN and GitHub, and community contributions are welcome.

54 ENVIRONMENTAL SCIENCES↗

Photodriven electron-transfer dynamics in a series of heteroleptic Cu(I)–anthraquinone dyads

Solar fuels catalysis is a promising route to efficiently harvesting, storing, and utilizing abundant solar energy. To achieve this promise, however, molecular systems must be designed with sustainable components that can balance numerous photophysical and chemical processes. To that end, here we report on the structural and photophysical characterization of a series of Cu(I)–anthraquinone-based electron donor–acceptor dyads. The dyads utilized a heteroleptic Cu(I) bis-diimine architecture with a copper(I) bis-phenanthroline chromophore donor and anthraquinone electron acceptor. We characterized the structures of the complexes using x-ray crystallography and density functional theory calculations and the photophysical properties via resonance Raman and optical transient absorption spectroscopy. The calculations and resonance Raman spectroscopy revealed that excitation of the Cu(I) metal-to-ligand charge-transfer (MLCT) transition transfers the electron to a delocalized ligand orbital. The optical transient absorption spectroscopy demonstrated that each dyad formed the oxidized copper–reduced anthraquinone charge-separated state. Unlike most Cu(I) bis-phenanthroline complexes where increasingly bulky substituents on the phenanthroline ligands lead to longer MLCT excited-state lifetimes, here, we observe a decrease in the long-lived charge-separated state lifetime with increasing steric bulk. The charge-separated state lifetimes were best explained in the context of electron-transfer theory rather than with the energy gap law, which is typical for MLCT excited states, despite the complete conjugation between the phenanthroline and anthraquinone moieties.

14 SOLAR ENERGY↗

Exploring the mass surface near the rare-earth abundance peak via precision mass measurements at JYFLTRAP

The JYFLTRAP double Penning trap at the Ion Guide Isotope Separator On-Line facility has been used to measure the atomic masses of 13 neutron-rich rare-earth isotopes. Eight of the nuclides, 161 Pm, 163 Sm, 164,165 Eu, 167 Gd, and 165,167,168 Tb, were measured for the first time. The systematics of the mass surface has been studied via one- and two-neutron separation energies as well as neutron pairing-gap and shell-gap energies. The proton-neutron pairing strength has also been investigated. Furthermore, the impact of the new mass values on the astrophysical rapid neutron capture process has been studied. The calculated abundance distribution results in a better agreement with the solar abundance pattern near the top of the rare-earth abundance peak at around A≈165.

74 ATOMIC AND MOLECULAR PHYSICS↗

Data-Centric Development of Lignin Structure–Solubility Relationships in Deep Eutectic Solvents Using Molecular Simulations

Lignin is a natural source of aromatic chemicals with significant potential as an abundant, renewable feedstock for value-added products. Deep eutectic solvents (DES)–solvents composed of a hydrogen bond donor (HBD) and acceptor (HBA) in varying ratios–have emerged as a highly tunable class of solvents for lignin solubilization. However, the variety of possible DES compositions and limited molecular-scale understanding of lignin solubility makes solvent selection a challenge without laborious trial-and-error experimentation. To address these challenges, we use classical molecular dynamics (MD) simulations to study the interactions of lignin model compounds with various DES–water systems. Quantitative parameters (descriptors) were calculated by postprocessing the MD results and used to train a regression model that predicts experimentally determined solubilities of lignin model compounds. This approach revealed that the most important descriptors of solubility are the system temperature, solute hydrophilicity, and metrics quantifying hydrogen bonding. Maximizing the interactions between solute–HBD (hydrophobic group), water–HBD (hydrophilic group), and water–HBA molecules led to the highest model compound solubility. Our results support a hydrotropic mechanism in which extensive DES–water hydrogen bonding and favorable HBD interactions with the solute promote high solubility. We applied the regression model derived using model compounds to predict the solubility of representative lignin oligomers. The model predicted lignin oligomers’ solubilities in good agreement with experiments, indicating that the simulations of model compounds can be extended to predict the solubility of larger lignin compounds across a range of solvent compositions and temperatures. Furthermore, these findings provide new molecular-scale insight into lignin solubilization mechanisms and a new method for computationally screening potential solvent systems for lignin valorization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Inelastic rate coefficients for collisions of C 4 H - with H 2

Carbon-chain anions were recently detected in the interstellar medium. These very reactive species are used as tracers of the physical and chemical conditions in a variety of astrophysical environments. However, the local thermodynamic equilibrium conditions are generally not fulfilled in these environments. Therefore, collisional as well as radiative rates are needed to accurately model the observed emission lines. We determine in this work the state-to-state rate coefficients of C 4 H- in collision with both ortho- and para-H 2 . A new ab initio 4D potential energy surface was computed using explicitly correlated coupled-cluster procedures. This surface was then employed to determine rotational excitation and de-excitation cross-sections and rate coefficients for the first 21 rotational levels (up to rotational level j 1 = 20) using the close-coupling method, while the coupled-state approximation was used to extend the calculations up to j 1 = 30. State-to-state rate coefficients were obtained for the temperature range 2–100K. The differences between the ortho- and para-H 2 rate coefficients are found to be small.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigating the Effect of Water on the Mechanical Properties of Cellulose from Multiscale Molecular Dynamics Simulations

Classical molecular dynamics (MD) simulations provide insight into the structure and physicochemical properties of materials with atomic resolution. However, the length and time scales accessible to atomistic MD are orders of magnitude smaller than many relevant processes such as the response of a bulk material to experimentally accessible strain rates, which presents challenges when comparing models to experimental measurements. Bottom-up coarse-graining provides a means for systematically mapping atomistic information to lower resolution models to increase the length and time scales achievable by simulation. Cellulose is an abundant carbohydrate biopolymer with applications to many fields of research, such as materials science and renewable energy, due to its desirable mechanical properties and viability for conversion into biofuel. The effect of moisture content on the Young's modulus of cellulose is of special interest due to its native environment often being in the hydrated secondary plant cell wall and the grinding energy requirements for biomass feedstock preprocessing. The current work investigates the effects of water solvent on the Young's modulus of cellulose calculated from coarse-grained MD mechanical stress simulations. The coarse-grained model was parametrized from atomistic MD calculations of cellulose-cellulose potentials of mean force using umbrella sampling techniques under vacuum and solvated conditions. The Young's moduli of the coarse-grained cellulose assemblies parametrized from cellulose in vacuum or solvated in water were computed via mechanical stress simulations to highlight the importance of capturing solvent interactions for modeling the mechanical behavior of cellulose.

BASIC BIOLOGICAL SCIENCES,RADIATION PROTECTION AND↗

The RADIOSTAR Project

Radioactive nuclei are the key to understanding the circumstances of the birth of our Sun because meteoritic analysis has proven that many of them were present at that time. Their origin, however, has been so far elusive. The ERC-CoG-2016 RADIOSTAR project is dedicated to investigating the production of radioactive nuclei by nuclear reactions inside stars, their evolution in the Milky Way Galaxy, and their presence in molecular clouds. So far, we have discovered that: (i) radioactive nuclei produced by s l o w ( 107 Pd and 182 Hf) and r a p i d ( 129 I and 247 Cm) neutron captures originated from stellar sources —asymptotic giant branch (AGB) stars and compact binary mergers, respectively—within the galactic environment that predated the formation of the molecular cloud where the Sun was born; (ii) the time that elapsed from the birth of the cloud to the birth of the Sun was of the order of 10 7 years, and (iii) the abundances of the very short-lived nuclei 26 Al, 36 Cl, and 41 Ca can be explained by massive star winds in single or binary systems, if these winds directly polluted the early Solar System. Our current and future work, as required to finalise the picture of the origin of radioactive nuclei in the Solar System, involves studying the possible origin of radioactive nuclei in the early Solar System from core-collapse supernovae, investigating the production of 107 Pd in massive star winds, modelling the transport and mixing of radioactive nuclei in the galactic and molecular cloud medium, and calculating the galactic chemical evolution of 53 Mn and 60 Fe and of the p -process isotopes 92 Nb and 146 Sm.

59 BASIC BIOLOGICAL SCIENCES↗

Thermal conductivity of CaSiO 3 perovskite at lower mantle conditions

Thermal conductivity (κ) of mantle minerals is key to understanding dynamics in the deep Earth. It controls the style of mantle convection and the timescale of cooling both the mantle and the core. Cubic CaSiO 3 perovskite (CaPv) is the third most abundant mineral in the lower mantle (7 vol %). Despite its importance, no theoretical or experimental estimate of CaPv's κ is available. Theoretical investigations of its properties are challenging because of its strong anharmonicity. Experimental measurements at relevant pressures and temperatures are equally challenging. Here we present ab initio results for CaPv's κ obtained using the phonon quasiparticle approach to address its strong anharmonicity. We also offer experimental measurements of κ up to 67 GPa and 1950 K. Predictions and measurements are in good agreement and reveal a surprisingly large κ for cubic CaPv that can be explained on the basis of its high crystal structure symmetry. Despite its relatively low abundance, CaPv's κ increases the lower mantle κ by ~10 %, if accounted for. κ of mantle regions enriched in subducted crustal materials will be more strongly impacted.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗