Complex organic molecules in protoplanetary disks

Context. Protoplanetary disks are vital objects in star and planet formation, possessing all the material, gas and dust, which may form a planetary system orbiting the new star. Small, simple molecules have traditionally been detected in protoplanetary disks; however, in the ALMA era, we expect the molecular inventory of protoplanetary disks to significantly increase.

Aims. We investigate the synthesis of complex organic molecules (COMs) in protoplanetary disks to put constraints on the achievable chemical complexity and to predict species and transitions which may be observable with ALMA.

Methods. We have coupled a 2D steady-state physical model of a protoplanetary disk around a typical T Tauri star with a large gas-grain chemical network including COMs. We compare the resulting column densities with those derived from observations and perform ray-tracing calculations to predict line spectra. We compare the synthesised line intensities with current observations and determine those COMs which may be observable in nearby objects. We also compare the predicted grain-surface abundances with those derived from cometary comae observations.

Results. We find COMs are efficiently formed in the disk midplane via grain-surface chemical reactions, reaching peak grain-surface fractional abundances similar to 10(-6)-10(-4) that of the H nuclei number density. COMs formed on grain surfaces are returned to the gas phase via non-thermal desorption; however, gas-phase species reach lower fractional abundances than their grain-surface equivalents, similar to 10(-12)-10(-7). Including the irradiation of grain mantle material helps build further complexity in the ice through the replenishment of grain-surface radicals which take part in further grain-surface reactions. There is reasonable agreement with several line transitions of H2CO observed towards T Tauri star-disk systems. There is poor agreement with HC3(N) lines observed towards LkCa 15 and GO Tau and we discuss possible explanations for these discrepancies. The synthesised line intensities for CH3OH are consistent with upper limits determined towards all sources. Our models suggest CH3OH should be readily observable in nearby protoplanetary disks with ALMA; however, detection of more complex species may prove challenging, even with ALMA "Full Science" capabilities. Our grain-surface abundances are consistent with those derived from cometary comae observations providing additional evidence for the hypothesis that comets (and other planetesimals) formed via the coagulation of icy grains in the Sun's natal disk.

Complex organic molecules along the accretion flow in isolated and externally irradiated protoplanetary disks

(Abridged) The birth environment of the Sun will have influenced the conditions in the pre-solar nebula, including the attainable chemical complexity, important for prebiotic chemistry. The formation and distribution of complex organic molecules (COMs) in a disk around a T Tauri star is investigated for two scenarios: (i) an isolated disk, and (ii) a disk irradiated externally by a nearby massive star. The chemistry is calculated along the accretion flow from the outer disk inwards using a comprehensive network. Two simulations are performed, one beginning with complex ices and one with simple ices only. For the isolated disk, COMs are transported without major alteration into the inner disk where they thermally desorb into the gas reaching an abundance representative of the initial assumed ice abundance. For simple ices, COMs efficiently form on grain surfaces under the conditions in the outer disk. Gas-phase COMs are released into the molecular layer via photodesorption. For the irradiated disk, complex ices are also transported inwards; however, they undergo thermal processing caused by the warmer conditions in the irradiated disk which tends to reduce their abundance along the accretion flow. For simple ices, grain-surface chemistry cannot synthesise COMs in the outer disk because the necessary grain-surface radicals, which tend to be particularly volatile, are not sufficiently abundant on the grain surfaces. Gas-phase COMs are formed in the inner region of the irradiated disk via gas-phase chemistry induced by the desorption of strongly bound molecules such as methanol; hence, the abundances are not representative of the initial molecular abundances injected into the outer disk. These results suggest that the composition of comets formed in isolated disks may differ from those formed in externally irradiated disks with the latter composed of more simple ices.

The nucleus of Comet 67P/Churyumov-Gerasimenko. A new shape model and thermophysical analysis

