Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Calculations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,351 records · Page 75

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↗

Calculations of anisotropic magnetic properties using spin-orbit energy variations

We analyze several methods of obtaining the accurate relativistic total energy (TE) variations using traditional perturbation theories (PTs) and proposed coupling constant integration (CCI) methods. For this purpose, we perform benchmark calculations within the density functional theory taking the spin-orbit coupling (SOC) and its derivative as a perturbation. The TE change due to SOC addition obtained from both PTs and CCI is shown to reach the accuracy of fully self-consistent TE calculations. Similar accuracy is also obtained even for the magnetocrystalline anisotropy energy (MAE). The real advantage of the proposed methods is to use PTs and CCI methods in those electronic structure methods where accurate total energies currently cannot be obtained with required accuracy. Correspondingly, we demonstrate the applicability of suggested methods for calculations of MAE in different magnetic materials using a dynamic mean-field method. All suggested PTs and CCI methods also provide convenient site, orbital, and spin decompositions of the TE variation, creating a powerful way to analyze microscopic physics in strongly correlated materials.

36 MATERIALS SCIENCE↗

First-principles calculations of metal surfaces. I. Slab-consistent bulk reference for convergent surface properties

The first-principles computation of the surfaces of metals is typically accomplished through slab calculations of finite thickness. The extraction of a convergent surface formation energy from slab calculations is dependent upon defining an appropriate bulk reference energy. I describe a method for an independently computed, slab-consistent bulk reference that leads to convergent surface formation energies from slab calculations that also provides realistic uncertainties for the magnitude of unavoidable nonlinear divergence in the surface formation energy with slab thickness. In conclusion, the accuracy is demonstrated on relaxed, unreconstructed low-index aluminum surfaces with slabs with up to 35 layers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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↗

Analytic calculation of the vison gap in the Kitaev spin liquid

Although the ground-state energy of the Kitaev spin liquid can be calculated exactly, the associated vison gap energy has to date only been calculated numerically from finite size diagonalization. Here we show that the phase shift for scattering Majorana fermions off a single bond-flip can be calculated analytically, leading to a closed-form expression for the vison gap energy Δ = 0.2633 J . Generalizations of our approach can be applied to Kitaev spin liquids on more complex lattices such as the three dimensional hyper-octagonal lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

First-principles calculations of phonons and Raman spectra in monoclinic CsSnCl 3

Density functional perturbation theory is used here to calculate the phonons at the zone center for monoclinic CsSnCl 3 . Furthermore, we report the calculated normal mode frequencies classified according to irreducible representations for both infrared-active and Raman-active modes. We also report the dielectric constants and Born-effective charges and investigate the anisotropy of the nonanalyticities near the zone center. The first-principles calculated Raman-active phonon frequencies and simulated Raman spectra are compared with experimental results from literature. We find that aside from the first-order allowed Raman modes of A g and B g symmetries, forbidden LO Raman modes appear prominently in the measured spectra. Several of these modes are close to or coincide with allowed modes and alter their intensity, while others appear as separate otherwise unexplained peaks. With this interpretation, good agreement is obtained with experimental mode frequencies to within a few cm -1 .

36 MATERIALS SCIENCE↗

Orbital-free extension to Kohn-Sham density functional theory equation of state calculations: Application to silicon dioxide

The liquid regime equation of state of silicon dioxide SiO 2 is calculated via quantum molecular dynamics in the density range of 5 to 15 g/cc and with temperatures from 0.5 to 100 eV, including the α-quartz and stishovite phase Hugoniot curves. Below 8 eV calculations are based on Kohn-Sham density functional theory (DFT), and above 8 eV a new orbital-free DFT formulation, presented here, based on matching Kohn-Sham DFT calculations is employed. Recent experimental shock data are found to be in very good agreement with the current results. Finally both experimental and simulation data are used in constructing a new liquid regime equation of state table for SiO 2 .

36 MATERIALS SCIENCE↗

Ab initiomany-body Green's function calculations of optical properties of LiF at high pressures

We present density functional theory (DFT) + quasiparticle self-energy (G 0 W 0 )+ Bethe-Salpeter calculations of the real and imaginary parts of the long-wavelength dielectric function of LiF between ambient pressure and P=5 Mbars. While the optical absorption spectrum is predicted to show dramatic pressure-dependent features above the optical gap, the index of refraction well below the gap is shown to exhibit the same trends as that seen in both DFT calculations and experiment: a roughly linear increase with density. This increase does not result from a decrease in the band gap, but rather follows from the increase in oscillator strength which counteracts a smaller increase in band gap with P. Our calculations also suggest that the index of refraction (for visible and near-UV light) of the higher-T B 2 phase should be quite close to that of the B 1 (ambient crystalline) phase. These findings may be of interest to researchers who use LiF as a window material in dynamic compression experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dielectric function and plasmons in graphene: A self-consistent-field calculation within a Markovian master equation formalism

