Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional calculation”

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

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

At least 163 records · Page 9

Discovery of single-atom alloy catalysts for CO 2 -to-methanol reaction by density functional theory calculations

The transformations of CO 2 molecules into valuable products are of increasing interest due to the negative impact of anthropogenic CO 2 emissions on global warming. The CO 2 -to-methanol hydrogenation is an economically profitable reaction of carbon fixation, but it still steps away from widespread industrialization because of the lack of efficient and selective catalysts. Recently, single-atom alloy (SAA) catalysts have been developed to work remarkably in CO 2 hydrogenation reactions. Doping isolated single atoms into metallic catalyst can dramatically alter the catalytic performance of the host. Here, we have performed a screening discovery on Ru and 6 RuX (X = Fe, Co, Ni, Cu, Ir and Pt) SAAs using density functional theory (DFT) computations. We considered 13 possible elementary reactions in 4 possible reaction pathways on Ru and all RuX surfaces. In the computed mechanisms, we found that the formation of *H 2 COOH and *HCOO intermediates plays a critical role in determining catalysts’ activities. Doping Co and Pt isolated single atoms into Ru surface can thermodynamically and kinetically facilitate these intermediates formation processes, eventually promoting the production of methanol. The combination of weak binding and enhanced charge redistribution on RuCo and RuPt surfaces give them improved catalytic activities over pure Ru. This work will ultimately facilitate the discovery and development of SAAs for CO 2 to methanol, serving as guidance to experiments and theoreticians alike.

25 ENERGY STORAGE↗

Magnetization distribution in Cu 0.6 Mn 2.4 Ge 2 ferromagnet from polarized and non-polarized neutron powder diffraction aided by density-functional theory calculations

The crystal structure and magnetic properties of Cu 0.6 Mn 2.4 Ge 2 have been re-investigated by a combination of extensive magnetic measurements and neutron scattering experiments, aided by electronic structure calculations. The material is found to be a soft ferromagnet with the ordering temperature T C = 316 K. The magnetocaloric effect evaluated from field-dependent magnetization isotherms is equal to 1.2 J/(kg·K) and 2.5 J/(kg·K) under the maximum applied magnetic field of 2 T and 5 T, respectively. The compound crystallizes in the hexagonal space group P6 3 /mmc. A complex structural disorder necessitated testing of several disorder models against the results of non-polarized and polarized neutron scattering experiments and magnetization measurements. Simulations at the density-functional theory level were also performed to identify the most robust solution that properly described the data observed. The final magnetic structure model reveals non-equal magnetic moments on the Mn1 and Mn2 atoms (2.29(9) µ B and 2.7(1) µ B , respectively) and the presence of vacancies and minor Cu substitution defects in both Mn sites. The work demonstrates how the non-polarized and polarized neutron scattering methods can be combined with electronic structure calculations to establish the microscopic structure of magnetic materials with complex crystallographic disorder.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Search for giant magnetic anisotropy in transition-metal dimers on defected hexagonal boron nitride sheet

Structural and magnetic properties of many transition-metal dimers embedded in a defected hexagonal boron nitride monolayer are investigated through density functional calculations to search for systems with magnetic anisotropy energies (MAEs) larger than 30meV. In particular, Ir–Ir@Dh–BN is found to have both large MAE (∼126 meV) and high structural stability against dissociation and diffusion, and it hence can serve as magnetic unit in spintronics and quantum computing devices. This giant MAE mainly results from the spin orbit coupling and the magnetization of the upper Ir atom, which is in a rather isolated environment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Local order of orthorhombic weberite-type Y 3 TaO 7 as determined by neutron total scattering and density functional theory calculations

Weberite-type oxides are a family of oxides with A 3 BO 7 stoichiometry that can adopt different space groups depending on chemical composition. There is a discrepancy in previous studies as to whether Y 3 TaO 7 is the orthorhombic C222 1 or the Ccmm space group. Here, we describe the short- and long-range structural properties of weberite-type Y 3 TaO 7 using neutron total scattering data, which has a high sensitivity to the oxygen sublattice. Simultaneous analysis of both short- and long-range structural data via conventional Rietveld and small box modelling demonstrates that Y 3 TaO 7 is best modeled as C222 1 . While the Ccmm describes equally well the long-range structure, pair distribution function analysis revealed that the local atomic configuration can be best modelled as C222 1 . This is corroborated by first-principles calculations that confirm the energetic preference of the C222 1 over Ccmm. Further, neutron total scattering data are reported for the first time for the Y 3 TaO 7 weberite-type ceramics.

36 MATERIALS SCIENCE↗

Surface stability and the selection rules of substrate orientation for optimal growth of epitaxial II-VI semiconductors

