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 73 records · Page 4

Prediction of molecular infrared band and line intensities

The general problem of predicting infrared intensities is examined. The problem of predicting line intensities is discussed with attention focussed specifically on the problem of predicting vibrational strengths (or total band intensities). Two methods are considered: predictions based on ab initio quantum mechanical calculations, and predictions based on the transfer of atomic polar tensor parameters from small model molecules to the molecule of interest. Comparisons of predicted spectra with experimental spectra are shown for CH3OH, CH3OCH3, C2H5OH, and HOC2. It is suggested that the most important application of this work may be to the prediction of intensities for free radicals, ions, and reactive intermediates so that the infrared spectrum can be used in the quantitative analysis for these interesting species.

Person, W. B.↗

Theoretical study of NH2 - Potential curves, transition moments, and photodissociation cross sections

Photodissociation cross sections from the ground state of NH2 have been calculated using a pseudodiatomic model. The potential curves needed in these calculations, functions of one NH bond length, were obtained by ab initio MCSCF calculations on the five lowest doublet states of NH2 using a contracted Gaussian basis set of double zeta quality augmented by polarization and Rydberg functions. Transition dipole moments between the ground and excited states were evaluated using MCSCF wave functions.

Saxon, R. P.↗

The nature of the metal-CO interaction and bonding

The adsorption of CO on metal surfaces is represented by molecular orbital cluster models of CO at an on top site and adsorbed normal to the surface carbon end down. Ab initio SCF and MCSCF calculations are performed for several clusters. The new constrained space orbital variation CSOV approach is used to analyze the bonding and to compare CO adsorption on Al, representative of sp metals, with that on Cu, representative of transition metals. There is a large repulsion between the superposed free CO and metal charge distributions which is considerably smaller for Cu than for Al because there are fewer valence sigma electrons for Cu than for Al. The CSOV analysis shows that the metal to CO pi donation is much more important than the CO to metal sigma donation. It is also shown that for Cu, the d pi contribution to the metal pi donation is larger than the valence 4p pi contribution. The d pi donation is compared between Fe, Ni, and Cu and this donation and the metal-CO interaction are found to be different in the order Fe greater than Ni greater than Cu.

Bagus, P. S.↗

Scandium and Chromium in the Strontium Filament in the Homunculus of eta Carinae

We continue a systematic study of chemical abundances of the Strontium Filament found in the ejecta of eta Carinae. To this end we interpret the emission spectrum of Sc II and Cr II using multilevel non-LTE models of these systems. Since the atomic data for these ions was previously unavailable, we carry out ab initio calculations of radiative transition rates and electron impact excitation rate coefficients. The observed spectrum is emitted from a mostly neutral region with electron density of the order of 10(exp 7) cm (exp -3) and a temperature between 6000 and 7000 K. These conditions are consistent with our previous diagnostics from [Ni II], [Ti II], amd [Sr II]. The observed spectrum indicates an abundance of Sc relative Ni that more than 40 times the solar values, while the Cr/Ni abundance ratio is roughly solar. Various scenarios of depletion and dust destruction are suggested to explain such abnormal abundances.

Gull, T.R.↗

Macromolecular recognition: Structural aspects of the origin of the genetic system

Theoretical simulation of prebiotic chemical processes is an invaluable tool for probing the phenomenon of the evolution of life. Using computational and modeling techniques and guided by analogies from present day systems, we seek to understand the emergence of the genetic apparatus, enzymatic catalysis and protein synthesis under prebiotic conditions. Modeling of the ancestral aminoacyl-tRNA-synthetases (aRS) may provide important clues to the emergence of the genetic code and the protein synthetic machinery. The minimal structural requirements for the catalysis of tRNA aminoacylation are being explored. A formation of an aminoacyl adenylate was studied in the framework of ab initio molecular orbital theory. The role of individual residues in the vicinity of the TyrRS active site was examined, and the effect of all possible amino acids substitutions near the active site was examined. A formation of aminoacyl tRNA was studied by the molecular modeling system SYBYL with the high resolution crystallographic structures of the present day tRNA, aRS's complexes. The ultimate goal is to propose a simple RNA segment that is small enough to be build in the primordial chemical environment but maintains the specificity and catalytic activity of the contemporary RNA enzyme. To understand the mechanism of ribozyme catalyzed reactions, ab initio and semi-empirical (ZINDO) programs were used to investigate the reaction path of transphosphorylation. A special emphasis was placed on the possible catalytic and structural roles played by the coordinated magnesium cation. Both the inline and adjacent mechanisms of transphosphorylation were studied. The structural characteristics of the target helices, particularly a possible role for the G-T pair, is also studied by a molecular dynamics (MD) simulation technique.

