Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles 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 1,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

First-principles study of the surface properties of uranium carbides

Uranium carbides have attracted renewed interest as advanced nuclear fuels for Generation IV reactors. As an important property required for gas bubble modeling in nuclear fuels, the surface energy of uranium carbides is scarce in literature. In this work, we study the surface properties of uranium carbides by first-principles density functional theory calculations. Surface orientations with maximum Miller index up to 3, 2 and 2 are investigated for UC, U 2 C 3 and α-UC 2 , respectively. By studying the effects of surface termination and chemical potential on surface energy, we identify the factors that determines the surface stability. From the calculated surface energies, the surface properties of uranium carbide single crystals are obtained from Wulff construction, including equilibrium morphology, dominant surface orientation and area weighted surface energy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Pressure-induced creation and annihilation of Weyl points in T d - Mo 0.5 W 0.5 Te 2 and 1 T " - Mo 0.5 W 0.5 Te 2

By means of first-principles density-functional theory calculations, we investigate the role of hydrostatic pressure in the electronic structure of the T d (Pmn2 1 ) and 1T" (Pm) phases of the Weyl semimetal Mo 0.5 W 0.5 Te 2 , which is a promising material for phase-change memory technology and superconductivity. We particularly focus on changes occurring in the distribution of the gapless Weyl points (WPs) within the 0 to 45 GPa pressure range. Here we further investigate the structural phase transition and lattice dynamics of the T d and 1T" phases within the aforementioned pressure range. Our calculations suggest that both the T d and 1T" phases of Mo 0.5 W 0.5 Te 2 host four WPs in their full Brillouin zone at zero pressure. The total number of WPs increases to 44 (36) with increasing pressure via pair creation up to 20 (15) GPa for the T d (1T") phase, and beyond this pressure pair annihilation of WPs starts occurring, leaving only 16 WPs at 45 GPa in both phases. The enthalpy versus pressure data reveal that the 1T" phase is more favorable below the critical pressure of 7.5 GPa; however, beyond this critical pressure the T d phase becomes enthalpically favorable. We also provide the calculated x-ray diffraction spectra along with the calculated Raman- and infrared-active phonon frequencies to facilitate the experimental identification of the studied phases

36 MATERIALS SCIENCE↗
DOE OSTI2023

Gate-Tunable Phonon Magnetic Moment in Bilayer Graphene

Here we develop a first-principles quantum scheme to calculate the phonon magnetic moment in solids. As a showcase example, we apply our method to study gated bilayer graphene, a material with strong covalent bonds. According to the classical theory based on the Born effective charge, the phonon magnetic moment in this system should vanish, yet our quantum mechanical calculations find significant phonon magnetic moments. Furthermore, the magnetic moment is highly tunable by changing the gate voltage. Our results firmly establish the necessity of the quantum mechanical treatment, and identify small-gap covalent materials as a promising platform for studying tunable phonon magnetic moment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Suppressed Fluctuations as the Origin of the Static Magnetic Order in Strained Sr 2 RuO 4

Combining first-principles density-functional calculations and Moriya's self-consistent renormalization theory, we explain the recently reported counterintuitive appearance of an ordered magnetic state in uniaxially strained Sr 2 RuO 4 beyond the Lifshitz transition. We show that strain weakens the quantum spin fluctuations, which destroy the static order, more strongly than the tendency to magnetism. Furthermore, a different rate of decrease of the spin fluctuations vs magnetic stabilization energy promotes the onset of a static magnetic order beyond a critical strain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

S 8 ring content in liquid sulfur determined from x-ray diffraction