The surface structures of ionic zinc-blende CdTe (001), (110), (111), and (211) surfaces are systematically studied by first-principles density functional calculations. Based on the surface structures and surface energies, we identify the detrimental twinning appearing in molecular beam epitaxy (MBE) growth of II-VI compounds as the (111) lamellar twin boundaries. To avoid the appearance of twinning in MBE growth, we propose the following selection rules for choosing optimal substrate orientations: (1) the surface should be nonpolar so that there is no large surface reconstructions that could act as a nucleation center and promote the formation of twins; (2) the surface structure should have low symmetry so that there are no multiple equivalent directions for growth. These straightforward rules, in consistent with experimental observations, provide guidelines for selecting proper substrates for high-quality MBE growth of II-VI compounds.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

In situ investigations of laser and thermally modified As{sub 2}S{sub 3} nanolayers: Synchrotron radiation photoelectron spectroscopy and density functional theory calculations

As-deposited, annealed, and in situ As{sub 2}S{sub 3} nanolayers, illuminated by blue (405 nm) and red (650 nm) laser light, were studied using synchrotron radiation photoelectron spectroscopy and DFT electronic structure calculations. Changes in composition and local atomic coordination occurring in the irradiated region of As{sub 2}S{sub 3} films were monitored by analysis of As 3d and S 2p core levels. These studies show that both the thermal treatment and the red laser illumination of As{sub 2}S{sub 3} nanolayers decrease the concentration of homopolar (As-As and S-S) bonds. From the other hand, an increasing concentration of As-rich structural units (s.u.) with a homopolar As-As bond was observed under in situ blue laser illumination of As{sub 2}S{sub 3} nanolayers. Molecular orbital energies were calculated for different As- and S-centered s.u. and used for the interpretation of the core levels and valence band spectra. The surface local structure of the As{sub 2}S{sub 3} nanolayers and its photoinduced transformation are discussed in detail.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Zero-Kelvin Compression Isotherms of the Elements 1 ≤ Z ≤ 92 to 100 GPa

Most of the chemical elements have now been compressed close to or above 100 GPa (1 Mbar) pressure in diamond-anvil cells and the pressure–volume room-temperature isotherms have been measured. We collect these data and use simple lattice-dynamics models to reduce the isotherms to 0 K. We have extended the published work by making new diamond-anvil-cell measurements on Cr and Rh, and by conducting density-functional calculations on the elements Po, At, Rn, Fr, Ra, and Ac. The 0 K data are tabulated for all elements 1 ≤ Z ≤ 92 and 0 ≤ P ≤ 100 GPa. These data are useful for generating wide-range equation of state models and for studying the stability of chemical compounds at high pressure (“Megabar chemistry”). The tables presented here are intended to be reference thermodynamic tables for use in high-pressure research. Further experimental and theoretical work will be needed to extend the tables to higher pressure and to improve accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular-Level Insights into the Reactivity of Siloxane-Based Electrolytes at a Lithium-Metal Anode

A molecular-level understanding of the reactions that occur at the lithium-metal anode/electrolyte interphase is essential to improve the performance of Li–O2 batteries. Experimental and computational techniques are applied to explore the reactivity of tri(ethylene glycol)-substituted trimethylsilane (1NM3), a siloxane-based ether electrolyte, at the lithium-metal anode. In situ/ex situ X-ray diffraction and Fourier-transform infrared spectroscopy studies provide evidence of the formation of lithium hydroxide and lithium carbonates at the anode upon gradual degradation of the metallic lithium anode and the solvent molecules in the presence of oxygen. Density functional calculations performed to obtain a mechanistic understanding of the reductive decomposition of 1NM3 indicate that the decomposition does not require any apparent barrier to produce lithium hydroxide and lithium carbonates when the reduced 1NM3 solvent molecules interact with the oxygen crossing over from the cathode. This study indicates that degradation may be more significant in the case of the 1NM3 solvent, compared to linear ethers such as tetraglyme or dioxalone, because of its relatively high electron affinity. Also, both protection of the lithium metal and prevention of oxygen crossover to the anode are essential for minimizing electrolyte and anode decomposition.

Environmental Molecular Sciences Laboratory↗

Understanding chemical reactions of CO{sub 2} and its isoelectronic molecules with 1-butyl-3-methylimidazolium acetate by changing the nature of the cation: The case of CS{sub 2} in 1-butyl-1-methylpyrrolidinium acetate studied by NMR spectroscopy and density functional theory calculations