Rein, Robert↗

Summary of Research/Publications

Summary of research/publications include:(1) Comment on broadening of water microwave lines by collisions with helium atoms; (2) Calculations of ion-molecule deuterium fractionation reactions involving HD; (3) Ab initio predictions on the rotational spectra of carbon-chain carbene molecules; (4) Theoretical IR spectra of ionized naphthalene; (5) Improved collisional excitation rates for interstellar water; (6) Calculations on the competition between association and reaction for C3H+ + H2; (7) Theoretical infrared spectra of some model polycyclic aromatic hydrocarbons: effect of ionization; (8) Calculations concerning interstellar isomeric abundance ratios for C3H and C3H2; (9) New calculations on the ion-molecule processes C2H2+ + H2 C2H3+ + H and C2H2+ + H2 C2H4+; (10) Anisotropic rigid rotor potential energy function for H2O-H2; (11) A correlated ab initio study of linear carbon-chain radicals CnH (n=2-7); (12) Ab initio characterization of MgCCH, MgCCH+, and MgC2 and pathways to their formation in the interstellar medium; (13) Why HOC+ is detectable in interstellar clouds: The rate of the reaction between HOC+ and H2; (14) A correlated ab initio study of the X 2A 1 and A 2E states of MgCH3; (15) On the stability of interstellar carbon clusters: The rate of the reaction between C3 and O; and (16) The rate of the reaction between CN and C2H2 at interstellar temperatures.

Source record↗

Experimental and Ab Initio Studies of the HDO Absorption Spectrum in the 13165-13500 1/cm Spectral Region

The HDO absorption spectrum has been recorded in the 13165 - 13500 cm(exp-1) spectral region by Intracavity Laser Absorption Spectroscopy. The spectrum (615 lines), dominated by the 2n2 + 3n3 and n1+3n3 bands was assigned and modeled leading to the derivation of 196 accurate energy levels of the (103) and (023) vibrational states. Finally, 150 of these levels have been reproduced by an effective Hamiltonian involving two vibrational dark states interacting with the (023) and ( 103) bright states. The rms deviation achieved by variation of 28 parameters is 0.05-1 cm, compared to an averaged experimental uncertainty of 0.007-1 cm, indicating the limit of validity of the effective Hamiltonian approach for HDO at high vibrational excitation. The predictions of previous ab initio calculations of the HDO spectrum were extensively used in the assignment process. The particular spectral region under consideration has been used to test and discuss the improvements of new ab initio calculations recently performed on the basis of the same potential energy surface but with an improved dipole moment surface. The improvements concern both the energy levels and the line intensities. In particular, the strong hybrid character of the n1+3n3 band is very well accounted for by the the new ab initio calculations.

Schwenke, David↗

Intermolecular potential for thermal H2O-He collisions

Theoretical potentials for rotational excitation of H2O by He were constructed via several methods, all of which start with a large basis set SCF interaction. The semiempirical Hartree-Fock with damped dispersion model adds a damped long-range attraction with parameters adjusted to fit experimental total differential cross sections. Purely ab initio potentials add correlation energies obtained via perturbation theory (MP2 and MP4) or a variational method (ICF1). Scattering calculations were performed on all surfaces to compare wih available beam scattering and pressure broadening data and to assess sensitivity of state-to-state rates to uncertainties in the potential. From comparison with the limited experimental data, the ICF1 surface appears to be marginally better than the MP4 surface. Thermal rates calculated from this surface should be accurate to better than 50 percent, at least for the larger, more important rates.

Palma, Amedeo↗

Hierarchical screening for Li-based solid electrolytes using fast, interpretable machine-learned potentials

Li-based solid-state electrolyte materials enable safer, all-solid-state batteries but the computational search for candidates with favorable stability and Li-ion conductivity is challenging due to the size of the search space and the cost of evaluating transport properties with ab initio methods. The prohibitive cost of high-throughput screening with DFT has lead to the development of surrogate models using geometric analysis, empirical potentials, and descriptors for ionic transport. Here, I will discuss a hierarchical screening approach for identifying promising materials using a combination of density functional theory, bond-valence methods, and machine learning potentials generated with the Ultra-Fast Force Fields (UF3) framework. We show how the inexpensive bond-valence method can be used to guide the generation of training samples for machine learning, in addition to filtering candidates. Finally, we apply the hierarchical workflow to screen for ionic conductivity across a database of Li-containing compounds.

Materials discovery↗

Hierarchical Screening for Li-Based Solid Electrolytes Using Fast, Interpretable Machine-Learned Potentials

