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

Hydrogen-Bonded Organic Frameworks: A Rising Class of Porous Molecular Materials

Hydrogen-bonded organic frameworks (HOFs) are a class of porous molecular materials that rely on the assembly of organic building blocks by means of hydrogen-bonding interactions to form two-dimensional (2D) and three-dimensional (3D) crystalline networks. The reversible nature of the hydrogen-bond formation endows HOFs with the attributes of solution processability and simple regeneration. High-quality single crystals of HOFs can be grown easily for unambiguous superstructure determination by single-crystal X-ray diffraction, which is crucial for the elucidation of superstructure–property relationships. During the past decade, considerable progress has been achieved in realizing stable HOFs with permanent porosities by focusing on the design of molecular building blocks in order to introduce rigidity, auxiliary [π···π] interactions, and interpenetration of their frameworks to sustain the extended networks. The applications of HOFs are far-reaching, spanning catalysis, energy, and biomedical products as well as the storage and separation of fine chemicals. In this paper, we, first of all, provide an overview of the chronological development of HOFs, starting from the seminal work by Marsh and Duchamp in 1969 on the crystal superstructure of the hydrogen-bonded networks of trimesic acid. We identify the development of novel hydrogen-bonding motifs such as diaminotriazine (DTA), the introduction of the concept of molecular tectonics, and the establishment of permanent porosity in HOFs as being some of the milestones, which incentivized the current burgeoning research endeavors on developing HOFs as multifunctional materials. This Account is focused primarily on surveying the strategies for constructing porous 3D HOFs based on organic building blocks with peripheral carboxyl groups. These strategies are presented in the following categories: (1) the polycatenation of 2D networks by trigonal building blocks to form global 3D frameworks, (2) the utilization of building blocks with 3D geometries—tetrahedral and trigonal prismatic—that are predisposed to form 3D networks, and (3) the docking by shape-fitting of geometrically labile building blocks. We emphasize how the molecular geometry of the building blocks plays an important role in modulating the superstructures of extended frameworks so as to address specific applications. Recognizing that the in silico design of HOFs is the ultimate goal of researchers in this field, we also discuss the recent advances in superstructure prediction that lead to the formation of porous supramolecular crystals and assess the complications in implementing computational methods for HOFs with complex superstructures. We hope this Account will inspire the development of new supramolecular designs and creative approaches to crystal engineering that aid and abet the assembly of multifunctional HOFs with customizable properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Molecular Dynamics Investigation of Nanoscale Hydrophobicity of Polymer Surfaces: What Makes Water Wet?

The wettability of a polymer surface–related to its hydrophobicity or tendency to repel water–can be crucial for determining its utility, such as for a coating or a purification membrane. While wettability is commonly associated with the macroscopic measurement of a contact angle between surface, water, and air, the molecular physics that underlie these macroscopic observations are not fully known, and anticipating relative behavior of different polymers is challenging. To address this gap in molecular-level understanding, we use molecular dynamics simulations to investigate and contrast interactions of water with six chemically distinct polymers: polytetrafluoroethylene, polyethylene, polyvinyl chloride, poly(methyl methacrylate), Nylon 66, and polyvinyl alcohol. We show that several prospective quantitative metrics for hydrophobicity agree well with experimental contact angles. Moreover, the behavior of water in proximity to these polymer surfaces can be distinguished with analysis of interfacial water dynamics, extent of hydrogen bonding, and molecular orientation–even when macroscopic measures of hydrophobicity are similar. The predominant factor dictating wettability is found to be the extent of hydrogen bonding between polymer and water, but the precise manifestation of hydrogen bonding and its impact on surface-water structure varies. In the absence of hydrogen bonding, other molecular interactions and polymer mechanics control hydrophobic ordering. Furthermore, these results provide new insights into how polymer chemistry specifically impacts water-polymer interactions and translates to surface hydrophobicity. Such factors may facilitate the design or processing of polymer surfaces to achieve targeted wetting behavior, and presented analyses can be useful in studying the interfacial physics of other systems.

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

Dramatic effect of electrode type on tunnel junction based molecular spintronic devices

