Abstract
Titan, the largest moon of Saturn, is one of the key planetary objects in the field of exobiology. Its dense, nitrogen-rich atmosphere is the site of important organic chemistry. This paper focuses on the organic aerosols produced in Titan's atmosphere that play an important role in atmospheric and surface processes and in organic chemistry as it applies to exobiological interests. To produce reliable laboratory analogues of these aerosols, we developed, tested, and optimized a device for the synthesis of clean tholins. The potential chemical evolution of Titan aerosols at Titan's surface has been studied, in particular, the possible interaction between aerosols and putative ammonia-water cryomagma. Modeling of the formation of Saturn's atmosphere has permitted the characterization of a composition of salts in the subsurface ocean and cryolava. From this new and original chemical composition, a laboratory study of several hydrolyses of tholins was carried out. The results obtained show the formation of many organic compounds, among them, species identified only in the presence of salts. In addition, a list of potential precursors of these compounds was established, which could provide a database for research of the chemical composition of tholins and/or aerosols of Titan. Key Words: Titan tholins—Titan aerosols—Hydrolysis—Carbonates—Titan's surface. Astrobiology 17, 8–26.
1. Introduction
T
Titan's atmospheric composition is driven by the coupled chemistry of N2 and CH4 (Raulin et al., 2012, and references therein). The complex organic chemistry occurring in this atmosphere is initiated in the ionosphere and involves large negative ions (Waite et al., 2007). The dissociation of methane induced by photons with wavelengths lower than 150 nm provides CH, CH2, and CH3 radicals. The dissociation of dinitrogen is generated by photons with wavelengths lower than 100 nm and by energetic particles from Saturn's magnetosphere. Then the resulting products chemically evolve, coagulate, and/or condense. This evolution yields more complex compounds up to the production of macromolecular organic particles, forming Titan's haze. Those complex organic aerosols provide the characteristic orange-brown color of Titan in the visible spectral range. They also play an important role in radiative transfer models (McKay et al., 1989, 1991), as well as in the composition and properties of the atmosphere of Titan. And they may interact chemically with geological features present on the surface and affect the chemical composition of the satellite's exterior. With the Cassini-Huygens mission, in particular the Synthetic Aperture Radar (SAR) and the Visual and Infrared Mapping Spectrometer (VIMS) images, the surface of Titan has been revealed to be geologically complex and diversified (Elachi et al., 2005; McCord et al., 2006). Although several large plains remain unknown (Lopes et al., 2010), many geological features present on Earth have been observed on Titan's surface, such as volcanoes, craters, sand dunes, lakes, seas, and mountains. These geological structures are constantly covered by solid organic aerosols deposited by precipitation or dry deposition. For those reasons, it is important to understand how Titan's aerosols chemically evolve once deposited at Titan's surface in order to determine whether astrobiologically interesting molecules could be formed.
In this paper, we focus on the likely cryovolcanic area potentially identified on Titan's surface (Barnes et al., 2005, 2006; Lopes et al., 2007, 2010, 2013; Nelson et al., 2009; Soderblom et al., 2009; Wall et al., 2009). Cryovolcanism refers to the eruption of materials of low density and low melting point, which constitute part of Titan's interior (as the plausible subsurface ocean). The presence of cryovolcanic processes has been suggested by Lorenz (1993) and Lorenz and Lunine (1996). It has been tentatively observed after Cassini's arrival to the Saturn system. In particular, the SAR and VIMS data have provided images from the surface that highlight potential cryovolcanic activity. With VIMS data, it is possible to detect surface changes that could indicate a currently active cryovolcanism on Titan. But the interpretation of those data is still controversial since the SAR images alone lead to a different conclusion, compared to those derived from both VIMS and SAR data. Lopes et al. (2007) observed signs of cryovolcanism flow in the Ganesa region on Titan via SAR data processing. Also, it has been pointed out that cryovolcanism could be a minor geological process with sporadic activity (Tobie et al., 2006; Mitri et al., 2008). Indeed, SAR data show that the cryovolcanism area would only represent a tiny part of the surface: roughly 1.2% between 90°N and 90°S (an area corresponding to 20% of the total titanian surface). Nevertheless, plains cover a large part of the surface of Titan, but their nature is still unclear. If their origin is cryovolcanic, it could suggest that cryovolcanism would have been involved at one moment in the history of Titan (Lopes et al., 2010, 2013). Lopes et al. (2010, 2013) proposed that most of the cryovolcanic features are located from 30°S to 60°N (Table 1). Within this location, two regions seem to be good candidates for cryovolcanism, in particular because of the presence of caldera-like structures and potential surface changes, Hotei Regio and Tui Regio (Barnes et al., 2005, 2006; Nelson et al., 2009; Soderblom et al., 2009; Wall et al., 2009; Lopes et al., 2010, 2013), which do not exclude the assumption of the sedimentary alternative.
Lopes et al. (2007, 2010, 2013) also suggested that Ganesa Macula, which was first interpreted as a volcanic dome (Elachi et al., 2005; Lopes et al., 2007), could be a cryovolcanism region. Indeed, recent studies have shown that the features of Ganesa do not fit with a volcanic dome (Kirk et al., 2008; Stiles et al., 2009). Furthermore, by combining RADAR and VIMS data, Lopes et al. (2011) explained that Ganesa Macula could not be a cryovolcanic area. A similar interpretation has been made for the Winia Fluctus region. The Sotra Facula region has also been suggested as a cryovolcanic feature candidate via SAR (Lopes et al., 2010) and RADAR stereogrammetry (Kirk et al., 2008). This region could also be partially covered by dunes, which would imply that the cryovolcanic features are relatively young (Lopes et al., 2010). This region includes the sites named Mohini Fluctus, Sotra Patera, Doom Mons, and Erebor Mons. It has been suggested that these four areas have been shaped by cryovolcanic processes (Lopes et al., 2010). Rohe Fluctus and Ara Fluctus are also considered potential candidates for cryovolcanism since they show caldera-like depressions associated with flow features as is the case for Mohini Fluctus (Lopes et al., 2007, 2013). Wall et al. (2009) showed that both Western Xanadu and Hotei Arcus regions exhibit variable infrared brightness that could be evidence of active cryovolcanoes. Nevertheless, Soderblom et al. (2009) questioned the approach used to extract the reflectance data, and they take this hypothesis with caution. Sotin et al. (2005) interpreted that the region Tortola Facula (also known as Sotra Patera) from VIMS data could be a cryovolcanic area, while Sotin et al. (2010) and Hayes et al. (2008) interpreted this region to be “hummocky and mountainous terrain units.” Lopes et al. (2011) also concluded that this region would have other origins than cryovolcanism.
In summary, Ganesa Macula and Winia Fluctus could not be cryovolcanic areas. For Tui Regio and Western Xanadu, the cryovolcanic interpretation remains uncertain, while Sotra Facula and Hotei Regio could be the strongest candidates for cryovolcanic regions (Lopes et al., 2011, 2013). Nevertheless, the low resolution of the data sets is a real issue, so several cryovolcanic features may not have been identified yet.
Although its composition remains unknown, it is often suggested that the cryomagma could be a mixture of water-ammonia (Grasset et al., 2000; Tobie et al., 2005) and may also include methanol (Lopes et al., 2007). Mitri et al. (2008) made the assumption that the composition of the cryomagma could be connected to the composition of the internal ocean. Models predict the presence of a water-ammonia ocean located between a low-pressure water ice layer and a high-pressure water ice layer (Tobie et al., 2005, 2012). The possibility of such an internal ocean is supported by several observations such as the changing of Titan's rotation (Lorenz et al., 2008) or the uncommon Schumann resonance in Titan's atmosphere (Béghin et al., 2012). Mitri et al. (2008) showed that the putative water-ammonia cryomagma could slowly reach the surface through cracks at the base of the low-pressure water ice layer. The fraction of ammonia usually considered in experimental simulation works was 10–25 wt % (Neish et al., 2009; Poch et al., 2012). However, a recent study shows that the fraction of ammonia in the subsurface ocean could be lower, around 2–3 wt % if taking into account the actual mass of nitrogen in Titan's atmosphere (Tobie et al., 2012). Additionally, studies have provided new highlights on the bulk composition of Titan for various gas species. Indeed, Saturn's observed atmosphere enrichment constrains the composition of the planetesimals present in the feeding zone of Saturn. The enrichment in volatiles in Saturn's atmosphere has been reproduced by assuming (Hersant et al., 2004, 2008) that (i) carbon is mostly in the form of CH4 clathrate hydrates and condensed CO2; (ii) nitrogen is mostly in the form of NH3 hydrates; (iii) sulfur is mostly in the form of H2S clathrate hydrates. In the present study, we assumed that those gas species have been trapped in the likely internal ocean. By taking into account the plausible acid-base properties of a water-ammonia ocean, a new probable composition of the cryomagma that could potentially interact with Titan's deposited aerosols has been determined (see Section 2.2).
In this work, hydrolyses were performed on Titan aerosol analogues with 5 wt % ammonia aqueous solution with and without salts (Na2CO3 and NaHS) in order to simulate the new probable composition of the cryomagma described above. Previous experiments have tested different kinds of hydrolyses: (i) acidic pH hydrolysis with HCl 6 mol L−1 (Khare et al., 1986; Nguyen, 2007; Taniuchi et al., 2013), HCl 2 mol L−1 (Nguyen, 2007); (ii) hydrolysis at pH 7 (Nguyen, 2007); (iii) basic pH hydrolysis with 25 wt % NH3 (Ramirez et al., 2010; Poch et al., 2012), with 13 wt % NH3 (Neish et al., 2009; Ramirez et al., 2010) or less (Ramirez et al., 2010). The experimental conditions of these works are summarized below in Section 3.2. Amino acids, alcohols, carboxylic acids, urea, and nucleobases have been identified as hydrolysis products. Only two studies have reported quantitative data on the production of these organics of astrobiological interest (Ramirez et al., 2010; Poch et al., 2012).
The experimental device used to synthesize representative free-oxygen analogues of Titan's aerosols (named Titan tholins) is detailed, as well as the calculation of the molar concentration of each minor species that could constitute the likely subsurface ocean. The detection/identification/quantification methods of hydrolysis products are described and the results related to the low-temperature alkaline pH hydrolysis of oxygen-free Titan tholins are presented and discussed.
2. Experimental
2.1. Plasma experimental setup
To synthesize representative analogues of Titan aerosols, the experimental setup (developed at LISA, Créteil, France) named PLASMA has been used (Coll et al., 1997). This device allows for simulation of the dissociation of the main gas molecules present in Titan's atmosphere by the establishment of a cold plasma in a gas mixture of N2 and CH4.
As shown in Fig. 1, the reactor is connected to a gas cylinder composed of a mixture of N2/CH4 (98/2 ratio), which is representative of the mean gas composition of Titan's atmosphere. This gas mixture continuously flows through the reactor, which allows for maintenance of the N2/CH4 ratio constant during the synthesis. A flowmeter/restrictor (Bronkhorst F-201CV-100-AAD-11-V) placed at the cylinder exit and upstream of the reactor maintains a constant gas flow of 20 sccm (1 sccm = 7.4 × 10−7 mol s−1: standard cubic centimeters per minute) and a pressure of 2.6 mbar. A dry pump (Edwards nXDS6i, maximum pumping speed = 6.2 m3 h−1) is used to maintain the low pressure in the system.

