Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SOLVATION”

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 1,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Low Adsorption of Magnetite Nanoparticles with Uniform Polyelectrolyte Coatings in Concentrated Brine on Model Silica and Sandstone

In subsurface imaging and oil recovery where temperatures and salinities are high, it is challenging to design polymer-coated nanoparticles with low retention (high mobility) in porous rock. Herein, the grafting of poly(2-acrylamido-2-methyl-1-propanesulfonic acid-$co$-acrylic acid) (poly(AMPS-$co$-AA)) on magnetic iron oxide nanoparticles was sufficiently uniform to achieve low adsorption on model colloidal silica and crushed Berea sandstone in highly concentrated API brine (8% NaCl and 2% CaCl 2 by weight). The polymer shell was grafted via amide bonds to an aminosilica layer, which was grown on silica-coated magnetite nanoparticles. The particles were found to be stable against aggregation in American Petroleum Institute (API) brine at 90 °C for 24 h. For IO nanoparticles with ~23% polymer content, Langmuir adsorption capacities on colloidal silica and crushed Berea Sandstone in batch experiments were extremely low at only 0.07 and 0.09 mg of IO/m 2 , respectively. Furthermore, upon injection of a 2.5 mg/mL IO suspension in API brine in a column packed with crushed Berea sandstone, the dynamic adsorption of IO nanoparticles was only 0.05 ± 0.01 mg/m 2 , which is consistent with the batch experiment results. The uniformity and high concentration of solvated poly(AMPS-$co$-AA) chains on the IO surfaces provided electrosteric stabilization of the nanoparticle dispersions and also weakened the interactions of the nanoparticles with negatively charged silica and sandstone surfaces despite the very large salinities.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Machine Learning-Driven Solvent Screening for Biobased 2,3-Butanediol Extraction

Biobased 2,3-butanediol (2,3-BDO) is a valuable biomass-derived chemical due to its versatility in being transformed into a wide variety of products. However, the separation and purification of 2,3-BDO from fermentation broth remain a significant challenge owing to its high boiling point and hydrophilic nature. Herein, we developed a machine learning (ML)-based screening workflow that uses molecular calculations as training data and requires only a small number of experimental measurements for validation to identify alternative solvent candidates for the liquid–liquid extraction (LLE) of 2,3-BDO from aqueous solution. In particular, 130 density functional theory (DFT) calculations with the implicit solvation method not only built a correlation between the computational partition coefficient and the experimental distribution coefficient of 2,3-BDO but also parameterized an Extra-Trees ML model to screen the distribution coefficient for a wider range of 6717 organic solvents. The experimental measurements of only 24 solvents were needed to validate the computational results. A list of 50 prioritized solvents was proposed for 2,3-BDO LLE, and seven additional experimental measurements were conducted to further verify our selected solvents. The impact of the extraction temperature and solvent-to-feed ratio was also investigated for selected solvents in experiments. Furthermore, this work suggested alternative solvents for 2,3-BDO LLE and proposed a versatile workflow that requires fewer experiments and can be applied to a broader range of LLE studies.

Extraction↗
DOE OSTI2017

Hexavalent Actinide Extraction Using N , N -Dialkyl Amides

