Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic and molecular physics”

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 181 records · Page 10

Experimental Resonance Enhanced Multiphoton Ionization (REMPI) studies of small molecules

Resonance enhanced multiphoton ionization (REMPI) utilizes tunable dye lasers to ionize an atom or molecule by first preparing an excited state by multiphoton absorption and then ionizing that state before it can decay. This process is highly selective with respect to both the initial and resonant intermediate states of the target, and it can be extremely sensitive. In addition, the products of the REMPI process can be detected as needed by analyzing the resulting electrons, ions, fluorescence, or by additional REMPI. This points to a number of exciting opportunities for both basic and applied science. On the applied side, REMPI has great potential as an ultrasensitive, highly selective detector for trace, reactive, or transient species. On the basic side, REMPI affords an unprecedented means of exploring excited state physics and chemistry at the quantum-state-specific level. An overview of current studies of excited molecular states is given to illustrate the principles and prospects of REMPI.

Dehmer, J. L.↗

Contamination on LDEF: Sources, distribution, and history

An introduction to contamination effects observed on the Long Duration Exposure Facility (LDEF) is presented. The activities reported are part of Boeing's obligation to the LDEF Materials Special Investigation Group. The contamination films and particles had minimal influence on the thermal performance of the LDEF. Some specific areas did have large changes in optical properties. Films also interfered with recession rate determination by reacting with the oxygen or physically shielding underlying material. Generally, contaminant films lessen the measured recession rate relative to 'clean' surfaces. On orbit generation of particles may be an issue for sensitive optics. Deposition on lenses may lead to artifacts on photographic images or cause sensors to respond inappropriately. Particles in the line of sight of sensors can cause stray light to be scattered into sensors. Particles also represent a hazard for mechanisms in that they can physically block and/or increase friction or wear on moving surfaces. LDEF carried a rather complex mixture of samples and support hardware into orbit. The experiments were assembled under a variety of conditions and time constraints and stored for up to five years before launch. The structure itself was so large that it could not be baked after the interior was painted with chemglaze Z-306 polyurethane based black paint. Any analysis of the effects of molecular and particulate contamination must account for a complex array of sources, wide variation in processes over time, and extreme variation in environment from ground to launch to flight. Surface conditions at certain locations on LDEF were established by outgassing of molecular species from particular materials onto adjacent surfaces, followed by alteration of those species due to exposure to atomic oxygen and/or solar radiation.

Pippin, Gary↗

Interfacing Relativistic and Nonrelativistic Methods: A Systematic Sequence of Approximations

A systematic sequence of approximations for the introduction of relativistic effects into nonrelativistic molecular finite-basis set calculations is described. The theoretical basis for the approximations is the normalized elimination of the small component (ESC) within the matrix representation of the modified Dirac equation. The key features of the normalized method are the retention of the relativistic metric and the ability to define a single matrix U relating the pseudo-large and large component coefficient matrices. This matrix is used to define a modified set of one- and two-electron integrals which have the same appearance as the integrals of the Breit-Pauli Hamiltonian. The first approximation fixes the ratios of the large and pseudo-large components to their atomic values, producing an expansion in atomic 4-spinors. The second approximation defines a local fine-structure constant on each atomic centre, which has the physical value for centres considered to be relativistic and zero for nonrelativistic centres. In the latter case, the 4-spinors are the positive-energy kinetic al ly-balanced solutions of the Levy-Leblond equation, and the integrals involving pseudo-large component basis functions on these centres, are set to zero. Some results are presented for test systems to illustrate the various approximations.

Dyall, Ken↗

A Solar System Perspective on Laboratory Astrophysics

