Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum mechanical calculations”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Experimental Protocols for Studying Organic Non-aqueous Redox Flow Batteries

We report that Redox flow batteries (RFBs) are promising devices for grid-scale energy storage due to the decoupling of power and energy, which can be independently scaled by the electrode area and storage tank size, respectively. To date, only aqueous RFBs, such as the vanadium RFB, have been implemented commercially. Nevertheless, the limited energy densities and high-cost materials may preclude their wider market penetration. Organic non-aqueous redox-flow batteries (O-NRFBs), which utilize redox-active organic molecules (ROMs), have been offered as an attractive alternative. Many possible advantages include the use of earth-abundant elements (C, H, N, O, S, F), wherein the ROMs can be prepared from low-cost and sustainable materials. Additionally, a large variety of electroactive moieties are accessible as building blocks, providing a synthetic platform to tune the properties of ROMs through rational design. As a consequence, the development of novel ROMs has attracted researchers with diverse backgrounds, prompting remarkable innovations in the past decade. However, consistency in experimental protocols (e.g., electrochemical methods, cycling stability, experimental conditions) has not coincided with this uptick in research, leading to sometimes ambiguous and incomparable results. This is further convoluted by the complexity innate to battery development, such as cell design, detection, and characterization of reactions and active components. The O-NRFB application imposes stringent requirements on the physical or/and chemical processes involved in electrochemical cycling. Yet such considerations are often overlooked. Further, while quantum mechanical calculations provide a convenient tool for evaluating and selecting ROM candidates, this simulation is not always performed. Thus, in this Energy Focus, we detail a means to standardize experimental protocols for studies of O-NRFBs and suggest practices to facilitate fundamental understanding and development of ROMs.

25 ENERGY STORAGE↗
DOE OSTI2025

Surface-Controlled TiO 2 Nanocrystals with Catalytically Active Single-Site Co Incorporation for the Oxygen Evolution Reaction

The design of advanced electrocatalysts is often hindered by uncertainties in identifying and controlling the active surfaces and catalytic centers within heterogeneous materials. Here we present the synthesis of single-site Co catalysts, substitutionally doped into surface-controlled TiO 2 anatase nanocrystals, aimed at enhancing the oxygen evolution reaction (OER). Grand canonical quantum mechanics calculations reveal that the kinetics of the OER, following an adsorbate evolution mechanism, is markedly influenced by the coordination environment of Co. The simulations suggest significantly higher turnover frequencies when Co is doped into the (001) surface of TiO 2 compared to the (101) surface. Consistent with the computational findings, experimental results show that Co-doped TiO 2 (Co-TiO 2 ) nanoplates with selectively exposed {001} surfaces exhibit enhanced current densities and turnover frequencies compared to Co-TiO 2 nanobipyramids with {101} surfaces. This study highlights the synergy between theoretical calculations and precision synthesis in the development of more effective catalysts.

25 ENERGY STORAGE↗
DOE OSTI2022

Realistic magnetic thermodynamics by local quantization of a semiclassical Heisenberg model

Abstract Classical Monte Carlo simulation of the Heisenberg model poorly describes many thermodynamic phenomena due to its neglect of the quantum nature of spins. Alternatively, we discuss how to semiclassically approach the quantum problem and demonstrate a simple method for introducing a locally approximate form of spin quantization. While the procedure underestimates magnetic short-range order, our results suggest a simple correction for recovering realistic spin–spin correlations above the critical temperature. Moreover, ensemble fluctuations are found to provide reasonably accurate thermodynamics, largely reproducing quantum mechanically calculated heat capacities and experimental magnetometry for ferromagnetic Fe and antiferromagnetic RbMnF 3 . Extensions of the method are proposed to address remaining inaccuracies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Oxygen evolution reaction over catalytic single-site Co in a well-defined brookite TiO 2 nanorod surface

Efficient electrocatalysts for the oxygen evolution reaction (OER) are paramount to the development of electrochemical devices for clean energy and fuel conversion. However, the structural complexity of heterogeneous electrocatalysts makes it a great challenge to elucidate the surface catalytic sites and OER mechanisms. Here, we report that catalytic single-site Co in a well-defined brookite TiO 2 nanorod (210) surface (Co-TiO 2 ) presents turnover frequencies that are among the highest for Co-based heterogeneous catalysts reported to date, reaching 6.6 ± 1.2 and 181.4 ± 28 s -1 at 300 and 400 mV overpotentials, respectively. Based on grand canonical quantum mechanics calculations and the single-site Co atomic structure validated by in situ and ex situ spectroscopic probes, we have established a full description of the catalytic reaction kinetics for Co-TiO 2 as a function of applied potential, revealing an adsorbate evolution mechanism for the OER. The computationally predicted Tafel slope and turnover frequencies exhibit exceedingly good agreement with experiment.

