Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular properties”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A step in the direction of resolving the paradox of Perdew–Zunger self-interaction correction. II. Gauge consistency of the energy density at three levels of approximation

The Perdew-Zunger(PZ) self-interaction correction (SIC) was designed to correct the one-electron limit of any approximate density functional for the exchange-correlation (xc) energy, while yielding no correction to the exact functional. Unfortunately, it spoils the slowly-varying-in-space limits of the uncorrected approximate functionals, where those functionals are right by construction. The right limits can be restored by locally scaling down the energy density of the PZ SIC in many-electron regions, but then a spurious correction to the exact functional would be found unless the self-Hartree and exact self-xc terms of the PZ SIC energy density were expressed in the same gauge. Only the local density approximation satisfies the same-gauge condition for the energy density, which explains why the recent local-scaling SIC (LSIC) is found here to work excellently for atoms and molecules only with this basic approximation, and not with the more advanced generalized gradient approximations (GGAs) and meta-GGAs, which lose the Hartree gauge via simplifying integrations by parts. The transformation of energy density that achieves the Hartree gauge for the exact xc functional can also be applied to approximate functionals. Doing so leads to a simple scaled-down self-interaction (sdSIC) correction that is typically much more accurate than PZ SIC in tests for many molecular properties (including equilibrium bond lengths). The present work shows unambiguously that the largest errors of PZ SIC applied to standard functionals at three levels of approximation can be removed by restoring their correct slowlyvarying- density limits. Here, it also confirms the relevance of these limits to atoms and molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

An extensive ab initio study of the structures, vibrational spectra, quadratic force fields, and relative energetics of three isomers of Cl2O2

Quantum mechanical computational methods are employed for an ab initio investigation of: (1) the molecular properties of the lowest isomers of the ClO dimer; and (2) predicted molecular and thermochemical properties. Techniques employed include electron correlation and particularly singles and doubles coupled-cluster (CCSD) theory with or without perturbational estimates of the effects of connected triple excitations. The isomers ClOClO and ClClO2 are found to have higher energies than the ClOOCl isomer, and the theoretical vibrational frequencies of the isomers are well correlated with experimental data. Experimental values of the heat of formation for the isomers are also compared with calculations based on an isodesmic reaction with Cl2O, H2O, and HOOH.

Lee, Timothy J.↗
DOE OSTI2021

Equation-of-motion coupled-cluster method with double electron-attaching operators: Theory, implementation, and benchmarks

We report a production-level implementation of the equation-of-motion (EOM) coupled-cluster (CC) method with double electron-attaching (DEA) EOM operators of 2 p and 3 p 1 h types, EOM-DEA-CCSD. This ansatz, suitable for treating electronic structure patterns that can be described as two-electrons-in-many orbitals, represents a useful addition to the EOM-CC family of methods. We analyze the performance of EOM-DEA-CCSD for energy differences and molecular properties. By considering reduced quantities, such as state and transition one-particle density matrices, we compare EOM-DEA-CCSD wave functions with wave functions computed by other EOM-CCSD methods. The benchmarks illustrate that EOM-DEA-CCSD is capable of treating diradicals, bond-breaking, and some types of conical intersections.

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

“Freedom of design” in chemical compound space: towards rational in silico design of molecules with targeted quantum-mechanical properties

The rational design of molecules with targeted quantum-mechanical (QM) properties requires an advanced understanding of the structure–property/property–property relationships (SPR/PPR) that exist across chemical compound space (CCS). In this work, we analyze these fundamental relationships in the sector of CCS spanned by small (primarily organic) molecules using the recently developed QM7-X dataset, a systematic, extensive, and tightly converged collection of 42 QM properties corresponding to ≈4.2M equilibrium and non-equilibrium molecular structures containing up to seven heavy/non-hydrogen atoms (including C, N, O, S, and Cl). By characterizing and enumerating progressively more complex manifolds of molecular property space—the corresponding high-dimensional space defined by the properties of each molecule in this sector of CCS—our analysis reveals that one has a substantial degree of flexibility or “freedom of design” when searching for a single molecule with a desired pair of properties or a set of distinct molecules sharing an array of properties. To explore how this intrinsic flexibility manifests in the molecular design process, we used multi-objective optimization to search for molecules with simultaneously large polarizabilities and HOMO–LUMO gaps; analysis of the resulting Pareto fronts identified non-trivial paths through CCS consisting of sequential structural and/or compositional changes that yield molecules with optimal combinations of these properties.

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

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy–cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe–Salpeter methods, second-order Møller–Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE’s functionality, including excited-state methods, RPA and Green’s function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE’s current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE’s development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

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