A new class of molecular spintronic devices can be fabricated by chemically bonding magnetic molecular channels to the electrodes of a prefabricated tunnel junction with exposed side edges. Prior experimental studies showed that the cyanide-bridged octametallic molecular cluster, [(pzTp)FeIII(CN)3]4[NiII(L)]4¬[O3SCF3]4 [(pzTp) = tetra(pyrazol-1-yl)borate; L = 1-S(acetyl)tris(pyrazolyl)decane] molecule impact depended on the type of metallic electrodes used in the tunnel junction testbed. Experimental magnetization and transport studies showed a dramatic difference in molecule response on tunnel junctions with different combinations of metallic electrodes. Transport via paramagnetic molecular channels on a tunnel junction involving paramagnetic and ferromagnetic metal electrodes was dramatically different than the suppressed current state observed on tunnel junctions involving two ferromagnetic electrodes. We conducted theoretical studies to understand the experimental data and explore a wide range of electrode materials on tunnel junction-based molecular spintronics devices (TJMSD). Here, we report a Monte Carlo simulation study that focuses on understanding the effect of electrodes on the magnetic and physical properties of TJMSD. A 3D Heisenberg model of cross-junction-shaped TJMSD was used for the simulation study. We studied the effects of ferromagnetic, paramagnetic, and antiferromagnetic electrode materials. This study provides insights for designing and understanding futuristic molecular spintronics devices.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Molecular Simulation of Covalent Bond Dynamics in Liquid Silicon

Many atomic liquids can form transient covalent bonds reminiscent of those in the corresponding solid states. These directional interactions dictate many important properties of the liquid state, necessitating a quantitative, atomic-scale understanding of bonding in these complex systems. A prototypical example is liquid silicon, wherein transient covalent bonds give rise to local tetrahedral order and consequent nontrivial effects on liquid-state thermodynamics and dynamics. To further understand covalent bonding in liquid silicon, and similar liquids, we present an ab initio-simulation-based approach for quantifying the structure and dynamics of covalent bonds in condensed phases. Through the examination of structural correlations among silicon nuclei and maximally localized Wannier function centers, we develop a geometric criterion for covalent bonds in liquid Si. We use this to monitor the dynamics of transient covalent bonding in the liquid state and estimate a covalent bond lifetime. We compare covalent bond dynamics to other processes in liquid Si and similar liquids and suggest experiments to measure the covalent bond lifetime.

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

In situ molecular imaging of ion clusters reveals the acid gas capture capacity and mechanism of water-lean ionic liquids

Water-lean solvents are a promising technology for capturing acid gases like carbon dioxide (CO 2 ). In situ liquid time-of-flight secondary ionization mass spectroscopy (ToF-SIMS) is used to study a representative solvent N-(2-ethoxyethyl)-3-morpholinopropan-1-amine (2-EEMPA) with different CO 2 loadings to reveal the complex solvent structure upon CO 2 capture. Characteristic peaks of 2-EEMPA, such as m/z – 215 C 11 H 23 N 2 O 2 – (deprotonated 2-EEMPA) and m/z + 217 C 11 H 25 N 2 O 2 + (protonated 2-EEMPA), are detected due to acid gas uptake. Also, solvent molecules and carboxylate ion pairs, such as m/z – 259 C 12 H 23 N 2 O 4 – [(deprotonated 2-EEMPA∙∙∙CO 2 )] and m/z + 261 C 12 H 25 N 2 O 4 + (protonated 2-EEMPA∙∙∙CO 2 ), are observed. Interestingly, more than one CO 2 molecule can be captured per each solvent molecule as evidenced in SIMS mass spectra, for example, m/z – 321 C 13 H 25 N 2 O 7 – [(deprotonated 2-EEMPA)∙∙∙2CO 2 ∙∙∙H 2 O], m/z + 305 C 13 H 25 N 2 O 4 + [(protonated 2-EEMPA)∙∙∙2CO 2 ], m/z – 389 C 17 H 29 N 2 O 8 – [(deprotonated 2-EEMPA)∙∙∙3CO 2 ∙∙∙3CH 2 ], and m/z + 373 C 16 H 25 N 2 O 8 + [(protonated 2-EEMPA)∙∙∙3CO 2 ∙∙∙2C]. However, the monomer of 2-EEMPA and CO 2 seems to be most prevalent. Furthermore, solvent clusters are detected in loaded solvents, for instance m/z + 433 C 22 H 49 N 4 O 4 + [(2-EEMPA)2∙∙∙H] and m/z + 646 [(2-EEMPA) 3 -2H], while capturing CO 2 at different amounts. Relative abundance of cluster ions provides a semi-qualitative venue to assess the free energies of gas capture energetics, indicating the relative stability trend within the same solvent system, previously impossible. These observed ion clusters are verified with molecular modeling, where dimer, trimer, and cluster ions are validated for their presence either due to weak molecular interactions or hydrogen bonds. In situ molecular imaging of ionic liquids and molecular modeling reveals that the acid gas capture mechanism by ionic liquids includes both physical adsorption and chemical bonding with multiple reaction pathways, engaging cluster formation and alteration of solvent structures.

