Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular Langmuir Model”

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.

Nitric Acid Uptake on Subtropical Cirrus Cloud Particles

The redistribution of HNO3 via uptake and sedimentation by cirrus cloud particles is considered an important term in the upper tropospheric budget of reactive nitrogen. Numerous cirrus cloud encounters by the NASA WB-57F high-altitude research aircraft during CRYSTAL-FACE were accompanied by the observation of condensed-phase HNO3 with the NOAA chemical ionization mass spectrometer. The instrument measures HNO3 with two independent channels of detection connected to separate forward- and downward-facing inlets that allow a determination of the amount of HNO3 condensed on ice particles. Subtropical cirrus clouds, as indicated by the presence of ice particles, were observed coincident with condensed-phase HNO3 at temperatures of 197 K - 224 K and pressures of 122 hPa - 224 hPa. Maximum levels of condensed-phase HNO3 approached the gas-phase equivalent of 0.8 ppbv. Ice particle surface coverages as high as 1.4- 10(exp 14) molecules/sq cm were observed. A dissociative Langmuir adsorption model, when using an empirically derived HNO3 adsorption enthalpy of -11.0 kcal/mol, effectively describes the observed molecular coverages to within a factor of 5. The percentage of total HNO3 in the condensed phase ranged from near zero to 100% in the observed cirrus clouds. With volume-weighted mean particle diameters up to 700 pm and particle fall velocities up to 10 m/s, some observed clouds have significant potential to redistribute HNO3 in the upper troposphere.

Source record↗

Nitric Acid Uptake on Subtropical Cirrus Cloud Particles

The redistribution of HNO3 via uptake and sedimentation by cirrus cloud particles is considered an important term in the upper tropospheric budget of reactive nitrogen. Numerous cirrus cloud encounters by the NASA WB-57F high-altitude research aircraft during the Cirrus Regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE) were accompanied by the observation of condensed-phase HNO3 with the NOAA chemical ionization mass spectrometer. The instrument measures HNO3 with two independent channels of detection connected to separate forward and downward facing inlets that allow a determination of the amount of HNO3 condensed on ice particles. Subtropical cirrus clouds, as indicated by the presence of ice particles, were observed coincident with condensed-phase HNO3 at temperatures of 197-224 K and pressures of 122-224 hPa. Maximum levels of condensed-phase HNO3 approached the gas-phase equivalent of 0.8 ppbv. Ice particle surface coverages as high as 1.4 # 10(exp 14) molecules/ square cm were observed. A dissociative Langmuir adsorption model, when using an empirically derived HNO3 adsorption enthalpy of -11.0 kcal/mol, effectively describes the observed molecular coverages to within a factor of 5. The percentage of total HNO3 in the condensed phase ranged from near zero to 100% in the observed cirrus clouds. With volume-weighted mean particle diameters up to 700 ?m and particle fall velocities up to 10 m/s, some observed clouds have significant potential to redistribute HNO3 in the upper troposphere.

Popp, P. J.↗

A DSMC Surface Chemistry Model for Carbon-Based Ablators

A detailed molecular surface chemistry model for the DSMC (Direct Simulation Monte Carlo) method is proposed and implemented into the SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer) DSMC solver. Molchanova et al. constructed a molecular model for surface recombination in DSMC that includes different surface processes (adsorption, desoprtion, Eley-Rideal and Langmuir-Hinshelwood). All surface processes can be divided into two groups: surface mechanisms, which involve only the particle adsorbed by the surface (desorption and Langmuir-Hinshelwood), and impact mechanisms, which also involve gas-phase particles (adsorption, Eley-Rideal). Using a similar approach, the 14-reaction kinetic model of oxygen-carbon interaction suggested by Zhlukhtov and Abe, as well as more recent models by Alba et al., Poovathinghal et al., and a new model developed in the scope of this work, are implemented in SPARTA. The computational results for the different oxidation models are compared with experimental results from Murray et al. (oxidation of a vitreous carbon surface due to a hyperthermal beam of O and O2), with a particular focus on fluxes, angular and Time-Of-Flight distributions of scattered particles.

oxidation↗

An equation that describes material outgassing for contamination modeling

A generalization of the Clausius-Claperon equation for vapor pressure is made for an outgassing material. The expression is derived using Langmuir's equation for the outgassing rate of a material and using an empirical equation for the vapor pressure of a material as a function of its molecular weight and temperature. Also, outgassing rate equations are derived in terms of the vapor pressure of the outgassing material for three general geometries.