Multiple-objective Reinforcement Learning for Inverse Design and Identification

The aim of inverse chemical design is to develop new molecules with given optimized molecular properties or objectives. Recently, generative deep learning (DL) networks are considered as the state-of-the-art in inverse chemical design and have achieved early success in generating molecular structures with desired properties in the pharmaceutical and material chemistry fields. However, satisfying a large number (> 10 objectives) of molecular objectives is a limitation of current generative models. To improve the model’s ability to handle a large number of molecule design objectives, we developed a Reinforcement Learning (RL) based generative framework to optimize chemical molecule generation. Our use of Curriculum Learning (CL) to fine-tune the pre-trained generative network allowed the model to satisfy up to 21 objectives and increase the generative network’s robustness. The experiments show that the proposed multiple-objective RL-based generative model can correctly identify unknown molecules with an 83% to 100% success rate, compared to the baseline approach of 0%. Additionally, this proposed generative model is not limited to just chemistry research challenges; we anticipate that problems that utilize RL with multiple objectives will benefit from this framework.

reinforcement learning, generative model, inverse ↗
DOE OSTI2021

Three-Dimensional Modeling of the Explicit Chemistry of Organic Aerosols: Formation, Removal, and Comparison with Observations

The major goals of the project were to further our understanding of the formation of secondary organic aerosols (SOA) from major anthropogenic and biogenic precursors, and to evaluate the GECKO-A chemistry model in the context of three DOE-led field campaigns (GoAmazon, BBOP, and HI-SCALE) with the long-term objective of improving the treatment of organic gases and particles in three-dimensional regional and global models. The three field campaigns were selected to evaluate and improve specific aspects of GECKO-A: (i) GoAmazon to help quantify interactions of biogenic and anthropogenic emissions; (ii) BBOP to predict the UV-visible absorption properties of biomass burning aerosols (brown carbon) from explicit molecular properties; and (iii) HI-SCALE to evaluate new model components that include aqueous (cloud) processing of organic gases and aerosols.

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

Oxadiazole-based Heterocycles as Building Block for Material Property Control: Design, Synthesis, and Characterization of 3,4-Bis(3-(4-nitro-1,2,5-oxadiazol-3-yl)- 1,2,4-oxadiazol-5-yl)-1,2,5-oxadiazole (LLM-210)

This paper presents an approach to novel oxadiazole-based melt-castable energetic material, LLM-210. Used in pour-in process and 3D printing, melt-castable energetic materials (MCEM) must possess required physical and safety properties of industry-standard energetic materials, also a narrow range of melting point between 80-100 °C. Melting point is difficult to accurately predict, therefore, searching for MCEM is challenging. 3,4-Bis(5-(4-nitro-1,2,5-oxadiazol-3-yl)-1,2,4-oxadiazol3-yl)-1,2,5-oxadiazole (LLM-205), composed of five-oxadiazole rings allied via carbon-carbon bonds, was selected as new MCEM candidate. The target molecule was synthesized and characterized, showing a density of 1.807g/cm3 ; thermal decomposing temperature at 299 °C; and insensitive to external stimuli, but the melting point of 104 °C was out of the rang. Based on the molecular properties and ab initio calculations, 3,4-bis(3-(4-nitro-1,2,5-oxadiazol-3-yl)-1,2,4-oxadiazol-5-yl)- 1,2,5-oxadiazole (LLM-210), an isomer of LLM-205, was designed and synthesized from a different way. LLM-210 was characterized as MCEM, possessing melting point of 87 °C; density of 1.812 g/cm3 ; and exothermally decomposing temperature at 314 °C; being insensitive to external stimuli. In this study, the relationship of molecular structure to properties of LLM-205 and LLM-210 based on the results of Density Functional Theory (DFT) and X-ray crystallographic analysis is also briefly discussed.

Chemistry - Chemical explosives↗
NASA NTRS2023

Molecular Dynamics Simulations of Austenite-Martensite Interfaces in NiTi Shape Memory Alloys

