Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio modeling”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Ionization and thermal equilibrium models for O star winds based on time-independent radiation-driven wind theory

Ab initio ionization and thermal equilibrium models are calculated for the winds of O stars using the results of steady state radiation-driven wind theory to determine the input parameters. Self-consistent methods are used for the roles of H, He, and the most abundant heavy elements in both the statistical and the thermal equilibrium. The model grid was chosen to encompass all O spectral subtypes and the full range of luminosity classes. Results of earlier modeling of O star winds by Klein and Castor (1978) are reproduced and used to motivate improvements in the treatment of the hydrogen equilibrium. The wind temperature profile is revealed to be sensitive to gross changes in the heavy element abundances, but insensitive to other factors considered such as the mass-loss rate and velocity law. The reduced wind temperatures obtained in observing the luminosity dependence of the Si IV lambda 1397 wind absorption profile are shown to eliminate any prospect of explaining the observed O VI lambda 1036 line profiles in terms of time-independent radiation-driven wind theory.

Drew, J. E.↗
NASA NTRS1974

Virial expansion of the second layer in physical adsorption - An ab initio calculation for helium on argon crystal

A model of helium adsorption on an argon crystal is built up from the premise that local adsorption predominates in the first layer and nonlocal adsorption in the second. Application of the virial expansion theorem to the second layer gives a series in which the first term represents the motion of a single molecule in the external potential field and the second a two-body interaction under this field. The thermodynamic functions of the adsorbed phase are calculated ab initio, the gas-solid interaction potential being derived from lattice summation and the partition function from an appropriate choice of a site-spacing polynomial to describe the periodic potential. The mutual interaction of adsorbed molecules is calculated with a two-dimensional Lennard-Jones potential. The second virial coefficient is calculated and its dependence on temperature and choice of potential is studied. It is found that the second virial coefficient is very well approximated by a two-dimensional gas in free space. The adsorption isotherm, isosteric heat, and specific heat are obtained and compared with the results of Ross and Steele, giving excellent agreement.

Oh, B. K.↗
NASA NTRS1977

Comment on the relation between the nonadiabatic coupling and the complex intersection of potential energy curves

Simple relations are discussed that provide a correspondence between the complex intersection of two potential surfaces and the nonadiabatic coupling matrix element between those surfaces. These are key quantities in semiclassical and quantum mechanical theories of collision induced electronic transitions. Within the two state approximation, the complex intersection is shown to be directly related to the location and magnitude of the peak in the nonadiabatic coupling. Two cases are considered: the avoided crossing between two potential surfaces; and the spin orbit interaction due to a P-2 halogen atom. Comparisons are made between the results of the two-state model and the results of ab initio quantum chemical calculations.

Jaffe, R. L.↗
NASA NTRS1977

Comment on the relation between the nonadiabatic coupling and the complex intersection of potential energy curves

Simple relations are discussed that provide a correspondence between the complex intersection of two potential surfaces and the nonadiabatic coupling matrix element between those surfaces. These are key quantities in semiclassical and quantum-mechanical theories of collision-induced electronic transitions. Within the two-state approximation, the complex intersection is shown to be directly related to the location and magnitude of the peak in the nonadiabatic coupling. Two cases have been considered: (1) the avoided crossing between two potential surfaces, and (2) the spin-orbit interaction due to a 2P halogen atom. Comparisons are made between the results of the two-state model and the results of ab initio quantum chemical calculations.

Jaffe, R. L.↗
NASA NTRS1985

Dynamic Creep Buckling: Analysis of Shell Structures Subjected to Time-dependent Mechanical and Thermal Loading

The objective of the present research is to develop a general mathematical model and solution methodologies for analyzing the structural response of thin, metallic shell structures under large transient, cyclic, or static thermomechanical loads. Among the system responses associated with these loads and conditions are thermal buckling, creep buckling, and ratcheting. Thus geometric and material nonlinearities (of high order) can be anticipated and must be considered in developing the mathematical model. A complete, true ab-initio rate theory of kinematics and kinetics for continuum and curved thin structures, without any restriction on the magnitude of the strains or the deformations, was formulated. The time dependence and large strain behavior are incorporated through the introduction of the time rates of metric and curvature in two coordinate systems: fixed (spatial) and convected (material). The relations between the time derivative and the covariant derivative (gradient) were developed for curved space and motion, so the velocity components supply the connection between the equations of motion and the time rates of change of the metric and curvature tensors.

Simitses, G. J.↗
NASA NTRS1995

Determination of Quantum Chemistry Based Force Fields for Molecular Dynamics Simulations of Aromatic Polymers

