Free Access
Issue
A&A
Volume 629, September 2019
Article Number A28
Number of page(s) 6
Section Astrophysical processes
DOI https://doi.org/10.1051/0004-6361/201834035
Published online 30 August 2019

© ESO 2019

1. Introduction

Methanol (CH3OH) is one of the most abundant organic molecules in interstellar clouds (e.g., Öberg et al. 2010) and has even been detected in extragalactic objects (e.g., Saito et al. 2015, 2016). However, its formation in the gas phase is highly inefficient (Geppert et al. 2006). Observations with the Infrared Space Observatory (ISO) have shown methanol on ice layers of dust grains (Ehrenfreund & Schutte 2000). This led to the proposal of formation mechanisms on the ice mantles of dust grains through hydrogenation of CO (see, e.g. Gibb et al. 2004 and references therein). It is assumed that collisions between molecules and dust grains in dense and cold regions of interstellar clouds lead to the deposition of gas-phase molecules on the grain surfaces, forming ice mantles (see, e.g. Zahorecz et al. 2017 and references therein). This initiates a rich chemistry that leads to the formation of complex organic molecules (COMs), including methanol.

Bacmann & Faure (2016) studied eight star-forming regions with different physical and chemical conditions. They found a relationship between the abundances of species such as HCO, CH3O, H2CO, and CH3OH. Formaldehyde (H2CO) was considered as a precursor of CH3OH (Bacmann & Faure 2016) on the ices. Several laboratory experiments that were conducted to predict the formation of COMs using CH3OH-rich ices with UV irradiation have been performed (e.g., Öberg et al. 2009). As a consequence of these experiments, a detailed chemical network has been established to explain the formation of COMs such as dimethyl ether on the ices. In addition to the ice chemistry, several reactions in the gas phase could lead to the formation of COMs. For example, Harada et al. (2017) have found that dimethyl ether CH3OCH3 can be produced in the gas phase at times shorter than 5 × 104 yr, while grain surface processes dominate for longer times. Gas-phase processes such as radiative associations have also been proposed to produce COMs from small radicals such as CH3O and CH3 (Balucani et al. 2015).

The discovery of COMs in cold dark clouds such as methyl formate, dimethyl ether (Cernicharo et al. 2012; see also Öberg et al. 2010; Bacmann et al. 2012), and methanol (Coutens et al. 2017) represents a challenge for astrochemical models as desorption mechanisms from the ices are not very efficient in these cold objects, which, in addition, are well protected against UV radiation. As a consequence, the observed abundances of COMs in cold environments cannot be reproduced by current chemical models (Vasyunin et al. 2017). In this context, theoretical and modeling efforts have been made by Vasyunin & Herbst (2013), Reboussin et al. (2014), Balucani et al. (2015), Chang & Herbst (2016), and Ruaud et al. (2015) to study the formation and abundances of COMs. The abundances reported, for example, by Ruaud et al. (2015) in their gas-grain models underestimate the observed values in cold dense cores (e.g., B1-b and L1689b) by more than two orders of magnitude.

The observations of Cernicharo et al. (2012) toward the cold core B1-b show methyl mercaptan (CH3SH), methoxy (CH3O), formic acid (HCOOH), propynal (HCCCHO), acetaldehyde (CH3CHO), dimethyl ether (CH3OCH3), methyl formate (CH3OCOH), and formyl radical (HCO). According to Cernicharo, CH3O is one of the key precursors of COMs. The formation of this radical, either in the ices or in the gas phase, could therefore be the first step to initiate the chemistry of COMs (Öberg et al. 2009; Antiñolo et al. 2016). In the gas phase, methoxy could be formed from the photodissociation of methanol or by the reaction of methanol with the radical OH. To clarify the importance of methanol as a generating species, laboratory experiments using the CRESU technique have been performed (Antiñolo et al. 2016) to measure the reaction CH3OH + OH at temperatures as low as 22 K. These experiments have shown that the reaction rate increases considerably at low temperatures, with the methoxy CH3O and hydroxymethyl CH2OH radical as the most likely products. When methoxy is available in the gas phase, its radiative association with CH3 could produce dimethyl ether (Balucani et al. 2015), and its reaction with other species could lead to more complex species.

