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 271 records · Page 15

First-principles analysis of anharmonic nuclear motion and thermal transport in thermoelectric materials

We show a first-principles approach for analyzing anharmonic properties of lattice vibrations in solids. We firstly extract harmonic and anharmonic force constants from accurate first-principles calculations based on the density functional theory. Using the many-body perturbation theory of phonons, we then estimate the phonon scattering probability due to anharmonic phonon-phonon interactions. We show the validity of the approach by computing the lattice thermal conductivity of Si, a typical covalent semiconductor, and selected thermoelectric materials PbTe and Bi{sub 2}Te{sub 3} based on the Boltzmann transport equation. We also show that the phonon lifetime and the lattice thermal conductivity of the high-temperature phase of SrTiO{sub 3} can be estimated by employing the perturbation theory on top of the solution of the self-consistent phonon equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

First-principles calculations and thermodynamic modeling of the Yb-Ni binary system

A complete thermodynamic description of the Yb-Ni binary system is developed by means of the CALculation of PHAse Diagram (CALPHAD) method in combination with first-principles calculations based on density functional theory (DFT) and available experimental data in the literature. Finite temperature thermodynamic properties of the Yb-Ni intermetallic compounds are predicted using the quasi-harmonic approach, where first-principles phonon calculations are performed to calculate the lattice vibrational entropy. Here, the associate solution model is used to describe the complex thermodynamic behavior of the liquid phase. The calculated phase diagram agrees well the experimental phase equilibrium data in the literature. By the coupling of CALPHAD modeling with first-principles calculations, the present work provides a more thermodynamically accurate model of the Yb-Ni system when compared to previous models.

36 MATERIALS SCIENCE↗

Comparative Accuracies of Models for Drag Prediction During Geomagnetically Disturbed Periods: A First Principles Model Versus Empirical Models

We examine the accuracy of density prediction by the first principles model Thermosphere Ionosphere Electrodynamics General Circulation Model (TIEGCM) developed by the National Center for Atmospheric Research and compare it to the accuracy of three empirical models: Jacchia 71, the Naval Research Laboratory Mass Spectrometer Incoherent Scatter Extended 2000 (NRLMSIS), Jacchia 1971, and Jacchia-Bowman 2008. Comparisons are made for three large storms: the October 2003 storm, the March 2013 storm, and the March 2015 storm. To evaluate the accuracy of these models we use tracking data for nine space objects in low Earth orbit. Additionally, we evaluate the accuracy of the TIEGCM and NRLMSIS with data from high precision accelerometers on the Challenging Minisatellite Payload (CHAMP) and Gravity field and Circulation Explorer (GOCE) satellites. The goal is to assess the use of a first principles model as a potential tool for forecasting satellite drag during large magnetic storms. For the storms considered, we found the TIEGCM, JB2008, and NRLMSIS models to be substantially more accurate than the Jacchia 71 model. The accuracies of the TIEGCM and JB2008 models were similar, but overall, the TIEGCM was more accurate. We found smaller differences for TIEGCM versus CHAMP than for NRLMIS for the Halloween Storm, and smaller differences than results published for JB2008 and the assimilative model HASDM. The empirical models are at present more practical for operational purposes, but the TIEGCM, developed as a research model, with a greater focus on operational use offers the potential for improved utility during stressing conditions.

R. L. Walterscheid↗

Revisiting W–ZrC interfaces: A first principles study

We investigate W–ZrC interfaces using first-principles calculations based on the density functional theory. There have been theoretical and experimental studies exploring W–ZrC interfaces, however, the debate regarding the most stable interface continues to persist. In this study, we systematically simulated various W–ZrC interfaces merging W and ZrC surfaces with different orientations. Subsequently, we evaluated their stabilities and explained the corresponding stabilities in terms of the nature of bonding and charge-transfer processes at the interface. We find ZrC(111)–W(110) is the most stable interface with higher adhesive energy than the other interfaces. The additional stability associated with the ZrC(111)–W(110) results from significant interface reconstruction. Three layers of W and ZrC adjacent to the interface are involved in the charge-transfer process leading to stronger ionic bonds in ZrC(111)–W(110) as compared to the other potential candidate: ZrC(100)–W(100). The C and W atoms are found to be displaced from their symmetric position during the reconstruction process at the interface to facilitate stronger bonds with shorter W–C and W–Zr bonds in ZrC(111)–W(110) as compared to ZrC(100)–W(100). This leads to stronger covalent bonds in ZrC(111)–W(110) than that in ZrC(100)–W(100). Therefore, we conclude that the stronger covalent and ionic forces in ZrC(111)–W(110) than those in ZrC(100)–W(100) are responsible for making ZrC(111)–W(110) to be the most stable interface. In conclusion, this study addresses the long-standing question of the most stable W–ZrC interface and derives a number of implications for other W-transition metal carbide interfaces which are potential candidates for improving the mechanical properties of plasma facing materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

