Data-driven molecular modeling with the generalized Langevin equation
Not Available
SEARCH · Engineering Papers
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.
Not Available
This work explores the effect of co-mingled dopants, Co(II) and Cr(III), on Tc(IV) incorporation and retention in magnetite when heat treated to 625 or 700 °C. Key trends in Tc retention in the high temperature regime were identified using a combination of density-functional-theory based ab initio molecular dynamics (AIMD) simulations, and batch experiments including solid phase characterization techniques, e.g. X-ray absorption spectroscopy. A stabilizing effect on Tc(IV) was observed when the number of Tc and Cr atoms are equal or when the magnetite surface is oversaturated with Tc and Cr inclusions. Here, oversaturation is hypothesized to force Cr from the magnetite surface to form a Cr2O3 phase, which may act as a protective layer that prevents Tc release. With the addition of Co, Tc(IV) is stabilized via redox processes. The presence of Cr in low concentrations interferes with this redox stabilization and Cr is preferentially stabilized as opposed to Tc. As a result, using Co as a stabilizing dopant for Tc may be compromised in the presence of Cr. When the relative concentration of Tc, Cr and Co is the same, or more Co atoms are added to high Cr incorporated systems, the formation of the Cr2O3 phase may be suppressed. Although waste streams with co-mingled Tc and Cr potentially may benefit from a Co dopant, since the formation of a Cr2O3 passivation layer may protect incorporated Tc from being released, the relative concentration of the three elements will be a critical parameter for maximizing effectiveness of this strategy.
Here, homogeneous and heterogeneous catalysts have advantages and disadvantages. Molecular catalysts supported on heterogeneous materials have the potential to combine the advantages of both fields. Jabbour et al. report a detailed study of the 3D structure of an Ir-catalyst on silica.
Not Available
Noncovalent intermolecular interactions are very important in many research areas. Therefore, it is vital to understand the extent to which approximate density functionals give a proper description of noncovalent interactions. Previous research has demonstrated that some approximate density functionals can predict usefully accurate interaction energies for many noncovalent systems; however, most of that work is limited to small and moderate-sized molecules. Very recently though, accurate benchmarks have become available for some very large molecules. Here, the present work applies 21 approximate density functionals to compute the binding energies of seven large molecular systems that have a number of atoms ranging from 200 to 910. The results are judged by comparison to the recently published CIM-DLPNO-CCSD(T) results, which are assumed to provide a reliable benchmark. The five most accurate methods among those tested are found to be PW6B95-D4, PW6B95-D3(BJ), revM11, M06-L, and MN15.
Performing reliable computer simulations of elementary processes occurring at metal–water interfaces is pivotal for novel catalyst design in sustainable energy applications. Computational catalyst design hinges on the ability to reliably and efficiently compute the potential energy surface (PES) of the system. Here, due to the large system sizes needed for studying processes at liquid water–metal interfaces, these systems can currently not be described using density functional theory (DFT). In this work, we used a hybrid quantum mechanical, molecular mechanical, and machine learning potential for studying the adsorption behavior of phenol, atomic hydrogen, 2-butanol, and 2-butanone on the (0001) facet of Ru under reducing conditions when Ru is not oxidized. Specifically, we describe the adsorbate and the surrounding metal atoms at the DFT level of theory. Here, we also considered the electrostatic field effect of the water molecules on adsorbate–metal interactions. Next, for the water–water and water–adsorbate interactions, we used established classical force fields. Finally, for the water–Ru surface interaction, for which no reliable force fields have been published, we used Behler–Parrinello high-dimensional neural network potentials (HDNNPs). Employing this setup, we used our explicit solvation for metal surface (eSMS) approach to compute the aqueous-phase effect on the low-coverage adsorption of selected molecules and atoms on the (0001) facet of Ru. In agreement with previous experimental and computational studies of oxygenated molecules over transition metal facets, we found that liquid water destabilizes the tested adsorbates on Ru(0001). Interestingly, our findings indicate that adsorbates on Ru are less affected by the presence of an aqueous phase than on other transition metals (e.g., Pt), highlighting the necessity of experimental investigations of Ru-based catalytic systems in liquid water.
The carbonate minerals of Ca and Mg are abundant throughout the lithosphere and have recently garnered significant research interest as possible long-term carbon sinks in the sequestration of atmospheric carbon dioxide. Nonetheless, an understanding of the atomic-level processes comprising their mineralization remains limited. Furthermore, we characterize and contrast the mechanisms of contact ion-pair formation in aqueous Ca and Mg carbonate systems, which represents the most fundamental step leading to the formation of their mineral solids. Utilizing multilevel embedded correlated wavefunction-based ab initio molecular dynamics/quantum mechanics simulations, we characterize not only the dynamics of these processes but also factors arising from the electronic structure of the involved species, revealing further details of the fundamentally different mechanisms for the interconversion between the contact ion-pairs and solvent-shared ion-pairs of Ca versus Mg carbonate.
For this work, we have investigated the structure and dynamics of water (D 2 O) adsorbed on TiO 2 (110) for coverages between 0 and 1 monolayer with infrared reflection absorption spectroscopy (IRAS) and ab initio molecular dynamics (AIMD) simulations. For D 2 O coverages as low as 0.4 monolayers (ML) on a hydroxylated surface, IR spectra typical of hydrogen-bonded chains of water molecules are observed. However, for D 2 O coverages ≥ 0.3 ML, a sharp, high frequency peak is also observed in the p-polarized spectra that is red-shifted relative to the bridging hydroxyl peak. This new peak is not observed for water adsorbed on an oxidized surface. Based on the AIMD simulations and comparisons with previous IR spectra for TiO 2 nanoparticles, the new peak is assigned to terminal hydroxyl groups produced by dissociative adsorption of some of the water on TiO 2 (110). The simulations indicate that water dissociation is related to the presence of defect electrons in the system, but not due to direct interactions between adsorbed water and bridging hydroxyls.
Here, the present work demonstrates the generality of molecular scaling relationships toward achieving faster catalysis at lower overpotentials. In iron porphyrin-catalyzed oxygen reduction, independent tuning of the catalyst E 1/2 and buffer pK a resulted in improvement of both the maximum turnover frequency (TOF max ) and the effective overpotential (η eff ). Interestingly, the most optimal system is achieved with the most reducing catalyst and a weak buffered acid—a perhaps counterintuitive conclusion. More generally, this study illustrates a path to optimization of catalytic systems through their differing responses to changes in the intrinsic features of the system.
Reactive organic fluid - mineral interactions at elevated temperatures contribute to the evolution of planetary matter. One of the less studied but important transformations in this regard involves the reactions of formic acid with naturally occurring clays such as sodium montmorillonite. To advance a mechanistic understanding of these interactions, we use ReaxFF reactive molecular dynamics simulations in conjunction with infrared (IR) spectroscopy and X-ray scattering experiments to investigate the speciation behavior of water-formic acid mixtures on sodium montmorillonite interfaces at 473 K and 1 atm. Using a newly developed reactive forcefield, we show that the experimental IR spectra of unreacted and reacted mixture can be accurately reproduced by ReaxFF/MD. We further benchmark the simulation predictions of sodium carbonate and bicarbonate formation in the clay interlayers using Small and Wide-Angle X-ray Scattering measurements. Subsequently, leveraging the benchmarked forcefield, we interrogate the pathway of speciation reactions with emphasis on carbonate, formate, and hydroxide groups elucidating the energetics, transition states, intermediates, and preferred products. Further, we also delineate the differences in reactivities and catalytic effects of clay edges, facets, and interlayers owing to their local chemical environments, which have far reaching consequences in their speciation behavior. The experimental and simulation approaches described in this study and the transferable forcefields can be applied translationally to advance the science of clay-fluid interactions for several applications including subsurface fluid storage and recovery and clay-pollutant dynamics.
Methyl-terminated p-type silicon photoelectrodes selectively drive CO 2 reduction by a homogeneous [Ru(tpy)(Mebim-py)(NCCH 3 )] 2+ catalyst (tpy = 2,2′:6′,2″-terpyridine, Mebim-py = 1-methylbenzimidazol-2-ylidene-3-(2′-pyridine)). A 460 mV photovoltage is quantified for the photoelectrode. Under 1 sun illumination, this system achieves a Faradaic efficiency of 87% for CO at −1.7 V vs Fc +/0 , matching reports of the same catalyst at metallic electrodes operating at −2.1 V. When 5% water is introduced, the CH 3 -terminated Si photoelectrode remains stable, selectivity for CO is retained, and current density increases. Methyl termination suppresses the competitive hydrogen evolution observed for H-terminated Si photoelectrodes, which under the same conditions produce ca. 60% CO and 8% H 2 and have unstable performance. Furthermore, these results establish that a semiconductor photoelectrode can power a molecular CO 2 reduction catalyst without hydrogen evolution by the photoelectrode itself. Methyl termination of p-Si allows CO 2 reduction to kinetically outcompete proton reduction, revealing an important design principle for selective fuel formation.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The effect of molecular weight (M) on the fragmentation kinetics of micelles formed by 1,2-polybutadiene-block-poly(ethylene oxide) (PB-PEO) copolymers was studied in the ionic liquid 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide. A series of six samples, with total M ranging from 10 4 to 10 5 g mol –1 and nearly constant composition (f PEO ≈ 0.4), were examined; all six formed spherical micelles with PEO coronas. Nonequilibrium PB-PEO micelles were prepared by direct dissolution, a process that systematically produces nanoparticles with mean aggregation numbers more than twice the equilibrium values. When subjected to high temperature annealing (170 °C), the average micelle radius was found to decrease substantially, as determined by temperature-jump dynamic light scattering (T-jump DLS) and time-resolved small-angle X-ray scattering (TR-SAXS). Furthermore, the characteristic fragmentation times (τ) were found to increase strongly with increasing degree of polymerization N, as τ ~ N 1.8 . This result compares favorably with the prediction of a previously untested model.
The nature of electron transfer across metal oxide–water interfaces depends significantly on the band gap of the oxide and its band edge energies relative to the potentials of relevant aqueous redox couples. Here we focus on the water interface with MgO, a prototypical wide band gap oxide whose conduction band edge is close in energy to that of water. We investigate the behavior of an excess electron at and out of equilibrium near the interface using ab initio molecular dynamics based on hybrid density functional theory. Our simulations show that under equilibrium conditions the excess electron (donated by an Al impurity in MgO) localizes to a midgap defect state comparable in energy and shape to a hydrated electron in bulk water. To characterize the electron transfer from the conduction band of MgO to interfacial product states, we dope near-equilibrium configurations of the pristine MgO–water system with Al and run short trajectories of these instantaneously out-of-equilibrium systems. Here we observe two distinct products associated with the excess electron: a surface-localized electron (e surf –) and an aqueous hydrogen radical (H • ). The H • pathway exhibits a much higher activation barrier despite being more exoergic, making e surf – the kinetic product. Our characterization of the pathways on the basis of Marcus theory is consistent with the poor observed utility of MgO for water radiolysis. Moreover, we anticipate that the computational framework employed here will be broadly applicable to assessing electron transfer mechanisms at aqueous, photocatalytic interfaces.
Explore the source record for details and available documents.
Abstract Lung cancer is frequently caused by activating mutations in the epidermal growth factor receptor (EGFR). Allosteric EGFR inhibitors offer promise as the next generation of therapeutics, as they are unaffected by common ATP-site resistance mutations and synergize with the drug osimertinib. Here, we examine combinations of ATP-competitive and allosteric inhibitors to better understand the molecular basis for synergy. We identify a subset of irreversible EGFR inhibitors that display positive binding cooperativity and synergy with the allosteric inhibitor JBJ-04-125-02 in several EGFR variants. Structural analysis of these complexes reveals conformational changes occur mainly in the phosphate-binding loop (P-loop). Mutation of F723 in the P-loop reduces cooperative binding and synergy, supporting a mechanism in which F723-mediated contacts between the P-loop and the allosteric inhibitor are critical for synergy. These structural and mechanistic insights will aid in the identification and development of additional inhibitor combinations with potential clinical value.
The blood clotting protein fibrin has extraordinary elastomeric properties due to its hierarchical structure. SAXS is combined with computational molecular modeling, providing insight in fibrin elasticity and guidelines for designing new polymers.