NMR spectroscopy ({sup 1}H, {sup 13}C, {sup 15}N) shows that carbon disulfide reacts spontaneously with 1-butyl-1-methylpyrrolidinium acetate ([BmPyrro][Ac]) in the liquid phase. It is found that the acetate anions play an important role in conditioning chemical reactions with CS{sub 2} leading, via coupled complex reactions, to the degradation of this molecule to form thioacetate anion (CH{sub 3}COS{sup −}), CO{sub 2}, OCS, and trithiocarbonate (CS{sub 3}{sup 2−}). In marked contrast, the cation does not lead to the formation of any adducts allowing to conclude that, at most, its role consists in assisting indirectly these reactions. The choice of the [BmPyrro]{sup +} cation in the present study allows disentangling the role of the anion and the cation in the reactions. As a consequence, the ensemble of results already reported on CS{sub 2}-[Bmim][Ac] (1), OCS-[Bmim][Ac] (2), and CO{sub 2}-[Bmim][Ac] (3) systems can be consistently rationalized. It is argued that in system (1) both anion and cation play a role. The CS{sub 2} reacts with the acetate anion leading to the formation of CH{sub 3}COS{sup −}, CO{sub 2}, and OCS. After these reactions have proceeded the nascent CO{sub 2} and OCS interact with the cation to form imidazolium-carboxylate ([Bmim] CO{sub 2}) and imidazolium-thiocarboxylate ([Bmim] COS). The same scenario also applies to system (2). In contrast, in the CO{sub 2}-[Bmim] [Ac] system a concerted cooperative process between the cation, the anion, and the CO{sub 2} molecule takes place. A carbene issued from the cation reacts to form the [Bmim] CO{sub 2}, whereas the proton released by the ring interacts with the anion to produce acetic acid. In all these systems, the formation of adduct resulting from the reaction between the solute molecule and the carbene species originating from the cation is expected. However, this species was only observed in systems (2) and (3). The absence of such an adduct in system (1) has been theoretically investigated using DFT calculations. The values of the energetic barrier of the reactions show that the formation of [Bmim] CS{sub 2} is unfavoured and that the anion offers a competitive reactive channel via an oxygen-sulphur exchange mechanism with the solute in systems (1) and (2)

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reaction Mechanism for the Conversion of γ-Valerolactone (GVL) over a Ru Catalyst: A First-Principles Study

Conversion of biomass-derived γ-valerolactone (GVL) to valuable chemicals has been studied extensively, and understanding the reaction mechanism is very valuable for improving turnover rates and selectivities. Here, we report first-principles density functional calculations, through which we show in detail the reaction pathways of GVL conversion on a Ru(0001) surface, in good agreement with recent experimental results performed on supported Ru catalysts. We find that (i) GVL undergoes a ring-opening reaction rather easily, and (ii) the rate-limiting step toward the formation of 1,4-pentanediol (1,4-PDO) and 2-pentanol (2-PeOH) is the hydrogenation step. The high energy barrier for this step is caused by a strong interaction between Ru and the unsaturated acyl intermediate that is formed after opening the ring. Among all the primary products, formation of 2-butanol (2-BuOH) has the smallest activation barrier, while the slowest step is C–C bond cleavage in the decarbonylation step. We further show that the same acyl intermediate after ring opening of GVL can also be formed by dehydrogenation of 1,4-PDO with a moderate activation barrier, which suggests that both 2-PeOH and 2-BuOH can also be produced in secondary steps.

Bababrik, Reda M.↗

Interfacial Effects on Lithium Superoxide Disproportionation in Li-O 2 Batteries

During the cycling of Li-O-2 batteries the discharge process gives rise to dynamically evolving agglomerates composed of lithium-oxygen nanostructures; however, little is known about their composition. In this paper, we present results for a Li-O-2 battery based on an activated carbon cathode that indicate interfacial effects can suppress disproportionation of a LiO2 component in the discharge product. High-intensity X-ray diffraction and transmission electron microscopy measurements are first used to show that there is a LiO2 component along with Li2O2 in the discharge product. The stability of the discharge product was then probed by investigating the dependence of the charge potential and Raman intensity of the superoxide peak with time. The results indicate that the LiO2 component can be stable for possibly up to days when an electrolyte is left on the surface of the discharged cathode. Density functional calculations on amorphous LiO2 reveal that the disproportionation process will be slower at an electrolyte/LiO2 interface compared to a vacuum/LiO2 interface. The combined experimental and theoretical results provide new insight into how interfacial effects can stabilize LiO2 and suggest that these interfacial effects may play an important role in the charge and discharge chemistries of a Li-O-2 battery.

interface↗

LaMnO 3 Dopants for Efficient Thermochemical Water Splitting Identified by Density Functional Theory Calculations