First principle calculations of effective exchange integrals: Comparison between SR (BS) and MR computational results

First principle calculations of effective exchange integrals (J) in the Heisenberg model for diradical species were performed by both symmetry-adapted (SA) multi-reference (MR) and broken-symmetry (BS) single reference (SR) methods. Mukherjee-type (Mk) state specific (SS) MR coupled-cluster (CC) calculations by the use of natural orbital (NO) references of ROHF, UHF, UDFT and CASSCF solutions were carried out to elucidate J values for di- and poly-radical species. Spin-unrestricted Hartree Fock (UHF) based coupled-cluster (CC) computations were also performed to these species. Comparison between UHF-NO(UNO)-MkMRCC and BS UHF-CC computational results indicated that spin-contamination of UHF-CC solutions still remains at the SD level. In order to eliminate the spin contamination, approximate spin-projection (AP) scheme was applied for UCC, and the AP procedure indeed corrected the error to yield good agreement with MkMRCC in energy. The CC double with spin-unrestricted Brueckner's orbital (UBD) was furthermore employed for these species, showing that spin-contamination involved in UHF solutions is largely suppressed, and therefore AP scheme for UBCCD removed easily the rest of spin-contamination. We also performed spin-unrestricted pure- and hybrid-density functional theory (UDFT) calculations of diradical and polyradical species. Three different computational schemes for total spin angular momentums were examined for the AP correction of the hybrid (H) UDFT. HUDFT calculations followed by AP, HUDFT(AP), yielded the S-T gaps that were qualitatively in good agreement with those of MkMRCCSD, UHF-CC(AP) and UB-CC(AP). Thus a systematic comparison among MkMRCCSD, UCC(AP) UBD(AP) and UDFT(AP) was performed concerning with the first principle calculations of J values in di- and poly-radical species. It was found that BS (AP) methods reproduce MkMRCCSD results, indicating their applicability to large exchange coupled systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

GPU acceleration of the Locally Selfconsistent Multiple Scattering code for first principles calculation of the ground state and statistical physics of materials

The Locally Self-consistent Multiple Scattering (LSMS) code solves the first principles Density Functional theory Kohn–Sham equation for a wide range of materials with a special focus on metals, alloys and metallic nano-structures. It has traditionally exhibited near perfect scalability on massively parallel high performance computer architectures. In this paper, we present our efforts to exploit GPUs to accelerate the LSMS code to enable first principles calculations of O(100,000) atoms and statistical physics sampling of finite temperature properties. We reimplement the scattering matrix calculation for GPUs with a block matrix inversion algorithm that only uses accelerator memory. Finally, using the Cray XK7 system Titan at the Oak Ridge Leadership Computing Facility we achieve a sustained performance of 14.5PFlop/s and a speedup of 8.6 compared to the CPU only code.

36 MATERIALS SCIENCE↗

First-principles thermodynamic theory of Seebeck coefficients

To begin, thermoelectric effects, measured by the Seebeck coefficients, refer to the phenomena in which a temperature difference or gradient imposed across a thermoelectric material induces an electrical potential difference or gradient, and vice versa, enabling the direct conversion of thermal and electric energies. All existing first-principles calculations of Seebeck coefficients have been based on the Boltzmann kinetic transport theory. In this work, we present a fundamentally different method for the first-principles calculations of Seebeck coefficients without using any assumptions of the electron-scattering mechanism, being in contrast to the traditional theory by Cutler and Mott that shows the dependence of the Seebeck coefficient on the scattering mechanisms. It is shown that the Seebeck coefficient is a well-defined thermodynamic quantity that can be determined from the change in the chemical potential of electrons induced by the temperature change and thus can be computed solely based on the electronic density of states through first-principles calculations at different temperatures. The proposed approach is demonstrated using the prototype PbTe and SnSe thermoelectric materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Time-dependent first-principles approaches to PV materials

Computational scheme for designing photovoltaic (PV) materials is presented. First-principles electron dynamics of photo-excitation and subsequent electron-hole splitting is performed based on the time-dependent density functional theory. Photo-induced enhancement of dipole moment was observed in a polar crystal and a donor-acceptor molecular pair. These experiences will pave a way to design PV material from first-principles simulations.

36 MATERIALS SCIENCE↗

Active Learning A Neural Network Model For Gold Clusters & Bulk From Sparse First Principles Training Data

