Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DENSITY FUNCTIONAL METHOD”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Valence-to-core X-ray emission spectroscopy of vanadium oxide and lithiated vanadyl phosphate materials

In this work, we report valence-to-core (VTC) X-ray emission spectroscopy (XES) measurements and theoretical calculations of the electrochemical sequence ε-VOPO 4 , ε-LiVOPO 4 , ε-Li 2 VOPO 4 and the reference oxides V 2 O 3 , VO 2 , and V 2 O 5 . In our analysis of these results, we establish a framework for interrogating chemical bonding that is generally applicable to a wide range of systems, including complex, extended inorganic compounds. While this latter regime has garnered less focused application than, e.g. , metalloenzymes in many excellent catalysis studies, we show that the technique provides high utility in materials-focused energy storage research. Here, sensitivities to the local atomic structure and hybridization schemes are discussed in detail. Similarly, the effect of lithiation on oxidation, delocalization, and shifts in ligand valence energy levels are all readily apparent in the analyzed results. Finally, the TDDFT projections clearly reveal the directional dependencies of the valence band at each of the vanadium sites. Our results demonstrate laboratory-based X-ray spectroscopy instrumentation is a viable route for attaining well-resolved VTC-XES features for inorganic compounds of 3d transition metals, even for samples of limited quantity or suffering from sensitivity to the atmosphere. The experimental results are in good agreement with results produced by real-space Green's function and time-dependent density functional theory (TDDFT) methods, respectively. Hence, we propose that VTC-XES, when equipped with appropriate theoretical support, can be a valuable complement to X-ray absorption pre-edge features for more detailed characterization of a compound's electronic structure. We expect similar analyses will find application in a broad range of materials chemistry research and provide both fundamental and applied insights.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

A substituent-dependent deep-blue pyrene-based chemosensor for trace nitroaniline sensing

Nitroaniline (NA) is a kind of high toxicity and pathogenic amine derivative, and the trace monitoring of nitroaniline has become an urgent problem in terms of both the environment and security. Herein, two K-region multiple substituted pyrene-based luminogens with reasonable photoluminescence efficiency and blue emission (408–430 nm) have been synthesized and characterized. Their photophysical properties were carefully investigated in solution and in the solid state by experimental and density functional theory (DFT) methods. It is worth mentioning that both of these π-electron rich and extended π-conjugated materials were used for the detection of electron deficient nitroaromatics (NACs) based on a photo-induced electron transfer (PET) mechanism, and a sensitive response towards NA, especially for p-NA with a limit of detection (LOD) as low as 8.58 × 10 -9 M, was observed, with a fluorescence quenching efficiency up to 99.5%. In conclusion, these research results provided a low-cost and simple strategy for trace nitroaniline sensing.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Advanced thermal/environmental barrier coatings of high-entropy rare earth disilicates tuned by strong anharmonicity of Eu 2 Si 2 O 7

Advancing thermal/environmental barrier coating (TEBC) materials with integrated thermal-mechanical functions is paramount for safeguarding SiC-based ceramic matrix composites (CMCs) in high-efficiency gas turbines. Herein, we employ a synergistic approach, combining density functional theory (DFT) methods and combinatorial chemistry techniques, to design high-performance and low-cost RE 2 Si 2 O 7 (RE = rare earth elements) TEBC materials tailored for enhanced compatibility with SiC-based CMCs. Expanding on phase stability of alloying pure RE 2 Si 2 O 7 , the investigation extends to the mechanical and thermal properties of solid solution systems, including Er 1/2 Y 3/4 Yb 3/4 Si 2 O 7 , Gd 1/4 Er 1/4 Y 3/4 Yb 3/4 Si 2 O 7 , and Eu 1/4 Er 1/4 Y 3/4 Yb 3/4 Si 2 O 7 . The solid solution systems exhibit a major reduction in lattice thermal conductivity relative to their pure counterparts, achieving ultralow values of 0.25 to 0.39 W m −1 K −1 at 1500 K. Furthermore, the coefficients of thermal expansion (CTE) of these solid solutions are precisely tuned within the desired range for SiC (4.4 to 5.5 × 10 −6 K −1 ), while maintaining good mechanical properties. Here, in particular, the addition of Eu 2 Si 2 O 7 demonstrates to be an important variable to the tuning of CTE and lattice thermal conductivity by leveraging its strong anharmonicity, presenting a pioneering avenue for fine-tuning material properties. In summary, this research not only identifies promising TEBC materials with superior thermal properties, but also introduces a valuable computational material design methodology for the rapid discovery of complex materials for harsh environments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Hyperfine and quadrupole interactions for Dy isotopes in DyPc 2 molecules