Highly energetic particles could have a significant effect on the chemistry of interstellar ices (Shingledecker & Herbst 2018), especially in cold dense molecular cloud cores, where thermal energy is not sufficient to drive the normal diffusive chemistry of grain species (Ruaud et al. 2015). In order to provide theoretical support for the formation of COMs in ices through the reaction of methanol with other species, we have undertaken a systematic ab initio study of the system (CH3OH)10 ice-mantle plus OH, OH, and OH+. We here present the first part results of ab initio molecular dynamics (AIMD) simulations of the irradiation of an amorphous methanol-ice mantle that is impacted by OH molecules. The results of OH radical and OH+ cation will be published in a second article, in which we analyze the charge effect and COM formation processes as well. The OH radical, which is one of the most abundant molecules in the ISM, can easily form its anion by the OH + e → OH + 1.83 eV process, leading to the stable anion OH. Its rotational and vibrational spectrum has been measured in the laboratory (Liu & Oka 1986; Owrutsky et al. 1985), but it has not been detected in space so far. Even if OH has a low abundance in the ISM, it is detected in other astrophysical environments such as the coma of comets series (Chaizy et al. 1991). From the theoretical point of view, the OH anion is extremely interesting because it is isoelectronic with H2O and is a closed-shell species. It can therefore be considered as a prototype for the interaction of the water molecule in interstellar and cometary ices (Cordiner & Charnley 2014), and moreover, it can relatively easily be modeled by calculations based on first principles.

In our study of the reaction of methanol under the impact of OH projectile, we use the Born-Oppenheimer ab initio molecular dynamics to identify which outcomes emerge depending on conditions such as the kinetic energy of the OH projectile and the impact parameters of the collision. We also comment in detail on how the outcome could be linked to the formation of COMs. The remainder of this paper is organized as follows. In Sect. 2 (computational details) we describe the calculation model in more detail. The main outcomes of our simulations are presented in Sect. 3 and discussed in Sect. 4. We present our conclusions in Sect. 5.

2. Computational details

We simulated a methanol ice mantle as a cluster with ten units of methanol (hereafter (CH3OH)10-ice-mantle) in collision with OH as a projectile. According to Boyd & Boyd (2007), this cluster corresponds to the most stable isomer. We chose this cluster because its size is a good compromise between a good representation of a real dust ice mantle and the computational cost of the AIMD we used to simulate the impact of OH. The collision was set to occur in the ground state of the species. This is a good approximation because the lowest excitation energy of methanol (at least 6 eV) (Varela et al. 2015) and hydroxyl (at least 4 eV) (Schofield & Kjaergaard 2004) is a significant fraction of the kinetic energy of the projectile we used during the simulations (10 to 22 eV), which implies a low probability of nonadiabatic dynamics. Therefore, we used Born-Oppenheimer molecular dynamics (BOMD) at constant energy. This microcanonical ensemble approximation reproduces the conditions of the ISM because the low particle density inhibits the dissipation of thermal energy on timescales similar to the collision (< 400 fs).

Forces on atoms and total energy were calculated using DFT within the Kohn-Sham (KS) formalism. The KS equations are a set of mono-electronic Schrödinger equations with a mean field for electrons that can recover the density of the actual interacting system. Although the KS equations are formally exact, in practice, only approximations to the exchange-correlation potential of the mean-field are known. Here we used the long range-corrected hybrid functional of Head-Gordon et al., ωB97X-D (Helgaker et al. 1990; Uggerud & Helgaker 1992; Bolton et al. 1998). KS orbitals were expanded in a basis set of Gaussian functions centered in each atom; namely, the Pople 6-31+G(d,p) basis set (Rassolov et al. 2001). This ωB97X-D basis set is known to be flexible enough to describe reactive collisions (McBride et al. 2013). A Verlet-velocity algorithm with a time step of 0.5 fs was used to integrate the equations of motion. This time step guarantees that the energy is conserved within a few hundreds of kcal mol−1 during the whole simulation. All calculations were made using the electronic structure package Gaussian 09 (Frisch et al. 2009).

