The Exoplanet Population Observation Simulator. I - The Inner Edges of Planetary Systems

The Kepler survey provides a statistical census of planetary systems out to the habitable zone. Because most planets are non-transiting, orbital architectures are best estimated using simulated observations of ensemble populations. Here, we introduce EPOS, the Exoplanet Population Observation Simulator, to estimate the prevalence and orbital architectures of multi-planet systems based on the latest Kepler data release, DR25. We estimate that at least 42% of sun-like stars have nearly coplanar planetary systems with 7 or more exoplanets. The fraction of stars with at least one planet within 1 au could be as high as 100% depending on assumptions about the distribution of single transiting planets. We estimate an occurrence rate of planets in the habitable zone around sun-like stars of eta_earth=36+-14%. The innermost planets in multi-planet systems are clustered around an orbital period of 10 days (0.1 au), reminiscent of the protoplanetary disk inner edge or could be explained by a planet trap at that location. Only a small fraction of planetary systems have the innermost planet at long orbital periods, with fewer than ~8% and ~3% having no planet interior to the orbit of Mercury and Venus, respectively. These results reinforce the view that the solar system is not a typical planetary system, but an outlier among the distribution of known exoplanetary systems. We predict that at least half of the habitable zone exoplanets are accompanied by (non-transiting) planets at shorter orbital periods, hence knowledge of a close-in exoplanet could be used as a way to optimize the search for Earth-size planets in the Habitable Zone with future direct imaging missions.Comment: Accepted in AAS journals, code available on githu

Earths in Other Solar Systems N-body simulations: the Role of Orbital Damping in Reproducing the Kepler Planetary Systems

The population of exoplanetary systems detected by Kepler provides opportunities to refine our understanding of planet formation. Unraveling the conditions needed to produce the observed exoplanets will sallow us to make informed predictions as to where habitable worlds exist within the galaxy. In this paper, we examine using N-body simulations how the properties of planetary systems are determined during the final stages of assembly. While accretion is a chaotic process, trends in the ensemble properties of planetary systems provide a memory of the initial distribution of solid mass around a star prior to accretion. We also use EPOS, the Exoplanet Population Observation Simulator, to account for detection biases and show that different accretion scenarios can be distinguished from observations of the Kepler systems. We show that the period of the innermost planet, the ratio of orbital periods of adjacent planets, and masses of the planets are determined by the total mass and radial distribution of embryos and planetesimals at the beginning of accretion. In general, some amount of orbital damping, either via planetesimals or gas, during accretion is needed to match the whole population of exoplanets. Surprisingly, all simulated planetary systems have planets that are similar in size, showing that the "peas in a pod" pattern can be consistent with both a giant impact scenario and a planet migration scenario. The inclusion of material at distances larger than what Kepler observes has a profound impact on the observed planetary architectures, and thus on the formation and delivery of volatiles to possible habitable worlds.Comment: Resubmitted to ApJ. Planet formation models available online at http://eos-nexus.org/genesis-database

The Exoplanet Population Observation Simulator. II -- Population Synthesis in the Era of Kepler

The collection of planetary system properties derived from large surveys such as Kepler provides critical constraints on planet formation and evolution. These constraints can only be applied to planet formation models, however, if the observational biases and selection effects are properly accounted for. Here we show how epos, the Exoplanet Population Observation Simulator, can be used to constrain planet formation models by comparing the Bern planet population synthesis models to the Kepler exoplanetary systems. We compile a series of diagnostics, based on occurrence rates of different classes of planets and the architectures of multi-planet systems, that can be used as benchmarks for future and current modeling efforts. Overall, we find that a model with 100 seed planetary cores per protoplanetary disk provides a reasonable match to most diagnostics. Based on these diagnostics we identify physical properties and processes that would result in the Bern model more closely matching the known planetary systems. These are: moving the planet trap at the inner disk edge outward; increasing the formation efficiency of mini-Neptunes; and reducing the fraction of stars that form observable planets. We conclude with an outlook on the composition of planets in the habitable zone, and highlight that the majority of simulated planets smaller than 1.7 Earth radii have substantial hydrogen atmospheres. The software used in this paper is available online for public scrutiny at https://github.com/GijsMulders/eposComment: Accepted in Ap

The spectral energy distribution of galaxies at z > 2.5: Implications from the Herschel/SPIRE color-color diagram

