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 235 records · Page 13

Calculating electronic stopping power in materials from first principles

From the early models of electronic stopping power to the current first principles simulations, the techniques evolved to increase the range of validity and to reduce empiricism. Thanks to a combination of theoretical advances provided by Time Dependent Density Functional Theory and the development of numerical codes, it is currently possible to predict electronic stopping power for realistic materials by performing direct simulations of the electron excitation processes beyond linear response, and including electronic band structure effects. Electronic stopping power is an important quantity used to predict and understand the effects of particle radiation in matter. First principles calculations of electronic stopping power can be applied to any atomistic system, solids, liquids and alloys. This review aims here to help graduate level students and researchers immerse themselves into state-of-the-art techniques to computationally model and calculate electronic stopping power.

36 MATERIALS SCIENCE↗

Improved first-principles equation-of-state table of deuterium for high-energy-density applications

Here, we present a first-principles equation-of-state (EOS) table of deuterium aimed at improving the previously established first-principles equation-of-state table (FPEOS). The EOS table presented here, referred to as iFPEOS, introduces (1) a universal density functional theory (DFT) treatment of all density and temperature conditions, (2) a fully consistent treatment of exchange-correlation (XC) thermal effects across the entire range of temperatures covered, and (3) quantum treatment of ions. Based on ab initio molecular dynamics driven by thermal density functional theory, iFPEOS includes density points in the range 1×10 –3 $\leqslant$ ρ $\leqslant$ 1.6×10 3 g/cm 3 and temperature points in the range 800 K $\leqslant$ T $\leqslant$ 256 MK, thus covering the challenging warm dense matter (WDM) regime. For an improved description of the electronic structure, iFPEOS employs an advanced free-energy XC density functional with explicit temperature dependence, which is at the metageneralized gradient approximation level of DFT. We use the latest orbital-free free-energy density functional for the high-temperature regime where it shows excellent agreement with standard Mermin-Kohn-Sham DFT. For quantum treatment of ions we use path-integral molecular dynamics in order to take into account nuclear quantum effects. Results are compared to other EOS models and most recent experimental measurements of deuterium properties such as the molecular-toatomic fluid transition, the principal and reshock Hugoniot, and sound speed. We find that iFPEOS provides an improved agreement with experimental data compared to other first-principles EOS models in the WDM regime for pressures up to 200 GPa and temperatures up to 60 000 K. For higher pressures and temperatures, however, iFPEOS is in agreement with other models in predicting lower compressibility and higher sound speed along the Hugoniot, compared to experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Data on the comprehensive first-principles diffusion study of the aluminum-magnesium system

First-principles calculation of diffusion coefficients between Mg and Al is investigated comprehensively using density functional theory (DFT). The effect of different uncertainty sources arising from first principles calculations has been investigated systematically. These sources include the diffusion model, energetic, entropic and attempt frequency calculations. Variation in self and impurity diffusion coefficients of Mg and Al in stable phases are quantified using different DFT settings and compared with the experiments. Using the optimal DFT settings, diffusion coefficients in metastable phases of Al and Mg are predicted. The dataset refers to “An integrated experimental and computational study of diffusion and atomic mobility of the aluminum-magnesium system” [1].

36 MATERIALS SCIENCE↗

First principles prediction of the Al-Li phase diagram including configurational and vibrational entropic contributions

The whole Al-Li phase diagram is predicted from first principles calculations and statistical mechanics including the effect of configurational and vibrational entropy. The formation enthalpy of different configurations at different temperatures was accurately predicted by means of cluster expansions that were fitted from first principles calculations. The vibrational entropic contribution of each configuration was determined from the bond length vs. bond stiffness relationships for each type of bond and the Gibbs free energy of the different phases was obtained as a function of temperature from Monte Carlo simulations. The predicted phase diagram was in excellent agreement with the currently accepted experimental one in terms of the stable (AlLi, Al 2 Li 3 , AlLi 2 , Al 4 Li 9 ) and metastable (Al 3 Li) phases, of the phase boundaries between them and of the maximum stability temperature of line compounds. In addition, it provided accurate information about the gap between Al 3 Li and AlLi solvus lines. Finally, the influence of the vibrational entropy on the correct prediction of the phase diagram is discussed. Overall, the methodology shows that accurate phase diagrams of alloys of technological interest can be predicted from first principles calculations.

