Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional 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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Reversible phase transformations between Pb nanocrystals and a viscous liquid-like phase

Phase transformations have been a prominent topic of study for both fundamental and applied science. Solid-liquid reaction–induced phase transformations can be hard to characterize, and the transformation mechanisms are often not fully understood. Here, we report reversible phase transformations between a metal (Pb) nanocrystal and a viscous liquid-like phase unveiled by in situ liquid cell transmission electron microscopy. The reversible phase transformations are obtained by modulating the electron current density (between 1000 and 3000 electrons Å –2 s –1 ). The metal-organic viscous liquid-like phase exhibits short-range ordering with a preferred Pb-Pb distance of 0.5 nm. Assisted by density functional theory and molecular dynamics calculations, we show that the viscous liquid-like phase results from the reactions of Pb with the CH 3 O fragments from the triethylene glycol solution under electron beam irradiation. Such reversible phase transformations may find broad implementations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Multiple Slater Determinants and Strong Spin-Fluctuations as Key Ingredients of the Electronic Structure of Electron- and Hole-Doped Pb 10− x Cu x (PO 4 ) 6 O

LK-99, with chemical formula Pb 10−x Cu x (PO 4 ) 6 O, was recently reported to be a room-temperature superconductor. While this claim has met with little support in a flurry of ensuing work, a variety of calculations (mostly based on density-functional theory) have demonstrated that the system possesses some unusual characteristics in the electronic structure, in particular flat bands. We have established previously that within DFT, the system is insulating with many characteristics resembling the classic cuprates, provided the structure is not constrained to the P3(143) symmetry nominally assigned to it. Here we describe the basic electronic structure of LK-99 within self-consistent many-body perturbative approach, quasiparticle self-consistent GW (QSGW) approximation and their diagrammatic extensions. QSGW predicts that pristine LK-99 is indeed a Mott/charge transfer insulator, with a bandgap gap in excess of 3 eV, whether or not constrained to the P3(143) symmetry. When Pb 9 Cu(PO 4 ) 6 O is hole-doped, the valence bands modify only slightly, and a hole pocket appears. However, two solutions emerge: a high-moment solution with the Cu local moment aligned parallel to neighbors, and a low-moment solution with Cu aligned antiparallel to its environment. In the electron-doped case the conduction band structure changes significantly: states of mostly Pb character merge with the formerly dispersionless Cu d state, and high-spin and low spin solutions once again appear. Thus we conclude that with suitable doping, the ground state of the system is not adequately described by a band picture, and that strong correlations are likely. Irrespective of whether this system class hosts superconductivity or not, the transition of Pb 10 (PO 4 ) 6 O from being a band insulator to Pb 9 Cu(PO 4 ) 6 O, a Mott insulator, and multi-determinantal nature of doped Mott physics make this an extremely interesting case-study for strongly correlated many-body physics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1975

Thermal conductivity of comets

The thermal conductivity of two comet models is calculated. Both models assume the comets to be heterogeneous in composition. One model considers the comet to have a nucleus of water-ice mixed with basaltic and meteoritic material in the form of dust and agglomerated particles surrounded by a layer of water-frost. The second model assumes the frost and water-ice layers have sublimated, leaving a porous structure to some depth composed of solid basaltic and meteoritic material with residual gases. K values are calculated as a function of depth, density, temperature and porosity.

Kumar, G. N.↗
NASA NTRS1975

Satellite disintegration dynamics

The subject of satellite disintegration is examined in detail. Elements of the orbits of individual fragments, determined by DOD space surveillance systems, are used to accurately predict the time and place of fragmentation. Dual time independent and time dependent analyses are performed for simulated and real breakups. Methods of statistical mechanics are used to study the evolution of the fragment clouds. The fragments are treated as an ensemble of non-interacting particles. A solution of Liouville's equation is obtained which enables the spatial density to be calculated as a function of position, time and initial velocity distribution.

Dasenbrock, R. R.↗
NASA NTRS1979

New atomic data of O/2+/

New atomic data are tabulated for the ion O(2+). Collision strengths are calculated for several energies of the exciting electron. The populations of the levels of O(2+) are calculated as a function of electron density under conditions appropriate for the solar atmosphere. The available solar data are compared with theoretical predictions of relative line intensities.