Scheme of the experimental setup, named PLASMA, used for the synthesis of Titan tholins.
By using a high-voltage generator, an electric discharge (2 kV–23.2 mA) is applied on the gas mixture, which establishes a cold plasma that induces the dissociation of N2 and CH4 and other chemical reactions leading to the formation of tholins that settle on the reactor's wall. The U-shape of the reactor (Fig. 1) facilitates the introduction of the tholins into a Dewar vessel full of liquid nitrogen and allows the simulation to continue at a very low temperature. However, in the present study, the synthesis was conducted only at ambient temperature since the production yields of tholins, at low temperature, are too low to provide the quantity of tholins necessary to achieve this new and exploratory experimental study.
Previous studies have shown that tholins synthesized with the PLASMA experimental setup were contaminated by the atmospheric molecular oxygen of the laboratory (Poch et al., 2012). Consequently, the entire experimental device, with the exception of the N2/CH4 gas cylinder and the high-voltage generator, were put into a glove box to insulate the experiment from the terrestrial atmosphere. The glove box was connected to a N2 gas cylinder to maintain a high pressure of internal gas in it. Moreover, silica gel granules were placed in the glove box to capture the residual water vapor. After 3 days and nights of synthesis, the condensed tholins are collected by scratching the reactor's walls. This step occurred in the glove box as well.
2.2. Plausible subsurface ocean chemical composition: theoretical calculation
The chemical composition of the cryomagma is likely in close association with the chemical composition of the potential subsurface ocean. So to simulate the cryovolcanism area (assumed to be made of cryomagma), we determined the possible chemical composition of the internal ocean. Models of the formation of Titan suggest that the water ocean was initially in contact with the atmosphere (Fortes, 2000). We assumed that the composition of the present ocean is similar to that of the early ocean, and neglected the pressure effect on the pH and the solubility.
If we consider that this ocean is only composed of 5 wt % NH3, that is, 2.88 mol L−1, then its pH can be determined with the classical weak base formula:
With pKa = 9.2, it comes to pH = 11.83, corresponding to (H+) = 1.5 × 10−12 mol L−1.
Knowledge of the molar concentration of ammonia and the plausible acid-base properties of the water-ammonia ocean allow for calculation of the concentration of the other minor species, which include NH4 +, OH−, H2S, HS−, S2−, CO2(aq), HCO3 −, CO3 2−.
2.2.1. NH4 +/NH3
We assume a concentration of 5 wt % of ammonia. This is equivalent to the molar concentration:
Then with the acid dissociation constant Ka
With a pKa of 9.2 and a pH of 11.83, it comes to
2.2.2. H+/OH−
The concentration of H+ is directly deduced from the pH value
And following in the same way, from the water equilibrium, the concentration of OH− is obtained by using the expression
Then the value of the hydroxide concentration is
2.2.3. H2S/HS−/S2−
For carbon dioxide and hydrogen sulfide, another parameter was needed to determine their maximum concentration: their solubility in water for the likely ocean's temperature.
At the temperature of T = 262 K, the solubility of H2S is 8.7 g(H2S)/kg(H2O). So the maximum concentration for sulfurs in water is
The mass balance for sulfides gives the following expression:
We can also write
With
it comes to
Then
With Ka1 = 10−7 and Ka2 = 10−12.9 at pH = 11.83, we get
and
2.2.4. CO2(aq)/HCO3 −/CO3 2−
As for sulfides, we need the solubility of carbonates in water to determine the concentration of all their species.
For a temperature of T = 262 K, the solubility of CO2 is 4.2 g(CO2)/kg(H2O). So the maximum concentration of carbonates in water is
With
this expression becomes
With
then it becomes
Then
With K′a1 = 10−6.4 and K′a2 = 10−10.3 at pH = 11.83, we get
and
2.2.5. Na+
Since the presence of sodium salt (as sodium carbonate for example) in the ocean is now suspected (Mitri et al., 2014), we have chosen to use the Na+ ion as a counter ion to get the electroneutrality in the ocean. By applying the electroneutrality law to the mixture, we get the following expression:
Then, it comes to
Table 2 summarizes the possible chemical composition of the likely subsurface ocean. Hydrogenosulfide and carbonate anions appear to be among the most abundant minor species.
2.2.6. The influence of the ocean's temperature on the pH value
The previous calculated pH value is valid at an ambient temperature of 298 K. However, Titan's ocean temperature would be close to 262 K (lower limit of the ocean)/270 K (upper limit of the ocean) (Tobie et al., 2012).
To determine the influence of the ocean's temperature and the ionic strength on the pH value, the Davies approach (Butler, 1998) has been followed to estimate the ion activity coefficients for each temperature.
With
γi = ion activity coefficient
I = ionic strength
z i = charge of the ion i
A = parameter depending on the temperature (T) and the dielectric constant (ɛ)
= 1.82 × 106 × (ɛT)−3/2
Since liquid water is the solvent, we use the dielectric constant of water for the calculations, ignoring the presence of ammonia. Its value is given by the Handbook of Chemistry and Physics, 87th edition (Lide, 2006, p 6–133).
The ionic strength I is given by the following equation:
With this new calculation method, it is possible to avoid the initial assumption, which was that the ion activity coefficient γi = 1.
By using the calculated concentration from Table 2, it is possible to determine the potential ionic strength of the internal ocean.
It is now possible to obtain the average ion activity coefficient at different temperatures.
which can be obtained with the Davies approach as mentioned previously
and
Finally, it is possible to calculate accurately the pH value by taking into account the ocean temperature.
Table 3 displays the results concerning the ion activity coefficients and the corrected pH values. We can see that the pH in the ocean would more likely be 11.76. So it seems that the temperature and the presence of salts in the ocean do not impact strongly the pH value.
Indeed, the molar concentration of the minor species has been recalculated for a pH of 11.76 as shown in Table 4. Then, we show that the values obtained are not significantly different from those calculated at pH 11.83, especially for the most abundant minor species HS− and CO3 2−.
Those new data have been included in our hydrolysis experiments as detailed in the following section.
2.3. Hydrolysis conditions
To integrate this new plausible composition of the cryomagma in our hydrolysis experiments, four different kinds of hydrolysis were performed on 10 mg of oxygen-free tholins as indicated in Table 5.
All hydrolysis solutions included 5 wt % of ammonia in water. One has carbonate also (hydrolysis n°3); another one has sulfide (hydrolysis n°2). The last one (hydrolysis n°4) has both carbonate and sulfide in order to identify the contribution of each species in the production of organics. Each hydrolysis was studied at three different temperatures: 293, 279, and 253 K. For each kind of hydrolysis, reference solutions were placed in the same experimental conditions.
2.4. Gas chromatography–mass spectrometry
After 10 weeks of evolution at the three different temperatures mentioned above, the aqueous phase was analyzed by gas chromatography–mass spectrometry using a Thermo Scientific TRACE gas chromatograph (GC) coupled to a Thermo Scientific PolarisQ mass spectrometer (MS; trap ion detection mode). The PTV injector was heated at 250°C with a split ratio of 1/25. The GC column was a fused silica capillary Phenomenex Zebron Inferno ZB-5HT (L = 30 m × i.d. = 0.25 mm × dt = 0.25 μm). Its almost apolar stationary phase (5% phenyl and 95% dimethylpolysiloxane) was efficient to separate derivatized products, which are molecules of very low polarity. The GC temperature program started at 70°C for 5 min, then with a temperature gradient of 6°C/min it increased from 70°C to 350°C, and finally it stayed 5 min more at 350°C. The carrier gas used was helium (99.9999%) with a constant flow of 1 mL min−1. To protect the filament from derivatization reagents, the start time of scan began after the first 8 min, allowing the remaining derivatization agents to be fully eluted. The electronic impact source was set at 70 eV, and the mass ranged from 40 to 550 m/z.
2.5. Characterization of hydrolysis products: derivatization and extraction phase
The GCMS analytical device requires a derivatization step for decreasing the polarity and increasing the stability and volatility of most of the produced organics, allowing for the analysis of compounds of astrobiological interest. Indeed, amino acids, carboxylic acids, and nucleobases are very polar and of low volatility. As described in Fig. 2, N-methyl-N-(tert-butyldimethylsilyl)-trifluoroacetamide (MTBSTFA) with dimethylformamide (DMF) as a solvent has been used as a derivatization solution (Buch et al., 2006, 2009). A total of 30 μL of MTBSTFA and 10 μL of DMF were added to each sample (and to each reference solution), which had been dried before, under a nitrogen flow at 40°C. This mixture was heated at 75°C for 15 min to perform the derivatization reaction.

