What can we tell about exoplanet habitability if currently only the stellar properties, planet radius, and the incoming stellar flux are known? A planet is in the Habitable Zone (HZ) if it harbors liquid water on its surface. The HZ is traditionally conceived as a sharp region around stars because it is calculated for one planet with specific properties. Such an approach is limiting because the planets' atmospheric and geophysical properties, which influence the presence of liquid water on the surface, are currently unknown but expected to be diverse. A statistical HZ description is outlined which does not favor one planet type. Instead the stellar and planet properties are treated as random variables and a continuous range of planet scenarios are considered. Various probability density functions are assigned to each random variable, and a combination of Monte Carlo sampling and climate modeling is used to generate synthetic exoplanet populations with known surface climates. Then, the properties of the liquid water bearing subpopulation is analyzed. Given our current observational knowledge, the HZ takes the form of a weakly-constrained but smooth probability function. The HZ has an inner edge but a clear outer edge is not seen. Currently only optimistic upper limits can be derived for the potentially observable HZ occurrence rate. Finally, we illustrate through an example how future data on exoplanet atmospheres will help to narrow down the probability that an exoplanet harbors liquid water and we identify the greatest observational challenge in the way of finding a habitable exoplanet.

High-precision eclipse spectrophotometry of transiting terrestrial exoplanets represents a promising path for the first atmospheric characterizations of habitable worlds and the search for life outside our solar system. The detection of terrestrial planets transiting nearby late-type M-dwarfs could make this approach applicable within the next decade, with soon-to-come general facilities. In this context, we previously identified GJ 1214 as a high-priority target for a transit search, as the transit probability of a habitable planet orbiting this nearby M4.5 dwarf would be significantly enhanced by the transiting nature of GJ 1214 b, the super-Earth already known to orbit the star. Based on this observation, we have set up an ambitious high-precision photometric monitoring of GJ 1214 with the Spitzer Space Telescope to probe the inner part of its habitable zone in search of a transiting planet as small as Mars. We present here the results of this transit search. Unfortunately, we did not detect any other transiting planets. Assuming that GJ 1214 hosts a habitable planet larger than Mars that has an orbital period smaller than 20.9 days, our global analysis of the whole Spitzer dataset leads to an a posteriori no-transit probability of ~ 98%. Our analysis allows us to significantly improve the characterization of GJ 1214 b, to measure its occultation depth to be 70+-35 ppm at 4.5 microns, and to constrain it to be smaller than 205ppm (3-sigma upper limit) at 3.6 microns. In agreement with the many transmission measurements published so far for GJ 1214 b, these emission measurements are consistent with both a metal-rich and a cloudy hydrogen-rich atmosphere.

We present new visible and infrared observations of the hot Jupiter Kepler-7b to determine its atmospheric properties. Our analysis allows us to 1) refine Kepler-7b's relatively large geometric albedo of Ag=0.35+-0.02, 2) place upper limits on Kepler-7b thermal emission that remains undetected in both Spitzer bandpasses and 3) report a westward shift in the Kepler optical phase curve. We argue that Kepler-7b's visible flux cannot be due to thermal emission or Rayleigh scattering from H2 molecules. We therefore conclude that high altitude, optically reflective clouds located west from the substellar point are present in its atmosphere. We find that a silicate-based cloud composition is a possible candidate. Kepler-7b exhibits several properties that may make it particularly amenable to cloud formation in its upper atmosphere. These include a hot deep atmosphere that avoids a cloud cold trap, very low surface gravity to suppress cloud sedimentation, and a planetary equilibrium temperature in a range that allows for silicate clouds to potentially form in the visible atmosphere probed by Kepler. Our analysis does not only present evidence of optically thick clouds on Kepler-7b but also yields the first map of clouds in an exoplanet atmosphere.

