Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density functional theory simulations”

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

Molecular growth upon ionization of van der Waals clusters containing HCCH and HCN is a pathway to prebiotic molecules

The growth mechanisms of organic molecules in an ionizing environment such as the interstellar medium are not completely understood. Here we examine by means of ab initio molecular dynamics (AIMD) simulations and density functional theory (DFT) computations the possibility of bond formation and molecular growth upon ionization of van der Waals clusters of pure HCN clusters, and mixed clusters of HCN and HCCH, both of which are widespread in the interstellar medium. Ionization of van der Waals clusters can potentially lead to growth in low temperature and low-density environments. Our results show, that upon ionization of the pure HCN clusters, strongly bound stable structures are formed that contain N∙∙∙H bonds, and growth beyond pairwise HCN molecules is seen only in a small percentage of cases. In contrast, mixed clusters, where HCCH is preferentially ionized over HCN, can grow up to 3 or 4 units long with new carbon–carbon and carbon–nitrogen covalent bonds. Moreover, cyclic molecules formed, such as the radical cation of pyridine, which is a prebiotic molecule. The results presented here are significant as they provide a feasible pathway for molecular growth of small organic molecules containing both carbon and nitrogen in cold and relatively denser environments such as in dense molecular clouds but closer to the photo-dissociation regions, and protoplanetary disks. In the mechanism we propose, first, a neutral van der Waals cluster is formed. Once the cluster is formed it can undergo photoionization which leads to chemical reactivity without any reaction barrier.

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

Chemisorbed oxygen or surface oxides steer the selectivity in Pd electrocatalytic propene oxidation observed by operando Pd L-edge X-ray absorption spectroscopy

Controlled electrochemical oxidation of hydrocarbons to desired products is an attractive approach in catalysis. Here we study the electrochemical propene oxidation under operando conditions using Pd L-edge X-ray absorption spectroscopy (XAS) as a sensitive probe to elucidate surface processes occurring during catalysis. Together with ab initio multiple-scattering calculations, our XAS results enable assignment of characteristic changes of the Pd L-edge intensity and energy position in terms of a mechanistic understanding of the selective oxidation of propene. The results, supported by electrochemical density functional theory DFT simulations, show that in the potential range of 0.8–1.0 V vs. the reversible hydrogen electrode (RHE), selective oxidation of propene to acrolein and acrylic acid occurs on the metallic Pd surface. Here, these reactions are proposed to proceed via the Langmuir–Hinshelwood mechanism. In contrast, for the potential range of 1.1–1.3 V vs. RHE, selective oxidation of propene to propylene glycol takes place on a Pd oxide surface.

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

Data-driven design of novel halide perovskite alloys

The great tunability of the properties of halide perovskites presents new opportunities for optoelectronic applications as well as significant challenges associated with exploring combinatorial chemical spaces. Here, in this work, we develop a framework powered by high-throughput computations and machine learning for the design and prediction of mixed cation halide perovskite alloys. In a chemical space of ABX 3 perovskites with a selected set of options for A, B, and X species, pseudo-cubic structures with B-site mixing are simulated using density functional theory (DFT) and several properties are computed, including stability, lattice constant, band gap, vacancy formation energy, refractive index, and optical absorption spectrum, using both semi-local and hybrid functionals. Neural networks (NN) are used to train predictive models for every property using tabulated elemental properties of A, B, and X site species as descriptors. Starting from a DFT dataset of 229 points, we use the trained NN models to predict the structural, energetic, electronic, and optical properties of an extensive dataset of 17955 compounds, and perform high-throughput screening in terms of stability, band gap, and defect tolerance, to obtain 392 promising compounds that are ranked as potential absorbers according to their photovoltaic figure of merit. Compositional trends in the screened set of attractive mixed cation halide perovskites are revealed and additional computations are performed on selected compounds. The data-driven design framework developed here is promising for designing novel mixed compositions and can be extended to a wider perovskite chemical space in terms of A, B, and X species, different kinds of mixing at the A, B, or X sites, non-cubic phases, and other properties of interest.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elucidating correlated defects in metal organic frameworks using theory-guided inelastic neutron scattering spectroscopy