Scheme of the MTBSTFA derivatization reaction.
Before each sample injection, two blank analyses were performed. The first one was a real blank without any injection, and the second one was an injection with only the derivatization reagent. In this way, we were able to identify the GC peaks directly produced by the derivatization reagent and solvent.
To characterize the products, their respective standards were analyzed, and a database of derivatized molecules was established, as detailed in Table 6.
2.6. Influence of salt on the derivatization reaction
In the presence of salts, which includes all samples except for hydrolysis n°1, an additional step is needed before the derivatization in order to eliminate the salts, which could interfere with the derivatization. For that purpose, the sample was eluted from a cationic exchange resin. This step was performed with an AG 50W-X8 resin #142-1441 (Bio-rad Laboratories, Inc.) so that neutral molecules and anions do not interact with the resin. Once the sample was dried up, pure distilled water was added in order to eliminate the ammonia (to avoid converting the counter-ions in the resin from H+ to NH4 +). Indeed, the counter-ions in the resin (here H+) are substituted by ions present in the sample that have the same charge (here Na+). In that case, by removing Na+ (replaced by H+) in the sample solution, Na2CO3 does not form, so the derivatization can occur as usual.
The sample was introduced into the hydrated resin with a Pasteur pipette and collected at the exit: this was solute n°1. Next, to recover carboxylic acids and amino acids trapped by the resin, a solution of ammonia was added to the resin: solute n°2 was collected. Finally, for recovering poor water soluble carboxylic acids, an aqueous solution of ethanol was introduced into the resin: solute n°3 was collected. Those three solutes were then dried separately before being derivatized (as detailed above in this section). This step was efficient for hydrolysis containing only carbonates (i.e., hydrolysis n°3). However, we have been unable, to date, to perform this extraction step correctly on samples that contain sulfur (i.e., hydrolysis n°2 and n°4). To resolve this issue, we plan to use a UHPLC analytical technique.
2.7. Hydrolysis products' quantification and production yields
To quantify and determine the production yields of the compounds produced from the hydrolysis step, 4 μL of an internal standard, methyl-laurate at 7.9 × 10−2 mol L−1 (Fluka, 97%), was added to each derivatized solution before the injection in the GCMS. To achieve quantitative measurements, we carried out a calibration curve for each targeted compound. For each molecule, a calibration curve (see Table 7) was traced by using the ratio of the chromatographic peak area of the compound over the chromatographic peak area of the internal standard (methyl-laurate) as a function of the amount of the targeted compound. Then by using the calibration curve (Table 7) it was possible to determine the exact quantification of the molecules produced after hydrolysis. The production yields were obtained by calculating the ratio of the mass of the product to the mass of hydrolyzed tholins.
To obtain a faithful quantification, for each injection made for the calibration and for the sample, a minimum of three replicas was achieved. Those replicas allowed us to formulate accurate error bars for each value.
3. Results
The chromatogram acquired for nonhydrolyzed tholins after derivatization (as detailed in Section 2.5) is shown in Fig. 3. It can be seen that the chromatogram does not display any oxygenated molecules. This demonstrates that the identified products from the hydrolysis experiments are due to the hydrolysis of tholins and not to oxygen contamination (from the air) of the tholins and guarantees the reliability of the data presented in this paper.

