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,315 records · Page 73

Transition Metal Nitrides as Promising Catalyst Supports for Tuning CO/H 2 Syngas Production from Electrochemical CO 2 Reduction

The electrochemical carbon dioxide reduction reaction (CO 2 RR) to produce synthesis gas (syngas) with tunable CO/H 2 ratios has been studied by supporting Pd catalysts on transition metal nitride (TMN) substrates. Combining experimental measurements and density functional theory (DFT) calculations, Pd-modified niobium nitride (Pd/NbN) is found to generate much higher CO and H 2 partial current densities and greater CO Faradaic efficiency than Pd-modified vanadium nitride (Pd/VN) and commercial Pd/C catalysts. In-situ X-ray diffraction identifies the formation of PdH in Pd/NbN and Pd/C under CO 2 RR conditions, whereas the Pd in Pd/VN is not fully transformed into the active PdH phase. DFT calculations show that the stabilized *HOCO and weakened *CO intermediates on PdH/NbN are critical to achieving higher CO 2 RR activity. This work suggests that NbN is a promising substrate to modify Pd, resulting in an enhanced electrochemical conversion of CO 2 to syngas with a potential reduction in precious metal loading.

25 ENERGY STORAGE↗

The LISE package: solvers for static and time-dependent superfluid local density approximation equations in three dimensions

Nuclear implementation of the density functional theory (DFT) is at present the only microscopic framework applicable to the whole nuclear landscape. The extension of DFT to superfluid systems in the spirit of the Kohn-Sham approach, the superfluid local density approximation (SLDA) and its extension to time-dependent situations, time-dependent superfluid local density ap- proximation (TDSLDA), have been extensively used to describe various static and dynamical problems in nuclear physics, neutron star crust, and cold atom systems. In this paper, we present the codes that solve the static and time-dependent SLDA equations in three-dimensional coordinate space without any symmetry restriction. These codes are fully parallelized with the message passing interface (MPI) library and take advantage of graphic processing units (GPU) for accelerating execution. The dynamic codes have checkpoint/restart capabilities and for initial conditions one can use any generalized Slater determinant type of wave function. The code can describe a large number of physical problems: nuclear fission, collisions of heavy ions, the interaction of quantized vor- tices with nuclei in the nuclear star crust, excitation of superfluid fermion systems by time dependent external fields, quantum shock waves, domain wall generation and propagation, the dynamics of the Anderson-Bogoliubov-Higgs mode, dynamics of fragmented condensates, vortex rings dynamics, generation and dynamics of quantized vortices, their crossing and recombinations and the incipient phases of quantum turbulence.

Jin, Shi↗

Atomic Fingerprinting of Heteroatoms Using Noncontact Atomic Force Microscopy

We report immense strides have been made in increasing the resolution of scanning probe microscopy. Noncontact atomic force microscopy (n-AFM) now offers one the ability to characterize and visualize single molecules with subatomic resolution. Specifically, nc-AFM with a carbon monoxide (CO) functionalized tip has the ability to discriminate functional groups (–C≡C–, –CH 2 , –C$=$O, …), although discriminating atomic species often remains as an ongoing challenge. Here, real-space pseudopotentials constructed within density functional theory are employed to accurately simulate nc-AFM images of molecules containing heteroatoms (S, I, and N) within dibenzothiophene (DBT), 2-iodotriphenylene (ITP), acridine (ACR) and ferrous phthalocyanine (FePc). It is found that S and I atoms can be easily identified from C based on their unique features. For N atoms, a use of tip functionalization is proposed to effectively discriminate them from C atoms.

36 MATERIALS SCIENCE↗

Controlling extrapolations of nuclear properties with feature selection