Small metal clusters are of fundamental scientific interest and of tremendous significance in catalysis. These nanoscale clusters display diverse geometries and structural motifs depending on the cluster size; a knowledge of this size-dependent structural motifs and their dynamical evolution has been of longstanding interest. Given the high computational cost of first-principles calculations, molecular modeling and atomistic simulations such as molecular dynamics (MD) has proven to be an important complementary tool to aid this understanding. Classical MD typically employ predefined functional forms which limits their ability to capture such complex size-dependent structural and dynamical transformation. Neural Network (NN) based potentials represent flexible alternatives and in principle, well-trained NN potentials can provide high level of flexibility, transferability and accuracy on-par with the reference model used for training. A major challenge, however, is that NN models are interpolative and requires large quantities (similar to 10 4 or greater) of training data to ensure that the model adequately samples the energy landscape both near and far-from-equilibrium. A highly desirable goal is minimize the number of training data, especially if the underlying reference model is first-principles based and hence expensive. In this work, we introduce an active learning (AL) scheme that trains a NN model on-the-fly with minimal amount of first-principles based training data. Our AL workflow is initiated with a sparse training dataset (similar to 1 to 5 data points) and is updated on-the-fly via a Nested Ensemble Monte Carlo scheme that iteratively queries the energy landscape in regions of failure and updates the training pool to improve the network performance. Using a representative system of gold clusters, we demonstrate that our AL workflow can train a NN with similar to 500 total reference calculations. Using an extensive DFT test set of similar to 1100 configurations, we show that our AL-NN is able to accurately predict both the DFT energies and the forces for clusters of a myriad of different sizes. Our NN predictions are within 30 meV/atom and 40 meV/angstrom of the reference DFT calculations. Moreover, our AL-NN model also adequately captures the various size-dependent structural and dynamical properties of gold clusters in excellent agreement with DFT calculations and available experiments. We finally show that our AL-NN model also captures bulk properties reasonably well, even though they were not included in the training data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First-Principles Cost Estimation of a Sodium Fast Reactor Nuclear Plant

A multi-tiered cost analysis is performed to estimate full costs of a nuclear power plant (NPP) based on sodium-cooled fast reactor (SFR) technology. To address the lack of fully transparent cost estimations from past undertakings for NPPs, we have developed a detailed and first-principles-based cost estimate for a generalized SFR NPP. Our intent is to achieve a high degree of transparency with our cost assumptions and develop a cost model that is flexible and easily extendable to variations in NPP design and other nuclear reactor types. Furthermore, we strive to achieve a clear organization of costs and complete identification of key cost drivers based on first principles. To this end, the cost results of our analysis as given in Table 24 and Table 25 are organized and categorized into a code of accounts (COA) under development at Idaho National Laboratory (INL). Varying degrees of first-principles methods are employed, such as design for manufacture and assembly® (DFMA® ), to elucidate costs in all process levels of the plant equipment, buildings and site structures, personnel, and other miscellaneous but significant cost elements. These approaches have been successfully applied in past cost analysis projects and are designed to enable rapid and flexible cost estimation. Application of these techniques for evaluating NPP costs is similar in concept to the full, detailed estimation of construction and fabrication costs determined in a later stage of NPP development. Note that our approach tries to avoid use of other past analysis results and data such as those from the legacy Energy Economic Data Base (EEDB) Program, as these resources are based on historical NPP costs and thus may not be indicative of new reactor technologies or construction and fabrication/manufacturing techniques. However, we provide a comparison of our SFR NPP cost results in Table 90 against those included in the EEDB for a representative pressurized water reactor (PWR).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

First-principles predictions of temperature-dependent infrared dielectric function of polar materials by including four-phonon scattering and phonon frequency shift

Recently, first-principles calculations based on density functional theory have been widely used to predict the temperature-dependent infrared spectrum of polar materials, but the calculations are usually limited to the harmonic frequency (0 K) and three-phonon scattering damping for the zone-center infrared-active optical phonon modes, and fail to predict the high-temperature infrared optical properties of materials such as sapphire (α–Al 2 O 3 ), GaAs, TiO 2 , etc., due to the neglect of high-order phonon scattering damping and phonon frequency shift. In this work, we implemented first-principles calculations to predict the temperature-dependent infrared dielectric function of polar materials by including four-phonon scattering and phonon frequency shift. The temperature-dependent phonon damping by including three- and four-phonon scattering as well as the phonon frequency shift by including cubic and quartic anharmonicity and the thermal expansion effect are calculated based on anharmonic lattice dynamics method. The infrared dielectric function of α–Al 2 O 3 is parameterized, and then the temperature-dependent infrared optical reflectance is determined. We find that our predictions agree better with the experimental data than the previous density functional theory-based methods. Lastly, this work will help to effectively predict the thermal radiative properties of polar materials at elevated temperature, which is generally difficult to measure, and will enable predictive design of new materials for radiative applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Magnetic materials at finite temperatures: thermodynamics and combined spin and molecular dynamics derived from first principles calculations