Recently, efforts towards closing the nuclear fuel cycle have considered oxidizing americium to a hexavalent, linear-dioxo cation for co-recovery with hexavalent U, Np, and Pu from molar nitric acid using solvating extraction ligands such as tri-n-butyl phosphate or diamyl amylphosphonate. This work assesses solvent extraction recovery of sodium bismuthate oxidized americium by N , N -di-(2-ethylhexyl)butyramide (DEHBA), N , N -di-(2-ethylhexyl)isobutyramide (DEHiBA), and N , N -dihexyloctanamide (DHOA). Extraction efficiency between the monoamides was found to increase in the order of DEHiBA < DHOA < DEHBA. For all monoamides, oxidized americium extraction was less than 50% from 4 M HNO 3 and below. Extraction efficiency above 50% was obtained using concentrations of 5 M HNO 3 or higher. The DEHBA extractant provided the highest distribution value of 5.4 at 7 M HNO 3 . Distribution values were found to be stable for up to 45 seconds aqueous/organic phase mixing times and indicated decreased reduction of hexavalent americium relative to separations completed with organophosphorus extractants. The simultaneous co-extraction of U, Np, Pu, and Am was demonstrated using DEHiBA, and was found to decrease with increasing atomic number (D U > D Np > D Pu > D Am ). Interestingly, a break in recovery was observed where ligher actinides, U and Np, were better recovered relative to the heavier actinides, Pu and Am, in this study. This observation seems to be related to differences in the extracted metal complex for light actinides compared to heavy actinides. In conclusion, the metal:ligand stoichiometry was found to be mixed 1:1 and 1:2 species for U and Np while Pu and Am displayed 1:2 stoichiometry.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Molybdenum(VI) Coordination in Tributyl Phosphate Chloride Based System

Fundamental coordination chemistry of Mo(VI) at its macroconcentrations in solvent extraction systems is of great importance for industrial processes that require the purification or recovery of large concentrations of Mo. The coordination of Mo(VI) in tri-n-butyl phosphate (TBP) from solutions of hydrochloric acid and up to 0.3 M Mo was investigated using UV, FTIR, and P-31 NMR spectroscopies, as well as EXAFS. From these techniques we resolved near-neighbor atoms, speciation, structural information on the coordination environment, and thermodynamic parameters affiliated with the solvent extraction of Mo(VI) and HCl. The solvated extracted form of Mo(VI) as MoO2Cl2 center dot 2TBP was identified. High extraction yield of Mo at >5 M HCl concentration is driven by replacing HCl in the organic phase by Mo. The existence of additional organic Mo adducts is also discussed with the aid of density functional theory, however no evidence of dimeric or polymeric Mo species was found to be present in TBP.

Tkac, Peter↗
DOE OSTI2022

In Situ Visualization of Local Distortions in the High- T c Molecule-Intercalated Li x (C 5 H 5 N) y Fe 2– z Se 2 Superconductor

A time-resolved synchrotron X-ray total scattering study sheds light on the evolution of the different structural length scales involved during the intercalation of the layered iron–selenide host by organic molecular donors, aiming at the formation of the expanded-lattice Li x (C 5 H 5 N) y Fe 2– z Se 2 hybrid superconductor. The intercalates are found to crystallize in the tetragonal ThCr 2 Si 2 -type structure at the average level, however, with an enhanced interlayer iron–selenide spacing ( d = 16.2 Å) that accommodates the heterocyclic molecular spacers. Quantitative atomic pair distribution function (PDF) analysis at variable times suggests distorted FeSe 4 tetrahedral local environments that appear swollen with respect to those in the parent β-FeSe. Simultaneously acquired in situ synchrotron X-ray powder diffraction data disclose that secondary phases (α-Fe and Li 2 Se) grow significantly when a higher lithium concentration is used in the solvothermal reaction or when the solution is aged. These observations are in line with the strongly reducing character of the intercalation medium’s solvated electrons that mediate the defect chemistry of the expanded-lattice superconductor. In the latter, intralayer correlated local distortions indicate electron-donating aspects that reflect in somewhat enlarged Fe–Se bonds. Additionally, they also reveal a degree of relief of chemical pressure associated with a large distance between Fe and Se sheets (“taller” anion height) and a stretched Fe–Fe square planar topology. The elongation of the latter, derived from the in situ PDF study, speaks for a plausible increase in the Fe-site vacancy concentration. The evolution of the local structural parameters suggests an optimum reaction window where kinetically stabilized phases resemble the distortions of the edge-sharing Fe–Se tetrahedra, required for a high- T c in expanded-lattice iron-chalcogenides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