Metal organic frameworks (MOFs) that incorporate metal oxide cluster nodes, exemplified by UiO-66, have been widely studied, especially in terms of their deviations from the ideal, defect-free crystalline structures. Although defects such as missing linkers, missing nodes, and the presence of adventitious synthesis-derived node ligands (such as acetates and formates) have been proposed, their exact structures remain unknown. Previously, it was demonstrated that defects are correlated and span multiple unit cells. The highly specialized techniques used in these studies are not easily applicable to other MOFs. Thus, there is a need to develop new experimental and computational approaches to understand the structure and properties of defects in a wider variety of MOFs. Here, we show how low-frequency phonon modes measured by inelastic neutron scattering (INS) spectroscopy can be combined with density functional theory (DFT) simulations to provide unprecedented insights into the defect structure of UiO-66. We are able to identify and assign peaks in the fingerprint region (<100 cm –1 ) which correspond to phonon modes only present in certain defective topologies. Specifically, this analysis suggests that our sample of UiO-66 consists of predominantly defect-free fcu regions with smaller domains corresponding to a defective bcu topology with 4 and 2 acetate ligands bound to the Zr 6 O 8 nodes. Importantly, the INS/DFT approach provides detailed structural insights (e.g., relative positions and numbers of acetate ligands) that are not accessible with microscopy-based techniques. Furthermore, the quantitative agreement between DFT simulations and the experimental INS spectrum combined with the relative simplicity of sample preparation, suggests that this methodology may become part of the standard and preferred protocol for the characterization of MOFs, and, in particular, for elucidating the structure defects in these materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Screening binary alloys for electrochemical CO 2 reduction towards multi-carbon products

Electrochemical reduction of CO 2 (eCO 2 R) to high-value chemicals presents an attractive approach for utilizing CO 2 . Copper (Cu) is presently the only electrocatalyst that fulfills this purpose with notable activity, but selectivity remains a problem. To identify catalysts for eCO 2 R with high selectivity towards multicarbon (C 2(+) ) products, we explore binary systems composed of strongly and weakly CO binding metals alloyed with Cu, Fe, Co, Ni, and Pd. A total number of 142 alloys with two commonly studied configurations, L1 2 and L1 0 , are simulated with density functional theory (DFT). We leverage recent progress in the atomistic understanding of the eCO 2 R mechanism and use the binding energies of CO* and C* as descriptors when screening for C 2(+) selectivity. We evaluate the stability of the binary alloys by analyzing the formation energy of the clean alloy surfaces. Our theoretical screening identifies about 16 Cu-based alloys and 18 non-Cu based alloys with optimal C 2(+) selective properties for eCO 2 R. For the non-Cu based binary alloys, the p-block elements play an important role in tuning the C* and CO* adsorption energies. In terms of stability, most of the Cu-based systems alloyed with metals that exhibit strong CO* binding are unstable. Ni-based alloys are more stable than the Co-based alloys followed by the Fe-based alloys, and all the Pd-based alloys are stable. In general, the L1 0 structural Fe, Co, Ni, and Pd-based alloys are more stable than the corresponding L1 2 alloys. Our approach identifies materials known to have good C 2(+) selectivity, but it also proposes several other promising materials that have not previously been tested for eCO 2 R.

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

Nonideal mixing effects in warm dense matter studied with first-principles computer simulations

Here, we study nonideal mixing effects in the regime of warm dense matter (WDM) by computing the shock Hugoniot curves of BN, MgO, and MgSiO 3 . First, we derive these curves from the equations of state (EOS) of the fully interacting systems, which were obtained using a combination of path integral Monte Carlo calculations at high temperature and density functional molecular dynamics simulations at lower temperatures. We then use the ideal mixing approximation at constant pressure and temperature to rederive these Hugoniot curves from the EOS tables of the individual elements. We find that the linear mixing approximation works remarkably well at temperatures above ~2 × 10 5 K, where the shock compression ratio exceeds ~3.2. The shape of the Hugoniot curve of each compound is well reproduced. Regions of increased shock compression, which emerge because of the ionization of L and K shell electrons, are well represented, and the maximum compression ratio of the Hugoniot curves is reproduced with high precision. Some deviations are seen near the onset of the L shell ionization regime, where ionization equilibrium in the fully interacting system cannot be well reproduced by the ideal mixing approximation. This approximation also breaks down at lower temperatures, where chemical bonds play an increasingly important role. However, the results imply that the equilibrium properties of binary and ternary mixtures in the regime of WDM can be derived from the EOS tables of the individual elements. This significantly simplifies the characterization of binary and ternary mixtures in the WDM and plasma phases, which otherwise requires large numbers of more computationally expensive first-principles computer simulations.

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

