Search NASA⌕ Search

SEARCH · Search NASA

Results for “optimized energy density functionals”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The Triple Catalytic Action of Tertiary Nitrogen Catalysts in Recyclable Epoxy-Anhydride Thermosets

The thermosetting polymer matrix in fiber reinforced composites is an important component for energy related applications, such as the lightweighting of vehicles or their use in wind and waterpower turbine blades, due to their ability to provide superior adhesion, stiffness, and applicability to a wide range of manufacturing processes. Despite these benefits, today's thermosets are widely considered to be unrecyclable; thus, there is a large interest in redesigning these materials to be inherently recyclable so that energy intensive production of fibers and monomers can be circumvented, bolstering composite manufacture supply chains. Polyester covalent adaptable networks (PECANs) are one such promising alternative to the incumbent, nonrecyclable epoxy-amine thermosets. PECANs can be formed from the ring-opening co-polymerization (ROCOP) of epoxy-anhydride monomer mixtures and subsequent curing at mild temperatures to exhibit similar performance to conventional epoxies while also possessing unique dynamic chemistries along the ester-hydroxyl backbone that are capable of transesterification and thus reprocessability. While significant advancements have been made in formulating these materials for improved mechanical properties or optimizing solvolysis and reprocessing strategies, less attention has been placed on the impact of the residing amine catalyst used to generate the polyester network. In this work, we evaluated the triple-catalytic efficacy of 12 tertiary amines that act as a curing (bulk ROCOP), a transesterification (internal bond exchange), and a deconstruction (methanolysis) catalyst for PECAN thermosets. Specifically, we first distinguish between chain-growth and step-growth polymerization mechanisms for epoxy-amine and epoxy-anhydride mechanisms. We also utilized density functional theory (DFT) to estimate the basicity (pKb) of each catalyst. Of the tested catalysts, the ROCOP of the studied PECAN network can be completed between 95 and 247 min (at 80 degrees C), with variable gelation phenomena. Additionally, the stress relaxation (transesterification metric) efficiency of the tested PECAN networks with alternative embedded catalysts ranged from 95% to 15% reduction in stress after 5 h at 200 degrees C, and the depolymerization efficacy ranged from 2.5% to 9.8% deconstruction after 36 h at 130 degrees C. Overall, the nitrogen-based moieties were demonstrated to influence polymerization kinetics, catalyze the dynamic transesterification exchange mechanism, and aid in the solvolysis of the thermosets at end-of-life.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Hollow-structured Ni-N-C catalysts for highly selective CO 2 electroreduction

Atomically dispersed single-atom catalysts have emerged as promising non-precious catalyst alternatives to expensive Ag and Au catalysts for electrochemical CO 2 reduction reaction (CO 2 RR). In particular, nickel-nitrogen-carbon (Ni-N-C) catalysts have demonstrated a high faradaic efficiency (FE) toward CO formation at low overpotentials. Nonetheless, the exact nature of Ni active sites under CO 2 RR remains elusive and conventional Ni-N-C catalysts are limited by microporosity and low density of Ni single atoms, hindering performance in CO 2 electrolyzers. Here, we report the synthesis of hollow-structured Ni-N-C ( hs -Ni-N-C) catalysts via a post-synthesis modification (PSM) strategy using partial ligand exchange of 2-methylimidazole with 3-amino-1,2,4-triazole. This approach enables the formation of a hollow structure, resulting in more than a twofold increase in Ni atom density compared to regular Ni-N-C (r-Ni-N-C). In a zero-gap CO 2 electrolyzer, the optimized hs -Ni-N-C allows for achieving an FE CO of 97% at a current density of > 100 mA cm⁻ 2 , while maintaining high CO selectivity with stable performance over 100 h at 2.5 V. hs-Ni-N-C shows a more than sevenfold increase in the CO partial current density relative to r-Ni-N-C resulting from the combined effects of a higher density of Ni single-atom sites, improved kinetics, and lower transport resistance under the operating conditions, as indicated by electrochemical impedance spectra and distribution of relaxation times analysis. Operando high energy-resolution X-ray absorption spectroscopy (XAS) reveals that atop-bonded CO on Ni single sites induces dynamic transformations of the Ni–N coordination environment, leading to a symmetric coordination structure of hs -Ni-N-C. Under CO 2 RR, the catalysts undergo a more pronounced structural change and form a minor fraction of Ni nanoparticles. Density functional theory calculations are consistent with the XAS results and provide molecular insights showing that the interplay between protonation and CO adsorption leads to adsorbate-induced restructuring of the Ni single atom. This work demonstrates the synergistic role of hollow structure and high-density Ni atoms in governing CO 2 RR selectivity and provides mechanistic insights into the structural dynamics of single-atom catalysts under operating conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Enhancing the accuracy and generality of the Debye–Grüneisen Model: Optimizing the volume dependence for accurate predictions across varied compositions

