Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles”

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 325 records · Page 18

Competing stripe and magnetic phases in the cuprates from first principles

Realistic description of competing phases in complex quantum materials has proven extremely challenging. For example, much of the existing density-functional-theory-based first-principles framework fails in the cuprate superconductors. Various many-body approaches involve generic model Hamiltonians and do not account for the interplay between the spin, charge, and lattice degrees of freedom. Here, by deploying the recently constructed strongly constrained and appropriately normed (SCAN) density functional, we show how the landscape of competing stripe and magnetic phases can be addressed on a first-principles basis both in the parent insulator YBa 2 Cu 3 O 6 and the near-optimally doped YBa 2 Cu 3 O 7 as archetype cuprate compounds. In YBa 2 Cu 3 O 7 , we find many stripe phases that are nearly degenerate with the ground state and may give rise to the pseudogap state from which the high-temperature superconducting state emerges. We invoke no free parameters such as the Hubbard U , which has been the basis of much of the existing cuprate literature. Lattice degrees of freedom are found to be crucially important in stabilizing the various phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Strain field of the monovacancy in silicene: First-principles study

The in-plane strain fields of single-vacancy silicene with different monovacancy (MV) concentrations, as well as the corresponding electronic band structures, are investigated by using the first-principle calculations. Firstly the self-healing MV is found to be the most stable ground structure in silicene, which is different from the other 2D hexagonal honeycomb materials, e.g. graphene, h-BN. In the isolated MV center, the bonds along the pentagons are compressed, creating a compress field, and those close to the distorted hexagons are stretched, creating a stretch field. As the MV concentration increasing, the interacted compress field tends to corrugate the defected silicene, while the interacted stretch field impacts little on the low-buckled structure. Especially, the corrugation presents in those supercells with small MV concentration, just as the (4, 5), (4, 6), (4, 7), (4, 8) supercells. The corrugations approach zero at both low and high MV concentrations, and the (4, 6) supercell with a MV concentration of about 0.021, has a peak value of 3.23Å. The electronic calculations show that the linear dispersion at Γ point in pristine silicene is broken by the lower lattice symmetry of the self-healing MV reconstruction, which translates it into metal as well.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

First-principles studies on the equation-of-state, thermal-conductivity, and opacity of deuterium-tritium and polystyrene (CH) for inertial confinement fusion applications

Using first-principles (FP) methods, we have performed ab initio compute for the equation of state (EOS), thermal conductivity, and opacity of deuterium-tritium (DT) in a wide range of densities and temperatures for inertial confinement fusion (ICF) applications. These systematic investigations have recently been expanded to accurately compute the plasma properties of CH ablators under extreme conditions. In particular, the first-principles EOS and thermal-conductivity tables of CH are self-consistently built from such FP calculations, which are benchmarked by experimental measurements. When compared with the traditional models used for these plasma properties in hydrocodes, significant differences have been identified in the warm dense plasma regime. When these FP-calculated properties of DT and CH were used in our hydrodynamic simulations of ICF implosions, we found that the target performance in terms of neutron yield and energy gain can vary by a factor of 2 to 3, relative to traditional model simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Towards Accurate and Efficient Predictions of Martensitic Transition Temperatures for Shape Memory Alloys from First Principles

Recent rapid progresses in physics theory and computational power have made it possible to predict the martensitic transition temperatures (MTTs) in shape memory alloys (SMAs) from first principles [1-3]. In particular, rigorous while time-consuming thermodynamic integration has been employed to compute the anharmonic phonon free energies, which play a crucial role in determining martensitic phase transitions in SMAs. However, this approach has only been applied to simple binaries, and its accuracy is unsatisfying for certain SMAs such as the most commonly used NiTi. In this work, we report on several new developments to our method that bring first-principles theory and experiment much closer into agreement including the MTT of NiTi, and that improve the computational efficiency significantly. We have applied our refined approach to investigate the Ni0.5Ti0.5-xHfx and PdxNi0.5-xTi0.5 ternaries, and the predicted MTT for each composition is within 100K compared with experiment. We will address various techniques to overcome the difficulty encountered in studying ternaries. Our theoretical approach is expected to be a broadly applicable and predictive theory for designing complex SMAs with desirable properties. [1] J.B. Haskins, A.E. Thompson,and J.W. Lawson, Phys. Rev B 94, 214110 (2016). [2] J.B. HaskinsandJ.W. Lawson, J. App. Phys. 121, 205103 (2017). [3] J.B. Haskins, H. Malmir, S. J. Honrao, L. A. Sandoval, and J.W. Lawson, Acta Materialia 212, 116872 (2017).

