Search NASA⌕ Search

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.

At least 109 records · Page 6

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.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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.

2-dimensional systems↗

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.

36 MATERIALS SCIENCE↗

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.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

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.

36 MATERIALS SCIENCE↗

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.

Hayrapetyan, Aram [Yerevan Phys. Inst.]↗

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.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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.

Atomic and molecular spectra↗

Universal Solvent Viscosity Reduction via Hydrogen Bonding Disruptors

Liquid Ion Solutions LLC (DBA RoCo Global) in partnership with Carnegie Mellon University and Carbon Capture Scientific LLC, has performed lab-scale development and evaluation of novel additives that lower the viscosity of water-lean amine solvents for post-combustion carbon dioxide capture. This project focuses on developing additives that minimize the formation of long-range electrostatic and hydrogen bonding (HB) networks, decreasing the solvent viscosity, improving diffusion, and improving the process economics. The project objectives included: 1) performing computer simulation to understand the molecular interactions of the additive molecules in water-lean CO 2 capture solvents, 2) design and synthesis of HB disruptors additives, 3) performance testing with additive molecules on model amine solvents, and 4) demonstration of the effectiveness of the optimized additives in the presence of synthetic flue gas. To meet the abovementioned objectives, the project team utilized a holistic approach that combines molecular simulation, experimental testing, and economic analysis studies. The project team developed ab initio molecular model and then perform computer simulation to develop relationship between hydrogen bonding, viscosity, and performed quantitative analysis of additive on the viscosity of the solvent. The team completed computational comparative study on a range of organic functional groups such as ethers, esters, cyclic carbonates, alkanes, and ammonium salts for their effect on viscosity gaining key insights into molecular interactions and the impact of various functional groups and molecular shapes on viscosity. Assisted with molecular simulation insights, the project team conducted additive synthesis and testing, including a proof-of-concept study, additive screening, optimization, and synthetic flue gas testing. The experimental proof-of-concept study proved that the hydrogen bonding acceptors result in significant decrease of viscosities. Detailed additive screening (exploring various functionalities and molecular structures) has been performed. Several promising additives showed excellent reduction in viscosity (30-41%) at 5% additive loading, and over 50% viscosity reduction at 10% additive loading for the model solvents. The team also performed complex screening studies on additive loadings and mixing effect among additives using the design of experiments. Based on multiple screening experiments, one additive-solvent candidate was down-selected for synthetic flue gas testing. A 100-hour continuous absorption/desorption study was conducted under simulated flue gas using a lab-scale continuous capture and separation system. No degradation (for both solvent and additive) was observed based on the GC results of the solvent samples collected from the continuous study. The team conducted preliminary engineering analyses and cost-benefit analyses to quantify the potential economic benefits of the additive approach for solvent viscosity reduction. Based on the experimental data, CO 2 capture cost savings from the capital and operating cost savings are estimated at $\$$4.7/tonne and $\$$0.3/tonne CO 2 captured, respectively. Considering the additive cost, the net benefit is estimated to be between $\$$4.32~$\$$4.86/tonne CO 2 captured.

20 FOSSIL-FUELED POWER PLANTS↗

Uranium Carbide and Uranium-Metal Thermal Scattering Law Evaluations and Cross Sections [Slides]

New Uranium Carbide and Uranium Metal evaluations have been developed using the FLASSH platform. Ab initio LD/MD models converged electronic, magnetic, and dynamical properties and demonstrate good agreement with experimental data. Thermal scattering cross sections generate inelastic and elastic thermal neutron scattering libraries in FLASSH . New TSL evaluations occurred in the production of ENDF File 7 data for Uranium Carbide and Uranium Metal submitted to NNDC for ENDF/B-VIII.1 release.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Temperature-dependent elastic constants of thorium dioxide probed using time-domain Brillouin scattering

Herein we report the adiabatic elastic constants of single-crystal thorium dioxide over a temperature range of 77–350 K. Time-domain Brillouin scattering, an all-optical, non-contact picosecond ultrasonic technique, is used to generate and detect coherent acoustic phonons that propagate in the bulk perpendicular to the surface of the crystal. These coherent acoustic lattice vibrations have been monitored in two hydrothermally grown single-crystal thorium dioxide samples along the (100) and (311) crystallographic directions. The three independent elastic constants of the cubic crystal (C 11 , C 12 , and C 44 ) are determined from the measured bulk acoustic velocities. The longitudinal wave along the (100) orientation provided a direct measurement of C 11 . Measurement of C 44 and C 12 was achieved by enhancing the intensity of quasi-shear mode in a (311) oriented crystal by adjusting the polarization angle relative to the crystal axes. We find the magnitude of softening of the three elastic constants to be ~2.5% over the measured temperature range. Good agreement is found between the measured elastic constants with previously reported values at room temperature, and between the measured temperature-dependent bulk modulus with calculated values. We find that semi-empirical models capturing lattice anharmonicity adequately reproduce the observed trend. We also determine the acoustic Grüneisen anharmonicity parameter from the experimentally derived temperature-dependent bulk modulus and previously reported temperature-dependent values of volumetric thermal expansion coefficient and heat capacity. This work presents measurements of the temperature-dependent elasticity in single-crystal thorium dioxide at cryogenic temperature and provides a basis for testing ab initio theoretical models and evaluating the impact of anharmonicity on thermophysical properties.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Probing the 𝑗 dependence of angular distributions and 𝑁=20 shell rigidity via the 36 S (𝑝,𝑑) 35 S reaction