Context. Comet 67P/Churyumov-Gerasimenko is the target of the European Space Agency Rosetta spacecraft rendez-vous mission. Detailed physical characteristation of the comet before arrival is important for mission planning as well as providing a test bed for ground-based observing and data-analysis methods. Aims: To conduct a long-term observational programme to characterize the physical properties of the nucleus of the comet, via ground-based optical photometry, and to combine our new data with all available nucleus data from the literature. Methods: We applied aperture photometry techniques on our imaging data and combined the extracted rotational lightcurves with data from the literature. Optical lightcurve inversion techniques were applied to constrain the spin state of the nucleus and its broad shape. We performed a detailed surface thermal analysis with the shape model and optical photometry by incorporating both into the new Advanced Thermophysical Model (ATPM), along with all available Spitzer 8-24 μm thermal-IR flux measurements from the literature. Results: A convex triangular-facet shape model was determined with axial ratios b/a = 1.239 and c/a = 0.819. These values can vary by as much as 7% in each axis and still result in a statistically significant fit to the observational data. Our best spin state solution has Psid = 12.76137 ± 0.00006 h, and a rotational pole orientated at Ecliptic coordinates λ = 78°(±10°), β = + 58°(±10°). The nucleus phase darkening behaviour was measured and best characterized using the IAU HG system. Best fit parameters are: G = 0.11 ± 0.12 and HR(1,1,0) = 15.31 ± 0.07. Our shape model combined with the ATPM can satisfactorily reconcile all optical and thermal-IR data, with the fit to the Spitzer 24 μm data taken in February 2004 being exceptionally good. We derive a range of mutually-consistent physical parameters for each thermal-IR data set, including effective radius, geometric albedo, surface thermal inertia and roughness fraction. Conclusions: The overall nucleus dimensions are well constrained and strongly imply a broad nucleus shape more akin to comet 9P/Tempel 1, rather than the highly elongated or "bi-lobed" nuclei seen for comets 103P/Hartley 2 or 8P/Tuttle. The derived low thermal inertia of

Main-belt Comet P/2012 T1 (PANSTARRS)

We present initial results from observations and numerical analyses aimed at characterizing the main-belt comet P/2012 T1 (PANSTARRS). Optical monitoring observations were made between 2012 October and 2013 February using the University of Hawaii 2.2 m telescope, the Keck I telescope, the Baade and Clay Magellan telescopes, Faulkes Telescope South, the Perkins Telescope at Lowell Observatory, and the Southern Astrophysical Research Telescope. The object's intrinsic brightness approximately doubles from the time of its discovery in early October until mid-November and then decreases by ~60% between late December and early February, similar to photometric behavior exhibited by several other main-belt comets and unlike that exhibited by disrupted asteroid (596) Scheila. We also used Keck to conduct spectroscopic searches for CN emission as well as absorption at 0.7 μm that could indicate the presence of hydrated minerals, finding an upper limit CN production rate of Q CN <1.5 × 1023 mol s-1, from which we infer a water production rate of Q_H_2O100 Myr and is unlikely to be a recently implanted interloper from the outer solar system, while a search for potential asteroid family associations reveals that it is dynamically linked to the ~155 Myr old Lixiaohua asteroid family. Some of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California, and the National Aeronautics and Space Administration, and made possible by the generous financial support of the W. M. Keck Foundation, the Magellan Telescopes located at Las Campanas Observatory, Chile, and the Southern Astrophysical Research (SOAR) telescope, which is a joint project of the Ministério da Ciência, Tecnologia, e Inovação (MCTI) da República Federativa do Brasil, the U.S. National Optical Astronomy Observatory (NOAO), the University of North Carolina at Chapel Hill (UNC), and Michigan State University (MSU).

Thermal properties, sizes, and size distribution of Jupiter-family cometary nuclei

We present results from SEPPCoN, an on-going Survey of the Ensemble Physical Properties of Cometary Nuclei. In this report we discuss mid-infrared measurements of the thermal emission from 89 nuclei of Jupiter-family comets (JFCs). All data were obtained in 2006 and 2007 using imaging capabilities of the Spitzer Space Telescope. The comets were typically 4-5 AU from the Sun when observed and most showed only a point-source with little or no extended emission from dust. For those comets showing dust, we used image processing to photometrically extract the nuclei. For all 89 comets, we present new effective radii, and for 57 comets we present beaming parameters. Thus our survey provides the largest compilation of radiometrically-derived physical properties of nuclei to date. We have six main conclusions: (a) The average beaming parameter of the JFC population is 1.03 ± 0.11, consistent with unity; coupled with the large distance of the nuclei from the Sun, this indicates that most nuclei have Tempel 1-like thermal inertia. Only two of the 57 nuclei had outlying values (in a statistical sense) of infrared beaming. (b) The known JFC population is not complete even at 3 km radius, and even for comets that approach to ˜2 AU from the Sun and so ought to be more discoverable. Several recently-discovered comets in our survey have small perihelia and large (above ˜2 km) radii. (c) With our radii, we derive an independent estimate of the JFC nuclear cumulative size distribution (CSD), and we find that it has a power-law slope of around -1.9, with the exact value depending on the bounds in radius. (d) This power-law is close to that derived by others from visible-wavelength observations that assume a fixed geometric albedo, suggesting that there is no strong dependence of geometric albedo with radius. (e) The observed CSD shows a hint of structure with an excess of comets with radii 3-6 km. (f) Our CSD is consistent with the idea that the intrinsic size distribution of the JFC population is not a simple power-law and lacks many sub-kilometer objects.

