National Aeronautics and Space Administration

NASA Institute Annual Science Report 2016 Team Report: The Evolution of Prebiotic Chemical Complexity and the Organic Inventory of Protoplanetary Disk and Primordial Planets, NASA The Evolution of Prebiotic Chemical Complexity and the Organic Inventory of Protoplanetary Disk and Primordial Planets

Lead Institution: NASA Ames Research Center

Team Overview The Evolution of Prebiotic Chemical Complexity and the Organic Inventory of Protoplanetary Disks and Primordial Planets Team seeks a greater understanding of the chemical processes that occur at every stage in the evolution of organic chemical complexity, from quiescent regions of dense molecular clouds, through all stages of disk and planet formation, and ultimately to the materials that rain down on planets. The effort is an integrated, coherent program involving the interaction of a number of well-integrated research projects:

• Modeling and observations of protoplanetary disks • Modeling and observations of exoplanets • Laboratory studies of gas-grain chemistry • Laboratory studies of ice photolysis • Computational quantum chemistry

Principal Investigator: Scott Sandford These projects interact closely with each other so that each benefits from advances made in the others and helps guide future work. For example, the modeling of the chemistry that takes place in protostellar disks benefits from inputs provided by spectral, physical, and chemical properties of molecules determined by the laboratory and computational projects, but also provides guidance for key areas of future computational and laboratory work. Similarly, the computational studies can be used to help interpret laboratory results and extend them to additional materials or environments, while the lab results can provide confirmation of computational reaction paths.

Team Website: http://amesteam.arc.nasa.gov

NASA Astrobiology Institute 1 Annual Report 2016 2016 Executive Summary During 2016 our Team made significant progress on all aspects of our combined research. Some highlights are described briefly below. More details can be found in our individual project reports.

Disk Modeling - We have developed global protoplanetary disk evolution models that incorporate the formation of planetesimals in a 1+1D viscous evolution framework that combines the effects of photoevaporation, dust evolution, turbulence, and planetesimal formation by the streaming instability. This work shows that, as the gas is removed, regions with an over-density of dust grains are created that lead to the formation of planetesimals on relatively short timescales, i.e., planetesimal formation in the outer disk is a natural outcome of disk evolution. We have also begun incorporating gas-grain chemistry into the disk-modeling network to study transport of ices within disks.

Exoplanets - Our Team is at the leading edge of extrasolar planet discovery, and we are core participants in the long-running Lick-Carnegie Exoplanet Survey. Over the past year we assembled and published a compendium of 60,949 precision radial velocities for 1,624 nearby sun-like stars and red dwarfs. This data, all taken with the Keck-I Telescope over a period of more than 20 years, supports the detection of hundreds of exoplanets, including 60 new candidate worlds. Our planet detection efforts also extend to our own Solar System. Strong dynamical evidence points to the potential existence of “Planet Nine”, a ~10 Earth-mass world on an eccentric ~17,000-year orbit. Our work pinpoints the as-yet unseen object at a distance of 950 AU, with a sky position near RA=3 hours, Dec=0 degrees, and apparent brightness V~23.

Laboratory Studies - Significant progress has been made in several areas of our laboratory research. An apparatus for studying gas-grain chemistry has been developed and constructed. We have also completed a preliminary study of the abiotic production of sugars and sugar derivatives during the photolysis of simple ice mixtures. We have also completed studies that have demonstrated that all of the nucleobases used by modern life, as well as a host of related compounds, can also be produced by ice photolysis (Fig. 1).

Sample TIC

Mass = 279 Da

Adenine Standard

12.6 12.65 12.7 12.75 12.8 12.85 12.9 Time (min) Fig. 1. Left panel: Total-ion chromatogram (TIC) of the residue produced from a UV-irradiated

H2O:NH3:purine ice (top), single-ion chromatogram (SIC) of the same residue for mass 279 Da (middle), and SIC of the adenine standard (mass 279 Da). Right panel: Mass spectrum of the peak identified as adenine

in the residue produced from a UV-irradiated H2O:NH3:purine ice (top), an adenine standard (middle), and 15 shifted mass spectrum of adenine in a residue produced from an irradiated ice containing NH3.

