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 343 records · Page 19

The use of strain to tailor electronic thermoelectric transport properties: A first principles study of 2H-phase CuAlO 2

Using first principles calculations, the use of strain to adjust electronic transport and the resultant thermoelectric (TE) properties is discussed using 2H phase CuAlO 2 as a test case. Transparent oxide materials, such as CuAlO 2 , a p-­type transparent conducting oxide (TCO), have recently been studied for high temperature thermoelectric power generators and coolers for waste heat. Given TCO materials with relative ease of fabrication, low cost of materials, and non-­toxicity, the ability to tailor them to specific temperature ranges, power needs, and size requirements, through the use of strain opens an interesting avenue. We find that strain can have a significant effect on these properties, in some cases detrimental and in others beneficial, including the potential for n-­type power factors larger than the highest p-type case. The physical reasons for this behavior are explained in the terms of the thermoelectric transport distribution and the Landauer distribution of modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Type-II Dirac cones and electron-phonon interaction in monolayer biphenylene from first-principles calculations

In this work, we report a first-principles investigation of electronic structure, topological bands, and electron-phonon interactions in metallic biphenylene sheets. Biphenylene is a recently synthesized sp 2 -bonded carbon allotrope. We find coupling of electrons at the Fermi surface to very high frequency carbon-derived phonons, analogous to superconducting MgB 2 . This leads to low-temperature weak coupling superconductivity due to an unusual combination of exceptionally large logarithmically averaged phonon frequency ω log =1369 K and moderate electron-phonon coupling. The electronic structure shows a two-band Fermi surface dominated by C p z orbitals and a pair of type-II tilted Dirac cones along the Γ-Y line at the Brillouin zone boundary. Berry curvature and edge-state calculations show that monolayer biphenylene is a two-dimensional $\mathbb{Z}$ 2 topological material. Thus, monolayer biphenylene is predicted to be a topological superconductor based on C p orbitals and high-frequency phonons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Thermodynamically informed priors for uncertainty propagation in first-principles statistical mechanics

Here, this work demonstrates how first-principles statistical mechanics approaches within a Bayesian framework can quantify and propagate uncertainties to downstream thermodynamic calculations. To address the issue of Bayesian prior selection, knowledge of 0 K ground states in the material system of interest is incorporated into the prior. The effectiveness of this framework is shown by creating a phase diagram for the fcc zirconium nitride system, including confidence intervals on order-disorder transition temperatures.

Bayesian methods↗

First-Principles Calculations, Experimental Study, and Thermodynamic Modeling of the Al-Co-Cr System

The phase relations and thermodynamic properties of the condensed Al-Co-Cr ternary alloy system are investigated using first-principles calculations based on density functional theory (DFT) and phase-equilibria experiments that led to X-ray diffraction (XRD) and electron probe micro-analysis (EPMA) measurements. A thermodynamic description is developed by means of the calculations of phase diagrams (CALPHAD) method using experimental and computational data from the present work and the literature. Emphasis is placed on modeling the bcc-A2, B2, fcc-γ, and tetragonal-σ phases in the temperature range of 1173 to 1623 K. Liquid, bcc-A2 and fcc-γ phases are modeled using substitutional solution descriptions. First-principles special quasirandom structures (SQS) calculations predict a large bcc-A2 (disordered)/B2 (ordered) miscibility gap, in agreement with experiments. A partitioning model is then used for the A2/B2 phase to effectively describe the order-disorder transitions. The critically assessed thermodynamic description describes all phase equilibria data well. A2/B2 transitions are also shown to agree well with previous experimental findings.

36 MATERIALS SCIENCE↗

Multilayered vacuum-arc nanocomposite TiN/ZrN coatings before and after annealing: Structure, properties, first-principles calculations