Heslin, T. M.↗

Development of an Effective Finite-rate Oxidation Model for NuSil-coated Charred Carbon Preform Ablators

Recently, a detailed effective finite-rate surface chemistry model was developed for the oxidation of FiberForm [1] using the molecular beam experimental data of Poovathingal et al., [2]. FiberForm is the major building block of the thermal protection system (TPS) material Phenolic Impregnated Carbon Ablator (PICA), commonly used by NASA. The surface chemistry model consists of detailed surface reaction mechanisms such as adsorption, desorption, and several types of Langmuir-Hinshelwood (LH) reactions to characterize the oxygen-carbon interactions at the surface. This model provides excellent agreement with the experimental data for oxidation product compositions and corresponding translational energy distributions. Further, an effective oxidation model was constructed that captures the equivalent interaction of oxygen inside the microstructure (including multiple surface collisions) via a modified reactivity for a smooth wall boundary condition. This enables the use of this model directly into Computational Fluid Dynamics (CFD) codes and Material Response (MR) codes to accurately simulate the gas–surface interactions within FiberForm without using the detailed micro-structure. However, this effective model is valid only for the virgin FiberForm. As the TPS material undergoes ablation, the preform carbon burns and turns into char. In addition, PICA is also coated with a protective silicon coating called NuSil for the purpose of mitigating the spread of phenolic dust, and limit contamination during clean room operations. Fig. 1 shows a X-ray microtomography image of a charred NuSil-coated PICA. The NuSil layer introduces species containing silicon into the product mixture. Further, the reactivity of the carbon within the char layer is different from the virgin FiberForm as shown in Fig. 2. Recently a new set of molecular beam experiments were performed on this NuSil-coated TPS material [3] with the same type of oxygen beam used in the previous experiments. Using the latest experimental data, the previously developed effective model will be extended to account for the charred carbon as well as the NuSil coating. Finally, this new effective model with three phases – preform carbon, char, and NuSil; will be compared and validated against the experimental product compositions.

K Swaminathan Gopalan↗

Mercury ion thruster research, 1977

The measured ion beam divergence characteristics of two and three-grid, multiaperture accelerator systems are presented. The effects of perveance, geometry, net-to-total accelerating voltage, discharge voltage and propellant are examined. The applicability of a model describing doubly-charged ion densities in mercury thrusters is demonstrated for an 8-cm diameter thruster. The results of detailed Langmuir probing of the interior of an operating cathode are given and used to determine the ionization fraction as a function of position upstream of the cathode orifice. A mathematical model of discharge chamber electron diffusion and collection processes is presented along with scaling laws useful in estimating performance of large diameter and/or high specific impluse thrusters. A model describing the production of ionized molecular nitrogen in ion thrusters is included.

Wilbur, P. J.↗

Catalytic oxidation for treatment of ECLSS and PMMS waste streams

Catalytic oxidation was added to the baseline multifiltration technology for use on the Space Station Freedom in order to convert low-molecular weight organic waste components such as alcohols, aldehydes, ketones, amides, and thiocarbamides to CO2 at low temperature (121 C), thereby reducing the total organic carbon (TOC) to below 500 ppb. The rate of reaction for the catalytic oxidation of aqueous organics to CO2 and water depends primarily upon the catalyst, temperature, and concentration of reactants. This paper describes a kinetic study conducted to determine the impact of each of these parameters upon the reaction rate. The results indicate that a classic kinetic model, the Langmuir-Hinshelwood rate equation for heterogeneous catalysis, can accurately represent the functional dependencies of this rate.

Akse, James R.↗

Catalytic methods using molecular oxygen for treatment of PMMS and ECLSS waste streams, volume 2

