Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT 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 973 records · Page 54

Layer-dependent spin-resolved electronic structure of ferromagnetic triple-layered ruthenate Sr4Ru3O10

High-resolution angle- and spin-resolved photoemission spectroscopy (ARPES) of the triple-layered ruthenate Sr4Ru3O10 reveals features of the electronic structure that extend our understanding of the layered strontium ruthenates. The spectra near the Fermi energy are very different from the nonmagnetic analogues Sr2RuO4 and Sr3Ru2O7 with distinct Fermi surfaces for wide electronlike minority spin bands around the zone center and narrow holelike majority spin Fermi surface contours around the zone corners. The most dramatic results are two narrow spectral peaks ∼30 meV below the Fermi level, a spin-minority holelike band at the Brillouin zone center, and a spin-majority saddle-band van Hove singularity at the zone edge, which exhibits almost 100% spin polarization at low temperature, and a strong temperature dependent coherence-incoherence crossover attributed to Hund metal correlations. Quantitative comparison of the ARPES to spin-polarized density functional theory (DFT) calculations identify the specific antibonding and nonbonding orbital origins of the narrow bands, with a prediction of different spatial localization in the central and outer layers. This is shown to be consistent with experimental ARPES multizone matrix element intensity variations, and implicates outer-layer-specific control of the in-plane metamagnetism. The renormalization of the bands relative to the mean-field DFT, the demonstration of spin-polarized oxygen bands, and of spin-minority and spin-majority band-crossing hybridization provide a more complete picture of the magnetism which displays aspects of both delocalized and local moment behavior.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

FINETUNA: fine-tuning accelerated molecular simulations

Abstract Progress towards the energy breakthroughs needed to combat climate change can be significantly accelerated through the efficient simulation of atomistic systems. However, simulation techniques based on first principles, such as density functional theory (DFT), are limited in their practical use due to their high computational expense. Machine learning approaches have the potential to approximate DFT in a computationally efficient manner, which could dramatically increase the impact of computational simulations on real-world problems. However, they are limited by their accuracy and the cost of generating labeled data. Here, we present an online active learning framework for accelerating the simulation of atomic systems efficiently and accurately by incorporating prior physical information learned by large-scale pre-trained graph neural network models from the Open Catalyst Project. Accelerating these simulations enables useful data to be generated more cheaply, allowing better models to be trained and more atomistic systems to be screened. We also present a method of comparing local optimization techniques on the basis of both their speed and accuracy. Experiments on 30 benchmark adsorbate-catalyst systems show that our method of transfer learning to incorporate prior information from pre-trained models accelerates simulations by reducing the number of DFT calculations by 91%, while meeting an accuracy threshold of 0.02 eV 93% of the time. Finally, we demonstrate a technique for leveraging the interactive functionality built in to Vienna ab initio Simulation Package (VASP) to efficiently compute single point calculations within our online active learning framework without the significant startup costs. This allows VASP to work in tandem with our framework while requiring 75% fewer self-consistent cycles than conventional single point calculations. The online active learning implementation, and examples using the VASP interactive code, are available in the open source FINETUNA package on Github.

97 MATHEMATICS AND COMPUTING↗

CO Adsorption on Au Nanoparticles Grown on Hexagonal Boron Nitride/Rh(111)

CO adsorption on size-controlled Au nanoparticles grown on an $h$-BN/Rh(111) nanomesh surface has been examined to probe their potential catalytic properties. A combination of high-resolution electron energy loss spectroscopy (HREELS), temperature-programmed desorption (TPD), and density functional theory (DFT) calculations demonstrate that the CO adsorption strength depends heavily on the Au deposition coverage and particle morphology. Particles resulting from low Au coverages deposited at the liquid nitrogen temperature exhibit significantly enhanced CO binding relative to bulk crystalline Au. The resulting CO TPD spectra, and the significantly red-shifted C–O stretching frequency and negative charging of the Au nanoparticles as evidenced by HREELS and DFT, all correspond to those reported for catalytically active Au nanoparticles grown on reactive metal oxides, even though the h-BN/Rh(111) surface is free of carbon, oxygen, or defects. DFT modeling further suggests that the enhanced CO adsorption occurs at highly undercoordinated Au atoms on the perimeters of the nanoparticles. As the Au coverage is increased, the CO adsorption energy and C–O stretching frequency converge toward values associated with bulk Au. Annealing to 600 K results in bulk-like CO adsorption characteristics at any Au coverage. These results suggest that the $h$-BN/Rh(111) surface represents a potential platform for evaluating the role of Au vs support in low-temperature CO oxidation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Kinetic control in the synthesis of metastable polymorphs: Bixbyite-to-Rh{sub 2}O{sub 3}(II)-to-corundum transition in In{sub 2}O{sub 3}