25 ENERGY STORAGE↗
DOE OSTI2025

Using phage display for rational engineering of a higher-affinity humanized 3’ phosphohistidine-specific antibody

Abstract Histidine phosphorylation is a non-canonical post-translational modification (PTM), with 1-phosphohistidine (1-pHis) and 3-phosphohistidine (3-pHis) isoforms, that is understudied due to a lack of robust reagents, including high-affinity pHis-specific antibodies. Engineering pHis antibodies is challenging due to the labile nature of its phosphoramidate (P-N) bond. We developed a strategy for in vitro engineering of antibodies for the detection of native 3-pHis targets, in which the rabbit SC44-8 anti-3-pTza mAb is humanized into a scaffold (hSC44) that is suitable for phage display. Six unique Fab phage-displayed hSC44 scaffold libraries were screened for antibodies that bound 3-pHis with higher affinity and had specificity for 3-pHis versus 3-pTza. hSC44.20N32F L , the best engineered antibody, has ~10-fold higher affinity for 3-pHis than parental hSC44. Eleven new Fab structures, including the first antibody-pHis peptide structures, together with structural and quantum mechanical calculations, provided molecular insights into 3-pHis and 3-pTza discrimination by hSC44.20N32F L and the increased affinity obtained through engineering. We demonstrated the utility of these high-affinity 3-pHis-specific antibodies for the recognition of pHis proteins in mammalian cells by immunoblotting and immunofluorescence staining. Our work describes a general method for engineering labile PTM-specific antibodies and provides novel antibodies for investigating the role of 3-pHis in cell biology.

Martyn, Gregory D.↗
DOE OSTI2022

Molecular-level understanding of metal ion retention in clay-rich materials

Clay minerals retain or adsorb metal ions in the Earth’s critical zone. Rocks, sediments and soils rich in clay minerals can concentrate rare earth elements (REEs) in ion adsorption-type deposits (IADs) and are similarly effective at metallic contaminant remediation. However, the molecular-scale chemical and physical mechanisms of metal ion retention remain only partly understood. Here, we describe the nature, location and energy requirements of metal retention at clay mineral surfaces. Retention originates mainly from electrostatic interactions during cation exchange at low pH and chemical bonding in surface complexation and precipitation at neutral and high pH. Surface complexation can induce surface redox reactions and precipitation mechanisms including neoformation of clay mineral layered structures. In IADs, outer-sphere adsorption is the major retention mechanism of REE ions. By contrast, the use of clay minerals in pollution control relies on various mechanisms that can coexist, including cation exchange, surface complexation and nucleation growth. To more effectively leverage clay mineral–metal interactions in resource recovery and contaminant remediation, complex mechanisms such as surface precipitation and redox reactions must be better understood; for instance, by utilizing advances in quantum mechanical calculations, close combination between synchrotron and simulation techniques, and upscaling of molecular-level information in macroscopic thermokinetic predictive models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Thermodynamic consequences of Tyr to Trp mutations in the cation–π-mediated binding of trimethyllysine by the HP1 chromodomain

Evolution has converged on cation–p interactions for recognition of quaternary alkyl ammonium groups such as trimethyllysine (Kme3). While computational modelling indicates that Trp provides the strongest cation–p interaction of the native aromatic amino acids, there is limited corroborative data from measurements within proteins. Herein we investigate a Tyr to Trp mutation in the binding pocket of the HP1 chromodomain, a reader protein that recognizes Kme3. Binding studies demonstrate that the Trpmediated cation–p interaction is about -5 kcal mol -1 stronger, and the Y24W crystal structure shows that the mutation is not perturbing. Quantum mechanical calculations indicate that greater enthalpic binding is predominantly due to increased cation–π interactions. NMR studies indicate that differences in the unbound state of the Y24W mutation lead to enthalpy–entropy compensation. These results provide direct experimental quantification of Trp versus Tyr in a cation–π interaction and afford insight into the conservation of aromatic cage residues in Kme3 reader domains.

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

Development and application of ReaxFF methodology for understanding the chemical dynamics of metal carbonates in aqueous solutions

