Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional theory calculations”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Collaborative Simulation Grid: Multiscale Quantum-Mechanical/Classical Atomistic Simulations on Distributed PC Clusters in the US and Japan

A multidisciplinary, collaborative simulation has been performed on a Grid of geographically distributed PC clusters. The multiscale simulation approach seamlessly combines i) atomistic simulation backed on the molecular dynamics (MD) method and ii) quantum mechanical (QM) calculation based on the density functional theory (DFT), so that accurate but less scalable computations are performed only where they are needed. The multiscale MD/QM simulation code has been Grid-enabled using i) a modular, additive hybridization scheme, ii) multiple QM clustering, and iii) computation/communication overlapping. The Gridified MD/QM simulation code has been used to study environmental effects of water molecules on fracture in silicon. A preliminary run of the code has achieved a parallel efficiency of 94% on 25 PCs distributed over 3 PC clusters in the US and Japan, and a larger test involving 154 processors on 5 distributed PC clusters is in progress.

Kikuchi, Hideaki↗
NASA NTRS2003

Scalability of a Low-Cost Multi-Teraflop Linux Cluster for High-End Classical Atomistic and Quantum Mechanical Simulations

Scalability of a low-cost, Intel Xeon-based, multi-Teraflop Linux cluster is tested for two high-end scientific applications: Classical atomistic simulation based on the molecular dynamics method and quantum mechanical calculation based on the density functional theory. These scalable parallel applications use space-time multiresolution algorithms and feature computational-space decomposition, wavelet-based adaptive load balancing, and spacefilling-curve-based data compression for scalable I/O. Comparative performance tests are performed on a 1,024-processor Linux cluster and a conventional higher-end parallel supercomputer, 1,184-processor IBM SP4. The results show that the performance of the Linux cluster is comparable to that of the SP4. We also study various effects, such as the sharing of memory and L2 cache among processors, on the performance.

Kikuchi, Hideaki↗
NASA NTRS2018

Application of Ab Initio Methods in the Development of Advanced Technical Ceramics

A significant need exists to develop materials not only capable of providing desired electronic and mechanical properties but also survival at extreme temperatures during service and device fabrication, such as co-firing. Computational methods offers an efficient and systematic manner to design new materials and guide their development. As an example computational -based material approaches can be used to determine the suitability of a given materials as a practical thermoelectric for energy harvesting. In this presentation examples are given of applications to relevant technical ceramics such as thermoelectrics, dielectrics, and magnetic systems. The calculation were based on density functional theory and carried out with norm conserving and projector augmented wave (PAW) methods using commercial codes Materials Studio ( Biovia, Inc) with the Cambridge Serial Total Energy Package (CASTEP) and MedeA (Materials Design Inc.) utilizing the Vienna Ab-initio Simulation Package (VASP) as the respective computational engines. This study makes predictions of relevant technical properties of ceramic materials.

thermoelectric↗
NASA NTRS1995

Theoretical Infrared Spectra for Polycyclic Aromatic Hydrocarbon Neutrals, Cations and Anions

Calculations are carried out using density functional theory (DFT) to determine the harmonic frequencies and intensities of the neutrals and cations of thirteen polycyclic aromatic hydrocarbons (PAHs) up to the size of ovalene. Calculations are also carried out for a few PAH anions. The DFT harmonic frequencies, when uniformly scaled by the factor of 0.958 to account primarily for anharmonicity, agree with the matrix isolation fundamentals to within an average error of about 10 per centimeter. Electron correlation is found to significantly reduce the intensities of many of the cation harmonics, bringing them into much better agreement with the available experimental data. While the theoretical infrared spectra agree well with the experimental data for the neutral systems and for many of the cations, there are notable discrepancies with the experimental matrix isolation data for some PAH cations that are difficult to explain in terms of limitations in the calculations. In agreement with previous theoretical work, the present calculations show that the relative intensities for the astronomical unidentified infrared (UIR) bands agree reasonably well with those for a distribution of polycyclic aromatic hydrocarbon (PAH) cations, but not with a distribution of PAH neutrals. We also observe that the infrared spectra of highly symmetrical cations such as coronene agree much better with astronomical observations than do those of, for example, the polyacenes such as tetracene and pentacene. The total integrated intensities for the neutral species are found to increase linearly with size, while the total integrated intensities are much larger for the cations and scale more nearly quadratically with size. We conclude that emission from moderate-sized highly symmetric PAH cations such as coronene and larger could account for the UIR bands.