In all collisions the oxygen atom of OH pointed to the center of mass of the (CH3OH)10-ice-mantle while the initial kinetic energy of the (CH3OH)10-ice-mantle was set to the equivalent of 30 K (Dawes et al. 2007). To find the range of energy of the projectile in which a rich chemistry occurs, we preliminary scanned impact energies from 5 to 40 eV. This range corresponds to an estimation of the limits of an elastic collision of the projectile with the cluster and an ultrafast penetration of the projectile with no time for reactions to occur (Kohanoff & Artacho 2008).

Impacts below 10 eV result in adsorption of OH on the surface of the (CH3OH)10-ice-mantle, while impacts with kinetic energy above 22 eV rapidly evaporate the mantle and breaking the methanol ice mantle structure. Consequently, we focused only on collisions with initial OH kinetic energy of 10, 12, 15, 18, 20, and 22 eV, where the richest chemistry takes place. For each impact energy, 24 trajectories with a total physical time of 400 fs were simulated. These trajectories start with the velocity of H–O axis of OH projectile pointing toward the center of mass of the (CH3OH)10-ice-mantle. The 24 starting points were evenly distributed over a sphere, as depicted in Fig. 1, to obtain diversity in the impact factors.

thumbnail Fig. 1.

Schematic representation of the initial conditions used for the molecular dynamics of the collision of OH and a methanol ice-mantle. Each blue dot represents an initial position of the OH projectile, whose O–H axis points toward the center of mass of the methanol cluster (black circle) with its oxygen atom facing the cluster. The 24 points on the surface are evenly distributed such that the solid angle subtended by three closest-neighbor points on the surface and the center of the sphere is π 5 $ \frac{\pi}{5} $ sr.

Each trajectory was visually inspected to search for well-formed intermediates or products. By well formed, we mean molecules that keep their structure for five or more times the period of the important vibrational mode involved in the reaction, that is, at least 60 fs. For instance, if an intermediate is formed through a hydrogen transfer to oxygen, then this newly formed bond should oscillate at least five times before the molecule follows its way along a new reaction path. Finally, when needed, the local population analysis of Mulliken (Mulliken 1955) was used to verify the charge of species involved in the reaction pathways. To simplify the nomenclature in Table 1, we indicate CH3OH + OH → products instead of (CH3OH)10 + OH → products + (CH3OH)x.

Table 1.

Main products obtained from (CH3OH)10-ice-mantle in collision with a OH projectile with kinetic energy of 10–22 eV.

3. Results

In the set of simulations using OH as a projectile, different COMs are formed, such as methylene glycol (CH2(OH)2) and radicals such as OCH2OH, which have not been detected in the ISM. Some of the other molecules produced here have already been detected in the ISM, for example, the methyl radical ( CH 3 $ {\rm{C}}{{\rm{H}}_3}^ \bullet $, methoxy radical CH3O) (Cernicharo et al. 2012), dimethyl ether (CH3OCH3) (Snyder et al. 1974), formaldehyde (H2CO) (Snyder et al. 1969), and formaldehyde protonated (H2COH+) (Ohishi et al. 1996).

As mentioned above, we used projectile kinetic energy ranges from 10 eV to 22 eV. Below 10 eV, the (CH3OH)10-ice-mantle only warms up, which reduces the number of hydrogen bonds, and deforms the structure, and ejects some methanol molecules into the gas phase (desorption). Above 22 eV, all molecules are ejected into the gas phase, which destroys the methanol ice mantle and breaks the CH3OH molecules.

Figure 2 shows a general overview of the BOMD results obtained between (CH3OH)10-ice-mantle and OH projectile. We analyze the results of the formation and destruction processes observed in the simulations below at impact energies of 10 and 12 eV in Sect. 3.1, followed by reactions at impact energies of 15 to 22 eV in Sect. 3.1.

thumbnail Fig. 2.

Chemical network of COMs obtained after bombardment of the (CH3OH)10-ice-mantle with OH using impact energies of 10, 12, 15, 18, 20, and 22 eV. The (CH3OH)10-ice-mantle is represented as CH3OH hereafter, and radicals are not indicated by dots for clarity.

3.1. Reactions at impact energies of 10 and 12 eV

Impacts1 with projectile energy of 10 eV always yield one main product that is listed as reaction (1) in Table 1. Figure 3 plots the potential energy as a function of time since impact in reaction (1). After 40 fs since the collision, a potential energy barrier is observed that is associated with the formation of a methoxy CH3O. After a time of 100 fs, no new chemical bonds are created or destroyed. Thus, only a primary process was observed.