Catalytic oxidation has proven to be an effective addition to the baseline sorption, ion exchange water reclamation technology which will be used on Space Station Freedom (SSF). Low molecular weight, polar organics such as alcohols, aldehydes, ketones, amides, and thiocarbamides which are poorly removed by the baseline multifiltration (MF) technology can be oxidized to carbon dioxide at low temperature (121 C). The catalytic oxidation process by itself can reduce the Total Organic Carbon (TOC) to below 500 ppb for solutions designed to model these waste waters. Individual challenges by selected contaminants have shown only moderate selectivity towards particular organic species. The combined technology is applicable to the more complex waste water generated in the Process Materials Management System (PMMS) and Environmental Control and Life Support System (ECLSS) aboard SSF. During the phase 3 Core Module Integrated Facility (CMIF) water recovery tests at NASA MSFC, real hygiene waste water and humidity condensate were processed to meet potable specifications by the combined technology. A kinetic study of catalytic oxidation demonstrates that the Langmuir-Hinshelwood rate equation for heterogeneous catalysts accurately represent the kinetic behavior. From this relationship, activation energy and rate constants for acetone were determined.

Akse, James R.↗

Implementation of Active Sites in DSMC to Capture Pitting of Oxidizing Carbon Materials

In this work we demonstrate a newly developed capability to capture pitting of carbon fibers in DSMC simulations, specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. State-of-the-art reactive surface models in DSMC compute collision dependent carbon consumption rates (usually through desorption of CO) based on a set of surface reactions that has been derived from molecular beam experiments. The reactivity on each carbon surface element is constant in those models, such that the carbon surface recedes uniformly as a result of ablation. However, it is well known that in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. These defective sites have a much higher reactivity than the average sites (2-3 orders of magnitude) and are first to react during ablation leading to its removal. This causes all the neighboring atoms to be defective and increase their reactivity, thus leading to the localized carbon removal around these ”active” sites. In this manner, these highly reactive defective sites serve as nucleation sites for the formation and growth of etch pits with potentially detrimental effects on the structural integrity. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various reaction mechanisms such as adsorption, desorption, Eley-Rideal (ER) and Langmuir-Hinshelwood (LH) mechanisms. Within this framework, we have implemented the capability of a single surface having multiple site sets with different reactivities. Using this feature, we can simulate the presence of active sites on carbon surfaces, whose reactivity is much greater than an average site as a result of defects. We have implemented the active site fraction as a property of surface elements within SPARTA, which is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface, and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we are able to capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗

Flow Characteristics of UHV Nitrogen Arcjet

Flow characteristics in the nozzle of a nitrogen arcjet have been simulated by a model. Electrons, ions, and neutral atoms and molecules are represented as separate species. Dissociation, ionization, and recombination are treated as separate reactions. Thermal non-equilibrium is represented by a two-temperature model. Energy input to plasma from the arc is modeled as a source determined by local current and electrical conductivity. Momentum and energy losses and recombinations at the nozzle wall are included as source/sink determined by fluxes to the wall. The results show essentially frozen gas-phase reactions and thermal nonequilibrium due to Joule heating of electrons and in the expansion part of the nozzle. Wall interactions have strong effects on the results, indicating that they play important roles in the flow due to very high surface area relative to the volume of the plasma. For example, plasma velocity decreases from 6000 m per s without wall interactions to 3000 m per s with modest amount of wall interactions. These results also show reasonable agreement with optical emission measurements, which confirms that the arcjet plasma is far from LTE. The spectra suggest nitrogen dissociation levels of 0.3% - 9%, depending on nitrogen flow rate and arc plasma current, which also determine the relative amounts of excited atomic and molecular nitrogen. Langmuir probe studies of the source show that electron and ion fluxes increase with increasing power, and that the ion energy distribution shifts to lower energies. Typical ion fluxes were on the order of 4E-9 A/cm**2 with a maximum ion kinetic energy of 3.5eV. The median electron energy was 1eV, with a maximum of less than 4eV.

Chang, C. H.↗

Implementation of active sites to capture pitting of oxidizing carbon materials in DSMC.

In this work we demonstrate a newly developed capability to capture pitting of carbon fibers in DSMC simulations, specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code [1]. State-of-the-art reactive surface models in DSMC compute collision dependent carbon consumption rates (usually through desorption of CO) based on a set of surface reactions that has been derived from molecular beam experiments [2]. The reactivity on each carbon surface element is constant in those models, such that the carbon surface recedes uniformly as a result of ablation. However, it is well known that in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale [3]. These defective sites have a much higher reactivity than the average sites (2-3 orders of magnitude) and are first to react during ablation leading to its removal. This causes all the neighboring atoms to be defective and increase their reactivity, thus leading to the localized carbon removal around these ”active” sites (as shown in Fig. 1). In this manner, these highly reactive defective sites serve as nucleation sites for the formation and growth of etch pits with potentially detrimental effects on the structural integrity. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various reaction mechanisms such as adsorption, desorption, Eley-Rideal (ER) and Langmuir-Hinshelwood (LH) mechanisms [4]. Within this framework, we have implemented the capability of a single surface having multiple site sets with different reactivities. Using this feature, we can simulate the presence of active sites on carbon surfaces, whose reactivity is much greater than an average site as a result of defects. We have implemented the active site fraction as a property of surface elements within SPARTA, which is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface, and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we are able to capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗

Diagnostic for Plasma Enhanced Chemical Vapor Deposition and Etch Systems

In order to meet NASA's requirements for the rapid development and validation of future generation electronic devices as well as associated materials and processes, enabling technologies ion the processing of semiconductor materials arising from understanding etch chemistries are being developed through a research collaboration between Stanford University and NASA-Ames Research Center, Although a great deal of laboratory-scale research has been performed on many of materials processing plasmas, little is known about the gas-phase and surface chemical reactions that are critical in many etch and deposition processes, and how these reactions are influenced by the variation in operating conditions. In addition, many plasma-based processes suffer from stability and reliability problems leading to a compromise in performance and a potentially increased cost for the semiconductor manufacturing industry. Such a lack of understanding has hindered the development of process models that can aid in the scaling and improvement of plasma etch and deposition systems. The research described involves the study of plasmas used in semiconductor processes. An inductively coupled plasma (ICP) source in place of the standard upper electrode assembly of the Gaseous Electronics Conference (GEC) radio-frequency (RF) Reference Cell is used to investigate the discharge characteristics and chemistries. This ICP source generates plasmas with higher electron densities (approximately 10(exp 12)/cu cm) and lower operating pressures (approximately 7 mTorr) than obtainable with the original parallel-plate version of the GEC Cell. This expanded operating regime is more relevant to new generations of industrial plasma systems being used by the microelectronics industry. The motivation for this study is to develop an understanding of the physical phenomena involved in plasma processing and to measure much needed fundamental parameters, such as gas-phase and surface reaction rates. species concentration, temperature, ion energy distribution, and electron number density. A wide variety of diagnostic techniques are under development through this consortium grant to measure these parameters. including molecular beam mass spectrometry (MBMS). Fourier transform infrared (FTIR) spectroscopy, broadband ultraviolet (UV) absorption spectroscopy, a compensated Langmuir probe. Additional diagnostics. Such as microwave interferometry and microwave absorption for measurements of plasma density and radical concentrations are also planned.

Cappelli, Mark A.↗

Products of Dissociative Recombination in the Ionosphere

SRI International undertook a novel experimental measurement of the product states formed by dissociative recombination (DR) of O2(+), NO(+), and N2(+) as a function of both electron energy and reactant ion vibrational level. For these measurements we used a recently developed experimental technique for measuring dissociation product distributions that allows both the branching ratios to be accurately determined and the electronic and rovibrational state composition of the reactant ions to be specified. DR is the dominant electron loss mechanism in all regions of the ionosphere. In this process, electron attachment to the molecular ion produces an unstable neutral molecule that rapidly dissociates. For a molecular ion such as O2(+), the dissociation recombination reaction is (1) O2(+) + e yields O + O + W. The atomic products of this reaction, in this case two oxygen atoms, can be produced in a variety of excited states and with a variety of kinetic energies, as represented by W in Eq. (1). These atoms are not only active in the neutral chemistry of the ionosphere, but are also especially important because their optical emissions are often used to infer in situ concentrations of the parent molecular ion and ambient electron densities. Many laboratory measurements have been made of DR reaction rates under a wide range of electron temperatures, but very little is known about the actual distributions among the final states of the atomic products. This lack of knowledge seriously limits the validity and effectiveness of efforts to model both natural and man-made ionospheric disturbances. Bates recently identified major deficiencies in the currently accepted branching ratios for O2(+) as they relate to blue and green line emission measurements in the nocturnal F-region. During our two-year effort, we partially satisfied our ambitious goals. We constructed and operated a variable pressure, electron-impact ion source and a high pressure, hollow-cathode discharge ion source for O2(+), NO(+), and N2(+) beams. Translational spectroscopy of the products of dissociative charge transfer in Cs vapor was used to accurately assay the composition of the O2(+) and NO(+) beams and to develop a methodology for the vibrationally controlled preparation of the ground state ion beams. Attempts to assay the N2(+) beam revealed a novel two-electron process in the charge transfer reactions. A coaxial electron gun for the DR measurements was constructed following an extensive numerical design of the fields. Tests of the gun, however, found substantial perturbations of the magnetic fields by the soft iron (CMI-C) assembly containing the Langmuir probe that locates the electron beam. Hydrogen annealing of the iron failed to eliminate the field perturbations, necessitating the removal of the probe assembly. During this work on the coaxial electron gun, we discovered that predissociated high Rydberg states of O2 could be produced by subjecting the molecules to a sudden perturbation by an electromagnetic field. This technique allowed a measurement of the product branching to the atomic limits for the lowest seven vibrational levels of O2(+).