Langhoff, Stephen R.↗
NASA NTRS1999

Effective Forces Between Colloidal Particles

Colloidal suspensions have proven to be excellent model systems for the study of condensed matter and its phase behavior. Many of the properties of colloidal suspensions can be investigated with a systematic variation of the characteristics of the systems and, in addition, the energy, length and time scales associated with them allow for experimental probing of otherwise inaccessible regimes. The latter property also makes colloidal systems vulnerable to external influences such as gravity. Experiments performed in micro-ravity by Chaikin and Russell have been invaluable in extracting the true behavior of the systems without an external field. Weitz and Pusey intend to use mixtures of colloidal particles with additives such as polymers to induce aggregation and form weak, tenuous, highly disordered fractal structures that would be stable in the absence of gravitational forces. When dispersed in a polarizable medium, colloidal particles can ionize, emitting counterions into the solution. The standard interaction potential in these charged colloidal suspensions was first obtained by Derjaguin, Landau, Verwey and Overbeek. The DLVO potential is obtained in the mean-field linearized Poisson-Boltzmann approximation and thus has limited applicability. For more precise calculations, we have used ab initio density functional theory. In our model, colloidal particles are charged hard spheres, the counterions are described by a continuum density field and the solvent is treated as a homogeneous medium with a specified dielectric constant. We calculate the effective forces between charged colloidal particles by integrating over the solvent and counterion degrees of freedom, taking into account the direct interactions between the particles as well as particle-counterion, counterion-counterion Coulomb, counterion entropic and correlation contributions. We obtain the effective interaction potential between charged colloidal particles in different configurations. We evaluate two- and three-body forces in the bulk as well as study the influence of soft walls. We qualitatively explain the effects of the walls on the forces and demonstrate that many-body effects are negligible in our system. With adjustments in the parameters, the DLVO pair-potential can describe the results quantitatively. Besides electrostatic interactions, entropic depletion effects that arise from (hard-core) exclusion play an important role in determining the behavior of multi-component colloidal suspensions. A standard theory for depletion forces is due to Asakura and Oosawa and is based on the ideal gas approximation. To go beyond this approximation, we have studied entropic forces in molecular dynamics simulations of systems of hard spheres (the effects of the solvent have been ignored). The effective depletion forces for these systems can be found either from equilibrium distribution functions or from direct momentum transfer calculations. Our results obtained by either method show qualitative differences from the Asakura-Oosawa forces, indicating a longer range, higher value at contact and most importantly a more complicated structure, comprising of several maxima and minima. Our calculations include the determination of effective forces between two spheres, a hard sphere and a wall, and the behavior of a hard sphere near a step-edge and a corner. We also demonstrate that such entropic forces do not necessarily satisfy pairwise additivity.

Tehver, Riina↗
NASA NTRS2001

Closed-Shell Polycyclic Aromatic Hydrocarbon Cations: A New Category of Interstellar Polycyclic Aromatic Hydrocarbons

Density functional theory has been employed to calculate the harmonic frequencies and intensities of a range of polycyclic aromatic hydrocarbon (PAH) cations that explore both size and electronic structure effects of the infrared spectroscopic properties of these species. The sample extends the size range of PAH species considered to more than 50 carbon atoms and includes several representatives from each of two heretofore unexplored categories of PAH cations: (1) fully benzenoid PAH cations whose carbon skeleton is composed of an odd number of carbon atoms (C(sub odd) PAHs); and (2) protonated PAH cations (HPAH+). Unlike the radical electronic structures of the PAH cations that have been the subject of previous theoretical and experimental work, the species in these two classes have a closed-shell electronic configuration. The calculated spectra of circumcoronene, C54H18 in both neutral and (radical) cationic form are also reported and compared with those of the other species. Overall, the C(sub odd) PAHs spectra are dominated by strong CC stretching modes near 1600 cm(exp -1) and display spectra that are remarkably insensitive to molecular size. The HPAH+ species evince a more complex spectrum consistent with the added contributions of aliphatic modes and their generally lower symmetry. Finally, for both classes of closed-shell cations, the intensity of the aromatic CH stretching modes is found to increase with molecular size far out of proportion with the number of CH groups, approaching a value more typical of neutral PAHs for the largest species studied.

