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

Influence of solution phase environmental heterogeneity and fluctuations on vibronic spectra: Perylene diimide molecular chromophore complexes in solution

Ensembles of ab initio parameterized Frenkel-exciton model Hamiltonians for different perylene diimide dimer systems are used, together with various dissipative quantum dynamics approaches, to study the influence of the solvation environment and fluctuations in chromophore relative orientation and packing on the vibronic spectra of two different dimer systems: a π-stacked dimer in aqueous solution in which the relative chromophore geometry is strongly confined by a phosphate bridge and a side-by-side dimer in dichloromethane involving a more flexible alkyne bridge that allows quasi-free rotation of the chromophores relative to one another. These entirely first-principles calculations are found to accurately reproduce the main features of the experimental absorption spectra, providing a detailed mechanistic understanding of how the structural fluctuations and environmental interactions influence the vibronic dynamics and spectroscopy of solutions of these multi-chromophore complexes.

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

High-resolution infrared spectroscopy of supersonically cooled singlet carbenes: Bromomethylene (HCBr) in the CH stretch region

First high-resolution spectra of cold (~35 K) singlet bromomethylene HCBr in the CH stretching (v 1 ) region from 2770 to 2850 cm-1 are reported using near quantum shot-noise limited laser absorption methods in a slit jet supersonic discharge expansion source. Three rovibrational bands are identified at high S/N (20:1 - 40:1) and rotationally assigned to i) the CH stretch fundamental (v 1 ) band $\tilde{X}$(1,0,0) ← $\tilde{X}$ (0,0,0) and ii) vibrational hot bands ($\tilde{X}$(1,1,0) ← $\tilde{X}$(0,1,0) and $\tilde{X}$(1,0,1) ← $\tilde{X}$(0,0,1)) arising from vibrationally excited HCBr populated in the discharge with single quanta in either the H–C–Br bend (v 2 ) or C–Br stretch (v 3 ) modes. Precision rotational constants are reported for a total of six states, with an experimentally determined CH stretch vibrational frequency (2799.38 cm -1 ) in good agreement with previous low-resolution fluorescence studies. Detailed analysis of the fundamental v 1 band highlights the presence of perturbations in the $\tilde{X}$(1,0,0) level, which we tentatively attribute to arise from the nearby triplet state $\tilde{a}$(0,0,1) through spin-orbit interaction or the multiple quanta $\tilde{X}$(0,2,1) singlet state via c-type Coriolis coupling. Reduced-Doppler resolution (60 MHz) in the slit-jet IR spectrometer permits clear observation of nuclear spin hyperfine structure, with experimental line shapes well reproduced by nuclear quadrupole/spin-rotation coupling constants from microwave studies. Lastly, the a-type to b-type transition intensity ratio for the fundamental CH stretch band is notably larger than predicted by a bond-dipole model, which from high level ab initio quantum calculations (CCSD(T)/PVQZ) can be attributed to vibrationally induced “charge-sloshing” of electron density along the polar C–Br bond.

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

Γ valley transition metal dichalcogenide moiré bands

Significance At small twist angles, moiré patterns form in transition metal dichalcogenide bilayers. Using sophisticated ab initio calculations, we develop a model able to describe the band structure of these systems. We found that the moiré pattern acts as a long-wavelength modulating potential, with a sixfold rotational symmetry not present in the underlying lattice. The modulation potential reduces hole band widths and localizes electrons or holes around honeycomb or kagome lattice sites.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Visualizing the interplay of Dirac mass gap and magnetism at nanoscale in intrinsic magnetic topological insulators