Highlights: • Multilayer TiN/ZrN films were deposited using sequential vacuum-arc deposition of Ti and Zr targets in nitrogen atmosphere; • Several samples were annealed in air at the temperature 700 °C; • First-principles calculations of TiN(111), ZrN(111) structures and TiN(111)/ZrN(111) multilayer were carried out; • All deposited samples were highly polycrystalline with quite large 20–25 nm crystals; • The nanohardness and elastic modulus of non-annealed coatings reached 42 GPa and 348 GPa, respectively; • Deposited coatings demonstrate good wear and oxidation resistance; - Abstract: Nanoscale multilayered TiN/ZrN films were deposited using sequential vacuum-arc deposition of Ti and Zr targets in a nitrogen atmosphere. Studies of film's properties were carried out using various modern methods of analysis, such as XRD, STEM, HRTEM, SIMS combined with results of nanoindentation and tribological tests. To interpret the mechanical properties of the deposited multilayer films first-principles calculations of TiN(111), ZrN(111) structures and TiN(111)/ZrN(111) multilayer were carried out. To study the influence of thermal annealing, several samples were annealed in air at the temperature 700 °C. All deposited samples were highly polycrystalline with quite large 20–25 nm crystals. The crystalline planes were very ordinated and demonstrated an excellent coordinated growth. The nanohardness and elastic modulus of non-annealed coatings reached 42 GPa and 348 GPa, respectively. Annealing in air at the temperature 700 °C led to partial oxidation of the multilayered coatings, however hardness of the non-oxidized part of the coatings remained as high, as for initial coatings. All deposited coatings demonstrate good wear resistance.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

First-principles calculation of principal Hugoniot and K-shell X-ray absorption spectra for warm dense KCl

Principal Hugoniot and K-shell X-ray absorption spectra of warm dense KCl are calculated using the first-principles molecular dynamics (FPMD) method. Evolution of electronic structures as well as the influence of the approximate description of ionization on pressure (caused by the underestimation of the energy gap between conduction bands and valence bands) in the first-principles method are illustrated by the calculation. It is shown that approximate description of ionization in FPMD has small influence on Hugoniot pressure due to mutual compensation of electronic kinetic pressure and virial pressure. The calculation of X-ray absorption spectra shows that the band gap of KCl persists after the pressure ionization of the 3p electrons of Cl and K taking place at lower energy, which provides a detailed understanding to the evolution of electronic structures of warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First-principles study of the T center in silicon

The T center in silicon is a well-known carbon-based color center that has been recently considered for quantum technology applications. Using first-principles computations, we show that the excited state is formed by a defect-bound exciton made of a localized defect state occupied by an electron to which a hole is bound. The localized state is of strong carbon p character and reminiscent of the localization of the unpaired electron in the ethyl radical molecule. The radiative lifetime for the defect-bound exciton is calculated to be on the order of microseconds, much longer than for other quantum defects such as the nitrogen vacancy center in diamond and in agreement with experiments. The longer lifetime is associated with the small transition dipole moment as a result of the very different nature of the localized and delocalized states forming the defect-bound exciton. Finally, we use first-principles calculations to assess the stability of the T center. Here, we find the T center to be stable against decomposition into simpler defects when keeping the stoichiometry fixed. However, we identify that the T center is easily prone to (de)hydrogenation and so requires very precise annealing conditions (temperature and atmosphere) to be efficiently formed.

36 MATERIALS SCIENCE↗

Long-range quadrupole electron-phonon interaction from first principles

Lattice vibrations in materials induce perturbations on the electron dynamics in the form of long-range (dipole and quadrupole) and short-range (octopole and higher) potentials. The dipole Fröhlich term can be included in current first-principles electron-phonon (e-ph) calculations and is present only in polar materials. The quadrupole e-ph interaction is present in both polar and nonpolar materials, but currently it cannot be computed from first principles. In this work, we show an approach to compute the quadrupole e-ph interaction and include it in ab initio calculations of e-ph matrix elements. The accuracy of the approach is demonstrated by comparing with direct density functional perturbation theory calculations. We apply our method to silicon as a case of a nonpolar semiconductor and tetragonal PbTiO 3 as a case of a polar piezoelectric material. In both materials we find that the quadrupole term strongly impacts the e-ph matrix elements. Analysis of e-ph interactions for different phonon modes reveals that the quadrupole term mainly affects optical modes in silicon and acoustic modes in PbTiO 3 , although the quadrupole term is needed for all modes to achieve quantitative accuracy. The effect of the quadrupole e-ph interaction on electron scattering processes and transport is shown to be important. Our approach enables accurate studies of e-ph interactions in broad classes of nonpolar, polar, and piezoelectric materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling Electron Injection at Semiconductor–Molecule Interfaces using First-Principles Dynamics Simulation: Effects of Nonadiabatic Coupling, Self-energy, and Surface Models