Nuclear spin levels play a vital role in understanding magnetization dynamics and implementation and control of quantum bits in lanthanide-based single-molecule magnets. We investigate the hyperfine and nuclear quadrupole interactions for 161 Dy and 163 Dy nuclei in anionic DyPc 2 (Pc=phthalocyanine) single-molecule magnets, using multiconfigurational ab initio methods (beyond density-functional theory) including spin–orbit interaction. The two isotopes of Dy are chosen because the others have zero nuclear spin. Both isotopes have the nuclear spin I =5/2, although the magnitude and sign of the nuclear magnetic moment differ from each other. The large energy gap between the electronic ground and first-excited Kramers doublets, allows us to map the microscopic hyperfine and quadrupole interaction Hamiltonian onto an effective Hamiltonian with an electronic pseudo-spin that corresponds to the ground Kramers doublet. Our ab initio calculations illustrate that the coupling between the nuclear spin and electronic orbital angular momentum contributes the most to the hyperfine interaction and that both the hyperfine and nuclear quadrupole interactions for 161 Dy and 163 Dy nuclei are much smaller than those for the 159 Tb nucleus in TbPc 2 single-molecule magnets. The calculated separations of the electronic-nuclear levels are comparable to experimental data reported for 163 DyPc 2 . We demonstrate that hyperfine interaction for the Dy Kramers ion leads to tunnel splitting (or quantum tunneling of magnetization) at zero field. This effect does not occur for TbPc 2 single-molecule magnets. The magnetic field values of the avoided level crossings for 161 DyPc 2 and 163 DyPc 2 are found to be noticeably different, which can be observed from the experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Density functional study of the variants of inter-Coulombic decay resonances in the photoionization of Cl@C 60

Inter-Coulombic decay (ICD) resonances in the photoionization of Cl@C 60 endofullerene molecule are calculated using a perturbative density functional theory (DFT) method. This is the first ICD study of an open shell atom in a fullerene cage. Here, three classes of resonances are probed: (i) Cl inner vacancies decaying through C 60 outer continua, (ii) C 60 inner vacancies decaying through Cl outer continua, and (iii) inner vacancies of either system decaying through the continua of Cl-C 60 hybrid levels, the hybrid Auger-ICD resonances. Comparisons with Ar@C 60 results reveal that the properties of hybrid Auger-ICD resonances are affected by the extent of level hybridization.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Orbital-selective Mott phase and non-Fermi liquid in FePS 3

The layered metal phosphorous trisulfide FePS 3 is reported to be a Mott insulator at ambient conditions and to undergo structural and insulator-metal phase transitions under pressure. However, the character of the resulting metallic states has not been understood clearly so far. Here, we theoretically study the phase transitions of FePS 3 using first-principles methods based on density functional theory and embedded dynamical mean field theory. We find that the Mott transition in FePS 3 can be orbital selective, with t 2g states undergoing a correlation-induced insulator-to-metal transition while e g states remain gapped. We show that this orbital-selective Mott phase, which occurs only when nonhydrostatic pressure is used, is a bad metal (or non-Fermi liquid) with large fluctuating moments due to Hund's coupling. Further application of pressure increases the crystal-field splitting and converts the system to a conventional Fermi liquid with low-spin configurations dominant. Finally, our results show that FePS 3 is an example of a system that realizes an orbital-selective Mott phase, allowing tuning between correlated and uncorrelated metallic properties in an accessible pressure range (≤ 18 GPa).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

In-Silico Analysis of High Refractive Index Materials Through Principles of Materials Design

The intent of the paper is to use specific principles of Materials Design that were developed and applied in the electronics industry for enabling understanding and design of improved high refractive index materials. Further, by combining first-principle based ab-initio, semiempirical interatomic potential methods, and machine learning approaches in conjunction with experimental data, we identified specific determinants of high refractive index materials, which can be critically applied for informing materials design and accelerating discovery. Specifically, it was demonstrated that chalcogenides and perovskites as bulk materials can exhibit higher refractive indices with appropriate engineering of specific aspects of the materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Radioactive decay of $\mathrm {{}^{90}Sr}$ in cement: a non-equilibrium first-principles investigation

