On the spectral features of dangling bonds in CH4/H2O amorphous ice mixtures.

Dangling bond (DB) bands in IR spectra, above 3600 cm-1, are a source of information on structural properties of amorphous water ice, and especially on ice mixtures of water and other frozen gases. We deal in this paper with the spectroscopic behavior of DB bands of CH4/H2O mixtures. We use ab initio methodology to predict theoretical results which are compared with experimental results. Our model mixtures are created by inserting a variable number of molecules of either species into a cell of appropriate size to reach an initial density of 1 g cm-3, which can be modified by including an empty space at the top, to simulate pores. The cell is taken as a unit cell for a solid state calculation The structure of the mixture is optimized and the IR spectrum is calculated for the converged geometry. We find two different kinds of dangling bonds, in which the O-H stretching responsible for this mode is directed either to an empty space of a pore or towards a nearby CH4 molecule, with which some interaction takes place. The spectral characteristics of these two DB types are clearly different, and follow satisfactorily the pattern observed in experimental spectra. Estimated band strengths for these DB bands are given for the first time.

Thermodynamic, Raman Spectroscopic, and UV-Visible Optical Characterization of the Deltic, Squaric, and Croconic Cyclic Oxocarbon Acids.

A precise and complete thermodynamic, Raman spectroscopic, and ultraviolet-visible (UV-vis) optical characterization of the deltic, squaric, and croconic cyclic oxocarbon acids is obtained using theoretical solid-state methods employing very demanding calculation parameters. The computed fundamental thermodynamic properties include the isobaric specific heat, the entropy, the enthalpy, and the Gibbs free energy as a function of temperature. The calculated specific heats at 298.15 K of the deltic, squaric, and croconic acids are 89.7, 111.2, and 133.2 J mol-1 K-1, respectively, and the corresponding entropies are 98.3, 117.3, and 136.5 J mol-1 K-1. The only value of these properties known from experimental measurements is the specific heat of the squaric acid, which differs from the computed value at 315 K by about 4.9%. The calculated values of the thermodynamic properties are then used to determine the thermodynamic properties of formation of these materials in terms of the elements. As an application of the calculated thermodynamic properties of formation, the Gibbs free energies of reaction and associated reaction constants are evaluated for the reactions of thermal decomposition and complete combustion of the squaric and croconic acids and the reaction of interconversion between them. The only available experimental values of these properties, namely, the enthalpies of combustion of squaric and croconic acids at room temperature, are reproduced theoretically with high accuracy. The Raman spectra of these materials are also computed using density functional perturbation theory. The analysis of the theoretical Raman spectra of these materials points out to significant differences with respect to their usual empirical assignment. Therefore, the Raman spectra of these materials are fully reassigned. Finally, the ultraviolet-visible (UV-vis) optical properties of the deltic, squaric, and croconic acids are computed. The UV-vis absorption spectrum of the croconic acid in the spectral region 225-425 nm and the UV absorption spectrum of the squaric acid in the region 200-350 nm, which had previously been measured experimentally, are well reproduced. The corresponding spectrum for the deltic acid and the reflectivity, optical conductivity, dielectric, refractive index, and loss optical functions of the three materials, which had never been published as far as we know, are reported as a function of the wavelength of incident radiation in the range 200-750 nm. The origin of the peaks in the absorption spectra, which had not been analyzed so far, is unveiled here by examining the interband electronic transitions in these materials.

Prediction of the near-IR spectra of ices by ab initio molecular dynamics.

A method to predict the near-infrared spectra of amorphous solids by means of ab initio molecular dynamics is presented. These solids can simulate molecular ices. To test the method, mixtures of methane, water and nitrogen are generated as amorphous samples of various concentrations. The full theoretical treatment includes as a first step, the optimization of their geometrical structure for a range of densities, after which, the most stable systems are taken as initial structures for molecular dynamics, performed at 200 K in trajectories of 4 ps duration with a 0.2 fs time step. All the dynamics are carried out using the first principles method, solving the quantum problem for the electrons using density-functional theory (DFT), and integrating the DFT forces, following the Born-Oppenheimer dynamics. After the dynamics, near-IR spectra are predicted by the Fourier transform of the macroscopic polarization autocorrelation function. The calculated spectra are compared with the experimental spectra of ice mixtures of CH4 and H2O recorded in our laboratory, and with some spectra recorded by the New Horizons mission on Pluto.