The S 8 ring content in liquid sulfur has been determined from the 4.47 Å peak in the x-ray pair distribution function over the temperature range −55 °C –330 °C. This peak represents the characteristic 3rd and 4th nearest atom-atom distances across an S 8 -ring. Three procedures were used to calculate the S 8 -ring fraction; scaling to crystalline β-S, a linear combination fit and a direct integration method. All methods gave similar values, indicating that approximately one third of the S 8 -rings break up upon melting into a distribution of short chains with an average length of 5 ± 2 atoms. Our results contradict established models of low temperature sulfur that indicate below the λ-transition the liquid comprises almost entirely of S8-rings. In addition, we find the S 8 -ring fraction increases to 76 ± 10% at −30 °C in the supercooled liquid, and decreases to 23% ± 3% at 330 °C in the high temperature melt. First principle molecular dynamics calculations confirm the onset of dynamic ring distortion via broad bond angle distributions, associated with the interstitial region between the first and second peaks in the pair distribution function. Our results are consistent with the presence of a propensity of short chains in the low temperature liquid, that are required for the bond switching mechanism responsible for the polymerization process to occur at the λ-transition temperature, T λ .

lambda transition↗
DOE OSTI2023

Thermodynamic modeling of the Nb-Ni system with uncertainty quantification using PyCalphad and ESPEI

Here, the Nb–Ni system is remodeled with uncertainty quantification (UQ) using software tools of PyCalphad and ESPEI (the Extensible, Self-optimizing Phase Equilibria Infrastructure) with the presently implemented capability of modeling site fraction based on Wyckoff positions. The five- and three-sublattice models are used to model the topologically close pack (TCP) μ-Nb 7 Ni 6 and δ-NbNi 3 phases according to their Wyckoff positions. The inputs for CALPHAD-based thermodynamic modeling include the thermochemical data as a function of temperature predicted by first-principles and phonon calculations based on density functional theory (DFT), ab initio molecular dynamics (AIMD) simulations, together with phase equilibrium and site fraction data in the literature. In addition to phase diagram and thermodynamic properties, the CALPHAD-based predictions of site fractions of Nb in μ-Nb 7 Ni 6 agree well with experimental data. Furthermore, the UQ estimation using the Markov Chain Monte Carlo (MCMC) method as implemented in ESPEI is applied to study the uncertainty of site fraction in μ-Nb 7 Ni 6 and enthalpy of mixing (ΔH mix ) in liquid.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Out-of-plane polarization and topological magnetic vortices in multiferroic CrPSe 3

Two-dimensional multiferroic materials are ideal systems for exploring new coupling mechanisms between different ferroic orders and producing novel quantum phenomena with potential applications. Here we employ first-principles density functional theory calculations to discover intrinsic ferroelectric and antiferroelectric phases of Cr⁢P⁢Se 3 , which show ferromagnetic order and compete with the centrosymmetric phase with an antiferromagnetic order. Our analysis shows that the electrical dipoles of such type-I multiferroic phases come from the out-of-plane displacements of phosphorus ions due to the stereochemically active lone pairs. The coupling between polar and magnetic orders creates the opportunity for tuning the magnetic ground state by switching from the centrosymmetric to the ferroelectric phase using an out-of-plane electric field. In ferroelectric and antiferroelectric phases, the combination of easy-plane anisotropy and Dzyaloshinskii-Moriya interactions indicates that they can host topological magnetic vortices such as meron pairs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electronic structure and lattice dynamics of 1⁢𝑇−VSe 2 : Origin of the three-dimensional charge density wave

To characterize in detail the charge density wave (CDW) transition of 1⁢𝑇−VSe 2 , its electronic structure and lattice dynamics are comprehensively studied by means of x-ray diffraction, muon spectroscopy, angle resolved photoemission (ARPES), diffuse and inelastic x-ray scattering, and state-of-the-art first-principles density functional theory calculations. Resonant elastic x-ray scattering does not show any resonant enhancement at either V or Se, indicating that the CDW peak at the 𝐾 edges describes a purely structural modulation of the electronic ordering. ARPES experiments identify (i) a pseudogap at 𝑇 > 𝑇⁢ CDW , which leads to a depletion of the density of states in the ML-M'L' plane at 𝑇 < 𝑇⁢ CDW , and (ii) anomalies in the electronic dispersion reflecting a sizable impact of phonons on it. A diffuse scattering precursor, characteristic of soft phonons, is observed at room temperature (RT) and leads to the full collapse of the low-energy phonon (𝜔 1 ) with propagation vector (0.25 0 −0.3) r.l.u. Here, we show that the frequency and linewidth of this mode are anisotropic in momentum space, reflecting the momentum dependence of the electron-phonon interaction (EPI), hence demonstrating that the origin of the CDW is, to a much larger extent, due to the momentum dependent EPI with a small contribution from nesting. The pressure dependence of the 𝜔 1 soft mode remains nearly constant up to 13 GPa at RT, with only a modest softening before the transition to the high-pressure monoclinic C2/m phase. The wide set of experimental data is well captured by our state-of-the art first-principles anharmonic calculations with the inclusion of van der Waals corrections in the exchange-correlation functional. The comprehensive description of the electronic and dynamical properties of VSe 2 reported here adds important pieces of information to the understanding of the electronic modulations in the family of transition-metal dichalcogenides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Thermal transport in Yb-based 1-2-20 materials

