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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Thermal transport in nanoporous holey silicon membranes investigated with optically induced transient thermal gratings

In this study, we use transient thermal gratings—a non-contact, laser-based thermal metrology technique with intrinsically high accuracy—to investigate room-temperature phonon-mediated thermal transport in two nanoporous holey silicon membranes with limiting dimensions of 120 nm and 250 nm, respectively. We compare the experimental results with ab initio calculations of phonon-mediated thermal transport according to the phonon Boltzmann transport equation (BTE) using two different computational techniques. We find that the calculations conducted within the Casimir framework, i.e., based on the BTE with the bulk phonon dispersion and diffuse scattering from surfaces, are in quantitative agreement with the experimental data and thus conclude that this framework is adequate for describing phonon-mediated thermal transport in silicon nanostructures with feature sizes of the order of 100 nm.

Duncan, Ryan A. (ORCID:0000000328318099)↗
DOE OSTI2020

Finite-temperature coupled cluster: Efficient implementation and application to prototypical systems

In this work, we discuss the theory and implementation of the finite temperature coupled cluster singles and doubles (FT-CCSD) method including the equations necessary for an efficient implementation of response properties. Numerical aspects of the method including the truncation of the orbital space and integration of the amplitude equations are tested on some simple systems, and we provide some guidelines for applying the method in practice. The method is then applied to the 1D Hubbard model, the uniform electron gas (UEG) at warm, dense conditions, and some simple materials. The performance of model systems at high temperatures is encouraging: for the one-dimensional Hubbard model, FT-CCSD provides a qualitatively accurate description of finite-temperature correlation effects even at U = 8, and it allows for the computation of systematically improvable exchange–correlation energies of the warm, dense UEG over a wide range of conditions. We highlight the obstacles that remain in using the method for realistic ab initio calculations on materials

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

A coupled cluster framework for electrons and phonons

Here, we describe a coupled cluster framework for coupled systems of electrons and harmonic phonons. Neutral and charged excitations are accessed via the equation-of-motion version of the theory. Benchmarks on the Hubbard–Holstein model allow us to assess the strengths and weaknesses of different coupled cluster approximations, which generally perform well for weak to moderate coupling. Finally, we report progress toward an implementation for ab initio calculations on solids and present some preliminary results on finite-size models of diamond with a linear electron–phonon coupling. We also report the implementation of electron–phonon coupling matrix elements from crystalline Gaussian type orbitals within the PySCF program package.

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

Ab-initio Cu alloy design for high-gradient accelerating structures

Operation of normal conducting accelerator structures at high accelerating gradients is beneficial for many accelerator applications in basic science, industry, medicine, and National Security. RF breakdown is the major factor that limits the achievable accelerating gradients. Previous experiments on copper (Cu) have demonstrated that RF breakdown probability can be significantly decreased by hardening the material and alloying Cu with solutes such as silver (Ag). In this paper, we propose a figure-of-merit (FOM) that characterizes the ability of Cu alloys to withstand high-gradients. The FOM represents a trade-off between hardening through solid solution strengthening and the additional thermal stress induced by incremental RF pulse heating resulting from changes in electronic properties induced by alloying. We performed high-throughput ab initio calculations and computed the FOM for a large number of binary Cu alloys. Several promising candidate alloys for high-gradient accelerating structures were identified, such as CuAg, CuCd, CuHg, CuAu, CuIn, and CuMg. CuAg alloys have previously exhibited low RF breakdown rates in experiments. The results provide guidance for selecting alloys for the future high-gradient normal conducting accelerating structures operating at very high gradients.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ab initio investigation of elastic properties of dilute Cu alloys for high-gradient accelerating structures