Planetary science deals with a wide variety of natural materials in a wide variety of environments. These materials include metals, minerals, ices, gases, plasmas, and organic chemicals. In addition, the newly defined discipline of astrobiology introduces biological materials to planetary science. The environments range from the interiors of planets with megapascal pressures to planetary magnetospheres, encompassing planetary mantles, surfaces, atmospheres, and ionospheres. The interplanetary environment includes magnetic and electrical fields, plasma, and dust. In order to understand planetary processes over these vast ranges, the properties of materials must be known, and most of the necessary information comes from the laboratory. Observations of the bodies and materials in the Solar System are accomplished over the full range of the electromagnetic spectrum by remote sensing from Earth or spacecraft. Comets exemplify this; molecular and atomic identifications are made from the hard ultraviolet to radio wavelengths, while X-rays are emitted as comets interact with the solar wind. Gamma rays from the surfaces of the Moon and asteroids are diagnostic of the mineral and ice content of those bodies; eventually, gamma rays will also be observed by probes to comets. A number of planetary materials are available in the laboratory for extensive Study: rocks from the Moon, Mars, several asteroids, as well as dust from comets (and perhaps the Kuiper Belt) are closely studied at every level, including atomic (isotopic). Even pre-solar interstellar grains isolated from meteorites are scrutinized for composition and crystalline structure. Beyond the materials themselves, various agents and processes have altered them over the 4.6-Gy age of the Solar System. Solar radiation, solar wind particles, trapped magnetospheric particles, cosmic rays, and micrometeoroid impacts have produced chemical, physical, and morphological changes in the atmospheres and on the surfaces of all planetary bodies. These processes are not well understood, so studies in a laboratory setting are especially needed.

Cruikshank, Dale P.↗

J-OCTA Applications to in Space Manufacturing and Nanofluidic Technologies

We consider examples of J-OCTA software application to physics-based analysis of in-space manufacturing technologies. Firstly, we consider an atomistic model of thermal welding at the polymer-polymer interface for applications to 3D manufacturing by following diffusion of polymer chains. We show that each component of the polymer blend has its own characteristic time of diffusion and that both strain–stress and shear viscosity curves agree with experimental data. Next, nonlinear shrinkage of the metal part during manufacturing by bound metal deposition under microgravity, is considered using a multi-scale physics-based approach that spans from atomistic dynamics to full-part shrinkage. At the smallest scale we use atomistic dynamics to simulated sintering of a few Ti6Al4V nanoparticles, study sintering mechanisms and estimate important parameters including grain boundary width and diffusion rates of the atoms at the surface that are used at the other scales of modelling. Finally, we analyse selective ionic conduction through an artificial nanopore in a graphene sheet motivated by applications to desalination and energy harvesting. We build distributions of the number and orientation of water molecules in the ion’s hydration shells providing novel insight into the ion-water-carbon interaction at the nanopore.

molecular dynamics↗

Physical conditions in photodissociation regions: Application to galactic nuclei

Infrared and sub-millimeter observations are used in a simple procedure to determine average physical properties of the neutral interstellar medium in Galactic photodissociation regions as well as in ensembles of clouds which exist in the nuclei of luminous infrared galaxies. The relevant observations include the Infrared Astronomy Satellite (IRAS) infrared continuum measurements, infrared spectroscopy of the fine-structure lines of SiII 35 microns, OI 63 microns, and CII 158 microns, and the 2.6 mm CO (J=1-0) rotational transition. The diagnostic capabilities of the OI 145 microns line is also addressed. Researchers attribute these emission lines as well as the continuum to the atomic/molecular photodissociation region on the surfaces of molecular clouds which are illuminated by strong ultraviolet fields. They use the theoretical photodissociation region models of Tielens and Hollenbach (1985, Ap. J., 291, 722) to construct simple diagrams which utilize line ratios and line to continuum ratios to determine the average gas density n, the average incident far-ultraviolet flux G sub o, and the temperature of the atomic gas T.

Wolfire, M. G.↗

The Materials Chemistry of Atomic Oxygen with Applications to Anisotropic Etching of Submicron Structures in Microelectronics and the Surface Chemistry Engineering of Porous Solids