thumbnail Fig. 3.

Potential energy as a function of time since impact for all BOMD simulations with projectile at 10 eV. After a time of 100 fs no new chemical bonds are created or destroyed.

The OH-group in the (CH3OH)10-ice-mantle is not reactive because of the hydrogen bridges. Thus, it is not easy to extract a proton from it. After the impact with OH, the substrate structure absorbs part of the kinetic energy to break the hydrogen bridges, which causes the hydroxyl group to react more easily. The breaking of these hydrogen bonds makes a hydroxyl available and tends to form CH3O instead of CH2OH. This behavior is identical in all simulations at 10 eV.

Impacts with a kinetic energy of 12 eV also most often lead to the formation of CH3O through reaction (1). However, three additional chemical pathways (reactions (2)−(4)) can also occur (see Table 1).

The CH3 ⋅ and OH ⋅ species are formed in reaction (2) through direct impact of the hydroxyl ion on the C–O bond in methanol. They can escape from the mantle with no further reactions.

Reaction (3) in Table 1 produces the formation of methylene glycol (or diol) CH2(OH)2, which has not been detected in the ISM. At higher impact energies, this species can undergo a secondary process, the dehydration to produce formaldehyde (H2CO) and water; both are abundant in the ISM (Bacmann & Faure 2016). This reaction requires the elimination of a proton via nucleophilic substitution (SN2 mechanism). The process takes about 50 fs.

Reaction (4) is a secondary process that starts after the products obtained from reaction (3) are stabilized. This mechanism exhibits the formation of H2. The BOMD simulation shows that the two hydrogen atoms hover over the grain surface until structural conditions are favorable for H2 association.

3.2. Reactions at an impact energy of 15 to 22 eV

When the kinetic energy of the projectile is 15 eV or higher, many more reaction paths are possible. The reactions occur far from equilibrium; thus, the energy profile gives little information on reaction barriers. The outcomes of the reactions are listed in Table 1. The last column lists the impact energies in which the reactions were produced. The first four reactions are the same as in the 12 eV kinetic energy case. The primary outcome between 15 eV and 22 eV is still CH3O. The second most probable product formed is methylene glycol (CH2(OH)2) at energies between 15 eV and 22 eV, as illustrated by reaction (3).

Using an impact energy of 15 eV on the (CH3OH)10-ice-mantle corresponds to an average energy of 1.36 eV per molecule in the system including the projectile. This energy is sufficient to produce a variety of reactions. We find 18 reaction pathways and report them in Table 1. Reactions (1)−(4) are the same as were obtained with a 12 eV projectile, but their formation time decreases because the impact energy is higher at 15 eV.

The hydroxymethyl radical CH2OH was obtained with impacts energies of 18 and 20 eV (reaction (12)). These products have a common origin: methylene glycol CH2(OH)2. This precursor can also lead through a secondary process (reactions (6) and (7)) to a new product the OCH2OH anion, which has not been identified in the ISM. In the first case (reaction (6)), the diol reacts with the methoxy CH3O to yield OCH2OH. In reaction (7), the diol reacts with H to form OCH2OH + H2.

Reaction (8) shows the formation of dimethyl ether CH3OCH3 after impacts of 15 eV. The mantle absorbs most of the kinetic energy as thermal energy. In this environment, the CH3O and CH 3 $ {\rm{C}}{{\rm{H}}_3}^ \bullet $ radicals are transient species with lifetimes between 45 and 164 fs. After this time, radicals trapped in the (CH3OH)10-ice-mantle diffuse and recombine to form CH3OCH3 + H2O + OH, which are then ejected into the gas phase. Thus, secondary processes may contribute to the formation of COMs in astrophysical environments. At an impact energy of 18 eV, reaction (11) is favored instead.