Vibrational spectra and physico-chemical properties of astrophysical analogs.

We undertake in this paper a theoretical study based on DFT methodology of amorphous solids formed by methane, water and nitrogen in a ratio of 1 : 3 : 3. By varying the size of the cell containing this mixture of molecules, we study the effect of the corresponding cell volume and density on the predicted IR spectra, in particular on the hydrogen bond modes. Also the relative stability of the structures as a function of the density is studied. We have enclosed a large density range, from a very low value that simulates a gas-phase mixture, to values corresponding to solids under fairly high internal stress, with an intermediate range that could be expected to cover the values of mixtures at astronomical conditions. The variation of the energy at constant temperature with the volume of the unit cell fits well to a Morse function, which allows finding an equation of state for the material in the range of volumes studied here.

Simulations and spectra of water in CO matrices.

Models for the inclusion of water molecules in carbon monoxide matrices are developed using density functional theory applied to amorphous solid systems. The models cover a large range of systems for smaller or larger CO matrices with different water content, consisting of either individual H2O molecules or small clusters linked by H-bonds. The vibrational spectra of the samples are predicted at the minimum of their potential energy surface. The spectra allow instances where the water molecules remain isolated or form aggregates to be discerned, and they also provide an indication of the strength of the H-bonding, when present. The calculations support recent experimental observations that linked IR bands at 3707 cm-1 and 3617 cm-1 to the presence of unbound water molecules in water-poor CO/H2O mixed ices. Assignment of some observed bands to water dimers or trimers is suggested as well. The residual static pressure in fixed-volume simulation cells is also calculated.

Structure and infrared spectra of hydrocarbon interstellar dust analogs.

A theoretical study of the structure and mid infrared (IR) spectra of interstellar hydrocarbon dust analogs is presented, based on DFT calculations of amorphous solids. The basic molecular structures for these solids are taken from two competing literature models. The first model considers small aromatic units linked by aliphatic chains. The second one assumes a polyaromatic core with hydrogen and methyl substituents at the edges. The calculated spectra are in reasonably good agreement with those of aliphatic-rich and graphitic-rich samples of hydrogenated amorphous carbon (HAC) generated in our laboratory. The theoretical analysis allows the assignment of the main vibrations in the HAC spectra and shows that there is a large degree of mode mixing. The calculated spectra show a marked dependence on the density of the model solids, which evinces the strong influence of the environment on the strengths of the vibrational modes. The present results indicate that the current procedure of estimating the hydrogen and graphitic content of HAC samples through the decomposition of IR features into vibrational modes of individual functional groups is problematic owing to the mentioned mode mixing and to the difficulty of assigning reliable and unique band strengths to the various molecular vibrations. Current band strengths from the literature might overestimate polyaromatic structures. Comparison with astronomical observations suggests that the average structure of carbonaceous dust in the diffuse interstellar medium lies probably in between those of the two models considered, though closer to the more aliphatic structure.

Theoretical model of the interaction of glycine with hydrogenated amorphous carbon (HAC).

A theoretical model of hydrogenated amorphous carbon (HAC) is developed and applied to study the interaction of glycine with HAC surfaces at astronomical temperatures. Two models with different H content are tried for the HAC surface. The theory is applied at the Density Functional Theory (DFT) level, including a semiempirical dispersion correlation potential, d-DFT or Grimme DFT-D2. The level of theory is tested on glycine adsorption on a Si(001) surface. Crystalline glycine is also studied in its two stable phases, α and ß, and the metastable γ phase. For the adsorption on Si or HAC surfaces, molecular glycine is introduced in the neutral and zwitterionic forms, and the most stable configurations are searched. All theoretical predictions are checked against experimental observations. HAC films are prepared by plasma enhanced vapor deposition at room temperature. Glycine is deposited at 20 K into a high vacuum, cold temperature chamber, to simulate astronomical conditions. Adsorption takes place through the acidic group COO(-) and when several glycine molecules are present, they form H-bond chains among them. Comparison between experiments and predictions suggests that a possible way to improve the theoretical model would require the introduction of aliphatic chains or a polycyclic aromatic core. The lack of previous models to study the interaction of amino-acids with HAC surfaces provides a motivation for this work.