Cement is an inexpensive and relatively easily manageable material that is used as a last barrier for nuclear waste disposal. Under these conditions, the cement is in contact with low radiation doses, but there is a distinct possibility of being contaminated with radioactive products. Of particular concern is the medium lived half-life product 90 Sr (28.8 years) due to its ability to replace Ca. 90 Sr undergoes β-decay to 90 Y which, in turn, β-decays to stable 90 Zr. In this work, we discuss systematically the chain of non-equilibrium processes that result as a consequence of β-decay events in cement. We first use density functional-based methods to study the consequences of the sudden increase of the nuclear charge from Z to Z+1, a possible induced ionization and the perturbation of the surrounding electronic charge. Secondly, we use molecular dynamics simulations to study the recoil of the daughter nucleus. Finally, we discuss the damage caused by the ionization cascade produced during the propagation of the β-electron and the resulting chemical and structural perturbation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Machine learning assisted prediction of tungsten heavy alloy plasma facing component performance for fusion energy applications

Tungsten and tungsten heavy alloys (WHAs), known for their remarkably high hardness, durability, and corrosion resistance, play a critical role in the thriving development of nuclear fusion reactors in recent years. However, the exploration in tungsten alloys for the nuclear-related applications has been limited by the difficulty of manufacturing and the complexity of experiments to reproduce the environment of nuclear reaction. Therefore, this project aims to utilize nanoscale simulation methods such as density functional theory (DFT) and molecular dynamics (MD) with the help of machine learning techniques to not only understand the mechanisms of tungsten alloys but also allow us to computationally predict their mechanical behaviors under extreme environments. One critical problem of the application of WHAs in nuclear reactors is the surface melting. In the current design of the SPARC reactor, the WHA, W97Ni2.1Fe0.9 or W97NiFe, is chosen to be the first wall components to confine the plasma where the particles are fiercely moving and colliding into each other to create nuclear fusion reaction. This process will generate extremely high heat flux onto these WHA tiles, leaving high surface temperature that could possibly melt the surface of the WHA tiles, As illustrated in Fig. 1(a). a laser experiment previously done illustrates that a rough surface damage would be made after the surface melting where the matrix area mainly composed of nickel and iron as shown in Fig. 1(b), will first melt and then leave vacancies between these tungsten grains. Unfortunately, these kinds of roughness on the first-wall components could be deadly to the plasma inside a Tokmak reactor because the heat that is supposed to dissipate at a designed ratio through the tiles may in turn be excessively absorbed and accumulated on any uneven area of the surface, which will eventually make the whole nuclear reaction fail. In this project, we will introduce a machine learning potential, Allegro, based on DFT calculation and then build a MD model for W-Ni-Fe alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First-Principles Studies of Tritium Species Diffusivity Across the Interface of Nickel-Plated Zircaloy-4

In this study, we will calculate the 3H2 and 3H2O diffusion mechanism in the Ni coated layer and the dominant 3H species dissolution into the Zircaloy-4 getter. We will optimize the Ni/Zr interface and introduce Sn and O impurities. To achieve the objectives, we will perform the investigation by using spin-polarized density functional theory (DFT) methods to capture the effects of Ni magnetic properties on 3H diffusion. We will introduce NiO/Ni(OH) on the Ni side closer to the interface. Then, the 3H diffusion pathways and barriers from Ni side to the zircalory-4 will be explored. We will investigate how NiO/Ni(OH) affects the 3H diffusion across the interface.

Paudel, Hari P.↗
DOE OSTI2024

Quantum Monte Carlo Calculations of Chemical Binding and Reactions