36 MATERIALS SCIENCE↗

First-principles calculations of thermal, electrical, and thermoelectric transport properties of semiconductors

The transport properties of semiconductors are key to the performance of many solid-state devices (transistors, data storage, thermoelectric cooling and power generation devices, etc). An understanding of the transport details can lead to material designs with better performances. In recent years simulation tools based on first-principles calculations have been greatly improved, being able to obtain the fundamental ground-state properties of materials (such as band structure and phonon dispersion) accurately. Accordingly, methods have been developed to calculate the transport properties based on an ab initio approach. We focus on the thermal, electrical, and thermoelectric transport properties of semiconductors, which represent the basic transport characteristics of the two degrees of freedom in solids—electronic and lattice degrees of freedom. Starting from the coupled electron-phonon Boltzmann transport equations, we illustrate different scattering mechanisms that change the transport features and review the first-principles approaches that solve the transport equations. We then present the first-principles results on the thermal and electrical transport properties of semiconductors. The discussions are grouped based on different scattering mechanisms including phonon-phonon scattering, phonon scattering by equilibrium electrons, carrier scattering by equilibrium phonons, carrier scattering by polar optical phonons, scatterings due to impurities, alloying and doping, and the phonon drag effect. We show how the first-principles methods allow one to investigate transport properties with unprecedented detail and also offer new insights into the electron and phonon transport. The current status of the simulation is mentioned when appropriate and some of the future directions are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Solution-based thermodynamic modeling of the Ni-Al-Mo system using first-principles calculations

A solution-based thermodynamic description of the ternary Ni–Al–Mo system is developed here, incorporating first-principles calculations and reported modeling of the binary Ni–Al, Ni–Mo and Al–Mo systems. To search for the configurations with the lowest energies of the N phase, the Alloy Theoretic Automated Toolkit (ATAT) was employed and combined with VASP. The liquid, bcc and γ-fcc phases are modeled as random atomic solutions, and the γ'-Ni3Al phase is modeled by describing the ordering within the fcc structure using two sublattices, summarized as (Al,Mo,Ni)0.75(Al,Mo,Ni)0.25. Thus, γ-fcc and γ'-Ni3Al are modeled with a single Gibbs free energy function with appropriate treatment of the chemical ordering contribution. In addition, notable improvements are the following: first, the ternary effects of Mo and Al in the B2-NiAl and D0a-Ni3Mo phases, respectively, are considered; second, the N-NiAl8Mo3 phase is described as a solid solution using a three-sublattice model; third, the X-Ni14Al75Mo11 phase is treated as a stoichiometric compound. Model parameters are evaluated using first-principles calculations of zero-Kelvin formation enthalpies and reported experimental data. In comparison with the enthalpies of formation for the compounds ψ-AlMo, θ-Al8Mo3 and B2-NiAl, the first-principles results indicate that the N-NiAl8Mo3 phase, which is stable at high temperatures, decomposes into other phases at low temperature. Resulting phase equilibria are summarized in the form of isothermal sections and liquidus projections. To clearly identify the relationship between the γ-fcc and γ'-Ni3Al phases in the ternary Ni–Al–Mo system, the specific γ-fcc and γ'-Ni3Al phase fields are plotted in x(Al)–x(Mo)–T space for a temperature range 1200–1800 K.

36 MATERIALS SCIENCE↗

Effects of interlayer screening and temperature on dielectric functions of graphene by first-principles

The dielectric functions of few-layer graphene and the related temperature dependence are investigated from the atomic scale using first-principles calculations. Compared with ellipsometry experiments in the spectral range of 190–2500 nm, the normalized optical constants of mono-layer graphene demonstrate good agreement and further validate first-principles calculations. To interpret dielectric function of mono-layer graphene, the electronic band structure and density of states are analyzed. By comparing dielectric functions of mono-, bi-, and tri-layer graphene, it shows that interlayer screening strengthens intraband transition and greatly enhances the absorption peak located around 1 eV. The strengthened optical absorption is intrinsically caused by the increasing electron states near the Fermi level. To investigate temperature effect, the first-principles calculations and lattice dynamics are combined. The lattice vibration enhances parallel optical absorption peak around 1 eV and induces redshift. Moreover, it is observed that the van der Waals force plays a key role in keeping the interlayer distance stable during dynamics simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