The unique properties of shape memory alloys (SMAs) arise from a reversible martensitic transformation. The nucleation and migration of austenite-martensite interfaces are the key to understanding the SMA properties. Molecular dynamics (MD) simulations can provide important atomic-scale information about these aspects of the transformation, but their time scales prevent the interface formation under near-equilibrium conditions relevant to experiment. We present a new MD methodology which allows for the natural formation of energetically preferred austenite-martensite interfaces under near-equilibrium conditions. Our simulation demonstrates that the interfaces in NiTi are semi-coherent, composed of a series of terrace planes and structural disconnections, and they migrate rapidly through single crystals with only a small thermodynamic driving force. In bi-crystals and polycrystals, the migration of these same interfaces is significantly impeded by grain boundaries and stored elastic energy. This behavior can result in SMA hysteresis via several mechanisms associated with nucleation and non-elastic strain accommodation.

Gabriel Plummer↗
NASA NTRS2024

Effect of the Liquid-Vacuum Transition on the Relative Abundances of Amino and Fatty Acids Sought as Biosignatures on Icy Ocean Worlds

Interpreting measurements by robotic missions searching for signs of life in material erupted by cryovolcanic ocean worlds (e.g., Enceladus, Europa, Ceres) requires linking the composition of this material to that of its subsurface liquid source. Irrespective of the properties of –and processes along– the path of ejection, material erupted on airless worlds undergoes a transition from liquid to vacuum. Here, we investigate experimentally the effect of a direct transition from liquid to vacuum environments on the relative abundances of amino and carboxylic (fatty) acids, a metric used to distinguish between biological and abiotic sources for these compounds. Amino acids were dissolved in parent solutions prepared with 1 wt.% NaCl and pH 9–10 to match properties inferred from Enceladus plume grains. Compositional analysis of the dry solids resulting from the injection of this solution into vacuum indicates a < 50% change in amino acid abundances relative to glycine. Injection of two fatty acids, phenylacetic acid (soluble and undersaturated in the parent solution) and palmitic acid (insoluble and supersaturated) resulted in a change of less than a factor of 7 in their relative abundances. At this bulk scale (all grains aggregated together), proportions of amino and fatty acids were sufficiently conserved to allow distinction between end-member example biological and abiotic sources. We did not find strong correlations between the relative enrichment or depletion of amino or fatty acids with molecular properties such as molecular mass, hydrophobicity, functionality, or charge, other than high under-or supersaturation. Most of the organic residue was deposited along the experimental injection path, which suggests that cryovolcanic conduit walls on icy ocean worlds may be enriched in organic material relative to material having undergone the liquid-vacuum transition.

Amino acids↗
DOE OSTI2024

Unified understanding of the impact of semiflexibility, concentration, and molecular weight on macromolecular-scale ring diffusion

Conformationally fluctuating, globally compact macromolecules such as polymeric rings, single-chain nanoparticles, microgels, and many-arm stars display complex dynamic behaviors due to their rich topological structure and intermolecular organization. Synthetic rings are hybrid objects with conformations that display both ideal random walk and compact globular features, which can serve as models of genomic DNA. To date, emphasis has been placed on the effect of ring molecular weight on their unusual behaviors. Here, we combine simulations and a microscopic force-level theory to build a unified understanding for how key aspects of ring dynamics depend on different tunable molecular properties including backbone rigidity, monomer concentration, degree of traditional entanglement, and molecular weight. Our large-scale molecular dynamics simulations of ring melts with very different backbone stiffnesses reveal unanticipated behaviors which agree well with our generalized theory. This includes a universal master curve for center-of-mass diffusion constants as a function of molecular weight scaled by a chemistry and thermodynamic state-dependent critical molecular weight that generalizes the concept of an entanglement cross-over for linear chains. The key physics is how backbone rigidity and monomer concentration induced changes of the entanglement length, interring packing, degree of interpenetration, and liquid compressibility slow down space-time dynamic-force correlations on macromolecular scales. A power law decay of the center-of-mass diffusion constant with inverse molecular weight squared is the first consequence, followed by an ultraslow activated hopping transport regime. Our results set the stage to address slow dynamics and kinetic arrest in different families of compact synthetic and biological polymeric systems.

Science & Technology - Other Topics↗
DOE OSTI2023

Gaussian approximation potentials for accurate thermal properties of two-dimensional materials