Li-based solid-state electrolyte materials enable safer, all-solid-state batteries but the computational search for candidates with favorable stability and Li-ion conductivity is challenging due to the size of the search space and the cost of evaluating transport properties with ab initio methods. The prohibitive cost of high-throughput screening with DFT has lead to the development of surrogate models using geometric analysis, empirical potentials, and descriptors for ionic transport. Here, I will discuss a hierarchical screening approach for identifying promising materials using a combination of density functional theory, bond-valence methods, and machine learning potentials generated with the Ultra-Fast Force Fields (UF3) framework. We show how the inexpensive bond-valence method can be used to guide the generation of training samples for machine learning, in addition to filtering candidates.

Materials discovery↗

Ab initio study of the alkali and alkaline-earth monohydroxides

A systematic study of the structures and dissociation energies of all the alkali and alkaline-earth monohydroxides is conducted. A theoretical model for determining accurate dissociation energies of ionic molecules is discussed. The obtained theoretical structures and dissociation energies of the alkali and alkaline-earth monohydroxides, respectively, are compared with experimental data. It is found that the theoretical studies of the bending potentials of BeOH, MgOH, and CaOH reveal the different admixture of covalent character in these systems. The BeOH molecule with the largest degree of covalent character is found to be bent (theta equals 147 deg). The MgOH is also linear. The theoretical dissociation energies for the alkali and akaline-earth hydroxides are thought to be accurate to 0.1 eV.

Bauschlicher, C. W., Jr.↗

Realistic Modeling of Multi-Scale MHD Dynamics of the Solar Atmosphere

Realistic 3D radiative MHD simulations open new perspectives for understanding the turbulent dynamics of the solar surface, its coupling to the atmosphere, and the physical mechanisms of generation and transport of non-thermal energy. Traditionally, plasma eruptions and wave phenomena in the solar atmosphere are modeled by prescribing artificial driving mechanisms using magnetic or gas pressure forces that might arise from magnetic field emergence or reconnection instabilities. In contrast, our 'ab initio' simulations provide a realistic description of solar dynamics naturally driven by solar energy flow. By simulating the upper convection zone and the solar atmosphere, we can investigate in detail the physical processes of turbulent magnetoconvection, generation and amplification of magnetic fields, excitation of MHD waves, and plasma eruptions. We present recent simulation results of the multi-scale dynamics of quiet-Sun regions, and energetic effects in the atmosphere and compare with observations. For the comparisons we calculate synthetic spectro-polarimetric data to model observational data of SDO, Hinode, and New Solar Telescope.

Realistic Modeling↗

Downfolding Complex Materials Problems Onto Model Hamiltonians for Quantum Computers

Simulating the properties of quantum materials is expected to be one of the exciting applications for quantum computers and where we hope to see advantages over classical hardware. The complexity of ab initio Hamiltonians describing the physics of application-relevant materials places them beyond the realm of possibility for solution on near-term hardware with a limited number of qubits. Various Hamiltonian approximations, including Hamiltonian downfolding, offers a possibility towards simulating complex materials on near-term hardware. This is accomplished by approximating the relevant physics of a given material through their representation by simpler model Hamiltonians, such as the Hubbard Hamiltonian or extensions of it. Here we employ a well-defined first-principles methodology for deriving downfolded multi-band extended Hubbard Hamiltonians of materials, capturing strong electronic correlation and electron-phonon coupling, based on the formalism of Wannier functions and the calculation of the screened Coulomb interaction. We demonstrate for a variety of systems that quantum simulation of these downfolded Hamiltonians reproduces key properties, thus establishing downfolding as a promising route to achieve near-term simulation of application-relevant systems on quantum hardware.

Antonios Markos Alvertis↗

Downfolding complex materials problems onto model Hamiltonians for quantum computers

Simulating the properties of quantum materials is expected to be one of the exciting applications for quantum computers and where we hope to see advantages over classical hardware. The complexity of ab initio Hamiltonians describing the physics of application-relevant materials places them beyond the realm of possibility for solution on near-term hardware with a limited number of qubits. Various Hamiltonian approximations, including Hamiltonian downfolding, offers a possibility towards simulating complex materials on near-term hardware. This is accomplished by approximating the relevant physics of a given material through their representation by simpler model Hamiltonians, such as the Hubbard Hamiltonian or extensions of it. Here we employ a well-defined first-principles methodology for deriving downfolded multi-band extended Hubbard Hamiltonians of materials, capturing strong electronic correlation and electron-phonon coupling, based on the formalism of Wannier functions and the calculation of the screened Coulomb interaction. We demonstrate for a variety of systems that quantum simulation of these downfolded Hamiltonians reproduces key properties, thus establishing downfolding as a promising route to achieve near-term simulation of application-relevant systems on quantum hardware.