The elastic properties and solid solution strengthening of several Cu alloys that may be promising for high gradient applications were studied using ab initio calculations based on density functional theory. Specifically, the temperature and solute concentration dependence of elastic properties including bulk and shear moduli, thermal expansion coefficient, and heat capacity of these alloys were calculated. We found that the bulk and shear moduli are linearly dependent on the solute concentration in the dilute regime and change significantly for different Cu alloys. The thermal expansion coefficient and heat capacity are comparatively less sensitive to solute concentrations. Solid solution strengthening is observed to be linearly correlated with the change of lattice parameter of Cu. This work provides an in-depth understanding of the elastic properties of Cu alloys, which is key for predicting their performance at high gradient in normal conducting accelerating structures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Acceptor and compensating donor doping of single crystalline SnO (001) films grown by molecular beam epitaxy and its perspectives for optoelectronics and gas-sensing

(La and Ga)-doped tin monoxide [stannous oxide, tin (II) oxide, SnO] thin films were grown by plasma-assisted and suboxide molecular beam epitaxy with dopant concentrations ranging from ≈ 5 × 10 18 to 2 × 10 21 cm -3 . In this concentration range, the incorporation of Ga into SnO was limited by the formation of secondary phases observed at 1.2 × 10 21 cm -3 Ga, while the incorporation of La showed a lower solubility limit. Transport measurements on the doped samples reveal that Ga acts as an acceptor and La as a compensating donor. While Ga doping led to an increase in the hole concentration from 1 × 10 18 -1 × 10 19 cm -3 for unintentionally doped (UID) SnO up to 5 × 10 19 cm -3 , La-concentrations well in excess of the UID acceptor concentration resulted in semi-insulating films without detectable n-type conductivity. Ab initio calculations qualitatively agree with our dopant assignment of Ga and La and further predict In Sn to act as an acceptor as well as Al Sn and B Sn as donors. These results show the possibilities of controlling the hole concentration in p-type SnO, which can be useful for a range of optoelectronic and gas-sensing applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Synergy of semiempirical models and machine learning in computational chemistry

Catalyzed by enormous success in the industrial sector, many research programs have been exploring data-driven, machine learning approaches. Performance can be poor when the model is extrapolated to new regions of chemical space, e.g., new bonding types, new many-body interactions. Another important limitation is the spatial locality assumption in model architecture, and this limitation cannot be overcome with larger or more diverse datasets. The outlined challenges are primarily associated with the lack of electronic structure information in surrogate models such as interatomic potentials. Given the fast development of machine learning and computational chemistry methods, we expect some limitations of surrogate models to be addressed in the near future; nevertheless spatial locality assumption will likely remain a limiting factor for their transferability. Here, we suggest focusing on an equally important effort—design of physics-informed models that leverage the domain knowledge and employ machine learning only as a corrective tool. In the context of material science, we will focus on semi-empirical quantum mechanics, using machine learning to predict corrections to the reduced-order Hamiltonian model parameters. The resulting models are broadly applicable, retain the speed of semiempirical chemistry, and frequently achieve accuracy on par with much more expensive ab initio calculations. These early results indicate that future work, in which machine learning and quantum chemistry methods are developed jointly, may provide the best of all worlds for chemistry applications that demand both high accuracy and high numerical efficiency.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Influence of fourth-order vibrational corrections on semi-experimental (reSE) structures of linear molecules

Semi-experimental structures (reSE) are derived from experimental ground state rotational constants combined with theoretical vibrational corrections. They permit a meaningful comparison with equilibrium structures based on high-level ab initio calculations. Typically, the vibrational corrections are evaluated with second-order vibrational perturbation theory (VPT2). The amount of error introduced by this approximation is generally thought to be small; however, it has not been thoroughly quantified. Herein, we assess the accuracy of theoretical vibrational corrections by extending the treatment to fourth order (VPT4) for a series of small linear molecules. Typical corrections to bond distances are on the order of 10−5 Å. Larger corrections, nearly 0.0002 Å, are obtained to the bond lengths of NCCN and CNCN. A borderline case is CCCO, which will likely require variational computations for a satisfactory answer. Treatment of vibrational effects beyond VPT2 will thus be important when one wishes to know bond distances confidently to four decimal places (10−4 Å). Certain molecules with shallow bending potentials, e.g., HOC+, are not amenable to a VPT2 description and are not improved by VPT4.

