Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional”

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 1,639 records · Page 91

Frustration and Atomic Ordering in a Monolayer Semiconductor Alloy

Frustrated interactions can lead to short-range ordering arising from incompatible interactions of fundamental physical quantities with the underlying lattice. The simplest example is the triangular lattice of spins with antiferromagnetic interactions, where the nearest-neighbor spin-spin interactions cannot simultaneously be energy minimized. In this work, we show that engineering frustrated interactions is a possible route for controlling structural and electronic phenomena in semiconductor alloys. Using aberration-corrected scanning transmission electron microscopy in conjunction with density functional theory calculations, we demonstrate atomic ordering in a two-dimensional semiconductor alloy as a result of the competition between geometrical constraints and nearest-neighbor interactions. Statistical analyses uncover the presence of short-range ordering in the lattice. In addition, we show how the induced ordering can be used as another degree of freedom to considerably modify the band gap of monolayer semiconductor alloys.

2-dimensional systems↗

Using simultaneous x-ray diffraction and velocity interferometry to determine material strength in shock-compressed diamond

We determine the strength of laser shock-compressed polycrystalline diamond at stresses above the Hugoniot elastic limit using a technique combining x-ray diffraction from the Linac Coherent Light Source with velocity interferometry. X-ray diffraction is used to measure lattice strains, and velocity interferometry is used to infer shock and particle velocities. These measurements, combined with density-dependent elastic constants calculated using density functional theory, enable determination of material strength above the Hugoniot elastic limit. Furthermore, our results indicate that diamond retains approximately 20 GPa of strength at longitudinal stresses of 150–300 GPa under shock compression.

36 MATERIALS SCIENCE↗

Binding of noble gas atoms by superhalogens

Because of their closed shells, noble gas (Ng) atoms (Ng = Ne, Ar, Kr, and Xe) seldom take part in chemical reactions, yet finding such mechanisms not only is of scientific interest but also has practical significance. Following a recent work by Mayer et al. [Proc. Natl. Acad. Sci. U. S. A. 116, 8167–8172 (2019)] on the room temperature binding of Ar to a superelectrophilic boron site embedded in a negative ion complex, B 12 (CN) 11 - , we have systematically studied the effect of cluster size and terminal ligands on the interaction of Ng by focusing on B 12 X 11 (Ng) (X = H, CN, and BO) and B 12 X 10 (Ng) 2 (X = CN and BO) whose stabilities are governed by the Wade–Mingos rule and on C 5 BX 5 (Ng) (X = H, F, and CN) and C 4 B 2 (CN) 4 (Ng) 2 whose stabilities are governed by the Huckel’s aromaticity rule. Our conclusion, based on density functional theory, is that both the cluster size and the terminal ligands matter—the interaction between the cluster and the Ng atoms becomes stronger with increasing cluster size and the electron affinity of the terminal ligands. Our studies also led to a counter-intuitive finding—removing multiple terminal ligands can enable electrophilic centers to bind multiple Ng atoms simultaneously without compromising their binding strength.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Candidate structure for the H 2 -PRE phase of solid hydrogen

Experimental progress finally reached the metallic solid hydrogen phase, which was predicted by Wigner and Huntington over 80 years ago. However, the different structures in the phase diagram are still being debated due to the difficulty of diffraction experiments for high-pressured hydrogen. The determination of crystal structures under extreme condition is both of the basic condensed matter physics, and in planetary science: the behavior of giant gaseous planets (e.g., Jupiter and Saturn) strongly depends on the properties of inner high-pressured hydrogen. This work describes possible structures appearing under high pressures of 400~600 GPa. We applied a structural search using particle swarm optimization with density functional theory (DFT) to propose several candidate structures. For these structures, we performed fixed-node diffusion Monte Carlo simulations combined with DFT zero-point energy corrections to confirm their relative stability. We found P2 1 /c-8 as a promising candidate structure for the H 2 -PRE phase. P2 1 /c-8 is predicted the most stable at 400 and 500 GPa. P2 1 /c-8 reproduces qualitatively the IR spectrum peaks observed in the H 2 -PRE phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ab Initio Transport Calculations: From Normal to Superconducting Current