An investigation of the N=20 36 S nucleus has been performed through a detailed study of the 3 6S (p,d)3 5S neutron-removal reaction, employing a 66-MeV proton beam at iThemba Laboratory for Accelerator Based Sciences and an innovative target design. Using the high-resolution K = 600 magnetic spectrometer, 98 states in 3 5S were identified up to 16-MeV excitation energy, including 47 previously unobserved states. Angular distributions and spectroscopic factors, including isobaric analog-state contributions, were extracted for 81 levels. A pronounced j dependence in the angular distributions of ℓ=2 states provides refined insights into the spin-orbit splitting. Finite-range adiabatic distorted-wave approximation calculations qualitatively reproduce the observed j dependence. Comparisons of the measured 1d 5/2 spectroscopic strength distribution with large-scale shell-model and ab initio calculations show good agreement overall, and the robustness of the N=20 shell closure in 36 S is confirmed when comparing the relatively low fp orbital occupancies in 40 Ca and 3 6S across the Fermi surface. This study underscores the utility of neutron-removal reactions in probing nuclear structure and the Fermi surface of sd nuclei and beyond. The findings advance our understanding of shell evolution and offer constraining data for theoretical models.

20 ≤ A ≤ 38↗

Initial decomposition mechanisms of 2,4,6-triamino-1,3,5-trinitrobenzene (TATB) and their kinetic isotope effect

2,4,6-triamino-1,3,5-trinitrobenzene (TATB) is an insensitive High Explosive (HE) that is widely studied to better understand the physical properties of safety and sensitivity of HE. A dominant initial decomposition mechanism of TATB is believed to be a dehydration reaction that forms mono- and di-furazans, although other mechanisms have been reported. In this work, seven initial decomposition mechanisms were modeled with ab initio simulations to calculate its free energy barriers, decomposition rates, and kinetic isotope effects. The energy barrier for mono-benzofurazan mechanisms was found to be high, >61 kcal/mol in the gas phase; however, the reaction energy can decrease significantly in a disordered condensed state. The predicted kinetic isotope effect ratio of the furazan mechanism was found to be k h /k d ≈ 1.41 at 60K, in agreement with the experiment. The NO 2 scission mechanism was found to be an entropy-driven mechanism because the free energy barrier decreased significantly with temperature, making it the most energetically favorable mechanism at high temperatures in the gas phase. So the results provide a better understanding of the atomistic decomposition mechanisms of TATB and may be useful for improving models of safety and sensitivity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Water‐mediated Decomposition Pathways of ε‐Hexanitrohexaazaisowurtzitane

The dominant initial decomposition mechanism of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (HNIW) is understood from the literature to be the N─NO 2 bond scission unimolecular mechanism. Here we have found a water-mediated HONO release initial decomposition mechanism energetically competitive with the N─NO 2 bond scission mechanism modeled using ab initio nudged elastic band simulations. The activation energy of the water-mediated HONO mechanism was calculated to be 138.9 kJ/mol with the Perdew, Burke, and Ernzerhof (PBE) functional (177.8 with PBE0), while the double NO 2 release mechanism is 171.9 kJ/mol (203.3 with PBE0). Branching secondary, tertiary, and quaternary decomposition steps were also discovered including oxidation of HONO or H 2 O causing ring-opening that leads to C─N bond breaking and N 2 O or NO release. The reaction rate of HONO oxidation is much faster than HONO release, making the HONO concentration low. The Helmholtz free energy barrier of the double NO 2 release mechanism is lower than the water-mediated HONO release above 438 K (165°C) due to the vibrational contribution to the free energy. This helps explain why the literature reports more NO 2 release at higher temperatures. The present study reveals a physical mechanism for how water can catalyze the decomposition of HNIW at low enough temperatures, and reveals secondary, tertiary, and quaternary reaction mechanisms leading to NO x release. In conclusion, as water is realistically always present, the reported barriers are important for building kinetic models needed to ensure the safety and reliability of HNIW.