Bhatia, A. K.↗
NASA NTRS1980

New atomic data for Fe/+19/

New atomic data are tabulated for Fe(+19). The electron collision strengths are calculated using the 'distorted wave' approximation and results are presented for several energies of the exciting electron. The populations of the ground levels are calculated as a function of electron density under conditions appropriate to solar flares and tokamak plasmas. Theoretical intensity ratios are tabulated for the UV lines.

Bhatia, A. K.↗
NASA NTRS1983

Theoretical studies of chromospheres and winds in cool stars

The formation of spectral lines in expanding spherical atmospheres was determined in a physically realistic way, taking into account multilevel atomic processes, partial frequency redistribution, and other non-LTE transfer effects that affect the formation of optically thick lines. The formation of MgII and Ca II circumstellar absorption lines in late type giants and supergiants is investigated. The radiative cooling rate as a function of density and temperature was calculated from the results of plane parallel chromospheric models and these results were used to approximate the radiative cooling in an extended wind. The run of temperature was calculated along with the density and velocity profiles. The most important prediction of these models is that a warm zone in the wind must exist as a result of the wave heating. Within this zone, the Ca II and Mg II atoms can be ionized to Ca III and Mg III, so that the gas is transparent in the resonance transitions.

Dupree, A.↗
NASA NTRS1985

Planetary Interior Modeling and Tectonic Implications

A technique is described for estimating spectral admittance functions using Pioneer Venus gravity and topography data. These admittance functions provide a convenient means to carry out systematic geophysical studies over much of the surface of Venus with a variety of interior density models. The admittance functions are calculated in the observation space of line-of-sight (LOS) gravity. Both closed and open system petrological models are considered for the Tharsis region of Mars. An analytic theory for isostatic compensation on one-plate planet is applied, including membrane stresses in the lithosphere, self gravitation, and rotational ellipticity. Crucial to this stress modeling and also to the petrological modeling is the observation that the earliest fracturing seen in the Tharsis region is associated with isostatic stresses. The radial fractures that extend far from Tharsis are associated with an additional and/or a completely different mechanism.

Phillips, R. J.↗
NASA NTRS1991

The electronic Hamiltonian for cuprates

A realistic many-body Hamiltonian for the cuprate superconductors should include both copper d and oxygen p states, hopping matrix elements between them, and Coulomb energies, both on-site and inter-site. We have developed a novel computational scheme for deriving the relevant parameters ab initio from a constrained occupation local density functional. The scheme includes numerical calculation of appropriate Wannier functions for the copper and oxygen states. Explicit parameter values are given for La2CuO4. These parameters are generally consistent with other estimates and with the observed superexchange energy. Secondly, we address whether this complicated multi-band Hamiltonian can be reduced to a simpler one with fewer basis states per unit cell. We propose a mapping onto a new two-band effective Hamiltonian with one copper d and one oxygen p derived state per unit cell. This mapping takes into account the large oxygen-oxygen hopping given by the ab initio calculations.

Annett, James F.↗
NASA NTRS2001

Resonance Scattering of Fe XVII X-ray and EUV Lines

Over the years a number of calculations have been carried out to derive intensities of various X-ray and EUV lines in Fe XVII to compare with observed spectra. The predicted intensities have not agreed with solar observations, particularly for the line at 1.5.02 Angstroms; resonance scattering has been suggested as the source for much of the disagreement. The atomic data calculated earlier used seven configurations having n=3 orbitals and the scattering calculations were carried out only for incident energies above the threshold of the highest fine-structure level. These calculations have now been extended to thirteen configurations having n=4 orbitals and the scattering calculations are carried out below as well as above the threshold of the highest fine structure level. These improved calculations of Fe XVII change the intensity ratios compared to those obtained earlier, bringing the optically thin F(15.02)/F(16.78) ratio and several other ratios closer to the observed values. However, some disagreement with the solar observations still persists, even thought the agreement of the presently calculated optically thin F(15.02)/F(15.26) ratio with the experimental results of Brown et al. (1998) and Laming et al. (2000) has improved. Some of the remaining discrepancy is still thought to be the effect of opacity, which is consistent with expected physical conditions for solar sources. EUV intensity ratios are also calculated and compared with observations. Level populations and intensity ratios are calculated, as a function of column density of Fe XVII, in the slab and cylindrical geometries. As found previously, the predicted intensities for the resonance lines at 15.02 and 15.26 Angstroms exhibit initial increases in flux relative to the forbidden line at 17.10 Angstroms and the resonance line at 16.78 Angstroms as optical thickness increases. The same behavior is predicted for the lines at 12.262 and 12.122 Angstroms. Predicted intensities for some of the allowed EUV lines are also affected by opacity.

