Limb darkening is a fundamental ingredient for interpreting observations of planetary transits, eclipsing binaries, optical/infrared interferometry and microlensing events. However, this modeling traditionally represents limb darkening by a simple law having one or two coefficients that have been derived from plane-parallel model stellar atmospheres, which has been done by many researchers. More recently, researchers have gone beyond plane-parallel models and considered other geometries. We previously studied the limb-darkening coefficients from spherically symmetric and plane-parallel model stellar atmospheres for cool giant and supergiant stars, and in this investigation we apply the same techniques to FGK dwarf stars. We present limb-darkening coefficients, gravity-darkening coefficients and interferometric angular diameter corrections from Atlas and SAtlas model stellar atmospheres. We find that sphericity is important even for dwarf model atmospheres, leading to significant differences in the predicted coefficients.

In my thesis I investigate the occurrence of mid-infrared excess around stars and their evolution. Since the launch of the first infrared satellite, IRAS, we have known that a large fraction of stars exhibit significant levels of infrared emission above their predicted photospheric level. Resolved optical and infrared images have revealed the majority of these excesses to arise from circumstellar disk structures, made up of distributions of planetesimals, rocks, and dust. These structures are descriptively called debris disks. The first part of my thesis analyzes the Spitzer Space Telescope Observations of δ Velorum. The 24 μm Spitzer images revealed a bow shock structure in front of the star. My analysis showed that this is a result of the star’s high speed interaction with the surrounding interstellar medium. We place this observation and model in context of debris disk detections and the origin of λ Boötis stars. The second part of my thesis summarizes our observational results on the open cluster Praesepe. Using 24 μm data, I investigated the fraction of stars with mid-infrared excess, likely to have debris disks. I also assembled all results from previous debris disk studies and followed the evolution of the fraction of stars with debris disks. The majority of debris disks systems are evolved, few hundred million or a Gyr old. Since the dissipation timescale for the emitting dust particles is less than the age of these systems, they have to be constantly replenished through collisional grinding of the larger bodies. The last two chapters of my thesis is a theoretical analysis of the collisional cascade in debris disks, the process that produces the constant level of dust particles detected. I introduce a numerical model that takes into account all types of destructive collisions in the systems and solves the full scattering equation. I show results of comparisons between my and other published models and extensive verification tests of my model. I also analyze the evolution of the particle size distribution as a function of the variables in my model and show that the model itself is quite robust against most variations.

In this thesis, we examine the role of accretion disks in the formation of stellar systems, focusing on young massive disks which regulate the flow of material from the parent molecular core down to the star. We study the evolution of disks with high infall rates that develop strong gravitational instabilities. We begin in chapter 1 with a review of the observations and theory which underpin models for the earliest phases of star formation and provide a brief review of basic accretion disk physics, and the numerical methods which we employ. In chapter 2 we outline the current models of binary and multiple star formation, and review their successes and shortcomings from a theoretical and observational perspective. In chapter 3 we begin with a relatively simple analytic model for disks around young, very massive stars, showing that instability in these disks may be responsible for the higher multiplicity fraction of massive stars, and perhaps the upper mass to which they grow. We extend these models in chapter 4 to explore the properties of disks and the formation of binary companions across a broad range of stellar masses. In particular, we model the role of global and local mechanisms for angular momentum transport in regulating the relative masses of disks and stars. We follow the evolution of these disks throughout the main accretion phase of the system, and predict the trajectory of disks through parameter space. We follow up on the predictions made in our analytic models with a series of high resolution, global numerical experiments in chapter 5. Here we propose and test a new parameterization for describing rapidly accreting, gravitationally unstable disks. We find that disk properties and system multiplicity can be mapped out well in this parameter space. Finally, in chapter 6, we address whether our studies of unstable disks are relevant to recently detected massive planets on wide orbits around their central stars.