First-Principles Petascale Simulations for Predicting Deflagration to Detonation Transition in Hydrogen-Oxygen Mixtures

Hydrogen has emerged as an important fuel across a range of industries as a means of achieving energy independence and to reduce emissions. DDT and the resulting detonation waves in hydrogen-oxygen can have especially catastrophic consequences in a variety of industrial and energy producing settings related to hydrogen. First-principles numerical simulations of flame acceleration and DDT are required for an in-depth understanding of the phenomena and facilitating design of safe hydrogen systems. The goals of this project were (1) to develop first-principles petascale reactive flow Navier-Stokes simulation code for predicting gaseous high-speed combustion and detonation (HSCD) phenomena and (2) demonstrate feasibility of first-principles simulations of rapid flame acceleration and deflagration-to-detonation transition (DDT) in stoichiometric hydrogen-oxygen mixture (2H 2 + O 2 ). The goals of the project have been accomplished. We have developed a novel numerical simulation code, named HSCD, for performing first-principles direct numerical simulations of high-speed hydrogen combustion. We carried out a series of validating numerical simulations of inert and reactive shock reflection experiments in shock tubes. We then performed a pilot numerical simulation of flame acceleration in a long pipe. The simulation showed the transition of the rapidly accelerating flame into a detonation. The DDT simulations were performed using BG/Q Mira at the Argonne National Laboratory, currently the fourth fastest super-computer in the world.

08 HYDROGEN↗

First-principles free energies by hybrid thermodynamic integration for phase equilibria and fission product solubility in molten salts

The overarching goal of this research project is quantitative microscopic modeling of phase diagrams and solubility of actinides and lanthanide-fission products in molten salts for molten salt reactors. Application of molten salts as coolants and liquid fuels for next-generation reactors requires quantitative understanding of phase diagrams, equations of state and transport properties of binary and ternary molten salt mixtures, as well as solubility of fuel, fission products and reactor-wall corrosion products, but these remain beyond the reach of current first-principles computational techniques. We will address this fundamental capability gap using a novel hybrid framework for thermodynamic integration to directly predict absolute free energies of solid and molten salts at any specified conditions of composition, temperature and pressure. This framework will combine a hierarchy of first-principles, molecular dynamics and continuum techniques to deliver first-principles-level accuracy at a fraction of the computational cost of conventional ab initio molecular dynamics approaches. In particular, this technique is expected to provide first-principles accuracy, while off-loading most of the computational work on much less expensive classical molecular dynamics and continuum techniques.

36 MATERIALS SCIENCE↗

First principles modeling of microscopic scintillation mechanisms using a Many-Body Green’s function approach (Final Report)

In 2010, we started the first lifecycle of a project on first principles modeling of microscopic unit processes that occur during scintillation in order to bring about a capability of formulating physics-based reaction-diffusion models that can predict scintillator light yield and explain their non-proportional response to the energy of incident ionizing radiation, with the hope that this will accelerate materials discovery that can mitigate this problem. During the ensuing years we focused on investigating from first principles the elementary quantum mechanical processes leading to polaron formation and exciton self-trapping in scintillators and identify the microscopic pathways leading to their transport. The effort to derive the master equation for carrier and exciton transport in realistic scintillators materials from first principles constitutes a challenge that requires an extensive background in various areas of computational physics. Previous theoretical investigations involved significant approximations and were limited in their scope. Most notably, the standard workhorse of the computational material scientist, density-functional theory, fails in describing the localized nature of self-trapped carriers as well as the electronic structure of rare-earth dopants. Advances in both algorithms and computational resources however, have rendered a comprehensive and accurate approach feasible. This has been amply demonstrated by the results produced by us during the first lifecycle of this project, which impressed upon the reviewers of the Independent Project Review: “The Review Team agreed, and commented that ’amazing and excellent progress had been made to date.’ The researchers tried and did improve the technology to advance the state of the art and show that it was adequate for the project.” In the following, we detail the progress achieved during each year of the project. We conclude by a final section where all our main achievements are summarized.