Neutral atomic oxygen is the most abundant component of the ionospheric plasma in the low Earth orbit environment (LEO; 200 to 700 kilometers altitude) and can produce significant degradation of some spacecraft materials. In order to produce a more complete understanding of the materials chemistry of atomic oxygen, the chemistry and physics of O-atom interactions with materials were determined in three radically different environments: (1) The Space Shuttle cargo bay in low Earth orbit (the EOIM-3 space flight experiment), (2) a high-velocity neutral atom beam system (HVAB) at Los Alamos National Laboratory (LANL), and (3) a microwave-plasma flowing-discharge system at JSC. The Space Shuttle and the high velocity atom beam systems produce atom-surface collision energies ranging from 0.1 to 7 eV (hyperthermal atoms) under high-vacuum conditions, while the flowing discharge system produces a 0.065 eV surface collision energy at a total pressure of 2 Torr. Data obtained in the three different O-atom environments referred to above show that the rate of O-atom reaction with polymeric materials is strongly dependent on atom kinetic energy, obeying a reactive scattering law which suggests that atom kinetic energy is directly available for overcoming activation barriers in the reaction. General relationships between polymer reactivity with O atoms and polymer composition and molecular structure have been determined. In addition, vacuum ultraviolet photochemical effects have been shown to dominate the reaction of O atoms with fluorocarbon polymers. Finally, studies of the materials chemistry of O atoms have produced results which may be of interest to technologists outside the aerospace industry. Atomic oxygen 'spin-off' or 'dual use' technologies in the areas of anisotropic etching in microelectronic materials and device processing, as well as surface chemistry engineering of porous solid materials are described.

Koontz, Steve L.↗

Angular motion of a PAH molecule in interstellar environment

Polycyclic aromatic hydrocarbon (PAH) molecules have recently been proposed as an important and hitherto undetected component of the Interstellar Medium (ISM). The theory was based on an explanation of the Unidentified IR Emission Bands by Leger et al. It has already led to a verified prediction on extended galactic and extragalactic emissions measured by IRAS, or by a recent balloon borne experiment. The physics that rules the motion of such molecules in the ISM was studied, taking into account their coupling with the ambient gas, the radiation field (absorption and emission) and the static magnetic field. This is important for many implications of the PAH theory such as the radio emission by these molecules or the expected polarization of their IR emission. A reflection nebulae is considered where the situation is rather well known. Every day life of a mean PAH molecule in such a region is as follows: every 3 hrs a UV photon is absorbed heating the molecule to a thousand degs; the temperature decay due to cooling by IR emission follows then within a few seconds. A collision with a molecule of gas occurs typically once a week, while an H atom is ejected or captured at the same rate. A typical cooling cycle after a heat impulse is given. The PAH molecules studied as representative of the family has typically 50 atoms, a radius of 4.5 A, is circular and has a molecular mass of M = 300; its permanent dipole moment is 3 Debye.

Rouan, D.↗

A Comparison of EAST Shock-Tube Radiation Measurements with a New Air Radiation Model

This paper presents a comparison between the recent EAST shock tube radiation measurements (Grinstead et al., AIAA 2008-1244) and the HARA radiation model. The equilibrium and nonequilibrium radiation measurements are studied for conditions relevant to lunar-return shock-layers; specifically shock velocities ranging from 9 to 11 kilometers per second at initial pressures of 0.1 and 0.3 Torr. The simulated shock-tube flow is assumed one-dimensional and is calculated using the LAURA code, while a detailed nonequilibrium radiation prediction is obtained in an uncoupled manner from the HARA code. The measured and predicted intensities are separated into several spectral ranges to isolate significant spectral features, mainly strong atomic line multiplets. The equations and physical data required for the prediction of these strong atomic lines are reviewed and their uncertainties identified. The 700-1020 nm wavelength range, which accounts for roughly 30% of the radiative flux to a peak-heating lunar return shock-layer, is studied in detail and the measurements and predictions are shown to agree within 15% in equilibrium. The plus or minus 1.5% uncertainty on the measured shock velocity is shown to cause up to a plus or minus 30% difference in the predicted radiation. This band of predictions contains the measured values in almost all cases. For the highly nonequilibrium 0.1 Torr cases, the nonequilibrium radiation peaks are under-predicted by about half. This under-prediction is considered acceptable when compared to the order-of-magnitude over-prediction obtained using a Boltzmann population of electronic states. The reasonable comparison in the nonequilibrium regions provides validation for both the non-Boltzmann modeling in HARA and the thermochemical nonequilibrium modeling in LAURA. The N2 (+)(1-) and N2(2+) molecular band systems are studied in the 290 480 nm wavelength range for both equilibrium and nonequilibrium regimes. The non-Boltzmann rate models for these systems, which have significant uncertainties, are tuned to improve the comparison with measurements.