Applying the Bogoliubov-de Gennes equations with density-functional theory, it is possible to formulate first principles description of current-phase relationships in superconducting/normal (magnetic)/superconducting tri-layers. Such structures are the basis for the superconducting analog of Magnetoresistive random access memory devices (JMRAM). In a recent paper1 we presented results from the first attempt to formulate such a theory, applied to the Nb/Ni/Nb trilayers. In the present work we provide computational details, explaining how to construct key ingredient (scattering matrices SN ) in a framework of linear muffin-tin orbitals (LMTO).

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗

Multipole responses in fissioning nuclei and their uncertainties

Electromagnetic multipole responses are key inputs to model the structure, decay, and reactions of atomic nuclei. With the introduction of the finite amplitude method (FAM), large-scale calculations of the nuclear linear response in heavy deformed nuclei have become possible. This work provides a detailed study of multipole responses in actinide nuclei with Skyrme energy density functionals. We quantify both systematic and statistical uncertainties induced by the functional parametrization in FAM calculations. We also extend the FAM formalism to perform blocking calculations with the equal filling approximation for odd-mass and odd-odd nuclei, and analyze the impact of blocking configurations on the response. In conclusion, by examining the entire plutonium isotopic chain from the proton to the neutron drip line, we find a large variability of the response with the neutron number and study how it correlates with the deformation of the nuclear ground state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Electron Thermal Conductivity of Pu and Zr Substituted Gamma-Uranium

Uranium alloys are attractive recycled nuclear fuels because of their high thermal conductivity (k) and fissile density; however, the effects of alloying elements on k remain unclear. Here, the electron thermal conductivity (k_e) of U-Pu-Zr compositions are calculated using density functional theory. The electronic structure is evaluated to understand the effects of plutonium (Pu) and zirconium (Zr) substitution on the k_e of ?-U. Alloys of up to 37.5 at. % Pu and 37.5 at. % Zr are examined. Two methods are applied to calculate k_e; we find that the accuracy of each method depends on the electronic and mass similarities between the solute and solvent atoms. Specifically, when the solute atom is similar in electronic structure and mass, the method that applies the electron relaxation time of ?-U is best, while if the elements are dissimilar, a mixed method that mixes several parameters associated with k_e from each element in the alloy is best. The introduction of all alloying elements decreases k_e; however, in binary compounds, Pu and Zr have different effects. Pu generally flattens the electronic bands but compensates for this deleterious effect by increasing electron density near the Fermi level. Zr flattens the electronic bands more severely without adding electron density near the Fermi level. Therefore, Zr decreases the k_e more than Pu in binary compounds. In ternary compounds, the difference between Pu and Zr is minimal due to the phononic change from the large mass change of Zr substitution, even at 12.5 at. %. Thus, we predict that higher loadings of Pu, and potentially other actinides, can be added to U-Pu-Zr compositions for faster recycling of spent fuel with without sacrificing k. We also note that these k_e calculation methods can be applied to non-fuel alloys that require k_e predictions, such as cladding, heat exchanger, and structural materials.

36 MATERIALS SCIENCE↗

Adsorption and activation of CO 2 on Pt/CeOx/TiO 2 (110): Role of the Pt-CeO x interface