Cosby, Philip↗

Surfactant Facilitated Spreading of Aqueous Drops on Hydrophobic Surfaces

Microgravity technologies often require aqueous phases to spread over nonwetting hydrophobic solid surfaces. Surfactants facilitate the wetting of water on hydrophobic surfaces by adsorbing on the water/air and hydrophobic solid/water interfaces and lowering the surface tensions of these interfaces. The tension reductions decrease the contact angle, which increases the equilibrium wetted area. Hydrocarbon surfactants; (i.e., amphiphiles with a hydrophobic moiety consisting of an extended chain of (aliphatic) methylene -CH2- groups attached to a large polar group to give aqueous solubility) are capable of reducing the contact angles on surfaces which are not very hydrophobic, but do not reduce significantly the contact angles of the very hydrophobic surfaces such as parafilm, polyethylene or self assembled monolayers. Trisiloxane surfactants (amphiphiles with a hydrophobe consisting of methyl groups linked to a trisiloxane backbone in the form of a disk ((CH3)3-Si-O-Si-O-Si(CH3)3) and an extended ethoxylate (-(OCH2CH2)a-) polar group in the form of a chain with four or eight units) can significantly reduce the contact angle of water on a very hydrophobic surface and cause rapid and complete (or nearly complete) spreading (termed superspreading). The overall goal of the research described in this proposal is to establish and verify a theory for how trisiloxanes cause superspreading, and then use this knowledge as a guide to developing more general hydrocarbon based surfactant systems which superspread. We propose that the trisiloxane surfactants superspread because their structure allows them to strongly lower the high hydrophobic solid/aqueous tension when they adsorb to the solid surface. When the siloxane adsorbs, the hydrophobic disk parts of the molecule adsorb onto the surface removing the surface water. Since the cross-sectional area of the disk is larger than that of the extended ethoxylate chain, the disks can form a space-filling mat on the surface which removes a significant amount of the surface water. In this presentation, we report the results of measurements of the molecular packing and rates of kinetic exchange of the trisiloxane surfactants at the air/water interface in order to confirm our picture of trisiloxane packing, and provide additional insight into the superspreading process. We used the pendant bubble technique as a Langmuir trough to measure the trisiloxane equation of state which relates the tension to the surface concentration. From these measurements we obtain accurate values for the maximum packing density. We find that trisiloxanes with 4 and 8 ethoxylate groups have the same maximum packing concentration, indicating that the maximum packing is controlled by the cross section of the head group. For trisiloxanes with larger than eight ethoxylates, the maximum packing increases with ethoxylate number, indicating that the disposition of the ethoxylate chain (i.e., its effective size) is controlling. This supports our picture of superspreading: The superspreading ability of trisiloxanes decreases considerably for trisiloxanes with larger than eight ethoxylates; the packing measurements indicate that for the higher ethoxylate number trisiloxanes, the compact nonpolar head groups are pushed apart by the ethoxylate chain. They leave spaces of surface water on adsorption and do not lower the solid tension as much as their lower chain analogues. Finally the report measurements of the dynamic tension reduction accompanying the adsorption of trisiloxanes onto an initially clean interface using the pendant bubble technique, and we obtain from these relaxations, the equation of state and a mass transfer model, the rate constants for kinetic exchange. We find that the rate constants for desorption of trisiloxanes are generally much slower than for analogous aliphatic polyethoxylate surfactants with identical ethoxylate chain lengths. When an aqueous drop of a superspreader solution is placed on a hydrophobic surface and begins to spread, the lower desorption rates allows the tension at the drop center to remain reduced relative to the tension of the expanding periphery, thus strengthening Marangoni forces which can assist the spreading. Marangoni forces can be especially significant in the case of superspreaders because their maximum equilibrium reduction in tension is several dynes/cm lower than for aliphatic surfactants.