A new ReaxFF reactive force field has been developed for metal carbonate systems including Na + , Ca 2+ , and Mg 2+ cations and the CO 3 2– anion. This force field is fully transferable with previous ReaxFF water and water/electrolyte descriptions. The Me–O–C (Me = metal) three-body valence angle parameters and Me–C non-reactive parameters of the force field have been optimized against quantum mechanical calculations including equations of state, heats of formation, heats of reaction, angle distortions and vibrational frequencies. The new metal carbonate force field has been validated using molecular dynamics simulations to study the solvation and reactivity of metal and carbonate ions in water at 300 K and 700 K. The coordination radius and self-diffusion coefficient show good consistency with existing experimental and simulation results. The angular distribution analysis explains the structural preference of carbonate ions to form carbonates and bicarbonates, where Na + predominantly forms carbonates due to weaker angular strain, while Ca 2+ and Mg 2+ prefer to form bicarbonate monodentate in nature. Residence time distribution analyses on different systems reveal the role of ions in accelerating and decelerating the dynamics of water and carbonate ions under different thermodynamic conditions. The formation and dissolution of bicarbonates and carbonates in solution were explored on the basis of the protonation capability in different systems. Here, the nucleation phenomenon of metal carbonates at ambient and supercritical conditions is explained from the perspective of cluster formation over time: Ca 2+ ions can form prenucleation clusters at ambient temperature but show saturation with increasing temperature, whereas Na + and Mg 2+ ions show a rapid increase in cluster size and amount upon increasing time and temperature.

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

Impact of lanthanide ion complexation and temperature on the chemical reactivity of N , N , N ', N '-tetraoctyl diglycolamide (TODGA) with the dodecane radical cation

The impact of trivalent lanthanide ion complexation and temperature on the chemical reactivity of N,N,N',N'-tetraoctyl diglycolamide (TODGA) with the n-dodecane radical cation (RH˙+) has been measured by electron pulse radiolysis and evaluated by quantum mechanical calculations. Additionally, Arrhenius parameters were determined for the reaction of the non-complexed TODGA ligand with the RH˙ + from 10–40 °C, giving the activation energy (E a = 17.43 ± 1.64 kJ mol –1 ) and pre-exponential factor (A = (2.36 ± 0.05) × 10 13 M –1 s –1 ). The complexation of Nd(III), Gd(III), and Yb(III) ions by TODGA yielded [LnIII(TODGA)3(NO3)3] complexes that exhibited significantly increased reactivity (up to 9.3× faster) with the RH˙ + , relative to the non-complexed ligand: k([Ln III (TODGA) 3 (NO 3 ) 3 ] + RH˙ + ) = (8.99 ± 0.93) × 10 10 , (2.88 ± 0.40) × 10 10 , and (1.53 ± 0.34) × 10 10 M –1 s –1 , for Nd(III), Gd(III), and Yb(III) ions, respectively. The rate coefficient enhancement measured for these complexes exhibited a dependence on atomic number, decreasing as the lanthanide series was traversed. Preliminary reaction free energy calculations—based on a model [Ln III (TOGDA)] 3+ complex system—indicate that both electron/hole and proton transfer reactions are energetically unfavorable for complexed TODGA. Furthermore, complementary average local ionization energy calculations showed that the most reactive region of model N,N,N',N'-tetraethyl diglycolamide (TEDGA) complexes, [Ln III (TEGDA) 3 (NO 3 ) 3 ], toward electrophilic attack is for the coordinated nitrate (NO 3 – ) counter anions. Furthermore, it is possible that radical reactions with the complexed NO 3 – counter anions dominate the differences in rates seen for the [Ln III (TODGA) 3 (NO 3 ) 3 ] complexes, and are likely responsible for the reported radioprotection in the presence of TODGA complexes.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Porosity modeling in a TiNbTaZrMo high-entropy alloy for biomedical applications

High-entropy alloys (HEAs) have attracted great attention for many biomedical applications. However, the nature of interatomic interactions in this class of complex multicomponent alloys is not fully understood. We report, for the first time, the results of theoretical modeling for porosity in a large biocompatible HEA TiNbTaZrMo using an atomistic supercell of 1024 atoms that provides new insights and understanding. Our results demonstrated the deficiency of using the valence electron count, quantification of large lattice distortion, validation of mechanical properties with available experimental data to reduce Young's modulus. We utilized the novel concepts of the total bond order density (TBOD) and partial bond order density (PBOD) via ab initio quantum mechanical calculations as an effective theoretical means to chart a road map for the rational design of complex multicomponent HEAs for biomedical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effect of f -element complexation on the radiolysis of 2-ethylhexylphosphonic acid mono-2-ethylhexyl ester (HEH[EHP])