NASA Astrobiology Institute 2 NASA AmesAnnual Research Report Center 2016 Computational Quantum Chemistry - In conjunction with the laboratory studies of nucleobase production, we investigated the formation of nucleobases in irradiated astrophysical ices containing and purine using computational quantum chemistry. This work demonstrated that the formation of nucleobases is energetically and kinetically favorable given the presence of one or several water molecules, i.e., it can occur in the solid state, but not in the gas phase (Fig. 2 and 3). This work also explains why some nucleobases () are made efficiently, while others () are not. In another study, electronic structure calculations were performed on + + + + C4H3 , C6H3 and C6H5 . C6H5 was found to be very stable and may be a good nucleation center for Fig. 2. Diagram illustrating the three reaction pathways from pyrimidine 1 to the growth of larger polycyclic uracil 3 and thymine 6 we studied using quantum chemical computations. aromatic hydrocarbons (PAHs), one of the most abundant families of molecules in space. We also continued our work with the Dutch Astrochemistry Network (DAN) on the computation of infrared emission from PAHs.

Team Members

Partha Bera Uma Gorti Martin Head-Gordon Gregory Laughlin Timothy Lee Andrew Mattioda Michel Nuevo Xander Tielens Fig. 3. Reaction scheme indicating the amino and hydroxyl group substitutions on the purine molecule and pathways to the nucleobases adenine and guanine Gustavo Cruz-Diaz studied with quantum chemical computations. Christopher Materese Tamar Stein

NASA Astrobiology Institute 3 Annual Report 2016 Project Reports

Planetesimal Formation and Transport of Ices We have begun incorporating gas grain chemistry in Protoplanetary Disks into the disk modeling network to study transport We developed global protoplanetary disk of ices. To that end, a gas-grain chemical network evolution models that incorporate the formation of was added to augment our current gas-phase planetesimals in a 1+1D viscous evolution framework. and photo-reaction chemistry in disk models. We combined the effects of photoevaporation, dust Dust grains continually collide and fragment or evolution, turbulence, and planetesimal formation coagulate (depending on local conditions in the by the streaming instability to study how dust disk) throughout the disk as it evolves, leading to forms planetesimals in an evolving disk. As the gas their size redistribution. Molecules preferentially is removed by photoevaporation, dust is retained freeze on small grains, which have greater total in the disk since flow gas densities are too low to cross-sectional area. Ices moreover change the local entrain dust. A steady increase in the amount of dust fragmentation threshold for destructive collisions relative to gas triggers the streaming instability due compared to bare grains, as ices are more “sticky” and to the back-reaction of the dust with the gas. Regions facilitate coagulation. We have constructed simple with an over-density of dust grains are created and semi-analytical formulations to include the effects lead to the formation of planetesimals on relatively of gas-grain chemistry and freeze-out to study the

short timescales. Our models produce ~76-153 distribution of CO and H2O ice as the disk evolves. Earth masses of planetesimals. Thus, planetesimal Parameter studies involving a few species at a time formation in the outer disk is found to be a natural are being conducted to gain a better understanding outcome of disk evolution. of ice chemistry in disks.

Fig. 4. The total mass in solids in different size ranges (color-coded) evolves with time as the disk disperses, and is shown in the figure above. The streaming instability can efficiently convert dust (<1 cm) into planetesimals (in blue), and in the model here, about 76 Earth masses of dust is converted during the disk lifetime of about 3 million years.

NASAThe Evolution Astrobiology of Prebiotic Institute Chemical Complexity and the Organic 4 NASA Ames Research Center Inventory of Protoplanetary Disks and Primordial Planets Exoplanet Studies We have continued our radial velocity surveys to In order to aid its detection, and to assess its likely detect and characterize extrasolar planets. The properties, we constructed atmospheric models that highlight of the past year’s effort has been the predict its reflection and emission spectra in the publication of the full catalog of high-precision optical and the infrared. We found that it should have velocities obtained over the past two decades with a very pure H-He atmosphere in which methane and the Keck Telescope (Butler, Vogt, Laughlin et al. other volatiles are fully condensed into deep-layer 2017). This compendium of 60,949 precision Doppler cloud. It should thus have a high albedo (and a light- measurements is accompanied by detailed spectral blue overall appearance) characterized by Rayleigh indices, tabulated host star properties for 1,624 scattering. (Fortney, Marley, Laughlin et al. 2016). nearby sun-like stars and red dwarfs, and a catalog of Through a Monte-Carlo survey, we established that sixty new candidate exoplanets, including a low-mass if it exists, Planet Nine is very likely (P~98%) in low- world orbiting Lalande 21185, the fourth nearest order mean-motion resonances with the most distant main-sequence stellar system to the Sun. Kuiper Belt Objects, including Sedna and 2012 VP 113. A large-scale supercomputer-based optimization We have also focused research on worlds much closer then allowed us to pin down its likely sky position to home. The past year has seen a flurry of activity (see accompanying Figure and Millholland and associated with Planet Nine, a proposed Solar System Laughlin 2016). planet with a period of order 17,000 years. If it exists, this world provides a novel and unexpected chance to advance the scientific goals of our NAI consortium.