The formation mechanism of CH3OCH3 can be divided into three steps:

  • (1)

    a homolytic rupture of CH3OH to produce CH 3 $ {\rm{C}}{{\rm{H}}_3}^ \bullet $ + OH;

  • (2)

    a proton abstraction from CH3OH to form water and a CH3O radical;

  • (3)

    radical-radical reaction to produce CH3OCH3. These steps are illustrated in Fig. 4.

    thumbnail Fig. 4.

    Formation of dimethyl ether CH3OCH3 at 15 eV following the chemical pathways as a function of time. a: initial conditions. In color the projectile and methanol molecules that actively react. The remaning methanol molecules are represented by grey lines. They absorve some of the collision energy even though they do not break or form new chemical bonds. b: radical formation at 87.5 fs. c: radical-radical association at 164 fs. The colored OH and H2O molecules escape the system.

At an impact energy of 18 eV, we found the HCOH isomer of H2CO. Its formation is a key factor of stabilization of its precursors CH3O + CH2OH to form two stable chemical species: CH3OH and HCOH in a “chair” structure. The HCOH appears in reactions (13) and (14) (see Table 1) with impact energies of 18 to 22 eV. The formation mechanism for this molecule can be described in the following three steps:

  • (1)

    diol formation: CH2(OH)2 + H;

  • (2)

    fragmentation: CH2(OH)2 → CH2OH + H2O + CH 3 $ {\rm{C}}{{\rm{H}}_3}^ \bullet $;

  • (3)

    proton abstraction: CH2OH → HCOH + H2O + CH3OH ⋅

Finally, formaldehyde H2CO was formed through a secondary process (reaction (15)) after impact energies of 18, 20, and 22 eV. The formation mechanism can be described in three steps:

  • (1)

    deprotonation of methyl-group;

  • (2)

    proton abstraction from the OH-group;

  • (3)

    stabilization through carbon-oxygen double bonds to form H2CO + H2O + H.

4. Discussion

The simulation of reactive collisions of an energetic OH projectile with an amorphous (CH3OH)10-ice mantle produces a rich diversity of COMS, including methoxy CH3O, and hydroxymethyl CH2OH radicals and methylene glycol CH2(OH)2 as primary products. The CH3O radical is the most frequent product at all energies, and the only product found at 10 eV. In this way, our results are in agreement with the gas-phase results reported by Shannon et al. (2013), which also led to the CH3O radical as the main product. If the reactants are correctly oriented, the reaction is expected to proceed by abstraction of a hydrogen atom at either the methyl or the hydroxyl site to produce CH3O or CH2OH radicals, respectively, which in most cases produces the CH3O radical over the CH2OH radical. Nevertheless, Antiñolo et al. (2016) pointed out that the role of grain-surface cannot be neglected in the formation of the CH2OH radical.

At collision energies higher than 10 eV, these COMs can participate in secondary reactions that generate more COMs. Thus, through these processes, relevant isomers such as methylene glycol and dimethyl ether are chemically methanol dependent.

Dimethyl ether CH3OCH3 has not been detected as a product of CH3OH ice reactions in previous experimental studies (Gerakines et al. 1996; Bennett et al. 2007; Baud & Wouterloot 1980). However, in our simulations, CH3OCH3 is formed by secondary processes involving CH 3 + CH 3 O $ \mathrm{CH}_{3}{^{\bullet}} + \mathrm{CH}_{3}\mathrm{O}^{\bullet} $.

The diol CH2(OH)2 forms by a one-step concerted mechanism through the addition of OH and proton abstraction at impact energies above 15 eV. It is the second most frequent product observed.

Formaldehyde is created in a new chemical reaction route starting from methanol ice mantles. This can explain why the hydroxymethyl radical has not been detected in astronomical observations. We find that when the hydroxymethyl CH2OH is formed, the radical will undergo proton abstraction to produce H2CO.

These simulations provide new possible scenarios for the formation of COMs in interstellar environments. We showed that molecular projectiles with adequate kinetic energy impacting on cold ice-covered dust grains can produce rich COMs that are released into the gas phase. The required kinetic energies of 10−22 eV are too high to be representative of typical conditions in molecular clouds, with typical gas kinetic temperatures of 8−20 K (Bergin & Tafalla 2007). However, these projectiles can correctly represent conditions in localized regions such as photodissociation regions at the interface of HII regions and their parent molecular clouds, or in shocked regions due to the impact of high-velocity protostellar jets and outflows in the parent dense molecular cloud core. For example, specific locations in the protostellar outflow source L1157 are known to have a rich chemistry including abundant COMs (e.g., Bachiller et al. 2001; Arce et al. 2008). Similar processes are likely to be prevalent in regions of high-mass star formation, as shown recently by Guzmán et al. (2018). Guzmán et al. (2014) suggested that photodesorption may be an efficient mechanism to release complex molecules into the gas phase in photodissociation regions. We showed that moderate energy impacts with molecular projectiles can contribute not only to the desorption of COMS into the gas, but also be a significant factor in the creation of COMs themselves in ice mantles on dust grains.