Stability of carbonaceous dust analogues and glycine under UV irradiation and electron bombardment.

The effect of UV photon (120-200 nm) and electron (2 keV) irradiation of analogues of interstellar carbonaceous dust and of glycine were investigated by means of IR spectroscopy. Films of hydrogenated amorphous carbon (HAC), taken as dust analogues, were found to be stable under UV photon and electron bombardment. High fluences of photons and electrons, of the order of 10(19) cm(-2), were needed for a film depletion of a few percent. UV photons were energetically more effective than electrons for depletion and led to a certain dehydrogenation of the HAC samples, whereas electrons led seemingly to a gradual erosion with no appreciable changes in the hydrocarbon structure. The rates of change observed may be relevant over the lifetime of a diffuse cloud, but cannot account for the rapid changes in hydrocarbon IR bands during the evolution of some proto-planetary nebulae. Glycine samples under the same photon and electron fluxes decay at a much faster rate, but tend usually to an equilibrium value different from zero, especially at low temperatures. Reversible reactions re-forming glycine, or the build-up of less transparent products, could explain this behavior. CO2 and methylamine were identified as UV photoproducts. Electron irradiation led to a gradual disappearance of the glycine layers, also with formation of CO2. No other reaction products were clearly identified. The thicker glycine layers (a few hundred nm) were not wholly depleted, but a film of the order of the electron penetration depth (80 nm), was totally destroyed with an electron fluence of -1 x 10(18) cm(-2). A 60 nm ice layer on top of glycine provided only partial shielding from the 2 keV electrons. From an energetic point of view, 2 keV electrons are less efficient than UV photons and, according to literature data, much less efficient than MeV protons for the destruction of glycine. The use of keV electrons to simulate effects of cosmic rays on analogues of interstellar grains should be taken with care, due to the low penetration depths of electrons in many samples of interest.

The formation of carbamate ions in interstellar ice analogues.

Carbon dioxide and ammonia are two of the most abundant species in astrophysical media, where they can react in the solid phase under certain conditions. This contribution presents a study of this reaction both in the presence of water and for anhydrous samples. It is shown that after deposition at 15 K, the reaction can start by warming the deposit, and the process continues on up to a temperature of 220 K. Reaction products are studied using infrared spectroscopy and mass spectrometry. For anhydrous samples, a 2 : 1 stoichiometry mixture of NH3 : CO2 gives the highest yield of products. The reaction is favored when a small amount of water is present, which enables ammonia and carbon dioxide molecules to collide within the pores and channels of the amorphous water solid. Large concentration of water, on the other hand, hampers such collisions. The main reaction product is found to be ammonium carbamate, but also carbamic acid is formed, and, in the presence of water, ammonium bicarbonate is produced as well. Theoretical calculations are carried out to provide the basis for the assignment of the spectra. Some of the experiments presented in this contribution consist of the generation of a compact water ice matrix where the carbamate and ammonium ions are embedded. If such a system was found in astrophysical media, it is shown that the ammonium ion could not be detected, whereas two infrared features of the carbamate ion in the 1040 to 1115 cm(-1) (9 to 9.6 µm) region could enable the observation of this species.

Crystallization of CO2 ice and the absence of amorphous CO2 ice in space.

Carbon dioxide (CO2) is one of the most relevant and abundant species in astrophysical and atmospheric media. In particular, CO2 ice is present in several solar system bodies, as well as in interstellar and circumstellar ice mantles. The amount of CO2 in ice mantles and the presence of pure CO2 ice are significant indicators of the temperature history of dust in protostars. It is therefore important to know if CO2 is mixed with other molecules in the ice matrix or segregated and whether it is present in an amorphous or crystalline form. We apply a multidisciplinary approach involving IR spectroscopy in the laboratory, theoretical modeling of solid structures, and comparison with astronomical observations. We generate an unprecedented highly amorphous CO2 ice and study its crystallization both by thermal annealing and by slow accumulation of monolayers from the gas phase under an ultrahigh vacuum. Structural changes are followed by IR spectroscopy. We also devise theoretical models to reproduce different CO2 ice structures. We detect a preferential in-plane orientation of some vibrational modes of crystalline CO2. We identify the IR features of amorphous CO2 ice, and, in particular, we provide a theoretical explanation for a band at 2,328 cm(-1) that dominates the spectrum of the amorphous phase and disappears when the crystallization is complete. Our results allow us to rule out the presence of pure and amorphous CO2 ice in space based on the observations available so far, supporting our current view of the evolution of CO2 ice.