Acid gas capture↗
DOE OSTI2024

Atomistic simulations to reveal HIP-bonding mechanisms of Al6061/Al6061

Molecular dynamics simulations were employed to understand the diffusion bonding process during hot isostatic pressing (HIP) of Al6061/Al6061 alloy. Simulations of the HIP process reveal atomistic phenomena that are difficult or unlikely to be observed experimentally and provide useful insights into the mechanism of diffusion and bonding. Here, the results reveal that at the start of the HIP process, a massive incursion of oxygen atoms occurs from the pre-existing γ-Al 2 O 3 to the 6061 region across the interphase interface. These oxygen atoms interact with the enriched Mg atom layer present at the existing γ-Al 2 O 3 and 6061 matrix to form a secondary complex Mg 2 Al 2 O 5 phase. Diffusion calculations also show that transport of atoms due to the applied pressure is 4–5 orders of magnitude higher than would occur in the absence of HIP conditions. The Mg 2 Al 2 O 5 phase also provides efficient pathways for the rapid transport of Mg atoms. Because of the higher diffusion coefficients observed for Mg within the phase, Mg atoms can move more swiftly compared to their diffusion within other phases such as γ-Al 2 O 3 . This accelerated mobility facilitates the rapid movement of Mg atoms across the interface, leading to changes in the local composition and the potential growth of the Mg 2 Al 2 O 5 phase.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A statistical mechanics framework for polymer chain scission, based on the concepts of distorted bond potential and asymptotic matching

To design increasingly tough, resilient, and fatigue-resistant elastomers and hydrogels, the relationship between controllable network parameters at the molecular level (bond type, non-uniform chain length, entanglement density, etc.) to macroscopic quantities that govern damage and failure must be established. Many of the most successful constitutive models for elastomers have been rooted in statistical mechanical treatments of polymer chains. Typically, such constitutive models have used variants of the freely jointed chain model with rigid links. However, since the free energy state of a polymer chain is dominated by enthalpic bond distortion effects as the chain approaches its rupture point, bond extensibility ought to be accounted for if the model is intended to capture chain rupture. To that end, a new bond potential is supplemented to the freely jointed chain model (as derived in the u FJC framework of Buche and Silberstein (2021) and Buche et al. (2022)), which we have extended to yield a tractable, closed-form model of single chain behavior that should be amenable to continuum-level constitutive model development. Inspired by the asymptotically matched u FJC model response in both the low/intermediate chain force and high chain force regimes, a simple, quasi-polynomial bond potential energy function is derived. This bond potential exhibits harmonic behavior near the equilibrium state and anharmonic behavior for large bond stretches tending to a characteristic energy plateau (akin to the Lennard-Jones and Morse bond potentials). Using this bond potential, approximate yet highly-accurate analytical functions for bond stretch and chain force dependent upon chain stretch are established. Then, using this polymer chain model, a stochastic thermal fluctuation-driven chain rupture framework is developed. This framework is based upon a force-modified tilted bond potential that accounts for distortional bond potential energy, allowing for the derivation and subsequent calculation of the dissipated chain scission energy. Here, the cases of rate-dependent and rate-independent scission are accounted for throughout the rupture framework. The impact of Kuhn segment number on chain rupture behavior is also investigated. The model is fit to single chain mechanical response data collected from atomic force microscopy tensile tests for validation and to glean deeper insight into the molecular physics taking place. Due to their analytical nature, this polymer chain model and the associated rupture framework can in the future be implemented in finite element models accounting for fracture and fatigue in polydisperse elastomer networks.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electrochemically Mediated Au–C(sp 2 ) Anchors for Molecular Electronics