The adsorption and dissociation of CO 2 on TiO 2 (110), CeO x /TiO 2 (110) and Pt/CeO x /TiO 2 (110) surfaces has been examined using Ambient Pressure X-ray Photoelectron Spectroscopy (AP-XPS). The substrates under study exhibited different degrees of complexity which were tested for the binding of the adsorbate and the cleavage of C-O bonds. The surfaces were prepared by depositing CeO x (0.1 ML) onto TiO 2 (110) to form a mixed oxide support, onto which Pt nanoparticles (0.2 ML) were deposited. This configuration yields a complex set of interfaces between metal and oxides and we have systematically titrated the active role of each component (Pt 4f, Ce 3d and Ti 2p regions) and the arising surface intermediates (C 1s and O 1s regions). CO 2 barely bonds to stoichiometric TiO 2 (110). It heals oxygen vacancies of this oxide surface (CO 2,gas → CO gas + O a ) and does not form stable carbonates. A stable carbonate was seen upon adsorption of CO 2 on CeO x /TiO 2 (110) and on this type of substrate the adsorbate also removed O vacancies leading to the oxidation of Ti 3+ and Ce 3+ sites. Pt nanoparticles dispersed on CeO x /TiO 2 (110) were highly effective for the binding and dissociation of CO 2 , with the formation of CO 3 , CO, C and CH x species on the Pt/CeO x /TiO 2 (110) system. The results of theoretical calculations based on density-functional theory (DFT) show that Pt/CeOx/TiO2(110) binds CO2 much stronger than surfaces of bulk platinum {(111), (100), (110)} or other late transition metals. On a Pt-CeO x interface, the molecule adsorbs with a bent configuration (~ 130° O-C-O bond angle) and with a substantial elongation (~ 0.1 Å) of the C-O bonds, facilitating its transformation into high value chemicals.

36 MATERIALS SCIENCE↗

Theoretical Insights on the Fragmentation of Cellulosic Oligomers to Form Hydroxyacetone and Hydroxyacetaldehyde

The presence of heavy unknown oligomeric sugar products in bio-oil is evidenced in experimental results reported in the literature. In this paper, we study the fragmentation reactions yielding acetol and glycolaldehyde from oligomeric sugars following previous work on dehydration reactions to propose structures of these oligomers. Acetol and glycolaldehyde are primary products of cellulose fast pyrolysis but the fragmentation reaction mechanism of these compounds from oligomers merits further study. The density functional theory (DFT) approach was employed to study this reaction. Results revealed that acetol and glycolaldehyde fragments are favorably removed from the non-reducing end based on their thermodynamic stabilities. Theoretical FTIR and NMR spectra were calculated to aid in understanding the structures of the oligomeric sugars. Also, the thermodynamics and physical properties of these compounds were estimated using the Group Contribution Method (GCM). In conclusion, these properties are essential in the design of processing technologies and the upgrading of products.

10 SYNTHETIC FUELS↗

Møller–Plesset and Density-Fixed Adiabatic Connections for a Model Diatomic System at Different Correlation Regimes

In recent years, adiabatic connection (AC) interpolations developed within density functional theory (DFT) have been found to provide good performances in the calculation of interaction energies when used with Hartree–Fock (HF) ingredients. The physical and mathematical reasons for such unanticipated performance have been clarified, to some extent, by studying the strong-interaction limit of the Møller–Plesset (MP) AC. In this work, we calculate both the MP and the DFT AC integrand for the asymmetric Hubbard dimer, which allows for a systematic investigation of different correlation regimes by varying two simple parameters in the Hamiltonian: the external potential, Δv, and the interaction strength, U. Notably, we find that, while the DFT AC integrand appears to be convex in the full parameter space, the MP integrand may change curvature twice. Furthermore, we discuss different aspects of the second-order expansion of the correlation energy in each AC, and we demonstrate why the derivative of the λ-dependent density in the MP AC at λ = 0 (i.e., at the HF density) is zero in the model. Concerning the strong-interaction limit of both ACs in the Hubbard dimer setting, we show that the asymptotic value of the MP AC, W ∞ HF , is lower than (or equal to) its DFT analogue, W ∞ KS , if the two are compared at a given density, just like in real space. Furthermore, we also show that this is not always the case if the two quantities are compared at a given external potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

What's Left for a Computational Chemist To Do in the Age of Machine Learning?

Machine learning (ML) has become a central focus of the computational chemistry community. In this paper, I will first discuss my personal history in the field. Then I will provide a broader view of how this resurgence in ML interest echoes and advances upon earlier efforts. Although numerous changes have brought about this latest wave, one of the most significant is the increased accuracy and efficiency of low-cost methods (e.g., density functional theory or DFT) that have made it possible to generate large data sets for ML models. ML has also been used to bypass, guide, or improve DFT. The field of computational chemistry thus finds itself at a crossroads as ML both augments and supersedes traditional efforts. I will present what I believe the role of the computational chemist will be in this evolving landscape, with specific focus on my experience in the development of autonomous workflows in computational materials discovery for open-shell transition-metal chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Efficient machine-learning model for fast assessment of elastic properties of high-entropy alloys