An example for kinetic control of a solid-state phase transformation, in which the system evolves via the path with the lowest activation barrier rather than ending in the thermodynamically most favorable state, has been demonstrated. As a case study, the phase transitions of indium sesquioxide (In{sub 2}O{sub 3}) have been guided by theoretical calculations and followed in situ under high-pressure high-temperature conditions in multi-anvil assemblies. The corundum-type rh-In{sub 2}O{sub 3} has been synthesized from stable bixbyite-type c-In{sub 2}O{sub 3} in two steps: first generating orthorhombic Rh{sub 2}O{sub 3}-II-type o′-In{sub 2}O{sub 3} which is thermodynamically stable at 8.5 GPa/850 °C and, thereafter, exploiting the preferred kinetics in the subsequent transformation to the rh-In{sub 2}O{sub 3} during decompression. This synthesis strategy of rh-In{sub 2}O{sub 3} was confirmed ex situ in a toroid-type high-pressure apparatus at 8 GPa and 1100 °C. The pressure–temperature phase diagrams have been constructed and the stability fields of In{sub 2}O{sub 3} polymorphs and the crystallographic relationship between them have been discussed. - Graphical abstract: In situ energy-dispersive XRD patterns in multi-anvil assemblies show the sequence of phase transition c-In{sub 2}O{sub 3}→o′-In{sub 2}O{sub 3}→rh-In{sub 2}O{sub 3} under particular pressure and temperature conditions. The tick marks refer to the calculated Bragg positions of bixbyite-type (c-In{sub 2}O{sub 3}), Rh{sub 2}O{sub 3}-II-type (o–-In2O{sub 3}) and corundum-type (rh-In{sub 2}O{sub 3}). - Highlights: • The solid-state synthesis methods can be employed for obtaining metastable phases. • The phase transition of In{sub 2}O{sub 3} was guided by DFT calculations. • The phase transition of In{sub 2}O{sub 3} was followed in situ under HP–HT conditions. • Orthorhombic o′-In{sub 2}O{sub 3} polymorph was synthesized from c-In{sub 2}O{sub 3} at 8.5 GPa/850 °C. • Metastable rh-In{sub 2}O{sub 3} was obtained from o′-In{sub 2}O{sub 3} at 5.5 GPa during decompression.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impact of Hydration and Sulfonation on the Morphology and Ionic Conductivity of Sulfonated Poly(phenylene) Proton Exchange Membranes

Multiple computational and experimental techniques are used here to understand the nanoscale morphology and water/proton transport properties in a series of sulfonated Diels–Alder poly(phenylene) (SDAPP) membranes over a wide range of temperature, hydration, and sulfonation conditions. New synthetic methods allow us to sulfonate the SDAPP membranes to much higher ion exchange capacity levels than has been previously possible. Nanoscale phase separation between the hydrophobic polymer backbone and the hydrophilic water/sulfonic acid groups was observed for all membranes studied. We find good agreement between structure factors calculated from atomistic molecular dynamics (MD) simulations and those measured by X-ray scattering. With increasing hydration, the scattering ionomer peak in SDAPP is found to decrease in intensity. This intensity decrease is shown to be due to a reduction of scattering contrast between the water and polymer and is not indicative of any loss of nanoscale phase separation. Both MD simulations and density functional theory (DFT) calculations show that as hydration levels are increased, the nanostructure morphology in SDAPP evolves from isolated ionic domains to fully percolated water networks containing progressively weaker hydrogen bond strengths. The conductivity of the membranes is measured by electrical impedance spectroscopy and the equivalent proton conductivity calculated from pulsed-field-gradient (PFG) NMR diffusometry measurements of the hydration waters. Comparison of the measured and calculated conductivity reveals that in SDAPP the proton conduction mechanism evolves from being dominated by vehicular transport at low hydration and sulfonation levels to including a significant contribution from the Grötthuss mechanism (also known as structural diffusion) at higher hydration and sulfonation levels. The observed increase in conductivity reflects the impact that changing hydration and sulfonation have on the morphology and hydrogen bond network and ultimately on the membrane performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Importance of elastic finite-size effects: Neutral defects in ionic compounds

