Search NASA⌕ Search

SEARCH · Search NASA

Results for “MOLECULAR BONDING”

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

Actinide 2-metallabiphenylenes that satisfy Hückel’s rule

Aromaticity and antiaromaticity, as defined by Hückel’s rule, are key ideas in organic chemistry, and are both exemplified in biphenylene—a molecule that consists of two benzene rings joined by a four-membered ring at its core. Biphenylene analogues in which one of the benzene rings has been replaced by a different (4n + 2) π-electron system have so far been associated only with organic compounds. In addition, efforts to prepare a zirconabiphenylene compound resulted in the isolation of a bis(alkyne) zirconocene complex instead. Here we report the synthesis and characterization of, to our knowledge, the first 2-metallabiphenylene compounds. Single-crystal X-ray diffraction studies reveal that these complexes have nearly planar, 11-membered metallatricycles with metrical parameters that compare well with those reported for biphenylene. Nuclear magnetic resonance spectroscopy, in addition to nucleus-independent chemical shift calculations, provides evidence that these complexes contain an antiaromatic cyclobutadiene ring and an aromatic benzene ring. Finally, spectroscopic evidence, Kohn–Sham molecular orbital compositions and natural bond orbital calculations suggest covalency and delocalization of the uranium f 2 electrons with the carbon-containing ligand.

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

Ultrafast structural dynamics of UV photoexcited cis , cis -1,3-cyclooctadiene observed with time-resolved electron diffraction

Conjugated diene molecules are highly reactive upon photoexcitation and can relax through multiple reaction channels that depend on the position of the double bonds and the degree of molecular rigidity. Understanding the photoinduced dynamics of these molecules is crucial for establishing general rules governing the relaxation and product formation. Here, in this study, we investigate the femtosecond time-resolved photoinduced excited-state structural dynamics of cis,cis-1,3-cyclooctadiene, a large-flexible cyclic conjugated diene molecule, upon excitation with 200 nm using mega-electron-volt ultrafast electron diffraction and trajectory surface hopping dynamics simulations. We tracked the photoinduced structural changes from the Franck–Condon region through the conical intersection seam to the ground state. Our findings revealed a novel primary reaction coordinate involving ring distortion, where the ring stretches along one axis and compresses along the perpendicular axis. The nuclear wavepacket remains compact along this reaction coordinate until it reaches the conical intersection seam, and it rapidly spreads as it approaches the ground state, where multiple products are formed.

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

Temperature dependence of nuclear quantum effects on liquid water via artificial neural network model based on SCAN meta-GGA functional

We investigate the temperature dependence of nuclear quantum effects (NQEs) on structural and dynamic properties of liquid water by training a neural network force field using first-principles molecular dynamics (FPMD) based on the strongly constrained and appropriately normed meta-generalized gradient approximation exchange-correlation approximation. The FPMD simulation based on density functional theory has become a powerful computational approach for studying a wide range of condensed phase systems. However, its large computational cost makes it difficult to incorporate NQEs in the simulation and investigate temperature dependence of various properties. To circumvent this difficulty, we use an artificial neural network model and employ the thermostatted ring polymer MD approach for studying the temperature dependence of NQEs on various properties. The NQEs generally bring the radial distribution functions closer to the experimental measurements. Translational diffusivity and rotational dynamics of water molecules are both slowed down by the NQEs. The competing inter-molecular and intra-molecular quantum effects on hydrogen bonds, as discussed by Habershon, Markland, and Manolopoulos [J. Chem. Phys. 131(2), 024501 (2019)], can explain the observed temperature dependence of the NQEs on the dynamical properties in our simulation.

Yao, Yi (ORCID:0000000299576979)↗
DOE OSTI2025

Influence of counterion substitution on the properties of imidazolium-based ionic liquid clusters