Kumar, Nitin↗

Influence of Nucleation Mechanisms on the Radiative Properties of Deep Convective Clouds and Subvisible Cirrus in CRYSTAL/FACE

During the past few years we have conducted work on several different topics, as reflected by our publications. As one of the Co-Project scientists for The Cirrus Regional Study of Tropical Anvils and Cirrus Layers - Florida Area Cirrus Experiment (CRYSTAL FACE) we worked to help design the mission and then conduct it in the field. Another major activity during the past two years has been to pull together various groups to formulate plans for follow on missions to CRYSTAL FACE. We organized a workshop at the University of Colorado during the summer of 2003 to assess the best locations for future missions. Working with a group of about 10 scientists from around the country we prepared a science-planning document (Tropical Composition, Cloud and Climate Coupling Experiment (TC(sup 4)) that outlined the rationale, locations, strategy to accomplish the goals, and possible payloads for a set of three tropical missions. We also prepared background materials for various NRAs being prepared at NASA Headquarters for missions in Costa Rica, Darwin and Guam. In conjunction with the group at NASA Ames we have helped build a new numerical model for deep convection and have applied that model to simulate the CRYSTAL data. Our goal in particular has been to better understand how convection distributes water vapor isotopes. CRYSTAL observations of water isotopes are very different from those suggested by previous workers who assumed the isotopes would obey Rayleigh fractionation. The water isotope study has several implications. First it is a check on the realism of the deep convection model. Second, the isotopes are a measure of the precipitation removal in the atmosphere. Hence they provide a constraint on a parameter that is difficult to otherwise measure. Finally it has been suggested that isotopes may be the key to unraveling the water transport into the stratosphere and upper troposphere. Such transport is critical both for the radiation balance and for stratospheric chemistry. Ours is the first model that is able to treat this transport. Our initial results have just been submitted to Geophys. Res. Lett (Smith et al., 2005). Essentially we are able to explain the vertical profiles of isotopes in the tropical tropopause transition layer. We are also able to account for stratospheric humidity and isotope abundances with this model. We have also been heavily involved in trying to improve our understanding of nitric acid condensation on ice. Gao et al (2004) have shown that water supersaturations above ice occur when the atmosphere is supersaturated with respect to nitric acid trihydrate. As one of the co-authors of that work, we suggested the mechanism that may explain why this is occurring. Essentially, ice does not like to grow near unit supersaturation, but does so because the water molecules can find sites on the ice surface to attach themselves to before they fly off the ice surface. This phenomena was well known in the 1960s when it was a source of debate about whether condensation and evaporation coefficients for ice would be the same. Evaporation does not require any molecular orientation, while condensation does, so it was possible that the coefficients would differ. They don't differ because the water molecules rapidly move across the surface and find places to attach. Nitric acid may be occupying these preferred sites and therefore the water molecules can't find a desirable place to attach. We anticipate that this research will be the subject of laboratory work during the coming few years. Another possibility that has been suggested is that cubic ice is forming in clouds. We have measured the vapor pressure of cubic ice, and plan to publish that result in the next few months. We have also been working on additional aspects of the condensation of nitric acid on ice. With Y. Kondo we studied the condensation of NOy on ice using the SOLVE data. Gamblin et al. have continued this work. The CRYSTAL NOy and HNO, groups have shown th their data can be fit using standard Langmuir isotherms as suggested in some, but not all, laboratory studies. We have found in the SOLVE data set that this is not the case. Moreover some laboratory studies show there are important kinetic effects that may be occurring in the atmosphere limiting the transfer of nitric acid to the ice. The SOLVE data seem consistent with these studies. We are currently re-analyzing the CRYSTAL data to look for these kinetic effects. There are a number of implications of these studies. One of the more interesting is that the nitric acid coating on ice can be used as a cloud clock to determine how long the cloud parcel has been in existence. We have also been involved with several laboratory studies. We have worked to improve the database on ice optical constants, which are critical for remote sensing. We have also studied the ways in which ice nucleates on clays. We suspect now that the standard theories used for depositional ice nucleation are completely incorrect. Further work will be needed to develop a new theory.

Toon, Owen B.↗