Small system sizes are a well known source of error in DFT calculations, yet computational constraints frequently dictate the use of small supercells, often as small as 96 atoms in oxides and compound semiconductors. In ionic compounds, electrostatic finite size effects have been well characterised, but self-interaction of charge neutral defects is often discounted or assumed to follow an asymptotic behaviour and thus easily corrected with linear elastic theory. Here we show that elastic effect are also important in the description of defects in ionic compounds and can lead to qualitatively incorrect conclusions if inadequatly small supercells are used; moreover, the spurious self-interaction does not follow the behaviour predicted by linear elastic theory. Considering the exemplar cases of metal oxides with fluorite structure, we show that numerous previous studies, employing 96-atom supercells, misidentify the ground state structure of (charge neutral) Schottky defects. We show that the error is eliminated by employing larger cells (324, 768 and 1500 atoms), and careful analysis determines that elastic effects, not electrostatic, are responsible. The spurious self-interaction was also observed in non-oxide ionic compounds and irrespective of the computational method used, thereby resolving long standing discrepancies between DFT and force-field methods, previously attributed to the level of theory. The surprising magnitude of the elastic effects are a cautionary tale for defect calculations in ionic materials, particularly when employing computationally expensive methods (e.g. hybrid functionals) or when modelling large defect clusters. We propose two computationally practicable methods to test the magnitude of the elastic self-interaction in any ionic system. In commonly studies oxides, where electrostatic effects would be expected to be dominant, it is the elastic effects that dictate the need for larger supercells | greater than 96 atoms.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

First principles calculations in support of Pu aging: calculating the effects of lattice imperfections on thermodynamics

Plutonium (Pu) has, in theory, well defined crystal structures: its atoms are arranged in regular spatial patterns. But Pu is radioactive, and as its nuclei decay those regular spatial patterns are interrupted. The interruptions are lattice imperfections, which are known to affect how materials respond to their environment. To adequately model Pu, we need to know which lattice imperfections are present, how they interact with each other, and how they affect the material’s response to its environment. The work presented here aims to use density functional theory (DFT) calculations to begin to answer the latter, in particular, how individual lattice imperfections affect measurable effects including thermal expansion (the change in volume in response to a change in temperature), heat capacity (the amount of thermal energy needed to change a material’s temperature), and elastic moduli (a material’s resistance to applied stresses). Pu poses many computational challenges. The $\textit{f}$ electrons require special attention. Of all the elements, Pu has the largest number of electrons that must be included in the calculations. Thermal effects require calculating the phonons (the lattice vibrations), which for systems containing lattice imperfections demand large, complex computational cells - but computational resources limit the size and complexity. With careful restructuring of how the calculations are performed, all these challenges have been met to enable calculations that provide insight into how lattice imperfections affect Pu’s response to its environment. Reported here are the computational challenges and the advances developed to meet them, along with the first results showing the strong effect that one prototype lattice imperfection (an interstitial Pu atom in a delta-phase Pu lattice) has on Pu’s response to its environment.

36 MATERIALS SCIENCE↗

Spin Signature of the C 60 Fullerene Anion: A Combined X- and D-Band EPR and DFT Study