An Aluminum-Based Metal–Organic Cage for Cesium Capture

Metal–organic cages are a class of supramolecular structures that often require the careful selection of organic linkers and metal nodes. Of this class, few examples of metal–organic cages exist where the nodes are composed of main group metals. Herein, we have prepared an aluminum-based metal–organic cage, H 8 [Al 8 (pdc) 8 (OAc) 8 O 4 ] (Al-pdc-AA), using inexpensive and commercially available materials. The cage formation was achieved via solvothermal self-assembly of solvated aluminum and pyridine-dicarboxylic linkers in the presence of a capping agent, acetic acid. The obtained supramolecular structure was characterized by single-crystal X-ray diffraction (SCXRD), thermogravimetric analysis, and NMR spectroscopy. Based on crystal structure and computational analyses, the cage has a 3.7 Å diameter electron-rich cavity suitable for the binding of cations such as cesium (ionic radius of 1.69 Å). Here, the host–guest interactions were probed with 1 H and 133 Cs NMR spectroscopy in DMSO, where at low concentrations, Cs + binds to Al-pdc-AA in a 1:1 ratio. The binding site was identified from the crystal structure of CsH 7 [Al 8 (pdc) 8 (OAc) 8 O 4 ] (Cs + Al-pdc-AA), and a binding affinity of ~10 6 –10 7 M –1 was determined from NMR titration experiments. The Al-pdc-AA showed improved selectivity for cesium binding over alkali metal cations (Cs + > Rb + > K + >> Na + ~ Li + ). Collectively, the study reports a novel aluminum cage that can serve as a promising host for efficient and selective cesium removal.

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

Solution Structures of Europium Terpyridyl Complexes with Nitrate and Triflate Counterions in Acetonitrile

Lanthanide–ligand complexes are key components of technological applications, and their properties depend on their structures in the solution phase, which are challenging to resolve experimentally or computationally. The coordination structure of the Eu 3+ ion in different coordination environments in acetonitrile is examined using ab initio molecular dynamics (AIMD) simulations and extended X-ray absorption fine structure (EXAFS) spectroscopy. AIMD simulations are conducted for the solvated Eu 3+ ion in acetonitrile, both with or without a terpyridyl ligand, and in the presence of either triflate or nitrate counterions. EXAFS spectra are calculated directly from AIMD simulations and then compared to experimentally measured EXAFS spectra. In acetonitrile solution, both nitrate and triflate anions are shown to coordinate directly to the Eu 3+ ion forming either ten- or eight-coordinate solvent complexes where the counterions are binding as bidentate or monodentate structures, respectively. Coordination of a terpyridyl ligand to the Eu 3+ ion limits the available binding sites for the solvent and anions. In certain cases, the terpyridyl ligand excludes any solvent binding and limits the number of coordinated anions. The solution structure of the Eu-terpyridyl complex with nitrate counterions is shown to have a similar arrangement of Eu 3+ coordinating molecules as the crystal structure. Finally, this study illustrates how a combination of AIMD and EXAFS can be used to determine how ligands, solvent, and counterions coordinate with the lanthanide ions in solution.

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

Coordination Chemistry and Photoluminescence of Sm(II) Dibenzo-24-crown-8 Complexes

