Abstract
Anomalous microwave emission (AME) is detected in many astrophysical environments as a foreground feature typically peaking between 20–30 GHz and extending over a 10–60 GHz range. One of the leading candidates for the source of AME is small spinning dust grains. Such grains should be very small (approx. ≤1 nm diameter) in order for the rotational emission to fall within the observed frequency range. In addition, these nanosized grains should possess a significant dipole moment to account for the observed emissivities. These constraints have been shown to be compatible with spinning bare nanosilicate clusters, assuming that ∼1% of the total Si mass budget is held in these ultrasmall grains. Silicate dust can be hydroxylated by processing in the interstellar medium and is generally known to provide seeds for molecular water ice nucleation in denser regions. Herein, we use quantum chemical calculations to investigate how the dipole moment of Mg-rich pyroxenic (MgSiO3) nanoclusters is affected by both accretion of molecular water and dissociative hydration. Our work thus provides an indication of how the formation of water ice mantles is likely to affect the capacity of nanosilicates to generate AME.
Introduction
Silicate dust is ubiquitously found in a wide range of astrophysical environments (). Chemical nucleation and growth in the circumstellar envelopes of oxygen-rich evolved stars is an important source of silicate dust (; ). From infrared (IR) emission observations, it is inferred that most of these newly formed silicate grains are Mg-rich and have sizes of the order 0.1 μm (diameter) (). Interactions between the dust and the stellar radiation field cause the grains to be ejected into the interstellar medium (ISM). Supernova shocks in the ISM induce grain–grain collisions which shatter the dust into smaller components (). A typical shock is thought to result in 5–15% of the original grain mass ending up in the form of ultrasmall nanosilicates (<3 nm diameter). Nanograin production from larger grains is also predicted to occur close to sources of intense radiation (e.g., massive stars and supernovae) via rotational disruption of dust grains by radiative torques (). These models of dust processing are compatible with estimates that approximately 10% of the Si in the ISM could be held in ≤3 nm diameter nanosilicate grains () without violating observational constraints.
Via gas–grain collisions and interactions with radiation, small nanosilicate grains will be rotationally excited. Such spinning nanosilicate grains have been proposed as the main carrier of anomalous microwave emission (AME) – a microwave foreground feature which typically peaks between 20–30 GHz and extends over a 10–60 GHz range (; ; ). AME is observed in a number of astrophysical environments including molecular clouds () and the diffuse ISM (). As smaller grains can be spun up to larger frequencies than larger grains in any given astronomical environment, proportionally they emit more power, and thus the smallest grains are expected to dominate the spinning dust microwave spectrum. For a grain to emit in the microwave region associated with AME via rotational emission, its size should be of the order of ∼1 nm in diameter (; ). In principle, contributions to the AME from the spinning of such nanograins could occur from any type of nanograin if it is sufficiently abundant and has a large enough dipole moment (; ). Observational evidence has tended to rule out small spinning organic dust species (i.e., polycyclic aromatic hydrocarbons—PAHs) as the most important source of AME (). The atomistically detailed structures and physical properties of nanosilicates with Mg-rich olivinic (Mg2SiO4)N and pyroxenic (MgSiO3)N compositions with diameters ≤1 nm have recently been determined by quantum chemical modeling (). These accurate calculations clearly show how the ionic chemical bonding in nanosilicates intrinsically tends to lead to grains with significant dipoles. Based on those structures, it has been shown that the AME can be explained if only 1% of the Si in the ISM were held in ultrasmall (i.e., ≤1 nm diameter) nanosilicate grains (). Although the dipoles possessed by these bare nanosilicates are more than sufficient to account for the observed AME, it is not clear if these results can be extended to nanosilicates, which become hydroxylated by processing and/or upon which water ice forms.
Silicate grains emerging from the outflows of evolved stars are processed from ion irradiation (e.g., protons, He+) from supernovae shocks in the ISM. From laboratory studies of the effect of the solar wind on lunar surface silicate material, proton irradiation was found to cause two important intertwined chemical processes: 1) hydroxyl formation () and 2) Mg cation depletion by preferential sputtering (). Analysis of interplanetary dust has provided direct evidence of such grain processing (). Such processes would tend to transform Mg-rich nano-olivine dust grains into hydroxylated nano-pyroxene grains. Relatively elevated levels of pyroxene (with respect to the initial dominance of olivine) are consistent with the pyroxene-rich environments of protostars () and the estimated fraction of amorphous pyroxene grains in the ISM ().
In denser regions of the ISM such as molecular clouds, the ultraviolet radiation field is shielded and molecular water can accrete on silicate grain surfaces. This marks the beginning of the formation of water ice–dominated mantles on silicate grains. Water is one of the most abundant molecules in the universe in the solid phase (; ). Interstellar water ices are thought to present an amorphous structural state because of the similarities of their IR features with those of amorphous solid water measured in terrestrial laboratories (). Gas-phase reactions based on low-temperature ion-neutral and high-temperature neutral-neutral processes allow for the explaining of the abundance of water as a gaseous component, but they cannot justify the large abundance of water ice in interstellar and protostellar clouds (; ). Accordingly, reaction routes based on surface chemical processes have long been recognized to be the major channel for the formation of water ice (; ). Therefore, water ice, rather than being formed by successive adsorption of water molecules coming from the gas phase, is formed in situ on the surfaces of dust grains. Different reactive channels leading to the formation of interstellar water on the surfaces of grains are hydrogenation of atomic O (e.g., ), of molecular O2 (e.g., ), and even of O3 (e.g., ), among others. We refer to the work of ) for a complete list of postulated reactions leading to water formation on interstellar grains.
Herein, we perform a systematic quantum chemical study of the effect of both hydroxylation and molecular water interactions on ultrasmall (∼1 nm diameter) Mg-rich pyroxenic nanosilicates with a (MgSiO3)6 composition. We specifically focus on how the dipole moments of the original bare P-6 nanosilicates are affected by these chemical changes. If we focus on the dipole moments of the most energetically stable grains in a population, we find that both processes can significantly affect the dipole moment of these selected grains. However, when considering a more realistic population of nanosilicate grains, our results suggest that neither hydroxylation nor molecular hydration has a significant effect on the average magnitude of the dipole moment. Our results thus suggest that the capacity of nanosilicates to be the source of AME should not be significantly diminished by processing and/or initial water ice mantle formation.
Methods
The initial bare (MgSiO3)6 isomer structures were obtained in a previous work (). We chose to compare two distinct isomers: 1) the lowest energy (global minimum) cluster structure with a very low dipole and 2) a slightly energetically metastable cluster structure with a high dipole. These bare (MgSiO3)6 structures are shown in Figure 1.
FIGURE 1
In the case of hydroxylated bare nanosilicate clusters, we carried out global optimization calculations to find the 50 lower energy isomers of (MgSiO3)6 after dissociative interaction of one, two, and three water molecules, respectively. The global optimization calculations were split into two steps. First, we used the Monte Carlo basin hopping (MCBH) algorithm (), employing previously developed polarizable classical interatomic potentials (IPs) (), to explore the potential energy surface (PES), followed by the structural configurations of the hydroxylated clusters. Previous studies have shown that the MCBH method is an efficient method for exploring the low energy PES of nanoclusters and for finding the global minimum (). The standard MCBH algorithm uses random atomic displacements at each step to explore the PES. In addition, we also allowed that 0.5% of the steps attempt an Mg ↔ Si cation exchange move assist in this exploration. To run the MCBH calculations, we employed a previously developed cascade MCBH code (), which uses the General Utility Lattice Package (GULP) () to locally optimize the cluster structures at each step. To avoid problems with the local optimization of distorted structures containing polarized ions, we preoptimize each cluster structure with an IP with polarizable shells removed before optimizing using the IP with shells. We note that, in principle, an MCBH global optimization search only depends on the chemical composition of the system and should be independent of the initial hydroxylated cluster structure. In practice, the exploration of the PES can be somewhat biased by the initial structure and thus, to help ensure that our searches were not biased, for each addition of water, we ran the MCBH searches a total of 20 times with different starting structures. Each MCBH run consisted of 10,000 steps.
After this global optimization procedure, we selected 100 distinct lowest energy nanocluster structures for each degree of hydroxylation for further refinement. The selected cluster structures were optimized using density functional theory (DFT)–based calculations using the PBE functional () with a light-tier one numerical atom-centered orbital (NAO) basis set. Previous studies on oxide nanoparticles have shown that this basis set has a similar quality to valence triple-ζ plus polarization Gaussian-type Orbital basis sets (). These DFT-based calculations employed the Fritz Haber Institute Ab Initio Molecular Simulation package (FHI-AIMS) ().
In the case of the molecularly hydrated (MgSiO3)6 clusters, the initial positions of the physisorbed water molecules on both bare cluster isomers (Figure 1) were prepared using an electrostatic compatibility criterion. Specifically, the negatively charged oxygen atom of each water molecule was placed at a distance 2 Å above each positively charged magnesium cation of the silicate cluster. Subsequently, for each water + cluster combination, the whole system structure was locally optimized. All optimizations were carried out using DFT-based calculations using the PBE functional and a def2-TZVPP (i.e., valence triple-ζ plus polarization) basis set (). This process was repeated for sampling all the Mg sites and all distinct combinations of adsorption sites when using two and three water molecules. The ORCA code () was used for these calculations.
Results and Discussion
In the following subsections, we describe our findings for both hydroxylation and molecular interaction with water for bare (MgSiO3)6 nanosilicate clusters (see Figure 1). Note that these bare nanosilicates contain 30 atoms and have a typical maximum diameter of ∼1 nm. The hydrated (MgSiO3)6(H2O)n nanosilicates we consider have up to 39 atoms, which induces only a minor increase in the average diameter. Although the chemical compositions of the hydrated nanosilicate systems considered in each case (i.e., hydroxylation and molecular interactions) are identical, astronomically each relate to distinct processes and regions of the ISM. As such, we discuss the two considered modes of nanosilicate hydration separately with relation to their astronomical context.
Dissociative Hydration
As discussed above, hydroxylation of silicate dust can occur through high energy processing by supernovae shocks in the diffuse ISM. Hydroxylation is also the natural chemical result of dissociative chemisorption of water molecules on dust grain surfaces. However, as water molecules are quickly photodissociated by high-energy ultraviolet (UV) photons in the diffuse ISM, chemical routes to hydroxylated nanosilicate dust would more likely occur via reactions with other species (e.g., O and H atoms). We note that the explicit modeling of these reactions and how the resulting hydroxylation affects subsequent H2 formation reactions were investigated in previous works on a (MgSiO3)4 nanopyroxene cluster (; ; ). Generally, hydroxylation of silicate grain surfaces is the most plausible mechanism by which hydration of silicates can occur in the diffuse ISM.
We consider hydroxylated nanopyroxene clusters with chemical compositions (MgSiO3)6(H2O)n, with n = 1–3. After finding the low-energy minima of these systems using the MCBH global optimization method, although some of the higher energy isomers displayed physisorbed water molecules, most clusters exhibited 2n hydroxyls. First, this shows that hydroxylation is energetically more favorable than molecular water adsorption. The lowest energy (i.e., candidate global minima) hydroxylated isomers found for each n are shown in Figure 2. By only considering the lowest energy isomers, we can follow the evolution of the dipole moment from the bare (MgSiO3)6 with increasing hydroxylation (see Figure 2). Second, we note that the global minimum bare (MgSiO3)6 nanosilicate isomer has a negligible dipole moment, due to its high symmetry. However, all lowest energy hydroxylated nanosilicates have a significant dipole moment and correspondingly low symmetry structures. Based on this relatively small change in the degree of hydroxylation for the lowest energy isomers, the initial relatively large increase in the dipole moment going from a bare cluster to one with two hydroxyls (+5.1 Debye) is then dampened with further hydroxylation increase.
FIGURE 2
In the ISM, due to energetic processing, the population of nanosilicates is likely to exhibit a range of isomeric structures and compositions with a corresponding range of energetic stabilities. From an astronomical perspective, it is thus more relevant to extract average properties of grain populations. From this more general perspective, it has previously been shown that the high symmetry of the global minimum bare (MgSiO3)6 nanosilicate is not typical of the majority of isomers. In fact, the average dipole moment of (MgSiO3)6 isomers, considering over 100 lowest energy structures, is found to be 7.7 Debyes (). In the case of the hydroxylated nanosilicates considered herein, our global optimization searches also produce a large number of low-energy isomer structures from which to obtain the average dipole moments. In Figure 3, we show the relative energetic stability and corresponding dipole moments for approximately 100 distinct nanocluster isomers for each degree of hydroxylation (i.e., n = 1–3) considered. In Table 1, we report the average dipole moments of all hydroxylated (MgSiO3)6(H2O)n nanoclusters based on these data. It has been shown that for a dipolar nanocluster species to contribute to AME, it should possess a sufficiently large size-normalized dipole moment (β=μ / ) such that β > 0.3 (, ). In Table 1, we also reported values of β based on the average dipole moments on our sets of hydroxylated nanosilicate clusters. We note that although there is a very slight decrease in the average dipole moment with increasing hydroxylation, the value of β is always over three times the size required according to the above AME constraint. These data clearly show that in all cases, a population of hydroxylated (MgSiO3)6(H2O)n nanoclusters could significantly contribute to the observed AME.
FIGURE 3
TABLE 1
| Water addition process | Dipole moment of initial bare nanosilicates/β | Average dipole moment/β (+ 1 × H2O) | Average dipole moment/β (+ 2 × H2O) | Average dipole moment/β (+ 3 × H2O) |
|---|---|---|---|---|
| Hydroxylation (chemisorption) | 7.7a/1.41 | 6.71/1.17 | 6.60/1.1 | 6.22/1.04 |
| Molecular adsorption (physisorption) | 0.01b/0.002 | 2.87/0.50 | 3.32/0.55 | 3.30/0.55 |
| Molecular adsorption (physisorption) | 13.67c/2.50 | 14.40/2.51 | 11.97/2.00 | 12.65/2.12 |
Evolution of total dipole moments (μ) in Debyes of (MgSiO3)6 nanosilicates with increasing hydration via hydroxylation (chemisorption) and molecular water adsorption (physisorption). The average dipoles of the hydrated nanosilicates are based on the data obtained herein and shown in Figures 2–4. Values of β are calculated as β=μ / .
Molecular Hydration
Hydroxylation via shock-induced processing or suitable chemical reactions is initially very energetically favorable for bare grains, but, with increasing hydroxylation, the surface sites for the hydroxylation to take place quickly become occupied. As such, nanosized grains can only maintain a limited number of hydroxyls on their surfaces before they become too energetically favorable for any formed water to interact with the hydroxylated grain in a molecular physisorbed manner (). As mentioned above, the main mechanism through which water ice forms and grows on the core grains (e.g., silicates) is by cumulative water formation on the same grain surface, where it remains firmly attached. Here, to check how the presence of interacting water molecules affects the dipole moment of the bare nanograin, we assume that water formation has already taken place on the binding sites where the water molecules are molecularly adsorbed.
Following our results for hydroxylation, we also consider the initial stages of molecular hydration in (MgSiO3)6(H2O)n species for n = 1–3. Note that unlike in the energetic process of hydroxylation in the ISM, in protected molecular cloud environments, each molecular hydration is a low-energy event in which we assume that the bare nanocluster structure is unaffected. We take the two nanocluster structures in Figure 1 as seeds for studying the initial stages of molecular physisorption. First, we consider the global minimum low dipole nanocluster (see Figure 1A) which has a negligible dipole due to its very symmetric structure. In Figure 4, we show the lowest energy (MgSiO3)6(H2O)n structures that we obtained for the addition of 1–3 water molecules. Upon nucleation of the first water molecule, the symmetry of the cluster is broken, resulting in a dramatic increase in the dipole moment. Conversely, with two water molecules, the lowest energy structure is highly symmetric, leading to a very low dipole moment. Finally, with three water molecules, the symmetry is once again broken and the lowest energy system is found to exhibit a high dipole moment. As in the case of the hydroxylated clusters, in astronomical environments there is likely to be a range of isomers for each degree of hydration, and the average effect of molecular hydration is of more interest. In this case, due to the limited number of distinct binding sites, we have fewer possibilities for each water addition than in the case of hydroxylated nanocluster isomers. Specifically, for the nucleation of one water molecule, we find six distinct adsorption sites, and for two and three water molecules, we consider 15 distinct adsorption sites. The relative energies and dipole moments of all these systems are shown in Figure 5. The corresponding average dipole moments (see Table 1) are all close to 3 Debyes, showing that molecular water adsorption can induce an average increase in the dipole moment, even for bare nanoclusters with negligible dipole moments. The corresponding values of β are all close to 0.5, which is more than sufficient for such hydrated nanoclusters to contribute to AME.
FIGURE 4
FIGURE 5
In the case of molecular nucleation of a bare cluster with an initial large dipole moment (see Figure 1B), we show the lowest energy (MgSiO3)6(H2O)n structures that we obtained in our searches in Figure 6. For these specific low-energy structures, we see that the dipole moment stays high (>9 Debyes) for each water addition. More generally, we consider a similar range of hydrated structures as in the low dipole moment case to obtain averages. In Figure 7, we show the relative energies and dipole moments of all considered systems. Here, we also see that the average dipole moments (see also Table 1) for each degree of hydration also stay high. Correspondingly, the values of β obtained from these average dipole moment values are also high (>2 Debyes), meaning that they would be strong emitters of microwave radiation.
FIGURE 6
FIGURE 7
Conclusion
We use searches based on global optimization and electrostatic considerations to obtain a range of structures for hydroxylated and molecularly hydrated nanoclusters, respectively, with a (MgSiO3)6(H2O)n composition (n = 1–3). These hydrated nanopyroxenes could be relevant for nanosilicate grain processing in the ISM or initial water ice mantle formation in molecular clouds. We have used accurate quantum chemical DFT calculations to calculate the dipole moments of all hydrated pyroxene nanoclusters. We find that both hydroxylation and molecular hydration, irrespective of the initial dipole moment of the bare nanocluster, induce and/or maintain systems with dipole moments that are sufficiently high to significantly contribute to AME. These results thus strongly indicate that the validity of previous reports showing that nanosilicates are very plausible carriers of AME is not affected by a modest degree of nanosilicate hydration. Overall, this study broadens the range of potential nanosilicate species and types of astrophysical environments that could produce AME signatures.
Statements
Data availability statement
The data supporting the findings of this study are available on reasonable request from the corresponding authors.
Author contributions
The original idea for the study was proposed by SB. All authors participated in the initial conception and planning of the study and in the discussion of the results. JMG and SF performed the calculations. SB and AR participated in the writing and editing of the manuscript and coordinated the project.
Funding
This project has received funding from the European Union’s Horizon 2020 research and innovation program under the Marie Sklodowska-Curie grant agreement No. 811312 for the project “Astro-Chemical Origins” (ACO) and from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program grant agreement No. 865657 for the project “Quantum Chemistry on Interstellar Grains” (QUANTUMGRAIN). We acknowledge financial support from the Spanish Ministerio de Ciencia, Innovación y Universidades (projects CTQ 2017-89132-P, RTI 2018-095460-B-100, and MDM-2017-0767 via the Spanish Structures of Excellence María de Maeztu program) and the Generalitat de Catalunya (projects 2017SGR1323 and 2017SGR13).
Acknowledgments
AR is indebted to the “Ramón y Cajal” program. The Red Española de Supercomputación (RES) and Consorci de Serveis Universitaris de Catalunya (CSUC) are also acknowledged for the provision of supercomputing time.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
Ali-HaïmoudY. (2013). Spinning Dust Radiation: A Review of the Theory. Adv. Astron.2013, 462697. 10.1155/2013/462697
2
BlumV.GehrkeR.HankeF.HavuP.HavuV.RenX.et al (2009). Ab Initio Molecular Simulations with Numeric Atom-Centered Orbitals. Comput. Phys. Commun.180, 2175–2196. 10.1016/j.cpc.2009.06.022
3
BoogertA. C. A.PontoppidanK. M.KnezC.LahuisF.Kessler‐SilacciJ.van DishoeckE. F.et al (2008). The c2dSpitzerSpectroscopic Survey of Ices Around Low‐Mass Young Stellar Objects. I. H2O and the 5-8 μm Bands1,2. ApJ678, 985–1004. 10.1086/533425
4
BradleyJ. P. (1994). Chemically Anomalous, Preaccretionally Irradiated Grains in Interplanetary Dust from Comets. Science265, 925–929. 10.1126/science.265.5174.925
5
CukoA.MaciàA.CalatayudM.BromleyS. T. (2017). Global Optimisation of Hydroxylated Silica Clusters: A Cascade Monte Carlo Basin Hopping Approach. Comput. Theor. Chem.1102, 38–43. 10.1016/j.comptc.2016.12.030
6
DemykK.JonesA. P.DartoisE.CoxP.d´HendecourtL. (1999). The Chemical Composition of the Silicate Dust Around RAFGL7009S and IRAS 19110+1045. Astron. Astrophys.349, 267–275.
7
DickinsonC.Ali-HaïmoudY.BarrA.BattistelliE. S.BellA.BernsteinL.et al (2018). The State-Of-Play of Anomalous Microwave Emission (AME) Research. New Astron. Rev.80, 1–28. 10.1016/j.newar.2018.02.001
8
DraineB. T.LazarianA. (1998a). Electric Dipole Radiation from Spinning Dust Grains. ApJ508, 157–179. 10.1086/306387
9
DraineB. T.LazarianA. (1998b). Diffuse Galactic Emission from Spinning Dust Grains. Astrophys. J.494, L19–L22. 10.1086/311167
10
DulieuF.AmiaudL.CongiuE.FillionJ.-H.MatarE.MomeniA.et al (2010). Experimental Evidence for Water Formation on Interstellar Dust Grains by Hydrogen and Oxygen Atoms. Astron. Astrophys.512, A30. 10.1051/0004-6361/200912079
11
GaleJ. D.RohlA. L. (2003). The General Utility Lattice Program (GULP). Mol. Simul.29, 291–341. 10.1080/0892702031000104887
12
GibbE. L.WhittetD. C. B.BoogertA. C. A.TielensA. G. G. M. (2004). Interstellar Ice: The Infrared Space Observatory Legacy. Astrophys J. Suppl. S151, 35–73. 10.1086/381182
13
GobrechtD.CherchneffI.SarangiA.BromleyS. T.BromleyS. T. (2016). Dust Formation in the Oxygen-Rich AGB Star Ik Tauri. A&A585, A6. 10.1051/0004-6361/201425363
14
GoumansT. P. M.BromleyS. T. (2011). Hydrogen and Oxygen Adsorption on a Nanosilicate - a Quantum Chemical Study. Mon. Not. R. Astron. Soc.414, 1285–1291. 10.1111/j.1365-2966.2011.18463.x
15
GoumansT. P. M.BromleyS. T. (2012). Efficient Nucleation of Stardust Silicates via Heteromolecular Homogeneous Condensation. Mon. Not. R. Astron. Soc.420, 3344–3349. 10.1111/j.1365-2966.2011.20255.x
16
HenningT. (2010). Cosmic Silicates. Annu. Rev. Astron. Astrophys.48, 21–46. 10.1146/annurev-astro-081309-130815
17
HensleyB. S.DraineB. T. (2017). Modeling the Anomalous Microwave Emission with Spinning Nanoparticles: No PAHs Required. ApJ836, 179. 10.3847/1538-4357/aa5c37
18
HensleyB. S.DraineB. T.MeisnerA. M. (2016). A Case against Spinning Pahs as the Source of the Anomalous Microwave Emission. ApJ827, 45. 10.3847/0004-637x/827/1/45
19
HildebrandtS. R.ReboloR.Rubiño-MartínJ. A.WatsonR. A.GutiérrezC. M.HoylandR. J.et al (2007). COSMOSOMAS Observations of the Cosmic Microwave Background and Galactic Foregrounds at 11 GHz: Evidence for Anomalous Microwave Emission at High Galactic Latitude. Mon. Not. R. Astron. Soc.382, 594–608. 10.1111/j.1365-2966.2007.12380.x
20
HoangT.VinhN.-A.LanN. Q. (2016). Spinning Dust Emission from Ultra-small Silicates: Emissivity and Polarization Spectrum. ApJ824, 18. 10.3847/0004-637x/824/1/18
21
HoangT.TramL. N.LeeH.AhnS.-H. (2019). Rotational Disruption of Dust Grains by Radiative Torques in Strong Radiation Fields. Nat. Astron.3, 766–775. 10.1038/s41550-019-0763-6
22
HollenbachD.McKeeC. F. (1989). Molecule Formation and Infrared Emission in Fast Interstellar Shocks. III - Results for J Shocks in Molecular Clouds. ApJ342, 306. 10.1086/167595
23
IoppoloS.CuppenH. M.RomanzinC.van DishoeckE. F.LinnartzH. (2008). Laboratory Evidence for Efficient Water Formation in Interstellar Ices. Astrophys. J.686, 1474–1479. 10.1086/591506
24
JonesA. P.TielensA. G. G. M.HollenbachD. J. (1996). Grain Shattering in Shocks: The Interstellar Grain Size Distribution. ApJ469, 740. 10.1086/177823
25
KellerL. P.McKayD. S. (1997). The Nature and Origin of Rims on Lunar Soil Grains. Geochim. Cosmochim. Acta61, 2331–2341. 10.1016/s0016-7037(97)00085-9
26
KemperF.VriendW. J.TielensA. G. G. M. (2004). The Absence of Crystalline Silicates in the Diffuse Interstellar Medium. ApJ609, 826–837. 10.1086/421339
27
KerkeniB.BromleyS. T. (2013). Competing Mechanisms of Catalytic H2 Formation and Dissociation on Ultrasmall Silicate Nanocluster Dust Grains. Mon. Not. R. Astron. Soc.435, 1486–1492. 10.1093/mnras/stt1389
28
KerkeniB.Bacchus-MontabonelM.-C.BromleyS. T. (2017). How Hydroxylation Affects Hydrogen Adsorption and Formation on Nanosilicates. Mol. Astrophys.7, 1–8. 10.1016/j.molap.2017.04.001
29
KerkeniB.Bacchus-MontabonelM.-C.ShanX.BromleyS. T. (2019). Understanding H2 Formation on Hydroxylated Pyroxene Nanoclusters: Ab Initio Study of the Reaction Energetics and Kinetics. J. Phys. Chem. A.123, 9282–9291. 10.1021/acs.jpca.9b06713
30
Lamiel-GarciaO.KoK. C.LeeJ. Y.BromleyS. T.IllasF. (2017). When Anatase Nanoparticles Become Bulklike: Properties of Realistic TiO2 Nanoparticles in the 1-6 Nm Size Range from All Electron Relativistic Density Functional Theory Based Calculations. J. Chem. Theor. Comput13, 1785–1793. 10.1021/acs.jctc.7b00085
31
LiA.DraineB. T. (2001). On Ultrasmall Silicate Grains in the Diffuse Interstellar Medium. Astrophys. J. Lett.550, L213–L217. 10.1086/319640
32
Macià EscatllarA.BromleyS. T. (2020). Assessing the Viability of Silicate Nanoclusters as Carriers of the Anomalous Microwave Emission: a Quantum Mechanical Study. Astronom. Astrophys.634, A77. 10.1051/0004-6361/201936419
33
Macià EscatllarA.LazaukasT.WoodleyS. M.BromleyS. T. (2019). Structure and Properties of Nanosilicates with Olivine (Mg2SiO4) N and Pyroxene (MgSiO3) N Compositions. ACS Earth Sp. Chem.3, 2390–2403. 10.1021/acsearthspacechem.9b00139
34
MokraneH.ChaabouniH.AccollaM.CongiuE.DulieuF.ChehrouriM.et al (2009). Experimental Evidence for Water Formation via Ozone Hydrogenation on Dust Grains at 10 K. ApJ705, L195–L198. 10.1088/0004-637x/705/2/l195
35
MolpeceresG.RimolaA.CeccarelliC.KästnerJ.UgliengoP.MatéB. (2019). Silicate-mediated Interstellar Water Formation: A Theoretical Study. Mon. Not. R. Astron. Soc.482, 5389–5400. 10.1093/mnras/sty3024
36
NeeseF.WennmohsF.BeckerU.RiplingerC. (2020). The ORCA Quantum Chemistry Program Package. J. Chem. Phys.152, 224108. 10.1063/5.0004608
37
NorrisB. R. M.TuthillP. G.IrelandM. J.LacourS.ZijlstraA. A.LykouF.et al (2012). A Close Halo of Large Transparent Grains Around Extreme Red Giant Stars. Nature484, 220–222. 10.1038/nature10935
38
PerdewJ. P.BurkeK.ErnzerhofM. (1996). Generalized Gradient Approximation Made Simple. Phys. Rev. Lett.77, 3865–3868. 10.1103/physrevlett.77.3865
39
TielensA. G. G. M.HagenW. (1982). Model Calculations of the Molecular Composition of Interstellar Grain Mantles. Astron. Astrophys.114, 245
40
van DishoeckE. F.HerbstE.NeufeldD. A. (2013). Interstellar Water Chemistry: From Laboratory to Observations. Chem. Rev.113, 9043–9085. 10.1021/cr4003177
41
WalesD. J.DoyeJ. P. K. (1997). Global Optimization by Basin-Hopping and the Lowest Energy Structures of Lennard-Jones Clusters Containing up to 110 Atoms. J. Phys. Chem. A.101, 5111–5116. 10.1021/jp970984n
42
WatanabeN.KouchiA. (2008). Ice Surface Reactions: A Key to Chemical Evolution in Space. Prog. Surf. Sci.83, 439–489. 10.1016/j.progsurf.2008.10.001
43
WatsonR. A.ReboloR.Rubiño-MartínJ. A.HildebrandtS.GutiérrezC. M.Fernández-CerezoS.et al (2005). Detection of Anomalous Microwave Emission in the Perseus Molecular Cloud with the COSMOSOMAS Experiment. ApJ624, L89–L92. 10.1086/430519
44
WeigendF.AhlrichsR. (2005). Balanced Basis Sets of Split Valence, Triple Zeta Valence and Quadruple Zeta Valence Quality for H to Rn: Design and Assessment of Accuracy. Phys. Chem. Chem. Phys.7, 3297. 10.1039/b508541a
45
ZellerE. J.RoncaL. B.LevyP. W. (1966). Proton-induced Hydroxyl Formation on the Lunar Surface. J. Geophys. Res.71, 4855–4860. 10.1029/jz071i020p04855
Summary
Keywords
anomalous microwave emission, dust grains, nanosilicates, water ices, interstellar medium, grain processing, hydroxylation, density functional theory
Citation
Mariñoso Guiu J, Ferrero S, Macià Escatllar A, Rimola A and Bromley ST (2021) Does Processing or Formation of Water Ice Mantles Affect the Capacity of Nanosilicates to Be the Source of Anomalous Microwave Emission?. Front. Astron. Space Sci. 8:676548. doi: 10.3389/fspas.2021.676548
Received
05 March 2021
Accepted
30 April 2021
Published
20 May 2021
Volume
8 - 2021
Edited by
Luca Bizzocchi, Max Planck Institute for Extraterrestrial Physics, Germany
Reviewed by
Martin Robert Stewart McCoustra, Heriot-Watt University, United Kingdom
Sergio Ioppolo, Queen Mary University of London, United Kingdom
Updates
Copyright
© 2021 Mariñoso Guiu, Ferrero, Macià Escatllar, Rimola and Bromley.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Albert Rimola, albert.rimola@uab.cat; Stefan T. Bromley, s.bromley@ub.edu
This article was submitted to Astrochemistry, a section of the journal Frontiers in Astronomy and Space Sciences
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.