Terminal anchor groups play a key role in the stability and electronic properties of molecular junctions. Single molecule junctions typically consist of two preinstalled terminal anchors linking organic molecules to metal electrodes. Here, in this work, we show that p -terphenyl derivatives containing only a single terminal anchor show conductance features similar to junctions with two preinstalled terminal anchors. A set of p -terphenyl derivatives with one terminal anchor was prepared using automated chemical synthesis and characterized using single molecule electronics experiments, molecular dynamics (MD) simulations, bulk electrochemistry and spectroscopy, and nonequilibrium Green’s function-density functional theory (NEGF-DFT) calculations. Our results show that 4-amino- p -terphenyl (PPP) and related analogs exhibit a well-defined high conductance state that is diminished or absent in other p-terphenyl derivatives lacking a preinstalled amine terminal anchor or fluorine or methyl substitutions at the terminal para position. However, a low conductance state is observed in all amino- p -terphenyl derivatives with one preinstalled anchor due to molecular junctions formed by noncovalent dimeric π–π stacking interactions. The observed high conductance state diminishes upon the addition of reducing agents and is restored upon the addition of an oxidizing agent. Our results suggest that the high conductance state arises due to Au–C(sp 2 ) bond formation facilitated by a single electron oxidation event at the electrode surface. A series of control experiments with different anchor groups shows that primary amines play a key role in forming Au–C bonds for molecular junctions. Overall, these results suggest that Au–C bond formation gives rise to high conductance pathways in organic molecules containing only one preinstalled terminal anchor. Insights from this work can be leveraged in the design of molecular electronic devices, particularly in understanding the mechanisms of molecular binding and junction formation.

charge transport↗
DOE OSTI2021

Probing molecular vibrations by monochromated electron microscopy

Chemical bonds fundamentally determine molecular properties and are prevalently characterized by various spectroscopic means such as infrared and Raman spectroscopies. However, the spatial resolution of these conventional approaches is insufficient to reveal nanoscale features. Recently, monochromated electron energy-loss spectroscopy (EELS) in the transmission electron microscope achieved a groundbreaking energy resolution of a few millielectronvolts and enabled direct observation of molecular vibrational spectrum with unmatched spatial resolution. Vibrational EELS is widely applicable to both organic and inorganic matter in the solid state or liquid phase. In this study, we introduce recent advancements and key concepts of this method, compare with other spectroscopic techniques, and discuss future developments for potential applications in research fields centered on catalysts, polymers, and live cells.

Vibrational Spectroscopy↗
DOE OSTI2023

Scrutinizing formally Ni IV centers through the lenses of core spectroscopy, molecular orbital theory, and valence bond theory

Nickel K- and L 2,3 -edge X-ray absorption spectra (XAS) are discussed for 16 complexes and complex ions with nickel centers spanning a range of formal oxidation states from II to IV. K-edge XAS alone is shown to be an ambiguous metric of physical oxidation state for these Ni complexes. Meanwhile, L 2,3 -edge XAS reveals that the physical d-counts of the formally Ni IV compounds measured lie well above the d 6 count implied by the oxidation state formalism. The generality of this phenomenon is explored computationally by scrutinizing 8 additional complexes. The extreme case of NiF 6 2- is considered using high-level molecular orbital approaches as well as advanced valence bond methods. The emergent electronic structure picture reveals that even highly electronegative F-donors are incapable of supporting a physical d 6 Ni IV center. The reactivity of Ni IV complexes is then discussed, highlighting the dominant role of the ligands in this chemistry over that of the metal centers.

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

Structural determinants of protocadherin-15 mechanics and function in hearing and balance perception