Zhigang Wu↗

Towards Accurate and Efficient Predictions of Martensitic Transition Temperatures for Shape Memory Alloys from First Principles

Recent rapid progresses in physics theory and computational power have made it possible to predict the martensitic transition temperatures (MTTs) in shape memory alloys (SMAs) from first principles [1-3]. In particular, rigorous while time-consuming thermodynamic integration has been employed to compute the anharmonic phonon free energies, which play a crucial role in determining martensitic phase transitions in SMAs. However, this approach has only been applied to simple binaries, and its accuracy is unsatisfying for certain SMAs such as the most commonly used NiTi. In this work, we report on several new developments to our method that bring first-principles theory and experiment much closer into agreement including the MTT of NiTi, and that improve the computational efficiency significantly. We have applied our refined approach to investigate the Ni0.5Ti0.5-xHfx and PdxNi0.5-xTi0.5 ternaries, and the predicted MTT for each composition is within 100K compared with experiment. We will address various techniques to overcome the difficulty encountered in studying ternaries. Our theoretical approach is expected to be a broadly applicable and predictive theory for designing complex SMAs with desirable properties.

Zhigang Wu↗

Photoluminescence spectra of point defects in semiconductors: Validation of first-principles calculations

Optically and magnetically active point defects in semiconductors are interesting platforms for the development of solid state quantum technologies. Their optical properties are usually probed by measuring photoluminescence spectra, which provide information on excitation energies and on the interaction of electrons with lattice vibrations. We present a combined computational and experimental study of photoluminescence spectra of defects in diamond and SiC, aimed at assessing the validity of theoretical and numerical approximations used in first-principles calculations, including the use of the Franck-Condon principle and the displaced harmonic oscillator approximation. We focus on prototypical examples of solid state qubits, the divacancy centers in SiC and the nitrogen-vacancy in diamond, and we report computed photoluminescence spectra as a function of temperature that are in very good agreement with the measured ones. As expected we find that the use of hybrid functionals leads to more accurate results than semilocal functionals. Interestingly our calculations show that constrained density functional theory (CDFT) and time-dependent hybrid DFT perform equally well in describing the excited state potential energy surface of triplet states; our findings indicate that CDFT, a relatively cheap computational approach, is sufficiently accurate for the calculations of photoluminescence spectra of the defects studied here. Finally, we find that only by correcting for finite-size effects and extrapolating to the dilute limit can one obtain a good agreement between theory and experiment. Our results provide a detailed validation protocol of first-principles calculations of photoluminescence spectra, necessary both for the interpretation of experiments and for robust predictions of the electronic properties of point defects in semiconductors.

36 MATERIALS SCIENCE↗

Magnetic ordering in TbMn{sub 0.5}Cr{sub 0.5}O{sub 3} studied by neutron diffraction and first-principles calculations

The structure and magnetic ordering of bulk TbMn{sub 0.5}Cr{sub 0.5}O{sub 3} was revealed through bulk magnetization and neutron diffraction measurements, and first-principles calculations, respectively. G-type antiferromagnetic ordering of Mn{sup 3+} and Cr{sup 3+} moments was observed in the neutron diffraction data below Néel temperature T{sub N} ∼ 84 K. In addition, below ∼40 K, short-range magnetic ordering was identified correlating to a ferromagnetic component due to the canting of the moments along the c-axis. The spin configuration is consistent with the first-principles calculations. The magnetic structure revealed in the present TbMn{sub 0.5}Cr{sub 0.5}O{sub 3} sample is distinct from that observed for both end members TbMnO{sub 3} and TbCrO{sub 3}.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Martensitic Phase Transition in Complex NiTi-Based Shape Memory Alloys from First-Principles Calculations