Predictions of nuclear properties far from measured data are inherently inaccurate because of uncertainties in our knowledge of nuclear forces and in our treatment of quantum many-body effects in strongly-interacting systems. While the model bias can be directly calculated when experimental data is available, only an estimate can be made in the absence of such measurements. Current approaches to compute the estimated bias quickly lose predictive power when their input variables are taken far from the training region, resulting in uncontrolled uncertainties in applications such as nucleosynthesis simulations. In this letter, we present a novel technique to identify the input variables of machine learning algorithms that can provide robust estimates of model bias. Our process is based on selecting input variables, or features, based on their probability distribution functions across the entire nuclear chart. We illustrate our approach on the problem of quantifying the model bias in nuclear binding energies calculated with Density Functional Theory (DFT). We prove that building model biases with only Z and N as features leads to highly unreliable extrapolations. Conversely, we show that proper feature selection can systematically improve theoretical predictions without increasing uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hybrid interatomic potential for Sn

To design materials for extreme applications, it is important to understand and predict phase transitions and their influence on material properties under high pressures and temperatures. Atomistic modeling can be a useful tool to assess these behaviors. However, this can be difficult due to the lack of fidelity of the interatomic potentials in reproducing this high pressure and temperature extreme behavior. Here, in this work, a hybrid EAM-R—which is the combination of embedded atom method (EAM) and rapid artificial neural network potential—for Tin (Sn) is described which is capable of accurately modeling the complex sequence of phase transitions between different metallic polymorphs as a function of pressure. This hybrid approach ensures that a basic empirical potential like EAM is used as a lower energy bound. By using the final activation function, the neural network contribution to energy must be positive, assuring stability over the whole configuration space. This implementation has the capacity to reproduce density functional theory results at 6 orders of magnitude slower than a pair potential for molecular dynamics simulation, including elastic and plastic characteristics and relative energies of each phase. Using calculations of the Gibbs free energy, it is demonstrated that the potential precisely predicts the experimentally observed phase changes at temperatures and pressures across the whole phase diagram. At 10.2 GPa, the present potential predicts a first-order phase transition between body-centered tetragonal (BCT) β-Sn and another polymorph of BCT-Sn. This structure transforms into body-centered cubic near the experimentally reported value at 33 GPa. Thus, the Sn potential developed in this paper can be used to study complex deformation mechanisms under extreme conditions of high pressure and strain rates unlike existing potentials. Moreover, the framework developed in this paper can be extended for different material systems with complex phase diagrams.

36 MATERIALS SCIENCE↗

Green’s Function Perspective on the Nonlinear Density Response of Quantum Many-Body Systems

We derive equations of motion for higher order density response functions using the theory of thermodynamic Green’s functions. We also derive expressions for the higher order generalized dielectric functions and polarization functions. Moreover, we relate higher order response functions and higher order collision integrals within the Martin–Schwinger hierarchy. We expect our results to be highly relevant to the study of a variety of quantum many-body systems such as matter under extreme temperatures, densities, and pressures.

Density Functional Theory↗

Efficient Langevin dynamics for “noisy” forces

Efficient Boltzmann-sampling using first-principles methods is challenging for extended systems due to the steep scaling of electronic structure methods with the system size. Additionally, stochastic approaches provide a gentler system-size dependency at the cost of introducing “noisy” forces, which could limit the efficiency of the sampling. When the forces are deterministic, the first-order Langevin dynamics (FOLD) offers efficient sampling by combining a well-chosen preconditioning matrix S with a time-step-bias-mitigating propagator [G. Mazzola and S. Sorella, Phys. Rev. Lett. 118, 015703 (2017)]. However, when forces are noisy, S is set equal to the force-covariance matrix, a procedure that severely limits the efficiency and the stability of the sampling. Here, we develop a new, general, optimal, and stable sampling approach for FOLD under noisy forces. We apply it for silicon nanocrystals treated with stochastic density functional theory and show efficiency improvements by an order-of-magnitude.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Weak magnetic field-dependent photoluminescence properties of lead bromide perovskites