Johnston, Christopher O.↗

The Diversity of Refractory Organic Material in Comets

Organic matter exists in comets (most notably in 81P/Wild 2 [Stardust], 67P/Churyomov- Gerasimenko (67P/C-G) [Rosetta], chondritic porous IDPS, and UCAMMs) and in primitive carbonaceous chondrites that likely retain some chemistry that reflects an origin in the prenatal cold molecular cloud (Alexander+2017). Heavy isotopic enrichments, 15N/14N and possibly D/H, signify preserved molecular cloud organics. In the cold outer disk, if grains are lofted above the disk mid-plane then organics likely experience significant UV processing (Ciesla+2012). In remote sensing of comet comae, organics in the dust are considered refractory or semirefractory. Semi-refractory organics have limited comae lifetimes and produce distributed sources of molecules (H2CO and CO). Rosetta's close passes of 67P/C-G's nucleus (10-15 km) reveals a distributed source of glycine, methyl amine and ethylamine (Altwegg+2016). Cometary samples and primitive meteorites have two types of organic matter: (1) acid-insoluble organic matter (IOM), which is a macromolecular polymer with a mixture of aromatic and aliphatic moieties, and (2) labile, soluble organics, which includes the amino acids, such as glycine (Stardust, Elsila+2009). Meteoritic IOM is robust, withstanding experimental temperatures of 1200 K (Dobrica+2011, Cody+2008). Nanoglobules are a type of IOM; they have a distinct physical structure, but often share the same chemistry as the other IOM from the same meteorite. Moderate-sized PAHs (20 C-atoms) are detected in Stardust samples (Clemett+2010). Refractory organic IOM is ubiquitous yet has a great diversity of abundances between cometary samples. IOM is in primitive chondrites, 67P/C-G (Rosetta), 81P/Wild 2 (Stardust), 1P/Halley, 26P/Grigg-Skjellerup, UCAMMs, anhydrous IDPs, and in chondritic porous IDPs (CP IDPs) and larger cluster IDPs (e.g., Fray+2016, Fomenkova+94, Busemann+ 09, Dobrica+2011, Dobrica+2012). 81P's refractory organic matter is of two types (De Gregorio+2011): nanoglobules of highly aromatic refractory organic matter and polyaromatic carbonyl-containing organic matter, which is similar to IOM in primitive meteorites and IDPs. Fray+2017 estimate that 50% of carbon in 67P/C-G is in IOM. 67P/C-G's organics appear to lack the soluble organic matter, aliphatic carbon, amino acids, and PAHs (Fray+2016). Other notable aspects of the diversity in IOM in cometary samples are the ranges of atomic ratios of N/C, O/C, and H/C, and the range of isotopic enrichments of 15N/14N and D/H. Aqueous and thermal processing on asteroids changes the balance of soluble to insoluble organics, and may be important for diversifying the range of OM delivered to Earth.

Wooden, Diane↗

Vibrational and Electronic Energy Transfer and Dissociation of Diatomic Molecules by Electron Collisions