Chemistry↗
DOE OSTI2024

Deprotonated sulfamic acid and its homodimers: Does sulfamic acid adopt zwitterion during cluster growth?

Here we present a joint experimental and computational study on the geometric and electronic structures of deprotonated sulfamic acid (SA) clusters [(SA) n –H] – (n = 1, 2) employing negative ion photoelectron spectroscopy and high-level ab initio calculations. The photoelectron spectra provide the vertical/adiabatic detachment energy (VDE/ADE) of the sulfamate anion (SM – ) H 2 N•SO 3 – at 4.85 ± 0.05 and 4.58 ± 0.08 eV, respectively, and the VDE and ADE of the SM – •SA dimer at 6.41 ± 0.05 and 5.87 ± 0.08 eV, respectively. The significantly increased electron binding energies of the dimer confirm the enhanced electronic stability upon the addition of one SA molecule. The CCSD(T)-predicted VDEs/ADEs agree excellently with the experimental data, confirming the identified structures as the most stable ones. Two types of dimer isomers possessing different hydrogen bonding (HB) motifs are identified, corresponding to SM – binding to a zwitterionic SA (SM – •SA z ) and a canonical SA (SM – •SA c ), respectively. Two N–H$\cdots$O HBs and one superior O–H$\cdots$O HB are formed in the lowest-lying SM – •SA c , while SM – •SA z has three moderate N–H$\cdots$O HBs, with the former being 4.71 kcal/mol more stable. Further theoretical analyses reveal that the binding strength advantage of SM – •SA c over SM – •SA z arises from its significant contributions of orbital interactions between fragments, illustrating that sulfamate strongly interacts with its parent SA acid and preferably chooses the canonical SA in the subsequent cluster formations. Given the prominent presence of SA, this study provides the first evidence that the canonical dimer model of sulfamic acid should exist as a superior configuration during cluster growth.

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

Cryogenic temperatures promote the pressure-induced polymorphic transition in CoCrFeMnNi high entropy alloy

Pressure-induced polymorphism has recently been demonstrated in several high entropy alloys. This offers a new window into the much-debated issue of phase selection and stability in these systems. Here, we examine the effect of cryogenic temperatures on the pressure-induced transition from face centered cubic to hexagonal close-packed structures of the prototype CoCrFeMnNi (Cantor) alloy. We observe a reduction in the critical pressure for the onset of the polymorphic transition as the temperature decreases, confirming the progressive stabilization of the hexagonal phase with decreasing temperature previously predicted by ab initio calculations accounting for magnetic interactions. We argue that in situ high-pressure experiments at cryogenic temperatures, which suppress time-dependent transformation triggered at higher temperatures, present a unique opportunity to significantly improve our understanding of these complex alloys.

alloys↗
DOE OSTI2025

Metal-organic chemical vapor deposition of MgGeN2 films on GaN and sapphire

MgGeN2 films were synthesized using metal-organic chemical vapor deposition on GaN/c-sapphire templates and c-plane sapphire substrates. Energy-dispersive x-ray spectroscopy was used to estimate the cation composition ratios. To mitigate magnesium evaporation, the films were grown at pyrometer temperature 745 °C with a wafer rotation speed of 1000 rpm. Growth rates were determined by fitting energy-dispersive x-ray spectroscopy spectra to film thicknesses using NIST DTSA-II software. The thickness estimates determined by this method were consistent with scanning transmission electron microscopy measurements done for selected samples. Scanning electron microscopy images revealed faceted surfaces indicative of a tendency toward three-dimensional growth. X-ray diffraction spectra confirmed that the films were highly crystalline and exhibited preferential orientation in alignment with the substrate. Atomic force microscopy measurements show that film thicknesses are consistent across samples grown on both GaN templates and sapphire substrates, with typical roughnesses around 10 nm. Transmittance spectra of films grown on double-side-polished sapphire substrates yielded band gaps of 4.28 ± 0.06 eV for samples exhibiting close-to-ideal stoichiometry. Comparison of the measured spectra with ab initio calculations are in good agreement both near the bandgap and at higher energies where excitation is into higher-lying bands. These findings provide insight into the growth and characterization of MgGeN2, contributing to the development of this material for potential applications in optoelectronics and power electronics.

