Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DENSITY FUNCTIONAL METHOD”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The Two-Energies Turbulence Scheme Coupled to the Assumed PDF Method

An update of the two-energy turbulence scheme is presented, the 2TE + APDF scheme. The original version of the two-energy scheme is able to successfully model shallow convection without the need of an additional parameterization for non-local fluxes. However, the performance of the two-energy scheme is worse in stratocumulus cases, where it tends to overestimate the erosion of the stable layers. We have identified the causes: the non-local stability parameter does not consider local stratification, the scheme lacks an internal parameter that could distinguish between a shallow convection regime and a stratocumulus regime, and it uses an inflexible turbulence length scale formulation. To alleviate this problem, we propose several modifications: an update of the stability parameter, a modified computation of the turbulence length scale, and the introduction of the entropy potential temperature to distinguish between a shallow convection and a stratocumulus regime. In addition, the two-energy scheme is coupled to a simplified assumed probability density function method in order to achieve a more universal representation of the cloudy regimes. The updated turbulence scheme is evaluated for several idealized cases and one selected real case in the ICOsahedral Nonhydrostatic (ICON) modeling framework. The results show that the updated scheme corrects the overmixing problem in the stratocumulus cases. The performance of the updated scheme is comparable to the operational setup, and can be thus used instead of the operational turbulence and shallow convection scheme in ICON. Additionally, the updated scheme improves the coupling with dynamics, which is beneficial for the modeling of coherent flow structures in the atmospheric boundary layer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Thermal effects on the electronic properties of sodium electride under high pressures

Sodium is a simple metal at ambient conditions, while it transits to an electride phase at pressures above ~160 GPa along the room temperature isotherm. Here we explore the thermal effects on the electronic properties of the $hP4$ phase of sodium along the $ρ$ = 5.872 g/cm 3 isochore. We quantitatively classify this phase as an insulator based on the criterion of nearsightedness of the one-particle density matrix. Ab initio calculations suggest that the band gap of this insulator decreases with increasing temperature along the isochore, primarily because of ionic distortions; culminating in an insulator-to-metal transition upon melting at $T_m \approx 2100$ $\hbox{K}$. This transition is accompanied by residual electronic localization (in d orbitals) in the form of dynamic electron bubbles and a change in hybridization form $p–d$ to $s–p$ upon melting. This transition is explored by tracking the electronic and electro-optical properties along the isochore under consideration.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Xenon Trapping in Metal-Supported Silica Nanocages

Xenon (Xe) is a valuable and scarce noble gas used in various applications, including lighting, electronics, and anesthetics, among many others. It is also a volatile byproduct of the nuclear fission of uranium. In this work, a novel material architecture consisting of silicate nanocages in contact with a metal surface and an approach for trapping single Xe atoms in these cages is presented. The trapping is done at low Xe pressures and temperatures between 400 and 600 K, and the process is monitored in situ using synchrotron-based ambient pressure X-ray photoelectron spectroscopy. Release of the Xe from the cages occurs only when heating to temperatures above 750 K. A model that explains the experimental trapping kinetics is proposed and tested using Monte Carlo methods. Density functional theory calculations show activation energies for Xe exiting the cages consistent with experiments. This work can have significant implications in various fields, including Xe production, nuclear power, nuclear waste remediation, and nonproliferation of nuclear weapons. The results are also expected to apply to argon, krypton, and radon, opening an even more comprehensive range of applications.

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

The influence of temperature on the strain-hardening behavior of Fe-22/25/28Mn-3Al-3Si TRIP/TWIP steels

In this work, the influence of temperature and stacking fault energy (SFE) on the strain-hardening behavior and critical resolved shear stress for twinning was investigated for three Fe–22/25/28Mn–3Al–3Si wt.% transformation- and twinning-induced plasticity (TRIP/TWIP) steels. The SFEs were calculated by two different methods, density functional theory and statistical thermodynamic modeling. The dislocation structure, observed at low levels of plastic deformation, transitions from “planar” to “wavy” dislocation glide with an increase in temperature, Mn content, and/or SFE. The change in dislocation glide mechanisms from planar to wavy reduces the strain hardening rate, in part due to fewer planar obstacles and greater cross slip activity. Additionally, the alloys exhibit a large decrease in strength and ductility with increasing temperature from 25 to 200 °C, attributed to a substantial reduction in the thermally activated component of the flow stress, predominate suppression of TRIP and TWIP, and a significant increase in the critical resolved shear stress for mechanical twinning. Interestingly, the increase in SFE with temperature had a rather minor influence on the critical resolved shear stress for mechanical twinning, and other temperature dependent factors which likely play a more dominant role are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Quantum-Chemically Informed Machine Learning: Prediction of Energies of Organic Molecules with 10 to 14 Non-hydrogen Atoms