A systematic study of the impact on the chemical reactivity of the oxidising n-dodecane radical cation (RH˙ + ) with f-element complexed 2-ethylhexylphosphonic acid mono-2-ethylhexyl ester (HEH[EHP]) has been undertaken utilizing time-resolved electron pulse radiolysis/transient absorption spectroscopy and high-level quantum mechanical calculations. Lanthanide ion complexed species, [Ln((HEH[EHP]) 2 ) 3 ], exhibited vastly increased reactivity (over 10× faster) in comparison to the non-complexed ligand in n-dodecane solvent, whose rate coefficient was k = (4.66 ± 0.22) × 10 9 M -1 s -1 . Similar reactivity enhancement was also observed for the corresponding americium ion complex, k = (5.58 ± 0.30) × 10 10 M -1 s -1 . The vastly increased reactivity of these f-element complexes was not due to simple increased diffusion-control of these reactions; rather, enhanced hole transfer mechanisms for the complexes were calculated to become energetically more favourable. Interestingly, the observed reactivity trend with lanthanide ion size was not linear; instead, the rate coefficients showed an initial increase (Lu to Yb) followed by a decrease (Tm to Ho), followed by another increase (Dy to La). This behaviour was excellently predicted by the calculated reaction volumes of these complexes. In conclusion, complementary cobalt-60 gamma irradiations for select lanthanide complexes demonstrated that the measured kinetic differences translated to increased ligand degradation at steady-state timescales, affording ~38% increase in ligand loss of a 1:1 [La((HEH[EHP]) 2 ) 3 ]: HEH[EHP] ratio system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Atomistic mechanisms for catalytic transformations of NO to NH3, N2O, and N2 by Pd

The industrial pollutant NO is a potential threat to the environment and to human health. Thus, selective catalytic reduction of NO into harmless N2, NH3, and/or N2O gas is of great interest. Among many catalysts, metal Pd has been demonstrated to be most efficient for selectivity of reducing NO to N2. However, the reduction mechanism of NO on Pd, especially the route of N−N bond formation, remains unclear, impeding the development of new, improved catalysts. We report here the elementary reaction steps in the reaction pathway of reducing NO to NH3, N2O, and N2, based on density functional theory (DFT)-based quantum mechanics calculations. We show that the formation of N2O proceeds through an Eley-Rideal (E−R) reaction pathway that couples one adsorbed NO* with one non−adsorbed NO from the solvent or gas phase. This reaction requires high NO* surface coverage, leading first to the formation of the trans-(NO)2* intermediate with a low N−N coupling barrier (0.58 eV). Notably, trans-(NO)2* will continue to react with NO in the solvent to form N2O, that has not been reported. With the consumption of NO and the formation of N2O* in the solvent, the Langmuir-Hinshelwood (L-H) mechanism will dominate at this time, and N2O* will be reduced by hydrogenation at a low chemical barrier (0.42 eV) to form N2. In contrast, NH3 is completely formed by the L-H reaction, which has a higher chemical barrier (0.87 eV). Our predicted E-R reaction has not previously been reported, but it explains some existing experimental observations. In addition, we examine how catalyst activity might be improved by doping a single metal atom (M) at the NO* adsorption site to form M/Pd and show its influence on the barrier for forming the N−N bond to provide control over the product distribution.

Physics↗
DOE OSTI2020

Predicting molecular dipole moments by combining atomic partial charges and atomic dipoles

The molecular dipole moment ( μ ) is a central quantity in chemistry. It is essential in predicting infrared and sum-frequency generation spectra as well as induction and long-range electrostatic interactions. Furthermore, it can be extracted directly—via the ground state electron density—from high-level quantum mechanical calculations, making it an ideal target for machine learning (ML). Here, we choose to represent this quantity with a physically inspired ML model that captures two distinct physical effects: local atomic polarization is captured within the symmetry-adapted Gaussian process regression framework which assigns a (vector) dipole moment to each atom, while the movement of charge across the entire molecule is captured by assigning a partial (scalar) charge to each atom. The resulting “MuML” models are fitted together to reproduce molecular μ computed using high-level coupled-cluster theory and density functional theory (DFT) on the QM7b dataset, achieving more accurate results due to the physics-based combination of these complementary terms. The combined model shows excellent transferability when applied to a showcase dataset of larger and more complex molecules, approaching the accuracy of DFT at a small fraction of the computational cost. We also demonstrate that the uncertainty in the predictions can be estimated reliably using a calibrated committee model. The ultimate performance of the models—and the optimal weighting of their combination—depends, however, on the details of the system at hand, with the scalar model being clearly superior when describing large molecules whose dipole is almost entirely generated by charge separation. These observations point to the importance of simultaneously accounting for the local and non-local effects that contribute to μ ; furthermore, they define a challenging task to benchmark future models, particularly those aimed at the description of condensed phases.

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