Stars of all ages and evolutionary stages are seen to be surrounded by discs of matieral. during the formation of a stellar system the stars are orbited by a massive protoplanetary disc composed of interstellar gas and dust, in which planet formation occurs. Betewwen 1 and 10 Myr the protoplanetary disc disperses, leaving behind the newly formed system of planets and smaller bodies. The remaining material which has not formed into planets is referred to as a debris disc. Even though the interstellar dust grains from the protoplanetary disc have long been removed from the system, debris discs can contain large quantities of dust due to collisions between larger bodies and cometary activity. such dust can be detected by its thermal emission. This thesis focuses on observational studies at far-infrared and sub-millimetre wavelengths of debris discs and the late stages of protoplanetary disc evolution. An overview of surveys for debris discs performed to date is presented, highlighting the limitations and statistical biases. the motivation, design and sample selection for two large surveys for debris discs around nearby stars, with the Hershel space observatory and the SCUBA-2 sub-millimetre camera on the James Clerk Maxwell Telescope, are described. The combination of a uniform obstevational strategy, longer wavelengths than previous surveys, and a large, clearly chosen sample - unbiased by stellar properties - will allow robust statistical conclusions of how the incidence and properties of debris discs depend on system parameters such as stellar mass, age, metallicity, binarity and the presence of planets. As a precursor to the Hershel and SCUBA-2 surveys, a volume-limited ample of 130 A type star systems was surveyed using observations at 24 and 70 μm, which were required to determine the presence of emission from dust, were predicted by fitting model flux distributions to optical and near-infrared photometry. Debris discs were detected around 46 systems, 12 of which including the system with the largest dust mass - are new discoveries. This survey adds to the results of previous studies which show that debris disc incidence is not correleated with host star metallicity despite the wll known giant planet - metallicity correlation, This is in accordance with what is predicted from the core accretion theory of planet formation. The most signigicant result from this survey is that, contrary to results reported in a previous work, debris discs are oberall less common around binary stars. Further investigation shows that systems with separations of ~3-150 AU are especially deficient of debris, while closer binaries and the primaries of wider binaries show debris detection rates consisten with those for single stars. A sample of circumstellar discs around 29 young stellar systems with ages of 5-30 Myr were observed with the LABOCA sub-millimetre instrument on the APEX telescope at 870μm, to provide disc masses or mass upper limits in support of a large Hershel programme. These targets included the η Chamaeleontis cluster and four bright Herbig Ae/Be stars which have not previously been observed at this wavelength. All but the Herbig Ae/Be stars were not detected, and 3σ dust mass upper limits of ~ 0.1-3 M are determined, with corresponding total disc masses of ~0.03-1Mjup. These mass limits indicate that there is insufficient remaining material in these discs to form gas giant planets, and add to the prevailing view that protoplanetary discs typically disperse within 10 Myr and that gas giant planet formation must be completed before this time. A search for cold dust emimission from two of the Solar System's nearest neighbours - α Centauri AB and ε Indi - was also performed with LABOCA. In both cases no debris disc emission was detected. A bright resolved feature was detected near α Centauri AB, nowever, follow-up observations at a second epoch, two years after the initial observations, showed that the feature is not co-moving with the stars. It is argued that the feature is most likely a pre-stellar core. The stars α Centauri A and B are detected, which is one of only very few detections of main sequence stellar photospheres at sub-millimetre wavelengths.

We present high signal-to-noise ratio, precise Y JH photometry and Y band (0.957-1.120 mu m) spectroscopy of HD 1160 B, a young substellar companion discovered from the Gemini NICI Planet Finding Campaign using the Subaru Coronagraphic Extreme Adaptive Optics instrument and the Gemini Planet Imager. HD 1160 B has typical mid-M dwarf-like infrared colors and a spectral type of M5.5(-0.5)(+1.0), where the blue edge of our Y band spectrum rules out earlier spectral types. Atmospheric modeling suggests HD 1160 B has an effective temperature of 3000-3100 K, a surface gravity of log g - 4-4.5, a radius of. 1.55 +/- 0.10 R-J, and a luminosity of log L/L circle dot - 2.76 +/- 0.05. Neither the primary's Hertzspring-Russell diagram position nor atmospheric modeling of HD 1160 B show evidence for a subsolar metallicity. Interpretation of the HD 1160 B spectroscopy depends on which stellar system components are used to estimate the age. Considering HD 1160 A, B and C jointly, we derive an age of 80-125 Myr, implying that HD 1160 B straddles the hydrogen-burning limit (70-90 M-J) If we consider HD 1160 A alone, younger ages (20-125 Myr) and a brown dwarf-like mass (35-90 M-J) are possible. Interferometric measurements of the primary, a precise Gaia parallax, and moderate-resolution spectroscopy can better constrain the system's age and how HD 1160 B fits within the context of (sub) stellar evolution.