5. Conclusions

Chemical reactions on icy grain mantles that subsequently eject products from the surface into the gas phase play a key role in the synthesis of COMs in astrophysical sources. We simulated the detailed interaction of an amorphous cold methanol decamer ((CH3OH)10) bombarded by OH projectiles with kinetic energies of 10−22 eV. We found that many complex organic molecules are produced in the collision, including CH3O, CH2(OH)2, CH2O, OCH2OH, CH3OCH3, CH2OH, HCOH, and H2CO.

The reaction CH3OH + OH → CH3O + H2O is most frequently observed in this work and thus the most probable to occur every time that methanol reacts with the hydroxyl anion. The projectile kinetic energy is crucial in breaking the hydrogen bridge and the subsequent modification of the ice-mantle structure. Different chemical products were observed as a function of the kinetic energy on impact. At high impact energies, the formation of CH2(OH)2 was found to be more probable than CH3O. At higher impact energies, at least 17 reaction pathways are found, as listed in Table 1.

We confirmed that H2CO is a product of the reaction of methanol with a hydroxyl anion. This is consistent with the earlier suggestion of a chemical relation between the fragments of CO, HCO, H2CO, CH3O, and CH3OH through an hydrogenation process (Bacmann & Faure 2016, and references therein). Furthermore, we showed that methylene glycol CH2(OH)2 is an outcome of primary process and that HCOH and OCH2OH are formed in secondary processes involving CH2(OH)2. We find a new mechanism to produce H2CO every time CH2OH is formed.

Given the projectile kinetic energies required for these reactions to occur, we suggest that these processes are likely relevant in the production of COMs in photodissociation regions and in shocks produced by high-velocity jets and outflows from young stellar objects.


1

Figure 1 shows the impact configurations.

Acknowledgments

This research was supported by FONDECYT-Chile grant 11140770 and CONICYT International Networks for young researchers Grant REDI170243. NI gratefully acknowledges support of Vicerectoría de Investigación y Postgrado (VRIP). DM acknowledges support from CONICYT project Basal AFB170002. CC and PF acknowledge support form Centers Of Excellence With Basal/Conicyt Financing, Grant FB0807 and FONDECYT-Chile grant 1181121, and 1180623. JG acknowledges support from Fondecyt postdoctoral fellowship 3170768. We thank Desmond MacLeod-Carey, Paola Caselli, Gil Jannes, and Christopher Shingledecker for their helpful comments and suggestions.