We combined descriptor-based analytical models for stiffness-matrix and elastic-moduli with mean-field methods to accelerate assessment of technologically useful properties of high-entropy alloys, such as strength and ductility. Model training for elastic properties uses Sure-Independence Screening (SIS) and Sparsifying Operator (SO) method yielding an optimal analytical model, constructed with meaningful atomic features to predict target properties. Computationally inexpensive analytical descriptors were trained using a database of elastic properties determined from density functional theory for binary and ternary subsets of Nb-Mo-Ta-W-V refractory alloys. The optimal Elastic-SISSO models, extracted from an exponentially large feature space, give an extremely accurate prediction of target properties, similar to or better than other models, with some verified from existing experiments. Here we also show that electronegativity variance and elastic-moduli can directly predict trends in ductility and yield strength of refractory HEAs, and reveals promising alloy concentration regions.

36 MATERIALS SCIENCE↗

4 d element induced improvement of structural disorder and development of weakly reentrant spin-glass behavior in NiRuMnSn

The pursuit of efficient spin polarization in quaternary Heusler alloys with the general formula XX'YZ (where X, X', and Y are transition metals and Z is a p-block element) has been a subject of significant scientific interest. While previous studies showed that isoelectronic substitution of a 4d element in place of a 3d element in quaternary Heusler alloys improves the half-metallic ferromagnetic characteristics, our research here, on the quaternary Heusler alloy NiRuMnSn, reveals a strikingly different scenario. In this study, we present a detailed structural analysis of the material using x-ray absorption fine structure and neutron diffraction techniques which confirms the formation of a single-phase compound with 50:50 site disorder between Ni and Ru atoms at 4c and 4d sites. Contrary to expectations, our density functional theory calculations suggest a considerable decrease in spin polarization even in the ordered structure. Additionally, we report on the compound's exceptional behavior, displaying a rare reentrant spin-glass property below ~ 60 K, a unique and intriguing feature for quaternary Heusler-type compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Powder X-ray diffraction of altrenogest, C 21 H 26 O 2