We propose to use Kepler in 2-wheel mode to conduct a detailed search for Earth-sized planets orbiting ultra-cool stars and brown dwarfs (spectral types from M7 to L3). This population of objects presents several advantages for exoplanet surveys. First, ultra-cool stars and brown dwarfs are small and thus result in favorable planet-to-star area ratios. Second, because of their low effective temperature, the inner edge of their habitable zone is extremely close (2 to 3 days only). Third, our targets are bright at infrared wavelengths, which will enable detailed follow-up studies. Our program therefore represents a unique opportunity to find a transiting Earth-size exoplanet for which atmospheric features (including biosignatures) could be detected with near-to-come facilities such as JWST. Such exoplanet has not been discovered yet. Kepler in 2-wheel mode provides the required stability and photometric precision to make this survey successful. Our initial target sample includes 60 ultra-cool stars and brown dwarfs from which we expect to detect at least one transiting planet. We propose to monitor each source for 4 days, resulting in a total program duration of ~240 days.

We explore the minimum distance from a host star where an exoplanet could potentially be habitable in order not to discard close-in rocky exoplanets for follow-up observations. We find that the inner edge of the Habitable Zone for hot desert worlds can be as close as 0.38 AU around a solar-like star, if the greenhouse effect is reduced ($\sim$ 1% relative humidity) and the surface albedo is increased. We consider a wide range of atmospheric and planetary parameters such as the mixing ratios of greenhouse gases (water vapor and CO$_2$), surface albedo, pressure and gravity. Intermediate surface pressure ($\sim$1-10 bars) is necessary to limit water loss and to simultaneously sustain an active water cycle. We additionally find that the water loss timescale is influenced by the atmospheric CO$_2$ level, because it indirectly influences the stratospheric water mixing ratio. If the CO$_2$ mixing ratio of dry planets at the inner edge is smaller than 10$^{-4}$, the water loss timescale is $\sim$1 billion years, which is considered here too short for life to evolve. We also show that the expected transmission spectra of hot desert worlds are similar to an Earth-like planet. Therefore, an instrument designed to identify biosignature gases in an Earth-like atmosphere can also identify similarly abundant gases in the atmospheres of dry planets. Our inner edge limit is closer to the host star than previous estimates. As a consequence, the occurrence rate of potentially habitable planets is larger than previously thought.

Context. Transition disks typically appear in resolved millimeter observations as giant dust rings surrounding their young host stars. More accurate observations with ALMA have shown several of these rings to be in fact asymmetric: they have lopsided shapes. It has been speculated that these rings act as dust traps, which would make them important laboratories for studying planet formation. It has been shown that an elongated giant vortex produced in a disk with a strong viscosity jump strikingly resembles the observed asymmetric rings. Aims. We aim to study a similar behavior for a disk in which a giant planet is embedded. However, a giant planet can induce two kinds of asymmetries: (1) a giant vortex, and (2) an eccentric disk. We studied under which conditions each of these can appear, and how one can observationally distinguish between them. This is important because only a vortex can trap particles both radially and azimuthally, while the eccentric ring can only trap particles in radial direction. Methods. We used the FARGO code to conduct the hydro-simulations. We set up a disk with an embedded giant planet and took a radial grid spanning from 0.1 to 7 times the planet semi-major axis. We ran the simulations with various viscosity values and planet masses for 1000 planet orbits to allow a fully developed vortex or disk eccentricity. Afterwards, we compared the dust distribution in a vortex-holding disk with an eccentric disk using dust simulations. Results. We find that vorticity and eccentricity are distinguishable by looking at the azimuthal contrast of the dust density. While vortices, as particle traps, produce very pronounced azimuthal asymmetries, eccentric features are not able to accumulate millimeter dust particles in azimuthal direction, and therefore the asymmetries are expected to be modest.