Antonios M. Alvertis↗

Analysis of the far infrared H2-He spectrum

Previous measurements of the far infrared absorption due to H2-He collisions at the temperatures of 77, 195, and 292 K are analyzed. The spectra are fitted by a semiempirical line shape representing the isotropic induced overlap component and combined anisotropic quadrupolar and overlap components. The experimental spectral moments are evaluated and compared with theory for several induced-dipole and potential models. From the isotropic contribution, the range and strength of the induced dipole is evaluated and compared with the results of ab initio calculations. Fitting parameters are obtained with physically plausible temperature dependences which allow simple and accurate representation of the spectra and of their moments at temperatures different from those of the measurements.

Cohen, E. R.↗

The effect of ionization on the infrared absorption spectra of PAHs: A preliminary report

The emission lines observed in many interstellar IR sources at 3.28, 6.2, 7.7, 8.7, and 11.3 microns are theorized to originate from polycyclic aromatic hydrocarbons (PAHs). These assignments are based on analyses of lab IR spectra of neutral PAHs. However, it is likely that in the interstellar medium that PAHs are ionized, i.e., are positively charged. Besides, as pointed out by Allamandola et al., although the IR emission band spectrum resembles what one might expect from a mixture of PAHs, it does not match in details such as frequency, band profile, or relative intensities predicted from the absorption spectra of any known PAH molecule. One source of more information to test the PAH theory is ab initio molecular orbital theory. It can be used to compute, from first principles, the geometries, vibrational frequencies, and vibrational intensities for model PAH compounds which are difficult to study in the lab. The Gaussian 86 computer program was used to determine the effect of ionization on the infrared absorption spectra of several small PAHs: naphthalene and anthracene. A preliminary report is presented of the results of these calculations.

Defrees, Doug J.↗

Electrostatic Properties of Polymers Subjected to Atmospheric Pressure Plasma Treatment; Correlation of Experimental Results with Atomistic Modeling

this study, PE, PTFE, PS and PMMA were exposed to a He+O2, APGD and pre and post treatment surface chemistries were analyzed by X-ray photoelectron spectroscopy and contact angle measurements. Semi-empirical and ab-initio calculations were performed to correlate the experimental results with sonic plausible molecular and electronic structure features of the oxidation process. For the PE and PS, significant surface oxidation showing C-O, C=O, and O-C=O bonding, and a decrease in the surface contact angles was observed. For the PTFE and PM MA, little change in the surface composition was observed. The molecular modeling calculations were performed on single and multiple oligomers and showed regardless of oxidation mechanism, e.g. -OH, =O or a combination thereof, experimentally observed levels of surface oxidation were unlikely to lead to a significant change in the electronic structure of PE and PS, and that the increased hydrophilic properties are the primary reason for the observed changes in its electrostatic behavior. Calculations for PTFE and PMMA argue strongly against significant oxidation of those materials, as confirmed by the XPS results.

Trigwell, S.↗

How Very Massive Metal-Free Stars Start Cosmological Reionization

The initial conditions and relevant physics for the formation of the earliest galaxies are well specified in the concordance cosmology. Using ab initio cosmological Eulerian adaptive mesh refinement radiation hydrodynamical calculations, we discuss how very massive stars start the process of cosmological reionization. The models include nonequilibrium primordial gas chemistry and cooling processes and accurate radiation transport in the case B approximation using adaptively ray-traced photon packages, retaining the time derivative in the transport equation. Supernova feedback is modeled by thermal explosions triggered at parsec scales. All calculations resolve the local Jeans length by at least 16 grid cells at all times and as such cover a spatial dynamic range of approx.10(exp 6). These first sources of reionization are highly intermittent and anisotropic and first photoionize the small-scale voids surrounding the halos they form in, rather than the dense filaments they are embedded in. As the merging objects form larger, dwarf-sized galaxies, the escape fraction of UV radiation decreases and the H II regions only break out on some sides of the galaxies, making them even more anisotropic. In three cases, SN blast waves induce star formation in overdense regions that were formed earlier from ionization front instabilities. These stars form tens of parsecs away from the center of their parent DM halo. Approximately five ionizing photons are needed per sustained ionization when star formation in 10(exp 6) stellar Mass halos is dominant in the calculation. As the halos become larger than approx.10(exp 7) Stellar Mass, the ionizing photon escape fraction decreases, which in turn increases the number of photons per ionization to 15-50, in calculations with stellar feedback only. Radiative feedback decreases clumping factors by 25% when compared to simulations without star formation and increases the average temperature of ionized gas to values between 3000 and 10,000 K.

Wise, John H.↗