Fullerenes attract much attention in various scientific fields, but their electronic properties are still not completely understood. Here we report on a combined EPR and DFT study of the fullerene anion C 60 – in solid glassy environment. DFT calculations were used to characterize its electronic structure through spin density distribution and magnetic resonance parameters. The electron spin density is not uniformly distributed throughout the C 60 – cage but shows a pattern similar to PC 61 BM – . EPR spectroscopy reveals a rhombic g-tensor sensitive to the environment in the frozen glassy solutions, which can be rationalized by deformation of the fullerenes along lowfrequency vibrational modes upon cooling. DFT modeling confirms that these deformations lead to variation in the C 60 – g values. Furthermore, the decrease in g-tensor anisotropy with sample annealing is related to the lessening of g-tensor strain upon temperature relaxation of the most distorted sites in the glassy state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High reactivity of mesoporous CeO 2 to dissociate chemical warfare agent sarin

We perform spectroscopic surface characterization on mesoporous ceria (CeO 2 ) under exposure of sarin gas, also known as GB. We find mesoporous CeO 2 displays extraordinary reactivity towards GB dissociation under high vacuum. The spectroscopic results combined with density functional theory (DFT) calculations suggest that the cleavage of the P–F and P–OR bonds are both possible pathways for GB to dissociate on the CeO 2 surfaces. The modeling reveals differences in dissociation pathways of GB on the crystallographic (111)- and (110)-surfaces of ceria. Here we also report some dissociation of GB on mesoporous CeO 2 under ambient conditions. The finding that mesoporous CeO 2 can degrade GB under mild conditions is encouraging and implies the strong potential of ceria-based materials to be applied in combating chemical warfare agents related to the GB series.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atom Probe Tomography Analysis of Ag Doping in 2D Layered Material (PbSe) 5 (Bi 2 Se 3 ) 3

Impurity doping in two-dimensional (2D) materials can provide a route to tuning electronic properties, so it is important to be able to determine the distribution of dopant atoms within and between layers. Here we report the totnographic mapping of dopants in layered 2D materials with atomic sensitivity and subnanometer spatial resolution using atom, probe tomography (APT). Also, APT analysis shows that Ag dopes both Bi 2 Se 3 and PbSe layers in (PbSe) 5 (Bi 2 Se 3 ) 3 , and correlations :in the position of Ag atoms suggest a pairing across neighboring Bi 2 Se 3 and PbSe layers. Finally, density functional theory (DFT) calculations confirm the favorability of substitutional-doping for both Pb and Bi and provide insights into the,observed spatial correlations in dopant locations.

2D materials↗

Electronic, structural and magnetic studies of niobium borides of group 8 transition metals, Nb{sub 2}MB{sub 2} (M=Fe, Ru, Os) from first principles calculations

The Nb{sub 2}FeB{sub 2} phase (U{sub 3}Si{sub 2}-type, space group P4/mbm, no. 127) is known for almost 50 years, but until now its magnetic properties have not been investigated. While the synthesis of Nb{sub 2}OsB{sub 2} (space group P4/mnc, no. 128, a twofold superstructure of U{sub 3}Si{sub 2}-type) with distorted Nb-layers and Os{sub 2}-dumbbells was recently achieved, “Nb{sub 2}RuB{sub 2}” is still not synthesized and its crystal structure is yet to be revealed. Our first principles density functional theory (DFT) calculations have confirmed not only the experimental structures of Nb{sub 2}FeB{sub 2} and Nb{sub 2}OsB{sub 2}, but also predict “Nb{sub 2}RuB{sub 2}” to crystalize with the Nb{sub 2}OsB{sub 2} structure type. According to chemical bonding analysis, the homoatomic B–B interactions are optimized and very strong, but relatively strong heteroatomic M–B, B–Nb and M–Nb bonds (M=Fe, Ru, Os) are also found. These interactions, which together build a three-dimensional network, are mainly responsible for the structural stability of these ternary borides. The density-of-states at the Fermi level predicts metallic behavior, as expected, from metal-rich borides. Analysis of possible magnetic structures concluded preferred antiferromagnetic ordering for Nb{sub 2}FeB{sub 2}, originating from ferromagnetic interactions within iron chains and antiferromagnetic exchange interactions between them. -- Graphical abstract: Nb{sub 2}FeB{sub 2} (U{sub 3}Si{sub 2} structure type, space group P4/mbm, no. 127) is predicted to order antiferromagnetically, due to the presence of iron chains which show ferromagnetic interactions in the chains and antiferromagnetic interactions between them. “Nb{sub 2}RuB{sub 2}” is predicted to crystallize with the recently discovered Nb{sub 2}OsB{sub 2} twofold superstructure (space group P4/mnc, no. 128) of U{sub 3}Si{sub 2} structure type. The building of ruthenium dumbbells instead of chains along [001] is found to be responsible for the stabilization of this superstructure. Highlights: • Nb{sub 2}FeB{sub 2} is predicted to order antiferromagnetically. • Ferromagnetic interactions found in iron chains and antiferromagnetic ones between them. • Unknown “Nb{sub 2}RuB{sub 2}” predicted to crystallize with a twofold U{sub 3}Si{sub 2} superstructure. • Puckering of Nb-layer and Ru-dumbbell formation responsible for superstructure occurrence.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Quantum Mechanical Study of Quartz (101) Interfacial Boundaries