Context. HR4796A is surrounded by a debris disc, observed in scattered light as an inclined ring with a high surface brightness. Past observations have raised several questions. First, a strong brightness asymmetry detected in polarised reflected light has recently challenged our understanding of scattering by the dust particles in this system. Secondly, the morphology of the ring strongly suggests the presence of planets, although no planets have been detected to date. Aims. We aim here at measuring with high accuracy the morphology and photometry of the ring in scattered light, in order to derive the phase function of the dust and constrain its near-infrared spectral properties. We also want to constrain the presence of planets and set improved constraints on the origin of the observed ring morphology. Methods. We obtained high-angular resolution coronagraphic images of the circumstellar environment around HR4796A with VLT/SPHERE during the commissioning of the instrument in May 2014 and during guaranteed-time observations in February 2015. The observations reveal for the first time the entire ring of dust, including the semi-minor axis that was previously hidden either behind the coronagraphic spot or in the speckle noise. Results. We determine empirically the scattering phase function of the dust in the H band from 13.6 degrees to 166.6 degrees. It shows a prominent peak of forward scattering, never detected before, for scattering angles below 30 degrees. We analyse the reflectance spectra of the disc from the 0.95 mu m to 1.6 mu m, confirming the red colour of the dust, and derive detection limits on the presence of planetary mass objects. Conclusions. We confirm which side of the disc is inclined towards the Earth. The analysis of the phase function, especially below 45 degrees, suggests that the dust population is dominated by particles much larger than the observation wavelength, of about 20 mu m. Compact Mie grains of this size are incompatible with the spectral energy distribution of the disc, however the observed rise in scattering efficiency beyond 50 degrees points towards aggregates which could reconcile both observables. We do not detect companions orbiting the star, but our high-contrast observations provide the most stringent constraints yet on the presence of planets responsible for the morphology of the dust.

The NASA K2 mission uses photometry to find planets transiting stars of various types. M dwarfs are of high interest since they host more short-period planets than any other type of main-sequence star and transiting planets around M dwarfs have deeper transits compared to other main-sequence stars. In this paper, we present stellar parameters from K and M dwarfs hosting transiting planet candidates discovered by our team. Using the SOFI spectrograph on the European Southern Observatory's New Technology Telescope, we obtained R approximate to 1000 J-, H-, and K-band (0.95-2.52 mu m) spectra of 34 late-type K2 planet and candidate planet host systems and 12 bright K4-M5 dwarfs with interferometrically measured radii and effective temperatures. Out of our 34 late-type K2 targets, we identify 27 of these stars as M dwarfs. We measure equivalent widths of spectral features, derive calibration relations using stars with interferometric measurements, and estimate stellar radii, effective temperatures, masses, and luminosities for the K2 planet hosts. Our calibrations provide radii and temperatures with median uncertainties of 0.059 R-circle dot (16.09%) and 160 K (4.33%), respectively. We then reassess the radii and equilibrium temperatures of known and candidate planets based on our spectroscopically derived stellar parameters. Since a planet's radius and equilibrium temperature depend on the parameters of its host star, our study provides more precise planetary parameters for planets and candidates orbiting late-type stars observed with K2. We find a median planet radius and an equilibrium temperature of approximately 3 R-circle plus and 500 K, respectively, with several systems (K2-18b and K2-72e) receiving near-Earth-like levels of incident irradiation.