High-fidelity quantum-chemical calculations can provide accurate predictions of molecular energies, but their high computational costs limit their utility, especially for larger molecules. We have shown in previous work that machine learning models trained on high-level quantum-chemical calculations (G4MP2) for organic molecules with one to nine non-hydrogen atoms can provide accurate predictions for other molecules of comparable size at much lower costs. Here we demonstrate that such models can also be used to effectively predict energies of molecules larger than those in the training set. To implement this strategy, we first established a set of 191 molecules with 10–14 non-hydrogen atoms having reliable experimental enthalpies of formation. We then assessed the accuracy of computed G4MP2 enthalpies of formation for these 191 molecules. The error in the G4MP2 results was somewhat larger than that for smaller molecules, and the reason for this increase is discussed. Two density functional methods, B3LYP and ωB97X-D, were also used on this set of molecules, with ωB97X-D found to perform better than B3LYP at predicting energies. The G4MP2 energies for the 191 molecules were then predicted using these two functionals with two machine learning methods, the FCHL-Δ and SchNet-Δ models, with the learning done on calculated energies of the one to nine non-hydrogen atom molecules. The better-performing model, FCHL-Δ, gave atomization energies of the 191 organic molecules with 10–14 non-hydrogen atoms within 0.4 kcal/mol of their G4MP2 energies. Thus, this work demonstrates that quantum-chemically informed machine learning can be used to successfully predict the energies of large organic molecules whose size is beyond that in the training set.

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

Crystal Structure of Colloidally Prepared Metastable Ag 2 Se Nanocrystals

Structural polymorphism is known for many bulk materials; however, on the nanoscale metastable polymorphs tend to form more readily than in the bulk, and with more structural variety. One such metastable polymorph observed for colloidal Ag 2 Se nanocrystals has traditionally been referred to as the “tetragonal” phase. While there are reports on the chemistry and properties of this metastable polymorph, its crystal structure, and therefore electronic structure, has yet to be determined. We report that an anti-PbCl 2 -like structure type (space group P2 1 /n) more accurately describes the powder X-ray diffraction and X-ray total scattering patterns of colloidal Ag 2 Se nanocrystals prepared by several different methods. Density functional theory (DFT) calculations indicate that this anti-PbCl 2 -like Ag 2 Se polymorph is a dynamically stable, narrow-band-gap semiconductor. Here, the anti-PbCl 2 -like structure of Ag 2 Se is a low-lying metastable polymorph at 5–25 meV/atom above the ground state, depending on the exchange-correlation functional used.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

The efficacy of Lewis affinity scale metrics to represent solvent interactions with reagent salts in all-inorganic metal halide perovskite solutions

Solvents employed in the solution processing of metal halide perovskites are known to play a key role in defining the morphology and properties of the resulting thin film, and thus the performance of perovskite solar cell devices. Accurate metrics are needed that are capable of differentiating among candidates, finding solvents that adequately solubilize the various precursor species in solution and facilitate the nucleation and growth of these materials. Existing metrics such as the unsaturated Mayer bond order (UMBO) and the Gutmann donor number (DN) have been tested for lead iodide perovskite systems; but there has yet to be a comprehensive study on their transferability to lead-free perovskite solutions. Here, we use ab initio methods (density functional theory) and regression analysis tools to study the usefulness of DN and BF 3 affinity scales in this regard. We compared the relative effectiveness of these scales to describe interactions between solvents and BXn perovskite salts of lead (Pb 2+ ), tin (Sn 2+ and Sn 4+ ), germanium (Ge 2+ ), bismuth (Bi 3+ ), and antimony (Sb 3+ and Sb 5+ ). The DN proved to be a better representation than the BF 3 of such interactions, reflecting the closer similarity of these species to the “parent” SbCl 5 Lewis acid than to BF 3 . In addition, we have uncovered the usefulness of the lithium cation affinity metric (LCA) to describe the strength of interactions between solvents and A-site cations (e.g. Na + , K + , Rb + and Cs + ) in all-inorganic metal halide perovskite solutions. We find that the coordination strengths of solvents towards species in all-inorganic metal halide perovskite solutions are best described by two different metrics with distinct modes of action: DN differentiates among BX n salt complexes, and LCA among A-site cation species. This revelation can help guide the choice of solvent to optimize processing conditions. It also emphasizes the importance of selecting solvents whose DN and LCA optimize coordination to key Lewis acid species in all-inorganic perovskite solutions.

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

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Perturbation theory and thermal transport in mass-disordered alloys: Insights from Green's function methods