Hu, Chenxi (ORCID:0000000236299660)↗
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 OSTI2023

PBE-GGA predicts the B8↔B2 phase boundary of FeO at Earth’s core conditions

FeO is a crucial component of the Earth’s core, and its thermodynamic properties are essential to developing more accurate core models. It is also a notorious correlated insulator in the NaCl-type (B1) phase at ambient conditions. It undergoes two polymorphic transitions at 300 K before it becomes metallic in the NiAs-type (B8) structure at ~100 GPa. Although its phase diagram is not fully mapped, it is well established that the B8 phase transforms to the CsCl-type (B2) phase at core pressures and temperatures. Here, we report a successful ab initio calculation of the B8↔B2 phase boundary in FeO at Earth’s core pressures. We show that fully anharmonic free energies computed with the Perdew–Burke–Ernzerhof-generalized gradient approximation coupled with thermal electronic excitations reproduce the experimental phase boundary within uncertainties at P > 255 GPa, including the largely negative Clapeyron slope of –52 MPa/K. This study validates the applicability of a standard density functional theory functional to FeO under Earth’s core conditions and demonstrates the theoretical framework that enables complex predictive studies of this region.

58 GEOSCIENCES↗
DOE OSTI2023

Direct observation of the collective modes of the charge density wave in the kagome metal CsV 3 Sb 5

A recently discovered group of kagome metals AV 3 Sb 5 (A = K, Rb, Cs) exhibit a variety of intertwined unconventional electronic phases, which emerge from a puzzling charge density wave phase. Understanding of this charge-ordered parent phase is crucial for deciphering the entire phase diagram. However, the mechanism of the charge density wave is still controversial, and its primary source of fluctuations—the collective modes—has not been experimentally observed. Here, we use ultrashort laser pulses to melt the charge order in CsV 3 Sb 5 and record the resulting dynamics using femtosecond angle-resolved photoemission. We resolve the melting time of the charge order and directly observe its amplitude mode, imposing a fundamental limit for the fastest possible lattice rearrangement time. These observations together with ab initio calculations provide clear evidence for a structural rather than electronic mechanism of the charge density wave. Our findings pave the way for a better understanding of the unconventional phases hosted on the kagome lattice.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Extended Rice–Thomson analysis and atomistic simulations revealing grain boundary effects on fracture in refractory high-entropy alloys

Significance This work serves to extend the fundamental ductile vs. brittle fracture theory, specifically the Rice–Thomson criterion, by introducing a grain boundary ahead of an initiating crack which propagates at an oblique angle to impinge the boundary. Atomistic fracture simulations on two refractory complex concentrated alloys, the brittle NbMoTaW and the ductile Nb 45 Ta 25 Ti 15 Hf 15 , demonstrate qualitative correspondence with the extended Rice–Thomson criterion and experimental observations. Abstract Understanding how grain boundaries mediate fracture remains a critical challenge in designing ductile, high-performance refractory alloys. Here, we extend the Rice–Thomson criterion to account for the angle between cracks and the impinging grain boundaries (GBs), capturing the competition between intergranular fracture and dislocation-mediated plasticity. Using machine learning interatomic potentials, we performed molecular statics simulations to probe fracture mechanisms in nanocrystalline NbMoTaW and Nb 45 Ta 25 Ti 15 Hf 15 , each with two different grain sizes, revealing trends consistent with experimental observations and the extended Rice model. Comparison with averaged R-curves for bulk samples demonstrates that GBs enhance ductility in Nb 45 Ta 25 Ti 15 Hf 15 in both grain sizes investigated. In contrast, GBs only locally improve fracture resistance in NbMoTaW when cracks are temporarily pinned at GBs inclined at high angles from the crack, but generally promote brittle intergranular fracture. These contrasting behaviors are attributed to differences in GB cohesion, reflecting clear alloying trends that align with ab-initio calculations and trends observed experimentally. Our results bridge classical fracture theory, atomistic simulations, and experimental observations, providing a comprehensive understanding of the fracture mechanisms in nanocrystalline refractory complex concentrated alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2020