Recent rapid progresses in physics theory and computational power have made it possible to accurately predict the phase transitions and martensitic transition temperatures (MTTs) in shape memory alloys (SMAs) from first principles. However, previously theory and calculations [1-4] were applied only to study highly ordered stoichiometric binary alloys such as NiTi, PdTi and NiHf. Here we report on our recent first-principles investigations [5,6] on Ni 0.5 Ti 0.5-x Hf x and Pd x Ni 0.5-x Ti 0.5 ternaries and off-stoichiometric NiTi, and the predicted martensitic phase transitions in these complex SMAs are in good agreement with experimental findings. In particular, the calculated MTTs for all these compositions are within 100 K compared with the corresponding measured data, and our results also reveal the origin of the striking asymmetry in MTT of the off-stoichiometric NiTi near equiatomic compositions. We will address various techniques to overcome the difficulty encountered in studying ternaries and off-stoichiometic binaries associated with disorder and/or much lowered symmetry. Our theoretical approach is expected to be a broadly applicable and predictive theory for designing complex SMAs with desirable properties. References: [1] J. B. Haskins, A. E. Thompson, and J. W. Lawson, Phys. Rev B 94 , 214110 (2016). [2] J. B. Haskins and J. W. Lawson, J. App. Phys. 121 , 205103 (2017). [3] J. B. Haskins, H. Malmir, S. J. Honrao, L. A. Sandoval, and J. W. Lawson, Acta Materialia 212 , 116872 (2017). [4] Z. Wu, J. W. Lawson, and O. Benafan, Phys. Rev. B 106 , L140102 (2022). [5] Z. Wu, H. Malmir, O. Benafan, and J. W. Lawson, Acta Materialia 261 , 119362 (2023). [6] Z. Wu, J. W. Lawson, and O. Benafan, Phys. Rev. B 108 , L140103 (2023).

Zhigang Wu↗

Trajectory sampling and finite-size effects in first-principles stopping power calculations

Abstract Real-time time-dependent density functional theory (TDDFT) is presently the most accurate available method for computing electronic stopping powers from first principles. However, obtaining application-relevant results often involves either costly averages over multiple calculations or ad hoc selection of a representative ion trajectory. We consider a broadly applicable, quantitative metric for evaluating and optimizing trajectories in this context. This methodology enables rigorous analysis of the failure modes of various common trajectory choices in crystalline materials. Although randomly selecting trajectories is common practice in stopping power calculations in solids, we show that nearly 30% of random trajectories in an FCC aluminum crystal will not representatively sample the material over the time and length scales feasibly simulated with TDDFT, and unrepresentative choices incur errors of up to 60%. We also show that finite-size effects depend on ion trajectory via “ouroboros” effects beyond the prevailing plasmon-based interpretation, and we propose a cost-reducing scheme to obtain converged results even when expensive core-electron contributions preclude large supercells. This work helps to mitigate poorly controlled approximations in first-principles stopping power calculations, allowing 1–2 order of magnitude cost reductions for obtaining representatively averaged and converged results.

36 MATERIALS SCIENCE↗

First-principles prediction of the softening of the silicon shock Hugoniot curve

Here, whock compression of silicon (Si) under extremely high pressures (>100 Mbar) was investigated by using two first-principles methods of orbital-free molecular dynamics (OFMD) and path integral Monte Carlo (PIMC). While pressures from the two methods agree very well, PIMC predicts a second compression maximum because of 1s electron ionization that is absent in OFMD calculations since Thomas–Fermi-based theories lack inner shell structure. The Kohn–Sham density functional theory is used to calculate the equation of state (EOS) of warm dense silicon for low-pressure loadings (P < 100 Mbar). Combining these first-principles EOS results, the principal Hugoniot curve of silicon for pressures varying from 0.80 Mbar to above ~10 Gbar was derived. We find that silicon is ~20% or more softer than what was predicted by EOS models based on the chemical picture of matter. Existing experimental data (P ≈ 1–2 Mbar) seem to indicate this softening behavior of Si, which calls for future strong-shock experiments (P > 10 Mbar) to benchmark our results.

36 MATERIALS SCIENCE↗

Validating first-principles phonon lifetimes via inelastic neutron scattering