At high altitudes and velocities equal to or greater than the geosynchronous return velocity (10 kilometers per second), the shock layer of a hypersonic flight will be in thermochemical nonequilibrium and partially ionized. The amount of ionization is determined by the velocity. For a trans atmospheric flight of 10 kilometers per second and at an altitude of 80 kilometers, a maximum of 1% ionization is expected. At a velocity of 12 - 17 kilometer per second, such as a Mars return mission, up to 30% of the atoms and molecules in the flow field will be ionized. Under those circumstances, electrons play an important role in determining the internal states of atoms and molecules in the flow field and hence the amount of radiative heat load and the distance it takes for the flow field to re-establish equilibrium. Electron collisions provide an effective means of transferring energy even when the electron number density is as low as 1%. Because the mass of an electron is 12,760 times smaller than the reduced mass of N2, its average speed, and hence its average collision frequency, is more than 100 times larger. Even in the slightly ionized regime with only 1% electrons, the frequency of electron-molecule collisions is equal to or larger than that of molecule-molecule collisions, an important consideration in the low density part of the atmosphere. Three electron-molecule collision processes relevant to hypersonic flows will be considered: (1) vibrational excitation/de-excitation of a diatomic molecule by electron impact, (2) electronic excitation/de-excitation, and (3) dissociative recombination in electron-diatomic ion collisions. A review of available data, both theory and experiment, will be given. Particular attention will be paid to tailoring the molecular physics to the condition of hypersonic flows. For example, the high rotational temperatures in a hypersonic flow field means that most experimental data carried out under room temperatures are not applicable. Also, the average electron temperature is expected to be between 10,000 and 20,000 K. Thus only data for low energy electrons are relevant to the model.

Huo, Winifred M.↗

Photoevaporation of Clumps in Photodissociation Regions

We present the results of an investigation of the effects of Far Ultraviolet (FUV) radiation (6.0eV < hv < 13.6eV) from hot early type OB stars on clumps in star-forming molecular clouds. Clumps in FUV-illuminated regions (or photodissociation regions or PDRs) undergo external heating and photodissociation as they are exposed to the FUV field, resulting in a loss of cold, molecular lump mass as it is converted to warm atomic gas. The heating, if rapid, creates strong photoevaporative mass flows off the clump surfaces, and drives shocks into the clumps, compressing them to high densities. The clumps lose mass on relatively short timescales. The evolution of an individual clump is found to be sensitive to three dimensionless parameters: Nc0, the ratio of the initial column density of the clump to the column N(0) approx. 10(exp 21) cm(exp -2) of a warm FUV-heated surface region; upsilon, the ratio of the sound speed in the heated surface to that in the cold clump material: and t(FUV)t(c), the ratio of the "turn-on time" t(FUV) of the heating flux on a clump to its initial sound crossing-time t(c). The evolution also depends on whether a confining interclump medium exists, or whether the interclump region has negligible pressure, as is the case for turbulence-generated clumps. In this paper, we use spherical 1-D numerical hydrodynamic models as well as approximate analytical models to study the dependence of clump photoevaporation on the physical parameters of the clump, and to derive the dynamical evolution, mass loss rates and photoevaporative timescales of a clump for a variety of astrophysical situations. Turbulent clumps evolve so that their column densities are equal to a critical value determined by the local FUV field, and typically have short photo evaporation timescales, approx. 10(exp 4-5) years for a 1 M(solar mass) clump in a typical star-forming region (Nc0 = 10, upsilon = 10). Clumps with insufficient magnetic pressure support, and in strong FUV fields may be driven to collapse by the compressional effect of converging shock waves. We also estimate the rocket effect on photoevaporating clumps and find that it is significant only for the smallest clumps, with sizes much less than the extent of the PDR itself. Clumps that are confined by all interclump medium may either get completely photoevaporated, or may preserve a shielded core with a warm, dissociated, protective shell that, absorbs the incident FUV flux. We compare our results with observations of some well studied PDRs: the Orion Bar, M17SW NGC 2023 and the Rosette Nebula. The data are consistent with both interpretations of clump origin. turbulence and pressure confinement, with a slight indication for favouring the turbulent model for clumps over pressure-confined

Gorti, Uma↗

Data Base for CFD Validation