36 MATERIALS SCIENCE↗

Deriving models for the Kitaev spin-liquid candidate material α - RuCl 3 from first principles

We use the constrained random phase approximation to derive from first principles the Ru-$t_{2_g}$ Wannier-function-based model for the Kitaev spin-liquid candidate material $α$-RuCl 3 . We find the nonlocal Coulomb repulsion to be sizable compared to the local one. In addition we obtain the contribution to the Hamiltonian from the spin-orbit coupling and find it to also contain non-negligible nonlocal terms. We invoke strong-coupling perturbation theory to investigate the influence of these nonlocal elements of the Coulomb repulsion and the spin-orbit coupling on the magnetic interactions. We find that the nonlocal Coulomb repulsions cause a strong enhancement of the magnetic interactions, which deviate from experimental fits reported in the literature. Our results contribute to the understanding and design of quantum spin-liquid materials via first-principles calculations.

36 MATERIALS SCIENCE↗

First-principles residual resistivity using a locally self-consistent multiple scattering method

The locally self-consistent multiple scattering (LSMS) method can perform efficient first-principles calculations of systems with a large number of atoms. Here, in this paper, we combine the Kubo-Greenwood equation with LSMS, enabling us to calculate the first-principles residual resistivity of large systems. This has been implemented in the open-source code LSMS . We apply this method to selected pure elements and binary random alloys. The results compare well with experiment, and with values obtained from a first-principles effective medium technique (the Korringa-Kohn-Rostoker coherent potential approximation). We discuss future applications of this method to complex systems where other methods are not applicable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Spin-state transition induced half metallicity in a cobaltate from first principles

Half metal is a promising spintronic material. Here, we explore, using first principles calculations, a spin-state transition induced half metallicity in a layered cobaltate via a physical or chemical pressure. Our exemplary first principles study shows that the layered cobaltate Sr{sub 2}CoO{sub 3}F would undergo a transition, under a pressure of 5.4 GPa, from a high-spin antiferromagnetic insulator to an intermediate-spin ferromagnetic half-metal. The former phase is associated with a superexchange in a Mott insulator, and the latter one is due to a broad band formation and a kinetic energy gain of the partially occupied e{sub g} orbital. Note that the above transition could also be induced by a chemical pressure via doping in (Sr{sub 1−x}Ca{sub x}){sub 2}CoO{sub 3}F (x > 0.3). This work suggests that a cobaltate would be of a particular interest if stabilized into an intermediate-spin state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Phase Stability for the Pd-Si System. First-Principles, Experiments, and Solution-Based Modeling

Relative stabilities of the compounds in the binary Pd-Si system were assessed using first-principles calculations and experimental methods. Calculations of lattice parameters and enthalpy of formation indicate that Pd 5 Si-μ, Pd 9 Si 2 -α, Pd 3 Si-β, Pd 2 Si-γ, and PdSi-δ are the stable phases at 0 K (-273 °C). X-ray diffraction analyses (XRD) and electron probe microanalysis (EPMA) of the as-solidified and heat-treated samples support the computational findings, except that the PdSi-δ phase was not observed at low temperature. Considering both experimental data and first-principles results, the compounds Pd 5 Si-μ, Pd 9 Si 2 -α, Pd 3 Si-β, and Pd 2 Si-γ are treated as stable phases down to 0 K (-273 °C), while the PdSi-δ is treated as being stable over a limited range, exhibiting a lower bound. Using these findings, a comprehensive solution-based thermodynamic model is formulated for the Pd-Si system, permitting phase diagram calculation. Moreover, the liquid phase is described using a three-species association model and other phases are treated as solid solutions, where a random substitutional model is adopted for Pd-fcc and Si-dia, and a two-sublattice model is employed for Pd 5 Si-μ, Pd 9 Si 2 -α, Pd 3 Si-β, Pd 2 Si-γ, and PdSi-δ. Model parameters are fitted using available experimental data and first-principles data, and the resulting phase diagram is reported over the full range of compositions.

36 MATERIALS SCIENCE↗

Magnetic Behaviors of 3d Transition Metal-Doped Silicane: a First-Principle Study