Chromatogram of a derivatized sample of nonhydrolyzed tholins. Positions of the chromatographic peaks of compounds identified by their standards are shown by arrows. The 30 m × 0.25 mm × 0.25 μm Zebron Inferno ZB-5HT column, operated in the split mode (1:25), was programmed at 6°C min−1 from 70°C to 300°C with an inlet flow of helium of 1 mL min−1.
3.1. Hydrolysis with and without carbonates
Figure 4 shows an example of a chromatogram obtained by gas chromatography–mass spectrometry analysis of tholin hydrolysis. The chromatogram displayed is related to the hydrolysis of free-oxygen tholins in a water-ammonia mixture at 293 K after derivatization and reports the hydrolysis products identified. Moreover, Fig. 5 presents the whole set of data obtained for the two different experiments: hydrolysis n°1 and hydrolysis n°3. As presented in Table 5, hydrolysis n°1 resulted from the water-ammonia (5 wt %) mixture alone, whereas hydrolysis n°3 includes carbonates. Figure 5 provides the production yields for each hydrolysis product displayed as histograms. The production yields have been multiplied by 106 to fit with this figure. The color code refers to the evolution temperature: the highest temperature is in black (293 K), the middle one in gray (279 K), and the lowest in white (253 K).

Chromatogram of a derivatized sample of hydrolyzed tholins in water-ammonia (5 wt % NH3—hydrolysis n°1) mixture at 293 K after derivatization. Positions of the chromatographic peaks of compounds identified by their standards are shown by arrows. The 30 m × 0.25 mm × 0.25 μm Zebron Inferno ZB-5HT column, operated in the split mode (1:25), was programmed at 6°C min−1 from 70°C to 300°C with an inlet flow of helium of 1 mL min−1.