We present a unified approach to describe the combined behavior of the atomic and magnetic degrees of freedom in magnetic materials. Using Monte Carlo simulations directly combined with first principles the Curie temperature can be obtained ab initio in good agreement with experimental values. The large scale constrained first principles calculations have been used to construct effective potentials for both the atomic and magnetic degrees of freedom that allow the unified study of influence of phonon-magnon coupling on the thermodynamics and dynamics of magnetic systems. The MC calculations predict the specific heat of iron in near perfect agreement with experimental results from 300K to above Tc and allow the identification of the importance of the magnon-phonon interaction at the phase-transition. Further Molecular Dynamics and Spin Dynamics calculations elucidate the dynamics of this coupling and open the potential for quantitative and predictive descriptions of dynamic structure factors in magnetic materials using first principles-derived simulations.

36 MATERIALS SCIENCE↗

Stability of vacancy-type defect clusters in Ni based on first-principles and molecular dynamics simulations

Using first-principles calculations based on density-functional theory, the energetics of different vacancy-type defects, including voids, stacking fault tetrahedra (SFT) and vacancy loops, in Ni are investigated. It is found that voids are more stable than SFT at 0 K, which is also the case after taking into account the volumetric strains. By carrying out ab initio molecular dynamics simulations at temperatures up to 1000 K, direct transformations from vacancy loops and voids into SFT are observed. Our results suggest the importance of temperature effects in determining thermodynamic stability of vacancy clusters in face-centered cubic metals.

36 MATERIALS SCIENCE↗

Magnetic properties of transition metal doped AlN nanosheet: First-principle studies

We carry out our first-principles calculations within density functional theory to study the 3d transition metal (TM) doped AlN nanosheets. The calculated results indicate that a stoichiometric AlN nanosheet is graphene-like structure and nonmagnetic. The TM impurities can induce magnetic moments, localized mainly on the 3d TM atoms and neighboring N atoms. Our calculated results of TM-doped nanosheet systems indicate a strong interaction between 3d orbit of TM atom and the 2p orbit of N atoms. In addition, the Mn- and Ni-doped AlN nanosheet with half-metal characters seems to be good candidates for spintronic applications. When substituting two Al atoms, the relative energies of the states between ferromagnetic and antiferromagnetic coupling are investigated sufficiently. The exchange coupling of Co- and Ni-doped AlN nanosheets exhibits a transformation with different distances of two TM atoms and that of Cr-, Mn-, and Fe-doped AlN nanosheets is not changed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electronic structures and optical spectra of BaO from first principles

We present the results of first-principles study for the electronic structure and optical absorption spectrum of the alkaline-earth metal oxide BaO. The quasiparticle band structure is evaluated within the Hedin's GW approximation [Phys. Rev. 139, A796 (1965)]. Thereafter, the electron-hole interaction is taken into consideration and the Bethe-Salpeter equation for the electron-hole two-particle Green function is solved. The calculated quasiparticle band gap of BaO is 4.1 eV, which is in good agreement with the experimental result. The calculated optical absorption spectrum of BaO is also in agreement with the experimental data. In particular, the calculated excitation energy for the lowest exciton peak in the optical absorption spectrum of BaO reproduces very well the corresponding experimental result.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Phonon-Assisted Auger-Meitner Recombination in Silicon from First Principles

Here, we present a consistent first-principles methodology to study both direct and phonon-assisted Auger-Meitner recombination (AMR) in indirect-gap semiconductors that we apply to investigate the microscopic origin of AMR processes in silicon. Our results are in excellent agreement with experimental measurements and show that phonon-assisted contributions dominate the recombination rate in both n-type and p-type silicon, demonstrating the critical role of phonons in enabling AMR. We also decompose the overall rates into contributions from specific phonons and electronic valleys to further elucidate the microscopic origins of AMR. Our results highlight potential pathways to modify the AMR rate in silicon via strain engineering.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multigrid based First-Principles Molecular Dynamics

MGmol ls a First-Principles Molecular Dynamics code. It relies on the Born-Oppenheimer approximation and models the electronic structure using Density Functional Theory, either LDA or PBE. Norm-conserving pseudopotentials are used to model atomic cores.

Fattebert, Jean-Luc↗