We present model atmospheres for an Earth-like planet orbiting the entire grid of main sequence FGK stars with effective temperatures ranging from Teff = 4250K to Teff = 7000K in 250K intervals. We model the remotely detectable spectra of Earth-like planets for clear and cloudy atmospheres at the 1AU equivalent distance from the VIS to IR (0.4 {\mu}m - 20 {\mu}m) to compare detectability of features in different wavelength ranges in accordance with JWST and future design concepts to characterize exo-Earths. We also explore the effect of the stellar UV levels as well as spectral energy distribution on a terrestrial atmosphere concentrating on detectable atmospheric features that indicate habitability on Earth, namely: H2O, O3, CH4, N2O and CH3Cl. The increase in UV dominates changes of O3, OH, CH4, N2O and CH3Cl whereas the increase in stellar temperature dominates changes in H2O. The overall effect as stellar effective temperatures and corresponding UV increase, is a lower surface temperature of the planet due to a bigger part of the stellar flux being reflected at short wavelengths, as well as increased photolysis. Earth-like atmospheric models show more O3 and OH but less stratospheric CH4, N2O, CH3Cl and tropospheric H2O (but more stratospheric H2O) with increasing effective temperature of Main Sequence stars. The corresponding spectral features on the other hand show different detectability depending on the wavelength observed. We concentrate on directly imaged planets here as framework to interpret future lightcurves, direct imaging and secondary eclipse measurements of atmospheres of terrestrial planets in the HZ at varying orbital positions.

One significant difference between the atmospheres of stars and exoplanets is the presence of condensed particles (clouds or hazes) in the atmosphere of the latter. The main goal of this paper is to develop a self-consistent microphysical cloud model for 1D atmospheric codes, which can reproduce some observed properties of Earth, such as the average albedo, surface temperature, and global energy budget. The cloud model is designed to be computationally efficient, simple to implement, and applicable for a wide range of atmospheric parameters for planets in the habitable zone. We use a 1D, cloud-free, radiative-convective, and photochemical equilibrium code originally developed by Kasting, Pavlov, Segura, and collaborators as basis for our cloudy atmosphere model. The cloud model is based on models used by the meteorology community for Earth's clouds. The free parameters of the model are the relative humidity and number density of condensation nuclei, and the precipitation efficiency. In a 1D model, the cloud coverage cannot be self-consistently determined, thus we treat it as a free parameter. We apply this model to Earth (aerosol number density 100 cm^-3, relative humidity 77 %, liquid cloud fraction 40 %, and ice cloud fraction 25 %) and find that a precipitation efficiency of 0.8 is needed to reproduce the albedo, average surface temperature and global energy budget of Earth. We perform simulations to determine how the albedo and the climate of a planet is influenced by the free parameters of the cloud model. We find that the planetary climate is most sensitive to changes in the liquid water cloud fraction and precipitation efficiency. The advantage of our cloud model is that the cloud height and the droplet sizes are self-consistently calculated, both of which influence the climate and albedo of exoplanets.

We investigate dust growth due to settling in a 1D vertical column of a protoplanetary disk. It is known from the observed 10 micron feature in disk SEDs, that small micron-sized grains are present at the disk atmosphere throughout the lifetime of the disk. We hope to explain such questions as what process can keep the disk atmospheres dusty for the lifetime of the disk and how does the particle properties change as a function of height above the midplane. We use a Monte Carlo code to follow the mass and porosity evolution of the particles in time. The used collision model is based on laboratory experiments performed on dust aggregates. As the experiments cannot cover all possible collision scenarios, the largest uncertainty of our model is the necessary extrapolations we had to perform. We simultaneously solve for the particle growth and motion. Particles can move vertically due to settling and turbulent mixing. We assume that the vertical profile of the gas density is fixed in time and only the solid component evolves. We find that the used collision model strongly influences the masses and sizes of the particles. The laboratory experiment based collision model greatly reduces the particle sizes compared to models that assume sticking at all collision velocities. We find that a turbulence parameter of alpha = 10^-2 is needed to keep the dust atmospheres dusty, but such strong turbulence can produce only small particles at the midplane which is not favorable for planetesimal formation models. We also see that the particles are larger at the midplane and smaller at the upper layers of the disk. At 3-4 pressure scale heights micron-sized particles are produced. These particle sizes are needed to explain the 10 micron feature of disk SEDs. Turbulence may therefore help to keep small dust particles in the disk atmosphere.