In intrinsic magnetic topological insulators, Dirac surface-state gaps are prerequisites for quantum anomalous Hall and axion insulating states. Unambiguous experimental identification of these gaps has proved to be a challenge, however. Here, we use molecular beam epitaxy to grow intrinsic MnBi 2 Te 4 thin films. Using scanning tunneling microscopy/spectroscopy, we directly visualize the Dirac mass gap and its disappearance below and above the magnetic order temperature. We further reveal the interplay of Dirac mass gaps and local magnetic defects. We find that, in high defect regions, the Dirac mass gap collapses. Ab initio and coupled Dirac cone model calculations provide insight into the microscopic origin of the correlation between defect density and spatial gap variations. Finally, this work provides unambiguous identification of the Dirac mass gap in MnBi 2 Te 4 and, by revealing the microscopic origin of its gap variation, establishes a material design principle for realizing exotic states in intrinsic magnetic topological insulators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Atomic structure of an FeCrMoCBY metallic glass revealed by high energy x-ray diffraction

Amorphous bulk metallic glasses with the composition Fe 48 Cr 15 Mo 14 C 15 B 6 Y 2 have been of interest due to their special mechanical and electronic properties, including corrosion resistance, high yield-strength, large elasticity, catalytic performance, and soft ferromagnetism. Here, we apply a reverse Monte Carlo technique to unravel the atomic structure of these glasses. The pair-distribution functions for various atomic pairs are computed based on the high-energy x-ray diffraction data we have taken from an amorphous sample. Monte Carlo cycles are used to move the atomic positions until the model reproduces the experimental pair-distribution function. The resulting fitted model is consistent with our ab initio simulations of the metallic glass. Our study contributes to the understanding of functional properties of Fe-based bulk metallic glasses driven by disorder effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Intertwined charge, spin, and pairing orders in doped iron ladders

Motivated by recent experimental progress on iron-based ladder compounds, we study the doped two-orbital Hubbard model for the two-leg ladder BaFe2S3. The model is constructed by using ab initio hopping parameters and the ground state properties are investigated using the density matrix renormalization group method. We show that the (π,0) magnetic ordering at half filling, with ferromagnetic rungs and antiferromagnetic legs, becomes incommensurate upon hole doping. Moreover, depending on the strength of the Hubbard U coupling, other magnetic patterns, such as (0,π), are also stabilized. We found that the binding energy for two holes becomes negative for intermediate Hubbard interaction strength, indicating hole pairing. Due to the crystal-field split among orbitals, the holes primarily reside in one orbital, with the other one remaining half filled. This resembles orbital selective Mott states. The formation of tight hole pairs continues with increasing hole density, as long as the magnetic order remains antiferromagnetic in one direction. The study of pair-pair correlations indicates the dominance of the intraorbital spin-singlet channel, as opposed to other pairing channels. Although in a range of hole doping pairing correlations decay slowly, our results can also be interpreted as corresponding to a charge density wave made of pairs, a precursor of eventual superconductivity after interladder couplings are included. Such a scenario of intertwined orders has been extensively discussed before in the cuprates, and our results suggest a similar physics could exist in ladder iron-based superconductors. Finally, we also show that a robust Hund's coupling is needed for pairing to occur.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Layer-dependent optically induced spin polarization in InSe

Optical control of spin in semiconductors has been pioneered using nanostructures of III-V and II-VI semi-conductors, but the emergence of two-dimensional van der Waals materials offers an alternative low-dimensional platform for spintronic phenomena. Indium selenide (InSe), a group-III monochalcogenide van der Waals material, has shown promise for optoelectronics due to its high electron mobility, tunable direct band gap, and quantum transport. In addition to these confirmed properties, there are predictions of spin-dependent optical selection rules, suggesting the potential for all-optical excitation and control of spin in a two-dimensional layered material. Despite these predictions, layer-dependent optical spin phenomena in InSe have yet to be explored. Here, we present measurements of layer-dependent optical spin dynamics in few-layer and bulk InSe. Polarized photoluminescence reveals layer-dependent optical orientation of spin, thereby demonstrating the optical selection rules in few-layer InSe. Spin dynamics are also studied in many-layer InSe using time-resolved Kerr rotation spectroscopy. By applying out-of-plane and in-plane static magnetic fields for polarized emission measurements and Kerr measurements, respectively, the g factor for InSe was extracted. Further investigations are done by calculating precession values using a k · p model, which is supported by ab initio density functional theory. Comparison of predicted precession rates with experimental measurements highlights the importance of excitonic effects in InSe for understanding spin dynamics. In conclusion, optical orientation of spin is an important prerequisite for optospintronic phenomena and devices, and these first demonstrations of layer-dependent optical excitation of spins in InSe lay the foundation for combining layer-dependent spin properties with advantageous electronic properties found in this material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Using spin alignment of inelastically excited nuclei in fast beams to assign spins: The spectroscopy of 13 O as a test case