In this work, we have introduced an optimized Debye-Grüneisen model that revolutionizes the determination of the Debye temperature and Grüneisen parameters. Unlike conventional methods, our model requires only the 0 K energy volume data for a material as input, eliminating the need to determine the bulk modulus and its pressure derivative, which often pose challenges due to numerical uncertainties. This unique feature sets our model apart from existing approaches and streamlines the process, enabling accurate predictions of thermal expansion behavior across various materials. To demonstrate its effectiveness, we showcase its excellent agreement with measured coefficients of thermal expansion (CTE) for the nickel-cobalt-chromium-aluminum-yttrium (Ni-Co-Cr-Al-Y) bond-coating system. Additionally, we apply our approach by conducting a high-throughput search for potential bond-coating materials among 90,000 compositions within the aluminum-cobalt-chromium-iron-nickel (Al-Co-Cr-Fe-Ni) system. From this extensive search, four compositions are synthesized, and the measured CTE values agree very well with theoretical predictions, hence validating our approach. In conclusion, the current optimized Debye-Grüneisen model combined with Density Functional Theory (DFT)-based thermodynamic database enables reliable and efficient high-throughput calculations of CTE of of a material without expensive phonon calculations.

Bond coating materials↗
DOE OSTI2022

Tunable metal hydroxide–organic frameworks for catalysing oxygen evolution

The oxygen evolution reaction is central to making chemicals and energy carriers using electrons. Combining the great tunability of enzymatic systems with known oxide-based catalysts can create breakthrough opportunities to achieve both high activity and stability. Here we report a series of metal hydroxide–organic frameworks (MHOFs) synthesized by transforming layered hydroxides into two-dimensional sheets crosslinked using aromatic carboxylate linkers. MHOFs act as a tunable catalytic platform for the oxygen evolution reaction, where the π–π interactions between adjacent stacked linkers dictate stability, while the nature of transition metals in the hydroxides modulates catalytic activity. Substituting Ni-based MHOFs with acidic cations or electron-withdrawing linkers enhances oxygen evolution reaction activity by over three orders of magnitude per metal site, with Fe substitution achieving a mass activity of 80A $g^{-1}_{catalyst}$ 1 at 0.3 V overpotential for 20 h. Density functional theory calculations correlate the enhanced oxygen evolution reaction activity with the MHOF-based modulation of Ni redox and the optimized binding of oxygenated intermediates.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Revealing the Nature of Binary-Phase on Structural Stability of Sodium Layered Oxide Cathodes

The emergence of layered sodium transition metal oxides featuring a multiphase structure presents a promising approach for cathode materials in sodium-ion batteries, showcasing notably improved energy storage capacity. However, the advancement of cathodes with multiphase structures faces obstacles due to the limited understanding of the integrated structural effects. Herein, the integrated structural effects by an in-depth structure-chemistry analysis in the developed layered cathode system Na x Cu 0.1 Co 0.1 Ni 0.25 Mn 0.4 Ti 0.15 O 2 with purposely designed P2/O3 phase integration, are comprehended. The results affirm that integrated phase ratio plays a pivotal role in electrochemical/structural stability, particularly at high voltage and with the incorporation of anionic redox. In contrast to previous reports advocating solely for the enhanced electrochemical performance in biphasic structures, it is demonstrated that an inappropriate composite structure is more destructive than a single-phase design. The in situ X-ray diffraction results, coupled with density functional theory computations further confirm that the biphasic structure with P2:O3 = 4:6 shows suppressed irreversible phase transition at high desodiated states and thus exhibits optimized electrochemical performance. Finally, these fundamental discoveries provide clues to the design of high-performance layered oxide cathodes for next-generation SIBs.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Computational Study of the Resonance Enhancement of Raman Signals of Ligands Adsorbed to CdSe Clusters through Photoexcitation of the Cluster