The flow behind the shock wave formed around objects which fly at hypervelocity behaves differently from that of a perfect gas. Molecules become vibrationally excited, dissociated, and ionized. The hot gas may emit or absorb radiation. When the atoms produced by dissociation reach the wall surface, chemical reactions, including recombination, may occur. The thermochemical phenomena of vibration, dissociation, ionization, surface chemical reaction, and radiation are referred to commonly as high-temperature real-gas phenomena. The phenomena cause changes in the dynamic behavior of the flow and the surface pressure and heat transfer distribution around the object. The character of a real gas is described by the internal degrees of freedom and state of constituent molecules; nitrogen and oxygen for air. The internal energy states, rotation, vibration and electronic, of the molecules are excited and, in the limit, the molecular bonds are exceeded and the gas dissociated into atomic and, possibly, ionic constituents. The process of energy transfer causing excitation, dissociation and recombination is a rate process controlled by particle collisions. Binary, two-body, collisions are sufficient to cause internal excitation, dissociation and ionization while three-body collisions are required to recombine the particles into molecular constituents. If the rates of energy transfer are fast with respect to the local fluid dynamic time scale the gas is in, or nearly in, equilibrium. If the energy transfer rates are very slow the gas can be described as frozen. In all other instances, wherein any of the energy exchange rates are comparable to the local fluid time scale, the gas will be thermally or chemically reacting and out of equilibrium. Real gas thermochemical nonequilibrium processes are important in the determination of aerodynamic heating; both convective (including wall catalytic effects) and radiative heating. To illustrate this we consider the hypervelocity flow over a bluff body typical of an atmospheric entry vehicle or an aerospace transfer vehicle (ASTV). The qualitative aspects of a hypersonic flow field over a bluff body are discussed in two parts, forebody and afterbody, with attention to which particular physical effects must be included in an analysis. This will indicate what type of numerical modeling will be adequate in each region of the flow. A bluff forebody flow field is dominated by the presence of the strong bow shock wave and the consequent heating, and chemical reaction of the gas. At high altitude hypersonic flight conditions the thermal excitation and chemical reaction of the gas occur slowly enough that a significant portion of the flow field is in a state of thermochemical nonequilibrium. A second important effect is the presence of the thick boundary layer along the forebody surface. In this region there are large thermal and chemical species gradients due to the interaction of the gas with the wall. Also at high altitudes the shock wave and the boundary layer may become so thick that they merge; in this case the entire shock layer is dominated by viscous effects.

Deiwert, George S.↗

Far Ultraviolet Astronomy

The Far Ultraviolet Spectroscopic Explorer (FUSE) is studying a wide range of astronomical problems in the 905-1187 Angstrom wavelength region through the use of high resolution spectroscopy. The FUSE bandpass forms a nearly optimal complement to the spectral coverage provided by the Hubble Space Telescope (HST), which extends down to approximately 1170 Angstroms. The photoionization threshold of atomic hydrogen (911 Angstroms) sets a natural short-wavelength limit for the FUV. FUSE was launched in June 1999 from Cape Canaveral, Florida, on a Delta II rocket into a 768 km circular orbit. Scientific observations started later that year. This spectral region is extremely rich in spectral diagnostics of astrophysical gases over a wide range of temperatures (100 K to over 10 million K). Important strong spectral lines in this wavelength range include those of neutral hydrogen, deuterium, nitrogen, oxygen, and argon (H I, D I, N I, O I, and Ar I), molecular hydrogen (H2), five-times ionized oxygen (O VI), and several ionization states of sulfur (S III - S VI). These elements are essential for understanding the origin and evolution of the chemical elements, the formation of stars and our Solar System, and the structure of galaxies, including our Milky Way. FUSE is one of NASA's Explorer missions and a cooperative project of NASA and the space agencies of Canada and France. These missions are smaller, more scientifically focused missions than the larger observatories, like Hubble and Chandra. FUSE was designed, built and operated for NASA by the Department of Physics and Astronomy at Johns Hopkins University. Hundreds of astronomers world-wide are using FUSE for a wide range of scientific research. Some of the important scientific discoveries from the first two years of the mission are described.

Sonneborn, George↗

Eta Carinae & the Homunculus: Far Infrared Sub-millimeter Spectral Lines detected with the Herschel Space Observatory

