Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio calculation”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

HTESP (High-throughput electronic structure package): A package for high-throughput ab initio calculations

High-throughput ab initio calculations are the indispensable parts of data-driven discovery of new materials with desirable properties, as reflected in the establishment of several online material databases. The accumulation of extensive theoretical data through computations enables data-driven discovery by constructing machine learning and artificial intelligence models to predict novel compounds and forecast their properties. Efficient usage and extraction of data from these existing online material databases can accelerate the next stage materials discovery that targets different and more advanced properties, such as electron–phonon coupling for phonon-mediated superconductivity. However, extracting data from these databases, generating tailored input files for different ab initio calculations, performing such calculations, and analyzing new results can be demanding tasks. Here, in this work, we introduce a software package named “HTESP” (High-Throughput Electronic Structure Package) written in Python and Bash languages, which automates the entire workflow including data extraction, input file generation, calculation submission, result collection and plotting. Our HTESP will help speed up future computational materials discovery processes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Ab initio calculations of scattering and fusion reactions

We review recent advancements achieved within the ab initio no-core shell model with continuum (NCSMC), a unified first-principle (or, ab initio) approach to nuclear bound and continuum states. We highlight its application to reactions of astrophysical interests, including the 7 Be(p, γ) 8 B solar fusion and α-α scattering, and to the structure of the exotic 12 Be nucleus. Furthermore, our findings provide a solid foundation for integrating ab initio calculations with experimental measurements, allowing for more precise evaluations of key thermonuclear reaction rates at astrophysical energies, and offer new insights into the factors that contribute to the emergence of clustering.