Due to their unique physiochemical properties that may be tailored for specific purposes, ionic liquids (ILs) have been investigated for various applications, including chemical separations, catalysis, energy storage, and space propulsion. The different cations and anions comprising ILs may be selected to optimize a range of desired properties, such as thermal stability, ionic conductivity, and volatility, leading to the designation of certain ILs as designer “green” solvents. The effect of counterions on the properties of ILs is of both fundamental scientific interest and technological importance. Herein, we report a systematic experimental and theoretical investigation of the size, charge, stability toward dissociation, and geometric/electronic structure of 1-ethyl-3-methyl imidazolium (EMIM)-based IL clusters containing two different atomic counterions (i.e., bromide [Br − ] and iodide [I − ]). This work extends our studies of EMIM + cations with atomic chloride (Cl − ) and molecular tetrafluoroborate (BF 4 − ) anions reported previously by Baxter et al. [Chem. Mater. 34, 2612 (2022)] and Zhang et al . [J. Phys. Chem. Lett. 11, 6844 (2020)], respectively. Distributions of anionic IL clusters were generated in the gas phase using electrospray ionization and characterized by high mass resolution mass spectrometry, energy-resolved collision-induced dissociation, and negative ion photoelectron spectroscopy experiments. The experimental results reveal anion-dependent trends in the size distribution, relative abundance, ionic charge state, stability toward dissociation, and electron binding energies of the IL clusters. Complementary global optimization theory provides molecular-level insights into the bonding and electronic structure of a selected subset of clusters, including their low energy structures and electrostatic potential maps, and how these fundamental characteristics are influenced by anion substitution. Collectively, our findings demonstrate how the fundamental properties of ILs, which determine their suitability for many applications, may be tuned by substituting counterions. These observations are critical in the sub-nanometer cluster size regime where phenomena do not scale predictably to the bulk phase, and each atom counts toward determining behavior.

cluster↗
DOE OSTI2020

Materials Data on Cl2 by Materials Project

Cl2 is Cubic alpha N2-like structured and crystallizes in the tetragonal P4_2/ncm space group. The structure is zero-dimensional and consists of eight molecular chlorine molecules. Cl is bonded in a single-bond geometry to one Cl atom. The Cl–Cl bond length is 1.99 Å.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Materials Data on Cl2 by Materials Project

Cl2 is Cubic alpha N2-like structured and crystallizes in the orthorhombic Cmce space group. The structure is zero-dimensional and consists of four molecular chlorine molecules. Cl is bonded in a single-bond geometry to one Cl atom. The Cl–Cl bond length is 2.01 Å.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Ordering of Trotterization: Impact on Errors in Quantum Simulation of Electronic Structure

Trotter–Suzuki decompositions are frequently used in the quantum simulation of quantum chemistry. They transform the evolution operator into a form implementable on a quantum device, while incurring an error—the Trotter error. The Trotter error can be made arbitrarily small by increasing the Trotter number. However, this increases the length of the quantum circuits required, which may be impractical. It is therefore desirable to find methods of reducing the Trotter error through alternate means. The Trotter error is dependent on the order in which individual term unitaries are applied. Due to the factorial growth in the number of possible orderings with respect to the number of terms, finding an optimal strategy for ordering Trotter sequences is difficult. In this paper, we propose three ordering strategies, and assess their impact on the Trotter error incurred. Initially, we exhaustively examine the possible orderings for molecular hydrogen in a STO-3G basis. We demonstrate how the optimal ordering scheme depends on the compatibility graph of the Hamiltonian, and show how it varies with increasing bond length. We then use 44 molecular Hamiltonians to evaluate two strategies based on coloring their incompatibility graphs, while considering the properties of the obtained colorings. We find that the Trotter error for most systems involving heavy atoms, using a reference magnitude ordering, is less than 1 kcal/mol. Relative to this, the difference between ordering schemes can be substantial, being approximately on the order of millihartrees. The coloring-based ordering schemes are reasonably promising—particularly for systems involving heavy atoms—however further work is required to increase dependence on the magnitude of terms. Finally, we consider ordering strategies based on the norm of the Trotter error operator, including an iterative method for generating the new error operator terms added upon insertion of a term into an ordered Hamiltonian.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Direct functionalization of C-H bonds by electrophilic anions