References

  1. Antiñolo, M., Agúndez, M., Jiménez, E., et al. 2016, ApJ, 823, 25 [NASA ADS] [CrossRef] [Google Scholar]
  2. Arce, H. G., Santiago-García, J., Jørgensen, J. K., Tafalla, M., & Bachiller, R. 2008, ApJ, 681, L21 [NASA ADS] [CrossRef] [Google Scholar]
  3. Bachiller, R., Pérez Gutiérrez, M., Kumar, M. S. N., & Tafalla, M. 2001, A&A, 372, 899 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Bacmann, A., & Faure, A. 2016, A&A, 587, A130 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Bacmann, A., Taquet, V., Faure, A., Kahane, C., & Ceccarelli, C. 2012, A&A, 541, L12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Balucani, N., Ceccarelli, C., & Taquet, V. 2015, MNRAS, 449, L16 [Google Scholar]
  7. Baud, B., & Wouterloot, J. G. A. 1980, A&A, 90, 297 [NASA ADS] [Google Scholar]
  8. Bennett, C., Chen, S.-H., Sun, B.-J., Chang, A., & Kaiser, R. 2007, ApJ, 660, 1588 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bergin, E. A., & Tafalla, M. 2007, ARA&A, 45, 339 [NASA ADS] [CrossRef] [Google Scholar]
  10. Bolton, K., Hase, W. L., & Peslherbe, G. H. 1998, Modern Methods for Multidimensional Dynamics Computations in Chemistry (Singapore: World Scientific), 143 [CrossRef] [Google Scholar]
  11. Boyd, S. L., & Boyd, R. J. 2007, J. Chem. Theory Comput., 3, 54 [CrossRef] [Google Scholar]
  12. Cernicharo, J., Marcelino, N., Roueff, E., et al. 2012, ApJ, 759, L43 [NASA ADS] [CrossRef] [Google Scholar]
  13. Chaizy, P., Reme, H., Sauvaud, J., et al. 1991, Nature, 349, 393 [NASA ADS] [CrossRef] [Google Scholar]
  14. Chang, Q., & Herbst, E. 2016, ApJ, 819, 145 [NASA ADS] [CrossRef] [Google Scholar]
  15. Cordiner, M. A., & Charnley, S. B. 2014, Meteorit. Planet. Sci., 49, 21 [NASA ADS] [CrossRef] [Google Scholar]
  16. Coutens, A., Rawlings, J., Viti, S., & Williams, D. A. 2017, MNRAS, 467, 737 [NASA ADS] [Google Scholar]
  17. Dawes, A., Hunniford, A., Holtom, P. D., et al. 2007, Phys. Chem. Chem. Phys., 9, 2886 [NASA ADS] [CrossRef] [Google Scholar]
  18. Ehrenfreund, P., & Schutte, W. A. 2000, Adv. Space Res., 25, 2177 [NASA ADS] [CrossRef] [Google Scholar]
  19. Frisch, M. J., Trucks, G. W., Schlegel, H. B., et al. 2009, Gaussian 09 Revision A.1 (Wallingford, CT: Gaussian Inc.) [Google Scholar]
  20. Geppert, W. D., Hamberg, M., Thomas, R. D., et al. 2006, Faraday Discuss., 133, 177 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  21. Gerakines, P., Schutte, W., & Ehrenfreund, P. 1996, A&A, 312, 289 [NASA ADS] [Google Scholar]
  22. Gibb, E., Whittet, D., Boogert, A., & Tielens, A. 2004, ApJS, 151, 35 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  23. Guzmán, V. V., Pety, J., Gratier, P., et al. 2014, Faraday Discuss., 168, 103 [NASA ADS] [CrossRef] [Google Scholar]
  24. Guzmán, A. E., Guzmán, V. V., Garay, G., Bronfman, L., & Hechenleitner, F. 2018, ApJS, 236, 45 [NASA ADS] [CrossRef] [Google Scholar]
  25. Harada, N., Hasegawa, Y., Aikawa, Y., et al. 2017, ApJ, 837, 78 [NASA ADS] [CrossRef] [Google Scholar]
  26. Helgaker, T., Uggerud, E., & Jensen, H. J. A. 1990, Chem. Phys. Lett, 173, 145 [NASA ADS] [CrossRef] [Google Scholar]
  27. Kohanoff, J., & Artacho, E. 2008, AIP Conf. Proc., 1080, 78 [NASA ADS] [CrossRef] [Google Scholar]
  28. Liu, D.-J., & Oka, T. 1986, J. Chem. Phys., 84, 2426 [NASA ADS] [CrossRef] [Google Scholar]
  29. McBride, E. J., Millar, T. J., & Kohanoff, J. J. 2013, J. Phys. Chem. A, 117, 9666 [Google Scholar]
  30. Mulliken, R. S. 1955, J. Chem. Phys., 23, 1833 [NASA ADS] [CrossRef] [Google Scholar]
  31. Öberg, K. I., Garrod, R. T., van Dishoeck, E. F., & Linnartz, H. 2009, A&A, 504, 891 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Öberg, K. I., Bottinelli, S., Jørgensen, J. K., & van Dishoeck, E. F. 2010, ApJ, 716, 825 [NASA ADS] [CrossRef] [Google Scholar]
  33. Ohishi, M., Ishikawa, S.-I., Amano, T., et al. 1996, ApJ, 471, L61 [NASA ADS] [CrossRef] [Google Scholar]
  34. Owrutsky, J. C., Rosenbaum, N. H., Tack, L. M., & Saykally, R. J. 1985, J. Chem. Phys., 83, 5338 [NASA ADS] [CrossRef] [Google Scholar]
  35. Rassolov, V. A., Ratner, M. A., Pople, J. A., Redfern, P. C., & Curtiss, L. A. 2001, J. Comput. Chem., 22, 976 [CrossRef] [Google Scholar]
  36. Reboussin, L., Wakelam, V., Guilloteau, S., & Hersant, F. 2014, MNRAS, 440, 3557 [NASA ADS] [CrossRef] [Google Scholar]
  37. Ruaud, M., Loison, J. C., Hickson, K. M., et al. 2015, MNRAS, 447, 4004 [NASA ADS] [CrossRef] [Google Scholar]
  38. Saito, T., Iono, D., Yun, M. S., et al. 2015, ApJ, 803, 60 [NASA ADS] [CrossRef] [Google Scholar]
  39. Saito, T., Iono, D., Espada, D., et al. 2016, ApJ, 834, 6 [NASA ADS] [CrossRef] [Google Scholar]
  40. Schofield, D. P., & Kjaergaard, H. G. 2004, J. Chem. Phys., 120, 6930 [NASA ADS] [CrossRef] [Google Scholar]
  41. Shannon, R. J., Blitz, M. A., Goddard, A., & Heard, D. E. 2013, Nat. Chem., 5, 745 [NASA ADS] [CrossRef] [Google Scholar]
  42. Shingledecker, C. N., & Herbst, E. 2018, Phys. Chem. Chem. Phys., 20, 5359 [Google Scholar]
  43. Snyder, L. E., Buhl, D., Zuckerman, B., & Palmer, P. 1969, Phys. Rev. Lett., 22, 679 [NASA ADS] [CrossRef] [Google Scholar]
  44. Snyder, L., Buhl, D., Schwartz, P., et al. 1974, ApJ, 191, L79 [NASA ADS] [CrossRef] [Google Scholar]
  45. Uggerud, E., & Helgaker, T. 1992, J. Am. Chem. Soc., 114, 4265 [CrossRef] [Google Scholar]
  46. Varela, K., Hargreaves, L. R., Ralphs, K., et al. 2015, J. Phys. B: At. Mol. Opt. Phys., 48, 115208 [CrossRef] [Google Scholar]
  47. Vasyunin, A., & Herbst, E. 2013, ApJ, 769, 34 [Google Scholar]
  48. Vasyunin, A. I., Caselli, P., Dulieu, F., & Jiménez-Serra, I. 2017, ApJ, 842, 33 [Google Scholar]
  49. Zahorecz, S., Jimenez-Serra, I., Testi, L., et al. 2017, A&A, 602, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Tables