Bhatia, A. K.↗
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 NTRS2013

Very Large Eddy Simulations of a Jet-A Spray Reacting Flow in a Single Element LDI Injector With and Without Invoking an Eulerian Scalar DWFDF Method

This paper presents the very large eddy simulations (VLES) of a Jet-A spray reacting flow in a single element lean direct injection (LDI) injector by using the National Combustion Code (NCC) with and without invoking the Eulerian scalar DWFDF method, in which DWFDF is defined as the density weighted time filtered fine grained probability density function. The flow field is calculated by using the time filtered compressible Navier-Stokes equations (TFNS) with nonlinear subscale turbulence models, and when the Eulerian scalar DWFDF method is invoked, the energy and species mass fractions are calculated by solving the equation of DWFDF. A nonlinear subscale model for closing the convection term of the Eulerian scalar DWFDF equation is used and will be briefly described in this paper. Detailed comparisons between the results and available experimental data are carried out. Some positive findings of invoking the Eulerian scalar DWFDF method in both improving the simulation quality and maintaining economic computing cost are observed.

Shih, Tsan-Hsing↗
DOE OSTI2024

BEAST DB: Grand-Canonical Database of Electrocatalyst Properties

We present BEAST DB, an open-source database comprised of ab initio electrochemical data computed using grand-canonical density functional theory in implicit solvent at consistent calculation parameters. The database contains over 20,000 surface calculations and covers a broad set of heterogeneous catalyst materials and electrochemical reactions. Calculations were performed at self-consistent fixed potential as well as constant charge to facilitate comparisons to the computational hydrogen electrode. This article presents common use cases of the database to rationalize trends in catalyst activity, screen catalyst material spaces, understand elementary mechanistic steps, analyze the electronic structure, and train machine learning models to predict higher fidelity properties. Users can interact graphically with the database by querying for individual calculations to gain a granular understanding of reaction steps or by querying for an entire reaction pathway on a given material using an interactive reaction pathway tool. BEAST DB will be periodically updated, with planned future updates to include advanced electronic structure data, surface speciation studies, and greater reaction coverage.

database↗
NASA NTRS1975

Short-term time variability of Cygnus X-1

Autocorrelation functions, power-density spectra, and cross-correlation functions were calculated for 71 observations of Cyg X-1 by the Uhuru satellite. The average autocorrelation function and power-density spectrum were found to be independent of energy and consistent with a shot-noise source with a characteristic time constant of 0.45 s. The data in different energy bands were found to be 100 percent correlated.

Weisskopf, M. C.↗
DOE OSTI2025

Development of a multiphase equation of state for gallium with experiments and ab initio free-energy calculations

Here, we present a five-phase equation of state (EOS) for elemental gallium (Ga) that is developed using both experimental data and new theoretical predictions. Four experimentally observed solid phases (Ga-I, Ga-II, Ga-III, and Ga-IV) and one liquid phase are included. To improve our understanding of the thermal behavior of Ga and its phase boundaries under compression, we have performed ab initio density functional theory (DFT) free-energy calculations for the Ga-III and liquid phases, which enables us to determine the melt temperature to be 2214 ± 100 K at 110 GPa, extending significantly beyond the existing experimental melt data, which are limited to 25 GPa. In order to best describe the electron-thermal contribution, which dominates the liquid free energy at high temperatures, we have carried out averaged-atom-in-jellium DFT calculations to cover the entire temperature and density range of the EOS. The resulting multiphase Ga EOS is able to accurately reproduce a diverse variety of data, including known phase boundaries, the principal Hugoniot, low-pressure liquid isobars, and diamond-anvil-cell isotherm measurements at high pressures. It agrees more closely with key experimentally measured properties than other Ga EOS models targeted for high-pressure applications.

Wu, Christine J. [Lawrence Livermore National Labo↗