We introduce a method for calculating the dielectric function of nanostructures with an arbitrary band dispersion and Bloch wave functions. The linear response of a dissipative electronic system to an external electromagnetic field is calculated by a self-consistent-field approach within a Markovian master equation formalism (SCF-MMEF) coupled with full-wave electromagnetic equations. The SCF-MMEF accurately accounts for several concurrent scattering mechanisms. The method captures interband electron-hole-pair generation, as well as the interband and intraband electron scattering with phonons and impurities. We employ the SCF-MMEF to calculate the dielectric function, complex conductivity, and loss function for supported graphene. From the loss-function maximum, we obtain plasmon dispersion and propagation length for different substrate types [nonpolar diamondlike carbon (DLC) and polar SiO 2 and hBN], impurity densities, carrier densities, and temperatures. Plasmons on the two polar substrates are suppressed below the highest surface phonon energy, while the spectrum is broad on the nonpolar DLC. Plasmon propagation lengths are comparable on polar and nonpolar substrates and are on the order of tens of nanometers, considerably shorter than previously reported. As a result, they improve with fewer impurities, at lower temperatures, and at higher carrier densities.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Anharmonic renormalization of the dispersion of flexural modes in graphene using atomistic calculations

In this work, we investigate through an atomistic method the effects of anharmonicity on the dispersion of flexural modes of graphene. Using a calculation based on ensemble averages of correlations among displacements and forces, we calculate the temperature-dependent frequencies for a semiempirical potential for graphene. We find that the dispersion relation of the flexural modes of graphene is renormalized by anharmonic coupling to other modes. Our calculations confirm that the anharmonic continuum results of Mariani and von Oppen [Phys. Rev. Lett. 100, 076801 (2008)] hold in detail for small wave number and at low temperatures. We examine the deviation from the continuum result outside of that range.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Atomistic calculations of dislocation core energy in aluminium

A robust molecular-dynamics simulation method for calculating dislocation core energies has been developed. This method has unique advantages: It does not require artificial boundary conditions, is applicable for mixed dislocations, and can yield converged results regardless of the atomistic system size. Utilizing a high-fidelity bond order potential, we have applied this method in aluminium to calculate the dislocation core energy as a function of the angle β between the dislocation line and the Burgers vector. These calculations show that, for the face-centered-cubic aluminium explored, the dislocation core energy follows the same functional dependence on β as the dislocation elastic energy: E c = Asin 2 β + Bcos 2 β, and this dependence is independent of temperature between 100 and 300 K. By further analyzing the energetics of an extended dislocation core, we elucidate the relationship between the core energy and the core radius of a perfect versus an extended dislocation. As a result, with our methodology, the dislocation core energy can accurately be accounted for in models of dislocation-mediated plasticity.

36 MATERIALS SCIENCE↗

Comment on “Comparative study of ab initio nonradiative recombination rate calculations under different formalisms”

In a recent article [Phys. Rev. B 91, 205315 (2015)] Shi, Xu, and Wang presented a comparison of several formalisms to calculate nonradiative recombination rates and concluded the “one-dimensional (1D) quantum formula” that was used by Alkauskas et al. [Phys. Rev. B 90, 075202 (2014)] is insufficient to accurately describe nonradiative capture rates. Our analysis of the results of Shi, Xu, and Wang indicates that their conclusions about the 1D quantum formula are unfounded and stem from an error in their calculations. Further, our own calculations demonstrate that the 1D quantum formula approach yields reliable and accurate results for nonradiative recombination rates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Using first-principles calculations to screen for fragile magnetism: Case study of LaCrGe 3 and LaCrSb 3

In this study, we present a coupled experimental/theoretical investigation of pressure effect on the ferromagnetism of LaCrGe 3 and LaCrSb 3 compounds. The magnetic, electronic, elastic, and mechanical properties of LaCrGe 3 and LaCrSb 3 at ambient condition are studied by first-principles density-functional theory calculations. The pressure dependences of the magnetic properties of LaCrGe 3 and LaCrSb 3 are also investigated. The ferromagnetism in LaCrGe 3 is rather fragile, with a ferro- to paramagnetic transition at a relatively small pressure (around 7 GPa from our calculations, and 2 GPa in experiments). The key parameter controlling the magnetic properties of LaCrGe 3 is found to be the proximity of the peak of Cr density of states to the Fermi level, a proximity that is strongly correlated with the distance between Cr atoms along the c axis, suggesting that there would be a simple way to suppress magnetism in systems with one-dimensional arrangement of magnetic atoms. By contrast, the ferromagnetism in LaCrSb 3 is not fragile. In conclusion, our theoretical results are consistent with our experimental results and demonstrate the feasibility of using first-principles calculations to aid experimental explorations in screening for materials with fragile magnetism.