Proton transfer and autoionization in HNO3·HCl·(H2O)n particles.

The structure and spectroscopic properties of clusters of HNO(3)·HCl·(H(2)O)(n), with n = 1 to 6, have been calculated at the MP2/aug-cc-pVDZ level of theory. Altogether 22 different clusters have been found as stable structures, with minima in their potential energy surfaces. The clusters can be grouped in families with the same number of water molecules, and with close aggregation energies within each family. The addition of each new water molecule increments the aggregation energy of the clusters by a nearly constant value of 76.2 ± 0.1 Hartree. The proton transfer parameter and the coordination number of HNO(3) and HCl in each cluster have been evaluated, and the wavenumber shifts for the X(-)-H(+) vibration from the corresponding mode in the isolated molecules have also been predicted. These values allow classification of the acidic species in the clusters into three types, characterized by the strength of the hydrogen bond and the degree of ionization. A correspondence is found between the coordination number of HNO(3) and the magnitude of the X(-)-H(+) vibrational shift.

An infrared study of solid glycine in environments of astrophysical relevance.

The conversion from neutral to zwitterionic glycine is studied using infrared spectroscopy from the point of view of the interactions of this molecule with polar (water) and non-polar (CO(2), CH(4)) surroundings. Such environments could be found on astronomical or astrophysical matter. The samples are prepared by vapour-deposition on a cold substrate (25 K), and then heated up to sublimation temperatures of the co-deposited species. At 25 K, the neutral species is favoured over the zwitterionic form in non-polar environments, whereas for pure glycine, or in glycine/water mixtures, the dominant species is the latter. The conversion is easily followed by the weakening of two infrared bands in the mid-IR region, associated to the neutral structure. Theoretical calculations are performed on crystalline glycine and on molecular glycine, both isolated and surrounded by water. Spectra predicted from these calculations are in reasonable agreement with the experimental spectra, and provide a basis to the assignments. Different spectral features are suggested as probes for the presence of glycine in astrophysical media, depending on its form (neutral or zwitterionic), their temperature and composition.

Interaction of CH4 and H2O in ice mixtures.

Ice mixtures of methane and water are investigated by means of IR spectroscopy in the 14-60 K range. The spectroscopic research is focused on the symmetry-forbidden nu(1) band of CH(4) and the dangling bond bands of water. The nu(1) band is visible in the spectra of the mixtures, revealing a distorted methane structure which co-exists with the normal crystalline methane. The water dangling bond bands are found to increase their intensity and appear at red-shifted frequency when distorted methane is present. Methane adsorbed on water micropores or trapped inside the amorphous solid water structure is assumed to be responsible for these effects. CH(4) mobility in water ice depends on the deposition method used to prepare the samples and on the temperature. After warming the samples to 60 K, above the methane sublimation point, a fraction of CH(4) is retained in the water ice. An adsorption isotherm analysis is performed yielding the estimation of the desorption energy of CH(4) on H(2)O amorphous surfaces.

Theoretical study on hydrogen-bond effects in IR spectra of high- and low-temperature phases of nitric acid dihydrate.

The low- and high-temperature phases (alpha and beta, respectively) of solid nitric acid dihydrate (NAD) are studied in depth by DFT methods. Each phase contains two types of complex structures (H(3)O(+)) x (H(2)O), designated A and B, with different hydrogen-bonding (HB) characteristics. The theoretical study reveals that type A complexes are weakly bound and could be described as (H(3)O)(+) and H(2)O aggregates, with decoupled vibrational modes, whereas in type B structures the proton is situated close to the centre of the O...O bond and induces strong vibrational coupling. The proton-transfer mode is predicted at quite different wavenumbers in each complex, which provides an important differentiating spectral feature, together with splitting of some bands in beta-NAD. Theoretical spectra are estimated by using two GGA parameterizations, namely, PBE and BLYP. The potential-energy surface for each type of HB in NAD is also studied, as is the spectral influence of displacement of the shared H atom along the O-O bond. The results are compared to literature infrared spectra recorded by different techniques, namely, transmission and reflection-absorption, with both normal and tilted incident radiation. This work provides a thorough assignment of the observed spectra, and predictions for some spectra not yet available. The usefulness of high-level theoretical calculations as performed herein to discriminate between two phases of a solid crystal is thus evidenced.