The auxiliary field quantum Monte Carlo method developed by the PIs has been shown to provide the most accurate description of strongly correlated electronic systems, from molecules to solids. Unlike other explicitly many‐body approaches, the quantum Monte Carlo method scales as a low order polynomial of systems size, similar to mean‐field methods such as density functional theory. However, the auxiliary field quantum Monte Carlo algorithm is significantly more expensive than traditional density functional calculations. This creates a bottleneck for applications to extended systems, such as large molecules and solids. One principal objective of this proposal was to develop new auxiliary field quantum Monte Carlo computational strategies to achieve improved scaling with system size, using downfolding and localization schemes, without sacrificing the predictive power of the calculations. A second goal is to extend the reach of auxiliary field quantum Monte Carlo to calculate excited states. This final report summarizes what has been achieved during the course the project toward these goals.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Augmenting Molecular Graphs with Geometries via Machine Learning Interatomic Potentials

Accurate molecular property predictions require 3D geometries, which are typically obtained using expensive methods such as density functional theory (DFT). Here, we attempt to obtain molecular geometries by relying solely on machine learning interatomic potential (MLIP) models. To this end, we first curate a large-scale molecular relaxation dataset comprising 3.5 million molecules and 300 million snapshots. Then MLIP pre-trained models are trained with supervised learning to predict energy and forces given 3D molecular structures. Once trained, we show that the pre-trained models can be used in different ways to obtain geometries either explicitly or implicitly. First, it can be used to obtain approximate low-energy 3D geometries via geometry optimization. While these geometries do not consistently reach DFT-level chemical accuracy or convergence, they can still improve downstream performance compared to non-relaxed structures. To mitigate potential biases and enhance downstream predictions, we introduce geometry fine-tuning based on the relaxed 3D geometries. Second, the pre-trained models can be directly fine-tuned for property prediction when ground truth 3D geometries are available. Our results demonstrate that MLIP pre-trained models trained on relaxation data can learn transferable molecular representations to improve downstream molecular property prediction and can provide practically valuable but approximate molecular geometries that benefit property predictions. Our code is publicly available at: https://github.com/divelab/AIRS/.