Excited electron transfer across semiconductor–molecule heterogeneous interfaces is central to various future electronic and optoelectronic devices. At the same time, first-principles modeling of such dynamical processes remains as a great challenge in theoretical chemistry and condensed matter physics for developing better understanding at the molecular scale. Excited electron transfer from a molecule to semiconductor surface is a particularly difficult case to model accurately, because the initial state of such an electron injection process often lies deep within the dense manifold of the conduction band states in the semiconductor. Nonadiabatic couplings and energy level alignments at such interfaces as well as the finite size error of the surface model all play important roles in numerical modeling of electron injection via first-principles theory. Here, using representative interfaces between a well-defined hydrogen-terminated Si(111) surface and series of covalently adsorbed conjugated molecules, we investigate the extent to which these theoretical and numerical considerations influence the description of electron injection at the semiconductor–molecule interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Double-Atom Catalysts Featuring Inverse Sandwich Structure for CO 2 Reduction Reaction: A Synergetic First-Principles and Machine Learning Investigation

Electrocatalytic CO 2 reduction reactions (CO 2 RR) based on scalable and highly efficient catalysis provide an attractive strategy for reducing CO 2 emissions. Here in this work, we combined first-principles density functional theory (DFT) and machine learning (ML) to comprehensively explore the potential of double-atom catalysts (DACs) featuring an inverse sandwich structure anchored on defective graphene (gra) to catalyze CO 2 RR to generate C 1 products. We started with five homonuclear M 2 ⊥gra (M = Co, Ni, Rh, Ir, and Pt), followed by 127 heteronuclear MM'⊥gra (M = Co, Ni, Rh, Ir, and Pt, M' = Sc–Au). Stable DACs were screened by evaluating their binding energy, formation energy, and dissolution potential of metal atoms, as well as conducting first-principles molecular dynamics simulations with and without solvent water molecules. Based on DFT calculations, Rh 2 ⊥gra DAC was found to outperform the other four homonuclear DACs and the Rh-based single- and double-atom catalysts of noninverse sandwich structures. Out of the 127 heteronuclear DACs, 14 were found to be stable and have good catalytic performance. An ML approach was adopted to correlate key factors with the activity and stability of the DACs, including the sum of radii of metal and ligand atoms (d M–M' , d M–C , and d M'–C ), the sum and difference of electronegativity of two metal atoms (P M + P M' , P M – P M '), the sum and difference of first ionization energy of two metal atoms (I M + I M' , I M – I M '), the sum and difference of electron affinity of two metal atoms (A M + A M' , A M – A M '), and the number of d-electrons of the two metal atoms (Nd). The obtained ML models were further used to predict 154 potential electrocatalysts out of 784 possible DACs featuring the same inverse sandwich configuration. Overall, this work not only identified promising CO 2 RR DACs featuring the reported inverse sandwich structure but also provided insights into key atomic characteristics associated with high CO 2 RR activity.

30 DIRECT ENERGY CONVERSION↗

First-principles equations of state and structures of liquid metals in multi-megabar conditions

Liquid metals at extreme pressures and temperatures are widely interested in the high-pressure community. Based on density functional theory molecular dynamics, we conduct first-principles investigations on the equation of state ( EOS) and structures of four metals (Cu, Fe, Pb, and Sn) at 1.5–5 megabar conditions and 5 × 10 3 –4 × 10 4 K. Our first-principles EOS data enable evaluating the performance of four EOS models in predicting Hugoniot densities and temperatures of the four systems. We find the melting temperature of Cu is 1000–2000 K higher and shows a similar Clapeyron slope, in comparison to those of Fe. Our structure, coordination number, and diffusivity analysis indicates all the four liquid metals form similar simple close-packed structures. Our results set theoretical benchmarks for EOS development and structures of metals in their liquid states and under dynamic compression.

36 MATERIALS SCIENCE↗

First-Principles Quantum Transport Modeling of Spin-Transfer and Spin-Orbit Torques in Magnetic Multilayers