The strong spin–orbit coupling (SOC) in lead halide perovskites, when inversion symmetry is lifted, has provided opportunities for investigating the Rashba effect in these systems. Moreover, the strong orbital moment, which, in turn, impacts the spin-pair in singlet and triplet electronic states, plays a significant role in enhancing the optoelectronic properties in the presence of external magnetic fields in lead halide perovskites. Here, we investigate the effect of weak magnetic fields (<1 T) on the photoluminescence (PL) properties of CsPbBr 3 nanocrystals with and without Ruddlesden–Popper (RP) faults and single crystals of CH 3 NH 3 PbBr 3 . Along with an enhancement in the PL intensity as a function of an external magnetic field, which is observed in both lead bromide perovskites, the PL emission red-shifts in CsPbBr 3 nanocrystals. Density-functional theory calculations of the electronic band-edge in CsPbBr 3 show almost no change in the energy gap as a function of the external magnetic field. The experimental results, thus, suggest the role of mixing of the triplet and singlet excitonic states under weak magnetic fields. This is further deduced from an enhancement in PL lifetimes as a function of the field in CsPbBr 3 ⁠. In CH 3 NH 3 PbBr 3⁠ , an increase in PL intensity is observed under weak magnetic fields; however, no changes in the peak energy or PL lifetimes are observed. The internal magnetic fields due to SOC are characterized for all three samples and found to be the highest for CsPbBr 3 nanocrystals with RP faults.

36 MATERIALS SCIENCE↗

Mechanisms of water oxidation on heterogeneous catalyst surfaces

Water oxidation, an essential step in photosynthesis, has attracted intense research attention. Understanding the reaction pathways at the electrocatalyst/water interface is of great importance for the development of water oxidation catalysts. How the water is oxidized on the electrocatalyst surface by the positive charges is still an open question. This review summarizes current advances in studies on surface chemistry within the context of water oxidation, including the intermediates, reaction mechanisms, and their influences on the reaction kinetics. The Tafel analyses of some electrocatalysts and the rate-laws relative to charge consumption rates are also presented. Moreover, how the multiple charge transfer relies on the intermediate coverage and the accumulated charge numbers is outlined. Finally, the intermediates and rate-determining steps on some water oxidation catalysts are discussed based on density functional theories.

14 SOLAR ENERGY↗

Systematic Study of Solid-State U(VI) Photoreactivity: Long-Lived Radicalization and Electron Transfer in Uranyl Tetrachloride

Reported are the syntheses, structural characterizations and luminescence properties of three novel [UO 2 Cl 4 ] 2- bearing compounds containing substituted 1,1’-dialkyl-4,4’-bipyridinum dications (i.e. viologens). These compounds undergo photoinduced luminescence quenching upon exposure to UV radiation. Further, this reactivity is concurrent with two phenomena: radicalization of the uranyl tetrachloride anion and photoelectron transfer to the viologen which constitutes the formal transfer of one electron from the [UO 2 Cl 4 ] 2- to the viologen species. This behavior is elucidated using electron paramagnetic resonance (EPR) spectroscopy and further probed through a series of characterization and computational techniques including Rehm-Weller analysis, time-dependent density functional theory (TD-DFT), and density of states (DOS). This work provides a systematic study of the photoreactivity of the uranyl unit in the solid state, an under-described aspect of fundamental uranyl chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural and Spectroscopic Properties of Butanediol-Modified Boehmite Materials

Glycoboehmite (GB) materials are synthesized by a solvothermal reaction to form layered aluminum oxyhydroxide (boehmite) modified by intercalated butanediol molecules. These hybrid materials offer a platform to design materials with potentially novel sorption, wetting, and catalytic properties. Several synthetic methods have been used, resulting in different structural and spectroscopic properties, but atomistic detail is needed to determine the interlayer structure to explore the synthetic control of GB materials. Here, in this study, we use classical molecular dynamics (MD) simulations to compare the structural properties of GB interlayers containing chemisorbed butanediol molecules as a function of diol loading. Accompanying quantum (density functional theory, DFT) static calculations and MD simulations are used to validate the classical model and compute the infrared spectra of various models. Classical MD results reveal the existence of two unique interlayer environments at higher butanediol loading, corresponding to smaller (cross-linked) and expanded interlayers. DFT-computed infrared spectra reveal the sensitivity of the aluminol O–H stretch frequencies to the interlayer environment, consistent with the spectrum of the synthesized material. Insight from these simulations will aid in the characterization of the newly synthesized GB materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

In Silico Design of Methyl-Driven Overhauser Dynamic Nuclear Polarization Agents

