Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional calculation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35

Unconventional Anomalous Hall Effect in Hexagonal Polar Magnet Y 3 Co 8 Sn 4

A rare realization of unconventional anomalous Hall effect (UAHE) both below and above the magnetic transition temperature (T C ) in a hexagonal noncentrosymmetric magnet Y 3 Co 8 Sn 4 , using a combined experimental and ab initio calculations. Occurrence of such UAHE is mainly attributed to the reciprocal ($\mathscr{KS}$) topology (i.e., the presence of topological Weyl points at/near the Fermi level), along with some contribution from the topological magnetic texture (at low temperatures), as inferred from the measured field-dependent ac susceptibility. The effect of UAHE on the measured transport behavior however evolves differently with temperature above and below T C = 53 K, suggesting different physical mechanism responsible in the two phases. A unique planar ferrimagnetic ordering is found to be the most stable state with ab-plane as the easy plane below T C . The simulated net magnetization and the moment per Co atom agrees fairly well with the experimentally measured values. A reasonably large anomalous Hall conductivity (AHC) is also observed in both the phases (above and below T C ) of the present compound, which is again not so ubiquitous. In conclusion, the results underscore the family of R 3 Co 8 Sn 4 (R = rare earth) polar magnets as a compelling backdrop for exploring the synergy of magnetism and topological non-trivial electronic states, pivotal for spintronic applications.

Density functional theory calculations↗

Thermodynamics and kinetics of H adsorption and intercalation for graphene on 6 H -SiC(0001) from first-principles calculations

Previous experimental observations for H intercalation under graphene on SiC surfaces motivate the clarification of configuration stabilities and kinetic processes related to intercalation. From first-principles density-functional-theory calculations, we analyze H adsorption and intercalation for graphene on a 6H-SiC(0001) surface, where the system includes two single-atom-thick graphene layers: the top-layer graphene (TLG) and the underling buffer-layer graphene (BLG) above the terminal Si layer. Our chemical potential analysis shows that in the low-H coverage regime (described by a single H atom within a sufficiently large supercell), intercalation into the gallery between TLG and BLG or into the gallery underneath BLG is more favorable thermodynamically than adsorption on top of TLG. However, intercalation into the gallery between TLG and BLG is most favorable. We obtain energy barriers of about 1.3 and 2.3 eV for a H atom diffusing on and under TLG, respectively. From an additional analysis of the energy landscape in the vicinity of a step on the TLG, we assess how readily one guest H atom on the TLG terrace can directly penetrate the TLG into the gallery between TLG and BLG versus crossing a TLG step to access the gallery. Finally, we also perform density functional theory calculations for higher H coverages revealing a shift in favorability to intercalation of H underneath BLG and characterizing the variation with H coverage in interlayer spacings.

36 MATERIALS SCIENCE↗

Ab initio electron-phonon interactions using atomic orbital wave functions

The interaction among electrons and lattice vibrations determines key physical properties of materials, including their electrical and heat transport, excited electron dynamics, phase transitions, and superconductivity. We present an ab initio method that employs atomic orbital (AO) wave functions to compute the electron-phonon (e-ph) interactions in materials and interpolate the e-ph coupling matrix elements to fine Brillouin zone grids. We illustrate the numerical implementation of such AO-based e-ph calculations, and benchmark them against direct density functional theory calculations and Wannier function (WF) interpolation. The prime advantages of AOs over WFs for e-ph calculations are outlined. Since AOs are fixed basis functions associated with the atoms, they circumvent the need to generate a material-specific localized basis set with a trial-and-error approach, as is needed in WFs. Therefore, AOs are ideal to compute e-ph interactions in chemically and structurally complex materials for which WFs are challenging to generate, and are also promising for high-throughput materials discovery. Although our findings focus on AOs, the formalism we present generalizes e-ph calculations to arbitrary localized basis sets, with WFs recovered as a special case.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Phase stability and property evolution of biphasic Ti-Ni-Sn alloys for use in thermoelectric applications