Fu, Cong [Texas A & M Univ., College Station, TX (↗
DOE OSTI2021

Radiation-Induced Oxidation Reactions of 2-Selenouracil in Aqueous Solutions: Comparison with Sulfur Analog of Uracil

One-electron oxidation of 2-selenouracil (2-SeU) by hydroxyl (●OH) and azide (●N3) radicals leads to various primary reactive intermediates. Their optical absorption spectra and kinetic characteristics were studied by pulse radiolysis with UV-vis spectrophotometric and conductivity detection and by the density functional theory (DFT) method. The transient absorption spectra recorded in the reactions of ●OH with 2-SeU are dominated by an absorption band with an λmax = 440 nm, the intensity of which depends on the concentration of 2-SeU and pH. Based on the combination of conductometric and DFT studies, the transient absorption band observed both at low and high concentrations of 2-SeU was assigned to the dimeric 2c-3e Se-Se-bonded radical in neutral form (2●). The dimeric radical (2●) is formed in the reaction of a selenyl-type radical (6●) with 2-SeU, and both radicals are in equilibrium with Keq = 1.3 × 104 M−1 at pH 4 (below the pKa of 2-SeU). Similar equilibrium with Keq = 4.4 × 103 M−1 was determined for pH 10 (above the pKa of 2-SeU), which admittedly involves the same radical (6●) but with a dimeric 2c-3e Se-Se bonded radical in anionic form (2●−). In turn, at the lowest concentration of 2-SeU (0.05 mM) and pH 10, the transient absorption spectrum is dominated by an absorption band with an λmax = 390 nm, which was assigned to the ●OH adduct to the double bond at C5 carbon atom (3●) based on DFT calculations. Similar spectral and kinetic features were also observed during the ●N3-induced oxidation of 2-SeU. In principle, our results mostly revealed similarities in one-electron oxidation pathways of 2-SeU and 2-thiouracil (2-TU). The major difference concerns the stability of dimeric radicals with a 2c-3e chalcogen-chalcogen bond in favor of 2-SeU.

2-selenouracil↗
DOE OSTI2023

Self-consistent implementation of locally scaled self-interaction-correction method

Recently proposed local self-interaction correction (LSIC) method is a one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange–correlation and Coulomb energy densities. The LSIC method is exact for the one-electron densities, also recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to the well-known Perdew–Zunger SIC (PZSIC) method as a special case. This article presents the self-consistent implementation of the LSIC method using the ratio of Weizsäcker and Kohn–Sham kinetic energy densities as an iso-orbital indicator. The atomic forces as well as the forces on the Fermi-Löwdin orbitals are also implemented for the LSIC energy functional. Results show that LSIC with the simplest local spin density functional predicts atomization energies of the AE6 dataset better than some of the most widely used generalized-gradient-approximation (GGA) functional and barrier heights of the BH6 database better than some of the most widely used hybrid functionals. The LSIC method [a mean absolute error (MAE) of 0.008 Å] predicts bond lengths of a small set of molecules better than the PZSIC-LSDA (MAE 0.042 Å) and LSDA (0.011 Å). Furthermore, this work shows that accurate results can be obtained from the simplest density functional by removing the self-interaction-errors using an appropriately designed SIC method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

How the Size and Density of Charge-Transfer Excitons Depend on Heterojunction’s Architecture

We have characterized the size, intensity, density, and distribution of charge-transfer (CT) excitons as a function of the acceptor–donor architecture of prototypical organic interfaces. This characterization was done by computational analysis of 17 models of varying numbers, positions, and orientations of the donor and acceptor molecules. The models’ building blocks were phenyl-C61-butyric acid methyl ester (PCBM) fullerene acceptors and dual-band donor polymers composed of thiophene, benzothiadiazole, and benzotriazole subunits. The electronic structure of the donor–acceptor complexes was computed with the time-dependent long-range-corrected density-functional tight-binding method and analyzed with the fragment-based one-electron transition density matrix. In all models, the complexes with edge-on orientation have denser spectra of low-energy CT states lying below the absorption bands compared to the complexes with face-on orientation. This CT-state distribution in edge-on complexes provides a gate to efficiently populate cold CT excitons. Moreover, the cold CT excitons have a higher degree of charge separation in the edge-on than in the face-on complexes. The CT amount and the CT exciton size generally increase with the energy of the CT states, although the electron remains localized on a single molecule in cold CT states. Delocalization over two PCBM molecules was observed for high-energy CT states. The exciton size also depends on the orientation. Larger excitons are produced by the delocalization of the electrons perpendicularly to the interface. When the delocalization is parallel, the smaller electron–hole distances yield moderately sized CT excitons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Time-dependent density functional theory with the orthogonal projector augmented wave method

The projector augmented wave (PAW) method of Blöchl linearly maps smooth pseudo wavefunctions to the highly oscillatory all-electron DFT orbitals. Compared to norm-conserving pseudopotentials (NCPP), PAW has the advantage of lower kinetic energy cutoffs and larger grid spacing at the cost of having to solve for non-orthogonal wavefunctions. We earlier developed orthogonal PAW (OPAW) to allow the use of PAW when orthogonal wavefunctions are required. In OPAW, the pseudo wavefunctions are transformed through the efficient application of powers of the PAW overlap operator with essentially no extra cost compared to NCPP methods. Previously, we applied OPAW to DFT. Here, we take the first step to make OPAW viable for post-DFT methods by implementing it in real-time time-dependent (TD) DFT. Using fourth-order Runge–Kutta for the time-propagation, we compare calculations of absorption spectra for various organic and biological molecules and show that very large grid spacings are sufficient, 0.6–0.7 bohr in OPAW-TDDFT rather than the 0.4–0.5 bohr used in traditional NCPP-TDDFT calculations. This reduces the memory and propagation costs by around a factor of 3. Our method would be directly applicable to any post-DFT methods that require time-dependent propagations such as the GW approximation and the Bethe–Salpeter equation.

Chemistry↗
DOE OSTI2026

A physics informed model for the prediction of surface energies and the shape of metal particles

Surface energies of metal-based systems are important for determining the Wulff-constructed shapes of metal nanoparticles and understanding the stability. We have developed a physics informed method to predict the total energy of metal-based systems across a wide range of configurations. Our method has been tested against density functional theory (DFT) calculations for late transition metals. This method enables on-the-fly surface energy predictions based on the fundamental of materials property and allows for the Wulff construction of metal particles for a random number of elemental atoms and without the need for DFT calculations. By making a division between atoms in the different layers of the model system, we can improve the accuracy of the model, suggesting a dissimilarity between the electronic structure due to an alternating compression and expansion of atomic layers. We find that our model accurately and effectively provides valuable insights into the distribution and stability of nanoparticles.

Mandal, Shyama Charan (ORCID:0000000245882874)↗