The thesis deals with the first stage of planet formation, namely dust coagulation from micron to millimeter sizes in circumstellar disks. For the first time, we collect and compile the recent laboratory experiments on dust aggregates into a collision model that can be implemented into dust coagulation models. We put this model into a Monte Carlo code that uses representative particles to simulate dust evolution. Simulations are performed using three different disk models in a local box (0D) located at 1 AU distance from the central star. We find that the dust evolution does not follow the previously assumed growth-fragmentation cycle, but growth is halted by bouncing before the fragmentation regime is reached. We call this the bouncing barrier which is an additional obstacle during the already complex formation process of planetesimals. The absence of the growth-fragmentation cycle and the halted growth has two important consequences for planet formation. 1) It is observed that disk atmospheres are dusty throughout their lifetime. Previous models concluded that the small, continuously produced fragments can keep the disk atmospheres dusty. We however show that small fragments are not produced because bouncing prevents fragmentation. 2) As particles do not reach the fragmentation barrier, their sizes are smaller compared to the sizes reached in previous dust models. We decided to investigate point 1) in more detail. A vertical column of a disk (1D) is modeled including the sedimentation of the particles. We find that already intermediate levels of turbulence can prevent particles settling to the midplane.

We examine the first phase of the core accretion model, namely the dust growth/fragmentation in binary systems. In our model, a gas and dust disk is present around the primary star and is perturbed by the secondary. We study the effects of a secondary with/without eccentricity on the dust population to determine what sizes the aggregates can reach and how that compares to the dust population in disks around single stars. We find that the secondary star has two effects on the dust population. 1.) The disk is truncated due to the presence of the secondary star and the maximum mass of the particles is decreased in the lowered gas densities. This effect is dominant in the outer disk. 2.) The perturbation of the secondary pumps up the eccentricity of the gas disk, which in turn increases the relative velocity between the dust and the gas. Therefore the maximum particle sizes are further decreased. The second effect of the secondary influences the entire disk. Coagulation is efficiently reduced even at the very inner parts of the disk. The average mass of the particles is reduced by four orders of magnitude (as a consequence, the stopping time is reduced by one order of magnitude) in disks around binary systems compared to dust in disks around single stars.

Context. Observations at sub-millimeter and mm wavelengths will in the near future be able to resolve the radial dependence of the mm spectral slope in circumstellar disks with a resolution of around a few AU at the distance of the closest star-forming regions. Aims. We aim to constrain physical models of grain growth and fragmentation by a large sample of (sub-)mm observations of disks around pre-main sequence stars in the Taurus-Auriga and Ophiuchus star-forming regions. Methods. State-of-the-art coagulation/fragmentation and disk-structure codes are coupled to produce steady-state grain size distributions and to predict the spectral slopes at (sub-)mm wavelengths. Results. This work presents the first calculations predicting the mm spectral slope based on a physical model of grain growth. Our models can quite naturally reproduce the observed mm-slopes, but a simultaneous match to the observed range of flux levels can only be reached by a reduction of the dust mass by a factor of a few up to about 30 while keeping the gas mass of the disk the same. This dust reduction can either be due to radial drift at a reduced rate or during an earlier evolutionary time (otherwise the predicted fluxes would become too low) or due to efficient conversion of dust into larger, unseen bodies.

The sticking of micron sized dust particles due to surface forces in circumstellar disks is the first stage in the production of asteroids and planets. The key ingredients that drive this process are the relative velocity between the dust particles in this environment and the complex physics of dust aggregate collisions. Here we present the results of a collision model, which is based on laboratory experiments of these aggregates. We investigate the maximum aggregate size and mass that can be reached by coagulation in protoplanetary disks. We model the growth of dust aggregates at 1 AU at the midplane at three different gas densities. We find that the evolution of the dust does not follow the previously assumed growth-fragmentation cycles. Catastrophic fragmentation hardly occurs in the three disk models. Furthermore we see long lived, quasi-steady states in the distribution function of the aggregates due to bouncing. We explore how the mass and the porosity change upon varying the turbulence parameter and by varying the critical mass ratio of dust particles. Particles reach Stokes numbers of roughly 10^-4 during the simulations. The particle growth is stopped by bouncing rather than fragmentation in these models. The final Stokes number of the aggregates is rather insensitive to the variations of the gas density and the strength of turbulence. The maximum mass of the particles is limited to approximately 1 gram (chondrule-sized particles). Planetesimal formation can proceed via the turbulent concentration of these aerodynamically size-sorted chondrule-sized particles.