High intensity x-ray interaction with a model bio-molecule system: double-core-hole states and fragmentation of formamide

Formamide, a simple model bio-molecule (HCONH 2 ), is irradiated with high intensity, ultrashort pulses from an x-ray free electron laser. Ionic fragments resulting from photoionization and subsequent decay processes are recorded, as well as the electronic signature of the different inner shell ionization events that can take place during the x-ray pulses. Additionally, the formation of double-core-hole states, where a second inner shell electron is removed before the first core hole has been refilled is observed in the electron spectra, recorded at all three sites (C, N, O) of the molecule. The individual ionization potentials are compared with results of ab initio calculations at different level of theory. Based on our results, future opportunities for advanced studies of inner-shell-induced electronic and nuclear dynamics are explored.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Observed metallization of hydrogen interpreted as a band structure effect

A recent experimental study of the metallization of hydrogen tracked the direct band gap and vibron frequency via infrared measurements up to ~425 GPa (Loubeyre et al (2020 Nature 577 631). Above this pressure, the direct gap has a discontinuous drop to below the minimum experimentally accessible energy (~0.1 eV). The authors suggested that this observation is caused by a structural phase transition between the C2/c-24 molecular phase to another molecular phase such as Cmca-12. Here, through ab initio calculations of pressure dependent vibron frequency and direct band gap, we find that the experimental data is consistent with the C2/c-24 phase up to 425 GPa, and suggest that this consistency extends beyond that pressure. Specifically, we find that qualitative changes in the band structure of the C2/c-24 phase lead to a discontinuous drop of the direct band gap, which can explain the observed drop without a structural transition. In conclusion, this alternative scenario, which naturally explains the absence of hysteresis in the measurements, will hopefully motivate further experimental studies to ascertain the structure of the phase above the high pressure 'phase transition'.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

The effect of Cu additions in FePt–BN–SiO 2 heat-assisted magnetic recording media

Structural and chemical order impact magnetic properties of solids, which are governed by spin–orbit coupling and exchange interaction. The ordered L1 0 phase of FePt is a key material to heat-assisted magnetic recording; to enable high storage density, a solid understanding is needed of structural and chemical disorder at small length scales, as well as associated modifications of the electronic band structure. Here, we investigate the effect of boron and copper additions (≲6 mol% Cu) on structural and magnetic properties of L1 0 FePt granular media. Two copper-driven mechanisms, although competing, can lead to improvements in both structural and magnetic properties. In particular, the Cu substitution on the Fe-site leads to a degradation of magnetic properties due to the delocalized electron orbitals originating from a larger Cu d -orbital occupancy. At the same time, Cu substitution leads to an enhanced crystallographic order and consequently magneto-crystalline anisotropy, which offsets the former effect to a large extent. Our study is based on magnetometry, x-ray absorption spectroscopy, ab-initio calculations and a phenomenological theory of disordered FePt granular media. We do not observe a sizable modification to Fe moments and electronic configuration; Cu reveals two different resonances associated with the presence and absence of Cu–B bonds that vary with total Cu concentration.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