Abstract Lattice thermal-transport properties of heavy-fermion YbT 2 Zn 20 -based (T = Co, Rh, and Ir) 1-2-20 compounds are calculated with the single-mode relaxation-time approximation and the full solution of the linearized phonon Boltzmann transport equation from first-principles anharmonic phonon calculations. We predict low lattice thermal conductivity κ L with a maximum value of ∼ 24.16 · K at 85.14 K, 137.29 · K at 9.67 K, and 23.55 · K at 67.50 K and ∼ 12.02, 40.04, and 10.30 · K at room temperature for YbCo 2 Zn 20 , YbRh 2 Zn 20 , and YbIr 2 Zn 20 , respectively. Based on the analysis of the cumulative κ L as a function of the phonon mean free path and the frequency-dependent joint density of states at various temperatures, we attribute the low κ L to a rattling mode with an average characteristic rattling frequency that enhanced phonon scattering processes, which reduced the phonon mean free paths, suppressed the phonon lifetime, and enhanced the probability of three-phonon scattering events. The predicted low κ L , especially in YbCo 2 Zn 20 and YbIr 2 Zn 20 , makes them promising candidate materials for thermoelectric applications and thermal management.

Physics↗
DOE OSTI2025

Gluon unpolarized, polarized, and transversity GPDs from lattice QCD: Lorentz-covariant parametrization

We identify the matrix elements necessary to determine the leading-twist gluon generalized parton distributions (GPDs) H g , E g , H ˜ g , E ˜ g , H g T , E g T , H ˜ g T , E ˜ g T in lattice QCD calculations. We present a method to achieve a Lorentz-covariant parameterization of the matrix elements in terms of a linearly independent basis of tensor structures. This parameterization is crucial for projecting lattice QCD matrix elements onto light cone distributions. For the first time, we determine the corresponding components that project onto the linear combinations of invariant amplitudes, which reduce to the different gluon GPDs in the light cone limit and enable their separation in a lattice QCD calculation for spin-0 and spin- 1 2 hadrons. Hence, this work lays the foundation for the numerical determination of the gluon GPDs from first-principle lattice QCD calculations, directly advancing our understanding of the mass and spin structures and mechanical properties of the nucleon, as well as the physics underlying deeply virtual Compton scattering and deeply virtual meson production in a range of experimental processes. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

The Local Electronic Structure of Supercritical CO 2 from X-ray Raman Spectroscopy and Atomistic-Scale Modeling

Supercritical CO 2 is encountered in several technical and natural systems related to biology, geophysics, and engineering. While the structure of gaseous CO 2 has been studied extensively, the properties of supercritical CO 2 , particularly close to the critical point, are not well-known. Here, in this work, we combine X-ray Raman spectroscopy, molecular dynamics simulations, and first-principles density functional theory (DFT) calculations to characterize the local electronic structure of supercritical CO 2 at conditions around the critical point. The X-ray Raman oxygen K-edge spectra manifest systematic trends associated with the phase change of CO 2 and the intermolecular distance. Extensive first-principles DFT calculations rationalize these observations on the basis of the 4sσ Rydberg state hybridization. X-ray Raman spectroscopy is found to be a sensitive tool for characterizing electronic properties of CO 2 under challenging experimental conditions and is demonstrated to be a unique probe for studying the electronic structure of supercritical fluids.

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