Excited states in O 13 were investigated using inelastic scattering of an E / A = 69.5 MeV O 13 beam off of a Be 9 target. The excited states were identified in the invariant-mass spectra of the decay products. Both single-proton and sequential two-proton decays of the excited states were examined. Furthermore, for a number of the excited states, the protons were emitted with strong anisotropy where emissions transverse to the beam axis are favored. The measured proton-decay angular distributions were compared to predictions from distorted-wave Born-approximation calculations of the spin alignment which was shown to be largely independent of the excitation mechanism. The deduced O 13 level scheme is compared to ab initio no-core shell model with continuum predictions. The lowest-energy excited states decay isotropically consistent with predictions of strong proton 1 s 1 / 2 structure. Above these states in the level scheme, we observed a number of higher-spin states not predicted within the model. Possibly these are associated with rotational bands built on deformed cluster configurations predicted by antisymmetrized molecular dynamics calculations. The spin alignment mechanism is shown to be useful for making spin assignments and may have widespread use.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Muon capture on Li 6 , C 12 , and O 16 from ab initio nuclear theory

Muon capture on nuclei is one of the most promising probes of the nuclear electroweak current driving the yet-hypothetical neutrinoless double-beta (0νββ) decay. Both processes involve vector and axial-vector currents at finite momentum transfer, q ~ 100 MeV, as well as the induced pseu doscalar and weak-magnetism currents. Comparing measured muon-capture rates with reliable ab initio nuclear-theory predictions could help us validate these currents. To this end, we compute partial muon-capture rates for 6 Li, 12 C and 16 O, feeding the ground and excited states in 6 He, 12 B and 16 N, using ab initio no-core shell model with two- and three-nucleon chiral interactions. Here, we remove the spurious center-of-mass motion by introducing translationally invariant operators and approximate the effect of hadronic two-body currents by Fermi-gas model. We solve the bound-muon wave function from the Dirac wave equations in the Coulomb field created by a finite nucleus. We find that the computed rates to the low-lying states in the final nuclei are in good agreement with the measured counterparts. We highlight sensitivity of some of the transitions to the sub-leading three-nucleon interaction terms. We also compare summed rates to several tens of final states with the measured total capture rates and note that we slightly underestimate the total rate with this simple approach due to limited range of excitation energies

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Ab initio investigation of the Li 7 ( p , e + e - ) Be 8 process and the X17 boson

Observations of anomalies in the electron-positron angular correlations in high-energy decays in 4 He, 8 Be, and 12 C have been reported recently by the ATOMKI collaboration. These could be explained by the creation and subsequent decay of a new boson with a mass of ≈ 17MeV. Theoretical understanding of pair creation in the proton capture reactions used in these experiments is important for the interpretation of the anomalies. We apply the ab initio no-core shell model with continuum (NCSMC) to the proton capture on 7 Li. The NCSMC describes both bound and unbound states in light nuclei in a unified way with chiral two- and three-nucleon interactions as the only input. We investigate the structure of 8 Be, the p+ 7 Li elastic scattering, the 7 Li(p,y)⁢ 8 Be cross section, and the internal pair creation 7 Li ⁢(p,e + ⁢e - ) 8 Be. Here we discuss the impact of a proper treatment of the initial scattering state on the electron-positron angular correlation spectrum and compare our results to available ATOMKI data sets. Finally, we calculate 7 Li ⁢(p,X)⁢ 8 Be cross sections for several proposed models of the hypothetical X17 particle.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Elastic Antiproton-Nucleus Scattering from Chiral Forces