Three Sm(II) dibenzo-24-crown-8 (db24c8) complexes were synthesized in anhydrous, air-free conditions via the reaction of SmI 2 with db24c8 and tetrabutylammonium tetraphenylborate ([TBA][BPh 4 ]; where needed) in acetonitrile (CH 3 CN), dimethoxyethane (DME), and tetrahydrofuran (THF) to yield [Sm(db24c8)(CH 3 CN) 2 ][BPh 4 ][I]·CH 3 CN, [Sm(db24c8)(DME)]I 2 , and [Sm(db24c8)(THF) 2 ]I 2 , respectively. In each case, a 10-coordinate, staggered dodecahedral (2:6:2) environment is formed around the Sm 2+ center that is completed by either two solvent molecules (CH 3 CN or THF) or one bidentate solvent molecule (DME). Inner-sphere solvent molecules can be excluded by reacting SmI 2 with db24c8 in 1:3 THF:toluene to yield Sm(db24c8)I 2 . Here, this molecule features a distorted, eight-coordinate, hexagonal pyramidal Sm 2+ metal center, where the coordinated db24c8 molecule shows a torsion angle unexpectedly close to the 180° antiperiplanar arrangement and two uncoordinated db24c8 oxygen atoms. Solution UV–vis–NIR measurements demonstrate that Sm 2+ is a good size match for the cavity of various db24c8 conformations and that Eu 2+ and Yb 2+ exhibit competition between acetonitrile solvation and the Eu 2+ and Yb 2+ /db24c8 complexes in solution. During excitation by 546 nm light, both [Sm(db24c8)(DME)]I 2 and [Sm(db24c8)(THF) 2 ]I 2 exhibit mixed 5d → 4f and 4f → 4f emission at 20 °C and exclusively 4f → 4f at −180 °C, whereas Sm(db24c8)I 2 only shows 5d → 4f emission regardless of temperature. Photoluminescence from [Sm(db24c8)(CH 3 CN) 2 ][BPh 4 ][I]·CH 3 CN is quenched.

Cations↗
DOE OSTI2015

Predicting Stability Constants for Uranyl Complexes Using Density Functional Theory

The ability to predict the equilibrium constants for the formation of 1:1 uranyl:ligand complexes (log K 1 values) provides the essential foundation for the rational design of ligands with enhanced uranyl affinity and selectivity. We also use density functional theory (B3LYP) and the IEFPCM continuum solvation model to compute aqueous stability constants for UO 2 2+ complexes with 18 donor ligands. Theoretical calculations permit reasonably good estimates of relative binding strengths, while the absolute log K 1 values are significantly overestimated. Accurate predictions of the absolute log K 1 values (root mean square deviation from experiment < 1.0 for log K 1 values ranging from 0 to 16.8) can be obtained by fitting the experimental data for two groups of mono and divalent negative oxygen donor ligands. The utility of correlations is demonstrated for amidoxime and imide dioxime ligands, providing a useful means of screening for new ligands with strong chelate capability to uranyl.

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

Homoleptic An(IV) Dimethyl Sulfoxide Complexes of Th–Pu Enable Insight into Solution Speciation and Redox

A family of 8- and 9- coordinate homoleptic dimethyl sulfoxide An(IV) complexes is used to evaluate the solvation properties of tetravalent actinide ions in dimethyl sulfoxide (DMSO) vs aqueous media. These compounds are prepared from aqueous mixtures of DMSO and yield crystalline solids, whose solid-state electronic absorption spectra are compared to solution phase spectra, providing insight into the solution speciation. The thermodynamics of the equilibria between the 8- and 9-coordinate compounds were determined experimentally and computationally, showing a trend that correlates with ionic radius. Cyclic voltammetry data for the +IV/+III couples of Np and Pu in the DMSO electrolyte indicate that the +IV ions are significantly stabilized compared to aqueous media due to the more donating nature of DMSO compared to water. As a result, computational studies of these systems indicate that further reactivity may take place upon reduction to the +III oxidation state.

Actinides↗
DOE OSTI2026

Influence of Nd(II) and Nd(III) Ions on the Speciation and Kinetic Dynamics of Radiolytic Transients in Molten LiCl-KCl Eutectic Salt Mixtures