Lowest-order quantum perturbation theory (Fermi's golden rule) for phonon-disorder scattering has been used to predict thermal conductivities in several semiconducting alloys with surprising success given its underlying hypothesis of weak and dilute disorder. In this paper, we explain how this is possible by focusing on the case of maximally mass-disordered Mg 2 Si 1-x Sn x . In this work, we use a Chebyshev polynomials Green's function method, which allows a full treatment of disorder on very large systems (tens of millions of atoms) to probe individual phonon linewidths and frequency-resolved thermal transport. We demonstrate that the success of perturbation theory originates from the specific form of mass disorder terms in the phonon Green's function and from the interplay between anharmonic and disorder scattering.

36 MATERIALS SCIENCE↗
DOE OSTI2023

First-principles derivation of magnetic interactions in the triangular quantum spin liquid candidates KYb C h 2 ( C h = S , Se , Te ) and A YbSe 2 ( A = Na , Rb )

The AYb Ch 2 (A=alkali metal, Ch=chalcogen) delafossites are a family of crystals which have a triangular lattice of effective S=1/2 moments on their Yb atoms. They hold great promise for the realization of the triangular quantum spin liquid state because of their defect-minimized growth and the ability to interchange chemical constituencies among the family. Here we use ab initio computations to evaluate the exchange couplings of four realized (and one theoretical) rhombohedral delafossite structures and examine the influence of chemical substitution on promoting the development of a quantum spin liquid state. We find good agreement with experiment regarding the antiferromagnetic nearest-neighbor exchange J 1 , but our calculations underestimate J 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

A density functional theory based comparative study of hybrid photoemissions from Cl@C 60 , Br@C 60 and I@C 60

Photoionization from atom-C 60 hybrid levels in halogen endofullerene molecules, Cl@C 60 , Br@C 60 and I@C 60 , are calculated using a linear response density functional method. Both the ordinary electron-configuration where the open shell halogen is at the center of C 60 and the stable configuration after the atom receives an electron from C 60 to form a closed shell anion are considered. Similar ground state hybridization is found for all three systems while, in general, a slight weakening of the effect is noticed after the electron transfer. At lower photon energies, cross sections of the outer hybrid levels attain identical shapes from enhancements driven by the C 60 plasmon resonances, while the higher energy emissions remain distinguishable from the differences in atomic responses. These results further show near insensitivity to the choice of a configuration. The inner hybrid cross sections in general exhibit similar overall structures, although differ in details between molecules. Furthermore, for these states the results significantly differ before and after the electron transfer – a feature that can be useful to experimentally determine the real configuration of the molecules via photoelectron spectroscopy.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Active Learning Framework

Machine learning (ML) of interatomic potentials show great promise to accelerate scientific simulation, e.g., by emulating expensive computations at a high accuracy but much reduced computational cost. Training datasets are calculated from computationally expensive ab initio quantum mechanics methods, density functional theory (DFT). Trained on this data, an ML model can be very successful in predicting energy and forces for new atomic configurations. A critical factor is the quality and diversity of the training dataset. Thus, a highly automated approach to dataset construction based on active learning framework is designed suitable for material physics. The active learning scheme begins with fully randomized atomic configurations. Then, many Molecular Dynamics (MD) trajectories are simulated using current ML potentials, where each MD trajectory is initialized to a random disordered configuration. The temperature is varied in order to diversify the sampled configuration during these simulations. The variance of predictions for eight neural networks within an ensemble is analyzed to determine whether the model is operating as expected. This helps in determining whether collecting more data would be helpful to the model by checking the ensemble variance is greater than the threshold. In this case, the MD trajectory is terminated and the final atomic configuration is placed on a queue (SQL database) for DFT calculations and added to training dataset. Periodically, ML model is retrained to the updated training model. This Active Learning loop is iterated until the cost of MD simulations becomes prohibitively expensive. The MD simulations will hopefully be sufficiently robust to support nucleation after many active learning iterations. In this sense, active learning scheme must automatically discover the important low energy and nonequilibrium physics.

Nebgen, Benjamin↗
DOE OSTI2024

Nuclear Computational Low Energy Initiative (NUCLEI)

The NUCLEI project, as defined by the scope of work, developed, implemented and run codes for large-scale computations of many topics in low-energy nuclear physics. Physics studied include the properties of nuclei and nuclear decays, nuclear structure and reactions, and the properties of nuclear matter. The computational techniques used include Quantum Monte Carlo, Configuration Interaction, Coupled Cluster, and Density Functional methods. The research program emphasized areas of high interest to current and possible future DOE nuclear physics facilities, including ATLAS and FRIB (nuclear structure and reactions, and nuclear astrophysics), TJNAF (neutron distributions in nuclei, few body systems, and electroweak processes), NIF (thermonuclear reactions), MAJORANA and FNPB (neutrinoless double-beta decay and physics beyond the Standard Model), and LANSCE (fission studies).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

On the Crystal Structure of Colloidally Prepared Metastable Ag2Se Nanocrystals

Structural polymorphism is known for many bulk materials; however, on the nanoscale metastable polymorphs tend to form more readily than in the bulk, and with more structural variety. One such metastable polymorph observed for colloidal Ag2Se nanocrystals has traditionally been referred to as the “tetragonal” phase of Ag2Se. While there are reports on the chemistry and properties of this metastable polymorph, its crystal structure, and therefore electronic structure, has yet to be determined. We report that an anti-PbCl2-like structure type (space group P21/n) accurately describes the powder X-ray diffraction and X-ray total scattering patterns of colloidal Ag2Se nanocrystals prepared by several different methods. Density functional theory (DFT) calculations indicate that the anti-PbCl2-like Ag2Se polymorph is a dynamically stable, narrow-band gap semiconductor. DFT results reveal a dense theoretical Ag2Se phase space with many low-energy polymorphs, which helps explain the large number of polymorphs reported in the literature. Analysis and calculation data are stored in the zip archive. The `ag2se-calcs.aiida.` contains the provenance of the calculations and can be imported into an AiiDA database instance. The antiPbCl2like_Ag2Se_laboratory.cif file is the Reitveld refined Ag2Se structure starting from the PbCl2 structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1996

Evaluation of joint probability density function models for turbulent nonpremixed combustion with complex chemistry

Two types of mixing sub-models are evaluated in connection with a joint-scalar probability density function method for turbulent nonpremixed combustion. Model calculations are made and compared to simulation results for homogeneously distributed methane-air reaction zones mixing and reacting in decaying turbulence within a two-dimensional enclosed domain. The comparison is arranged to ensure that both the simulation and model calculations a) make use of exactly the same chemical mechanism, b) do not involve non-unity Lewis number transport of species, and c) are free from radiation loss. The modified Curl mixing sub-model was found to provide superior predictive accuracy over the simple relaxation-to-mean submodel in the case studied. Accuracy to within 10-20% was found for global means of major species and temperature; however, nitric oxide prediction accuracy was lower and highly dependent on the choice of mixing sub-model. Both mixing submodels were found to produce non-physical mixing behavior for mixture fractions removed from the immediate reaction zone. A suggestion for a further modified Curl mixing sub-model is made in connection with earlier work done in the field.

