SEARCH · Engineering Papers
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.
A model of fluid drop behavior in clusters has been developed including the interactions induced by the drop proximity. The model is based upon the global conservation equations for the cluster coupled to isolated fluid drop equations previously developed.
Azirinyl cation is an aromatic cyclic molecule that is isoelectronic with cyclopropenylidene, c-C3H2, and c-C3H3+. Cyclopropenylidene has been shown to be ubiquitous, existing in many different astrophysical environments. Given the similar chemistry between C and N, and the relative abundances between C and N in astrophysical environments, it is expected that there should be aromatic ringed molecules that incorporate N in the ring, but as yet, no such molecule has been identified. To address this issue, the present study uses high levels of electronic structure theory to compute a highly accurate quartic force field (QFF) for azirinyl cation and its two lowest lying isomers, the cyanomethyl and isocyanomethyl cations. The theoretical approach uses the singles and doubles coupled-cluster method that includes a perturbative correction for connected triple excitations, CCSD(T), together with extrapolation to the one-particle basis set limit and corrections for scalar relativity and core-correlation. The QFF is then used in a second order vibrational perturbation theory analysis (VPT2) to compute the fundamental vibrational frequencies and rovibrational spectroscopic constants for all three C2NH2+ isomers. The reliability of the VPT2 vibrational frequencies is tested by comparison to vibrational configuration interaction (VCI) calculations and excellent agreement is found between the two approaches. Fundamental vibrational frequencies and rovibrational spectroscopic constants for all singly substituted 13C, 15N, and D isotopologues are also reported. It is expected that the highly accurate spectroscopic data reported herein will be useful in the identification of these cations in high resolution experimental or astronomical observations.
Identifying new nitrogenated hydrocarbon molecules in the interstellar medium (ISM) is challenging due to the lack of comprehensive spectroscopic data from experiments. In this computational work, we focus on investigating the structures, relative energies, spectroscopic constants, and energy landscape of the cyclic isomers of diazirine (c-CH2N2) using ab initio quantum chemical methods. Density functional theory (DFT) methods and coupled cluster theory with singles and doubles including perturbative triples [CCSD(T)] and CCSD(T) with the explicitly correlated F12 b correction [CCSD(T)-F12b] were employed for this purpose along with large correlation consistent cc-pVTZ, cc-pVQZ, and cc-pV5Z basis sets. Harmonic vibrational frequencies, infrared vibrational intensities, rotational constants, and dipole moments are reported. An harmonic vibrational fundamentals along with centrifugal distortion constants, and vibration-rotation interaction constants are also reported for all the cyclic isomers. The energies computed with the CCSD(T) and CCSD(T)-F12b methods were extrapolated to the one-particle complete basis set (CBS) limit following a three-point formula. At the CCSD(T)-F12b/CBS level of theory the 3,3H-diazirine (c-CH2N2) is the lowest energy cyclic isomer followed by 1,3H-diazirine, (E)-1,2H-Diazirine, and (Z)-1,2H-diazirine which are 20.1, 47.8, and 51.3kcal mol-1above the 3,3H-Diazirine, respectively. Accurate structures and spectroscopic constants that are reported here could be useful for future identification of these cyclic nitrogenated organic molecules in the interstellar medium or circumstellar disks.
Due to coherence time limitations, reducing the resources required to run quantum algorithms and simulate physical systems on a quantum computer is crucial. With regards to Hamiltonian simulation, a significant effort has focused on building efficient algorithms using various factorizations and truncations, typically derived from the Hamiltonian alone. We introduce a new paradigm for improving Hamiltonian simulation and reducing the cost of ground state problems based on ideas recently developed for classical chemistry simulations. The key idea is that one can find efficient ways to reduce resources needed by quantum algorithms by making use of two key pieces of information: the Hamiltonian operator and an approximate ground state wavefunction. We refer to our algorithm as the self-consistent quantum iteratively sparsified Hamiltonian (SQuISH). By performing our scheme iteratively, one can drive SQuISH to create an accurate wavefunction using a truncated, resource-efficient Hamiltonian. By utilizing this more compact Hamiltonian, our algorithm provides an approach to reduce the gate complexity of ground state calculations on quantum hardware. As proof of principle, we implement SQuISH using configuration interaction for small molecules and coupled cluster for larger systems. Through our combination of approaches, we demonstrate how it performs on a range of systems, the largest of which would require more than 200 qubits to run on quantum hardware.
In a high temperature combustion situation such as a gas turbine or rocket engine, combustion products include water vapor. This water vapor reacts with oxides in the hot-gas-path to form volatile hydroxides and oxyhydroxides. The oxides may be structural oxide components or protective oxides on alloys or non-oxide ceramics. In either case, this enhanced volatilization leads to degradation of components. Reliable thermochemical data for many metal hydroxides and oxyhydroxides has not been available. We review a database put together at NASA over many years based on experimental and ab initio computations for thermochemical data of Al, Si, Ti, Cr, Mn, Y, Zr, La, Ta, Gd, and Yb hydroxides and oxyhydroxides. Experimental methods for measuring thermodynamic quantities of metal hydroxides and oxyhydroxides are limited since the required high temperatures and high pressures preclude many conventional methods. Nonetheless, transpiration has proven a valuable tool in our laboratories for the Si-O-H, Cr-O-H, and Ti-O-H systems. Data are analyzed with second and third law methods. These data can be leveraged with quantum chemistry to provide reliable free energy functions for the third law calculation and comparison to the derived enthalpies of formation. Both wave function and DFT methods can be utilized for deriving the thermodynamic functions for the various metal hydroxides and oxyhydroxides. We have found that, in general, reasonable molecular geometries and spectroscopic data can be derived with a relatively low level of theory, i.e. DFT with the B3LYP functional. Accurate calculations of the enthalpy of formation require higher levels of theory, such as the Coupled Cluster, Singles, Doubles, and perturbative Triples CCSD(T) approach with complete basis sets (CBS). Also, corrections for hindered rotations and anharmonicity are necessary. For some symmetries, the anharmonic correction has proven to be a challenge and we have found that this is only treated properly by certain codes. The database lists the fHo(298), So(298), and Cp in the form Cp = A + BT + CT2 + D/T + E/T2, where A, B, C, D, and E are the fitting constants. Thus, it can be easily used in any of the common computational thermochemistry codes. To date, the database contains 51 gaseous species of the formula MOx(OH)y.
The ability to perform ab initio molecular dynamics simulations using potential energy surfaces provided by quantum computers would open the door to virtually exact dynamics for a variety of chemical and biochemical systems, with impacts on catalysis and biophysics. Nonetheless, performing molecular dynamics on surfaces produced by quantum hardware has been hampered by the noisy energies typically produced by quantum computers and challenges associated with computing gradients and scaling to large systems interest. A recent set of advances in machine learning, known as transfer learning, provides a new path forward for molecular dynamics simulations on quantum hardware. Transfer learning offers a workaround, where one first trains models on larger, less accurate classical datasets and then refines them on smaller, more accurate quantum datasets. We explore this approach by training machine learning models to predict a molecule's potential energy based on its geometric structure using Behler-Parrinello neural networks. When successfully trained, the model enables energy gradient predictions necessary for dynamic simulations. To reduce the quantum resources needed, the model is initially trained with data derived from classical density functional theory and subsequently refined with a smaller dataset obtained from a variational quantum eigensolver optimization of the unitary coupled cluster ansatz. We show that this approach significantly reduces the size of the needed quantum training dataset while capturing the high accuracies needed within quantum chemistry simulations. The success of this two-step training method opens more opportunities to apply machine learning models on quantum data, a significant stride towards efficient quantum-classical hybrid computational models.
Matrix-Holzer analyses for predicting free vibration modes of clustered launch vehicle configurations including Titan IIIC and uncoupled Saturn I cases
We present performance results for three distributed computing environments using the three simulated CFD applications in the NAS Parallel Benchmark suite. These environments are the DCF cluster, the LACE cluster, and an Intel iPSC/860 machine. The DCF is a prototypic cluster of loosely coupled SGI R3000 machines connected by Ethernet. The LACE cluster is a tightly coupled cluster of 32 IBM RS6000/560 machines connected by Ethernet as well as by either FDDI or an IBM Allnode switch. Results of several parallel algorithms for the three simulated applications are presented and analyzed based on the interplay between the communication requirements of an algorithm and the characteristics of the communication network of a distributed system.
The performance of a loosely and tightly-coupled workstation cluster is compared against a conventional vector supercomputer for the solution the Reynolds- averaged Navier-Stokes equations. The application geometries include a transonic airfoil, a tiltrotor wing/fuselage, and a wing/body/empennage/nacelle transport. Decomposition is of the manager-worker type, with solution of one grid zone per worker process coupled using the PVM message passing library. Task allocation is determined by grid size and processor speed, subject to available memory penalties. Each fluid zone is computed using an implicit diagonal scheme in an overset mesh framework, while relative body motion is accomplished using an additional worker process to re-establish grid communication.
We propose a numerical study designed to interpret the origin and evolution of galaxy properties revealed by space- and ground-based imaging and spectroscopical surveys. Our aim is to unravel the physical processes responsible for the development of different galaxy morphologies and for the establishment of scaling laws such as the Tully-Fisher relation for spirals and the Fundamental Plane of ellipticals. In particular, we plan to address the following major topics: (1) The morphology and observability of protogalaxies, and in particular the relationship between primordial galaxies and the z approximately 3 'Ly-break' systems identified in the Hubble Deep Field and in ground-based searches; (2) The origin of the disk and spheroidal components in galaxies, the timing and mode of their assembly, the corresponding evolution in galaxy morphologies and its sensitivity to cosmological parameters; (3) The origin and redshift evolution of the scaling laws that link the mass, luminosity size, stellar content, and metal abundances of galaxies of different morphological types. This investigation will use state-of-the-art N-body/gasdynamical codes to provide a spatially resolved description of the galaxy formation process in hierarchically clustering universes. Coupled with population synthesis techniques. our models can be used to provide synthetic 'observations' that can be compared directly with observations of galaxies both nearby and at cosmologically significant distances. This study will thus provide insight into the nature of protogalaxies and into the formation process of galaxies like our own Milky Way. It will also help us to assess the cosmological significance of these observations within the context of hierarchical theories of galaxy formation and will supply a theoretical context within which current and future observations can be interpreted.
We implemented the ADAPT-VQE algorithm into our recent classical sparse wavefunction circuit simulator to demonstrate that classical resources can (1) efficiently find a physically motivated compact wavefunction ansatz for further refinement on near-term quantum hardware and (2) benchmark expected results of VQE-based algorithms once the quantum hardware is available to study large-scale applications. In particular, we study the role of the ADAPT-VQE operator pool, molecular basis set selection, and variations such as TETRIS-ADAPT-VQE on the performance of our classical circuit simulator. This work demonstrates the promise of using classical resources to generate highly accurate wavefunctions that can be prepared on quantum hardware to initiate other quantum algorithms such as phase estimation. Our approach harnesses the power of high-performance computing resources with the more limited available quantum computers to map a path toward quantum advantage for electronic structure calculations in chemistry and materials science.
Global IMAGE/EUV observations have revealed complex changes in plasmaspheric structures as the plasmasphere responds to geomagnetic activity while remaining under varying degrees of influence by co-rotation, depending on the radial distance. The complex plasmaspheric dynamics, with different scales of variability, is clearly far from being well understood. There is now renewed interest in the plasmasphere due to its apparent connections with the development of the ring current and radiation belt, and loss of ionospheric plasmas. Early in the mission, the Cluster spacecraft only crossed the plasmapause (L - 4) occasionally and made measurements of the outer plasmasphere and plasmaspheric drainage plumes. The study by Darrouzet et al. [2006] provided detailed analyses of in situ Cluster observations and IMAGE EUV observations of three plasmaspheric plumes detected in April-June, 2002. Within the next couple of years, Cluster orbit will change, causing perigee to migrate to lower altitudes, and thus providing excellent opportunities to obtain more detailed measurements of the plasmasphere. In this paper, we report our analyses of the earlier Cluster-IMAGE events by incorporating the different perspectives provided by the IMAGE Radio Plasma Imager (RPI) observations. We will discuss our new understanding of the structure and dynamics of the Cluster-IMAGE events.
Global EUV imaging observations by the IMAGE satellite have revealed spectacularly complex changes in plasmaspheric structures as the plasmaspheric plasmas respond to geomagnetic activity while remaining under varying degrees of influence by co-rotation, depending on the radial distance. This complex plasmaspheric plasma dynamics, with its numerous scales of variability, is clearly far from being well understood. There is now renewed interest in studying the plasmasphere due to its apparent connections with the development of the ring current and radiation belt and loss of ionospheric plasmas. Earlier in the mission, the Cluster spacecraft only crossed the plasmapause (L - 4) occasionally and made measurements of the cold plasma in the plasmasphere and plasmaspheric drainage plumes. The study by Darrouzet et al. [2006], for example, provided detailed analyses of in situ Cluster observations and IMAGE EUV imaging observations of three plasmaspheric plumes detected in April-June, 2002. Within the next couple of years, Cluster orbit will change, causing perigee to migrate to lower altitudes, and thus providing excellent opportunities to obtain more detailed measurements of the plasmasphere. In this paper, we report our analyses of the earlier Cluster-IMAGE events by incorporating the different perspectives provided by the IMAGE Radio Plasma Imager (RPI) observations through the plasmasphere. We will describe our new understanding of the Cluster-IMAGE events and their implications on plasmaspheric dynamics.
Our analysis concerns macro and meso-scale aspects of coupling between the IMF and the magnetosphere-ionosphere system, as opposed to the microphysics of determining how electron gyrotropy is broken and merging actually occurs. We correlate observed behaviors at Cluster and at Polar with temporal variations in other regions, such as in the ionosphere as measured by SuperDARN. Addressing problems with simultaneous observations from diverse locations properly constrains our interpretations.
We have obtained deep charge coupled device (CCD)V and I images of a number of fields near the center of the Pleiades open cluster. We have also obtained imaging data for Praesepe, a very similar cluster in terms of distance and richness but nearly 10 times older than the Pleiades. Because brown dwarfs are predicted to become much fainter and cooler between Pleiades and Praesepe ages, this provides a powerful differential technique for placing constraints on the brown dwarf population in open clusters. Combined with our previously reported observations, we now have data for about 0.4 sq deg in the Pleiades, corresponding roughly to 5% of the area of that cluster. We have searched the new CCD frames for additional Pleiades brown dwarf candidates. Two possible candidates are present, the faintest of which has V approximately equal to 22.5, (V-I)(sub K) approximately equal to 4.6. Because we do not have proper motion data and the colors of these objects are not redder than the reddest known field stars, it is possible that some or all of our candidates are somewhat higher mass field stars rather than Pleiades-age brown dwarfs. Even if all six of the proposed brown dwarf candidates in our 0.4 sq deg field are Pleiades members, the relatively small number found suggests that low mass stars or brown dwarfs do not contribute significantly to the total mass of the cluster.
The geometries of the B(n)H(+) clusters for n=1-13 have been optimized at the B3LYP level of theory. Excluding B2H(+) all clusters are low-spin coupled. The structures appear to be related to the bare B(+)(n) clusters in two ways: 1) an H atom is added to the B(+)(n) cluster or 2) a BH is added to the B(+)(n-1) cluster. The B(+)(n)-H binding energies are computed at the B3LYP level of theory, and calibrated n using the CCSD(T) level of theory. The computed results tend to fall between the experimental best estimates and the experimental lower bounds. The B(+)(n)-H binding energies show an inverse correspondence with the B+(n-1)-B values) while the B(+)(n-1)-B and B(+)(n-1)-BH results tend to parallel each other.
We report the identification of five short-period variables near the center of the metal-rich globular cluster M71. Our observations consist of multiepoch VI charge coupled device (CCD) images centered on the cluster and covering a 6.3 min x 6.3 min field. Four of these variables are contact eclipsing binaries with periods between 0.35 and 0.41 days; one is a detached or semidetached eclipsing binary with a period of 0.56 days. Two of the variables were first identified as possible eclipsing binaries in an earlier survey by Hodder et al. (1992). We have used a variety of arguments to conclude that all five binary stars are probable members of M71, a result that is consistent with the low number (0.15) of short-period field binaries expected along this line of sight. Based on a simple model of how contact binaries evolve from initially detached binaries, we have determined a lower limit of 1.3% on the frequency of primordial binaries in M71 with initial orbital periods in the range 2.5 - 5 days. This implies that the overall primordial binary frequency, f, is 22(sup +26)(sub -12)% assuming df/d log P = const ( the 'flat' distribution), or f = 57(sup +15)(sub -8)% for df/d log P = 0.032 log P + const as observed for G-dwarf binaries in the solar neighborhood (the 'sloped' distribution). Both estimates of f correspond to binaries with initial periods shorter than 800 yr since any longer-period binaries would have been disrupted over the lifetime of the cluster. Our short-period binary frequency is in excellent agreement with the observed frequency of red-giant binaries observed in globulars if we adopt the flat distribution. For the sloped distribution, our results significantly overestimate the number of red-giant binaries. All of the short-period M71 binaries lie within 1 mag of the luminosity of the cluster turnoff in the color-magnitude diagram despite the fact we should have easily detected similar eclipsing binaries 2 - 2.5 mag fainter than this. We discuss the implications of this on our estimates of the binary frequency in M71 and on the formation of blue stragglers.
The age of the universe should be calculable by independent methods with similar results. Previous calculations using nucleochronometers, globular clusters and dynamical measurements coupled with Friedmann models and nucleosynthesis constraints have given different values of the age. A consistent age is reported, whose implications for the constituent mass density are very interesting and are affected by the existence of a third neutrino flavor, and by allowing the possibility that neutrinos may have a non-zero rest mass.