Bond formation between two of the most inert molecules, alkanes and [B12X12]2- (X = Cl, Br), is performed in a two step process. Fragmentation of [B12X12]2- in the gas phase generates highly reactive [B12X11]- ions which spontaneously react with alkanes. The reaction mechanism was investigated using tandem mass spectrometry and gas-phase vibrational spectroscopy combined with electronic structure calculations. This effort revealed the unprecedented and counterintuitive electrophilic substitution of a proton in an alkane by a negatively charged ion resulting in a B-C bond formation. The product is a dianionic [B12X11CnH2n+1]2- species, in which H+ is electrostatically bound to the dianion. High flux ion soft-landing was performed to co-deposit [B12X11]- and complex organic molecules (phthalates) in thin layers on surfaces. Molecular structure analysis of the product films revealed that C-H functionalization by [B12X11]- occured in the presence of other more reactive functional groups. This observation demonstrates the utility of highly reactive fragment ions for defined bond formation processes and may pave the way to the use of gas phase ion chemistry for chemical synthesis in the condensed phase.

Warneke, Jonas↗
DOE OSTI2024

Deep-Learning Interatomic Potential Connects Molecular Structural Ordering to the Macroscale Properties of Polyacrylonitrile

Polyacrylonitrile (PAN) is an important commercial polymer, bearing atactic stereochemistry resulting from nonselective radical polymerization. As such, an accurate, fundamental understanding of governing interactions among PAN molecular units is indispensable for advancing the design principles of final products at reduced processability costs. While ab initio molecular dynamics (AIMD) simulations can provide the necessary accuracy for treating key interactions in polar polymers, such as dipole–dipole interactions and hydrogen bonding, and analyzing their influence on the molecular orientation, their implementation is limited to small molecules only. Herein, we show that the neural network interatomic potentials (NNIPs) that are trained on the small-scale AIMD data (acquired for oligomers) can be efficiently employed to examine the structures and properties at large scales (polymers). NNIP provides critical insight into intra- and interchain hydrogen-bonding and dipolar correlations and accurately predicts the amorphous bulk PAN structure validated by modeling the experimental X-ray structure factor. Furthermore, the NNIP-predicted PAN properties, such as density and elastic modulus, are in good agreement with their experimental values. Overall, the trend in the elastic modulus is found to correlate strongly with the PAN structural orientations encoded in the Hermans orientation factor. In conclusion, this study enables the ability to predict the structure–property relations for PAN and analogues with sustainable ab initio accuracy across scales.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Au 3 -to-Ag 3 coordinate-covalent bonding and other supramolecular interactions with covalent bonding strength

An efficient strategy for designing charge-transfer complexes using coinage metal cyclic trinuclear complexes (CTCs) is described herein. Due to opposite quadrupolar electrostatic contributions from metal ions and ligand substituents, [Au(m-Pz-(i-C 3 H7) 2 )] 3 ∙[Ag(m-Tz-(n-C 3 F 7 ) 2 )] 3 (Pz = pyrazolate, Tz = triazolate) has been obtained and its structure verified by single crystal X-ray diffraction – representing the 1st crystallographically-verified M 3 @M' 3 stacked adduct of monovalent coinage metal CTCs. Abundant supramolecular interactions with aggregate covalent bonding strength arise from a combination of M–M' (Au / Ag), metal–π, π–π interactions and hydrogen bonding in this charge-transfer complex, according to density functional theory analyses, yielding a computed binding energy of 66 kcal mol -1 between the two trimer moieties – a large value for intermolecular interactions between adjacent d 10 centres (nearly doubling the value for a recently claimed Au(I) / Cu(I) polar-covalent bond: Proc. Natl. Acad. Sci. U.S.A., 2017, 114, E5042) – which becomes 87 kcal mol -1 with benzene stacking. Surprisingly, DFT analysis suggests that: (a) some other literature precedents should have attained a stacked M 3 @M' 3 product akin to the one herein, with similar or even higher binding energy; and (b) a high overall intertrimer bonding energy by inferior electrostatic assistance, underscoring genuine orbital overlap between M and M' frontier molecular orbitals in such polar-covalent M–M' bonds in this family of molecules. The Au / Ag bonding is reminiscent of classical Werner-type coordinate-covalent bonds such as H 3 N: → Ag in [Ag(NH 3 ) 2 ] + , as demonstrated herein quantitatively. Solidstate and molecular modeling illustrate electron flow from the p-basic gold trimer to the p-acidic silver trimer with augmented contributions from ligand-to-ligand’ (LL'CT) and metal-to-ligand (MLCT) charge transfer.

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