Hudgins, Douglas M.↗
NASA NTRS2000

A Theoretical Investigation of the Infrared Spectroscopic Properties of Closed-Shell Polycyclic Aromatic Hydrocarbon Cations

Density functional theory has been employed to calculate the harmonic frequencies and intensities of a range of PAH cations which explore both size and electronic structure effects on the infrared spectroscopic of these species. The sample extends the size range of PAH species considered to more than 50 carbon atoms and includes several representatives from each of two heretofore unexplored categories of PAH cations: (1) fully benzenoid PAH cations whose carbon skeleton is composed of an odd number of carbon atoms and (2) protonated PAH cations. Unlike the radical electronic structures of the PAH cations that have been the subject of previous theoretical and experimental work, the species in these two classes have a closed-shell electronic configuration. The calculated spectra of circumcoronene, C54H18, in both neutral and (radical) cationic form are also reported and compared to those of the other species. Closed-shell species are inherently less reactive than radical (or open-shell) cations and are known to play a role in combustion chemistry. Since interstellar PAHs are typically exposed to abundant atomic hydrogen and are thought to originate under pseudo-combustion conditions in carbon-rich circumstellar shells, such species may represent an important component of the interstellar PAH population. Furthermore, species larger than 50 carbon atoms are more representative of the size of the typical interstellar PAH. Overall, as has been the case for previous studies of PAH radical cations, the general pattern of band positions and intensities are consistent with that of the interstellar infrared emission spectrum. In addition, the spectra of closed-shell and open-shell cations are found to converge with increasing molecular size and are found to be relatively similar for species containing about 50 carbon atoms.

Hudgins, Douglas M.↗
NASA NTRS1994

The Molecular Structure of cis-FONO

The molecular structure of cis-FONO has been determined with the CCSD(T) correlation method using an spdf quality basis set. In agreement with previous coupled-cluster calculations but in disagreement with density functional theory, cis-FONO is found to exhibit normal bond distances. The quadratic and cubic force fields of cis-FONO have also been determined in order to evaluate the effect of vibrational averaging on the molecular geometry. Vibrational averaging is found to increase bond distances, as expected, but it does not affect the qualitative nature of the bonding. The CCSD(T)/spdf harmonic frequencies of cis-FONO support our previous assertion that a band observed at 1200 /cm is a combination band (upsilon(sub 3) + upsilon(sub 4)), and not a fundamental.

Lee, Timothy J.↗
NASA NTRS1996

Infrared Spectra of Perdeuterated Naphthalene, Phenanthrene, Chrysene, and Pyrene

Calculations are carried out using density functional theory (DFT) to determine the harmonic frequencies and intensities of perdeuterated naphthalene, phenanthrene, pyrene, and chrysene. We also report matrix- isolation spectra for these four species. The theoretical and experimental frequencies and relative intensities for the perdeuterated species are in generally good agreement. The effect of perdeuteration is to reduce the sum of the integrated intensities by a factor of about 1.75. This reduction occurs for all vibrational motions, except for the weak low frequency ring deformation modes. There is also a significant redistribution of the relative intensities between the out-of-plane C-D bands relative to those found for the out-of-plane C-H bands. The theoretical isotopic ratios provide an excellent diagnostic of the degree of C-H(C-D) involvement in the vibrational bands, allowing in most cases a clear distinction of the type of motion.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1997

Infrared Spectra of Substituted Polycyclic Aromatic Hydrocarbons

Calculations are carried out using density functional theory (DFT) to determine the harmonic frequencies and intensities of 1-methylanthracene, 9-methylanthracene, 9-cyanoanthracene, 2-aminoanthracene, acridine, and their positive ions. The theoretical data are compared with matrix-isolation spectra for these species also reported in this work. The theoretical and experimental frequencies and relative intensities for the neutral species are in generally good agreement, whereas the positive ion spectra are only in qualitative agreement. Relative to anthracene, we find that substitution of amethyl or CN for a hydrogen does not significantly affect the spectrum other than to add the characteristic methyl C-H stretch and C-N stretch (near 2200/cm), respectively. However, addition of NH2 dramatically affects the spectrum of the neutral. Not only are the NH2 modes themselves strong, but this electron withdrawing group induces sufficient partial charge on the ring to give the neutral molecule characteristics of the anthracene cation spectrum. The sum of the absolute intensities is about four times larger for 2-aminoanthracene than for 9-cyanoanthracene. Substituting nitrogen in the ring at the nine position (acridine) does not greatly alter the spectrum compared with anthracene.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1998