Density functional theory (DFT) calculations were performed on periodic and molecular models to explore the energetics of bonding modes at quartz–quartz interfaces, common grain contacts in sandstones. Four interface types were modeled: H 2 O-mediated, silanol H-bond mediated, O–O peroxide bonds, and siloxane bonds. Each type of interface may exist in the subsurface of the Earth, depending upon water activity, temperature, and pressure, as interfacial structures are stable within DFT molecular dynamics simulations up to 473 K. The results predict interfacial energies SiOSi > SiOOSi > SiOH-6H 2 O > 2SiOH, as expected based on similar bond strengths in other systems. The proposed SiOOSi interface is not generally accounted for in chemomechanical models but could be important in sandstones at depth in the Earth. The implications for intergrain fracturing and H2 generation are discussed.

Hydrogen↗

Adaptive DFT-based Interferometer Fringe Tracking

An automatic interferometer fringe tracking system has been developed, implemented, and tested at the Infrared Optical Telescope Array (IOTA) observatory at Mt. Hopkins, Arizona. The system can minimize the optical path differences (OPDs) for all three baselines of the Michelson stellar interferometer at IOTA. Based on sliding window discrete Fourier transform (DFT) calculations that were optimized for computational efficiency and robustness to atmospheric disturbances, the algorithm has also been tested extensively on off-line data. Implemented in ANSI C on the 266 MHz PowerPC processor running the VxWorks real-time operating system, the algorithm runs in approximately 2.0 milliseconds per scan (including all three interferograms), using the science camera and piezo scanners to measure and correct the OPDs. The adaptive DFT-based tracking algorithm should be applicable to other systems where there is a need to detect or track a signal with an approximately constant-frequency carrier pulse.

Wilson, Edward↗

Photoionization of Benzene and Small Polycyclic Aromatic Hydrocarbons in Ultraviolet-Processed Astrophysical Ices: A Computational Study

We employed density functional theory (DFT) calculations to model the photoionization behavior of benzene and small polycyclic aromatic hydrocarbons when they are embedded in a matrix of water ice in order to investigate issues raised by recent experimental work by Gudipati and Allamandola. The ionization energies of benzene, naphthalene, anthracene, and pyrene were found to be lowered by 1.5-2.1 eV in water ice. Low-lying vertical electronic excitation energies were computed with time-dependent DFT for both neutral and ionized species and are found in both cases to be remarkably unaffected by the ice matrix. Chemical behavior in ultraviolet-photoprocessed ices is also discussed, with a focus on electron recombination and pathways leading to phenol and analogous products.

Woon, D. E.↗

Density functional theory calculations of self- and Xe diffusion in U 3 Si 2