Fig. 5. Planet Nine orbit and likely current location. A 3D manipulable version of this figure is available at https://smillholland. github.io/P9_Orbit/

NASA Astrobiology Institute 5 Annual Report 2016 Laboratory Studies of Gas-Grain Chemistry During 2016, Team members Drs. Mattioda and Cruz-Diaz modified a Harrick Praying Mantis Diffuse Reflectance apparatus and corresponding low temperature cell to accommodate a flowing hydrogen discharge lamp (see Fig. 6 and 7). This was accomplished through the design and fabrication of a high vacuum adapter for which a new technology report is being submitted to NASA.

The high vacuum adapter is equipped with two gas inlet ports, permitting the deposition

of multiple gas (H2O, etc.) and organic (PAH, etc.) samples onto a mineral surface for the creation of astrobiologically interesting ice and mineral mixtures. The apparatus will allow us to study the adsorption/ desorption processes of Polycyclic Aromatic Hydrocarbons (PAHs) and ice mixtures onto mineral grain surfaces. We will also be able to Fig. 6. Harrick Praying Mantis low-temperature setup with the flowing hydrogen lamp adapter. Note the two gas inlet ports on the adapter, which permit the determine any catalytic properties the grains introduction of astrobiologically relevant gases. may exhibit when irradiated via high energy UV photons (primarily Lyman alpha). This year we have measured the infrared reflectance spectra of the hematite, magnetite, silicon dioxide, graphite as well as iron (II) and iron (III) oxides. We have also performed initial studies of benzene adsorption onto several mineral surfaces.

Dr. Ana Ferreira de Barros (Professor and Researcher at the Federal Center for Technological Education – CEFET/ RJ), Rio de Janeiro, Brazil, started a 1-year Visiting Researcher position in which she'll be working with Dr. Mattioda on the photochemistry of PAHs in water ices. This research exchange is in alignment with the agreements between NASA and the Brazilian Space Agency (AEB) to enhance collaborations between the agencies. Dr. Barros has completed the infrared measurements on the photochemistry of

the PAH Coronene (C24H12) under various Fig. 7. Close-up of the low temperature Harrick Praying Mantis cell. The enclosed PAH:Water ice ratios and has started the area allows the sample to be heated and cooled under vacuum, as well as for the investigation of the photochemistry of introduction of gas mixtures onto the solid sample. The photolysis of the sample with UV photons can be monitored in real time via in-situ infrared spectroscopy. the PAH Ovalene (C34H12) in water ice. The larger sizes of these PAHs are astrophysically relevant to the PAHs found in interstellar environments.

The Evolution of Prebiotic Chemical Complexity and the Organic 6 NASA Ames Research Center Inventory of Protoplanetary Disks and Primordial Planets Laboratory Studies of Ice Photochemistry and the Production of Sugar Derivatives and Nucleobases During the past year, we studied the formation of (i) Nucleobases can be subdivided into two major sugars and their derivatives, in particular those seen categories: the (uracil, , and in meteorites, and (ii) purine-based compounds, in thymine) and the purines (adenine and guanine). particular the purine-based nucleobases, during We earlier demonstrated that UV irradiation of ices the UV irradiation of low-temperature ices under containing pyrimidine leads to the production of astrophysical conditions. uracil, cytosine, thymine, and other pyrimidine derivatives. During the past year we studied Simultaneous deposition and UV irradiation of similar ices containing purine and demonstrated

H2O:CH3OH ice mixtures was found to produce com- the formation of adenine and guanine, as plex sugars and sugar derivatives with up to 6 carbon well as other functionalized purines including atoms, including ribose (the sugar of RNA; Fig. 8), hypoxanthine, isoguanine, several aminopurines, glucose, and some of their isomers. Contamination of and 2,6-diaminopurine. The relative abundance of our samples was ruled out using two methods: (i) the photoproducts appears to be controlled by three analysis of residues produced from ices whose metha- factors in decreasing order of importance: 1) the nol was isotopically labeled with 13C, and (ii) the use of number of functional group additions required to a chiral derivatization agent to separate enantiomers. form the product; 2) the type of functional group

We also conducted experiments in which H2O:CO2 added; and 3) the position where the addition takes

and H2O:CH3OH:CO2 ices were used to see how they place. Collectively, our results on pyrimidine and matched the high abundance of sugar acids in mete- purine demonstrate that all the biological nucleobases orites. Preliminary results show a slight increase in the can form under the same astrophysical conditions. abundances of sugar acids in these residues.