We perform a first-principle computation on the magnetic properties of 3d transition metal (TM)-doped silicane. All the substituted systems are energetically stable. Robust half-metallic properties are found in Ti-, V- and Mn-substituted silicane. At the dopant concentration of 12.5%, we evaluate that the Curie temperature of Ti- and Mn-substituted silicane is at 340 and 666 K respectively in the mean-field approximation. Consequently, room-temperature ferromagnetism is easily achievable in Ti- and Mn-substituted silicane. Our work demonstrates potential application of Ti and Mn-substituted silicane for room-temperature spintronic devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Structural, Elastic, Electronic, and Magnetic Properties of a New Full-Heusler Alloy Mn{sub 2}MgGe: First-Principles Calculations

We present the first-principles calculations of the structural, elastic, electronic, and magnetic properties for a new full-Heusler alloy Mn{sub 2}MgGe. Both L2{sub 1} and XA structures are considered for both nonmagnetic and ferromagnetic states. The results show that the XA structure in the ferromagnetic state is the energetically most favorable for the full-Heusler alloy, and exhibits ductile behavior, significant anisotropy, and robust half-metallicity. The total spin moment is 2.000 μ{sub B} per formula unit in equilibrium state, which follows the Slater-Pauling rule. The spin-up electrons are metallic, whereas the spin-down bands are semiconductor with a gap of 1.086 eV at the equilibrium lattice constant of 6.066 Å. Half-metallicity is maintained within the lattice constant range from 5.6 to 6.1 Å. Our results indicate that Mn{sub 2}MgGe is an interesting candidate in the area of spintronics.

36 MATERIALS SCIENCE↗

Investigation of the reaction mechanism of the hydrodeoxygenation of propionic acid over a Rh(1 1 1) surface: A first principles study

Microkinetic models based on first principles calculations have been used to study the vapor and liquid phase hydrodeoxygenation of propionic acid on a Rh(111) surface. Calculations suggest that both decarboxylation and decarbonylation do not occur at an appreciable rate under all reaction environments. Propanol and propionaldehyde are the main products on this surface and they are produced at similar rates in both vapor and liquid phase environments. While a condensed phase can shift the reaction rate, the dominant pathways and selectivity towards the various products are hardly affected. At 473 K, the turnover frequency is increased by about a factor of 1.5 in liquid water relative to the gas phase. In liquid 1,4-dioxane, the turnover frequency is also slightly increased relative to the gas phase. Here, given the uncertainty in Rh cavity radius in the liquid phase calculations, computations with different cavity radii have been performed. With larger Rh cavity radius, the promotional effect of the solvents on the turnover frequency becomes more significant, while practically no changes are observed for a smaller cavity radius.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interfacial and Alloying Effects on Activation of Ethanol from First-Principles

Here, we present a first-principles density-functional theory study of ethanol activation at oxide/Rh(111) interface and the alloying effect on mitigating carbon deposition, which are essential to direct ethanol fuel cell (DEFC) anode reaction and steam reforming of ethanol (SRE) reaction. Our calculated results show that charge can transfer from Rh(111) substrate to MO x chain (e.g., MoO 3 and MnO 2 ), or from MO x chain (e.g., MgO, SnO 2 , ZrO 2 , and TiO 2 ) to Rh(111) substrate. The OH-binding strength is increased exponentially with M δ+ charge ranging from 1.4 to 2.2, which renders MnO 2 /Rh(111) and MgO/Rh(111) interfaces weaker OH-binding, and thereby enhanced oxidizing functionality of OH* for promoting ethanol oxidation reaction (EOR) at DEFC anode. For efficient C–C bond breaking, a large number of Rh ensemble sizes are critically needed at the interface of MO x /Rh(111). We found that Rh 1 Au 3 near surface alloy has the weakest C* and CO* binding, followed by Rh 1 Cu 3 and Rh 1 Pd 3 near surface alloys, while Rh 1 Ir 3 and Rh 1 Ru 3 surface alloys have C* and CO* binding strength similar to that of pure Rh metal. The general implication of this study is that by engineering alloyed structure of weakened C* and CO* binding complemented with metal oxides of weakened OH-binding, high-performance DEFC anode or SRE catalysts can be identified.

36 MATERIALS SCIENCE↗