Overhauser effect (OE) dynamic nuclear polarization (DNP) has drawn attention owing to its enhanced performance at ultrahigh magnetic fields. The lack of design principles has, nevertheless, limited the development of OE polarizing agents when compared to those used for cross-effect DNP. Here, we measured the 19 F OE DNP performance of a series of CF 3 -functionalized Blatter-type radicals. Using density functional theory calculations, we accurately predict the methyl-driven OE performance, paving the way for computer-aided design of optimized OE DNP polarizing agents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nanoconfinement of High Hydrogen-to-Metal Ratio Lanthanum Hydrides in Functionalized Carbon Hosts

Metal hydrides with a high hydrogen content are important for materials-based hydrogen storage and high-temperature superconductivity. Nanoconfinement of metal hydrides in porous hosts is a promising strategy to tune the thermodynamic stability and control the hydrogen-to-metal ratio. However, lanthanum hydride (LaH x ) nanoconfinement in porous materials has been limited due to the challenges associated with isolating and stabilizing nanoparticles of La or La-hydrides. Here we successfully demonstrated the chemical reduction of La­(III) salts to La(0), and subsequent infiltration into pure CMK-3 and nitrogen-doped CMK-3 (NCMK-3) porous carbons. Transmission electron microscopy measurements revealed a uniform distribution of LaH x species within the carbon hosts, while X-ray absorption and photoelectron spectroscopy provided detailed information about the local chemical environment. Sieverts measurements indicate that LaH x @NCMK-3 could desorb up to 0.75 wt % hydrogen, which is higher than non-nitrogen-functionalized CMK-3 (0.43 wt % H). Density Functional Theory and ab initio molecular dynamics calculations indicate that host–guest interaction energies are favorable for porous carbon with nitrogen defects, which is supported by experimental evidence. Moreover, high-pressure synchrotron X-ray diffraction measurements were conducted using a diamond anvil cell up to 60 GPa and reveal that the nitrogen-functionalized nanoporous carbon host favors the formation of higher H:La ratios in the presence of ammonia borane compared to the pure CMK-3 host. This approach could serve as a suitable platform for developing nanoscale superconducting materials at lower pressures and temperatures compared to bulk.

Shivanna, Mohana↗

Computational and Experimental Mechanistic Insights into the Ethanol-to-Butanol Upgrading Reaction over MgO

The mechanism of ethanol upgrading to higher products is still under debate, especially regarding intermediate species and hydrogenation and dehydrogenation steps. In this work, we conducted a combined theoretical and experimental approach to contribute to this discussion. For such, detailed electronic structure density functional theory calculations (aiming at probing density of states, infrared spectra, geometric parameters, charge densities, and reaction energetics) and diffuse reflectance infrared Fourier transform spectroscopy experiments were carried out revealing the relevance of an appropriate combination of reactive surface sites to support the formation of several intermediates that are formed in the C-C coupling over MgO. The roles of Mg and O sites were also studied under an electronic perspective and different geometrical arrangements. We found that a kink configuration was the most adequate for ethanol to 1-butanol upgrading. Our calculations also gave us arguments to propose distinct reaction routes, whose mutual predominance would depend upon reaction temperature. At temperatures up to 573 K, the so-called β-route, which goes through scission of a Cβ-H bond and formation of an oxametallacycle-like intermediate, would dominate the coupling, whereas at higher temperatures, up to 673 K, a more usual Guerbet mechanism, via an aldol coupling step and then consecutive hydrogenations, would be expected. Here, the theoretical conclusions were followed by a careful experimental strategy using sequential experimental planning techniques in order to estimate accurate parameters with the lowest possible experimental load. Information from these different sources were coupled to develop a mathematical model for the rate of the ethanol upgrading reaction, using a Langmuir-Hinshelwood-Hougen-Watson approach. The developed and statistically validated model adequately described the experimental data at 673 K and 1.1 bar total pressure for ethanol partial pressures in the range from 0 to 20 kPa.

09 BIOMASS FUELS↗

Regioselective On-Surface Synthesis of [3]Triangulene Graphene Nanoribbons