While ceria is the standard material for two-step water splitting, perovskites are emerging as viable alternatives. In this work, based on the orthorhombic LaMnO 3 supercell, we substitute Li Na K Rb Mg Ca Sr Ba on the A-sites (La sites) and Al Ga In Mg Zn on the B-sites (Mn sites) at a concentration of 37.5%. The range of temperature and oxygen partial pressure at which each composition is stable is predicted. For compositions that are stable in relevant temperature and pressure ranges, the oxygen vacancy formation energies are determined for all of the oxygen vacancy site positions available in the computational supercell. Mg, Ca, Sr, and Ba A-site-substituted LaMnO 3 and Al and In B-site-substituted LaMnO 3 meet these two criteria for candidates in solar-thermal water splitting applications. Finally, oxygen vacancy formation energy can also be controlled by adjusting the doping strategy.

08 HYDROGEN↗

Insights into the physical and chemical properties of a cement-polymer composite developed for geothermal wellbore applications

To isolate injection and production zones from overlying formations and aquifers during geothermal operations, cement is placed in the annulus between well casing and formations. Wellbore cement can fracture due to chemical and physical stress. A polymer-cement (composite) with self-healing properties, stable at high temperatures, was recently developed by our group. Here we present data on the elemental distribution, mineral composition, internal microstructure, and tensile elasticity of this composite and how it compares to cement without polymer (control cement). Entrapped air porosity was similar between control cement (0.44% by volume) and composite (0.22% by volume). Total voids increased with polymer content in the cement (10% and 25% composite), and polymer was found to fill most of these voids. Polymer appears to act as a (cement hydration) retarder increasing the amount of unhydrated silicate minerals in the matrix. Tensile elasticity was improved for composites as compared to control cement, suggesting that the composite materials perform better under mechanical stress. Density-functional calculations reproduce the drop in Young’s modulus by about 60% observed experimentally when polymer is added to the cement. In addition, calculations indicate that the bonding interactions between the cement and polymer remain stable over the range of strain studied.

Rod, Kenton A. (ORCID:0000000318189094)↗

On the applicability of one- and many-electron quantum chemistry models for hydrated electron clusters

We evaluate the applicability of a hierarchy of quantum models in characterizing the binding energy of excess electrons to water clusters. In particular, we calculate the vertical detachment energy of an excess electron from water cluster anions with methods that include one-electron pseudopotential calculations, density functional theory (DFT) based calculations, and ab initio quantum chemistry using MP2 and eom-EA-CCSD levels of theory. The examined clusters range from the smallest cluster size (n = 2) up to nearly nanosize clusters with n = 1000 molecules. The examined cluster configurations are extracted from mixed quantum-classical molecular dynamics trajectories of cluster anions with n = 1000 water molecules using two different one-electron pseudopotenial models. We find that while MP2 calculations with large diffuse basis set provide a reasonable description for the hydrated electron system, DFT methods should be used with precaution and only after careful benchmarking. Strictly tested one-electron psudopotentials can still be considered as reasonable alternatives to DFT methods, especially in large systems. The results of quantum chemistry calculations performed on configurations, that represent possible excess electron binding motifs in the clusters, appear to be consistent with the results using a cavity structure preferring one-electron pseudopotential for the hydrated electron, while they are in sharp disagreement with the structural predictions of a non-cavity model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiple Redox Modes in the Reversible Lithiation of High-Capacity, Peierls-Distorted Vanadium Sulfide

Vanadium sulfide VS4 in the patronite mineral structure is a linear chain compound comprising vanadium atoms coordinated by disulfide anions [S2]2–. 51V NMR shows that the material, despite having V formally in the d1 configuration, is diamagnetic, suggesting potential dimerization through metal–metal bonding associated with a Peierls distortion of the linear chains. This is supported by density functional calculations, and is also consistent with the observed alternation in V–V distances of 2.8 and 3.2 Å along the chains. Partial lithiation results in reduction of the disulfide ions to sulfide S2–, via an internal redox process whereby an electron from V4+ is transferred to [S2]2– resulting in oxidation of V4+ to V5+ and reduction of the [S2]2– to S2– to form Li3VS4 containing tetrahedral [VS4]3– anions. On further lithiation this is followed by reduction of the V5+ in Li3VS4 to form Li3+xVS4 (x = 0.5–1), a mixed valent V4+/V5+ compound. Eventually reduction to Li2S plus elemental V occurs. Despite the complex redox processes involving both the cation and the anion occurring in this material, the system is found to be partially reversible between 0 and 3 V. The unusual redox processes in this system are elucidated using a suite of short-range characterization tools including 51V nuclear magnetic resonance spectroscopy (NMR), S K-edge X-ray absorption near edge spectroscopy (XANES), and pair distribution function (PDF) analysis of X-ray data.

Britto, Sylvia↗