Outgassing Behavior of C/2012 S1 (ISON) from 2011 September to 2013 June

We report photometric observations for comet C/2012 S1 (ISON) obtained during the time period immediately after discovery (r = 6.28 AU) until it moved into solar conjunction in mid-2013 June using the UH2.2 m, and Gemini North 8 m telescopes on Mauna Kea, the Lowell 1.8 m in Flagstaff, the Calar Alto 1.2 m telescope in Spain, the VYSOS-5 telescopes on Mauna Loa Hawaii and data from the CARA network. Additional pre-discovery data from the Pan STARRS1 survey extends the light curve back to 2011 September 30 (r = 9.4 AU). The images showed a similar tail morphology due to small micron sized particles throughout 2013. Observations at submillimeter wavelengths using the James Clerk Maxwell Telescope on 15 nights between 2013 March 9 (r = 4.52 AU) and June 16 (r = 3.35 AU) were used to search for CO and HCN rotation lines. No gas was detected, with upper limits for CO ranging between 3.5-4.5 × 1027 molecules s-1. Combined with published water production rate estimates we have generated ice sublimation models consistent with the photometric light curve. The inbound light curve is likely controlled by sublimation of CO2. At these distances water is not a strong contributor to the outgassing. We also infer that there was a long slow outburst of activity beginning in late 2011 peaking in mid-2013 January (r ~ 5 AU) at which point the activity decreased again through 2013 June. We suggest that this outburst was driven by CO injecting large water ice grains into the coma. Observations as the comet came out of solar conjunction seem to confirm our models.

Irradiation of oxygen and water ices at by 4 keV singly and doubly charged ions; sputtering and molecular synthesis

The icy surfaces of dust grains in the Interstellar Medium and those of comets, satellites and Kuiper Belt Objects are continuously exposed tophoton and charged particle irradiation. These energetic particles maysputter and induce chemical changes in the ices and the underlyingsurfaces.In the present work 258 nm thick O2 and H2O ices were deposited at 10 K with the thickness measured by a laser interferometer method. Asimple model fit to the reflected laser intensity as measured by aphotodiode detector enabled the refractive index of the ices to bedetermined. The ices were then irradiated with various singly and doublycharged ions such as He+, 13C+, N+, O+ , Ar+, 13C2+, N2+ and O2+ at 4keV. The decrease in ice thickness as a function of ion dose wasmonitored by a laser interferometer and the model used to determine thesputtering yield as shown in Figure 1.In the case of O2 ice thesputtering yields increased with increasing ion mass in good agreementwith a model calculation [Fama, J, Shi, R.A Baragiola, Surface Sci.,602, 156 (2007)]. In the case of O2 ice, O2+ has a significant lowersputtering yield when compared to O+. The sputtering yields for O2 icewere found to be at least 9 times larger compared to those for H2O ice.For H2O ice the sputter yields for C, N and O ions were found todecrease with increasing mass. Doubly charged C, N and O ions which werefound to have the same sputtering yield as the singly charged ionswithin the experimental errors. A preliminary TPD study was carried outusing a QMS to detect the desorbed species from water ice afterirradiation by 6 × 10^15 ions of 13C+ and 13C2+. The formation of13CO and 13CO2 was observed with the yield of 13CO almost of a factor of100 larger than of 13CO2. This is in contrast to our earlier work whereonly CO¬2 was observed.

Sublimation-Driven Activity in Main-Belt Comet 313p/Gibbs