Decoupling between Shockley partials and stacking faults strengthens multiprincipal element alloys

We report that mechanical properties are fundamental to structural materials, where dislocations play a decisive role in describing their mechanical behavior. Although the high-yield stresses of multiprincipal element alloys (MPEAs) have received extensive attention in the last decade, the relation between their mechanistic origins remains elusive. Our multiscale study of density functional theory, atomistic simulations, and high-resolution microscopy shows that the excellent mechanical properties of MPEAs have diverse origins. The strengthening effects through Shockley partials and stacking faults can be decoupled in MPEAs, breaking the conventional wisdom that low stacking fault energies are coupled with wide partial dislocations. This study clarifies the mechanistic origins for the strengthening effects, laying the foundation for physics-informed predictive models for materials design.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Plasmon spectroscopy for the determination of Ti 3 C 2 T x MXene few layer stacks architecture

Like many 2D materials, numerous properties of MXene multilayers, and especially the most popular one Ti 3 C 2 T x , have been shown to significantly depend on their architecture, i.e. the number of layers and interlayer distance. These structural parameters are thus key elements to be characterized for the analysis of MXene properties. Focusing on valence electron energy-loss spectroscopy (VEELS) as performed in a transmission electron microscope (TEM), and using density functional theory (DFT) simulations, we here analyze the layer dependent large changes in the VEEL spectra of Ti 3 C 2 T x multilayers as a probe of their total thickness, and emphasize the bulk plasmon energy sensitivity to interlayer distance. Together these findings allow to directly quantify the absolute number of layers in a Ti 3 C 2 T x stack up to ~10 nm thickness and give access to interlayer distance modifications with sub-angström sensitivity, evidencing VEELS as a powerful method for the characterization of MXene multilayers on the nanometer scale. We expect these results to be relevant for the study of structure/properties correlations in this class of materials, especially with the development of in situ or environmental TEM experiments.

Materials Science↗
DOE OSTI2022

All-electron, density-functional-based method for angle-resolved tunneling ionization in the adiabatic regime

Here we develop and test a method that integrates many-electron weak-field asymptotic theory (ME-WFAT) [O. I. Tolstikhin, L. B. Madsen, and T. Morishita, Phys. Rev. A 89, 013421 (2014)] in the integral representation (IR) into the density-functional-theory (DFT) framework. In particular, we present modifications of the integral formula in the IR ME-WFAT to incorporate the potential terms unique to DFT. By solving an adiabatic rate equation for the angle-resolved ionization yield in our DFT-based ME-WFAT method, we show that the results are in excellent agreement with those of real-time time-dependent density-functional-theory (TDDFT) simulations for NO, OCS, C⁢H 3 ⁢Br, and C⁢H 3 ⁢Cl interacting with one- and two-color laser fields with a fundamental wavelength of 800 nm. This agreement is significant because the WFAT calculations take only a small fraction of the time of full TDDFT calculations. These results suggest that in the wavelength region commonly used in strong-field experiments (800 nm and longer), our DFT-based WFAT treatment can be used to rapidly screen for the ionization properties of a large number of molecules as a function of alignment or orientation between the molecule and the strong field.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

First-principles coupled cluster theory of the electronic spectrum of transition metal dichalcogenides

The electronic properties of two-dimensional transition metal dichalcogenides (2D TMDs) have attracted much attention during the last decade. Here we show how a diagrammatic ab initio coupled cluster singles and doubles (CCSD) treatment paired with a careful thermodynamic limit extrapolation in two dimensions can be used to obtain converged band gaps for monolayer materials in the MoS 2 family. We find CCSD gaps to lie in the upper range of the spread of $\textit{GW}$ approximation based on density functional theory (DFT) simulations, and also find slightly higher effective hole masses compared to previous reports. We also investigate the ability of CCSD to describe trion states, finding a reasonable qualitative structure, but poor excitation energies due to the lack of screening of three-particle excitations in the effective Hamiltonian. Our study provides an independent high-level benchmark of the role of many-body effects in 2D TMDs and showcases the potential strengths and weaknesses of diagrammatic coupled cluster approaches for realistic materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Prediction of high strain polar phases in antiferroelectric PbZrO 3 from multiscale approach