The disk mass is among the most important input parameter for every planet formation model to determine the number and masses of the planets that can form. We present an ALMA 887 mu m survey of the disk population around objects from similar to 2 to 0.03 M-circle dot in the nearby similar to 2 Myr old Chamaeleon I star-forming region. We detect thermal dust emission from 66 out of 93 disks, spatially resolve 34 of them, and identify two disks with large dust cavities of about 45 au in radius. Assuming isothermal and optically thin emission, we convert the 887 mu m flux densities into dust disk masses, hereafter M-dust. We find that the M-dust-M* relation is steeper than linear and of the form M-dust proportional to (M*)(1.3-1.9), where the range in the power-law index reflects two extremes of the possible relation between the average dust temperature and stellar luminosity. By reanalyzing all millimeter data available for nearby regions in a self-consistent way, we show that the 1-3 Myr old regions of Taurus, Lupus, and Chamaeleon. I share the same M-dust-M* relation, while the 10 Myr old Upper. Sco association has a steeper relation. Theoretical models of grain growth, drift, and fragmentation reproduce this trend and suggest that disks are in the fragmentation-limited regime. In this regime millimeter grains will be located closer in around lower-mass stars, a prediction that can be tested with deeper and higher spatial resolution ALMA observations.

We present new astrometric measurements from our ongoing monitoring campaign of the HR 8799 directly imaged planetary system. These new data points were obtained with NIRC2 on the W.M. Keck II 10 m telescope between 2009 and 2014. In addition, we present updated astrometry from previously published observations in 2007 and 2008. All data were reduced using the SOSIE algorithm, which accounts for systematic biases present in previously published observations. This allows us to construct a self-consistent data set derived entirely from NIRC2 data alone. From this data set, we detect acceleration for two of the planets (HR 8799b and e) at >3 sigma. We also assess possible orbital parameters for each of the four planets independently. We find no statistically significant difference in the allowed inclinations of the planets. Fitting the astrometry while forcing coplanarity also returns chi(2) consistent to within 1 sigma of the best fit values, suggesting that if inclination offsets of less than or similar to 20 degrees are present, they are not detectable with current data. Our orbital fits also favor low eccentricities, consistent with predictions from dynamical modeling. We also find period distributions consistent to within 1 sigma with a 1:2:4:8 resonance between all planets. This analysis demonstrates the importance of minimizing astrometric systematics when fitting for solutions to highly undersampled orbits.

We analyze full-orbit phase curve observations of the transiting hot Jupiters WASP-19b and HAT-P-7b at 3.6 and 4.5 mu m, obtained using the Spitzer Space Telescope. For WASP-19b, we measure secondary eclipse depths of 0.485% +/- 0.024% and 0.584% +/- 0.029% at 3.6 and 4.5 mu m, which are consistent with a single blackbody with effective temperature 2372 +/- 60 K. The measured 3.6 and 4.5 mu m secondary eclipse depths for HAT-P-7b are 0.156% +/- 0.009% and 0.190% +/- 0.006%, which are well described by a single blackbody with effective temperature 2667 +/- 57 K. Comparing the phase curves to the predictions of one-dimensional and three-dimensional atmospheric models, we find that WASP-19b's dayside emission is consistent with a model atmosphere with no dayside thermal inversion and moderately efficient day-night circulation. We also detect an eastward-shifted hotspot, which suggests the presence of a superrotating equatorial jet. In contrast, HAT-P-7b's dayside emission suggests a dayside thermal inversion and relatively inefficient day-night circulation; no hotspot shift is detected. For both planets, these same models do not agree with the measured nightside emission. The discrepancies in the model-data comparisons for WASP-19b might be explained by high-altitude silicate clouds on the nightside and/or high atmospheric metallicity, while the very low 3.6 mu m nightside planetary brightness for HAT-P-7b may be indicative of an enhanced global C/O ratio. We compute Bond albedos of 0.38 +/- 0.06 and 0 (<0.08 at 1 sigma) for WASP-19b and HAT-P-7b, respectively. In the context of other planets with thermal phase curve measurements, we show that WASP-19b and HAT-P-7b fit the general trend of decreasing day-night heat recirculation with increasing irradiation.