Ab initio quantum chemistry calculations for model molecules can be used to parameterize force fields for molecular dynamics simulations of polymers. Emphasis in our research group is on using quantum chemistry-based force fields for molecular dynamics simulations of organic polymers in the melt and glassy states, but the methodology is applicable to simulations of small molecules, multicomponent systems and solutions. Special attention is paid to deriving reliable descriptions of the non-bonded and electrostatic interactions. Several procedures have been developed for deriving and calibrating these parameters. Our force fields for aromatic polyimide simulations will be described. In this application, the intermolecular interactions are the critical factor in determining many properties of the polymer (including its color).

Jaffe, Richard↗
NASA NTRS1991

Ab initio study of C + H3+ reactions

The reaction C + H3+ --> CH(+) + H2 is frequently used in models of dense interstellar cloud chemistry with the assumption that it is fast, i.e. there are no potential energy barriers inhibiting it. Ab initio molecular orbital study of the triplet CH3+ potential energy surface (triplet because the reactant carbon atom is a ground state triplet) supports this hypothesis. The reaction product is 3 pi CH+; the reaction is to exothermic even though the product is not in its electronic ground state. No path has been found on the potential energy surface for C + H3+ --> CH2(+) + H reaction.

NASA Discipline Exobiology↗
NASA NTRS2003

Surface Segregation Energies of BCC Binaries from Ab Initio and Quantum Approximate Calculations

We compare dilute-limit segregation energies for selected BCC transition metal binaries computed using ab initio and quantum approximate energy method. Ab initio calculations are carried out using the CASTEP plane-wave pseudopotential computer code, while quantum approximate results are computed using the Bozzolo-Ferrante-Smith (BFS) method with the most recent parameterization. Quantum approximate segregation energies are computed with and without atomistic relaxation. The ab initio calculations are performed without relaxation for the most part, but predicted relaxations from quantum approximate calculations are used in selected cases to compute approximate relaxed ab initio segregation energies. Results are discussed within the context of segregation models driven by strain and bond-breaking effects. We compare our results with other quantum approximate and ab initio theoretical work, and available experimental results.

Good, Brian S.↗
NASA NTRS2011

Ab initio Computations of the Electronic, Mechanical, and Thermal Properties of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2

Refractory materials such as metallic borides, often considered as ultra high temperature ceramics (UHTC), are characterized by high melting point, high hardness, and good chemical inertness. These materials have many applications which require high temperature materials that can operate with no or limited oxidation. Ab initio, first principles methods are the most accurate modeling approaches available and represent a parameter free description of the material based on the quantum mechanical equations. Using these methods, many of the intrinsic properties of these material can be obtained. We performed ab initio calculations based on density functional theory for the UHTC materials ZrB2 and HfB2. Computational results are presented for structural information (lattice constants, bond lengths, etc), electronic structure (bonding motifs, densities of states, band structure, etc), thermal quantities (phonon spectra, phonon densities of states, specific heat), as well as information about point defects such as vacancy and antisite formation energies.

Lawson, John W.↗
NASA NTRS1993

Applicability of PM3 to transphosphorylation reaction path: Toward designing a minimal ribozyme

A growing body of evidence shows that RNA can catalyze many of the reactions necessary both for replication of genetic material and the possible transition into the modern protein-based world. However, contemporary ribozymes are too large to have self-assembled from a prebiotic oligonucleotide pool. Still, it is likely that the major features of the earliest ribozymes have been preserved as molecular fossils in the catalytic RNA of today. Therefore, the search for a minimal ribozyme has been aimed at finding the necessary structural features of a modern ribozyme (Beaudry and Joyce, 1990). Both a three-dimensional model and quantum chemical calculations are required to quantitatively determine the effects of structural features of the ribozyme on the reaction it catalyzes. Using this model, quantum chemical calculations must be performed to determine quantitatively the effects of structural features on catalysis. Previous studies of the reaction path have been conducted at the ab initio level, but these methods are limited to small models due to enormous computational requirements. Semiempirical methods have been applied to large systems in the past; however, the accuracy of these methods depends largely on a simple model of the ribozyme-catalyzed reaction, or hydrolysis of phosphoric acid. We find that the results are qualitatively similar to ab initio results using large basis sets. Therefore, PM3 is suitable for studying the reaction path of the ribozyme-catalyzed reaction.

Manchester, John I.↗
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 NTRS1971

Interpretation of experimental differential elastic scattering cross section for H/+/ + Ne.

The interatomic V(r) found by potential-model calculations is compared to an ab initio calculation of the intermolecular potential for NeH+ due to Peyerimhoff (1965). The reason for a rather significant difference in the two results is discussed, and a new method due to Remler (1971) involving Regge poles is applied. The method can be used to calculate the differential cross section, where the starting point for such calculations is based on intuition in light of the classical deflection function.

Bobbio, S. M.↗
NASA NTRSDate not supplied

Molten Salt Electrolytes: Computational Analysis of Transport, Electrochemical Properties and Designing New Mixtures