The growth processes from protoplanetary dust to planetesimals are not fully understood. Laboratory experiments and theoretical models have shown that collisions among the dust aggregates can lead to sticking, bouncing, and fragmentation. However, no systematic study on the collisional outcome of protoplanetary dust has been performed so far so that a physical model of the dust evolution in protoplanetary disks is still missing. We intend to map the parameter space for the collisional interaction of arbitrarily porous dust aggregates. This parameter space encompasses the dust-aggregate masses, their porosities and the collision velocity. With such a complete mapping of the collisional outcomes of protoplanetary dust aggregates, it will be possible to follow the collisional evolution of dust in a protoplanetary disk environment. We use literature data, perform own laboratory experiments, and apply simple physical models to get a complete picture of the collisional interaction of protoplanetary dust aggregates. In our study, we found four different types of sticking, two types of bouncing, and three types of fragmentation as possible outcomes in collisions among protoplanetary dust aggregates. We distinguish between eight combinations of porosity and mass ratio. For each of these cases, we present a complete collision model for dust-aggregate masses between 10^-12 and 10^2 g and collision velocities in the range 10^-4 to 10^4 cm/s for arbitrary porosities. This model comprises the collisional outcome, the mass(es) of the resulting aggregate(s) and their porosities. We present the first complete collision model for protoplanetary dust. This collision model can be used for the determination of the dust-growth rate in protoplanetary disks.

As accretion in protoplanetary disks is enabled by turbulent viscosity, the border between active and inactive (dead) zones constitutes a location where there is an abrupt change in the accretion flow. The gas accumulation that ensues triggers the Rossby wave instability, that in turn saturates into anticyclonic vortices. It was suggested that the trapping of solids within them leads to a burst of planet formation on very short timescales. We perform two-dimensional global simulations of the dynamics of gas and solids in a non-magnetized thin protoplanetary disk with the Pencil Code. We use multiple particle species of radius 1, 10, 30, and 100 cm, solving for the particles' gravitational interaction by a particle-mesh method. The dead zone is modeled as a region of low viscosity. Adiabatic and locally isothermal equations of state are used. We find that the Rossby wave instability is triggered under a variety of conditions, thus making vortex formation a robust process. Inside the vortices, fast accumulation of solids occurs and the particles collapse into objects of planetary mass in timescales as short as five orbits. Because the drag force is size-dependent, aerodynamical sorting ensues within the vortical motion, and the first bound structures formed are composed primarily of similarly-sized particles. In addition to erosion due to ram pressure, we identify gas tides from the massive vortices as a disrupting agent of formed protoplanetary embryos. We also estimate the collisional velocity history of the particles that compose the most massive embryo by the end of the simulation, finding that the vast majority of them never experienced a collision with another particle at speeds faster than 1 m/s.

Context: There is increasing need for good algorithms for modeling the aggregation and fragmentation of solid particles (dust grains, dust aggregates, boulders) in various astrophysical settings, including protoplanetary disks, planetary- and sub-stellar atmospheres and dense molecular cloud cores. Here we describe a new algorithm that combines advantages of various standard methods into one. Aims: The aim is to develop a method that 1) can solve for aggregation and fragmentation, 2) can easily include the effect and evolution of grain properties such as compactness, composition, etc., and 3) can be built as a coagulation/fragmentation module into a hydrodynamics simulations. Methods: We develop a Monte-Carlo method in which we follow the 'life' of a limited number of representative particles. Each of these particles is associated with a certain fraction of the total dust mass and thereby represents a large number of true particles which all are assumed to have the same properties as their representative particle. Under the assumption that the total number of true particles vastly exceeds the number of representative particles, the chance of a representative particle colliding with another representative particle is negligibly small, and we therefore ignore this possibility. This now makes it possible to employ a statistical approach to the evolution of the representative particles. Results: The method reproduces the known analytic solutions of simplified coagulation kernels, and compares well to numerical results for Brownian motion using other methods. For reasonably well-behaved kernels it produces good results even for moderate number of swarms.