Quaglioni, Sofia [Lawrence Livermore National Labo↗
DOE OSTI2022

Converged ab initio calculations of heavy nuclei

We propose a novel storage scheme for three-nucleon (3N) interaction matrix elements relevant for the normal-ordered two-body approximation used extensively in ab initio calculations of atomic nuclei. This scheme reduces the required memory by approximately two orders of magnitude, which allows the generation of 3N interaction matrix elements with the standard truncation of E 3max =28, well beyond the previous limit of 18. We demonstrate that this is sufficient to obtain the ground-state energy of 132 Sn converged to within a few MeV with respect to the E 3max truncation. In addition, we study the asymptotic convergence behavior and perform extrapolations to the un-truncated limit. Finally, we investigate the impact of truncations made when evolving free-space 3N interactions with the similarity renormalization group. We find that the contribution of blocks with angular momentum J rel > 9/2 to the ground-state energy is dominated by a basis-truncation artifact, which vanishes in the large-space limit, so these computationally expensive components can be neglected. For the two sets of nuclear interactions employed in this work, the resulting binding energy of 132 Sn agrees with the experimental value within theoretical uncertainties. This work enables converged ab initio calculations of heavy nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Energy and mechanical properties predictions in Fe-Ni binary system by ab initio calculations

Here, the fcc random solid solutions in Fe-Ni binary system were investigated using ab initio calculations. Density functional theory (DFT) was applied along with supercells generated by the Special Quasi-Random Structure (SQS) method to predict the enthalpy of mixing and magnetic and elastic properties. The results were compared with relevant experimental and computational data produced by other groups. The effect of SQS supercell size and relaxation procedure on the predicted properties of the alloys were also studied. Our calculations show that 16-atoms SQS model provides adequate approach to the mixing energy of the alloy. However, calculation of mechanical properties requires larger supercell. Excellent agreement between experimental and calculated elastic properties was found in fcc region of Fe-Ni system.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Temperature-dependent phonon lifetimes and thermal conductivity of silicon by inelastic neutron scattering and ab initio calculations

In this paper, inelastic neutron scattering on a single crystal of silicon was performed at temperatures from 100 to 1500 K. These experimental data were reduced to obtain phonon spectral intensity at all wave vectors $\overset{→}{Q}$ and frequencies ω in the first Brillouin zone. Thermal broadenings of the phonon peaks were obtained by fitting and by calculating with an iterative ab initio method that uses thermal atom displacements on an ensemble of superlattices. Agreement between the calculated and experimental broadenings was good, with possible discrepancies at the highest temperatures. Distributions of phonon widths versus phonon energy had similar shapes for computation and experiment. These distributions grew with temperature but maintained similar shapes. Parameters from the ab initio calculations were used to obtain the thermal conductivity from the Boltzmann transport equation, which was in good agreement with experimental data. Despite the high group velocities of longitudinal acoustic phonons, their shorter lifetimes reduced their contribution to the thermal conductivity, which was dominated by transverse acoustic modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Ab Initio Calculation of the Alpha-Particle Monopole Transition Form Factor

Here we present a parameter-free ab initio calculation of the α-particle monopole transition form factor in the framework of nuclear lattice effective field theory. We use a minimal nuclear interaction that was previously used to reproduce the ground state properties of light nuclei, medium-mass nuclei, and neutron matter simultaneously with no more than a few percent error in the energies and charge radii. The results for the monopole transition form factor are in good agreement with recent precision data from Mainz.

4He↗
DOE OSTI2020

Ab initio calculations of the rate of carrier trapping and release at dopant sites in NaI: Tl beyond the harmonic approximation

We present ab initio calculations of electron capture and release coefficients as a function of temperature for Tl dopants in the widely used scintillator NaI, which is a soft ionic crystal. The modelled capture and release events occur by transitions mediated by multiple phonon absorption and emission between states with significantly different local geometries around the trapping site. We demonstrate that such transitions are not well-described by the normal harmonic approximation to the nuclear dynamics. We avoid the harmonic approximation by numerically solving the vibrational Schrödinger equation for the motion along/of a single effective phonon coordinate. The localized trapped state is not correctly described by semilocal density functionals, so we employ a hybrid functional tuned to reproduce the band gap and lattice constant of the host material. Our calculations, which combine plane-wave and cluster electronic structure methods in a novel way, are in reasonable agreement with available experimental data for detrapping and predict an unusual temperature dependence for trapping rates.

Prange, Micah P.↗
DOE OSTI2021

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.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Ab initio calculations of the carbon and oxygen isotopes: Energies, correlations, and superfluid pairing

We perform ab initio nuclear lattice calculations of the neutron-rich carbon and oxygen isotopes using high-fidelity chiral interactions. We find good agreement with the observed binding energies and compute correlations associated with each two-nucleon interaction channel. For the isospin T = 1 channels, we show that the dependence on T z provides a measure of the correlations among the extra neutrons in the neutron-rich nuclei. For the spin-singlet S-wave channel, we observe that any paired neutron interacts with the nuclear core as well as its neutron pair partner, while any unpaired neutron interacts primarily with only the nuclear core. For the other partial waves, the correlations among the extra neutrons grow more slowly and smoothly with the number of neutrons. These general patterns are observed in both the carbon and oxygen isotopes and may be universal features that appear in many neutron-rich nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Ab initio calculations of overlap integrals for μ → e conversion in nuclei

The rate for μ → e conversion in nuclei is set to provide the most stringent test of lepton-flavor symmetry and a window into physics beyond the Standard Model. However, to disentangle new lepton-flavor-violating interactions, in combination with information from μ → ey and μ → 3e, it is critical that uncertainties at each step of the analysis be controlled and fully quantified. In this regard, nuclear response functions related to the coupling to neutrons are notoriously problematic, since they are not directly constrained by experiment. We address these shortcomings by combining ab initio calculations with a recently improved determination of charge distributions from electron scattering by exploiting strong correlations among charge, point-proton, and point-neutron radii and densities. We present overlap integrals for 27 Al, 48 Ca, and 48 Ti including full covariance matrices, allowing, for the first time, for a comprehensive consideration of nuclear structure uncertainties in the interpretation of μ → e experiments.

ab initio calculations↗
DOE OSTI2021

Nuclear energy density functionals grounded in ab initio calculations

Here, we discuss the construction of a nuclear energy density functional (EDF) from ab initio computations and advocate the need for a methodical approach that is free from ad hoc assumptions. The equations of state (EoSs) of symmetric nuclear and pure neutron matter are computed using the chiral NNLO sat and the phenomenological AV4' + UIX c Hamiltonians as inputs to self-consistent Green's function (SCGF) and auxiliary field diffusion Monte Carlo (AFDMC) methods. We propose a convenient parametrization of the EoS as a function of the Fermi momentum and fit it on the SCGF and AFDMC calculations. We apply the ab initio based EDF to carry out an analysis of the binding energies and charge radii of different nuclei in the local density approximation. The NNLO sat -based EDF produces encouraging results, whereas the AV4' + UIX c -based one is farther from experiment. Possible explanations of these different behaviors are suggested, and the importance of gradient and spin-orbit terms is analyzed. Our paper paves the way for a practical and systematic way to merge ab initio nuclear theory and density functional theory, while shedding light on some critical aspects of this procedure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Anharmonic multi-phonon nonradiative transition: An ab initio calculation approach

Nonradiative carrier recombinations at deep centers in semiconductors are of great importance for both fundamental physics and device engineering. Here, we provide a revised analysis of Huang’s original nonradiative multi-phonon (NMP) theory with ab initio calculations. First, we confirmed at the first-principles level that Huang’s concise formula gives the same results as the matrix-based formula, and that Huang’s high-temperature formula provides an analytical expression for the coupling constant in Marcus theory. Secondly, we correct for anharmonic effects by taking into account local phonon-mode variations for different charge states of a defect. The corrected capture rates for defects in GaN and SiC agree well with experiments.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

High-Precision Spectroscopy of 20 O Benchmarking Ab Initio Calculations in Light Nuclei

The excited states of unstable 20 O were investigated via γ-ray spectroscopy following the 19 O(d, p) ⁢20 O reaction at 8 A⁢ MeV. By exploiting the Doppler shift attenuation method, the lifetimes of the $2^+_2$ and $3^+_1$ states were firmly established. From the γ-ray branching and E⁢2/M⁢1 mixing ratios for transitions deexciting the $2^+_2$ and $3^+_1$ states, the B⁡(E⁢2) and B⁡(M⁢1) were determined. Various chiral effective field theory Hamiltonians, describing the nuclear properties beyond ground states, along with a standard USDB interaction, were compared with the experimentally obtained data. Such a comparison for a large set of γ-ray transition probabilities with the valence space in medium similarity renormalization group ab initio calculations was performed for the first time in a nucleus far from stability. It was shown that the ab initio approaches using chiral effective field theory forces are challenged by detailed high-precision spectroscopic properties of nuclei. Finally, the reduced transition probabilities were found to be a very constraining test of the performance of the ab initio models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Ab initio calculations of third-order elastic coefficients

Third-order elasticity (TOE) theory predicts strain-induced changes in second-order elastic coefficients (SOECs) and can model elastic wave propagation in stressed media. Although third-order elastic tensors have been determined based on first principles in previous studies, their current definition is based on an expansion of thermodynamic energy in terms of the Lagrangian strain near the natural, or zero pressure, reference state. This definition is inconvenient for predictions of SOECs under significant initial stresses. Therefore, when TOE theory is necessary to study the strain dependence of elasticity, the seismological community has resorted to an empirical version of the theory. This study reviews the thermodynamic definition of the third-order elastic tensor and proposes using an “effective” third-order elastic tensor. An explicit expression for the effective third-order elastic tensor is given and verified. Additionally, we extend the ab initio approach to calculate third-order elastic tensors under finite pressure and apply it to two cubic systems, namely, NaCl and MgO. As applications and validations, we evaluate (a) strain-induced changes in SOECs and (b) pressure derivatives of SOECs based on ab initio calculations. Good agreement between third-order elasticity-based predictions and numerically calculated values confirms the validity of our theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2023

Ab initio calculation of carrier mobility in semiconductors including ionized-impurity scattering

The past decade has seen the emergence of ab initio computational methods for calculating phonon-limited carrier mobilities in semiconductors with predictive accuracy. More realistic calculations ought to take into account additional scattering mechanisms such as, for example, impurity and grain-boundary scattering. Here, in this paper, we investigate the effect of ionized-impurity scattering on the carrier mobility. We model the analytical impurity potential parameterized from first principles by a collection of randomly distributed Coulomb scattering centers, and we include this relaxation channel into the ab initio Boltzmann transport equation, as implemented in the EPW code. We demonstrate this methodology by considering silicon, silicon carbide, and gallium phosphide, for which detailed experimental data are available. Our calculations agree well with experiments over a broad range of temperatures and impurity concentrations. For each compound investigated here, we compare the relative importance of electron-phonon scattering and ionized-impurity scattering, and we critically assess the reliability of Matthiessen's rule. We also show that an accurate description of dielectric screening and carrier effective masses can improve quantitative agreement with experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Revisit the VEC criterion in high entropy alloys (HEAs) with high-throughput ab initio calculations: A case study with Al-Co-Cr-Fe-Ni system

Valence electron concentration (VEC) was treated as a useful parameter to predict the stability of solid solution phases. However, the available experimental data to support this criterion is far from enough. In the current study, the high-throughput ab initio modeling is applied to investigate the relative stability of FCC and BCC single crystals of the Al-Co-Cr-Fe-Ni high entropy alloys (HEAs) by using the special quasi-random structure (SQS) approach. Furthermore, the predictions start with pure elements of the Al-Co-Cr-Fe-Ni system and are continued with binaries, ternaries, and quaternary compositions, which come up with 180 compositions (360 structures). After that, the reliability of the VEC criterion is testified. The results show that the VEC criterion not only works for the stable structure but also works effectively for metastable structure when both FCC and BCC are not thermodynamic stable. However, it is found that the old VEC criterion proposed by Guo et al. fails to work effectively for compositions containing high concentrations of light-weight metals such as Al at VEC<5. To solve this problem, the present work proposes a new VEC rule to define the stability of FCC and BCC structures at the ground state. With the implementation of the new VEC rule, the effectiveness of the VEC rule (EVEC) of both FCC and BCC structures is enhanced, especially for pure elements and binary compositions, indicating that this rule does not only work effectively for multicomponent systems but also works for low-order systems.

36 MATERIALS SCIENCE↗