Two-dimensional materials (2DMs) continue to attract a lot of attention, particularly for their extreme flexibility and superior thermal properties. Molecular dynamics simulations are among the most powerful methods for computing these properties, but their reliability depends on the accuracy of interatomic interactions. While first principles approaches provide the most accurate description of interatomic forces, they are computationally expensive. In contrast, classical force fields are computationally efficient, but have limited accuracy in interatomic force description. Machine learning interatomic potentials, such as Gaussian Approximation Potentials, trained on density functional theory (DFT) calculations offer a compromise by providing both accurate estimation and computational efficiency. Here, in this work, we present a systematic procedure to develop Gaussian approximation potentials for selected 2DMs, graphene, buckled silicene, and h-XN (X = B, Al, and Ga, as binary compounds) structures. We validate our approach through calculations that require various levels of accuracy in interatomic interactions. The calculated phonon dispersion curves and lattice thermal conductivity, obtained through harmonic and anharmonic force constants (including fourth order) are in excellent agreement with DFT results. HIPHIVE calculations, in which the generated GAP potentials were used to compute higher-order force constants instead of DFT, demonstrated the first-principles level accuracy of the potentials for interatomic force description. Molecular dynamics simulations based on phonon density of states calculations, which agree closely with DFT-based calculations, also show the success of the generated potentials in high-temperature simulations.

2D Materials↗
DOE OSTI2022

Polar liquids at charged interfaces: A dipolar shell theory

The structure of polar liquids and electrolytic solutions, such as water and aqueous electrolytes, at interfaces underlies numerous phenomena in physics, chemistry, biology, and engineering. In this work, we develop a continuum theory that captures the essential features of dielectric screening by polar liquids at charged interfaces, including decaying spatial oscillations in charge and mass, starting from the molecular properties of the solvent. The theory predicts an anisotropic dielectric tensor of interfacial polar liquids previously studied in molecular dynamics simulations. We explore the effect of the interfacial polar liquid properties on the capacitance of the electrode/electrolyte interface and on hydration forces between two plane-parallel polarized surfaces. In the linear response approximation, we obtain simple formulas for the characteristic decay lengths of molecular and ionic profiles at the interface.

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

Prediction of Research/Motor Octane Number and Octane Sensitivity Using Artificial Neural Networks

Octane sensitivity (OS), dfined as the research octane number (RON) minus the motor octane number (MON) of a given fuel, has gained interest among researchers due to its apparent ffect on knocking conditions in internal combustion engines. Compounds with a high OS enable higher efficiencies, especially with respect to advanced compression ignition engines. RON/MON must be experimentally tested to determine OS; however, the experimental methods utilized require a substantial amount of time, a significant monetary investment, and specialized equipment. To this end, predictive computational models trained with existing experimental data and molecular properties would allow for the preemptive screening of compounds prior to performing these experiments. The present work proposes two methods for predicting the OS of a given compound: using artficial neural networks (ANNs) trained with quantitative structure-property relationship (QSPR) descriptors to predict RON and MON individually to compute OS from RON/MON predictions (derived octane sensitivity, dOS), and using an ANN trained with QSPR descriptors to directly predict OS. ANNs trained to predict RON and MON achieved test set root-mean-square errors (RMSEs) of 10.499 and 7.551 respectively. dOS calculations were found to have a test set RMSE of 6.432 while predicting OS directly resulted in a test set RMSE of 7.019, showing it is more bene cial to obtain OS from RON/MON predictions than predicting it directly. Furthermore, relationships between individual QSPR descriptors and RON/MON/OS are discussed, highlighting correlations between specfic molecular features and these properties.

09 BIOMASS FUELS↗
DOE OSTI2025

Probing specific ion effects at air-aqueous dibutyl phosphate interfaces using vibrational sum frequency generation spectroscopy

Molecular properties at air–liquid and liquid–liquid interface hold the key to many processes involving molecular transport across phase boundaries from aerosol formation to carbon cycling and material separation using solvent extraction techniques. Using dibutyl phosphate (DBP) as a representative for partially aqueous soluble surfactants, the specific ion effect (SIE) of the Hofmeister series cations Cs + , Na + , Li + , and Mg 2+ on the partition and interaction between surfactant molecules and water molecules in the air–aqueous interface are investigated using vibrational sum frequency generation spectroscopy and surface tension measurements. In the presence of 1 mM and 1M bulk aqueous phase ionic strength salt concentrations, fundamental qualitative relationships are observed for the salting out of DBP relative to bulk aqueous phase nitrate salt concentrations and the specific cations species. At 1 mM ionic strength, the interfacial charge and hence the interfacial potential modulates the electrostatic interactions; in particular, the counter cations partially screen the negatively charged interface induced by the DBP in a direct Hofmeister order. At 1M ionic strength, the electric field at the interface or interfacial potential is effectively neutralized, and the counter cations promote the partitioning of DBP to the interface depending on their specific interaction with the DBP head group and metal ion hydration properties. The present results lay a foundation to study SIEs of heavier metals on hydrophobic-aqueous DBP interfaces.

Louie, Christina [Washington State University, Pul↗