Development of a Reduced C3 NEQAIR Radiation Model Based on Ab-Initio Calculations
This work reports an updated model for C3 absorption based on ab initio data.
SEARCH · Search NASA
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.
This work reports an updated model for C3 absorption based on ab initio data.
Thermodynamic functions for Si(OH)4(g) and SiO(OH)2(g) have been measured using the transpiration method. A second law enthalpy of formation and entropy and a third law enthalpy of formation has been calculated for Si(OH)4. The results are in very good agreement with previous experimental measurements, ab-initio calculations, and estimates.
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.
The results of ab-initio calculations are reported for (1) small transition metal clusters and (2) potential energy surfaces for chemical reactions important in hydrogen combustion and high temperature air chemistry.
Dissociative recombination (DR) of molecular ions with electrons has important consequences in many areas of physical science. Ab-initio calculations coupled with resonant scattering theory and multichannel quantum defect studies have produced detailed results illuminating the role of ion vibrational excitation, the quantum yields of the DR products, and the role of Rydberg states. The theoretical and experimental results are discussed.
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.
We have resolved the long-standing problem consisting of 30%-50% theoretical underestimates of the band gaps of non-metallic materials. We describe the Bagayoko, Zhao, and Williams (BZW) method that rigorously circumvents the basis-set and variational effect presumed to be a cause of these underestimates. We present ab-initio, computational results that are in agreement with experiment for diamond (C), silicon (Si), silicon carbides (3C-SiC and 4H-SiC), and other semiconductors (GaN, BaTiO3, AlN, ZnSe, ZnO). We illustrate the predictive capability of the BZW method in the case of the newly discovered cubic phase of silicon nitride (c-Si3N4) and of selected carbon nanotabes [(10,0), and (8,4)]. Our conclusion underscores the inescapable need for the BZW method in ab-initio calculations that employ a basis set in a variational approach. Current nanoscale trends amplify this need. We estimate that the potential impact of applications of the BZW method in advancing our understanding of nonmetallic materials, in informing experiment, and particularly in guiding device design and fabrication is simply priceless.
It is well known that liquid oxygen has a sufficient paramagnetic susceptibility that a strong magnetic field gradient can lift it in the earth's gravitational field. The movement of liquid oxygen is vital to the space program since it one of the primary oxidizers used for propulsion. Transport of liquid oxygen (LOX) via direct interaction of the magnetic fields (B field) with the fluid is a current topic of research and development at Kennedy Space Center, FL. This method of transporting (i.e. pumping) LOX may have particular advantages on Mars and other reduced gravitational environments, namely safety and reliability. This paper will address transport of a magnetic fluid, LOX, via phased-pulsed electromagnets acting on the edge of the column of fluid. The authors have developed a physical model from first-principles for the motion of a magnetic fluid in a particular U-tube geometry subjected to a pulsed magnetic field from an arbitrary solenoidal electromagnet. Experimental data that have been collected from the analogous geometry correlate well to that of the ab-initio calculations.
Although not yet identified in the interstellar medium (ISM), N-heterocycles including nucleobases the information subunits of DNA and RNA are present in carbonaceous chondrites, which indicates that molecules of biological interest can be formed in non-terrestrial environments via abiotic pathways. Recent laboratory experiments and ab-initio calculations have already shown that the irradiation of pyrimidine in pure H2O ices leads to the formation of a suite of oxidized pyrimidine derivatives, including the nucleobase uracil. In the present work, NH3:pyrimidine and H2O:NH3:pyrimidine ice mixtures with different relative proportions were irradiated with UV photons under astrophysically relevant conditions. Liquid- and gas-chromatography analysis of the resulting organic residues has led to the detection of the nucleobases uracil and cytosine, as well as other species of prebiotic interest such as urea and small amino acids. The presence of these molecules in organic residues formed under abiotic conditions supports scenarios in which extraterrestrial organics that formed in space and were subsequently delivered to telluric planets via comets and meteorites could have contributed to the inventory of molecules that triggered the first biological reactions on their surfaces.
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.
We report the detection in the circumstellar envelope IRC+10216 of millimeter lines of the rare isotopomers (25)MgNC and (26)MgNC, as well as of a line at 234433 MHz, which could be the J= 7-6 transition of (26)AlF (an alternate, although less likely identified would be the J= 9-8 transition of NaF). The derived Mg-24:Mg-25:Mg-26 isotopic abundance ratios (78 : 11+/- 1 : 11 +/-1) are consistent with the solar system values (79.0:10.0:11.0), following Anders & Grevesse 1989). According to new calculations of evolutionary models of 3 solar mass and 5 solar mass asymptotic giant branch (AGB) stars, these ratios and the previously measured N, O and Si isotopic ratios imply that the central star had an initial mass 3 solar mass (less than or equal to M(sub *, ini) less than 5 solar mass and has already experienced many 3rd dredge-up events. From this, it can be predicted that the Al-26/Al-27 isotopics ratio lies between 0.01 and 0.08; in fact, the value derived in the case that U234433 arises from (26)AlF is Al-26/Al-27 = 0.04. The identification of the (25)MgNC and (26)MgNC lines was made possible by ab-initio quantum mechanical calculations of the molecule geometrical structure. It was confirmed through millimeter-wave laboratory measurements. The quantum mechanical calculations are briefly described and the laboratory results presented in some detail. The rotation constants B, D, H and the spin-rotation constant gamma of (25)MgNC and (26)MgNC are determined from a fit of laboratory and astronomical data.
We have developed computational tools, based on particle codes, for molecular dynamics (MD) simulation of carbon nanotubes (CNT) in aqueous environments. The interaction of CNTs with water is envisioned as a prototype for the design of engineering nano-devices, such as artificial sterocillia and molecular biosensors. Large scale simulations involving thousands of water molecules are possible due to our efficient parallel MD code that takes long range electrostatic interactions into account. Since CNTs can be considered as rolled up sheets of graphite, we expect the CNT-water interaction to be similar to the interaction of graphite with water. However, there are fundamental differences between considering graphite and CNTs, since the curvature of CNTs affects their chemical activity and also since capillary effects play an important role for both dynamic and static behaviour of materials inside CNTs. In recent studies Gordillo and Marti described the hydrogen bond structure as well as time dependent properties of water confined in CNTs. We are presenting results from the development of force fields describing the interaction of CNTs and water based on ab-initio quantum mechanical calculations. Furthermore, our results include both water flows external to CNTs and the behaviour of water nanodroplets inside heated CNTs. In the first case (external flows) the hydrophobic behaviour of CNTs is quantified and we analyze structural properties of water in the vicinity of CNTs with diagnostics such as hydrogen bond distribution, water dipole orientation and radial distribution functions. The presence of water leads to attractive forces between CNTs as a result of their hydrophobicity. Through extensive simulations we quantify these attractive forces in terms of the number and separation of the CNT. Results of our simulations involving arrays of CNTs indicate that these exhibit a hydrophobic behaviour that leads to self-organising structures capable of trapping water clusters. In the second case (internal flows) we study the behaviour of water droplets confined inside CNTs. Constant temperature simulations allow us to capture structural properties such as the contact angles and density profiles of the equilibrated drops. By heating and subsequently cooling of the CNT, we are able to measure the evaporation and the condensation rate of the entrapped water.
Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity solutions of environments containing combined continuum/rarified flow. The use of DSMC, a stochastic, particle-based method, is necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally very computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the used thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. While the transport models governing CFD and DSMC simulations are starkly different, they can effectively be linked via their use of reduced Chapman-Enskog collision integrals. In CFD, these integrals are typically stored as fitted polynomial expressions and used to directly compute gas transport properties via mixing rules or the full Chapman-Enskog formulation. In DSMC, they can be used to derive the collision parameters needed for the phenomenological collision cross-section models that govern particle interactions, via a Nelder-Mead optimization scheme. The goal of this work is to provide a unified DSMC transport database encompassing the vast majority of known gas species encountered during atmospheric entry, on Earth or any other Solar body. This goal is largely possible due to recently performed ab-initio quantum chemistry calculations. Combined with other high-fidelity literature sources, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 17000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute Variable Soft Sphere (VSS) collision model parameters for DSMC, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database will be able to be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below. Equilibrium mixture transport properties for a 19-species Titan atmospheric model [4] are computed using both fitted VSS parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches is less than 5% for the entire temperature range.
The Wigner-Eisenbud R matrix method has been combined with the frame transformation theory to study electron scattering from molecular systems. The R matrix, calculated at the boundary point of the molecular core radius, has been transformed to the space frame in order to continue the solution of the scattering equations in the outer region where rotational motion of the nuclei is taken into account. This procedure has been applied to a model calculation of thermal energy electron scattering from CO.
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.
While hydrocarbon combustion and gas-surface interactions have been simulated for a variety of applications in continuum solvers, such processes have been rarely investigated in Direct Simulation Monte Carlo (DSMC). Thus Variable Hard/Soft Sphere (VHS/VSS)collision parameters are not available for most molecules of interest in the decomposition of methane, a significant pyrolysis gas product and potential atmospheric species for Titan re-entries. As these properties are essential for accurate modeling of gas transport and thermochemical effects, a method must be devised to compute them for use in DSMC. In this study, the methodology outlined in Stephani et al is utilized to compute the VHS/VSS parameters by fitting the Ω(1,1)and Ω(2,2) collision integrals [1]. These collision integrals are typically computed from intermolecular potentials and are used in CFD simulations to calculate transport properties such as viscosity, thermal conductivity, and diffusion. This approach ensures consistency between micro-scale collisions and macro-scale transport properties. To produce the necessary collision integrals needed for the fitting routine, we utilize the method described by Laricchiuta et al, where collision integrals are computed based on the constituent particles’ polarizability and effective electron number [2].Preliminary computations of pure-species Ω(1,1)and Ω(2,2)collision integrals produced from the coupled approach of Stephani et al and Laricchiuta et al for methane show good agreement with those computed via available ab-initio Potential Energy Surface Computations, as seen in Figure 1. Above 300K, no more than 12% deviation is observed between the two results over the temperature range of the fitted PES data [3].Current results of this study include full collision-specificVSS transport properties for a complex-species mixture involving major byproducts of pure methane decomposition. Targeted final results will include a comprehensive list of VSS parameters relevant toDSMC simulations of Titan atmospheric entry, based on the reduced kinetic model of Savajano et al [4].
In this work we present a molecular level study of N2+N collisions, focusing on excitation of internal energy modes and non-equilibrium dissociation. The computation technique used here is the direct molecular simulation (DMS) method and the molecular interactions have been modeled using an ab−initio potential energy surface (PES) developed at NASA's Ames Research Center. We carried out vibrational excitation calculations between 5000K and 30000K and found that the characteristic vibrational excitation time for the N + N2 process was an order of magnitude lower than that predicted by the Millikan and White correlation. It is observed that during vibrational excitation the high energy tail of the vibrational energy distribution gets over populated first and the lower energy levels get populated as the system evolves. It is found that the non-equilibrium dissociation rate coefficients for the N + N2 process are larger than those for the N2 + N2 process. This is attributed to the non-equilibrium vibrational energy distributions for the N + N2 process being less depleted than that for the N2 +N2 process. For an isothermal simulation we find that the probability of dissociation goes as 1/T(sub tr) for molecules with internal energy (epsilon(sub int)) less than approximately 9.9eV, while for molecules with epsilon (sub int) greater than 9.9eV the dissociation probability was weakly dependent on translational temperature of the system. We compared non-equilibrium dissociation rate coefficients and characteristic vibrational excitation times obtained by using the ab-initio PES developed at NASA's Ames Research Center to those obtained by using an ab-initio PES developed at the University of Minnesota. Good agreement was found between the macroscopic properties and molecular level description of the system obtained by using the two PESs.