A unified approach for computing (i) spin-transfer torque in magnetic trilayers like spin valve and magnetic tunnel junction, where injected charge current flows perpendicularly to interfaces, and (ii) spin-orbit torque in magnetic bilayers of the type ferromagnet/spin-orbit-coupled material, where injected charge current flows parallel to the interface, was reviewed. The experimentally explored and technologically relevant spin-orbit-coupled materials include 5d heavy metals, topological insulators, Weyl semimetals, and transition metal dichalcogenides. This approach requires to construct the torque operator for a given Hamiltonian of the device and the steady-state nonequilibrium density matrix, where the latter is expressed in terms of the nonequilibrium Green’s functions and split into three contributions. Tracing these contributions with the torque operator automatically yields field-like and damping-like components of spin-transfer torque or spin-orbit torque vector, which is particularly advantageous for spin-orbit torque where the direction of these components depends on the unknown-in-advance orientation of the current-driven nonequilibrium spin density in the presence of spin-orbit coupling. Illustrative examples are provided by computing spin-transfer torque in a one-dimensional toy model of a magnetic tunnel junction and realistic Co/Cu/Co spin valve, both of which are described by first-principles Hamiltonians obtained from noncollinear density functional theory calculations, as well as by computing spin-orbit torque in a ferromagnetic layer described by a tight-binding Hamiltonian which includes spin-orbit proximity effect within ferromagnetic monolayers assumed to be generated by the adjacent monolayer transition metal dichalcogenide. In addition, it is shown here how spin-orbit proximity effect, quantified by computing (via first-principles retarded Green’s function) spectral functions and spin textures on monolayers of realistic ferromagnetic material like Co in contact with heavy metal or monolayer transition metal dichalcogenide, can be tailored to enhance the magnitude of spin-orbit torque. Furthermore, errors made in the calculation of spin-transfer torque are quantified when using Hamiltonian from collinear density functional theory, with rigidly rotated magnetic moments to create noncollinear magnetization configurations, instead of proper (but computationally more expensive) self-consistent Hamiltonian obtained from noncollinear density functional theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Thermodynamic modeling of the La-Te system aided by first-principles calculations

A complete thermodynamic description of the La-Te binary system is developed by means of CALculation of PHAse Diagram (CALPHAD) method in combination with available experimental data in the literature and the present first-principles calculations based on density functional theory. The intermetallic phases with homogeneity ranges, La3-xTe4 and LaTe2-x, are modeled using a two-sublattice (La,Va)3(Te)4 model and a three-sublattice (La)1(Te)1(Te,Va)1 model based on their structure features, respectively. The intermetallic phases, LaTe and LaTe3, are treated as stoichiometric compounds. The thermodynamic properties of the intermetallic compounds and their corresponding end members at finite temperatures are predicted using first-principles quasi-harmonic approach. The associate solution model is used to describe the short-range ordering behavior of the liquid phase. The calculated phase diagram agrees well with the available phase equilibrium data in the literature.

Chemistry↗

First-principles equation of state of polystyrene and its effect on inertial confinement fusion implosions

Obtaining an accurate equation of state (EOS) of polystyrene (CH) is crucial to reliably design inertial confinement fusion (ICF) capsules using CH/CH-based ablators. Thus, with first-principles calculations, we have investigated the extended EOS of CH over a wide range of plasma conditions (ρ = 0.1 to 100 g/cm 3 and T = 1,000 to 4,000,000 K). When compared with the widely used SESAME -EOS table, the first-principles equation of state (FPEOS) of CH has shown significant differences in the low-temperature regime, in which strong coupling and electron degeneracy play an essential role in determining plasma properties. Hydrodynamic simulations of cryogenic target implosions on OMEGA using the FPEOS table of CH have predicted ~5% reduction in implosion velocity and ~30% decrease in neutron yield in comparison with the usual SESAME simulations. This is attributed to the ~10% lower mass ablation rate of CH predicted by FPEOS. Simulations using CH-FPEOS show better agreement with measurements of Hugoniot temperature and scattered lights from ICF implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First principles modeling of the metal-electrolyte interface: A novel approach to the study of the electrochemical interface

The research objective of this proposal is the computational modeling of the metal-electrolyte interface purely from first principles. The accurate calculation of the electrostatic potential at electrically biased metal-electrolyte interfaces is a current challenge for periodic “ab-initio” simulations. It is also an essential requisite for predicting the correspondence between the macroscopic voltage and the microscopic interfacial charge distribution in electrochemical fuel cells. This interfacial charge distribution is the result of the chemical bonding between solute and metal atoms, and therefore cannot be accurately calculated with the use of semi-empirical classical force fields. The project aims to study in detail the structure and dynamics of aqueous electrolytes at metallic interfaces taking into account the effect of the electrode potential. Another side of the project is to produce an accurate method to simulate the water/metal interface. While both experimental and theoretical surface scientists have made a lot of progress on the understanding and characterization of both atomistic structures and reactions at the solid/vacuum interface, the theoretical description of electrochemical interfaces is still lacking behind. A reason for this is that a complete and accurate first principles description of both the liquid and the metal interfaces is still computationally too expensive and complex, since their characteristics are governed by the explicit atomic and electronic structure built at the interface as a response to environmental conditions. This project will characterize in detail how different theoretical levels of modeling describer the metal/water interface. In particular the role of van der Waals interactions will be carefully analyzed and prescriptions to perform accurate simulations will be produced.