The integration of low-energy states into bottom-up engineered graphene nanoribbons (GNRs) is a robust strategy for realizing materials with tailored electronic band structure for nanoelectronics. Low-energy zero-modes (ZMs) can be introduced into nanographenes (NGs) by creating an imbalance between the two sublattices of graphene. This phenomenon is exemplified by the family of [n]triangulenes (n ϵ $\mathbb{N}$). Here, we demonstrate the synthesis of [3]triangulene-GNRs, a regioregular one-dimensional (1D) chain of [3]triangulenes linked by five-membered rings. Hybridization between ZMs on adjacent [3]triangulenes leads to the emergence of a narrow band gap, E g,exp ~ 0.7 eV, and topological end states that are experimentally verified using scanning tunneling spectroscopy. Tight-binding and first-principles density functional theory calculations within the local density approximation corroborate our experimental observations. Our synthetic design takes advantage of a selective on-surface head-to-tail coupling of monomer building blocks enabling the regioselective synthesis of [3]triangulene-GNRs. Detailed ab initio theory provides insights into the mechanism of on-surface radical polymerization, revealing the pivotal role of Au-C bond formation/breakage in driving selectivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental determination of the valence band offsets of ZnGeN 2 and (ZnGe) 0.94 Ga 0.12 N 2 with GaN

A predicted type-II staggered band alignment with an approximately 1.4 eV valence band offset at the ZnGeN 2 /GaN heterointerface has inspired novel band-engineered III-N/ZnGeN 2 heterostructure-based device designs for applications in high performance optoelectronics. We report on the determination of the valence band offset between metalorganic chemical vapor deposition grown (ZnGe) 1- x Ga 2 x N 2 , for x = 0 and 0.06, and GaN using x-ray photoemission spectroscopy. The valence band of ZnGeN 2 was found to lie 1.45–1.65 eV above that of GaN. This result agrees well with the value predicted by first-principles density functional theory calculations using the local density approximation for the potential profile and quasiparticle self-consistent GW calculations of the band edge states relative to the potential. For (ZnGe) 0.94 Ga 0.12 N 2 the value was determined to be 1.29 eV, ~10%–20% lower than that of ZnGeN 2 . The experimental determination of the large band offset between ZnGeN 2 and GaN provides promising alternative solutions to address challenges faced with pure III-nitride-based structures and devices.

Physics↗

Investigating the effects of local environment on nitrogen vacancies in high-entropy metal nitrides

High-entropy metal nitrides are an important material class in a variety of applications, and the role of nitrogen vacancies is of great importance for understanding their stability and mechanical properties. Here, we study six different high-entropy nitrides with eight different metal species to build a predictive model of the nitrogen-vacancy formation energy. We construct sets of supercells that maximize the number of unique nitrogen environments for a given chemistry, and then use density-functional theory to calculate the energy density for all nitrogen sites, and the vacancy formation energies for the highest, lowest, and a median subset based on the energy densities. The energy density of nitrogen sites correlates with the vacancy formation energies, for binary, ternary, and high-entropy nitrides. A linear regression model predicts the vacancy formation energies using only the nearest-neighbor composition; across our eight metals, we find the largest vacancy formation energies next to Hf, then Zr, Ti, V, Cr, Ta, Nb, and the lowest near Mo. Additionally, we see that binary nitride data show qualitatively similar vacancy formation energy trends for high-entropy nitrides; however, the binary data alone are insufficient to predict the complex nitride behavior. Our model is both predictive and easily interpretable, and correlates with experimental data.

DeSilva, Charith R. [Univ. of Illinois at Urbana-C↗

Statistics of some atmospheric turbulence records relevant to aircraft response calculations

Methods for characterizing atmospheric turbulence are described. The methods illustrated include maximum likelihood estimation of the integral scale and intensity of records obeying the von Karman transverse power spectral form, constrained least-squares estimation of the parameters of a parametric representation of autocorrelation functions, estimation of the power spectra density of the instantaneous variance of a record with temporally fluctuating variance, and estimation of the probability density functions of various turbulence components. Descriptions of the computer programs used in the computations are given, and a full listing of these programs is included.

Mark, W. D.↗