Histograms presenting the production yields of detected organics by gas chromatography–mass spectrometry after free-oxygen tholin hydrolysis in a water-ammonia (5 wt %) mixture without (hydrolysis n°1—top) and with (hydrolysis n°3—bottom) carbonates. Yields have been multiplied by 106 for clarity.
For hydrolysis n°3, only data for two temperatures are reported. The data obtained for the lowest evolution temperature are not reported because they are not reliable due to technical issues during the gas chromatography–mass spectrometry analysis.
Table 8 sums up the production yield values of each hydrolysis product displayed in Fig. 5 but without multiplication factor in order to facilitate the reading of Fig. 5.
The black dashed box represents compounds only produced in basic media (NH3) without carbonates (hydrolysis n°1). The gray lined box represents compounds only produced in the presence of carbonates in a basic media (hydrolysis n°3).
ND = not detected. D = detected but no quantification possible.
Several astrobiologically interesting organics were characterized after free-oxygen tholin hydrolysis. The production of amino acids, carboxylic acids, and fatty acids is observed for both hydrolyses. For hydrolysis n°1, 12 products were identified: alanine, glycine, urea, α-aminobutyric acid, γ-aminobutyric acid, N-acetylglycine, β-alanine, oxalic acid, succinic acid, methylsuccinic acid, palmitic acid, and stearic acid. For hydrolysis n°3, 13 products, though several are different from the previous list, were also identified: alanine, glycine, aspartic acid, urea, α-aminobutyric acid, β-alanine, oxalic acid, succinic acid, methylsuccinic acid, palmitic acid, glycerol, homoserine, and crotonic acid. Moreover, we can see that all production yields are lower in the presence of carbonates (hydrolysis n°3), except for urea at 279 K and methylsuccinic acid for both temperatures. Indeed, for hydrolysis n°1, the production yields range from 10−1% to 10−3%, while for hydrolysis n°3 it goes from 1% to 10−5% for the compounds detected in both hydrolyses.
As can be seen in Fig. 5 and Table 8, some compounds are specific of a particular hydrolysis. Stearic acid, γ-aminobutyric acid, and N-acetylglycine are only produced in basic media (NH3) without carbonates (hydrolysis n°1—black dashed box in Table 8 and Fig. 6). While aspartic acid, glycerol, homoserine, and crotonic acid are only produced in the presence of carbonates in a basic media (hydrolysis n°3—gray lined box in Table 8 and Fig. 6).