Uranium silicide, U 3 Si 2 , has been proposed as an advanced nuclear fuel to be used in light water reactors (LWRs). Development of this alternative to the predominant current fuel, UO 2 , is motivated by enhanced accident tolerance as a result of higher thermal conductivity as well as improved fuel cycle economics through increased uranium density. In order to accurately model the fuel performance of U 3 Si 2 , the diffusion rate of point defects, which is related to self-diffusion, and of fission gas atoms must be determined. DFT calculations are used to predict the U and Si point defect concentrations, the corresponding self-diffusivities, the preferred Xe trap site and the Xe diffusivity. Effects of irradiation are not considered. A low defect formation energy and a high entropy for Si interstitials give rise to Si-rich non-stoichiometry at elevated temperatures. Both U and Si self-diffusion and Xe diffusion are anisotropic as a consequence of the tetragonal crystal structure of U 3 Si 2 . Si diffusion occurs by interstitial mechanisms in both the a-b plane and along the c axis, while the U c axis diffusion rate is controlled by a vacancy mechanism. Interstitial diffusion of U is very fast in the a-b plane of the U 3 Si 2 crystal structure. Xe atoms prefer to occupy U vacancy trap sites. The highest Xe diffusion rate occurs by a vacancy mechanism in both the a-b plane and along the c axis. The diffusion rate is similar in the a-b plane and along the c axis. U and Si self-diffusion and Xe diffusion are all faster in U3Si2 than intrinsic U and Xe diffusion in conventional UO 2 nuclear fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Boron Substituted Na 3 V 2 (P 1 -x B x O 4 ) 3 Cathode Materials with Enhanced Performance for Sodium-Ion Batteries

The development of excellent performance of Na-ion batteries remains great challenge owing to the poor stability and sluggish kinetics of cathode materials. Herein, B substituted Na 3 V 2 P 3–x B x O 12 (0 ≤ x ≤ 1) as stable cathode materials for Na-ion battery is presented. A combined experimental and theoretical investigations on Na 3 V 2 P 3–x B x O 12 (0 ≤ x ≤ 1) are undertaken to reveal the evolution of crystal and electronic structures and Na storage properties associated with various concentration of B. X-ray diffraction results indicate that the crystal structure of Na 3 V 2 P 3–x B x O 12 (0 ≤ x ≤ 1/3) consisted of rhombohedral Na 3 V 2 (PO 4 ) 3 with tiny shrinkage of crystal lattice. X-ray absorption spectra and the calculated crystal structures all suggest that the detailed local structural distortion of substituted materials originates from the slight reduction of V–O distances. Na 3 V 2 P 3-1/6 B 1/6 O 12 significantly enhances the structural stability and electrochemical performance, giving remarkable enhanced capacity of 100 and 70 mAh g -1 when the C-rate increases to 5 C and 10 C. Spin-polarized density functional theory (DFT) calculation reveals that, as compared with the pristine Na 3 V 2 (PO 4 ) 3 , the superior electrochemical performance of the substituted materials can be attributed to the emergence of new boundary states near the band gap, lower Na + diffusion energy barriers, and higher structure stability.

36 MATERIALS SCIENCE↗

A comprehensive diffusion mobility database comprising 23 elements for magnesium alloys

We report that reliable experimental diffusion coefficients of 10 key alloying elements in Mg obtained by the present authors together with experimental data in the literature enabled us to perform a systematic test of the reliability of diffusion coefficients obtained from DFT calculations. The computed activation energy values were found to be quite accurate (mostly within 0.2 eV) but the computed pre-factors were less reliable. Such insights allowed us to develop a practical and yet robust strategy to perform diffusion mobility assessments by adopting the computed activation energy while fitting only the pre-factor when available experimental data are limited to a narrow temperature range. The overall good agreement between the DFT data and experimental data also gave us the confidence to employ the computed data for those that were still missing or inaccessible from experimental measurements. A systematic assessment of both the measured and computed diffusion data in hcp Mg was performed using the above holistic approach to yield the most comprehensive open Mg mobility database to date, comprising 23 elements (Mg, Ag, Al, Be, Ca, Cd, Ce, Cu, Fe, Ga, Gd, In, La, Li, Mn, Nd, Ni, Pu, Sb, Sn, U, Y, Zn). This more reliable mobility database will contribute to future development of advanced Mg alloys. The holistic approach developed in this study will be very beneficial to the future establishment of reliable mobility databases for other alloy systems as well.

36 MATERIALS SCIENCE↗