Infrared spectra of substituted polycylic aromatic hydrocarbons

Calculations are carried out using density functional theory (DFT) to determine the harmonic frequencies and intensities of 1-methylanthracene, 9-methylanthracene, 9-cyanoanthracene, 2-aminoanthracene, acridine, and their positive ions. The theoretical data are compared with matrix-isolation spectra for these species also reported in this work. The theoretical and experimental frequencies and relative intensities for the neutral species are in generally good agreement, whereas the positive ion spectra are only in qualitative agreement. Relative to anthracene, we find that substitution of a methyl or CN for a hydrogen does not significantly affect the spectrum other than to add the characteristic methyl C-H and C triple bond N stretches near 2900 and 2200 cm-1, respectively. However, addition of NH2 dramatically affects the spectrum of the neutral. Not only are the NH2 modes themselves strong, but this electron-withdrawing group induces sufficient partial charge on the ring to give the neutral molecule spectra characteristics of the anthracene cation. The sum of the absolute intensities is about four times larger for 2-aminoanthracene than those for 9-cyanoanthracene. Substituting nitrogen in the ring at the nine position (acridine) does not greatly alter the spectrum compared with anthracene.

NASA Discipline Exobiology↗
NASA NTRSDate not supplied

Machine Learning Approaches for Rare-Earth Silicate Environmental Barrier Coating Thermochemical and Thermomechanical Property Predictions

Environmental barrier coatings (EBCs) are a necessary enabling technology for the transition from superalloys to silicon carbide (SiC) ceramic matrix composites (CMCs) in gas turbine engines for increased efficiency and decreased fuel costs. SiC-based CMCs are prone to oxidation-based degradation in the engine hot section, and rare-earth (RE) silicates are promising candidates for EBCs due to their close thermal expansion match to the composite substrate and oxidation resistance. However, the design of EBCs is hindered by the large chemical space of candidate materials and the difficulty in obtaining material properties for engineering optimization. This is especially difficult as research continues into mixed-cation or “high-entropy” RE silicates. First-principles computational methods such as density functional theory (DFT) are highly effective at calculating material properties to guide coating design but are limited by their computational cost. Atomistic simulations have the potential to both accelerate property calculations and expand the properties able to be calculated due to their lower computational compared to DFT. However, they require interatomic potentials (IAPs) specific to the material system of interest, and, to our knowledge, there are no suitable IAPs for RE silicates. Machine learning (ML) is a promising technique to accelerate material property predictions indirectly by generating IAPs for atomistic simulations or via direct prediction. In this work, we present two ML approaches to accelerate the calculation of RE silicate properties relevant to EBC design: 1) a ML-derived interatomic potential (IAP) for atomistic simulations of yttrium disilicate (Y2Si2O7) from DFT training data, and 2) a neural network (NN) model to directly predict thermochemical properties of RE silicates and oxides directly from easily obtainable unit cell parameters. Classical MD simulations using the IAP yield lattice properties and bond lengths in good agreement with both DFT and experimental results from x-ray diffraction. Thermodynamic properties calculated using the finite-displacement phonon method and quasi-harmonic approximation were orders of magnitude faster than DFT with good agreement to the DFT results. The IAP was also used to calculate properties such as coefficient of thermal expansion (CTE) that require large simulation supercells and are therefore difficult with DFT. The IAP correctly predicted the anisotropic nature of the CTE in three different phases of Y2Si2O7. The NN model predicts constant pressure heat capacity, Cp, orders of magnitude faster than DFT calculations, which can enable its use as a surrogate model for multiscale simulations. The two methods presented in this work demonstrate the utility of ML for accelerating the prediction of RE silicate properties, which can in turn accelerate EBC design and optimization.

machine learning↗
NASA NTRS2023