In the pyrochemical reprocessing of used nuclear fuel, understanding the effects of ionizing radiation on the chemical behavior of fission products in molten chloride eutectic salt is crucial. Here, this study investigates the impact of radiation on the chemical and kinetic speciation of Nd ions in molten LiCl-KCl eutectic, which are challenging to separate due to their multivalent oxidation states and comproportionation reactions. Electron pulse radiolysis techniques were used to determine transient absorption spectra, chemical kinetics, and activation parameters for the reaction of Nd ions with the primary products of molten LiCl-KCl eutectic radiolysis: the solvated electron (e S - ) and dichlorine radical anion (Cl 2 •- ). Initially, Nd(III) reacted rapidly with e S - , forming Nd(II) with a second-order rate coefficient of k = (4.54 ± 0.07) × 10 10 M -1 s -1 at 673 K. The resulting Nd(II) ions then reacted slowly with Cl 2 •- , regenerating Nd(III) with an estimated rate of k = (1.72 ± 0.04) × 10 10 M -1 s -1 at 673 K. Additionally, the comproportionation reaction of Nd(III) and Nd(s) to form Nd(II) was monitored chemically and electrochemically. This work suggests that the comproportionation reaction of Nd(III) with Nd(s) to form Nd(II) is unstable in molten LiCl-KCl eutectic.

36 - MATERIALS SCIENCE↗
DOE OSTI2020

Dynamics of Cation-Induced Conformational Changes in Nanometer-Sized Uranyl Peroxide Clusters

Conformational changes of the pyrophosphate (Pp)-functionalized uranyl peroxide nanocluster [(UO 2 ) 24 (O 2 ) 24 (P 2 O 7 ) 12 ] 48– ({U 24 Pp 12 }), dissolved as a Li/Na salt, can be induced by the titration of alkali cations into solution. The most symmetric conformer of the molecule has idealized octahedral ( O h ) molecular symmetry. One-dimensional 31 P NMR experiments provide direct evidence that both K + and Rb + ions trigger an O h -to- D 4 h conformational change within {U 24 Pp 12 }. Variable-temperature 31 P NMR experiments conducted on partially titrated {U 24 Pp 12 } systems show an effect on the rates; increased activation enthalpy and entropy for the D 4 h -to- O h transition is observed in the presence of Rb + compared to K + . Two-dimensional, exchange spectroscopy 31 P NMR revealed that magnetization transfer links chemically unique Pp bridges that are present in the D 4 h conformation and that this magnetization transfer occurs via a conformational rearrangement mechanism as the bridges interconvert between two symmetries. The interconversion is triggered by the departure and reentry of K (or Rb) cations out of and into the cavity of the cluster. This rearrangement allows Pp bridges to interconvert without the need to break bonds. Cs ions exhibit unique interactions with {U 24 Pp 12 } clusters and cause only minor changes in the solution 31 P NMR signatures, suggesting that O h symmetry is conserved. Single-crystal X-ray diffraction measurements reveal that the mixed Li/Na/Cs salt adopts D 2 h molecular symmetry, implying that while solvated, this cluster is in equilibrium with a more symmetric form. Thus, these results highlight the unusually flexible nature of the actinide-based {U 24 Pp 12 } and its sensitivity to countercations in solution.

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

Predicting Small Molecule Transfer Free Energies by Combining Molecular Dynamics Simulations and Deep Learning

Accurately predicting small molecule partitioning and hydrophobicity is critical in the drug discovery process. There are many heterogeneous chemical environments within a cell and entire human body. For example, drugs must be able to cross the hydrophobic cellular membrane to reach their intracellular targets, and hydrophobicity is an important driving force for drug–protein binding. Atomistic molecular dynamics (MD) simulations are routinely used to calculate free energies of small molecules binding to proteins, crossing lipid membranes, and solvation but are computationally expensive. Machine learning (ML) and empirical methods are also used throughout drug discovery but rely on experimental data, limiting the domain of applicability. We present atomistic MD simulations calculating 15,000 small molecule free energies of transfer from water to cyclohexane. This large data set is used to train ML models that predict the free energies of transfer. We show that a spatial graph neural network model achieves the highest accuracy, followed closely by a 3D-convolutional neural network, and shallow learning based on the chemical fingerprint is significantly less accurate. A mean absolute error of ~4 kJ/mol compared to the MD calculations was achieved for our best ML model. We also show that including data from the MD simulation improves the predictions, tests the transferability of each model to a diverse set of molecules, and show multitask learning improves the predictions. This work provides insight into the hydrophobicity of small molecules and ML cheminformatics modeling, and our data set will be useful for designing and testing future ML cheminformatics methods.

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