Phases of solid methanol.

The solid phases of methanol were investigated using IR spectroscopy and numerical calculations with the SIESTA method. Improved spectra are reported of amorphous methanol at 90 K, and in particular of the alpha and beta phases at 130 and 165 K, respectively, with assignments of bands not previously measured. The main features of the spectra of each phase are discussed and compared. A study of spectral changes with temperature leads to the conclusion that the metastable phase previously reported might be a mixture of the two known stable phases. Such a mixture could explain all spectral features observed in this investigation. The theoretical calculations provide reasonable agreement with the experimental data for most of the parameters, but predict H-bondings stronger than those observed. Differences between the spectra of the alpha and beta phases are predicted with similar characteristics to the experimental results.

Ices of CO2/H2O mixtures. Reflection-absorption IR spectroscopy and theoretical calculations.

Ice mixtures of CO2 and H2O are studied using Fourier transform reflection-absorption infrared (RAIR) spectroscopy. Mixtures are prepared by sequential deposition or co-deposition of the two components from the gas phase onto an Al plate kept at 87 K inside a low-pressure chamber. Two CO2 structures are found in most experiments: a crystalline form similar to pure CO2, which evaporates when warming at 105 K, and a noncrystalline species which remains embedded in amorphous water ice after warming. Significant spectral variations are found depending on the deposition method and the thickness of the solid. Features characteristic of the RAIR technique appear in the spectral regions of the normal modes of the bending and asymmetric stretching CO2 vibrations. Simulations using Fresnel theory and Mie scattering are carried out with acceptable agreement with the experimental spectra of solids of variable thickness, from approximately 1 microm to the limit of nanoparticles. Theoretical calculations of a pure CO2 crystal are performed. The relaxed geometry of the solid and its vibrational spectrum are determined and compared to the experimental results.

Orientation effects on nitric acid dihydrate films.

An investigation of orientation effects in films of nitric acid dihydrate (NAD) is presented, based on a systematic study of transmission and reflection-absorption infrared (RAIR) spectra of samples of varying thickness. The samples are prepared by vapor deposition on Ge (for transmission spectroscopy) and on Al substrates (for RAIR spectroscopy) at 175 K to produce crystalline alpha-NAD films. Transmission spectra were recorded at normal incidence, and RAIR spectra were recorded at a grazing angle of 75 degrees, with polarized radiation. The observed spectra are compared with predictions of a classical Fresnel model, to test the available optical indices of NAD, which are of great importance for the accurate interpretation of data from remote sensing measurements. Whereas the procedure yields satisfactory results for transmission and s-polarized RAIR spectra, it is found that the agreement is not acceptable for p-polarized RAIR spectra. An explanation is suggested in terms of a preferential alignment of the films, with the (10-1) crystallographic plane of the crystal situated parallel to the substrate. The infrared activity of a band at approximately 1170 cm(-1) is explained in terms of a preferential orientation of the crystal domains in the film.

On the use of wavelet filtering and correlation techniques in atmospheric condensed phase spectroscopy.

The application of wavelet filtering and analysis in spectroscopy is discussed in relation to the analysis of complex atmospheric spectra, where contributions from condensed phase particles and gas phase molecules are present in the form of broad-band features and narrow lines, respectively. The broad-band contribution can be extracted as the 'smooth signal' component of the wavelet transform, with a large reduction in the size of the corresponding data files. This procedure is applied to an investigation of the H2SO4 aerosol content of a series of atmospheric spectra measured in the ATMOS missions. The sulfate content of the smooth signal is analysed by means of correlation techniques, using a set of laboratory reference spectra of varying sulfuric acid concentration and temperature. Correlation density maps and correlation curves are used to select the most appropriate spectral zones for sulfate analysis and to assess the sulfate aerosol content in the atmosphere subsequent to the eruption of the Mount Pinatubo volcano.