We use the Herschel SPIRE color-color diagram to study the spectral energy distribution (SED) and the redshift estimation of high-z galaxies. We compiled a sample of 57 galaxies with spectroscopically confirmed redshifts and SPIRE detections in all three bands at $z=2.5-6.4$, and compared their average SPIRE colors with SED templates from local and high-z libraries. We find that local SEDs are inconsistent with high-z observations. The local calibrations of the parameters need to be adjusted to describe the average colors of high-z galaxies. For high-z libraries, the templates with an evolution from z=0 to 3 can well describe the average colors of the observations at high redshift. Using these templates, we defined color cuts to divide the SPIRE color-color diagram into different regions with different mean redshifts. We tested this method and two other color cut methods using a large sample of 783 Herschel-selected galaxies, and find that although these methods can separate the sample into populations with different mean redshifts, the dispersion of redshifts in each population is considerably large. Additional information is needed for better sampling.Comment: 17 pages, 14 figures, accepted for publication in A&

Studying the evolution of galaxies in compact groups over the past 3 Gyr - II. The importance of environment in the suppression of star formation

We present an in depth study on the evolution of galaxy properties in compact groups over the past 3 Gyr. We are using the largest multi-wavelength sample to-date, comprised 1770 groups (containing 7417 galaxies), in the redshift range of 0.01<z<0.23. To derive the physical properties of the galaxies we rely on ultraviolet (UV)-to-infrared spectral energy distribution modeling, using CIGALE. Our results suggest that during the 3 Gyr period covered by our sample, the star formation activity of galaxies in our groups has been substantially reduced (3-10 times). Moreover, their star formation histories as well as their UV-optical and mid-infrared colors are significantly different from those of field and cluster galaxies, indicating that compact group galaxies spend more time transitioning through the green valley. The morphological transformation from late-type spirals into early-type galaxies occurs in the mid-infrared transition zone rather than in the UV-optical green valley. We find evidence of shocks in the emission line ratios and gas velocity dispersions of the late-type galaxies located below the star forming main sequence. Our results suggest that in addition to gas stripping, turbulence and shocks might play an important role in suppressing the star formation in compact group galaxies.Comment: (Accepted for publication in MNRAS, date of submission November 18, 2015

Lava channel formation during the 2001 eruption on Mount Etna: evidence for mechanical erosion

We report the direct observation of a peculiar lava channel that was formed near the base of a parasitic cone during the 2001 eruption on Mount Etna. Erosive processes by flowing lava are commonly attributed to thermal erosion. However, field evidence strongly suggests that models of thermal erosion cannot explain the formation of this channel. Here, we put forward the idea that the essential erosion mechanism was abrasive wear. By applying a simple model from tribology we demonstrate that the available data agree favorably with our hypothesis. Consequently, we propose that erosional processes resembling the wear phenomena in glacial erosion are possible in a volcanic environment.Comment: accepted for publication in Physical Review Letter

Welcome to the Twilight Zone: The Mid-Infrared Properties of Poststarburst Galaxies

We investigate the optical and Wide-field Survey Explorer (WISE) colors of "E+A" identified post-starburst galaxies, including a deep analysis on 190 post-starbursts detected in the 2{\mu}m All Sky Survey Extended Source Catalog. The post-starburst galaxies appear in both the optical green valley and the WISE Infrared Transition Zone (IRTZ). Furthermore, we find that post-starbursts occupy a distinct region [3.4]-[4.6] vs. [4.6]-[12] WISE colors, enabling the identification of this class of transitioning galaxies through the use of broad-band photometric criteria alone. We have investigated possible causes for the WISE colors of post-starbursts by constructing a composite spectral energy distribution (SED), finding that mid-infrared (4-12{\mu}m) properties of post-starbursts are consistent with either 11.3{\mu}m polycyclic aromatic hydrocarbon emission, or Thermally Pulsating Asymptotic Giant Branch (TP-AGB) and post-AGB stars. The composite SED of extended post- starburst galaxies with 22{\mu}m emission detected with signal to noise >3 requires a hot dust component to produce their observed rising mid-infrared SED between 12 and 22{\mu}m. The composite SED of WISE 22{\mu}m non-detections (S/N<3), created by stacking 22{\mu}m images, is also flat, requiring a hot dust component. The most likely source of this mid-infrared emission of these E+A galaxies is a buried active galactic nucleus. The inferred upper limit to the Eddington ratios of post-starbursts are 1e-2 to 1e-4, with an average of 1e-3. This suggests that AGNs are not radiatively dominant in these systems. This could mean that including selections able to identify active galactic nuclei as part of a search for transitioning and post-starburst galaxies would create a more complete census of the transition pathways taken as a galaxy quenches its star formation.Comment: 13 pages, 11 figures, accepted for publication in the Astrophysical Journa

CO Depletion in Protoplanetary Disks: A Unified Picture Combining Physical Sequestration and Chemical Processing

The gas-phase CO abundance (relative to hydrogen) in protoplanetary disks decreases by up to 2 orders of magnitude from its ISM value ${\sim}10^{-4}$, even after accounting for freeze-out and photo-dissociation. Previous studies have shown that while local chemical processing of CO and the sequestration of CO ice on solids in the midplane can both contribute, neither of these processes appears capable of consistently reaching the observed depletion factors on the relevant timescale of $1{-}3\mathrm{~Myr}$. In this study, we model these processes simultaneously by including a compact chemical network (centered on carbon and oxygen) to 2D ($r+z$) simulations of the outer ($r>20\mathrm{~au}$) disk regions that include turbulent diffusion, pebble formation, and pebble dynamics. In general, we find that the CO/H$_2$ abundance is a complex function of time and location. Focusing on CO in the warm molecular layer, we find that only the most complete model (with chemistry and pebble evolution included) can reach depletion factors consistent with observations. In the absence of pressure traps, highly-efficient planetesimal formation, or high cosmic ray ionization rates, this model also predicts a resurgence of CO vapor interior to the CO snowline. We show the impact of physical and chemical processes on the elemental (C/O) and (C/H) ratios (in the gas and ice phases), discuss the use of CO as a disk mass tracer, and, finally, connect our predicted pebble ice compositions to those of pristine planetesimals as found in the Cold Classical Kuiper Belt and debris disks.Comment: Accepted for publication in The Astrophysical Journa

A hyper luminous starburst at z=4.72 magnified by a lensing galaxy pair at z=1.48

International audienceWe serendipitously discovered in the Herschel Reference Survey an extremely bright infrared source with S500 ∼ 120 mJy in the line of sight of the Virgo cluster which we name Red Virgo 4 (RV4). Based on IRAM/EMIR and IRAM/NOEMA detections of the CO(5−4), CO(4−3), and [CI] lines, RV4 is located at a redshift of 4.724, yielding a total observed infrared luminosity of 1.1 ± 0.6 × 1014 L⊙. At the position of the Herschel emission, three blobs are detected with the VLA at 10 cm. The CO(5−4) line detection of each blob confirms that they are at the same redshift with the same line width, indicating that they are multiple images of the same source. In Spitzer and deep optical observations, two sources, High-z Lens 1 (HL1) West and HL1 East, are detected at the center of the three VLA/NOEMA blobs. These two sources are placed at z = 1.48 with X-shooter spectra, suggesting that they could be merging and gravitationally lensing the emission of RV4. HL1 is the second most distant lens known to date in strong lensing systems. Constrained by the position of the three VLA/NOEMA blobs, the Einstein radius of the lensing system is 2.2″ ± 0.2 (20 kpc). The high redshift of HL1 and the large Einstein radius are highly unusual for a strong lensing system. In this paper, we present the insterstellar medium properties of the background source RV4. Different estimates of the gas depletion time yield low values suggesting that RV4 is a starburst galaxy. Among all high-z submillimeter galaxies, this source exhibits one of the lowest L[CI] to LIR ratios, 3.2 ± 0.9 × 10−6, suggesting an extremely short gas depletion time of only 14 ± 5 Myr. It also shows a relatively high L[CI] to LCO(4−3) ratio (0.7 ± 0.2) and low LCO(5−4) to LIR ratio (only ∼50% of the value expected for normal galaxies) hinting at low density of gas. Finally, we discuss the short depletion time of RV4. It can be explained by either a very high star formation efficiency, which is difficult to reconcile with major mergers simulations of high-z galaxies, or a rapid decrease of star formation, which would bias the estimate of the depletion time toward an artificially low value

Towards understanding the relation between the gas and the attenuation in galaxies at kpc scales

[abridged] Aims. The aim of the present paper is to provide new and more detailed relations at the kpc scale between the gas surface density and the face-on optical depth directly calibrated on galaxies, in order to compute the attenuation not only for semi-analytic models but also observationally as new and upcoming radio observatories are able to trace gas ever farther in the Universe. Methods. We have selected a sample of 4 nearby resolved galaxies and a sample of 27 unresolved galaxies from the Herschel Reference Survey and the Very Nearby Galaxies Survey, for which we have a large set of multi-wavelength data from the FUV to the FIR including metallicity gradients for resolved galaxies, along with radio HI and CO observations. For each pixel in resolved galaxies and for each galaxy in the unresolved sample, we compute the face-on optical depth from the attenuation determined with the CIGALE SED fitting code and an assumed geometry. We determine the gas surface density from HI and CO observations with a metallicity-dependent XCO factor. Results. We provide new, simple to use, relations to determine the face-on optical depth from the gas surface density, taking the metallicity into account, which proves to be crucial for a proper estimate. The method used to determine the gas surface density or the face-on optical depth has little impact on the relations except for galaxies that have an inclination over 50d. Finally, we provide detailed instructions on how to compute the attenuation practically from the gas surface density taking into account possible information on the metallicity. Conclusions. Examination of the influence of these new relations on simulated FUV and IR luminosity functions shows a clear impact compared to older oft-used relations, which in turn could affect the conclusions drawn from studies based on large scale cosmological simulations.Comment: 24 pages, 21 figures, accepted for publication in A&