A Quantum Mechanical Study of Quartz (101) Interfacial Boundaries

Density functional theory (DFT) calculations were performed on periodic and molecular models to explore the energetics of bonding modes at quartz–quartz interfaces, common grain contacts in sandstones. Four interface types were modeled: H 2 O-mediated, silanol H-bond mediated, O–O peroxide bonds, and siloxane bonds. Each type of interface may exist in the subsurface of the Earth, depending upon water activity, temperature, and pressure, as interfacial structures are stable within DFT molecular dynamics simulations up to 473 K. The results predict interfacial energies SiOSi > SiOOSi > SiOH-6H 2 O > 2SiOH, as expected based on similar bond strengths in other systems. The proposed SiOOSi interface is not generally accounted for in chemomechanical models but could be important in sandstones at depth in the Earth. The implications for intergrain fracturing and H2 generation are discussed.

Hydrogen↗
DOE OSTI2022

The Crystalline Behavior and Device Function of Nonfullerene Acceptors in Organic Solar Cells

The current research investigates the structure features and intermolecular interactions of nonfullerene acceptors (NFAs) in single crystal and thin films, as well as their solar cell applications. Guiding parameters and key intermolecular forces that lead to 2D brickwork or 3D web packing are identified. The atomic modification is shown as the key to induce hydrogen bonding or π–π stacking column, which results in different crystalline packing. The molecular assembly in thin film is initiated by hydrogen bonding and completed by π–π stacking reorganization. The packing energy is seen as a guiding parameter that dictates the NFA crystalline morphology in blended thin films. The crystalline packing motif is not directly related with device efficiency. However, the crystalline morphology is the key parameter to influence exciton/carrier dynamics and device performance. In conclusion, a broader picture on the scaling behavior of organic semiconductor crystals ranging from oligoacenes to NFAs is established.

14 SOLAR ENERGY↗
DOE OSTI2025

PySIDT: Subgraph Isomorphic Decision Trees for Molecular Property Prediction

