Introduction
Astronomical studies have shown several small bodies, presently located in or originating from the outer Solar System, having surfaces covered with material of extremely low geometric albedo that is either neutral (black) or dark red in colour (Cruikshank Reference Cruikshank1987). Two independent analyses from the earliest of such research initiatives, one using pyrolysis–gas chromatography–mass spectrometry (Py–GC–MS) on a fallen Orgueil (C1) chondrite (Bandurski & Nagy Reference Bandurski and Nagy1976) and a second investigating astronomical reflectance spectra of the Centaur 944 Hidalgo (Gradie & Veverka Reference Gradie and Veverka1980), proposed the possible presence of complex organic material akin to terrestrial kerogen, although of non-biogenic origin, on the surfaces of small Solar System objects. These observations were supported by the argument that complex organic solids could not only be synthesized on inner Solar System asteroids by Fischer–Tropsch-like reactions (Hayatsu et al. Reference Hayatsu, Studier and Anders1971; Studier et al. Reference Studier, Hayatsu and Anders1972; Hayatsu & Anders Reference Hayatsu and Anders1981), but also on interstellar grain surfaces (Sagan & Khare Reference Sagan and Khare1979). Overall the research carried out to date, in forms of laboratory simulations, astronomical observations and analyses of extraterrestrial material therefore, indicate that complex organics could possibly be the reservoirs of prebiotic species and they could be abundant throughout the Solar System (Kerridge Reference Kerridge1999; Kwok Reference Kwok2009).
The term ‘tholins’ refers to a wide class of brownish–blackish highly viscous organic residues formed from irradiation of cosmically relevant gases (Sagan & Khare Reference Sagan and Khare1979; Cable et al. Reference Cable, Hörst, Hodyss, Beauchamp, Smith and Willis2011). Tholins are synthesized with several gaseous precursors and under varied experimental conditions. These diverse materials and conditions generate distinct organic residues but with certain overarching properties. Such distinct tholins simulate macromolecular organic material on Saturn's moon Titan, icy moons such as Triton, primitive Earth, and Centaurs, comets and trans-Neptunian objects (TNOs). Tholins synthesized by tesla coil electrical discharge under conditions that simulate the N2 and methane (CH4) dominant atmosphere of Titan (Khare et al. Reference Khare1984) have been a common ingredient of mixtures, used to fit the astronomical near-infrared (NIR) reflectances of Centaurs and TNOs. Centaurs whose astronomical NIR spectra are modelled with mixtures containing Titan tholins include 5145 Pholus (Cruikshank et al. Reference Cruikshank1998), 10199 Chariklo (Dotto et al. Reference Dotto, Barucci, Leyrat, Romon, de Bergh and Licandro2003a), 52872 Okyrhoe (DeMeo et al. Reference DeMeo, Barucci, Merlin, Guilbert-Lepoutre, Alvarez-Candal, Delsanti, Fornasier and de Bergh2010), 8405 Asbolus (Barucci et al. Reference Barucci, de Bergh, Cuby, Le Bras, Schmitt and Romon2000), (55576) 2002 GB10 and (83982) 2002 GO9 (Dorresoundiram et al. Reference Dorresoundiram, Barucci, Tozzi, Poulet, Boehnhardt, de Bergh and Peixinho2005), and the TNOs include 90482 Orcus, (73480) 2002 PN34* (DeMeo et al. Reference DeMeo, Barucci, Merlin, Guilbert-Lepoutre, Alvarez-Candal, Delsanti, Fornasier and de Bergh2010), (47171) 1999 TC36 (Dotto et al. Reference Dotto, Barucci, Boehnhardt, Romon, Dorresoundiram, Peixinho, de Bergh and Lazzarin2003b) and (55565) 2002 AW197 (Dorresoundiram et al. Reference Dorresoundiram, Barucci, Tozzi, Poulet, Boehnhardt, de Bergh and Peixinho2005).
Despite the essentiality of Titan tholins in simulation mixtures, their laboratory synthesis conditions do not represent the environment on Centaurs and TNOs, particularly the volatile precursors, their compositions, temperature, pressure and radiation flux regimes. Scattered attempts to synthesize tholins approximately representing the organic material on icy Solar System bodies have been undertaken. These tholins were designated differently as Triton or ice tholins (Khare et al. Reference Khare, Thompson, Chyba, Arakawa and Sagan1989, Reference Khare, Thompson, Cheng, Chyba, Sagan, Arakawa, Meisse and Tuminello1993; McDonald et al. Reference McDonald, Thompson, Heinrich, Khare and Sagan1994, Reference McDonald, Whited, DeRuiter, Khare, Patnaik and Sagan1996); however, until now not much analytical progress has been made regarding these tholins. The single-most factor for this research limitation has been their low production amounts as compared with Titan tholins. Yet, along with Titan tholins, Triton tholins have also been a regular component in the geographical mixtures simulating the reflectance of several Centaurs and TNOs.
The motivation of our investigation is to comprehend the chemical composition of Titan tholins, especially the components that allow the spectral matching of Centaurs and TNOs. This consistent use makes Titan tholins a plausible analogue to the surface chemical compositions of Centaurs and TNOs; however this plausibility demanded investigations beyond the NIR studies applied generally. We elucidate the material properties that contribute to the extremely low geometric albedo (0.05–0.15) by investigating the chemical structure and composition of Titan tholins using multiple analytical techniques such as Py–GC–MS, laser desorption–time-of-flight–mass spectrometry (LD–TOF–MS), Raman spectroscopy and field-emission scanning electron microscopy (FE-SEM).
Materials and methods
Production of tholins
The Titan tholins (hereafter ‘tholins’) used in this study were synthesized at the NASA Ames Research Center, USA. Following the description by McKay (Reference McKay1996), a gas mixture of 10% CH4 in N2 was introduced in a 3 litres glass reactor at a gas flow rate of ~3 cm3 min–1. The pressure within the reactor was maintained at slightly above ambient. Three sets of tungsten electrodes, each with an arc-gap of 0.5 cm, were inserted into the reactor. They were energized via electrical discharge from a tesla coil, where each coil was cyclically switched on for 15 min. The entire apparatus was operated for 3 weeks at ambient temperature, which eventually led to accumulation of tholins on the inner walls of the reactor. The accumulated tholins were scraped carefully and packed in air-tight vials.
The elemental composition of tholins produced in these simulations has been reported to be 67.2% C, 5.79% H, 14.2% N and 12% O (McKay Reference McKay1996). Oxygen was not present in the initial gas mixture and is believed to be an oxidative contaminant present throughout the tholins chemical structure. On discarding O from the mass balance the stoichiometry calculated was C11H11N2 (McKay Reference McKay1996).
Py–GC–MS
Approximately 1 mg of raw tholins with 100 mg of pure SiO2 was crushed to powder. The powdered mixture was loaded on a platinum filament of the Pyrola® 2000 pyrolyser system. Small and approximately equal quantities of raw tholins samples were pyrolysed at several temperatures – 175, 275, 375, 475 and 575°C – consecutively. The filament was heated up to the set pyrolysis temperature for 2 s. The volatilized species were analysed using a Varian 3800 GC (Varian, Inc.) equipped with a Restek® Rtx-20 chromatographic column (30 m × 0.25 mm × 0.15 µm) and coupled to a Varian 4000 ion-trap MS (Varian, Inc.) with internal ionization. The column temperature was ramped up from 30°C at a rate of 15°C min−1 to reach to 250°C where it was maintained constant during the remainder GC run. The total duration of the GC run was 25 min. Further the chemical species separated by the GC were identified by the MS using the scan detection mode at a rate of 1.5 scans s−1 and detection range of m/z 35–1000.
LD–TOF–MS
Raw tholin powder was directly pressed onto a customized removable sample stub that was analysed by LD–TOF–MS (Bruker autoflex™ speed, Bruker Daltonics). Furthermore, 5 mg of raw tholin powder was dissolved in 100 µl n-hexane, of which ~3 µl n-hexane-extracted solution was pipetted onto a highly polished stainless steel multi-target plate for mass analysis. The LD–TOF–MS instrument was equipped with a frequency-tripled Nd:YAG (355 nm) commercial laser running at 50 Hz, with an energy around 200–250 µJ, a pulse width of ~3 ns, and laser beam size of ~200 × 200 µm2. A sum of 1500 laser shots (spectra) was collected with the laser rastering over the sample surface. Both positive and negative ion modes were available on this instrument. The tholins caused a strong background signal in the positive mode measurements. However, in the negative mode the results showed more identifiable peaks. Therefore, only the negative ion mode data were considered in this work.
Raman spectroscopy
The instrument used was a LabRAM-HR 800 UV Raman spectrometer (Horiba Jobin–Yvon) attached to an Olympus BX41 microscope equipped with an Olympus MPlane 100× objective.
The Raman spectrometer was equipped with a 488 nm-line Ar+ laser (106020B0S, Melles Griot, Carlsbad, CA, USA), with a laser power of 20 mW and a focal length of 800 mm. The laser was dispersed by a 600 litres mm−1 grating on a charge-coupled device detector with 1024 × 256 pixels, producing a spectral resolution of 0.43 cm−1. For internal calibration a silicon standard with a major peak at 520.4 cm–1 was used. The spectra recording and processing were done with the LabSpec™ version 5.19.17 (Horiba Jobin–Yvon, Villeneuve d'Ascq, France). Background subtraction was applied using manual line fitting. Band fitting for the two main bands was conducted with a Gauss–Lorentz function from the background-corrected spectrum.
FE-SEM
The FE-SEM used was a LEO 1530 Gemini (Carl Zeiss) instrument with an accelerating voltage of 15 kV. The raw tholins samples were mounted on carbon pads and then coated with 15 nm thickness of Pt/Au with a sputter coater. No wet brush was used; the resultant dust was removed with air.
Results and discussion
Py–GC–MS analyses
The tholins pyrolysate chromatogram at 575°C displays aromatic hydrocarbons (Fig. 1) that were identified by linear ion-trap MS based on their fragmentation pattern. Of the numerous chromatograms obtained at different pyrolysis temperatures, the chromatogram at 575 °C, owing to higher temperature, showed the highest diversity of polycyclic aromatic hydrocarbons (PAHs). Therefore, we will only discuss this chromatogram in the following. The GC separated a variety of nitrile species that are typical of Titan tholins pyrolysate such as 2-butenenitrile (McGuigan et al. Reference McGuigan, Waite, Imanaka and Sacks2006), benzonitrile (Coll et al. Reference Coll, Coscia, Gazeau, Guez and Raulin1998; McGuigan et al. Reference McGuigan, Waite, Imanaka and Sacks2006; Kawai et al. Reference Kawai, Jagota, Kaneko, Obayashi, Yoshimura, Khare, Deamer, McKay and Kobayashi2013) and 2-napthalenecarbonitrile (Sagan et al. Reference Sagan, Khare, Thompson, McDonald, Wing, Bada, Vo-Dinh and Arakawa1993). We identified several aromatic hydrocarbons with unfused and fused rings (see Fig. 1 and Table 1).

Fig. 1. Total ion current chromatogram of raw tholins pyrolysed at 575°C. Labelled species are volatilized aromatic nitrile, unfused and fused aromatic hydrocarbons.
Table 1. Aromatic hydrocarbons detected in raw tholins using Py–GC–MS.

SNR, signal-to-noise ratio. Bold values signifies Molecular ion peaks.
Imanaka et al. (Reference Imanaka, Khare, Elsila, Bakes, McKay, Cruikshank, Sugita, Matsui and Zare2004) pointed out the PAH abundance in certain types of tholins (Khare et al. Reference Khare1984; Sagan et al. Reference Sagan, Khare, Thompson, McDonald, Wing, Bada, Vo-Dinh and Arakawa1993; Ehrenfreund et al. Reference Ehrenfreund, Boon, Commandeur, Sagan, Thompson and Khare1995), whereas deficiency in other types were reported (Coll et al. Reference Coll, Coscia, Smith, Gazeau, Ramírez, Cernogara, Israël and Raulin1999). PAH-abundant Titan tholins were consistently produced in reactors with low pressures (13–26 Pa), relatively high CH4 concentration (>1%), and with tesla- and inductively coupled cold plasma discharges (Hodyss et al. Reference Hodyss, McDonald, Sarker, Smith, Beauchamp and Beauchamp2004; Trainer et al. Reference Trainer, Pavlov, Jimenez, McKay, Worsnop, Toon and Tolbert2004; Cable et al. Reference Cable, Hörst, Hodyss, Beauchamp, Smith and Willis2011). In contrast, PAH-deficient Titan tholins have been characteristically produced in reactors maintained at high pressures (~100 Pa), low temperatures (100–150 K) and with capacitively coupled cold plasma discharge (Ramírez et al. Reference Ramírez, Navarro-González, Coll and Raulin2005; Quirico et al. Reference Quirico2008).
The tholins in this study were produced at low pressure, with tesla discharge, and with continuous gas flow and high CH4 content in the starting mixture. These conditions have been shown to correlate with high PAH production (Trainer et al. Reference Trainer, Pavlov, Jimenez, McKay, Worsnop, Toon and Tolbert2004; Cable et al. Reference Cable, Hörst, Hodyss, Beauchamp, Smith and Willis2011). Trainer et al. (Reference Trainer, Pavlov, Jimenez, McKay, Worsnop, Toon and Tolbert2004) reported that particles produced from a mixture of 10% CH4 in N2 are consistent with a large fraction of aromatics, including specific m/z peaks likely due to PAHs. However, at lower concentrations of CH4 (<1%), the mass fraction of PAHs greatly diminishes and an aliphatic chemical formation pathway dominates. Benzene (C6H6) and toluene (C6H5CH3) identified in the pyrolysate are known precursors of PAHs, and their detection among the other volatilized PAHs direct that they were volatilized from a large PAH (LPAH) structural component existing within the greater chemical structure of tholins.
LD–TOF–MS analyses
We obtained mass spectra of both raw and n-hexane extracted tholins, and the latter was used as a blank measurement. The spectrum of n-hexane-extracted tholins was also compared with that of pure n-hexane which was used as solvent blank reference (Fig. 2). The comparison clearly distinguishes the mass peaks contributed by the organic solvent n-hexane (Fig. 2(a)) with those contributed by tholins (Fig. 2(b)).

Fig. 2. Negative ion LD–TOF mass spectra of (a) pure n-hexane (solvent blank) and (b) n-hexane-extracted tholins.
Figure 3 shows two individual m/z spectral segments of negative ion mode measurements of n-hexane-extracted tholins. The segment from m/z 20–100 of n-hexane-extracted tholins (Fig. 3(a)) displays regularly spaced peak clusters, each separated by ∆ m/z 12 corresponding to the mass of atomic carbon. The clusters were observed from m/z 24 to 96 and could be assigned to negatively ionized C2 to C8. Although the structural conformations and the formation mechanisms of the assigned species could not be construed, we postulate these peaks arise due to carbon vapours formed from LD of a highly carbonaceous tholins. Many ion species that we tentatively assigned to the observed peaks are known precursors to soot (Cochran Reference Cochran1987; Gerhardt et al. Reference Gerhardt, Löffler and Homann1987; Ebert Reference Ebert1990; Gingerich et al. Reference Gingerich, Finkbeiner and Schmude1994) and mass spectrometric studies have already demonstrated emission of negative ionized C2 to C8 on molecular sublimation of graphite (Honig Reference Honig1954). The species that could represent these peaks have also been observed astronomically. Diatomic carbon (C2) was detected in the jets of comet C/1989 X1 (Austin) (Suzuki et al. Reference Suzuki, Kurihara and Watanabe1990); triatomic carbon (C3) was detected in comets and circumstellar discs (Douglas Reference Douglas1951; Goebel et al. Reference Goebel, Bregman, Strecker, Witteborn and Erickson1978; Cernicharo et al. Reference Cernicharo, Goicoechea and Caux2000; Haddad et al. Reference Haddad, Zhao, Linnartz and Ubachs2013). Linear C4 was first identified and characterized in gas phase, when it was produced by ultraviolet laser vaporization of graphite (Heath & Saykally Reference Heath and Saykally1991). Both C4 and C5 are speculated in carbon-rich evolved stars (IRC+10216) and in the molecular cloud Sagittarius B2 (Bernath et al. Reference Bernath, Hinkle and Keady1989; Cernicharo et al. Reference Cernicharo, Goicoechea and Benilan2002). We also observed a strong peak at m/z 63 that could be tentatively assigned to C5H3 −. C5H3 − is an anticipated species in interstellar clouds (Herbst & Leung Reference Herbst and Leung1989) and acts as a precursor to fullerenes, PAH and soot (Bettens & Herbst Reference Bettens and Herbst1997; Richter & Howard Reference Richter and Howard2000; Hansen et al. Reference Hansen, Klippenstein, Miller, Wang, Cool, Law, Westmoreland, Kasper and Höinghaus2006).

Fig. 3. Negative ion LD–TOF mass spectra of raw tholins at lower m/z range. (a) Clusters of regularly spaced mass peaks (∆m/z 12) in the m/z range of 20–100. (b) Intense peak at m/z 474 that could be tentatively assigned to C38H18 −.
In the next segment of the n-hexane-extracted tholins mass spectrum (Fig. 3(b)) we observed a cluster with a strong peak at m/z 474 and minor peaks at 475 and 476. With confirmed reproducibility of these peaks, we speculate that the strongest peak (m/z 474) may be a potential LPAH with a tentative molecular formula C38H18 −. Semi-empirical studies state that C38H18 − is a conformationally stable LPAH (Pogodin & Agranat Reference Pogodin and Agranat2001) and has also been detected in high resolution mass spectral studies of coke synthesized from pyrolysis of benzene (Kousoku et al. Reference Kousoku, Ashida, Miyasato, Miyake and Miura2014). However, we were unable to comprehend the reason for the resolution of this peak cluster and the absence of other nearby peaks with comparable intensities.
Figure 4 shows two individual m/z spectral segments of negative ion mode measurements of raw tholins along with the entire spectrum (Fig. 4(c)). In the raw tholins mass spectrum (Fig. 4(a)) we observed peaks at m/z 696, 708, 924 and 936, which could be tentatively estimated for C58 −, C59 −, C77 − and C78 −, respectively. Furthermore, in the next segment of the same mass spectrum (Fig. 4(b)) mass peaks at m/z 1153, 1177 and 1417 could be tentatively estimated as average molecular masses for C96 −, C98 − and C118 −, respectively.

Fig. 4. Negative ion LD–TOF mass spectra of raw tholins at higher m/z range. (a) Segment of spectrum displaying mass peak clusters at around m/z 696, 708, 924 and 936 that could be tentatively assigned to C58 −, C59 −, C77 − and C78 −, respectively. (b) Segment displaying peak clusters at m/z 1153, 1177 and 1417 that could be tentatively assigned to C96 −, C98 − and C118 −, respectively. (c) Entire mass spectrum of raw tholins including the segments (a) and (b).
Raman spectroscopy analyses
The variations in the carbonaceous structure of tholins were studied using Raman spectroscopy. The Raman spectrum revealed two prominent features, the D band at 1357 cm–1 and G band at 1576 cm–1 (Fig. 5). These bands are a straightforward characteristic of disordered or defected graphitic material (Vidano & Fischbach Reference Vidano and Fischbach1978). The D band has an A1g symmetry in first-order scattering progression that arises due to defects in the sp 2 bonded carbon graphite structure (Ferrari & Robertson Reference Ferrari and Robertson2000). These defects are probably due to in plane substituted H and N heteroatoms and vacancies (Gupta & Saxena Reference Gupta and Saxena2009) within the tholins. The G band is due to tangential stretching of C–C bonds having an E2g symmetry.

Fig. 5. Raman spectrum (λ 488 nm) decomposition into three Gaussian curves (coloured red, yellow and green) for raw tholins.
The Ag(2) symmetry at around 1462 cm–1, known to be the most prominent Raman feature exhibited by C60 fullerene (Gupta & Saxena Reference Gupta and Saxena2009) is not clearly visible in our Raman spectrum (Fig. 5). Also, the LD–TOF–MS did not reveal a mass peak at m/z 720 that could be assigned to C60 (Fig. 4(a)). The peak at ~2950 cm–1 is identified as the D + D′ band. The broad shape of this peak is characteristic for defected graphite structures (Nemanich & Solin Reference Nemanich and Solin1977; Sato et al. Reference Sato, Kamo and Setaka1978).
The D + D′ band displays a shouldered peak at ~3250 cm−1, characteristic to ordered graphite (Nemanich & Solin Reference Nemanich and Solin1977; Sato et al. Reference Sato, Kamo and Setaka1978) identified as 2D′ band. Taken together, the D + D′ and 2D′ bands suggest that the tholins structure probably consists of both disordered and ordered graphite.
FE-SEM analyses
The FE-SEM images (Fig. 6) of raw tholins display two distinct carbonaceous structures – fractal-like aggregates of soot nanoparticles (Fig. 6(a)) and cauliflower-type graphitic globules (Fig. 6(b)). The nanoparticles monomers do not appear to be uniform in size; the largest are ~100 nm in diameter and they appear to be growing in a three-dimensional fractal manner forming aggregates. The loosely packed cauliflower-type graphitic globules have non-uniform diameters ranging from 100 to 1000 nm. These globules appear concentric and discontinuous; the latter probably indicating diversity in the structural conformations of the nucleation centres. Both the aggregates and the globules were possibly formed from different precursor LPAHs and with different fullerenes acting as nucleation centre. But, the exact synthesis mechanism of carbonized material from LPAHs is yet to be understood.

Fig. 6. High-resolution FE-SEM images of raw tholins displaying two distinct carbonaceous structures: (a) fractal-like soot particle aggregates and (b) cauliflower-like graphitic globules.
Our observations from various analytical techniques now enable us to conceive a ‘proof of concept’ that N2:CH4 (9:1) gaseous mixture, when subjected to tesla coil discharge, forms extensive PAH assemblages and, with longer duration of irradiation, eventually undergoes carbonization. The presence of H and N heteroatoms in the starting mixture and the oxidation contributed O atoms should be potentially forming defects in the graphite. However, certain segments of tholins have seemingly formed ordered graphite, indicating complete elimination of heteroatoms. The regularly spaced mass peaks from m/z 24 to 96 and the mass peak clusters observed at higher m/z hint at the presence of graphitic structure. LD of graphite is known to generate fullerenes and so the identification of native fullerenes within graphitic material by LD–TOF–MS have for long been argued (Kroto et al. Reference Kroto, Heath, O'Brien, Curl and Smalley1985; Hammond & Zare Reference Hammond and Zare2008).
In this study, however, we did not intend to determine the nativity of fullerenes within tholins. The occurrence of high molecular mass peaks in our LD–TOF–MS data (Fig. 4), which could be tentatively assigned to fullerenes, directs to the existence of graphite components within the larger tholins structure. Since we were unable to identify any Raman signatures indicative for fullerenes, we preliminarily infer that the tholins apparently do not contain innate fullerenes. Nevertheless, we do not completely rule out the possible role of ‘curved’ fullerenes as an intermediate transition structure between the ‘planar’ LPAH and ‘curved’ primary soot particles, the latter being observed in the FE-SEM. This assumption stands on a well-established research that PAHs contribute to the generation of soot (Faccinetto et al. Reference Faccinetto, Desgroux, Ziskind, Therssen and Focsa2011). Therefore, we stress on the need for further research in exploring the characteristics of any tholins innate fullerenes.
Carbonaceous material exhibits a variety of chemical bond hybridizations (sp, sp 2 and sp 3) that in turn cause characteristic optical features. In graphite, the C–C bonds are trigonal, sp 2 hybridized and form extensive planar structures with sixfold rings. The planar structure leads to complete π-electron delocalization thereby imparting the low albedo. The N heteroatom in disordered graphite of tholins could further assist the π-electron delocalization due to its lone pair of electrons (Cruikshank et al. Reference Cruikshank, Imanaka and Dalle Ore2005). The low geometric albedo of comets and other similar SSSBs could possibly be arising from graphitization of tholin-like complex organic material present on their surfaces.
The potential presence of graphitic tholins-like material on Centaurs, TNOs and comets is of great astrobiological significance. Such material possibly could have comprised a large fraction of the organic payload of the ~1.8 × 1023 g of cometary material delivered during the late heavy bombardment (~4.1–3.8 Ga) to the primitive Earth (Gomes et al. Reference Gomes, Levison, Tsiganis and Morbidelli2005). Tholin-like material could be a candidate source of prebiotic amino acids (Khare et al. Reference Khare, Sagan, Ogino, Nagy, Er, Schram and Arakawa1986) and graphite; the latter probably acting as a template for the synthesis of a primordial genetic architecture (Sowerby et al. Reference Sowerby, Cohn, Heckl and Holm2001).
Conclusions
In astronomical models of several Centaurs and TNOs, Titan tholins were presumed to simulate only the NIR signature while amorphous carbon within the same mixtures was assumed to simulate the low albedo. Based on our observations on Titan tholins, we postulate that the low geometric albedo must not necessarily be due to a separate source of amorphous carbon. Alternatively, it can also result from the gradual carbonization of the same complex organic material that contributes to the NIR signature. Although Titan tholins may not be entirely representative of organic material on Centaurs and TNOs, the overarching common properties of complex organic materials makes them close analogues. Therefore, the structural and chemical properties of complex organic materials on Centaurs and TNOs have to be further explored.
Acknowledgements
This paper is dedicated to the late Dr Bishun Khare, who passed away on August 20, 2013. Khare made noteworthy contribution to the global effort to understand tholins and encouraged initiating this study as well. We thank Dr Hiroshi Imanaka, David Beeler and Dr Seema Jagota for discussion of tholins plasma-discharge experimental setup. We thank Dr Martin Hilchenbach for the discussions on preliminary Raman analyses of tholins. C.G. acknowledges the International Max Planck Research School for the Doctoral and Postdoctoral Fellowships awarded by the Max Planck Institute for Solar System Research, Germany. T.G. was supported by an appointment to the NASA Postdoctoral Programme at the Goddard Space Flight Center, administered by Oak Ridge Associated Universities through a contract with NASA. The authors would like to thank the reviewers, Dr Rainer Oswald and Dr Tim Leefmann for their valuable remarks and suggestions to improve this paper.