The evolved massive binary star η Carinae underwent eruptive mass loss events that formed the complex bi-polar “Homunculus” nebula harboring tens of solar masses of unusually nitrogen-rich gas and dust. Despite expectations for the presence of a significant molecular component to the gas, detections have been observationally challenged by limited access to the far-infrared and the intense thermal continuum. A spectral survey of the atomic and rotational molecular transitions was carried out with theHerschel Space Observatory, revealing a rich spectrum of broad emission lines originating in the ejecta. Velocity profiles of selected PACS lines correlate well with known substructures: Hi in the central core; NH and weak [Cii] within the Homunculus; and [Nii] emissions in fast-moving structures external to the Homunculus.We have identified transitions from [Oi], Hi, and 18 separate light C- and O-bearing molecules including CO, CH, CH+, and OH, and a wide set of N-bearing molecules,NH, NH+, N2H+, NH2, NH3, HCN, HNC, CN, and N2H+. Half of these are new detections unprecedented for any early-type massive star environment. A very low ratio[12C/13C]≤4 is estimated from five molecules and their isotopologues. We demonstrate that non-LTE effects due to the strong continuum are significant. Abundance patterns are consistent with line formation in regions of carbon and oxygen depletions with nitrogen enhancements, reflecting an evolved state of the erupting star with efficient transport of CNO-processed material to the outer layers. The results offer many opportunities for further observational and theoretical investigations of the molecular chemistry under extreme physical and chemical conditions around massive stars in their final stages of evolution.

T. R. Gull↗

Modeling Materials: Design for Planetary Entry, Electric Aircraft, and Beyond

NASA missions push the limits of what is possible. The development of high-performance materials must keep pace with the agency's demanding, cutting-edge applications. Researchers at NASA's Ames Research Center are performing multiscale computational modeling to accelerate development times and further the design of next-generation aerospace materials. Multiscale modeling combines several computationally intensive techniques ranging from the atomic level to the macroscale, passing output from one level as input to the next level. These methods are applicable to a wide variety of materials systems. For example: (a) Ultra-high-temperature ceramics for hypersonic aircraft-we utilized the full range of multiscale modeling to characterize thermal protection materials for faster, safer air- and spacecraft, (b) Planetary entry heat shields for space vehicles-we computed thermal and mechanical properties of ablative composites by combining several methods, from atomistic simulations to macroscale computations, (c) Advanced batteries for electric aircraft-we performed large-scale molecular dynamics simulations of advanced electrolytes for ultra-high-energy capacity batteries to enable long-distance electric aircraft service; and (d) Shape-memory alloys for high-efficiency aircraft-we used high-fidelity electronic structure calculations to determine phase diagrams in shape-memory transformations. Advances in high-performance computing have been critical to the development of multiscale materials modeling. We used nearly one million processor hours on NASA's Pleiades supercomputer to characterize electrolytes with a fidelity that would be otherwise impossible. For this and other projects, Pleiades enables us to push the physics and accuracy of our calculations to new levels.

Supercomputing↗

Photochemistry of Saturn's Atmosphere: Hydrocarbon Chemistry and Comparisons with ISO Observations - 1

To investigate the details of hydrocarbon photochemistry on Saturn, we have developed a one-dimensional diurnally averaged model that couples hydrocarbon and oxygen photochemistry, molecular and eddy diffusion, radiative transfer, and condensation. The model results are compared with observations from the Infrared Space Observatory (ISO) to place tighter constraints on molecular abundances, to better define Saturn's eddy diffusion coefficient profile, and to identify important chemical schemes that control the abundances of the observable hydrocarbons in Saturn's upper atmosphere. From the ISO observations, we determine that the column 12 densities of CH3, CH3C2H, and C4H2 above 10 mbar are 4 (sup +2) (sub -1.5) x 10 (exp 13) cm (sup -2), (1.1 plus or minus 0.3) x 10 (exp 15) cm (exp -2), and (1.2 plus or minus 0.3) x 10 (exp 14) cm (sup -2), respectively. The observed ISO emission features also indicate C2H2 mixing ratios of 1.2 (sup +0.9) (sub -0.6) x 10 (exp -6) at 0.3 mbar and (2.7 plus or minus 0.8) x 10 (exp -7) at 1.4 mbar, and a C2H6 mixing ratio of (9 plus or minus 2.5) x 10 (exp -6) at 0.5 mbar. Upper limits are provided for C2H4, CH2CCH2, C3H8, and C6H2 sensitivity of the model results to variations in the eddy diffusion coefficient profile, the solar flux, the CH4 photolysis branching ratios, the atomic hydrogen influx, and key reaction rates are discussed in detail. We find that C4H2 and CH3C2H are particularly good tracers of important chemical processes and physical conditions in Saturn's upper atmosphere, and C2H6 is a good tracer of the eddy diffusion coefficient in Saturn's lower stratosphere. The eddy diffusion coefficient must be smaller than approximately 3 x 10 (exp 4) sq cm s (sup -1) at pressures greater than 1 mbar in order to reproduce the C2H6 abundance inferred from ISO observations. The eddy diffusion coefficients in the upper stratosphere could be constrained by observations of CH3 radicals if the low-temperature chemistry of CH3 were better understood. We also discuss the implications of our modeling for aerosol formation in Saturn's lower stratosphere-diacetylene, butane, and water condense between approximately 1 and 300 mbar in our model and will dominate stratospheric haze formation at nonauroral latitudes. Our photochemical models will be useful for planning observational sequences and for analyzing data from the upcoming Cassini mission.