This paper describes density-functional-theory-based computations of resonance Raman (RR) spectra of ligand molecules adsorbed to the surface of a Cd 16 Se 13 cluster. Signals from asymmetric vibrational modes of ligand binding groups, such as the asymmetric O–C–O stretching modes of carboxylates, are enhanced relative to the symmetric vibrational modes when the excitation energy is on-resonance with the excitonic energy of the cluster. Certain ligand molecules have frontier orbitals with the correct energies and symmetries to mix with the orbitals of the CdSe cluster, and as a result, the wave functions of the electron and the hole delocalize from the cluster onto the ligand molecules; experimentally, this delocalization results in a bathochromic shift of the band edge excitonic absorption. Increased excitonic delocalization results in greater vibronic coupling between the exciton and the ligand vibrations and, on average, preferential enhancements in the RR signals of those vibrations. This work suggests that the use of exciton-delocalizing ligands to optimize electronic coupling between neighboring CdSe nanoparticles may, at the same time, enhance the rates of nonradiative exciton decay by coupling the exciton and ligand vibrational modes.

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

Substituent and Heteroatom Effects on π–π Interactions: Evidence That Parallel-Displaced π-Stacking is Not Driven by Quadrupolar Electrostatics

Stacking interactions are a recurring motif in supramolecular chemistry and biochemistry, where a persistent theme is a preference for parallel-displaced aromatic rings rather than face-to-face π-stacking. This is usually explained in terms of quadrupole–quadrupole interactions between the arene moieties but that interpretation is inconsistent with accurate calculations, which reveal that the quadrupolar picture is qualitatively wrong. At typical π-stacking distances, quadrupolar electrostatics may differ in sign from an exact calculation based on charge densities of the interacting arenes. We apply symmetry-adapted perturbation theory to dimers composed of substituted benzene and various aromatic heterocycles, which display a wide range of electrostatic interactions, and we investigate the interplay of Pauli repulsion, dispersion, and electrostatics as it pertains to parallel-displaced π-stacking. Profiles of energy components along cofacial slip-stacking coordinates support a prominent role for the “van der Waals model” (dispersion in competition with Pauli repulsion), even for polar monomers where electrostatic interactions are significant. While electrostatic interactions are necessary to explain the optimal face-to-face π-stacking distance and to account for the relative orientation of one polar arene with respect to another, we find no evidence to support continued invocation of quadrupolar electrostatics as a basis for π-stacking. Our results suggest that a driving force for offset-stacking exists even in the absence of electrostatic interactions. Consequently, tuning electrostatics via functionalization does not guarantee that slip-stacking can be avoided. This has implications for rational design of soft materials and other supramolecular architectures.

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

Mixed Ionic Liquid Improves Electrolyte Dynamics in Supercapacitors

Well-tailored mixtures of distinct ionic liquids can act as optimal electrolytes that extend the operating electrochemical window and improve charge storage density in supercapacitors. Here, we explore two room-temperature ionic liquids, 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (EmimTFSI) and 1-ethyl-3-methylimidazolium tetrafluoroborate (EmimBF 4 ). We study their electric double-layer behavior in the neat state and as binary mixtures on the external surfaces of onion-like carbon electrodes using quasielastic neutron scattering (QENS) and classical density functional theory techniques. Computational results reveal that a mixture with 4:1 EmimTFSI/EmimBF 4 volume ratio displaces the larger [TFSI – ] anions with smaller [BF 4 – ] ions, leading to an excess adsorption of [Emim + ] cations near the electrode surface. These findings are corroborated by the manifestation of nonuniform ion diffusivity change, complementing the description of structural modifications with changing composition, from QENS measurements. In conclusion, molecular-level understanding of ion packing near electrodes provides insight for design of ionic liquid formulations that enhance the performance of electrochemical energy storage devices.

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

Site Identity and Importance in Cosubstituted Bixbyite In2O3