Thermoelectric properties and phase evolution have been studied in biphasic Ti–Ni–Sn materials containing full-Heusler TiNi 2 Sn embedded within half-Heusler thermoelectric TiNiSn. Materials, prepared by levitation induction melting followed by annealing, were of the nominal starting composition of TiNi 1+x Sn, with x between 0.00 and 0.25. Phases and microstructure were determined using synchrotron X-ray diffraction and optical and electron microscopy. The full-Heusler phase is observed to be semi-coherent with the half-Heusler majority phase. Differential thermal analysis was performed to determine melting temperatures of the end-member compounds. The thermal conductivity is reduced with the introduction of a dispersed, full-Heusler phase within the half-Heusler material. This leads to an increased thermoelectric figure of merit, ZT, from 0.35 for the stoichiometric compound to 0.44 for TiNi 1.15 Sn. Beyond x = 0.15 ZT decreases due to a rise in thermal conductivity. Density functional theory calculations using hybrid functionals were performed to determine band alignments between the half- and full-Heusler compounds, as well as comparative energies of formation. The hybrid functional band structure of TiNiSn is presented as well.

Douglas, Jason E.↗

Hard templating ultrathin polycrystalline hematite nanosheets: effect of nano-dimension on CO 2 to CO conversion via the reverse water-gas shift reaction

Understanding how nano-dimensionality impacts iron oxide based catalysis is central to a wide range of applications. Here, we focus on hematite nanosheets, nanowires and nanoparticles as applied to catalyze the reverse water gas shift (RWGS) probe reaction. We introduce a novel approach to synthesize ultrathin (4–7 nm) hematite nanosheets using copper oxide nanosheets as a hard template and propose a reaction mechanism based on density functional theory (DFT) calculations. Hematite nanowires and nanoparticles were also synthesized and characterized. H2 temperature programmed reduction (H2-TPR) and RWGS reactions were performed to glean insights into the mechanism of CO2 conversion to CO over the iron oxide nanomaterials and were compared to H2 binding energy calculations based on density functional theory. While the nanosheets did exhibit high CO2 conversion, 28% at 510 °C, we found that the iron oxide nanowires had the highest CO2 conversion, reaching 50% at 750 °C under atmospheric pressure. No products besides CO and H2O were detected.

Fishman, Zachary S.↗

d -Electron Heavy Fermion Behavior in a Near-Room-Temperature Polar Metallic Ferrimagnet: A Case of Mn 5 SiC

Polycrystalline Mn 5 SiC was synthesized by using a high-temperature solid-state method. Mn 5 SiC adopts a polar space group (Cmc2 1 ) with six crystallographic Mn sites confirmed by X-ray and neutron diffraction, transmission electron microscopy, and second harmonic generation experiments. The complex crystal structure features edge-sharing trigonal prisms and icosahedra, as well as face/edge-sharing pentagonal prisms. Magnetic measurements indicate ferrimagnetic ordering with a transition temperature of 284 K. The ferrimagnetic structure (magnetic space group Cm’c’2 1 ) was further identified by powder neutron diffraction, where collinear Mn spins align along the crystallographic c-axis. The refined magnetic moment for each crystallographic Mn site at 4 K is 1.8(2), −2.42(9), −1.72(8), 0.51(6), 0.50(4), and 1.7(2) μ B . Density functional theory calculations confirm both the metallic behavior and the ferrimagnetic structure observed experimentally and further provide insight into the observed Mn moment dependence across crystallographic sites. The resistivity and specific heat measurements and density functional theory calculations reveal a substantially large Kadowaki–Woods ratio of 5 × 10 –5 μΩ·cm/(mJ/mol) 2 and a many-body renormalization factor of 5.5, indicating the unusual heavy Fermion behavior in such an itinerant magnetic metal.