Moses, Julianne I.↗

High-Performance Computing Optimization for Aladyn – Adaptive Neural Network Molecular Dynamics Mini-Application

This report provides a description and performance evaluation of the optimization techniques for high performance computing (HPC) implementation of the open source Computational Materials mini-application Aladyn (https://github.com/nasa/aladyn). Aladyn is a basic molecular dynamics code written in FORTRAN 2003, which is designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to efficiently reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of the more expensive quantum mechanics-based calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory (DFT) method. While achieving orders of magnitude faster computational performance than DFT, the ANN-based approach was still very computationally demanding compared to the conventional approach of using empirically fitted energy functions. After its initial development, Aladyn was evaluated and optimized by experts at the NASA Advanced Supercomputing (NAS) division to exploit modern supercomputer architectures. The code has been optimized for execution on multicore central processing units (CPUs), including Intel® Skylake microarchitecture, and on graphic accelerators, such as Nvidia® V100 graphic processing units (GPUs), using Open Multi-Processing (OpenMP) and Open Accelerators (OpenACC) programming interfaces. The optimization achieved a speedup of 4.7 times the baseline version on CPU performance and an additional 2.4 times on CPU+GPU performance. Atomistic computer simulations are a fundamental tool in materials research to model material properties form physics-based first principles. Atomic interaction, governed by Quantum Mechanics (QM) require sophisticated and highly computationally demanding mathematical models to calculate [1]. Classical methods use approximate functional forms, empirically fitted through a set of variable parameters to emulate atomic energies as direct functions of atomic coordinates [2]. While empirical potentials are computationally much simpler, allowing simulations of large-scale systems of up to a trillion (1012) atoms [3], they are substantially less accurate compared to quantum calculations and applicable only to very specific atomic configurations or predefined crystallographic phases. A recently suggested approach is to use heuristic machine learning methods [4], such as those based on Adaptive Neural Networks (ANNs) to predict atomic energies, after being trained on a sufficiently large database of QM-calculated structures [5,6]. This approach reduces significantly the computational complexity, allowing for simulations of orders of magnitude larger systems compared to QM-based methods without compromising accuracy. Still, compared to classical methods using empirical energy functions, ANN methods remain two- to three orders of magnitude more computationally demanding. Hence, the computational cost of simulations, together with the need for extensive training of ANNs, still makes the practical implementation of ANN-based methods quite challenging. The purpose of the Aladyn mini-application software [7], available as open source at https://github.com/nasa/aladyn, is to be a testbed for exploring possible optimization strategies to develop highly scalable parallel algorithms for ANN-based atomistic simulations. Aladyn is aimed at utilizing the architecture of the high-end modern highperformance computing (HPC) hardware based on multicore central processing units (CPUs) equipped with graphic processing unit (GPU) accelerators. Specifically, the goal is to optimize the performance on a single HPC compute node, before implementing scaling to multi-node parallelization using message passing interface (MPI). At the same time, the open source code of Aladyn can serve as a training model for students and professors in academia.

Yamakov, Vesselin I.↗