PbZrO 3 is regarded as the first antiferroelectric and currently under intense reexamination. We propose a multiscale approach that combines classical and first principles density functional theory based simulations to explore polar phases in this material. Application of this approach to PbZrO 3 predicts three polar phases, none of which is the common $R3c$. The lowest energy polar phase - $Cc$ - is metastable at zero Kelvin but can be stabilized by application of electric field. Its structural and polarization response to the electric field is in good agreement with experimental data from the literature. The other two phases are $Ima2$ and $I4cm$ of orthorhombic and tetragonal symmetry, respectively, and predicted to stabilize at finite temperatures and under application of larger electric fields. These two phases exhibit very large strain which is of technological significance.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Phase Switching as the Origin of Large Piezoelectric Response in Organic-Inorganic Perovskites: A First-Principles Study

Piezoelectrics are critical functional components of many practical applications such as sensors, ultrasonic transducers, actuators, medical imaging, telecommunications. So far, the best performing piezoelectrics are ferroelectric ceramics, many of which are toxic, heavy, hard, and cost-ineffective. Recently, a groundbreaking discovery of extraordinarily large piezoelectric coefficients in the family of organic-inorganic perovskites gave a hope for a cheaper, environmentally friendly, inexpensive, light-weight and flexible alternative. However, the origin of such response in organic-inorganic ferroelectrics whose spontaneous polarization is an order of magnitude smaller than for inorganic counterparts remains unclear. In this study we employ first-principles simulations to predict that the mechanism associated with large piezoelectric constants is of extrinsic origin and associated with switching between the stable phase and previously overlooked energetically competitive metastable phase that can be stabilized by the external stress. Here, the phase switching changes the polarization direction and, therefore, gives origin to large piezoelectric response, similar to PbZr$_{1-x}$Ti$_x$O$_3$ near morphotropic phase boundary. Existence of such metastable phases is likely to manifest as the dynamical molecular disorder above the Curie temperature and, therefore, could be intrinsic to the entire family of organic-inorganic ferroelectrics with such disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Understanding the lattice thermal conductivity of SrTiO 3 from an ab initio perspective

Here, we present a detailed analysis of the structure dependence of the lattice thermal conductivity of SrTiO 3 . We have used both ab initio molecular dynamic simulations and density functional theory calculations to decouple the effect of different structural distortions on the thermal conductivity. We have identified two main mechanisms for tuning the thermal conductivity when a distortion is applied. First, the modification of the acoustic-modes energy dispersion when a change in the lattice parameters is imposed and second, the low energy polar modes. In particular and counterintuitively, we have found that an increase in the angle of the oxygen octahedral rotations increases the thermal conductivity due to its coupling to these polar modes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Atomistic simulations for investigation of substrate and salt effects on lipid in-source fragmentation in secondary ion mass spectrometry: A follow-up study

In-source fragmentation (ISF) poses a significant challenge in secondary ion mass spectrometry (SIMS). These fragment ions increase the spectral complexity and can lead to incorrect annotation of fragments as intact species. The presence of salt that is ubiquitous in biological samples can influence the fragmentation and ionization of analytes in a significant manner, but their influences on SIMS have not been well characterized. To elucidate the effect of substrates and salt on ISF in SIMS, we have employed experimental SIMS in combination with atomistic simulations of a sphingolipid on a gold surface with various NaCl concentrations as a model system. Our results revealed that a combination of bond dissociation energy and binding energy between N-palmitoyl-sphingomyelin and a gold surface is a good predictor of fragment ion intensities in the absence of salt. However, ion-fragment interactions play a significant role in determining fragment yields in the presence of salt. Additionally, the charge distribution on fragment species may be a major contributor to the varying effects of salt on fragmentation. This study demonstrates that atomistic modeling can help predict ionization potential when salts are present, providing insights for more accurate interpretations of complex biological spectra.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Thermodynamic phase transitions of nematic order in magnetic helices

A nematic phase lacks translation order but has orientational order. Nematic phases have been discovered in a variety of systems, including liquid crystals, correlated materials, and superconductors. Here, we report on a magnetic nematic phase, where the basis components are composed of magnetic helices. We directly probed the order parameters associated with the magnetic helices using resonant soft x-ray scattering and find two distinct nematic phases with complex spatiotemporal signatures. Using x-ray correlation spectroscopy, we find that near the phase boundary between the two nematic phases, fluctuations coexist on multiple disparate timescales. Our micromagnetic simulations and density functional theory calculations show that the fluctuations occur concomitantly with a reorientation of the magnetic helices, indicating spontaneous symmetry breaking and the emergence of additional degrees of freedom. Our results provide a framework for characterizing exotic phases that can be extended to a broad class of physical systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The Critical Stack Pressure to Alter Void Generation at Li/Solid-Electrolyte Interfaces during Stripping