36 MATERIALS SCIENCE↗

Structural stability and electronic properties of β-tetragonal boron: A first-principles study

It is known that elemental boron has five polymorphs: α- and β-rhombohedral, α- and β-tetragonal, and the high-pressure γ phase. β-tetragonal (β-t) boron was first discovered in 1960, but there have been only a few studies since then. We have thoroughly investigated, using first-principles calculations, the atomic and electronic structures of β-t boron, the details of which were not known previously. The difficulty of calculation arises from the fact that β-t boron has a large unit cell that contains between 184 and 196 atoms, with 12 partially-occupied interstitial sites. This makes the number of configurations of interstitial atoms too great to calculate them all. By introducing assumptions based on symmetry and preliminary calculations, the number of configurations to calculate can be greatly reduced. It was eventually found that β-t boron has the lowest total energy, with 192 atoms (8 interstitial atoms) in an orthorhombic lattice. The total energy per atom was between those of α- and β-rhombohedral boron. Another tetragonal structure with 192 atoms was found to have a very close energy. The valence bands were fully filled and the gaps were about 1.16 to 1.54 eV, making it comparable to that of β-rhombohedral boron. - Graphical abstract: Electronic density distribution for the lowest-energy configuration (N=192) viewed from the 〈1 0 0〉 direction. Left: isosurface (yellow) at d=0.09 electrons/a.u.{sup 3} Right: isosurface (orange) at d=0.12 electrons/a.u.{sup 3}. - Highlights: • β-tetragonal boron was thoroughly investigated using first-principles calculations. • The lowest energy structure contains 192 atoms in an orthorhombic lattice. • Another tetragonal structure with 192 atoms has a very close energy. • The total energy per atom is between those of α- and β-rhombohedral boron. • The band gap of the lowest energy structure is about 1.16 to 1.54 eV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Large-Scale First-Principles Molecular Dynamics Simulations with Electrostatic Embedding: Application to Acetylcholinesterase Catalysis

Enzymes are complicated solvated systems that typically require many atoms to simulate their function with any degree of accuracy. We have recently developed numerical techniques for large scale First-Principles molecular dynamics simulations and applied them to study the enzymatic reaction catalyzed by acetylcholinesterase. We carried out Density functional theory calculations for a quantum mechanical (QM) sub- system consisting of 612 atoms with an O(N) complexity finite-difference approach. The QM sub-system is embedded inside an external potential field representing the electrostatic effect due to the environment. We obtained finite temperature sampling by First-Principles molecular dynamics for the acylation reaction of acetylcholine catalyzed by acetylcholinesterase. Our calculations shows two energies barriers along the reaction coordinate for the enzyme catalyzed acylation of acetylcholine. In conclusion, the second barrier (8.5 kcal/mole) is rate-limiting for the acylation reaction and in good agreement with experiment.

59 BASIC BIOLOGICAL SCIENCES↗

Effect of Ni Dopant on Furan Activation over Mo 2 C Surface: Insights from First-Principles-Based Microkinetic Modeling

Here, the role of Ni dopant in a molybdenum carbide (Mo 2 C) catalyst for butadiene and propyne production from furan is investigated by using a Mo 2 C(001) surface model and first-principles-based microkinetic modeling (MKM). A reaction mechanism including 18 reaction steps is generated to compute the turnover frequency (TOF) of butadiene and propyne from furan gas. The butadiene production is dominant with the TOF being 3-6 orders of magnitude higher than the TOF of propyne on a bare Mo 2 C surface due to an active C-O bond scission reaction compared to C-C bond cleavage reaction. The exothermic binding of ring-opened furan (CHCHCHCHO*) inhibits the activity of a bare Mo 2 C surface for butadiene and propyne production. Calculations predict that Mo 2 C surface modification with a Ni dopant energetically destabilizes the adsorption of CHCHCHCHO*, which facilitates the formation of butadiene and propyne. In the presence of Ni dopants on a Mo 2 C(001) surface, the selectivity toward the butadiene formation is decreased compared to propyne due to relatively favorable C-C bond scission as opposed to the C-O bond scission. This first-principles-based study provides mechanistic and kinetic insights into the role of dopants in cost-efficient Mo 2 C catalysts during the vapor phase conversion of furan to gaseous products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