The bixbyite structure of In2O3 has two nonequivalent, 6-coordinate cation sites and, when Sn is doped into In2O3, the Sn prefers the “b-site” and produces a highly conductive material. When divalent/tetravalent cation pairs are cosubstituted into In2O3, however, the conductivity increases to a lesser extent and the site occupancy is less understood. We examine the site occupancy in the MgxIn2−2xSnxO3 and ZnxIn2−2xSnxO3 systems with high resolution X-ray and neutron diffraction and density functional theory computations, respectively. In these sample cases and those that are previously reported in the MxIn2−2xSnxO3 (M = Cu, Ni, or Zn) systems, the solubility limit is greater than 25%, ensuring that the b-site cannot be the exclusively preferred site as it is in Sn:In2O3. Prior to this saturation point, we report that the M2+ cation always has at least a partial occupancy on the d-site and the Sn4+ cation has at least a partial occupancy on the b-site. The energies of formation for these configurations are highly favored, and prefer that the divalent and tetravalent substitutes are adjacent in the crystal lattice, which suggests short range ordering. Diffuse reflectance and 4-point probe measurements of MgxIn2−xSnxO3 demonstrate that it can maintain an optical band gap >2.8 eV while surpassing 1000 S/cm in conductivity. Understanding how multiple constituents occupy the two nonequivalent cation sites can provide information on how to optimize cosubstituted systems to increase Sn solubility while maintaining its dopant nature, achieving maximum conductivity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Water limiting current measurements in anion exchange membrane fuel cells (AEMFCs); part 1: Water limiting current method development

We report water management is critical in optimizing anion exchange membrane fuel cell (AEMFCs) performance and durability. The role of membranes on determining cell water balance in AEMFCs has received limited study. A more thorough understanding of membrane water transport parameters is necessary for thorough cell optimization. In this study, a novel method for probing limiting current based on water flux across an AEMFC is presented for the first time. The water limiting current has been investigated as a function of relative humidity of anode and cathode streams, cathode backpressure and cell operating temperature. From these studies, it is clearly shown that water diffusion through the membrane from the anode is the most critical mechanism for achieving high current density performance and durability. The developed method can act as a novel screening technique that focuses on the impacts of water transport parameters and helps identify promising membrane materials based on their water transport properties. A critical but underappreciated facet of AEMFC performance and durability.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2013

Excitation of helium Rydberg states and doubly excited resonances in strong extreme ultraviolet fields: Full-dimensional quantum dynamics using exponentially tempered Gaussian basis sets

Recently optimized exponentially tempered Gaussian basis sets [P. R. Kapralova-Zdanska and J. Smydke, J. Chem. Phys. 138, 024105 (2013)] are employed in quantitative simulations of helium absorption cross-sections and two-photon excitation yields of doubly excited resonances. Linearly polarized half-infinite and Gaussian laser pulses at wavelengths 38–58 nm and large intensities up to 100 TW/cm{sup 2} are considered. The emphasis is laid on convergence of the results with respect to the quality of the Gaussian basis sets (typically limited by a number of partial waves, density, and spatial extent of the basis functions) as well as to the quality of the basis set of field-free states (typically limited by the maximum rotational quantum number and maximum excitation of the lower electron). Particular attention is paid to stability of the results with respect to varying complex scaling parameter. Moreover, the study of the dynamics is preceded by a thorough check of helium energies and oscillator strengths as they are obtained with the exponentially tempered Gaussian basis sets, being also compared with yet unpublished emission wavelengths measured in electric discharge experiments.

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

Sensitizing properties of luminescence centers on the emission of Er{sup 3+} in Si-rich SiO{sub 2} film

In this paper, we report on the luminescence-center (LC)-mediated excitation of Er{sup 3+} as a function of annealing temperature in Er-doped Si-rich SiO{sub 2} (SRO) films fabricated by electron beam evaporation. It is found that the annealing temperature has significant effects on the emission of Er{sup 3+} and the specific optical-active point-defects called LCs within Er-doped SRO films. Different luminescence centers generated by the evolution of microstructures during annealing process act as efficient sensitizers for Er{sup 3+} in the films when the annealing temperature is below 1100 °C. Moreover, the temperature dependence of the energy coupling between LCs and Er{sup 3+} demonstrates the effective phonon-mediated energy transfer process. In addition, when the annealing temperature reaches 1100 °C, the decreased density of activable erbium ions induced by the aggregation of Er will bring detrimental effects on the emission of Er{sup 3+}. It is demonstrated that an appropriate annealing process can be designed to achieve efficiently enhanced emissions from Er{sup 3+} ions by optimizing the density of LCs and the coupling between Er{sup 3+} and LCs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Benchmarking and validation of global model code for negative hydrogen ion sources