Influence of Cation Species on Thermal Expansion of Y2Si2O7–Gd2Si2O7 Solid Solutions

Mixtures of Y 2 Si 2 O 7 and Gd 2 Si 2 O 7 were synthesized by solid-state reaction at 1600°C and characterized via in situ x-ray diffraction (XRD) to determine their coefficients of thermal expansion (CTE). All solid solutions within the system exhibited the orthorhombic δ-RE 2 Si 2 O 7 (Pna2 1 ) structure. Thermal expansion measurements of Y 2 Si 2 O 7 and Gd 2 Si 2 O 7 correlated well with reported values in literature, and all synthesized solid solutions exhibited CTEs between Y 2 Si 2 O 7 and Gd 2 Si 2 O 7 . Generally, there was a slight decrease in CTE exhibited by the materials with increasing Gd 2 Si 2 O 7 content, with Gd 2 Si 2 O 7 having the lowest CTEs and Y 2 Si 2 O 7 the highest CTEs. The decrease in CTE was attributed to stronger bonds of Gd-O over Y-O, as determined by calculated crystal orbital Hamilton populations using density functional theory. However, such differences were very small and crystal structure was the dominating factor in CTE trends.

rare earth silicates↗
NASA NTRS1992

Electronic and total energy properties of ternary and quaternary semiconductor compounds, alloys, and superlattices: Theoretical study of Cu/graphite bonding

The goals of the research were to provide a fundamental science basis for why the bonding of Cu to graphite is weak, to critically evaluate the previous analysis of the wetting studies with particular regard to the values used for the surface energies of Cu and graphite, and to make recommendations for future experiments or other studies which could advance the understanding and solution of this technological problem. First principles electronic structure calculations were used to study the problem. These are based on density functional theory in the local density approximation and the use of the linear muffin-tin orbital band structure method. Calculations were performed for graphite monolayers, single crystal graphite with the hexagonal AB stacking, bulk Cu, Cu(111) surface, and Cu/graphite superlattices. The study is limited to the basal plane of graphite because this is the graphite plane exposed to Cu and graphite surface energies and combined with the measured contact angles to evaluate the experimental adhesion energy.

Lambrecht, Walter R. L.↗
NASA NTRS1996

On the Reaction of FNO2 with CH3, t-butyl, and C13H21

Theoretical studies are reported for the reaction of FN02 with the radicals CH3, t-butyl, and C13H21, which are templates for the radical site of a hydrogenated diamond surface. All structures axe fully optimized using density functional theory (DFT) based on the B3LYP functional. Calibration calculations axe performed for CH3 + FNO2 using the coupled cluster approach, the the internally contracted multireference configuration. interaction method, and second order perturbation theory based upon complete active space SCF (CASSCF) reference wave function. These calibration calculations support the B3LYP approach for the calculation of bond energies, but show the B3LYP barrier is too low. Combining the calibration calculations with the larger clusters yields our best estimate of a barrier of about 10 kcal/mol. for the reaction of FNO2 with a radical site on hydrogenated diamond.

Thuemmel, H. T.↗
NASA NTRS2002

Effects of Structural Deformation and Tube Chirality on Electronic Conductance of Carbon Nanotubes

A combination of large scale classical force-field (UFF), density functional theory (DFT), and tight-binding Green's function transport calculations is used to study the electronic properties of carbon nanotubes under the twist, bending, and atomic force microscope (AFM)-tip deformation. We found that in agreement with experiment a significant change in electronic conductance can be induced by AFM-tip deformation of metallic zigzag tubes and by twist deformation of armchair tubes. The effect is explained in terms of bandstructure change under deformation.

Svizhenko, Alexei↗
NASA NTRS1998

Modeling of Branched (L, T and Y) Carbon Nanotubes

Models for connecting two or three carbon nanotubes (CNT) using topological defects (i.e., pentagons and heptagons) are presented for the characterization of experimentally observed L, T and Y CNT junctions. The effects of the separation and orientation of the topological defects on the structures and energetics of these junctions are investigated using the nonlocal density function theory (DFT) and semi-empirical molecular orbital (AM1) calculations, and the Brenner empirical potential molecular mechanics simulations. The potential applications of L, Y and T CNT junctions in nanoelectronic devices are also discussed.

Han, Jie↗