36 MATERIALS SCIENCE↗

Ab initio no-core Gamow shell-model calculations of multineutron systems

The existence of multineutron systems has always been a debatable question. Indeed, both internucleon correlations and a large continuum coupling occur in these states. We then employ the ab initio no-core Gamow shell model to calculate the resonant energies and widths of the trineutron and tetraneutron systems with realistic interactions. In this work, our results indicate that trineutron and tetraneutron are both unbound and bear broad widths. The calculated energy and width of tetraneutron are also comparable with recent experimental data. Moreover, our calculations suggest that the energy of trineutron is lower than that of tetraneutron, while its resonance width is also narrower. This strongly suggests that trineutron is more likely to be experimentally observed than tetraneutron. We thus suggest experimentalists to search for trineutron at low energy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Benchmark calculations of pure neutron matter with realistic nucleon-nucleon interactions

Here, we report benchmark calculations of the energy per particle of pure neutron matter as a function of the baryon density using three independent many-body methods: Brueckner–Bethe–Goldstone, Fermi hypernetted chain/single-operator chain, and auxiliary-field diffusion Monte Carlo. Significant technical improvements are implemented in the latter two methods. The calculations are made for two distinct families of realistic coordinate-space nucleon-nucleon potentials fit to scattering data, including the standard Argonne v 18 interaction and two of its simplified versions, and four of the new Norfolk Δ -full chiral effective field theory potentials. Primarily because of the advancements in the auxiliary-field diffusion Monte Carlo, we observe good agreement among the three many-body techniques up to nuclear saturation density—the maximum difference in the energy per particle is within 1.5 MeV for all the potentials we consider. At higher densities, the divergences become more important, and are mainly driven by the Fermi hypernetted chain/single-operator calculations. We also study the connection between nucleon-nucleon scattering data and the energy per particle of pure neutron matter. Our results suggest that fitting to higher-energy nucleon-nucleon scattering helps reduce the spread of energies among the models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Calculation of the C 12 + C 12 sub-barrier fusion cross section in an imaginary-time-dependent mean field theory

The 12 C + 12 C sub-barrier fusion cross section is calculated within the framework of a time-dependent Hartree-Fock-based classical model using the Feynman path-integral method. The modified astrophysical S* factor is compared to direct and indirect experimental results. A good agreement with the direct data is found. In the lower-energy region where recent analyses of experimental data obtained with the Trojan horse method (THM) lead to contrasting results, the model predicts a nonresonant S* factor half-way between those results. Low-energy resonances revealed in the THM data are added to the calculation, and the relative reaction rate in the Gamow region is calculated. In particular, including 0 + resonances result in some agreement with the THM data. Here, the role of different resonances is discussed in detail, and their influence on the reaction rate at temperatures relevant to stellar evolution is investigated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantum many-body calculations using body-centered cubic lattices

It is often computationally advantageous to model space as a discrete set of points forming a lattice grid. This technique is particularly useful for computationally difficult problems such as quantum many-body systems. For reasons of simplicity and familiarity, nearly all quantum many-body calculations have been performed on simple cubic lattices. Since the removal of lattice artifacts is often an important concern, it would be useful to perform calculations using more than one lattice geometry. In this paper we show how to perform quantum many-body calculations using auxiliary-field Monte Carlo simulations on a three-dimensional body-centered cubic (BCC) lattice. As a benchmark test we compute the ground state energy of 33 spin-up and 33 spin-down neutrons in the unitary limit, which is an idealized limit where the interaction range is zero and scattering length is infinite. As a fraction of the free Fermi gas energy E FG , we find that the ground state energy is E 0 /E FG =0.369(2),0.371(2), using two different definitions of the finite-system energy ratio. This is in excellent agreement with recent results obtained on a cubic lattice [He et al., Phys. Rev. A 101, 063615 (2020)]. We find that the computational effort and performance on a BCC lattice is approximately the same as that for a cubic lattice with the same number of lattice points. We discuss how the lattice simulations with different geometries can be used to constrain the size of lattice artifacts in simulations of continuum quantum many-body systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ab initio calculation of the β decay from 11 Be to a 10 Be + $\textit{p}$ resonance

In this work, the exotic β -delayed proton emission is calculated in 11 Be from first principles using chiral two- and three-nucleon forces. To investigate the unexpectedly large branching ratio measured in Ayyad et al. [Phys. Rev. Lett. 123, 082501 (2019)], we calculate the proposed (1/2 + , 1/2) proton resonance in 11 B using the no-core shell model with continuum. This calculation helps to address whether this enhancement is caused by unknown dark decay modes or an unobserved proton resonance. We report a branching ratio of $b_p$ = (1.3 ± 0.5) × 10 –6 , suggesting that its unexpectedly large value is caused by an unobserved proton resonance in 11 B.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