The crystal structure of altrenogest has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional techniques. Altrenogest crystallizes in space group P2 1 2 1 2 1 (#19) with a = 7.286 916(16), b = 10.580 333(19), c = 22.266 08(7) Å, V = 1716.671(6) Å 3 , and Z = 4 at 295 K. Thermal expansion between 113 and 295 K is anisotropic. An O–H···O hydrogen bond links the molecules into chains along the c-axis. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

36 MATERIALS SCIENCE↗

Molecular dynamics study on the interfacial properties of mixtures of monomers of polyvinylpyrrolidone (PVP)-based battery binders on graphene and graphite surfaces

This study investigates the behavior of two different mixtures of monomers of polyvinylpyrrolidone (PVP)-based battery binders, polyvinylpyrrolidone:polyvinylidene difluoride (PVP:PVDF) and polyvinylpyrrolidone:polyacrylic acid (PVP:PAA), at graphene and graphite interfaces using classical molecular dynamics simulations. The aim is to identify the best performing monomer binder blend and carbon-based material for the design of battery-optimized energy devices. The PVP:PAA monomer binder blend and graphite are found to have the best interaction energies, densification upon adsorption, and more ordered structure. The adsorption of both monomer binder blends is strongly guided by the higher affinity of PVP and PAA monomeric molecules for the surfaces compared to PVDF. The structure of adsorbed layers of PVP:PVDF monomer binder blend on graphene and graphite develops more quickly than PVP:PAA, indicating faster kinetics. This study complements a previous density functional theory study recently reported by our group and contributes to a better understanding of the nanoscopic features of relevant interfacial regions involving mixtures of monomers of PVP-based battery binders and different carbon-based materials. In conclusion, the effect of a blend of commonly used monomer binders on carbon-based materials is essential for obtaining tightly bound anode and cathode active materials in lithium-ion batteries, which is crucial for designing battery-optimized energy devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding atom probe’s analytical performance for iron oxides using correlation histograms and ab initio calculations

Abstract Field evaporation from ionic or covalently bonded materials often leads to the emission of molecular ions. The metastability of these molecular ions, particularly under the influence of the intense electrostatic field (10 10 Vm −1 ), makes them prone to dissociation with or without an exchange of energy amongst them. These processes can affect the analytical performance of atom probe tomography (APT). For instance, neutral molecules formed through dissociation may not be detected at all or with a time of flight no longer related to their mass, causing their loss from the analysis. Here, we evaluated the changes in the measured composition of FeO, Fe 2 O 3 and Fe 3 O 4 across a wide range of analysis conditions. Possible dissociation reactions are predicted by density-functional theory calculations considering the spin states of the molecules. The energetically favoured reactions are traced on to the multi-hit ion correlation histograms, to confirm their existence within experiments, using an automated Python-based routine. The detected reactions are carefully analyzed to reflect upon the influence of these neutrals from dissociation reactions on the performance of APT for analysing iron oxides.

74 ATOMIC AND MOLECULAR PHYSICS↗

Quantum physical reality of polar-nonpolar oxide heterostructures

Conducting interfaces between polar and nonpolar insulating oxides, e.g., LaAlO 3 /SrTiO 3 , have generated interest for both fundamental physics and oxide-electronics applications. Current understanding is based on an amalgamation of a classical electrostatic model (polar catastrophe model) that was originally derived for semi-infinite solids and quantum density-functional-theory (DFT) results on ultrathin films. Here we report comprehensive DFT calculations that unveil a very different purely quantum physical reality. We show that, for ultrathin polar films, the interfacial dipole does not control the electrostatic potential in the polar film—the surface and interface play equal roles, and the absence or presence of centrosymmetry in the physical LAO film results in different, purely quantum mechanisms for the generation of a conducting interface, neither involving physical-charged transfer. Predictions are made that can be tested and can guide technology development.

36 MATERIALS SCIENCE↗

The Electron Thermal Conductivity of Pu and Zr Substituted $\mathcal{γ}$-U

Uranium alloys are attractive recycled nuclear fuels because of their high thermal conductivity (𝑘) and fissile density. Limited experimental studies of the 𝑘 of U-Pu-Zr alloys in the range of 15 to 20 wt% Pu and 6 to 15 wt% Zr indicate that increasing the content of either Zr or Pu tends to lower 𝑘. However, which element has the greater effect on 𝑘, and the associated mechanisms, remains unclear. Here, in this study, the electron thermal conductivity (𝑘 𝑒 ) of U-Pu-Zr compositions are calculated using density functional theory. The electronic structure is evaluated to understand the effects of plutonium (Pu) and zirconium (Zr) substitution on the 𝑘 𝑒 of 𝛾-U. Alloys of up to 37.5 at. % Pu and 37.5 at. % Zr are examined. Two methods are applied to calculate 𝑘 𝑒 ; we find that the accuracy of each method depends on the electronic and mass similarities between the solute and solvent atoms. Specifically, when the solute atom is similar in electronic structure and mass, the more accurate method is that which employs the electron relaxation time of 𝛾-U, while if the elements are dissimilar, a mixed method that mixes several parameters associated with JNW_S⁢3033426825100132 from each element in the alloy is best. The introduction of all alloying elements decreases 𝑘 𝑒 ; however, in binary compounds, Pu and Zr have different effects. Pu flattens the electronic bands but compensates for this deleterious effect by increasing electron density near the Fermi level. Zr flattens the electronic bands more severely without adding electron density near the Fermi level. Therefore, Zr decreases 𝑘 𝑒 more than Pu in binary compounds. In ternary compounds, the difference between Pu and Zr is minimal due to the phononic change from the large mass change of Zr substitution, even at 12.5 at. %. Thus, we predict that higher loadings of Pu, and potentially other actinides, can be added to U-Pu-Zr compositions for faster recycling of spent fuel without sacrificing 𝑘. We also note that these 𝑘 𝑒 calculation methods can be applied to non-fuel alloys that require 𝑘 𝑒 predictions, such as cladding, heat exchanger, and structural materials.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