SEARCH · Search NASA
Results for “ab initio modeling”
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.
Tuning Point Defects by Elastic Strain Modulates Nanoparticle Exsolution on Perovskite Oxides
Exsolution generates stable and catalytically active metal nanoparticles via phase precipitation out of a host oxide. An ability to control the size and dispersion of the exsolution particles is desirable for design of nanostructured (electro)catalysts. In this work, we demonstrate that tuning point defects by lattice strain affects both the thermodynamics and the kinetics of iron (Fe 0 ) exsolution on La 0.6 Sr 0.4 FeO 3 (LSF) thin film model. By combining in situ surface characterization and ab initio defect modeling, we show oxygen vacancy and Schottky defects to be the primary point defects formed upon Fe0 exsolution. Lattice strain tunes the formation energy, and thus the abundance of these defects, and alters the amount and size of the resulting exsolution particles. In addition, we find that the density of exsolved nanoparticles matches the concentration of oxygen vacancy pairs, thus pointing to the surface oxygen vacancy pairs as preferential nucleation sites for exsolution. The tensile-strained LSF with a facile formation of these critical point defects results in a higher Fe 0 metal concentration, a larger density of nanoparticles, and a reduced particle size at its surfaces. These results provide important mechanistic insights and highlight the role of point-defect engineering in designing nanostructured catalysts in energy and fuel conversion technologies.
Signatures of a liquid–liquid transition in an ab initio deep neural network model for water
Significance Water is central across much of the physical and biological sciences and exhibits physical properties that are qualitatively distinct from those of most other liquids. Understanding the microscopic basis of water’s peculiar properties remains an active area of research. One intriguing hypothesis is that liquid water can separate into metastable high- and low-density liquid phases at low temperatures and high pressures, and the existence of this liquid–liquid transition could explain many of water’s anomalous properties. We used state-of-the-art approaches in computational quantum chemistry, statistical mechanics, and machine learning and obtained evidence consistent with a liquid–liquid transition, supporting the argument for the existence of this phenomenon in real water.
From bulk to surface: Structure and dynamics of amorphous alumina from deep potential molecular dynamics
Understanding the atomic-scale structure and dynamics of amorphous oxide surfaces is essential for interpreting their chemical reactivity, mechanical stability, and interfacial behavior, yet direct experimental characterization remains challenging. We employ Deep Potential (DP) molecular dynamics to generate large-scale, ab initio -quality models of amorphous Al 2 O 3 bulk glasses and melt-quenched free surfaces, enabling a quantitative analysis of both structure and relaxation dynamics with statistical confidence inaccessible to direct ab initio simulation. The trained DP model reproduces experimental liquid and glass structure, captures the cooling-rate dependence of the bulk glass transition, and corrects systematic biases in the polyhedral populations predicted by widely used classical force fields. At the free surface, mass density recovers to bulk values over ~10 Å, while local coordination requires a slightly wider subsurface region to fully converge. The outermost layer is oxygen-enriched, exhibits altered polyhedral connectivity with contracted Al–O bonds, and hosts a broad population of under-coordinated motifs (notably AlO 3 and OAl 2 ) whose abundances are governed by glass stability. These under-coordinated surface motifs exhibit distinct vibrational signatures and occur as locally paired Lewis acid and Brønsted base sites consistent with bond-valence compensation, yet remain spatially dispersed rather than aggregating into extended clusters. Despite this pronounced structural heterogeneity, surface relaxation and the glass-transition temperature remain comparable to their bulk counterparts, suggesting that the disordered surface is kinetically stable once formed. Together, these results establish a molecular-level picture of amorphous alumina surfaces and demonstrate the capability of machine-learned potentials to resolve structure–property relationships in disordered oxide interfaces.
The Gamow shell model with realistic interactions: a theoretical framework for ab initio nuclear structure at drip-lines
Abstract Ab initio approaches are among the most advanced models to solve the nuclear many-body problem. In particular, the no-core–shell model and many-body perturbation theory have been recently extended to the Gamow shell model framework, where the harmonic oscillator basis is replaced by a basis bearing bound, resonance and scattering states, i.e. the Berggren basis. As continuum coupling is included at basis level and as configuration mixing takes care of inter-nucleon correlations, halo and resonance nuclei can be properly described with the Gamow shell model. The development of the no-core Gamow shell model and the introduction of the Q ¯ ˆ -box method in the Gamow shell model, as well as their first ab initio applications, will be reviewed in this paper. Peculiarities compared to models using harmonic oscillator bases will be shortly described. The current power and limitations of ab initio Gamow shell model will also be discussed, as well as its potential for future applications.
Ab initio calculations of material properties for modeling debris (FY2021 Progress Report)
This report summarizes the work done during the first 6 months of the project “Ab initio calculations of material properties for modeling debris” that is part of the Rapid Response Research (R3) Venture. The main focus for this fiscal year (FY) has been the evaluation of classical interatomic potentials for the description of the condensation of uranium oxides in the regimes relevant to fireball conditions. We also delineate how these functions will be used within the molecular dynamics and Monte Carlo frameworks to compute thermodynamic properties of these materials. These descriptions are to be utilized within continuum modeling formulations to generate a more accurate characterization of the homogeneous nucleation process.
Magnetoelastic coupling anisotropy in the Kitaev material α-RuCl 3
Here we show that the Kitaev material α-RuCl 3 is among the most prominent candidates to host a quantum spin-liquid state endowed with fractionalized excitations. Recent experimental and theoretical investigations have separately revealed the importance of both the magnetoelastic coupling and the magnetic anisotropy in dependence of the applied magnetic-field direction. In this combined theoretical and experimental research, we investigate the anisotropic magnetic and magnetoelastic properties for magnetic fields applied along the main crystallographic axes as well as for fields canted out of the honeycomb plane. We found that the magnetostriction anisotropy is unusually large compared to the anisotropy of the magnetization, which is related to the strong magnetoelastic $\tilde{Γ}'$-type coupling in our ab-initio derived model. We observed large, nonsymmetric magnetic anisotropy for magnetic fields canted out of the honeycomb ab plane in opposite directions, namely, toward the +c* or -c* axes, respectively. The observed directional anisotropy is explained by considering the relative orientation of the magnetic field with respect to the coaligned RuCl 6 octahedra. Magnetostriction measurements in canted fields support this nonsymmetric magnetic anisotropy; however, these experiments are affected by magnetic torque effects. Comparison of theoretical predictions with experimental findings allow us to recognize the significant contribution of torque effects in experimental setups where α-RuCl 3 is placed in canted magnetic fields.
Learning intermolecular forces at liquid–vapor interfaces
By adopting a perspective informed by contemporary liquid-state theory, we consider how to train an artificial neural network potential to describe inhomogeneous, disordered systems. Here, we find that neural network potentials based on local representations of atomic environments are capable of describing some properties of liquid-vapor interfaces but typically fail for properties that depend on unbalanced long-ranged interactions that build up in the presence of broken translation symmetry. These same interactions cancel in the translationally invariant bulk, allowing local neural network potentials to describe bulk properties correctly. By incorporating explicit models of the slowly varying long-ranged interactions and training neural networks only on the short-ranged components, we can arrive at potentials that robustly recover interfacial properties. We find that local neural network models can sometimes approximate a local molecular field potential to correct for the truncated interactions, but this behavior is variable and hard to learn. Generally, we find that models with explicit electrostatics are easier to train and have higher accuracy. We demonstrate this perspective in a simple model of an asymmetric dipolar fluid, where the exact long-ranged interaction is known, and in an ab initio water model, where it is approximated.
Correlated conformational dynamics of the human GluN1-GluN2A type N-methyl-D-aspartate (NMDA) receptor
N-Methyl-D-aspartate receptors (NMDARs) are glutamate-gated ion channels found in the nerve cell membranes. As a result of overexcitation of NMDARs, neuronal death occurs and may lead to diseases such as epilepsy, stroke, Alzheimer’s disease, and Parkinson’s disease. In this study, human GluN1- GluN2A type NMDAR structure is modeled based on the X-ray structure of the Xenopus laevis template and missing loops are added by ab-initio loop modeling. The final structure is chosen according to two different model assessment scores. To be able to observe the structural changes upon ligand binding, glycine and glutamate molecules are docked into the corresponding binding sites of the receptor. Subsequently, molecular dynamics simulations of 1.3 μs are performed for both apo and ligand-bound structures. Structural parameters, which have been considered to show functionally important changes in previous NMDAR studies, are monitored as conformational rulers to understand the dynamics of the conformational changes. Moreover, principal component analysis (PCA) is performed for the equilibrated part of the simulations. Additionally, from these analyses, the differences in between apo and ligand-bound simulations can be summarized as the following: The girdle right at the beginning of the pore loop, which connects M2 and M3 helices of the ion channel, partially opens. Ligands act like an adhesive for the ligand-binding domain (LBD) by keeping the bi-lobed structure together and consequently this is reflected to the overall dynamics of the protein as an increased correlation of the LBD with especially the amino-terminal domain (ATD) of the protein.
Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods
We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.
Moiré Kanamori-Hubbard model for transition metal dichalcogenide homobilayers
Ab initio and continuum model studies predicted that the Γ valley transition metal dichalcogenide (TMD) homobilayers could simulate the conventional multiorbital Hubbard model on the moiré honeycomb lattice. Here, in this work, we perform the Wannierization starting from the continuum model and show that a more general moiré Kanamori-Hubbard model emerges, beyond the extensively studied standard multiorbital Hubbard model, which can be used to investigate the many-body physics in the Γ valley TMD homobilayers. Using the unrestricted Hartree-Fock and Lanczos techniques, we study these half-filled multiorbital moiré bands. By constructing the phase diagrams we predict the presence of an antiferromagnetic state and in addition we found unexpected and dominant states, such as a S = 1 ferromagnetic insulator and a charge density wave state. Our theoretical predictions made using this model can be tested in future experiments on the Γ valley TMD homobilayers.
Stochastic shock waves as a candidate mechanism for the formation of the He I 10830-A line in cool giant stars
The question of whether strong shocks produced in time-dependent stochastic wave models can explain the formation of the He I 10830-A line in cool giant stars is examined. The research is based on the ab initio chromosphere model for Arcturus by Cuntz (1987), showing that a stochastic distribution of wave periods leads to overtaking and merging of shocks, which occasionally produce very strong shocks with temperatures larger than 40,000 K in the postshock regions. These temperatures can easily produce a significant population in the 2s 3S state by electron collisional excitation. The 10830-A line occurs in absorption when the densities in the shocks exceed 10 million/cu cm.
Confinement effects on the solvation structure of solvated alkaline metal cations in a single-digit 1T-MoS 2 nanochannel: A first-principles study
Confinement plays an important role in determining ion transport in porous materials, which, in turn, may influence the performance of many energy storage and desalination devices. In this work, we combined density functional theory (DFT) with an implicit solvation model and ab initio molecular dynamics (AIMD) to investigate the effects of nanoconfinement on several solvated alkaline metal cations in a single-digit 1T-MoS 2 nanochannel. Our DFT calculations with a solvation model indicated that cations with stronger hydration energy introduce a higher number of co-intercalated water molecules into the channel, consistent with early experimental observation obtained for MXene (2D transition metal carbide) channels. The predicted optimal water numbers for the cations were then used for AIMD simulations that explicitly include the effects of the solvent. When compared with the cations in bulk solution, our simulations showed that the hydration structure and coordination number (CN) of the solvated cations confined in the MoS 2 channel can be significantly altered. Further, we found that larger cations with weaker hydration energy (K + , Rb + , and Cs + ) exhibited a distinctive CN decrease under confinement, while smaller cations (Li + and Na + ) retained a similar hydration shell as in the bulk solution. More specifically, the hydration shell of large cations (K + , Rb + , and Cs + ) in MoS 2 showed similar features of the coordination angle to the bulk, which suggests the partially broken hydration shell with no geometry change under confinement. Our simulations provided insights into the change of the hydration structure of alkaline metal cations under confinement, which may have important implications on their transport in the 1T-MoS 2 channel.
Anderson impurity mechanism for a multi-level model in δ-Pu
Abstract Electronic correlations and spin–orbit interactions in plutonium create variations in the bonding behavior of each of its allotropes. In δ -Pu, the 5f electrons lie at the tipping point between itinerant and localized behavior which has made the use of mixed-level models successful in describing its mechanical properties. The mechanism for the emergence of a mixed-level model has not yet been understood. We use a series of density functional theory approximations to understand the interactions that create a mixed-level description of δ -Pu which leads to accurate physical properties. With the intersite interactions present in the hybrid functional, we show that a single 5f electron engages in orbital-selective bonding that can be understood with an Anderson impurity picture. The Anderson model gives us a mechanism to understand how the bonding in δ -Pu evolves as a function of the interactions in the material such that we obtain both the accuracy and physics of the multi-level models from ab initio theory.
Understanding palladium–tellurium cluster formation on WTe2: From a kinetically hindered distribution to thermodynamically controlled monodispersity
Abstract A fundamental understanding of the transition metal dichalcogenide (TMDC)–metal interface is critical for their utilization in a broad range of applications. We investigate how the deposition of palladium (Pd), as a model metal, on WTe2(001), leads to the assembly of Pd into clusters and nanoparticles. Using X-ray photoemission spectroscopy, scanning tunneling microscopy imaging, and ab initio simulations, we find that Pd nucleation is driven by the interaction with and the availability of mobile excess tellurium (Te) leading to the formation of Pd-Te clusters at room temperature. Surprisingly, the nucleation of Pd-Te clusters is not affected by intrinsic surface defects, even at elevated temperatures. Upon annealing, the Pd-Te nanoclusters adopt an identical nanostructure and are stable up to ∼523 K. Density functional theory calculations provide a foundation for our understanding of the mobility of Pd and Te atoms, preferential nucleation of Pd-Te clusters, and the origin of their annealing-induced monodispersity. These results highlight the role the excess chalcogenide atoms may play in the metal deposition process. More broadly, the discoveries of synthetic pathways yielding thermally robust monodispersed nanostructures on TMDCs are critical to the manufacturing of novel quantum and microelectronics devices and catalytically active nano-alloy centers.
Observation of Long-Range Collective Flow in O + O and Ne + Ne Collisions and Implications for Nuclear Structure Studies
The long-range collective flow of particles produced in oxygen-oxygen ( O + O ) and neon-neon ( Ne + Ne ) collisions is measured with the CMS detector at the CERN LHC. The data samples were collected at a center-of-mass energy per nucleon pair of 5.36 TeV, with integrated luminosities of 7 nb - 1 and 0.8 nb - 1 for O + O and Ne + Ne collisions, respectively. Two- and four-particle azimuthal correlations are measured over nearly five units of pseudorapidity. Significant elliptic ( v 2 ) and triangular ( v 3 ) flow harmonics are observed in both systems. The ratios of v n coefficients between Ne + Ne and O + O collisions suggest sensitivity to the intrinsic nuclear structure of the respective Ne and O ions. Hydrodynamic models with ab initio nuclear structure inputs qualitatively reproduce the collision-overlap dependence of both the v n values and the Ne + Ne to O + O ratios. These measurements provide new constraints when modeling the interplay of nuclear structure and collective dynamics in collisions of O 16 and Ne 20 ions, respectively.
Single-Step Mechanism for Regioselective Nitration of 9,10-BN-Naphthalene with Acetyl Nitrate in the Gas Phase
Here, the energetics of the regioselective mononitration of 9,10-BN-naphthalene with acetyl nitrate (H 3 C 2 NO 4 ) were modeled with ab initio simulations in the gas phase and an acetonitrile solvent. The single-electron-transfer (SET) nitration mechanism leading to a σ-complex and a single-step nitration mechanism were modeled. The energy barrier for the single-step mechanism was lower than that for the SET mechanism in the gas phase. However, the two are much more energetically competitive in the solvent. The σ-complex was found to be unstable in the gas phase owing to the interaction with the counterion. Using the single-step mechanism, the carbon site 1 nearest boron had the lowest activation energy for nitration of 22.6 kcal/mol, while site 3 had the second lowest barrier of 24.6 kcal/mol. Finally, details on the molecular structures at intermediate and transition states as well as charges in different configurations are discussed.
Mixed quantum-classical methods for polaron spectral functions
In this work, using two distinct semiclassical approaches—namely, the mean-field Ehrenfest method and the mapping approach to surface hopping—we investigate the spectral function of a single charge interacting with phonons on a lattice. This quantity is relevant for the description of angle-resolved photoemission experiments. Focusing on the one-dimensional Holstein model, we compare the performance of these approaches across a range of coupling strengths and lattice sizes, exposing the relative strengths and weaknesses of each. We demonstrate that these approaches can be efficiently applied with reasonable accuracy to ab initio polaron models. Furthermore, our work provides a route to the calculation of spectral properties in realistic electron–phonon-coupled systems in a computationally inexpensive manner with encouraging accuracy.