Benchmarking and validation are prerequisites for using simulation codes as predictive tools. In this work, we have developed a Global Model for Negative Hydrogen Ion Source (GMNHIS) and performed benchmarking of the GMNHIS against another independently developed code, Global Enhanced Vibrational Kinetic Model (GEVKM). This is the first study to propose a quite comprehensive benchmarking test of this kind for models of negative hydrogen ion sources (NHIS), and excellent agreements have been achieved for collisional energy loss per electron-ion pair created, electron number density, electron temperature, densities of H$^{+}_{3}$ and H$^{+}_{2}$ ions, and densities of H(n = 1–3) atoms. Very small discrepancies in number densities of H - ions and H + ions, as well as the vibrational distribution function of hydrogen molecules, can be attributed to the differences in the chemical reactions datasets. The GEVKM includes additional chemical reactions that are more important at high pressures. Moreover, we validated the GMNHIS against experimental data obtained in an electron cyclotron resonance discharge used for H - production. The model qualitatively (and even quantitatively for certain conditions) reproduces the experimental H - number density. The H - number density as a function of pressure first increases at pressures below 1.6 Pa and then saturates for higher pressures. This dependence was analyzed by evaluating contributions from different reaction pathways to the creation and loss of the H - ions. The developed codes can be used for predicting the H - production, improving the performance of NHIS, and ultimately optimizing the parameters of negative ion beams for fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Fluid-Kinetic Coupling: Advanced Discretizations for Simulations on Emerging Heterogeneous Architectures (LDRD FY20-0643)

Plasma physics simulations are vital for a host of Sandia mission concerns, for fundamental science, and for clean energy in the form of fusion power. Sandia's most mature plasma physics simulation capabilities come in the form of particle-in-cell (PIC) models and magnetohydrodynamics (MHD) models. MHD models for a plasma work well in denser plasma regimes when there is enough material that the plasma approximates a fluid. PIC models, on the other hand, work well in lower-density regimes, in which there is not too much to simulate; error in PIC scales as the square root of the number of particles, making high-accuracy simulations expensive. Real-world applications, however, almost always involve a transition region between the high-density regimes where MHD is appropriate, and the low-density regimes for PIC. In such a transition region, a direct discretization of Vlasov is appropriate. Such discretizations come with their own computational costs, however; the phase-space mesh for Vlasov can involve up to six dimensions (seven if time is included), and to apply appropriate homogeneous boundary conditions in velocity space requires meshing a substantial padding region to ensure that the distribution remains sufficiently close to zero at the velocity boundaries. Moreover, for collisional plasmas, the right-hand side of the Vlasov equation is a collision operator, which is non-local in velocity space, and which may dominate the cost of the Vlasov solver. The present LDRD project endeavors to develop modern, foundational tools for the development of continuum-kinetic Vlasov solvers, using the discontinuous Petrov-Galerkin (DPG) methodology, for discretization of Vlasov, and machine-learning (ML) models to enable efficient evaluation of collision operators. DPG affords several key advantages. First, it has a built-in, robust error indicator, allowing us to adapt the mesh in a very natural way, enabling a coarse velocity-space mesh near the homogeneous boundaries, and a fine mesh where the solution has fine features. Second, it is an inherently high-order, high-intensity method, requiring extra local computations to determine so-called optimal test functions, which makes it particularly suited to modern hardware in which floating-point throughput is increasing at a faster rate than memory bandwidth. Finally, DPG is a residual-minimizing method, which enables high-accuracy computation: in typical cases, the method delivers something very close to the $L^2$ projection of the exact solution. Meanwhile, the ML-based collision model we adopt affords a cost structure that scales as the square root of a standard direct evaluation. Moreover, we design our model to conserve mass, momentum, and energy by construction, and our approach to training is highly flexible, in that it can incorporate not only synthetic data from direct-simulation Monte Carlo (DSMC) codes, but also experimental data. We have developed two DPG formulations for Vlasov-Poisson: a time-marching, backward-Euler discretization and a space-time discretization. We have conducted a number of numerical experiments to verify the approach in a 1D1V setting. In this report, we detail these formulations and experiments. We also summarize some new theoretical results developed as part of this project (published as papers previously): some new analysis of DPG for the convection-reaction problem (of which the Vlasov equation is an instance), a new exponential integrator for DPG, and some numerical exploration of various DPG-based time-marching approaches to the heat equation. As part of this work, we have contributed extensively to the Camellia open-source library; we also describe the new capabilities and their usage. We have also developed a well-documented methodology for single-species collision operators, which we applied to argon and demonstrated with numerical experiments. We summarize those results here, as well as describing at a high level a design extending the methodology to multi-species operators. We have released a new open-source library, MLC, under a BSD license; we include a summary of its capabilities as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Electrochemical dynamics of imidazolium ionic liquids at graphene electrodes for energy storage applications