Computing Relative Binding Affinity of Ligands to Receptor: An Effective Hybrid Single-Dual-Topology Free-Energy Perturbation Approach in NAMD

An effective hybrid single-dual-topology protocol is designed for the calculation of relative binding affinities of small ligands to a receptor. The protocol was developed as an extension of the NAMD molecular dynamics program, which exclusively supports a dual-topology framework for relative alchemical free-energy perturbation (FEP) calculations. Here, the alchemical end states are represented as two separate molecules sharing a common substructure identified through maximum structural mapping. Within the substructure, an atom-to-atom correspondence is established, and each pair of corresponding atoms is holonomically constrained to share identical coordinates at all time throughout the simulation. The forces are projected and combined at each step for propagation. Following this formulation, a set of illustrative calculations of reliable experiment/simulation data, including relative solvation free energies of small molecules and relative binding affinities of drug compounds to proteins, are presented. To enhance sampling of the dual-topology region, the FEP calculations were carried out within a replica-exchange MD scheme supported by the multiple-copy algorithm module of NAMD, with periodically attempted swapping of the thermodynamic coupling parameter lambda between neighboring states. The findings are consistent with experiments and benchmarks reported in the literature, lending support to the validity of the current protocol. In summary, this hybrid single-dual-topology approach combines the conceptual simplicity of the dual-topology paradigm with the advantageous sampling efficiency of the single-topology method, making it an ideal strategy for high-throughput in silico drug design.

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

Predicting Excitation Energies of Twisted Intramolecular Charge-Transfer States with the Time-Dependent Density Functional Theory: Comparison with Experimental Measurements in the Gas Phase and Solvents Ranging from Hexanes to Acetonitrile

Electronically excited states characterized by intramolecular charge transfer play an essential role in many biological processes and optical devices. The ability to make quantitative ab initio predictions of the relative energetics involved is a challenging yet desirable goal, especially for large molecules in solution. In this work, we present a data set of 61 experimental measurements of absorption and emission processes, both in the gas phase and in solvents representing a broad range of polarities, which involve intramolecular charge transfer mediated by a nonzero, "twisted" dihedral angle between one or more donor and acceptor subunits. Among a variety of density functionals investigated within the framework of linear-response theory, the "optimally tuned" LRC-ωPBE functional, which utilizes a system-specific yet nonempirical procedure to specify the range-separation parameter, emerges as the preferred choice. For the entire set of excitation energies, involving changes in dipole moment ranging from 4 to >20 Debye, the mean signed and absolute errors are 0.02 and 0.18 eV, respectively (compared, e.g., to -0.30 and 0.30 for PBE0, 0.44 and 0.47 for LRC-ωPBEh, 0.83 and 0.83 for ωB97X-V). We analyze the performance of polarizable continuum solvation models available in Q-Chem that partition the solvent response into fast and slow time scales, and clear trends emerge when measurements corresponding to the four small 4-(dimethylamino)benzonitrile (DMABN)-like molecules and a charged species are excluded. We make the case that the large errors found only for small molecules in the gas phase and weak solvents cannot be expected to improve via the optimal tuning procedure, which enforces a condition that is exact only in the well-separated donor-acceptor limit, and present empirical evidence implicating the outsized importance for small donor-acceptor systems of relaxation effects that cannot be accounted for by the linear-response time-dependent density functional theory within the adiabatic approximation. Finally, we demonstrate the utility of the optimally tuned density functional approach by targeting the charge-transfer states of a large biomimetic model system for light-harvesting structures in Photosystem II.

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