Smith, N. S. A.↗
NASA NTRS2001

Endo-Fullerene and Doped Diamond Nanocrystallite Based Models of Qubits for Solid-State Quantum Computers

Models of encapsulated 1/2 nuclear spin H-1 and P-31 atoms in fullerene and diamond nanocrystallite, respectively, are proposed and examined with ab-initio local density functional method for possible applications as single quantum bits (qubits) in solid-state quantum computers. A H-1 atom encapsulated in a fully deuterated fullerene, C(sub 20)D(sub 20), forms the first model system and ab-initio calculation shows that H-1 atom is stable in atomic state at the center of the fullerene with a barrier of about 1 eV to escape. A P-31 atom positioned at the center of a diamond nanocrystallite is the second model system, and 3 1P atom is found to be stable at the substitutional site relative to interstitial sites by 15 eV, Vacancy formation energy is 6 eV in diamond so that substitutional P-31 atom will be stable against diffusion during the formation mechanisms within the nanocrystallite. The coupling between the nuclear spin and weakly bound (valance) donor electron coupling in both systems is found to be suitable for single qubit applications, where as the spatial distributions of (valance) donor electron wave functions are found to be preferentially spread along certain lattice directions facilitating two or more qubit applications. The feasibility of the fabrication pathways for both model solid-state qubit systems within practical quantum computers is discussed with in the context of our proposed solid-state qubits.

Park, Seongjun↗
DOE OSTI2023

Accelerating the density-functional tight-binding method using graphical processing units

Acceleration of the density-functional tight-binding (DFTB) method on single and multiple graphical processing units (GPUs) was accomplished using the MAGMA linear algebra library. Herein two major computational bottlenecks of DFTB ground-state calculations were addressed in our implementation: the Hamiltonian matrix diagonalization and the density matrix construction. The code was implemented and benchmarked on two different computer systems: (1) the SUMMIT IBM Power9 supercomputer at the Oak Ridge National Laboratory Leadership Computing Facility with 1–6 NVIDIA Volta V100 GPUs per computer node and (2) an in-house Intel Xeon computer with 1–2 NVIDIA Tesla P100 GPUs. The performance and parallel scalability were measured for three molecular models of 1-, 2-, and 3-dimensional chemical systems, represented by carbon nanotubes, covalent organic frameworks, and water clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