We present an observational and dynamical study of newly discovered main-belt comet 313P/Gibbs. We find that the object is clearly active both in observations obtained in 2014 and in precovery observations obtained in 2003 by the Sloan Digital Sky Survey, strongly suggestingthat its activity is sublimation-driven. This conclusion is supported by a photometric analysis showing an increase in the total brightness of the comet over the 2014 observing period, and dust modeling resultsshowing that the dust emission persists over at least three months during both active periods, where we find start dates for emission nolater than 2003 July 24 ± 10 for the 2003 active period and 2014 July 28 ± 10 for the 2014 active period. From serendipitous observations by the Subaru Telescope in 2004 when the object was apparently inactive, we estimate that the nucleus has an absolute R-band magnitude of HR = 17.1 ± 0.3, corresponding to aneffective nucleus radius of re ∼ 1.00 ± 0.15 km.The object’s faintness at that time means we cannot rule out the presence of activity, and so this computed radius should be consideredan upper limit. We find that 313P’s orbit is intrinsically chaotic, having a Lyapunov time of Tl = 12,000 yr and beinglocated near two three-body mean-motion resonances with Jupiter andSaturn, 11J-1S-5A and 10J+12S-7A, yet appears stable over >50 Myr in an apparent example of stable chaos. We furthermore find that 313P is the second main-belt comet, after P/2012 T1 (PANSTARRS), to belong tothe ∼155 Myr old Lixiaohua asteroid family.

Visible-Wavelength Survey of Jupiter-Family Cometary Nuclei as Part of SEPPCoN

We present observations of a statistically-significant number of Jupiter-family cometary nuclei as part of SEPPCoN (Survey of the Ensemble Physical Properties of Cometary Nuclei). We present preliminary results on distributions of albedos and shapes.

Physical Characterisation of Fast-Rotating Near-Earth Asteroids

We present VLT+VISIR mid-IR observations of fast-rotating near-Earth asteroids. Diameters and albedos are determined with thermal models. These NEAs may have unusual surface properties, e.g. from regolith transport/stripping due to the YORP effect.

Discovery of Main-Belt Comet P/2006 VW139 by Pan-STARRS1

We describe the discovery of comet-like activity in main-belt asteroid (300163) 2006 VW139 (later re-designated as Comet P/2006 VW139) by Pan-STARRS1. We also detail follow-up photometric, spectroscopic, and dynamical analyses of the object.

The Influence of Global-Selfheating on the Yarkovsky and YORP Effects

The Yarkovsky (orbital drift) and YORP (altering spin state) effects are important mechanisms governing the evolution of asteroids. We have included global-selfheating into a new model and demonstrate that it significantly affects YORP predictions.

Castalia — A Mission to a Main Belt Comet

Main Belt Comets (MBCs) are a newly identified class of solar system object, having asteroid-like orbits but sometimes exhibiting comet-like behavior. We present the case for a mission to an MBC, to be submitted to the European Space Agency.

Visible-wavelength Observations of Jupiter-family Comet Nuclei as Part of Seppcon.

We present the latest analysis and results from SEPPCoN (Survey of Ensemble Physical Properties of Cometary Nuclei). This on-going survey involves studying 100 JFCs - about 25% of the known population - at both mid-infrared and visible wave-lengths to constrain the distributions of sizes, shapes, spins, and albedos of this population. Having earlier reported results from measuring thermal emissions of our sample nuclei [1,2,3,4], we report here progress on the visible-wavelength observations that we have obtained at many ground-based facilities in Chile, Spain, and the United States. To date we have attempted observations of 91% of our sample of 100 JFCs, and at least 64 of those were successfully detected. In most cases the comets were at heliocentric distances between 3.0 and 6.5 AU so as to decrease the odds of a comet having a coma. Of the 64 detected comets, 48 were apparently bare, having no extended emission. Our datasets are further augmented by archival data and photometry from the NEAT program [5]. An important goal of SEPPCoN is to accumulate a large comprehensive set of high quality physical data on cometary nuclei in order to make accurate statistical comparisons with other minor-body populations such as Trojans, Centaurs, and Kuiper-belt objects. Information on the size, shape, spin-rate, albedo and color distributions is critical for understanding their origins and evolutionary processes affecting them.

Solar phase functions of cometary nuclei

Relatively few measurements of the solar phase function of cometary nuclei exist, despite the importance of this parameter in determining accurate sizes and its use in modeling surface properties. We make use of robotic telescopes and servicemode observing to monitor cometary nuclei over months at a time, combining intensive observations at a single epoch with regular short light-curve segments to efficiently account for brightness changes due to both nucleus rotation and changing solar phase angle. We present our latest results on comets 8P/Tuttle, 14P/Wolf, 67P/Churyumov- Gerasimenko and 110P/Hartley 3.