Determining Partial Atomic Charges for Liquid Water: Assessing Electronic Structure and Charge Models

Partial atomic charges provide an intuitive and efficient way to describe the charge distribution and the resulting intermolecular electrostatic interactions in liquid water. Many charge models exist and it is unclear which model provides the best assignment of partial atomic charges in response to the local molecular environment. In this work, we systematically scrutinize various electronic structure methods and charge models (Mulliken, natural population analysis, CHelpG, RESP, Hirshfeld, Iterative Hirshfeld, and Bader) by evaluating their performance in predicting the dipole moments of isolated water, water clusters, and liquid water as well as charge transfer in the water dimer and liquid water. Although none of the seven charge models is capable of fully capturing the dipole moment increase from isolated water (1.85 D) to liquid water (about 2.9 D), the Iterative Hirshfeld method performs best for liquid water, reproducing its experimental average molecular dipole moment, yielding a reasonable amount of intermolecular charge transfer, and showing modest sensitivity to the local water environment. The performance of the charge model is dependent on the choice of the density functional and the quantum treatment of the environment. The computed molecular dipole moment of water generally increases with the percentage of the exact Hartree–Fock exchange in the functional, whereas the amount of charge transfer between molecules decreases. For liquid water, including two full solvation shells of surrounding water molecules (within about 5.5 Å of the central water) in the quantum chemical calculation converges the charges of the central water molecule. Furthermore, our final pragmatic quantum chemical charge-assigning protocol for liquid water is the Iterative Hirshfeld method with M06-HF/aug-cc-pVDZ and a quantum region cutoff radius of 5.5 Å.

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

Forces from Stochastic Density Functional Theory under Nonorthogonal Atom-Centered Basis Sets

We develop a formalism for calculating forces on the nuclei within the linear-scaling stochastic density functional theory (sDFT) in a nonorthogonal atom-centered basis set representation (Fabian et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2019, 9, e1412, 10.1002/wcms.1412) and apply it to the Tryptophan Zipper 2 (Trp-zip2) peptide solvated in water. We use an embedded-fragment approach to reduce the statistical errors (fluctuation and systematic bias), where the entire peptide is the main fragment and the remaining 425 water molecules are grouped into small fragments. We analyze the magnitude of the statistical errors in the forces and find that the systematic bias is of the order of 0.065 eV/Å (~1.2 × 10 -3 E h / a 0 ) when 120 stochastic orbitals are used, independently of system size. This magnitude of bias is sufficiently small to ensure that the bond lengths estimated by stochastic DFT (within a Langevin molecular dynamics simulation) will deviate by less than 1% from those predicted by a deterministic calculation.

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

Accurate Free Energies for Complex Condensed-Phase Reactions Using an Artificial Neural Network Corrected DFTB/MM Methodology

Semiempirical methods like density functional tight-binding (DFTB) allow extensive phase space sampling, making it possible to generate free energy surfaces of complex reactions in condensed-phase environments. Such a high efficiency often comes at the cost of reduced accuracy, which may be improved by developing a specific reaction parametrization (SRP) for the particular molecular system. Thiol–disulfide exchange is a nucleophilic substitution reaction that occurs in a large class of proteins. Its proper description requires a high-level ab initio method, while DFT-GAA and hybrid functionals were shown to be inadequate, and so is DFTB due to its DFT-GGA descent. We develop an SRP for thiol–disulfide exchange based on an artificial neural network (ANN) implementation in the DFTB+ software and compare its performance to that of a standard SRP approach applied to DFTB. Furthermore, as an application, we use both new DFTB-SRP as components of a QM/MM scheme to investigate thiol–disulfide exchange in two molecular complexes: a solvated model system and a blood protein. Demonstrating the strengths of the methodology, highly accurate free energy surfaces are generated at a low cost, as the augmentation of DFTB with an ANN only adds a small computational overhead.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