Organics identified after the hydrolysis of free-oxygen tholins with or without carbonates. The black dashed box represents compounds only produced in basic media (NH3) without carbonates (hydrolysis n°1). The gray lined box represents compounds only produced in the presence of carbonates (hydrolysis n°3).
Table 8 and Fig. 5 show the influence of the temperature on the formation of the majority of the detected molecules. At low temperatures, glycine and succinic acid are not produced in the presence of carbonates, and on the contrary, methylsuccinic acid is not produced in the absence of carbonates. On the other hand, α-aminobutyric acid is not produced at low temperatures in the presence or absence of carbonates. This suggests that this compound cannot be formed at the surface of Titan considering its very low temperature (around 94 K).
In the absence of carbonates (hydrolysis n°1), for most of the compounds, there is a correlation between the hydrolysis temperature and the production yields of organics. For these compounds, the formation reaction is endothermic, endergonic, or both. On the contrary, in the presence of carbonates (hydrolysis n°3), for most compounds, when the hydrolysis temperature increases, the production yields of the hydrolysis products decrease. Thus in the presence of carbonates, for several compounds, the formation reaction becomes exothermic, exergonic, or both. Under the temperature range studied, alanine and β-alanine are an exception. Indeed, temperature does not seem to have any influence. Even though the temperature is not the only parameter impacting the rate of reaction, commonly when the evolution temperature increases, the rate of reaction increases, so the production yield increases. In the presence of carbonates, this is not the case, and we do not have an endothermic/endergonic reaction but an exothermic/exergonic one. This is why it is possible to assume that carbonates could involve another formation pathway.
3.2. Comparison with other works
Other teams have studied different hydrolysis conditions of Titan tholins: • Acidic pH hydrolysis with hydrochloric acid (HCl) at 6 mol L−1 (Khare et al., 1986; Nguyen, 2007; Taniuchi et al., 2013) and at 2 mol L−1 (Nguyen, 2007). • Neutral pH hydrolysis (pH = 7) (Nguyen, 2007). • Alkaline pH hydrolysis with water-ammonia 25 wt % (Ramirez et al., 2010; Poch et al., 2012), 13 wt % (Neish et al., 2009; Ramirez et al., 2010), or 6.25 and 3.125 wt % (Ramirez et al., 2010).
Table 9 sums up the qualitative identification results from both hydrolyses investigated in the present study and from other hydrolysis experiments cited above. These studies differ in their hydrolysis experimental conditions but also in the kind of Titan tholins used (produced by different experimental devices). However, in these works, the analytical techniques used were able to detect most of the hydrolysis products identified in the present study. Nevertheless, Neish et al. (2009) had a different purpose for their study that did not involve the identification of the hydrolysis products. This is why the results related to the work of Neish et al. (2009) are not reported in Table 9. Moreover, Poch et al. (2012) and Ramirez et al. (2010) conducted their studies on tholins synthesized with the experimental device PLASMA before its improvement to avoid air contamination, though they used the same general experimental conditions that were used in the present study. This is why the results related to these authors' works are not reported in Table 9 as well.
Most of the compounds identified in our study have also been detected in other studies as displayed in Table 9. The production of alanine, glycine, urea, α-aminobutyric acid, γ-aminobutyric acid, N-acetylglycine, β-alanine, aspartic acid, oxalic acid, succinic acid, and methylsuccinic acid does not seem to be specific of the experimental conditions carried out in the present study. On the contrary, palmitic acid, stearic acid, homoserine, crotonic acid, and glycerol were identified for the first time after Titan tholin hydrolysis. Concerning palmitic acid, the latter is produced with or without carbonates, so it appears that carbonates do not impact their formation, while the water-ammonia (5 wt %) mixture does. However, glycerol, homoserine, and crotonic acid are produced only in the presence of carbonates (hydrolysis n°3). So it seems that the main factor in the formation of these three compounds might be the presence of carbonates, which reveals the noticeable carbonate impact on hydrolysis product.
4. Discussion
4.1. Astrobiological consequences
Previous studies on the hydrolysis of Titan's tholins have shown that they can react efficiently in alkaline solutions, even at low temperature, to produce complex organic compounds. The incorporation of O atoms provided by water molecules into the tholins appears faster in the presence of ammonia (Neish et al., 2009), producing various molecules of biological interest, such as asparagine, aspartic acid, glutamine, and glutamic acid (Neish et al., 2010). Ramirez et al. (2010) showed that the hydrolysis rate strongly decreases with temperature but that even at low temperatures a large variety of amino acids are still detectable. Similar results were obtained by Poch et al. (2012), who observed that the main product is urea. All these works suggest that alkaline pH hydrolysis of Titan's tholins plays a catalytic role in the formation of astrobiologically interesting compounds. The present study, using a smaller and realistic ammonia concentration, shows that the diversity of organic products, including compounds of astrobiological interest, can be largely enhanced by the presence of carbonates in the starting aqueous solution.
In all our experiments, urea was the main hydrolysis product both in the presence and absence of carbonates. Moreover, several amino acids, carboxylic acids, and fatty acids were produced as well. Three compounds were produced only when carbonates were present in the aqueous ammonia solution—glycerol, crotonic acid, and homoserine. These compounds are of astrobiological importance because of their potential role in prebiotic chemistry. Glycerol is the starting molecule toward glycerol monoacyl derivatives, which would have been important for prebiotic membranes (Maurer et al., 2009). Crotonic acid is one of the simplest unsaturated fatty acids, which are components of lipids. Homoserine is one of the important nonprotein amino acids in prebiotic chemistry (Zaia et al., 2008).
4.2. Possible chemical pathways
The production of some of the identified compounds can be explained by the chemical composition of tholins. Indeed, understanding the chemical pathways responsible for the formation of hydrolysis product can give insights into the composition of tholins' soluble phase and, by extrapolation, of Titan aerosols' soluble phase. For example, urea could be synthesized from a carbodiimide group (-N=C=N-) via nucleophilic attack. It must be noticed that the presence of carbodiimide group in Titan tholins was suspected by Imanaka et al. (2004).
Concerning the production of amino acids in an alkaline pH environment, the pathway most often considered is the Strecker reaction. But it begins with an oxygenated molecule, an aldehyde (RCHO), which cannot be available in the starting clean and O-free tholins. Thus, we must look for non-oxygenated precursors. Another chemical pathway was proposed by Poch et al. (2012), which involves the hydrolysis of the nitrile group (CN) of an aminonitrile (RCH(NH2)CN) followed by the elimination of ammonia (NH3) as detailed in Fig. 7. Nevertheless, it is also assumed that tholins could be constituted of heteropolymers that would induce a similar, though more complex, chemical pathway concerning the amino acids' production (Thompson and Sagan, 1989) involving the reduction of an imide group followed by a hydrolysis step.