Fig. 8. Identification of racemic ribose in a residue produced from the UV irradiation of an 2H O:CH3OH (2:1) ice mixture at 12 K for 18 hours. The top plots show the GC-MS chromatogram of the residue in the retention time range of where ribose elutes (left), and the mass spectrum of the first peak (right). The bottom plots show the same plots of a DL-ribose standard where L is 1.6 times more abundant than D. This confirms the identification of ribose in the residue, present in racemic amounts. The peak overlapping with D-ribose in the chromatogram of the residue (top left) corresponds to D-arabinose, an isomer of ribose.

NASA Astrobiology Institute 7 AnnualAnnual ReportReport 20162016 Fig. 9. Significant minima and transition states on the C6H5+ potential energy surface. The transition states for the ring-closure reaction are reported in red, while other transition states are reported as curved pathways.

Computational Quantum Chemistry We investigated the formation of nucleobases in characterize the structures, relative energies, and irradiated astrophysical ices containing pyrimidine spectroscopic and physical properties of the low- and purine. Quantum chemistry computations energy isomers of the azirinyl cation, which could be indicate thymine can be produced by reactions present in the interstellar medium. between pyrimidine, and hydroxyl and methyl radicals with only limited efficiency due to the Our collaboration with the Dutch astrochemistry difficulty of methylation, and only in the presence groups continued with two publications on the of a water matrix. Work estimating the efficiency of anharmonic vibrational spectra of five non-linear PAH production of the purine-based nucleobases adenine molecules. Computed anharmonic spectra of the full and guanine in gas-phase only and condensed- set of vibrational frequencies were shown to be in phase environments revealed that the formation of good agreement with new high-resolution spectra nucleobases is energetically and kinetically favorable in the C-H stretch region and with lower-resolution provided the presence of one or several water spectra that spanned the entire vibrational manifold. molecules. We collaborated with colleagues at Georgia Southern Electronic structure calculations were performed University to publish predicted rovibrational + + + on C4H3 , C6H3 and C6H5 (Fig. 9). On the basis of the spectrum and spectroscopic constants of proton calculated minima and transition states and their bound molecular complexes for possible assignment + relative energies, the likely pathways to C4H3 and in high-resolution astronomical observations and on + C6H3 likely arise from unimolecular decomposition of the nature of bidirectionally expanded closed-shell + + hot C4H4 and C6H4 by H atom elimination. In contrast, PAH and PANH anions. + C6H5 was found to be very stable to fragmentation and elimination pathways, and may be a good nucleation center for the growth of larger polycyclic aromatic hydrocarbons in interstellar conditions. Ab initio quantum chemical methods were used to

NASA Astrobiology Institute 8 Annual Report 2016 Evolution of Prebiotic Chemical Complexity and the Organic Inventory of Protoplanetary Disks and Primordial Planets: 2016 Publications

Bandyopadhyay, B., Stein, T., Fang, Y., Kostko, O., White, A., Head-Gordon, M., and Ahmed, M. (2016). Probing Ionic Complexes of Ethylene and Acetylene with Vacuum Ultraviolet Radiation. Journal of Physical Chemistry A, 120, 5053-5064. DOI: 10.1021/acs. jpca.6b00107

Batygin, K., Bodenheimer, P. and Laughlin, G. (2016). “In Situ Formation and Dynamical Evolution of Hot Jupiter Systems” Astrophysical Journal, Volume 829, Issue 2, article id. 114, 17 pp. DOI: 10.3847/0004- 637X/829/2/114

Bera, P.P., Nuevo, M., Materese, C.K., Sandford, S.A., and Lee, T.J. (2016). Mechanisms for the Formation of Thymine Under Astrophysical Conditions and Implications for the Origin of Life. Journal of Chemical Physics, 144: 144308. DOI: 10.1063/1.4945745

Fortenberry, R.C., Francisco, J. S., and Lee, T.J. (2016). Towards the Astronomical Detection of the Proton-Bound Complex NN-HCO+: Implications for the Spectra of Protoplanetary Disks. Astrophysical Journal, 819, 141. DOI: 10.3847/0004- 637X/819/2/141

Fortenberry, R.C., Lee, T.J., and Francisco, J.S. (2016). Quantum Chemical Analysis of the CO−HNN+ Proton-Bound Complex. Journal of Physical Chemistry A, 120, 7745. DOI: 10.1021/acs. jpca.6b07515

Fortenberry, R.C., Moore, M.M., and Lee, T.J. (2016). Excited State Trends in Bidirectionally Expanded Closed-Shell PAH and PANH Anions. Journal of Physical Chemistry A, 120, 7327. DOI: 10.1021/acs. jpca.6b06654