The vertebrate inner ear, responsible for hearing and balance, is able to sense minute mechanical stimuli originating from an extraordinarily broad range of sound frequencies and intensities or from head movements. Integral to these processes is the tip-link protein complex, which conveys force to open the inner-ear transduction channels that mediate sensory perception. Protocadherin-15 and cadherin-23, two atypically large cadherins with 11 and 27 extracellular cadherin (EC) repeats, are involved in deafness and balance disorders and assemble as parallel homodimers that interact to form the tip link. Here we report the X-ray crystal structure of a protocadherin-15 + cadherin-23 heterotetrameric complex at 2.9-Å resolution, depicting a parallel homodimer of protocadherin-15 EC1-3 molecules forming an antiparallel complex with two cadherin-23 EC1-2 molecules. In addition, we report structures for 10 protocadherin-15 fragments used to build complete high-resolution models of the monomeric protocadherin-15 ectodomain. Molecular dynamics simulations and validated crystal contacts are used to propose models for the complete extracellular protocadherin-15 parallel homodimer and the tip-link bond. Steered molecular dynamics simulations of these models suggest conditions in which a structurally diverse and multimodal protocadherin-15 ectodomain can act as a stiff or soft gating spring. These results reveal the structural determinants of tip-link–mediated inner-ear sensory perception and elucidate protocadherin-15’s structural and adhesive properties relevant in disease.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Prospects for Forming C–N Bonds from Dinitrogen

The formation of C–N bonds from molecular dinitrogen (N 2 ) offers a synthetic route to value-added nitrogen-containing compounds without relying on prefunctionalized nitrogen sources. Here, this Perspective highlights recent advances in forming C–N bonds from N 2 with homogeneous transition metal complexes. After discussing how transition metal complexes activate N 2 through various coordination modes, we focus on the reactivity of reduced nitrogen intermediates with different kinds of carbon sources. Carbon electrophiles and nucleophiles enable C–N bond formation via insertion, substitution, or radical pathways. Cycloaddition reactions, particularly involving polarized N 2 ligands and unsaturated carbon electrophiles, offer routes to more complex products. We describe efforts to achieve catalytic turnover and emphasize the remaining obstacles to catalytic C–N bond construction.

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

Chemical bond effects in classical site density functional theory of inhomogeneous molecular liquids

Intra-molecular interactions or chemical bonds represent one of the main distinguishing characteristic of molecular fluids. Development of accurate and practical methods to treat these effects is one of the long standing problems in classical site density functional theory (SDFT). One particular instance when these issues become particularly severe is the case of classical interactions potentials with auxiliary sites or dummy atoms. In this situation current SDFT implementations, such as three-dimensional reference interaction site model (3D-RISM), lead to nonphysical results. We re-examine this issue in this work using our recent reformulation of SDFT. We put forward a simple practical solution to this problem, and illustrate its utility for the case of spherical solutes in diatomic liquids.

Chuev, Gennady N.↗
DOE OSTI2024

Toward sustainable production of N-containing products via nonthermal plasma-enhanced conversion of natural gas resources

Nonthermal plasmas can directly activate and cleave the strong chemical bonds in molecular nitrogen and methane to facilitate the transformation of these inherently stable molecules through the application of electrical energy. Here, we report a low temperature, atmospheric pressure, nonthermal plasma for the “one-pot” synthesis of olefins, alkynes, higher molecular weight hydrocarbons, ammonia, and nitrogen-containing liquids from a representative shale gas feed enriched with nitrogen. Furthermore, we reproducibly observe a wide range of valuable and synthetically challenging gas-phase and liquid-phase products containing C-N, C-C, and N-H bonding by controlling the N 2 concentration in the inlet feed stream. In nitrogen-lean regimes, hydrocarbon products dominate (e.g., ethylene, acetylene, etc.), while nitrogen-rich regimes promote incorporation of nitrogen into the products, leading to the formation of ammonia and liquid products containing a variety of functionalities (e.g., nitriles, amines, heterocycles). High resolution electrospray ionization mass spectrometry was used to measure molecular weights and identify the chemical formulas of the liquid products. Van Krevelen diagrams were created and showed many products with compositions around H/C = 2 and N/C = 0.5, indicating the potential importance of intermediate species with these ratios for liquid formation (e.g., CH 3 CN + H).

03 NATURAL GAS↗
DOE OSTI2026

Cosolvent-tuned interactions in ionic liquids: A vibrational and quantum-chemical study of ethylene glycol ratio effects