Phonon lifetimes are a key component of quasiparticle theories of transport; yet first-principles lifetimes are rarely directly compared with inelastic neutron scattering (INS) results. Existing comparisons show discrepancies even at temperatures where perturbation theory is expected to be reliable. In this paper, we demonstrate that the reciprocal space voxel (q voxel), which is the finite region in reciprocal space required in INS data analysis, must be explicitly accounted for within theory in order to draw a meaningful comparison. Here, we demonstrate accurate predictions of peak widths of the scattering function when accounting for the q voxel in CaF 2 and ThO 2 . Passing this test implies high fidelity of the phonon interactions and the approximations used to compute the Green's function, serving as a critical benchmark of theory and indicating that other material properties should be accurately predicted, which we demonstrate for thermal conductivity.

36 MATERIALS SCIENCE↗

First-principles theory of the Dirac semimetal Cd 3 As 2 under Zeeman magnetic field

Time-reversal broken Weyl semimetals have attracted much attention recently, but certain aspects of their behavior, including the evolution of their Fermi surface topology and anomalous Hall conductivity with Fermi-level position, have remained underexplored. A promising route to obtain such materials may be to start with a nonmagnetic Dirac semimetal and break time-reversal symmetry via magnetic doping or magnetic proximity. Here we explore this scenario in the case of the Dirac semimetal Cd 3 As 2 , based on first-principles density-functional calculations and subsequent low-energy modeling of Cd 3 As 2 in the presence of a Zeeman field applied along the symmetry axis. We clarify how each four-fold degenerate Dirac node splits into four Weyl nodes, two with chirality ±1 and two higher-order nodes with chirality ±2. Using a minimal k · p model Hamiltonian whose parameters are fit to the first-principles calculations, we detail the evolution of the Fermi surfaces and their Chern numbers as the Fermi energy is scanned across the region of the Weyl nodes at fixed Zeeman field. We also compute the intrinsic anomalous Hall conductivity as a function of Fermi-level position, finding a characteristic inverted-dome structure. Cd 3 As 2 is especially well suited to such a study because of its high mobility, but the qualitative behavior revealed here should be applicable to other Dirac semimetals as well.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electron-Phonon Scattering in the Presence of Soft Modes and Electron Mobility in SrTiO 3 Perovskite from First Principles

Structural phase transitions and soft phonon modes pose a longstanding challenge to computing electron-phonon (e-ph) interactions in strongly anharmonic crystals. Here we develop a first-principles approach to compute e-ph scattering and charge transport in materials with anharmonic lattice dynamics. Our approach employs renormalized phonons to compute the temperature-dependent e-ph coupling for all phonon modes, including the soft modes associated with ferroelectricity and phase transitions. We show that the electron mobility in cubic SrTiO 3 is controlled by scattering with longitudinal optical phonons at room temperature and with ferroelectric soft phonons below 200 K. Our calculations can accurately predict the temperature dependence of the electron mobility in SrTiO 3 between 150-300 K, and reveal the microscopic origin of its roughly T -3 trend. Our approach enables first-principles calculations of e-ph interactions and charge transport in broad classes of crystals with phase transitions and strongly anharmonic phonons.

74 ATOMIC AND MOLECULAR PHYSICS↗

First-principles study, fabrication and characterization of (Zr 0.25 Nb 0.25 Ti 0.25 V 0.25 )C high-entropy ceramics

The formation possibility of a new (Zr 0.25 Nb 0.25 Ti 0.25 V 0.25 )C high-entropy ceramics (ZHC-1) was first analyzed by the first-principles calculations and thermodynamical analysis and then it was successfully fabricated by hot pressing sintering technique. The first-principles calculation results showed that the mixing enthalpy of ZHC-1 was 5.526 kJ/mol and the mixing entropy of ZHC-1 was in the range of 0.693R–1.040R. The thermodynamical analysis results showed that ZHC-1 was thermodynamically stable above 959 K owing to its negative mixing Gibbs free energy. The experimental results showed that the as-prepared ZHC-1 (95.1% relative density) possessed a single rock-salt crystal structure, some interesting nano-plate-like structures, and high compositional uniformity from nanoscale to microscale. By taking advantage of these unique features, compared with the initial metal carbides (ZrC, NbC, TiC and VC), it showed a relatively low thermal conductivity of 15.3 ± 0.3 W/(m&middot;K) at room temperature, which was due to the presence of solid solution effects, nanoplates and porosity. Meanwhile, it exhibited the relatively high nanohardness of 30.3 ± 0.7 GPa and elastic modulus of 460.4 ± 19.2 GPa and the higher fracture toughness of 4.7 ± 0.5 MPam 1/2 , which were attributed to the solid solution strengthening mechanism and nanoplate pullout and microcrack deflection toughening mechanism.