Elastic scattering of antiprotons off He 4 , C 12 , and O 16 , 18 is described for the first time with a consistent microscopic approach based on the calculation of an optical potential (OP) describing the antiproton-target interaction. The OP is derived using the recent antiproton-nucleon ( p ¯ N ) chiral interaction to calculate the p ¯ N t matrix, while the target densities are computed with the ab initio no-core shell model using chiral interactions as well. Here, our results are in good agreement with the existing experimental data and the results computed at different chiral orders of the p ¯ N interaction display a well-defined convergence pattern.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Ab Initio Prediction of the 4 $\mathrm{H}$e($d, γ$) 6 $\mathrm{L}$i Big Bang Radiative Capture

We report the rate at which helium ( 4 He ) and deuterium (d) fuse together to produce lithium-6 ( 6 Li) and a γ ray, 4 $\mathrm{H}$e($d, γ$) 6 $\mathrm{L}$i, is a critical puzzle piece in resolving the discrepancy between big bang predictions and astronomical observations for the primordial abundance of 6 Li . The accurate determination of this radiative capture rate requires the quantitative and predictive description of the fusion probability across the big bang energy window (30 keV ≲ E ≲ 400 keV), where measurements are hindered by low counting rates. We present first-principle (or, ab initio) predictions of the 4 $\mathrm{H}$e($d, γ$) 6 $\mathrm{L}$i astrophysical S factor using validated nucleon-nucleon and three-nucleon interactions derived within the framework of chiral effective field theory. By employing the ab initio no-core shell model with continuum to describe 4 He-d scattering dynamics and bound 6 Li product on an equal footing, we accurately and consistently determine the contributions of the main electromagnetic transitions driving the radiative capture process. Our results reveal an enhancement of the capture probability below 100 keV owing to previously neglected magnetic dipole (M1) transitions and reduce by an average factor of 7 the uncertainty of the thermonuclear capture rate between 0.002 and 2 GK.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Electronic and structural properties of RbCe X 2 ( X 2 : O 2 , S 2 , SeS, Se 2 , TeSe, Te 2 )

We report triangular lattice delafossite compounds built from magnetic lanthanide ions are a topic of recent interest due to their frustrated magnetism and realization of quantum disordered magnetic ground states. Here we report the evolution of the structure and electronic ground states of RbCeX 2 compounds, built from a triangular lattice of Ce 3+ ions, upon varying their anion character (X 2 =O 2 , S 2 , SeS, Se 2 , TeSe, Te 2 ). This includes the discovery of a new member of this series, RbCeO 2 , that potentially realizes a quantum disordered ground state analogous to NaYbO 2 . Magnetization and susceptibility measurements reveal that all compounds manifest mean-field antiferromagnetic interactions and, with the exception of the oxide, possess signatures of magnetic correlations onset below 1 K. The crystalline electric field level scheme is explored via neutron scattering and ab initio calculations in order to model the intramultiplet splitting of the J=5/2 multiplet. In addition to the two excited doublets expected within the J=5/2 manifold, we observe one extra local mode present across the sample series. This added mode shifts downward in energy with increasing anion mass and decreasing crystal field strength, suggesting a long-lived anomalous mode endemic to anion motion about the Ce3 + sites.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

First-principles ionized-impurity scattering and charge transport in doped materials

Scattering of carriers with ionized impurities governs charge transport in doped semiconductors. However, electron interactions with ionized impurities cannot be fully described with quantitative first-principles calculations, so their understanding relies primarily on simplified models. Here we show an ab initio approach to compute the interactions between electrons and ionized impurities or other charged defects. It includes the short- and long-range electron-defect (e-d) interactions on equal footing and allows for efficient interpolation of the e-d matrix elements. Here we combine the e-d and electron-phonon interactions in the Boltzmann transport equation to compute the carrier mobilities in doped silicon over a wide range of temperature and doping concentrations, seamlessly spanning the defect- and phonon-limited transport regimes. The individual contributions of the defect- and phonon-scattering mechanisms to the carrier relaxation times and mean-free paths are analyzed. Our method provides a powerful tool to study electronic interactions in doped materials. It broadens the scope of first-principles transport calculations, enabling studies of a wide range of doped semiconductors and oxides with application to electronics, energy, and quantum technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2021