Electric aircraft propulsion has gained a traction over the last decade due to possible high-energy density electrochemistries with reasonable cycle life and identification of non-flammable chemistries that can guarantee much better safety. Commercial Li-ion batteries cannot reach such high capacities because of weight limitations that arise from the widely used intercalation electrodes. Also, use of organic electrolytes make them highly susceptible to fire on exposure to air and humidity. Currently, use of Lithium metal anode along with conversion electrochemistries such as Li-O2 and Li-S are being pursued to achieve such high energy densities. Safer inorganic solid-state electrolytes, recently discovered Water-in-Salt electrolytes, molten salt electrolytes are some of the alternatives being pursued for a safer/non-flammable battery. Of these safer alternatives, molten-salt electrolytes offer some attractive properties: liquid at operating temperatures resulting in better electrode-wetting, stable interface with Li-metal anodes and good ionic conductivity; their only drawback being high operating temperatures. In this talk, we will discuss some of the computational studies, driven via atomistic simulations, of the properties of molten-salt electrolytes including transport mechanism and interface stability. We will also discuss predicting electrochemical properties of these high-temperature electrolytes. Further, we will discuss a thermodynamic approach to predicting new molten-salt eutectics with lower melting points.

Molten Salt Electrolytes↗
NASA NTRS2001

Endo-Fullerene and Doped Diamond Nanocrystallite Based Models of Qubits for Solid-State Quantum Computers

Models of encapsulated 1/2 nuclear spin H-1 and P-31 atoms in fullerene and diamond nanocrystallite, respectively, are proposed and examined with ab-initio local density functional method for possible applications as single quantum bits (qubits) in solid-state quantum computers. A H-1 atom encapsulated in a fully deuterated fullerene, C(sub 20)D(sub 20), forms the first model system and ab-initio calculation shows that H-1 atom is stable in atomic state at the center of the fullerene with a barrier of about 1 eV to escape. A P-31 atom positioned at the center of a diamond nanocrystallite is the second model system, and 3 1P atom is found to be stable at the substitutional site relative to interstitial sites by 15 eV, Vacancy formation energy is 6 eV in diamond so that substitutional P-31 atom will be stable against diffusion during the formation mechanisms within the nanocrystallite. The coupling between the nuclear spin and weakly bound (valance) donor electron coupling in both systems is found to be suitable for single qubit applications, where as the spatial distributions of (valance) donor electron wave functions are found to be preferentially spread along certain lattice directions facilitating two or more qubit applications. The feasibility of the fabrication pathways for both model solid-state qubit systems within practical quantum computers is discussed with in the context of our proposed solid-state qubits.

Park, Seongjun↗
NASA NTRS1986

A model of the Schumann-Runge continuum of O2

Ab initio calculations of the B3Sigma sub u(-) potential energy curve of O2 are reported together with the transition moment connecting the B3Sigma sub u(-) and the ground electronic states. Photodissociation cross sections are presented over the wavelength region 127-152 nm. Small adjustments to the potential energy curve and the transition moment lead to a model of the absorption which agrees well with the experimental cross sections and which reproduces the structural features observed on the high-energy edge. The shoulder at 135.6 nm is attributed to the lowest 3Pi sub u state and the auxiliary maxima to transitions into the B3Sigma sub u(-) state.

Allison, A. C.↗
NASA NTRSDate not supplied

Ab-Initio Study of Stacking Fault Segregation Behavior in Ni-based Superalloys at High Temperatures

Ni-based superalloys demonstrate extraordinary creep properties compared to traditional metallic systems due to the formation of L12 (γ’) precipitates. Despite experimental evidence of alloying element segregation to the stacking faults, prior ab-initio efforts indicated repulsion of Nb from superlattice intrinsic stacking faults (SISFs) . Due to the importance of this element in key strengthening mechanisms, we developed a new technique to account for multi-component segregation behavior at temperature combining ab-initio molecular dynamics (AIMD) and Monte-Carlo (MC) simulations. Using this technique, we replicated experimentally observed segregation behavior of Nb, Co, and Ti to the SISF in L12 (Figure 1). Furthermore, we investigated the co-segregation behaviors of individual pairs of alloying elements to determine the interaction effects which lead to this segregation behavior in the full alloy system. By simulating the reduced-chemistry models using our approach, we reconciled previous ab-initio calculations demonstrating repulsion of Nb from the SISF and the experimentally observed segregation behavior.

Molecular Dynamics↗
NASA NTRS1987

Theoretical approaches to metal chemistry

Theoretical calculations on metals ranging from very accurate ab initio studies of diatomic and triatomic systems to model studies of larger clusters are presented. Recent improvements in the representation of one-particle and n-particle spaces have made possible the prediction that Al2 has a 3Pi(u) ground state, even though the 3Sigma(-)g state lies within 200/cm. Results suggest that cluster geometry varies dramatically with cluster size, and that rather large clusters are required before the bulk structure becomes optimal. Al cluster studies show that three-body terms are needed for quantitative agreement with ab initio calculations.

Bauschlicher, Charles W.↗