Table 1.

Main products obtained from (CH3OH)10-ice-mantle in collision with a OH projectile with kinetic energy of 10–22 eV.

All Figures

thumbnail Fig. 1.

Schematic representation of the initial conditions used for the molecular dynamics of the collision of OH and a methanol ice-mantle. Each blue dot represents an initial position of the OH projectile, whose O–H axis points toward the center of mass of the methanol cluster (black circle) with its oxygen atom facing the cluster. The 24 points on the surface are evenly distributed such that the solid angle subtended by three closest-neighbor points on the surface and the center of the sphere is π 5 $ \frac{\pi}{5} $ sr.

In the text
thumbnail Fig. 2.

Chemical network of COMs obtained after bombardment of the (CH3OH)10-ice-mantle with OH using impact energies of 10, 12, 15, 18, 20, and 22 eV. The (CH3OH)10-ice-mantle is represented as CH3OH hereafter, and radicals are not indicated by dots for clarity.

In the text
thumbnail Fig. 3.

Potential energy as a function of time since impact for all BOMD simulations with projectile at 10 eV. After a time of 100 fs no new chemical bonds are created or destroyed.

In the text
thumbnail Fig. 4.

Formation of dimethyl ether CH3OCH3 at 15 eV following the chemical pathways as a function of time. a: initial conditions. In color the projectile and methanol molecules that actively react. The remaning methanol molecules are represented by grey lines. They absorve some of the collision energy even though they do not break or form new chemical bonds. b: radical formation at 87.5 fs. c: radical-radical association at 164 fs. The colored OH and H2O molecules escape the system.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.