36 MATERIALS SCIENCE↗

Mechanistic Understanding of Uranyl Ion Complexation on Montmorillonite Edges: A Combined First-Principles Molecular Dynamics–Surface Complexation Modeling Approach

Systematic first-principles molecular dynamics (FPMD) simulations were carried out to study the structures, free energies, and acidity constants of UO 2 2+ surface complexes on montmorillonite in order to elucidate the surface complexation mechanisms of the uranyl ion (UO 2 2+ ) on clay mineral edges at the atomic scale. Four representative complexing sites were investigated, that is, ≡Al(OH) 2 and ≡AlOHSiO on the (010) surface and ≡AlOHOa and ≡SiOOa on the (110) surface. The results show that uranyl ions form bidentate complexes on these sites. We report all calculated binding free energies for these complexes are very similar. These bidentate complexes can be hydrolyzed, and their corresponding derived p K a values (around 5.0 and 9.0 for p K a1 and p K a2 , respectively) indicate that UO 2 (OH) + and UO 2 (OH) 2 surface groups are the dominant surface species in the environmental pH range. The OH groups of UO 2 (OH) 2 surface complexes can act as complexing sites for subsequent metals. Additional simulations showed that such multinuclear adsorption is feasible and can be important at high pH. Furthermore, FPMD simulation results served as input parameters for an electrostatic thermodynamic surface complexation model (SCM) that adequately reproduced adsorption data from the literature. Overall, this study provides an improved understanding of UO 2 2+ complexation on clay mineral edge surfaces.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Ethylene glycol reforming on Pt(111): first-principles microkinetic modeling in vapor and aqueous phases

First-principles, periodic density functional theory (DFT) calculations and mean-field microkinetic modeling have been used to investigate the decomposition of ethylene glycol for hydrogen production on Pt(111) in vapor and aqueous phases. All dehydrogenated species derived from ethylene glycol (C 2 H x O 2 , x = 0–6) and methanol (CH y O, y = 0–4), and all elementary C–C, C–H, and O–H bond breaking steps are included in the microkinetic model. Reaction path analysis in vapor phase indicates that both initial C–H and O–H dehydrogenation steps are kinetically relevant at all temperatures (470–530 K). Initial O–H bond cleavage is reversible at low temperatures but accounts for an increasingly dominant fraction of the total reaction flux at higher temperatures. C–C bond scission is observed only after significant dehydrogenation (x ≤ 3) and is more likely to happen in surface intermediates where one of the cleavage products is CO. Here, the process is highly selective to the production of H2 compared to methanol. For aqueous-phase model development, free energies of solvation were computed for all surface intermediates and transition states using a continuum solvation approach. Our aqueous-phase microkinetic model predicts a 0.4 eV lower apparent activation energy and an order of magnitude larger turnover frequencies. Initial C–H bond cleavage becomes more important but the general trends are similar to the vapor phase, suggesting that the reaction chemistry is similar in both vapor and aqueous phases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First principle study of manganese doped cadmium sulphide sheet

First-principle electronic structure calculations for cadmium sulphide (CdS) sheet in hexagonal phase, with Manganese substitution and addition, as well as including the Cd defects, are investigated. The lattice constants calculated for CdS sheet agrees fairly well with results reported for thin films experimentally. The calculations of total spin density of states and partial density of states in different cases shows substantial magnetic dipole moments acquired by the sheet. A magnetic dipole moment 5.00612 μ{sub B} and band gap of the order 1 eV are found when cadmium atom is replaced by Manganese. The magnetism acquired by the sheet makes it functionally important candidate in many applications.

36 MATERIALS SCIENCE↗