Fortenberry, R.C., Roueff, E., and Lee, T.J. (2016). Inclusion of 13C and D in Protonated Acetylene. Chemical Physics Letters, 650: 126. DOI: 10.1016/j. cplett.2016.02.068

NASA Astrobiology Institute 9 Annual Report 2016 Fortney, J. J., Marley, M. S., Laughlin, G., et al. (2016). Lieman-Sifry, J., Hughes, A. M., Carpenter, J. M., Gorti, The Hunt for Planet Nine: Atmosphere, Spectra, U., Hales, A., and Flaherty, K. M. (2016). Debris Disks Evolution and Detectability. Astrophysical Journal in the Scorpius-Centaurus OB Association Resolved Letters, Volume 824, Issue 2, article id. L25, 6 pp. DOI: by ALMA. The Astrophysical Journal, 828, 25. DOI: 10.3847/2041-8205/824/2/L25 10.3847/0004-637X/828/1/25

Gorti, U. (2016). Photoevaporation and Disk Dispersal. Mackie, C.J., Candian, A., Huang, X., Maltseva, E., Proceedings of IAU Symposium No. 314, Young Stars Petrignani, A., Oomens, J., Mattioda, A.L., Buma, and Planets Near the Sun, J.H. Kastner, B. Stelzer, S.A. W.J., Lee, T.J., and Tielens, A.G.G.M. (2016). The Metchev, eds., (Cambridge University Press), pp., anharmonic quartic force field infrared spectra of 153-158. DOI: 10.1017/S1743921315006420 five non-linear Polycyclic Aromatic Hydrocarbons: benz[a]anthracene, chrysene, phenanthrene, Kastner, J. H., Principe, D. A., Punzi, K., Stelzer, B., Gorti, pyrene, and triphenylene. Journal of Chemical U., Pascucci, I., and Argiroffi, C. (2016). M Stars in Physics, 145, 084313. DOI: 10.1063/1.4961438 the TW Hya Association: Stellar X-Rays and Disk Dissipation. The Astronomical Journal, 152, 3. DOI: Maltseva, E., Petrignani, A., Candian, A., Mackie, C. J., 10.3847/0004-6256/152/1/3 Huang, X., Lee, T.J., Tielens, A.G.G.M., Oomens, J., and Buma, W.J. (2016). High-resolution IR absorption Kastner, J. H., Qi, C., Gorti, U., Hily-Blant, P., Oberg, spectroscopy of polycyclic aromatic hydrocarbons K., Forveille, T., Andrews, S., and Wilner, D. (2016). in the 3µ region: Role of periphery. Astrophys. J. 831,

Rings of C2H in the Molecular Disks Orbiting TW 58. DOI: 10.3847/0004-637X/831/1/58 Hya and V4046 Sgr. Proceedings of IAU Symposium No. 314, Young Stars and Planets Near the Sun, J.H. Millholland, S., Wang, S. and Laughlin, G. (2016). On the Kastner, B. Stelzer, S.A. Metchev, eds., (Cambridge Detection of Non-transiting Hot Jupiters in Multiple- University Press), pp. 193-194. DOI: 10.1017/ planet Systems, Astrophysical Journal Letters, Volume S1743921315005979 823, Issue 1, article id. L7, 6 pp. DOI: 10.3847/2041- 8205/823/1/L7 Kokkila-Schumacher, S. I. L., Bera, P.P., and Lee, T.J. (2016). Characterization of the Azirinyl Cation and Peverati, R., Bera, P.P., Lee, T.J., and Head-Gordon, its Isomers. Journal of Physical Chemistry A, 120, M. (2016). Insights into Hydrocarbon Chain 1275. DOI: 10.1021/acs.jpca.5b12373 and Ring Formation in the Interstellar Medium: + Computational Study of the Isomers of C4H3 , + + Kospal, A., Abraham, P., Csengeri, T., Gorti, U., Henning, C6H3 , and C6H5 and their Formation Pathways. T., Moor, A., Semenov, D. A., Szucs, L., and Gusten, R. Astrophysical Journal, 830, 128. DOI: 10.3847/0004- (2016). Cold CO Gas in the Disk of the Young Eruptive 637X/830/2/128 Star EX Lup. The Astrophysical Journal, 821, L4. DOI: 10.3847/2041-8205/821/1/L4 Sandford, S.A., Engrand, C., and Rotundi, A. (2016). Organic Matter in Cosmic Dust. Elements 12, 185- 189. DOI: 10.2113/gselements.12.3.185

The Evolution of Prebiotic Chemical Complexity and the Organic 10 NASA Ames Research Center Inventory of Protoplanetary Disks and Primordial Planets