Accurate molecular property prediction is important across all fields of chemistry. Deep neural networks (DNNs) have become increasingly popular due to their ability to train automatically, avoiding the incredibly tedious process of constructing and extending traditional property estimation schemes. However, DNNs require large amounts of training data, are challenging to interpret, require large amounts of memory to load even during inference, and have severe difficulties incorporating qualitative chemical knowledge, which are often desired for molecular property prediction tasks. Here, in this study, we present PySIDT (https://github.com/zadorlab/PySIDT), a software for training and running inference on Subgraph Isomorphic Decision Trees (SIDTs). SIDTs are graph-based decision trees made of nodes associated with molecular substructures. Inference is done by descending target molecular structures down the decision tree to nodes with matching subgraph isomorphic substructures and making predictions based on the final (most specific) nodes matched. SIDTs scale down well to dataset sizes much smaller than is feasible for DNNs. As trees of molecular substructures, SIDTs are inherently readable and easy to visualize, making them easy to analyze. They are also straightforward to extend and retrain, facilitate uncertainty estimation, and enable easy integration of expert knowledge. We demonstrate the SIDT approach discussing its application to a diverse range of molecular prediction tasks: rate coefficient estimation, diffusion coefficient estimation, thermochemistry estimation, transition state bond stretch prediction, p K a prediction, stability of molecular structures, stability of surface structures, and prediction of surface lateral interaction energetics. Additionally, we demonstrate the power of the SIDT algorithms in two direct learning curve vanilla comparisons with the popular DNN-based software Chemprop and the popular gradient boosted trees-based software XGBoost on enthalpy of formation and rate coefficient prediction tasks. In particular, in the enthalpy of formation case, vanilla PySIDT is able to outperform vanilla Chemprop and XGBoost across the full range of training/validation set sizes out to 11,560 data points.

Johnson, Matthew Sean [Sandia National Laboratorie↗
DOE OSTI2023

Kβ X-ray Emission Spectroscopy of Cu(I)-Lytic Polysaccharide Monooxygenase: Direct Observation of the Frontier Molecular Orbital for H 2 O 2 Activation

Lytic polysaccharide monooxygenases (LPMOs) catalyze the degradation of recalcitrant carbohydrate polysaccharide substrates. These enzymes are characterized by a mononuclear Cu(I) active site with a three-coordinate T-shaped “His-brace” configuration including the N-terminal histidine and its amine group as ligands. Here, this study explicitly investigates the electronic structure of the d 10 Cu(I) active site in a LPMO using Kβ X-ray emission spectroscopy (XES). The lack of inversion symmetry in the His-brace site enables the 3d/p mixing required for intensity in the Kβ valence-to-core (VtC) XES spectrum of Cu(I)-LPMO. These Kβ XES data are correlated to density functional theory (DFT) calculations to define the bonding, and in particular, the frontier molecular orbital (FMO) of the Cu(I) site. These experimentally validated DFT calculations are used to evaluate the reaction coordinate for homolytic cleavage of the H 2 O 2 O–O bond and understand the contribution of this FMO to the low barrier of this reaction and how the geometric and electronic structure of the Cu(I)-LPMO site is activated for rapid reactivity with H 2 O 2 .

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

A linear surrogate for optimising functions of an orthogonal matrix with applications in wave function theory

The technique of surrogate optimisation is to use a simpler function to approximate a complex function that is time-consuming to evaluate. Here we show that the maximum of a special type of surrogate function f(U)=Tr(AU),UϵO(n) is at A T (AA T ) 1/2 , and that there is one and only one local maximum both in SO(n) and O(n)–SO(n). This function f(U) has been found to be useful in various aspects of electronic structure theory, including proving the Carlson-Keller theorem, and localising orbitals. As one other example, we apply it here to optimise the ground state of molecules using the Generalised Valence Bond wavefunction.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Distinguishing homolytic vs heterolytic bond dissociation of phenylsulfonium cations with localized active space methods

Modeling chemical reactions with quantum chemical methods is challenging when the electronic structure varies significantly throughout the reaction and when electronic excited states are involved. Multireference methods, such as complete active space self-consistent field (CASSCF), can handle these multiconfigurational situations. However, even if the size of the needed active space is affordable, in many cases, the active space does not change consistently from reactant to product, causing discontinuities in the potential energy surface. The localized active space SCF (LASSCF) is a cheaper alternative to CASSCF for strongly correlated systems with weakly correlated fragments. The method is used for the first time to study a chemical reaction, namely the bond dissociation of a mono-, di-, and triphenylsulfonium cation. LASSCF calculations generate smooth potential energy scans more easily than the corresponding, more computationally expensive CASSCF calculations while predicting similar bond dissociation energies. Furthermore, our calculations suggest a homolytic bond cleavage for di- and triphenylsulfonium and a heterolytic pathway for monophenylsulfonium.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