Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT 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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Achieving accurate energetics beyond (semi-)local density functional theory: Illustrated with transition metal disulfides, Cu 2 ZnSnS 4 , and Na 3 PS 4 related semiconductors

Understanding, designing, and processing functional metal sulfides are significant challenges partially because of the high temperatures and pressures and the number of secondary phases encountered in these complex systems. In particular, the underlying thermochemical properties are not well understood that would allow prediction of equilibrium conditions and driving forces. In addition, obtaining accurate values for the energetics of metal sulfide systems is far from complete, suggesting application of density functional theory (DFT) calculations. Here, the results of an examination of 69 phases by DFT using 12 exchange-correlation (X-C) functionals indicate that (i) the key source of error in predicting the Gibbs energy of a phase is the enthalpy calculated at 0 K rather than entropy at finite temperatures from phonon calculations and (ii) an improved prediction of the thermodynamic properties at 0 K relies on the selected nonlocal X-C functional such as the hybrid potential. Regarding metal sulfides, we conclude from the present DFT results that (1) the secondary phase Cu 2 ZnSn 3 S 8 , associated with the desired photovoltaic material Cu 2 ZnSnS 4 , is not stable at 0 K, but it becomes slightly stable with increasing temperature (i.e., >800 K), primarily due to the vibrational entropy, which makes it difficult to be detected in a typical thin-film growth process; (2) the hybrid X-C functional improves the predicted energetics for most of the layered transition-metal disulfides such as TiS 2 , MoS 2 , and WS 2 , but not for the nonlayered RuS 2 , OsS 2 , and IrS 2 as well as the layered PdS 2 ; and (3) the formation of the solid-state electrolyte Na3PS4 is thermodynamically favored. We further conclude that accurate energetics as a function of temperature for the materials of interest is feasible to be achieved beyond the semilocal DFT calculations with the key being enthalpy predicted at 0 K.

36 MATERIALS SCIENCE↗
DOE OSTI2026

HfTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-titanium (Hf-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ti. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2026

TiV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys titanium-vanadium (Ti-V). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Ti and V. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2025

HfZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-zirconium (Hf-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

HfTa_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-tantalum (Hf-Ta). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ta. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect BCC lattice sites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,040 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Sub-4 nm PtZn Intermetallic Nanoparticles for Enhanced Mass and Specific Activities in Catalytic Electrooxidation Reaction

Atomically ordered intermetallic nanoparticles (iNPs) have sparked considerable interest in fuel cell applications by virtue of their exceptional electronic and structural properties. However, the synthesis of small iNPs in a controllable manner remains a formidable challenge because of the high temperature generally required in the formation of intermetallic phases. Here we report a general method for the synthesis of PtZn. iNPs (3.2 +/- 0.4 nm) on multiwalled carbon nanotubes (MWNT) via a facile and capping agent free strategy using a sacrificial mesoporous silica (mSiO(2)) shell. The as-prepared PtZn iNPs exhibited ca. 10 times higher mass activity in both acidic and basic solution toward the methanol oxidation reaction (MOR) compared to larger PtZn iNPs synthesized on MWNT without the mSiO2 shell. Density functional theory (DFT) calculations predict that PtZn systems go through a "non-CO" pathway for MOR because of the stabilization of the OH* intermediate by Zn atoms, while a pure Pt system forms highly stable COH* and CO* intermediates, leading to catalyst deactivation. Experimental studies on the origin of the backward oxidation peak of MOR coincide well with DFT predictions. Moreover, the calculations demonstrate that MOR on smaller PtZn iNPs is energetically more favorable than larger iNPs, due to their high density of corner sites and lower-lying energetic pathway. Therefore, smaller PtZn iNPs not only increase the number but also enhance the activity of the active sites in MOR compared with larger ones. This work opens a new avenue for the synthesis of small iNPs with more undercoordinated and enhanced active sites for fuel cell applications.

Qi, Zhiyuan↗
DOE OSTI2019

Probing Surface Chemistry at an Atomic Level: Decomposition of 1-Propanethiol on GaP(001) (2 × 4) Investigated by STM, XPS, and DFT

The adsorption and decomposition mechanisms for 1-propanethiol on a Ga-rich GaP(001) (2 × 4) surface are investigated at an atomic level using scanning tunneling microscopy (STM), X-ray photoelectron spectroscopy, and density functional theory (DFT) calculations. Using a combination of experimental and theoretical tools, we probe the detailed structures and energetics of a series of reaction intermediates in the thermal decomposition pathway from 130 to 773 K. At 130 K, the propanethiolate adsorbates are observed at the edge gallium sites, with the thiolate–Ga bonding configuration maintained up to 473 K. Further decomposition produces two new surface features, Ga–S–Ga and P-propyl species at 573 K. Finally, S-induced (1 × 1) and (2 × 1) reconstructions are observed at 673–773 K, which are reportedly associated with arrays of surface Ga–S–Ga bonds and subsurface diffusion of S. To understand the observed site-selectivity on the hydrogen dissociation of the thiol molecule at 130 K, the two most likely dissociation pathways (Ga–P vs Ga–Ga dimer sites) are investigated using DFT Gibbs energy calculations. While the theory predicts the kinetic advantage for the dissociation reaction occurring on the Ga–P dimer (Lewis acid–base combination), we only observed dissociation products on the Ga–Ga dimer (Lewis acid). Here, the DFT calculations clarify that the reversible thiolate diffusion along the Ga dimer row prevents recombinative desorption, which is probable on the Ga–P dimer. Together with experimental and theoretical results, we suggest a thermal decomposition mechanism for the thiol molecule with atomic-level structural details.

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

First-principles investigation of lanthanides diffusion in HCP zirconium via vacancy-mediated transport

The diffusion of lanthanide fission products plays an important role in the growth of the fuel-cladding chemical interaction (FCCI) region in metallic fuels. The use of a Zr interdiffusion barrier may mitigate the transport of lanthanides from the fuel to the cladding, but the efficacy of such a liner is not yet known. In this paper, the stability and vacancy-mediated diffusion of La, Ce, Pr, and Nd in hexagonal close-packed (HCP) Zr is investigated via density functional theory (DFT) calculations and self-consistent mean field (SCMF) analysis. DFT is used to calculate the formation, binding, and migration energies of vacancies and vacancy-solute pairs. The DFT energetics are used in the KineCluE code to calculate the Onsager transport coefficients. La is found to be the fastest diffusing species in HCP Zr and experiences an almost isotropic diffusion behavior. The other three species (Ce, Pr, and Nd) demonstrate anisotropic diffusion where the diffusion in the basal planes is significantly faster than that along the c-axis. Further, the calculated lanthanide diffusivities in HCP Zr are fitted to an Arrhenius relation and the activation energies and prefactors are reported for the first time. Furthermore, the vacancy drag and the segregation tendencies were analyzed using the calculated off-diagonal transport coefficients. According to our vacancy-mediated diffusion model, lanthanides are expected to be enriched at vacancy sinks at low temperatures, while at high temperatures, lanthanides are depleted at sinks and will preferably diffuse into the bulk. The enrichment/depletion transition temperature depends on the diffusion direction (basal or axial) and hence will be controlled by the grain texture and orientation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2004

Theoretical Investigation of OCN(-) Charge Transfer Complexes in Condensed Phase Media: Spectroscopic Properties in Amorphous Ice

Density functional theory (DFT) calculations of cyanate (OCN(-)) charge-transfer complexes were performed to model the "XCN" feature observed in interstellar icy grain mantles. OCN(-) charge-transfer complexes were formed from precursor combinations of HNCO or HOCN with either NH3 or H2O. Three different solvation strategies for realistically modeling the ice matrix environment were explored, including (1) continuum solvation, (2) pure DFT cluster calculations, and (3) an ONIOM DFT/PM3 cluster calculation. The model complexes were evaluated by their ability to reproduce seven spectroscopic measurements associated with XCN: the band origin of the OCN(-) asymmetric stretching mode, shifts in that frequency due to isotopic substitutions of C, N, O, and H, plus two weak features. The continuum solvent field method produced results consistent with some of the experimental data but failed to account for other behavior due to its limited capacity to describe molecular interactions with solvent. DFT cluster calculations successfully reproduced the available spectroscopic measurements very well. In particular, the deuterium shift showed excellent agreement in complexes where OCN(-) was fully solvated. Detailed studies of representative complexes including from two to twelve water molecules allowed the exploration of various possible solvation structures and provided insights into solvation trends. Moreover, complexes arising from cyanic or isocyanic acid in pure water suggested an alternative mechanism for the formation of OCN(-) charge-transfer complexes without the need for a strong base such as NH3 to be present. An extended ONIOM (B3LYP/PM3) cluster calculation was also performed to assess the impact of a more realistic environment on HNCO dissociation in pure water.

Park, Jin-Young↗
DOE OSTI2018

Molecular complexes and solvation interactions in the reaction of quinone imines with thiols

This study aims to investigate the role of complexation between reagents and the role of solvation of reagents by solvents in the kinetics of chain reactions of quinone imines with thiols. The thermodynamic characteristics of the complexation of quinone imines with thiophenol in CCl{sub 4}, chlorobenzene, and ethanol, as well as of the complexation of quinone imines and thiophenol with these solvents were calculated by quantum chemical methods (DFT calculations at the PBE/cc-pVDZ level of theory) and in terms of the additive-multiplicative model. Both approaches give consistent results. The formation of molecular complexes in quinone imine–thiphenol systems is accompanied by a 10–30 kJ mol{sup –1} decrease in enthalpy and has only a slight effect on the reaction mechanism.

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

Effect of hydrogen on plasticity of $α$-Fe: A multi-scale assessment

A multi-scale study was carried out to quantify the effect of interstitial hydrogen concentration on plasticity in $α$-Fe. Here, in this work, the influence of hydrogen on the screw dislocation glide behavior was examined across several length-scales. The insights obtained were integrated to provide an accurate continuum description for the effect of hydrogen on the dislocation based plasticity in polycrystalline $α$-Fe. At the outset of this work, a new Fe—H interatomic potential was formulated that enhanced the atomistic estimation of the variation in dislocation glide behavior in presence of hydrogen. Next, the dislocation core reconstruction observed due to the addition of hydrogen using atomistic simulations was validated with the help of large-scale DFT calculations based on the DFT-FE framework. Several atomistic simulations were carried out to comprehensively quantify the effect of hydrogen on the non-Schmid behavior exhibited during the dislocation glide in $α$-Fe. Finally, crystal plasticity simulations were carried out to understand the effect of hydrogen on the meso-scale deformation behavior of polycrystalline $α$-Fe.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Anisotropic Hydrolysis Susceptibility in Deformed Polydimethylsiloxanes

Chemical reactions involving the polydimethylsiloxane (PDMS) backbone can induce significant network rearrangements and ultimately degrade macro-scale mechanical properties of silicone components. Using two levels of quantum chemical theory, we identify a possible electronic driver for chemical susceptibility in strained PDMS chains and explore the complicated interplay between hydrolytic chain scissioning reactions, mechanical deformations of the backbone, water attack vector, and chain mobility. Density functional theory (DFT) calculations reveal that susceptibility to hydrolysis varies significantly with the vector for water attacks on silicon backbone atoms, which matches strain-induced anisotropic changes in the backbone electronic structure. Efficient semiempirical density functional tight binding (DFTB) calculations are shown to reproduce DFT predictions for select reaction pathways and facilitate more exhaustive explorations of configuration space. We show that concerted strains of the backbone must occur over at least a few monomer units to significantly increase hydrolysis susceptibility. In addition, we observe that sustaining tension across multiple monomer lengths by constraining molecular degrees of freedom further enhances hydrolysis susceptibility, leading to barrierless scission reactions for less substantial backbone deformations than otherwise. We then compute chain scission probabilities as functions of the backbone degrees of freedom, revealing complicated configurational interdependencies that impact the likelihood for hydrolytic degradation. As a result, the trends identified in our study suggest simple physical descriptions for the synergistic coupling between local mechanical deformation and environmental moisture in hydrolytic degradation of silicones.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Robust Methodology to Elucidate Kinetics of Room Temperature Electrochemical Propane Adsorption on Platinum

Electrocatalytic activation of alkanes can further decarbonize chemical manufacturing by leveraging affordable renewable electricity and readily available shale gas reserves in the United States. Earlier works have identified the unique role of Pt in adsorbing and activating alkanes, like propane, at room temperature in acidic, aqueous electrolytes, revealing spontaneous formation of deeply dehydrogenated propane-derived surface species with an intact C 3 - backbone. Although an adsorption mechanism was hypothesized, it has not been explicitly investigated to date, preventing the quantification of kinetic rate parameters. A robust methodology to investigate and benchmark propane adsorption kinetics on Pt is critical for the rational design of electrocatalysts that exhibit higher selectivity toward desired partially oxidized products. Herein, we analyze an oxidative current transience that appears during the adsorption of propane on Pt in aqueous electrochemical conditions and develop a methodology that elucidates the adsorption mechanism and enables quantification of rate parameters such as order dependences and apparent activation barriers. This method yields an expected first-order dependence with respect to propane concentration at low coverage and reveals a second-order dependence with respect to the concentration of surface active sites. Additionally, the apparent activation barrier for propane adsorption was calculated using an Arrhenius analysis of the current transience under temperature control. The experimentally measured activation barrier of 35 kJ mol –1 is in excellent agreement with the theoretical barrier calculated by density functional theory (DFT). The kinetic analysis was extended, via the use of transition state theory, to extract entropy and enthalpy of activation, yielding consistent results with the proposed two-step adsorption mechanism and DFT calculations. These results demonstrate reliable quantification of kinetic parameters for electrocatalytic activation of C–H bonds in alkanes that can be employed for rational catalyst development for a versatile range of electrocatalytic conditions.

alkane activation↗
DOE OSTI2021

Computational catalyst discovery: Active classification through myopic multiscale sampling

We report the recent boom in computational chemistry has enabled several projects aimed at discovering useful materials or catalysts. We acknowledge and address two recurring issues in the field of computational catalyst discovery. First, calculating macro-scale catalyst properties is not straightforward when using ensembles of atomic-scale calculations [e.g., density functional theory (DFT)]. We attempt to address this issue by creating a multi-scale model that estimates bulk catalyst activity using adsorption energy predictions from both DFT and machine learning models. The second issue is that many catalyst discovery efforts seek to optimize catalyst properties, but optimization is an inherently exploitative objective that is in tension with the explorative nature of early-stage discovery projects. In other words, why invest so much time finding a “best” catalyst when it is likely to fail for some other, unforeseen problem? We address this issue by relaxing the catalyst discovery goal into a classification problem: “What is the set of catalysts that is worth testing experimentally?” Here, we present a catalyst discovery method called myopic multiscale sampling, which combines multiscale modeling with automated selection of DFT calculations. It is an active classification strategy that seeks to classify catalysts as “worth investigating” or “not worth investigating” experimentally. Our results show an ~7–16 times speedup in catalyst classification relative to random sampling. These results were based on offline simulations of our algorithm on two different datasets: a larger, synthesized dataset and a smaller, real dataset.

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

Phonon Spectrum of SrFe2As2 determined by multizone phonon refinement

The ferropnictidesuperconductors exhibit a sensitive interplay between the lattice and magnetic degrees of freedom, including a number of phonon modes that are much softer than predicted by nonmagnetic calculations using density functional theory (DFT). However, it is not known what effect, if any, the long-range magnetic order has on phonon frequencies above 23 meV, where several phonon branches are very closely spaced in energy and it is challenging to isolate them from each other. We measured these phonons using inelastic time-of-flight neutron scattering in 40 Brillouin zones, and developed a technique to determine their frequencies. We find this method capable of determining phonon energies to 0.1 meV accuracy, and that the DFT calculations using the experimental structure yield qualitatively correct energies and eigenvectors. We do not find any effect of the magnetic transition on these phonons.

Parshall, D↗
DOE OSTI2023

A Progress Report for a 44-month Aging Study on Nitroplasticizer (NP) Stability - Summary of DFT and FTIR Results

In the past few years, a 44-month long aging experiment was conducted under various aging conditions. The properties of aged nitroplasticizer (NP) samples were analyzed using Fourier transform infrared spectroscopy, Karl Fischer (KF) titration, thermogravimetric analysis (TGA), liquid chromatography tandem quadrupole time of flight mass spectrometry (LC/QTOF), and ion chromatography (IC). In this progress report, some FTIR results are documented and discussed. On the theoretical front, density functional theory (DFT) calculations were carried out. Some DFT results are discussed as well.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Breakdown products of gaseous polycyclic aromatic hydrocarbons investigated with infrared ion spectroscopy

We report on a common fragment ion formed during the electron-ionization-induced fragmentation of three different three-ring polycyclic aromatic hydrocarbons (PAHs), fluorene (C{sub 13}H{sub 10}), 9,10-dihydrophenanthrene (C{sub 14}H{sub 12}), and 9,10-dihydroanthracene (C{sub 14}H{sub 12}). The infrared spectra of the mass-isolated product ions with m/z = 165 were obtained in a Fourier transform ion cyclotron resonance mass spectrometer whose cell was placed inside the optical cavity of an infrared free-electron laser, thus providing the high photon fluence required for efficient infrared multiple-photon dissociation. The infrared spectra of the m/z = 165 species generated from the three different precursors were found to be similar, suggesting the formation of a single C{sub 13}H{sub 9}{sup +} isomer. Theoretical calculations using density functional theory (DFT) revealed the fragment's identity as the closed-shell fluorenyl cation. Decomposition pathways from each parent precursor to the fluorenyl ion are proposed on the basis of DFT calculations. The identification of a single fragmentation product from three different PAHs supports the notion of the existence of common decomposition pathways of PAHs in general and can aid in understanding the fragmentation chemistry of astronomical PAH species.

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

NbTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys niobium-titanium (Nb-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP), https://github.com/ORNL/ScalableWorkflow_VASP_Calculations. Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Nb and Ti. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