A comprehensive first-principles study of the effects of the exchange-correlation functional and magnetism on defect and diffusion properties of the CoCrNi medium-entropy alloy

The present work is a novel, systematic study of the effect of density functional theory input parameters on the vacancy formation energy (VFE), migration barrier for diffusion, and electronic structure for each element in the CoCrNi medium-entropy alloy (MEA). In particular, the novelties include: (1) calculating the aforementioned properties of Co, Cr, or Ni, in the CoCrNi MEA using magnetic and non-magnetic states, and two versions of the generalized gradient approximation: Perdew, Burke, and Ernzerhof (PBE) and the PBE version for solids (PBEsol), and (2) a detailed comparison of 0 K activation energy to experimental creep activation energies. First-principles calculations at 0 K are performed using the Vienna ab-initio simulation package. Special quasirandom structures (SQS) and Widom-type substitution are employed. For each element, Co, Cr, or Ni, non-magnetic calculations result in a higher VFE and larger range of calculated values for the configurations studied. The averaged migration barrier is the highest for Co in the CoCrNi for three of four sets of calculation parameters in the configurations studied. Finally, the results indicate that the average 0 K activation energy for diffusion makes up 70–80% of the experimental creep activation energy, depending on the exchange-correlation functional employed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Radiative capture rates at deep defects from electronic structure calculations

We present a methodology to calculate radiative carrier capture coefficients at deep defects in semiconductors and insulators from first principles. Electronic structure and lattice relaxations are accurately described with hybrid density functional theory. Calculations of capture coefficients provide an additional validation of the accuracy of these functionals in dealing with localized defect states. We also discuss the validity of the Condon approximation, showing that even in the event of large lattice relaxations the approximation is accurate. We test the method on GaAs:VGa-TeAs and GaN:CN, for which reliable experiments are available, and demonstrate very good agreement with measured capture coefficients.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Vacancy-driven variations in the phonon density of states of fast neutron irradiated nuclear graphite

Research examining radiation damage and its effects in graphite began in the 1940’s during the development of moderated nuclear reactors. Interest in this topic is expanding because of emerging applications associated with fullerenes and carbon nanostructures as well as its long-standing use as a fission moderator. In this work, we report the measurements of the full phonon density of states of irradiated nuclear graphite using inelastic neutron scattering experiments at room temperature. The samples were previously exposed to different levels of neutron doses and irradiation temperatures. The phonon density of states of perfect and defected graphite supercells including different configurations of vacancies and interstitials were calculated using first-principles direct method. Despite high neutron irradiation doses, the scattering measurements show that the induced damage is localized, and the layered structure is preserved. In conclusion, a comparison of the measured phonon densities of states of irradiated samples with those calculated of defected supercells indicates that the main changes observed are mainly attributed to formation of vacancies.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Exciton Enhanced Nonlinear Optical Responses in Monolayer h-BN and MoS 2 : Insight from First-Principles Exciton-State Coupling Formalism and Calculations

Excitons are vital in the photophysics of materials, especially in low-dimensional systems. The conceptual and quantitative understanding of excitonic effects in nonlinear optical (NLO) processes is more challenging compared to linear ones. Here, we present an ab initio approach to second-order NLO responses, incorporating excitonic effects, that employs an exciton-state coupling formalism and allows for a detailed analysis of the role of individual excitonic states. Taking monolayer h-BN and MoS 2 as two prototype 2D materials, we calculate their second harmonic generation (SHG) susceptibility and shift current conductivity tensor. We find strong excitonic enhancement in the NLO responses requires that the resonant excitons are not only optically bright themselves but also able to couple strongly to other bright excitons. Our results explain the occurrence of two strong peaks in the SHG of monolayer h-BN and why the A and B excitons of MoS 2 unexpectedly exhibit minimal excitonic enhancement in both SHG and shift current generation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Quasi-one-dimensional uniform spin- 1 2 Heisenberg antiferromagnet KNaCuP 2 O 7 probed by P 31 and Na 23 NMR