Electric double layer capacitors (EDLCs) are prominent energy storage systems that constitute the foundation of more reliable and sustainable energy infrastructures. Modern EDLCs often incorporate ionic liquids (ILs) as a key component in their electrolytes, leveraging the high electrochemical stability of ILs to enhance device performance. The performance and functionality of these capacitors also hinge on the interfacial behavior of ILs at electrode surfaces, which remain insufficiently understood. Here, we performed synchrotron infrared nanospectroscopy (SINS) in combination with density functional theory (DFT) calculations to investigate the electric double layers (EDLs) of three imidazolium-based ILs in a custom-designed graphene liquid cell. This approach revealed new insights into the dynamics of IL EDLs and the underpinning factors originating from the IL structures. Variations in anion size and structure were found to tune the ILs’ ability to form EDLs with compact and closely-correlated ion arrangements, which are critical for enhancing their capacitive performance. These findings highlight the intricate interactions between ions governed by their structures and charge behaviors, which underscores the opportunities for targeted design of IL-based electrolytes to optimize EDLC functionality.

Electric double layers↗
DOE OSTI2023

NbTaV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution ternary alloys niobium-tantalum-vanadium (Nb-Ta-V). The first-principle code that has been used to run the calculations is the closed-source Vienna Ab-Initio Simulation Package (VASP). 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 8. 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 10,500 randomized atomic structures over 105 chemical compositions. The calculations have been collected on NERSC-Perlmutter using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for ternary alloys is AXBYCZ, where A, B, and C refer to the constituents, and X, Y, and Z are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. 0.POSCAR: input file that defines the atomic structure of a system 3. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 4. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 6. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 7. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 8. OUTCARL: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 10. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 11. vasprun.xml: contains similar information to OUTCAR, but in an xml format. This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216.

36 MATERIALS SCIENCE↗
DOE OSTI2016

SU-C-204-03: DFT Calculations of the Stability of DOTA-Based-Radiopharmaceuticals

Purpose: Application of the density function theory (DFT) to investigate the structural stability of complexes applied in cancer therapy consisting of the 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid (DOTA) chelated to Ac225, Fr221, At217, Bi213, and Gd68 radio-nuclei. Methods: The possibility to deliver a toxic payload directly to tumor cells is a highly desirable aim in targeted alpha particle therapy. The estimation of bond stability between radioactive atoms and the DOTA chelating agent is the key element in understanding the foundations of this delivery process. Thus, we adapted the Vienna Ab-initio Simulation Package (VASP) with the projector-augmented wave method and a plane-wave basis set in order to study the stability and electronic properties of DOTA ligand chelated to radioactive isotopes. In order to count for the relativistic effect of radioactive isotopes we included Spin-Orbit Coupling (SOC) in the DFT calculations. Five DOTA complex structures were represented as unit cells, each containing 58 atoms. The energy optimization was performed for all structures prior to calculations of electronic properties. Binding energies, electron localization functions as well as bond lengths between atoms were estimated. Results: Calculated binding energies for DOTA-radioactive atom systems were −17.792, −5.784, −8.872, −13.305, −18.467 eV for Ac, Fr, At, Bi and Gd complexes respectively. The displacements of isotopes in DOTA cages were estimated from the variations in bond lengths, which were within 2.32–3.75 angstroms. The detailed representation of chemical bonding in all complexes was obtained with the Electron Localization Function (ELF). Conclusion: DOTA-Gd, DOTA-Ac and DOTA-Bi were the most stable structures in the group. Inclusion of SOC had a significant role in the improvement of DFT calculation accuracy for heavy radioactive atoms. Our approach is found to be proper for the investigation of structures with DOTA-based-radiopharmaceuticals and will enhance our understanding of processes occurring at subatomic levels.

60 APPLIED LIFE SCIENCES↗