The lithium stripping process generates vacancies, which may accumulate as voids and lead to uneven current distribution and dendrite growth in the following plating cycles. A stack pressure is typically required during stripping, but how to optimize the stack pressure is not clear. In this work, extremely lithiophilic Li/Li 2 O and lithiophobic Li/LiF interfaces were used to reveal the combining effect of interface interaction and stack pressure induced lithium creep on the stripping critical current density (CCD). A multiscale simulation scheme with Density Functional Theory (DFT), kinetic Monte Carlo (KMC) simulations, and an analytical model was developed. The analytical model predicted lithiophobic interfaces require a higher stack pressure than lithiophilic interfaces to reach the same CCD. The KMC simulations also showed higher stack pressure is needed at lithiophobic interfaces to accelerate Li vacancy diffusion into the bulk and maintain a flat surface. This stack pressure needs to be high enough to alter the Li forward-and-backward hopping barriers at the interface. This multiscale simulation scheme illustrates the importance to include the chemical-mechanical effects during Li stripping morphology evolution. It can be used to design ideal interlayer coating materials to maintain a flat Li surface during cycling.

25 ENERGY STORAGE↗
DOE OSTI2020

Computational and Experimental Characterization of Intermediate Amorphous Phases in Geological Materials

In the subsurface, MgO engineered barriers are employed at the Waste Isolation Pilot Plant (WIPP), a transuranic waste repository near Carlsbad, NM. During service, the MgO will be exposed to high concentration brine environments and may form stable intermediate phases that can alter the barriers effectiveness. Here, MgO was aged in water and three different brine solutions. X-ray diffraction (XRD) and 1 H nuclear magnetic resonance (NMR) analysis were performed to identify the formation of secondary phases. After aging, ~4% of the MgO was hydrated and fine-grained powders resulted in greater loss of crystallinity than hard granular grains. 1 H magic angle spinning (MAS) NMR spectra resolved minor phases not visible in XRD, indicating that diverse 1 H environments are present along with Mg(OH) 2 . Density functional theory (DFT) simulations for several proposed Mg-O-H, Mg-CI-O-H, and Na-O-H containing phases were performed to index peaks in the experimental 1 H MAS NMR spectra. While proposed intermediate crystal structures exhibited overlapping 1 H NMR peaks, Mg-O-H intermediates were attributed to the growth of the 1.0-0.0ppm peak while the Mg-CI-O-H structures contributed to the 2.5- 5.0ppm peak in the chloride containing brines. Overall, NMR analysis of aged MgO indicates the formation of a range of possible intermediate structures that cannot be resolved with XRD analysis alone.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Carbon Capture in Novel Porous Liquids

Direct air capture (DAC) of CO 2 is one of the negative emission technologies under development to limit the impacts of climate change. The dilute concentration of CO 2 in the atmosphere (~400 ppm) requires new materials for carbon capture with increased CO 2 selectivity that is not met with current materials. Porous liquids (PLs) are an emerging material that consist of a combination of solvents and porous hosts creating a liquid with permanent porosity. PLs have demonstrated excellent CO 2 selectivity, but the features that control how and why PLs selectively capture CO 2 is unknown. To elucidate these mechanisms, density functional theory (DFT) simulations were used to investigate two different PLs. The first is a ZIF-8 porous host in a water/glycol/2-methylimidazole solvent. The second is the CC13 porous organic cage with multiple bulky solvents. DFT simulations identified that in both systems, CO 2 preferentially bound in the pore window rather than in the internal pore space, identifying that the solvent-porous host interface controls the CO 2 selectivity. Additionally, SNL synthesized ZIF-8 based PL compositions. Evaluation of the long-term stability of the PL identified no change in the ZIF-8 crystallinity after multiple agitation cycles, identifying its potential for use in carbon capture systems. Through this project, SNL has developed a fundamental understanding of solvent-host interactions, as well as how and where CO 2 binds in PLs. Through these results, future efforts will focus not on how CO 2 behaves inside the pore, but on the porous host-solvent interface as the driving force for PL stability and CO 2 selectivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