A possible chemical pathway for the formation of amino acids via alkalino pH hydrolysis of Titan tholins. Adapted from Poch et al. (2012).
For carboxylic acids (RCOOH), it is possible to apply the chemical pathway proposed for amino acid formation (Fig. 7), replacing the aminonitrile by a nitrile (RCN). The alcohol group, which is present for example in homoserine, can be obtained from the hydrolysis of an alkene (C=C). However, the presence of a carbon-carbon double bond within a molecule with a nitrile group could induce the production of several different molecules with an alcohol group. This shows that these chemical schemes must be considered very carefully.
These different chemical pathways suggest potential precursors that lead to the formation of the whole hydrolysis products in an alkaline pH environment as displayed in Table 10. This list contains mostly nitriles. This supports the idea that tholins could be made of HCN polymers or oligomers (Raulin, 2005, and references therein). Moreover, some of these potential precursors have already been identified in Titan tholins synthesized from CH4/N2 (5/95) such as aminoacetonitrile and butanedinitrile (He and Smith, 2014).
The list in Table 10 can be used as a database of target molecules or at least functional groups that can be identified in Titan tholins and—by extrapolation—in Titan aerosols by looking at the potential precursors. In this way, we could get some leads about the chemical composition of these particles. This kind of work is already the subject of many studies that use various analytical instruments such as UV spectroscopy, IR spectroscopy, UV fluorescence, elementary analysis, mass spectroscopy, pyrolysis gas chromatography and gas chromatography–mass spectrometry, and nuclear magnetic resonance spectroscopy. However, this rarely leads to the identification of the full molecular structure, except for new studies carried out with nuclear magnetic resonance spectroscopy on tholins (e.g., He et al., 2012; He and Smith, 2013, 2014).
4.3. Extrapolation to Titan
These hydrolysis experiments were conducted at temperatures between 293 and 253 K, while Titan's surface temperature is around 94 K. Could the identified compounds in this study be produced under Titan's surface conditions? Indeed, we can imagine that the chemical kinetics of all reactions leading to the formation of the compounds produced during those experiments would be much slower at Titan's surface temperature, and the contribution of several reactions involved in their synthesis could be quite different. Nevertheless, although the results need to be considered with caution, Ramirez et al. (2010) conducted hydrolysis experiments on tholins at 96 K, which is close to Titan's surface temperature, and showed that some hydrolysis products can be detected at this low temperature. Consequently, it is reasonable to assume that some of the compounds of astrobiological interest characterized in this study could also be present at the surface of Titan if there is indeed cryovolcanic activity and by extension the presence of a water-ammonia mixture.
Thus, the list in Table 10 can also be used as a database of target molecules that can be searched for at the surface of Titan, in particular in potentially cryovolcanic areas by looking at the hydrolysis products. This is the case of Sotra Facula (Sotra Patera) and Hotei Regio (Lopes et al., 2011, 2013), which are among the best candidates for cryovolcanism activity locations on Titan. Nevertheless, it is important to evaluate whether these compounds have a lifetime long enough under the satellite's surface conditions to be detected. This is why it would be necessary to carry out new experiments in order to study their potential lifetime at the surface of Titan, for example, by putting in place experiments of degradation of organic material, which will take into account the actual knowledge on energetic flux that reaches the surface. One of those energetic sources reaching the surface that could be simulated in the laboratory is cosmic rays. Recently, it has been suggested that Titan's internal ocean could be free of ammonia. Such a different composition could impact the nature and the diversity of hydrolysis products since it has been shown that hydrolysis of tholins is more efficient in the presence of ammonia than it is in pure water (Neish et al., 2009; Ramirez et al., 2010; Poch et al., 2012).
Finally, as pointed out by C. McKay (personal communication), the approach that was followed in the present work would also be of interest for the processing of Titan's tholins in the high-alkaline medium of Enceladus' ocean (Glein et al., 2015).
Footnotes
Acknowledgments
This study has been supported by grants (U27 & U45) from the French Space Agency (CNES). F. Raulin thanks ESA for its support as IDS of the Cassini-Huygens mission. The authors wish to thank Bernard Frouin for his help in the preparation of this manuscript, and the two referees for their contribution in improving the manuscript.