Here we present the structural and magnetic properties of KNaCuP 2 O 7 investigated via x-ray diffraction, magnetization, specific heat, and P 31 and Na 23 NMR measurements and complementary electronic structure calculations. The temperature-dependent magnetic susceptibility and P 31 NMR shift could be modeled very well by the uniform spin- 1 2 Heisenberg antiferromagnetic chain model with a nearest-neighbor interaction J / k B ≃ 58.7 K . The corresponding mapping using first-principles electronic structure calculations leads to J DFT / k B ≃ 59 K with negligibly small interchain couplings, further confirming that the system is indeed a one-dimensional uniform spin- 1 2 Heisenberg antiferromagnet. The diverging trend of NMR spin-lattice relaxation rates ( 1 / T 1 31 and 1 / T 1 23 ) implies the onset of a magnetic long-range ordering at around T N ≃ 1 K . From the value of T N , the average interchain coupling is estimated to be J ' / k B ≃ 0.28 K . Moreover, the NMR spin-lattice relaxation rates show the dominant contributions from uniform ( q = 0 ) and staggered ( q = ± π / a ) spin fluctuations in the high- and low-temperature regimes, respectively, mimicking one-dimensionality of the spin lattice. We have also demonstrated that 1 / T 1 31 in high temperatures varies linearly with 1 / H , reflecting the effect of spin diffusion on the dynamic susceptibility. The temperature-dependent unit cell volume could be described well using the Debye approximation with a Debye temperature of Θ D ≃ 294 K , consistent with the heat capacity data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Measurement of the Sixth-Order Cumulant of Net-Proton Multiplicity Distributions in Au + Au Collisions at s NN = 27 , 54.4, and 200 GeV at RHIC

According to first principle Lattice QCD calculations, the transition from quark-gluon plasma to hadronic matter is a smooth crossover in the region µB ≤ Tc. As a result, higher-order cumulants and their ratios are predicted to be negative, C 6 /C 2 < 0, for example. In this paper, we report the first measurement of the midrapidity net-proton C 6 /C 2 from 27, 54.4 and 200 GeV Au+Au collisions at RHIC. The dependence on collision centrality and kinematic acceptance in (p T, y ) are analyzed. While for 27 and 54.4 GeV collisions the C 6 /C 2 values are close to zero within uncertainties, it is observed that for 200 GeV collisions, the C 6 /C 2 ratio becomes progressively negative from peripheral to central collisions. Transport model calculations without critical dynamics predict values around zero. These observations seem to favor a smooth crossover in the high energy nuclear collisions at RHIC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Clathrate BaNi 2 P 4 : An Interplay of Heat and Charge Transport Due to Strong Host–Guest Interactions

Heat and charge transport properties of the metallic unconventional clathrate BaNi 2 P 4 , hosting Ba cations in oversized Ni 8 P 16 cages, are investigated. A novel method of single-crystal growth was developed, yielding 2–3 mm sized crystals of BaNi 2 P 4 . We also developed a setup to accurately measure thermal conductivity and electrical resistivity of the synthesized single crystals in a wide temperature range avoiding crystal remounting. BaNi 2 P 4 has a metallic temperature dependence of its electrical resistivity (decreasing with decreasing temperature) and manifests an unconventional T 2 power law for 50 K < T < 300 K; below 50 K, the power-law exponent increases gradually such that below 10 K the power law is T 5 , a predicted but extremely rarely experimentally observed dependence for peculiar electron–phonon interactions. Furthermore, electronic band structure calculations, consistent with measurements of de Haas–van Alphen oscillations, show large band dispersions with significant contributions of Ba orbitals to states near the Fermi level, which is atypical for clathrates. The thermal properties of BaNi 2 P 4 were probed using a combination of variable-temperature single-crystal X-ray diffraction experiments, heat capacity measurements, first-principles phonon dispersion calculations, and inelastic neutron scattering measurements. BaNi 2 P 4 exhibits significant hybridization of the Ba-guest and Ni–P-framework vibrational modes, which may be enhanced via the detected split of the Ba position, which results in strong Ba-framework interactions.

36 MATERIALS SCIENCE↗