Ionic liquids (ILs) are attractive media for CO 2 capture but remain limited by viscosity and cost. Blending ILs with ethylene glycol (EG) is a practical route to mitigate these constraints, yet the molecular origins of cosolvent effects and their dependence on composition are not well resolved. We combine Fourier-transform infrared (FT-IR) spectroscopy with quantum-chemical (DFT) analysis to elucidate how the IL:EG molar ratio modulates intermolecular interactions and electronic structure. Computed vibrational frequencies enable mode assignment and deconvolution of overlapping bands, revealing systematic, ratio-dependent shifts and broadenings in (i) EG O–H stretching, (ii) cation and EG C–H stretchings (imidazolium C2–H, C4–H, C5–H, methyl and ethyl groups, -CH2 of EG), (iii) anion signature modes (e.g., CN motifs), and (iv) EG C–O and C–C stretchings, consistent with the redistribution of hydrogen-bonding networks. Molecular electrostatic potential (MESP) maps quantify attenuation of extreme potential regions with increasing EG, indicating progressive screening of cation–anion electrostatic interactions. Quantum Theory of Atoms in Molecules (QTAIM) identifies emergent bond critical points between EG and the IL ions, while Reduced Density Gradient–Noncovalent Interaction (RDG–NCI) analysis differentiates strong directional hydrogen bonds from dispersive contacts across compositions. Together, these results show that EG fraction controls a switch from predominantly ion–ion to mixed ion–EG coordination, altering local polarity and polarizability that underlie the observed FT-IR trends. The framework provides composition–structure–spectrum relationships that can guide rational selection of IL:EG ratios to balance favorable molecular interactions with practical performance targets in scalable CO 2 capture systems.

DAC↗
DOE OSTI2023

Connectivity between Static Field and Continuous Wave Field Effects on Excitation-Induced H 2 Activation

Due to the tremendous applications of the plasmon resonance excitation process, such as improvements in catalytic efficiency due to plasmonic enhancement and/or hot-electron processes, understanding the mechanism behind these processes has become a popular topic in recent years. In this work, we focus on unraveling the mechanism of excitation-induced H 2 activation using a simplified triangular Au 6 /Ag 6 cluster to investigate the effects of the electric field on electron redistribution and bond activation. We applied both static and continuous wave fields to investigate how these fields affect the systems. Geometrical changes (such as bond lengthening), molecular orbital reordering (affecting the relative energies of orbitals corresponding to hot-electron and charge-transfer excited states), and electronic charge redistribution between the cluster and the adsorbate occur upon application of a static electric field. To study H 2 activation, we apply Ehrenfest dynamics with real-time time-dependent density functional theory and examine how different excitation frequencies and polarizations affect bond activation. Moreover, electron-only dynamics are examined with real-time time-dependent density functional theory, and the time-dependent variations in the orbital populations and electronic transitions provide information about the excitation and relaxation processes of hot electrons with applied electric fields. Here, the static field results represent structures that can be accessed during the evolution of the systems when applying continuous wave fields. Through these studies, the effects of static and continuous wave field effects on plasmon-induced H 2 activation can be understood.

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

Minimal Experimental Bias on the Hydrogen Bond Greatly Improves Ab Initio Molecular Dynamics Simulations of Water

Experiment Directed Simulations (EDS) is a method within a class of techniques seeking to improve molecular simulations by minimally biasing the system Hamiltonian to reproduce certain experimental observables. In a previous application of EDS to ab initio molecular dynamics (AIMD) simulation based on electronic density functional theory (DFT), the AIMD simulations of water were biased to reproduce its experimentally derived solvation structure. In particular, by solely biasing the O-O pair correlation functions, other structural and dynamical properties that were not biased were improved. In this work, the hypothesis is tested that directly biasing the OH pair correlation (and hence the H-O∙∙∙H hydrogen bonding), will provide an even better improvement of DFT-based water properties in AIMD simulations. The logic behind this hypothesis is that for most electronic DFT descriptions of water the hydrogen bonding is known to be deficient due to anomalous charge transfer and over polarization in the DFT. Using recent advances to the EDS learning algorithm, we thus train a minimal bias on AIMD water that reproduces the O-H radial distribution function derived from the highly accurate MB-pol model of water. Finally, it is then confirmed that biasing the O-H pair correlation alone can lead to improved AIMD water properties, with structural and dynamical properties in even closer to experiment than the previous EDS-AIMD model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