Steele, Brad A. [Lawrence Livermore National Labor↗

Microsecond Isomer at the N = 20 Island of Shape Inversion Observed at FRIB

Excited-state spectroscopy from the first experiment at the Facility for Rare Isotope Beams (FRIB) is reported. A 24(2)-μs isomer was observed with the FRIB Decay Station initiator (FDSi) through a cascade of 224- and 401-keV γ rays in coincidence with 32 Na nuclei. This is the only known microsecond isomer (1 μs ≤T 1/2 <1 ms) in the region. This nucleus is at the heart of the N=20 island of shape inversion and is at the crossroads of the spherical shell-model, deformed shell-model, and ab initio theories. It can be represented as the coupling of a proton hole and neutron particle to 32 Mg, 32 Mg+π -1 +ν +1 . This odd-odd coupling and isomer formation provides a sensitive measure of the underlying shape degrees of freedom of 32 Mg, where the onset of spherical-to-deformed shape inversion begins with a low-lying deformed 2 + state at 885 keV and a low-lying shape-coexisting 0$^+_2$ state at 1058 keV. Here, we suggest two possible explanations for the 625-keV isomer in 32 Na: a 6 - spherical shape isomer that decays by E2 or a 0 + deformed spin isomer that decays by M2. The present results and calculations are most consistent with the latter, indicating that the low-lying states are dominated by deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Orientation Control in Epitaxial PdO Thin Films Grown on MgO (001) – Role of Oxygen Chemical Potential

Control of crystal orientations in thin films of functional materials allowsedictive tuning of their strain states, electronic properties, and surface chemical reactivity. Here, conditions for orientation control in epitaxial PdO films are investigated. Due to its tetragonal structure, PdO can form two orientational relationships with the MgO (001). It is shown that, under an oxygen-rich environment provided by oxygen-plasma-assisted molecular beam epitaxy, both (00l)- and (100)-oriented PdO domains form on MgO (001). Subsequent thermal annealing in a vacuum promotes film restructuring to a predominantly (100)-oriented PdO with improved crystallinity. Ab initio calculations reveal that the (001) orientation has lower strain energy but weaker interfacial interactions and serves as an oxygen vacancy sink, whereas the (100) orientation benefits from significantly stronger MgO─PdO bonding. Consequently (100)-oriented domains become favored under oxygen-poor conditions. A mechanism is proposed whereby vacuum annealing drives orientation transformation by generating oxygen vacancies that destabilize the (001) domains and promote (100) ordering. These findings deepen the understanding of how oxygen content impacts interfacial stability and reorganization, thereby offering a route to tune domain orientations in oxide thin films.

36 MATERIALS SCIENCE↗

2D Mg-Cu Intermetallic Compounds with Nontrivial Band Topology and Dirac Nodal Lines

Synthesis of vertical heterostructures that include atomically thin layers of materials with topologically nontrivial energy bands is desirable for exploring exotic quantum states. In this study, the authors report on atomic-layer-by-layer deposition of magnesium on copper(111) surface by molecular beam epitaxy, monitored in situ by low-energy electron microscopy and diffraction, and modeled by ab initio theory. It is found that a 2D MgCu 2 intermetallic compound forms during initial Mg deposition and persists till a full monolayer of Mg is formed. Deposition of additional Mg triggers a phase transition from the commensurate MgCu 2 to an incommensurate Mg 2 Cu layer and enables growth of the second layer of Mg. Ab initio calculations indicate non-trivial topology of the electron bands in the interfacial Mg 2 Cu layer and the existence of Dirac nodal lines near the Fermi level. The new 2D Mg 2 Cu material emerges as a promising platform to study new topological states of matter.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

High‐Capacity Splitting of Mono‐ and Dibranched Hexane Isomers by a Robust Zinc‐Based Metal–Organic Framework

Abstract High‐efficiency separation of C 6 alkanes, particularly the mono‐ and dibranched isomers by using porous solids, is of paramount significance in the petrochemical industry and, remains a daunting challenge. In this work, we report the complete separation of linear/monobranched hexanes from their dibranched isomers through selective size‐exclusion by a microporous MOF, Zn‐tcpt (H 3 tcpt=2,4,6‐tris(4‐carboxyphenoxy)‐1,3,5‐triazine), with a two‐fold interpenetrated structure of hms nets. Importantly, its adsorption capacity and selectivity are notably higher than those of the previously reported adsorbents that can split mono‐ and dibranched alkane isomers. Dynamic breakthrough measurements verify the excellent separation of C 6 alkane isomers by Zn‐tcpt, and the size‐exclusion based separation mechanism has been confirmed by ab initio materials modeling. The high‐efficiency separation of alkane isomers by Zn‐tcpt can be attributed to its optimal pore dimensions as well as high porosity.

Yu, Liang↗