PyCrystalField : software for calculation, analysis and fitting of crystal electric field Hamiltonians

PyCrystalField is a Python software package for calculating single-ion crystal electric field (CEF) Hamiltonians. This software can calculate a CEF Hamiltonian ab initio from a point charge model for any transition or rare earth ion in either the J basis or the LS basis, perform symmetry analysis to identify nonzero CEF parameters, calculate the energy spectrum and observables such as neutron spectrum and magnetization, and fit CEF Hamiltonians to any experimental data. Here, the theory, implementation and examples of its use are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Direct Air Reactive Capture and Conversion for Utility-Scale Energy Storage (Final Report)

This final report for FEW0277 summarizes the work performed over the project performance period of October 2021 – March 2025. This project was funded under the “Reactive Capture and Conversion R&D” lab call released in FY2021. The goal of the project was to develop dual-function materials and process for capturing CO 2 from the atmosphere and converting it into CH 4 . The work was organized into four parallel tracks in 1) direct air capture materials synthesis and characterization, 2) catalysts for CO 2 conversion, 3) mechanistic investigations via ab initio simulations, and 4) process modeling, technoeconomic analysis, and lifecycle assessment. The project was split into two budget periods. The first budget period focused on development of amine-based materials, due to their known performance for CO 2 direct air capture and their potential to act synergistically with metal catalysts to enable a low-temperature methanation pathway. The second budget period focused on development of alkali-based materials and a simulated-moving-bed process for high conversion catalytic reduction of captured CO 2 to CH 4 . All project milestones were completed during the project performance period and are summarized in this report. Our work resulted in publication of eight peer-reviewed manuscripts, one patent application, and numerous presentations given at domestic and international conferences and invited academic department seminars.

03 NATURAL GAS↗
DOE OSTI2022

$O(N)$ ab initio calculation scheme for large-scale moiré structures

Here we present a two-step method specifically tailored for band structure calculation of the small-angle moiré-pattern materials which contain tens of thousands of atoms in a unit cell. In the first step, the self-consistent field calculation for the ground state is performed with the O(N) Krylov subspace method implemented in openmx. Second, the crystal momentum-dependent Bloch Hamiltonian and overlap matrix are constructed from the results obtained in the first step and only a small number of eigenvalues near the Fermi energy are solved with shift-invert and Lanczos techniques. By systematically tuning two key parameters, the cutoff radius for electron hopping interaction and the dimension of the Krylov subspace, we obtained the band structures for both rigid and corrugated twisted bilayer graphene structures down to the first magic angle (θ = 1.08°) with high enough accuracy at affordable costs. The band structures are in good agreement with those from tight-binding models, continuum models, plane-wave pseudopotential based ab initio calculations, and experimental observations. This method is also shown to be efficient in twisted double-bilayer graphene and bilayer WSe 2 . We think this two-step method can play a crucial role in other twisted two-dimensional materials, especially those with much more complex band structure and where the effective model is hard to construct.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Protein C-GeM: A Coarse-Grained Electron Model for Fast and Accurate Protein Electrostatics Prediction

The electrostatic potential (ESP) is a powerful property for understanding and predicting electrostatic charge distributions that drive interactions between molecules. In this study, we compare various charge partitioning schemes including fitted charges, density-based quantum mechanical (QM) partitioning schemes, charge equilibration methods, and our recently introduced coarse-grained electron model, C-GeM, to describe the ESP for protein systems. Furthermore, when benchmarked against high quality density functional theory calculations of the ESP for tripeptides and the crambin protein, we find that the C-GeM model is of comparable accuracy to ab initio charge partitioning methods, but with orders of magnitude improvement in computational efficiency since it does not require either the electron density or the electrostatic potential as input.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