Messegee, Zachary T. [George Mason University, Fai↗

CHARGED POINT DEFECTS AND DEFECT COMPLEXES IN PHOSPHORENE AND STRUCTURAL PREDICTION OF TAU-11, A DENSITY FUNCTIONAL THEORY APPROACH

With the advent of powerful computers, first principle calculations have become a robust and reliable tool to complement experiments in the study of materials. First principle calculations such as density functional theory can be used to guide experiments to save time and effort when it comes to new and advanced materials discovery. Study of charge defects in materials and crystal structure prediction are such areas in which density functional theory calculations facilitate the prediction of defect formation energy, charge transition levels and elemental site occupations. Experimentally, studying charge defects requires inferring results from various techniques. Likewise, determining site occupancies for similar elements using experimental techniques is challenging and sometimes impossible. In this work, we study the charged intrinsic and extrinsic defects and defect complexes in phosphorene and determine the site occupancies of Al, Fe and Si in the τ11 phase using density functional theory. In the defect work, we calculate the defect formation energy, charge transition level, binding energies and Stokes shift for vacancy, dopant substitution and dopant-vacancy defect complexes in phosphorene. We found that vacancy defect in phosphorene becomes negatively charged in n-type doping and may passivate the dopants and reduce carrier concentration and mobility. For non-metal dopants in phosphorene, we predict that O, S and Si prefer to form dopant-vacancy complexes removing the vacancy defect states from the band gap. Mn dopant-vacancy defect complex exhibits possibility of switching between two magnetic spin states. Lastly, using density functional theory, we complement the results from neutron powder diffraction to determine the site occupancies of Al and Si in τ11 phase which was used to determine the alloying element to develop a high temperature and low density alloy.

36 MATERIALS SCIENCE↗

Phonon Olympics: Phonon property and lattice thermal conductivity benchmarking from open-source packages

Three widely used open-source packages for determining phonon properties and lattice thermal conductivities (ALAMODE, phono3py, and ShengBTE) are benchmarked by teams of expert users and the package developers. The phonons for Ge, RbBr, monolayer MoSe 2 , and AlN are modeled at zero temperature, and they scatter through three-phonon and phonon-isotope processes, with thermal conductivities obtained from the linearized Peierls–Boltzmann transport equation with input from density functional theory calculations. Over a wide range of temperatures, the thermal conductivities calculated by the teams fall within at most ±15% of their mean values for each of the four materials. The phonon frequencies, obtained from the harmonic force constants, do not show large differences between the calculations, indicating that the modal heat capacities and group velocities are not responsible for the thermal conductivity variations. It is the lifetimes associated with three-phonon scattering, obtained from the cubic force constants, that drive the variations. The many decisions required to calculate the cubic force constants (e.g., supercell size, atomic displacement, neighbor cutoff, and application of symmetries) make identification of the precise origin of the thermal conductivity variations challenging. The calculated thermal conductivities do not generally show agreement with experimental measurements, which is attributed to the limitations of the density functional theory calculations. Guidance for the development of best practices is provided, which will help to standardize protocols needed for building thermal conductivity databases. The results provide a baseline for future benchmarking of other packages and more advanced calculations.

McGaughey, Alan J. H. [Carnegie Mellon Univ., Pitt↗

Influence of Alloying Elements and Effect of Stress on Anisotropic Hydrogen Diffusion in Zr-Based Alloys Predicted by Accelerated Kinetic Monte Carlo Simulations

The formation of hydrides degrades the toughness of Zr-based claddings, and thus is one of the primary concerns regarding cladding integrity during fuel operation and used fuel storage. However, the formation mechanism of hydrides, particularly the nucleation mechanism, is yet to be fully understood. To explore the nucleation process, atomic scale studies including density functional theory, molecular dynamics and kinetic Monte Carlo have been used to investigate H diffusion, segregation, and the possible nucleation paths of hydrides. Here, density functional theory calculations are used to parameterize a kinetic Monte Carlo model, which can accurately predict the effective H diffusivity in various Zr-based claddings including Zircaloy, Zirlo and Zr-2.5Nb. Assisted by density functional theory calculations, a recently developed charge-optimized-many-body potential is further improved and used to investigate the possible nucleation paths of hydrides. It’s shown that the formation of face-centered hydrides starts with coherent zeta-phase, shearing of which in the basal plane changes the hcp stacking into fcc and leads to the formation of gamma-hydride nuclei. The shearing involves a negligible barrier and negligible strain accumulation in the matrix by adopting three equivalent shear partials. The results obtained at the atomic scale are consistent with previous theory and experiments, and improves the current understanding on how and how fast hydride forms in Zr based cladding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

First-principles high-pressure unreacted equation of state and heat of formation of crystal 2,6-diamino-3, 5-dinitropyrazine-1-oxide (LLM-105)

We report dispersion-corrected density functional theoretical calculations of the unreacted equation of state (EOS) of crystal 2,6-diamino-3, 5-dinitropyrazine-1-oxide (LLM-105) under hydrostatic compression of up to 45 GPa. Convergence tests for k-points sampling in the Brillouin zone show that a 3 × 1 × 2 mesh is required to reproduce the X-ray crystal structure at ambient conditions, and we confirm our finding with a separate supercell calculation. Our high-pressure EOS yields a bulk modulus of 19.2 GPa, and indicates a tendency towards anisotropic compression along the b lattice vector due to molecular orientations within the lattice. We find that the electronic energy band gap decreases from a semiconductor type of 1.3 eV at 0 GPa to quasi-metallic type of 0.6 eV at 45 GPa. The extensive intermolecular hydrogen bonds involving the oxide (–NO) and dioxide (–NO{sub 2}) interactions with the amine (–NH{sub 2}) group showed enhanced interactions with increasing pressure that should be discernible in the mid IR spectral region. We do not find evidence for structural phase transitions or chemically induced transformations within the pressure range of our study. The gas phase heat of formation is calculated at the G4 level of theory to be 22.48 kcal/mol, while we obtain 25.92 kcal/mol using the ccCA-PS3 method. Density functional theory calculations of the crystal and the gas phases provided an estimate for the heat of sublimation of 32.4 kcal/mol. We thus determine the room-temperature solid heat of formation of LLM-105 to be −9.9 or −6.5 kcal/mol based on the G4 or ccCA-PS3 methods, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Finite-size correction for slab supercell calculations of materials with spontaneous polarization

Abstract The repeated slab approach has become a de facto standard to accurately describe surface properties of materials by density functional theory calculations with periodic boundary conditions. For materials exhibiting spontaneous polarization, we show that the conventional scheme of passivation with pseudo hydrogen is unable to realize a charge-neutral surface. The presence of a net surface charge induces via Gauss’s law a macroscopic electric field through the slab and results in poor size convergence with respect to the thickness of the slab. We propose a modified passivation method that accounts for the effect of spontaneous polarization, describes the correct bulk limits and boosts convergence with respect to slab thickness. The robustness, reliability, and superior convergence of energetics and electronic structure achieved by the proposed method are demonstrated using the example of polar ZnO surfaces.

Chemistry↗

Constraints on magnetism and correlations in RuO 2 from lattice dynamics and Mössbauer spectroscopy

We provide experimental evidence for the absence of a magnetic moment in bulk RuO 2 , a candidate altermagnetic material, by using a combination of Mössbauer spectroscopy, nuclear forward scattering, inelastic X-ray and neutron scattering, and density functional theory calculations. Using complementary Mössbauer and nuclear forward scattering, we determine the Ru magnetic hyperfine splitting to be negligible. Inelastic X-ray and neutron scattering-derived lattice dynamics of RuO 2 are compared to density functional theory calculations of varying flavors. Comparisons among theory with experiments indicate that electronic correlations, rather than magnetic order, are key in describing the lattice dynamics.

Mössbauer spectroscopy↗

Analytic calculations of hyper-Raman spectra from density functional theory hyperpolarizability gradients

We present the first analytic calculations of the geometrical gradients of the first hyperpolarizability tensors at the density-functional theory (DFT) level. We use the analytically calculated hyperpolarizability gradients to explore the importance of electron correlation effects, as described by DFT, on hyper-Raman spectra. In particular, we calculate the hyper-Raman spectra of the all-trans and 11-cis isomers of retinal at the Hartree-Fock (HF) and density-functional levels of theory, also allowing us to explore the sensitivity of the hyper-Raman spectra on the geometrical characteristics of these structurally related molecules. We show that the HF results, using B3LYP-calculated vibrational frequencies and force fields, reproduce the experimental data for all-trans-retinal well, and that electron correlation effects are of minor importance for the hyper-Raman intensities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tuning the reactivity of carbon surfaces with oxygen-containing functional groups

Oxygen-containing carbons are promising supports and metal-free catalysts for many reactions. However, distinguishing the role of various oxygen functional groups and quantifying and tuning each functionality is still difficult. Here we investigate the role of Brønsted acidic oxygen-containing functional groups by synthesizing a diverse library of materials. By combining acid-catalyzed elimination probe chemistry, comprehensive surface characterizations, 15N isotopically labeled acetonitrile adsorption coupled with magic-angle spinning nuclear magnetic resonance, machine learning, and density-functional theory calculations, we demonstrate that phenolic is the main acid site in gas-phase chemistries and unexpectedly carboxylic groups are much less acidic than phenolic groups in the graphitized mesoporous carbon due to electron density delocalization induced by the aromatic rings of graphitic carbon. The methodology can identify acidic sites in oxygenated carbon materials in solid acid catalyst-driven chemistry.

Zhou, Jiahua↗

First-principle natural band alignment of GaN / dilute-As GaNAs alloy

Density functional theory (DFT) calculations with the local density approximation (LDA) functional are employed to investigate the band alignment of dilute-As GaNAs alloys with respect to the GaN alloy. Conduction and valence band positions of dilute-As GaNAs alloy with respect to the GaN alloy on an absolute energy scale are determined from the combination of bulk and surface DFT calculations. The resulting GaN / GaNAs conduction to valence band offset ratio is found as approximately 5:95. Our theoretical finding is in good agreement with experimental observation, indicating the upward movements of valence band at low-As content dilute-As GaNAs are mainly responsible for the drastic reduction of the GaN energy band gap. In addition, type-I band alignment of GaN / GaNAs is suggested as a reasonable approach for future device implementation with dilute-As GaNAs quantum well, and possible type-II quantum well active region can be formed by using InGaN / dilute-As GaNAs heterostructure.

36 MATERIALS SCIENCE↗

A Polar Magnetic and Insulating Double Corundum Oxide: Mn2MnSbO6 with Ordered Mn(II) and Mn(III) Ions

A new magnetic insulator Mn2MnSbO6 with a polar crystal structure and an ordered Mn2+ and Mn3+ arrangement was synthesized under a high pressure of 7.5 GPa and 1300 degrees C. The crystal structure of Mn2MnSbO6, investigated by synchrotron powder X-ray diffraction, was found to be isomorphous with that of Ni3TeO6-type, space group R3. The non-centrosymmetric structure was confirmed by the second-harmonic generation measurements. The X-ray absorption near-edge spectroscopy measurement confirmed the nominal oxidation states of Mn22+Mn3+SbO6. Magnetic measurements indicate that Mn2MnSbO6 orders antiferromagnetically below 44 K and undergoes a field-induced spin-flop transition at 5 K. First-principles calculations indicate an antiferromagnetic ground state with up/up/up/down/down/down (uuuddd) spin configuration of the six crystallographically unique Mn ions in the c-axis doubled magnetic structure. The density functional theory calculations also substantiate the experimentally observed charge ordering of the Mn2+/Mn3+ ions and the insulating behavior due to a bandgap of 0.52 eV. To the best of our knowledge, this is the first double corundum oxide containing Jahn-Teller active Mn3+ ions.

Feng, Hai L.↗

Ti{sub 3}CrCu{sub 4}: A possible 2-D ferromagnetic spin fluctuating system

Ti{sub 3}CrCu{sub 4} is a new ternary compound which crystallizes in the tetragonal Ti{sub 3}Pd{sub 5} structure type. The Cr atoms form square nets in the a-b plane (a = 3.124 Å) which are separated by an unusually large distance c = 11.228 Å along the tetragonal axis, thus forming a -2-D Cr-sublattice. The paramagnetic susceptibility is characterized by a low effective moment, μ{sub eff} = 1.1 μ{sub B}, a low paramagnetic Curie temperature θ{sub P} (below 7 K) and a temperature independent χ{sub 0} = 6.7 x 10{sup −4} emu/mol. The magnetization at 1.8 K increases rapidly with field nearly saturating to 0.2 μ{sub B}/f.u. The zero field heat capacity C/T shows an upturn below 7 K (∼190 mJ/mol K{sup 2} at ∼0.1K) which is suppressed in applied magnetic fields and interpreted as suggesting the presence of spin fluctuations. The resistivity at low temperatures shows non-Fermi liquid behavior. Overall, the experimental data thus reveal an unusual magnetic state in Ti{sub 3}CrCu{sub 4}, which likely has its origin in the layered nature of the Cr sub-lattice and ferromagnetic spin fluctuations. Density functional theoretical calculations reveal a sharp Cr density of states peak just above the Fermi level, indicating the propensity of Ti{sub 3}CrCu{sub 4} to become magnetic.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