Machine learning substitutional defect formation energies in ABO 3 perovskites

Perovskite oxides are a promising material platform for use in a wide range of technological applications including electronics, sensors, fuel cells, and catalysis. This is owing to the extraordinary tunability of their physical and chemical properties via defect engineering. The feasibility and the stability of a defect, such as a substitutional dopant, in the host lattice is usually obtained via experiments and/or through detailed quantum mechanical calculations. Both of these conventional routes are expensive and time consuming. An alternative is a data-driven machine learning (ML)-based approach. In this work, we have applied ML techniques to identify the factors that influence defect formation energy, which is an important measure of the stability of the defects, in perovskite oxides. Using 13 elemental properties as features and random forest regression, we demonstrate a systematic approach to down-selecting from the larger set of features to those that are important, establishing a framework for accurate predictions of the defect formation energy. We quantitatively show that the most important factors that control the dopant stability are the dopant ionic size, heat of formation, effective tolerance factor, and oxidation state. Our work reveals previously unknown correlations, chemical trends, and the interplay between stability and underlying chemistries. Hence, these results showcase the efficacy of ML tools in identifying and quantifying different feature-dependencies and provide a promising route toward dopant selection in the perovskites. Herein, we have developed a framework that itself is general and can be potentially applied to other material classes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Subspace recursive Fermi-operator expansion strategies for large-scale DFT eigenvalue problems on HPC architectures

Quantum mechanical calculations for material modeling using Kohn–Sham density functional theory (DFT) involve the solution of a nonlinear eigenvalue problem for N smallest eigenvector-eigenvalue pairs, with N proportional to the number of electrons in the material system. Here, these calculations are computationally demanding and have asymptotic cubic scaling complexity with the number of electrons. Large-scale matrix eigenvalue problems arising from the discretization of the Kohn–Sham DFT equations employing a systematically convergent basis traditionally rely on iterative orthogonal projection methods, which are shown to be computationally efficient and scalable on massively parallel computing architectures. However, as the size of the material system increases, these methods are known to incur dominant computational costs through the Rayleigh–Ritz projection step of the discretized Kohn–Sham Hamiltonian matrix and the subsequent subspace diagonalization of the projected matrix. This work explores the potential of polynomial expansion approaches based on recursive Fermi-operator expansion as an alternative to the subspace diagonalization of the projected Hamiltonian matrix to reduce the computational cost. Subsequently, we perform a detailed comparison of various recursive polynomial expansion approaches to the traditional approach of explicit diagonalization on both multi-node central processing unit and graphics processing unit architectures and assess their relative performance in terms of accuracy, computational efficiency, scaling behavior, and energy efficiency.

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

The structure of the high-affinity nickel-binding site in the Ni,Zn-HypA•UreE2 complex

Abstract The maturation pathway for the nickel-dependent enzyme urease utilizes the protein UreE as a metallochaperone to supply Ni(II) ions. In Helicobacter pylori urease maturation also requires HypA and HypB, accessory proteins that are commonly associated with hydrogenase maturation. Herein we report on the characterization of a protein complex formed between HypA and the UreE2 dimer. Nuclear magnetic resonance (NMR) coupled with molecular modelling show that the protein complex apo, Zn-HypA•UreE2, forms between the rigorously conserved Met-His-Glu (MHE motif) Ni-binding N-terminal sequence of HypA and the two conserved His102A and His102B located at the dimer interface of UreE2. This complex forms in the absence of Ni(II) and is supported by extensive protein contacts that include the use of the C-terminal sequences of UreE2 to form additional strands of β-sheet with the Ni-binding domain of HypA. The Ni-binding properties of apo, Zn-HypA•UreE2 and the component proteins were investigated by isothermal titration calorimetry using a global fitting strategy that included all of the relevant equilibria, and show that the Ni,Zn-HypA•UreE2 complex contains a single Ni(II)-binding site with a sub-nanomolar KD. The structural features of this novel Ni(II) site were elucidated using proteins produced with specifically deuterated amino acids, protein point mutations, and the analyses of X-ray absorption spectroscopy, hyperfine shifted NMR features, as well as molecular modeling coupled with quantum-mechanical calculations. The results show that the complex contains a six-coordinate, high-spin Ni(II) site with ligands provided by both component proteins.

Zambelli, Barbara (ORCID:0000000238760051)↗