Massive stars like company. However, low-mass companions have remained extremely difficult to detect at angular separations ($\rho$) smaller than 1" (~1000-3000 au considering typical distance to nearby massive stars) given the large brightness contrast between the companion and the central star. Constraints on the low-mass end of the companions mass-function for massive stars are however needed, for example to help distinguishing between various scenarios for the formation of massive stars. To obtain statistically significant constraint on the presence of low-mass companions beyond the typical detection limit of current surveys ($\Delta \mathrm{mag} \lesssim 5$ at $\rho \lesssim 1$"), we initiated a survey of O and Wolf-Rayet stars in the Carina region using the SPHERE coronagraphic instrument on the VLT. In this first paper, we aim to introduce the survey, to present the methodology and to demonstrate the capability of SPHERE for massive stars using the multiple system QZ~Car. High-contrast imaging techniques, such as angular- and spectral-differential imaging techniques as well as PSF-fitting, were applied to detect and measure the relative flux of companions in each spectral channel of the instrument. We detected 19 sources around the QZ~Car system with detection limits of 9~mag at $\rho > 200$~mas for IFS and as faint as 13~mag at $\rho$ > 1".8 for IRDIS (corresponding to sub-solar masses for potential companions). All but two are reported here for the first time. Based on this proof of concept, we showed that VLT/SPHERE allows us to reach the sub-solar mass regime of the companion mass function. This paves the way for this type of observation with a large sample of massive stars to provide novel constraints on the multiplicity of massive stars in a region of the parameter space that has remained inaccessible so far.

The streaming instability (SI) has been extensively studied in the linear and non-linear regimes as a mechanism to concentrate solids and trigger planetesimal formation in the midplane of protoplanetary discs. A related dust settling instability (DSI) applies to particles while settling towards the midplane. The DSI has previously been studied in the linear regime, with predictions that it could trigger particle clumping away from the midplane. This work presents a range of linear calculations and non-linear simulations, performed with FARGO3D, to assess conditions for DSI growth. We expand on previous linear analyses by including particle size distributions and performing a detailed study of the amount of background turbulence needed to stabilize the DSI. When including binned size distributions, the DSI often produces converged growth rates with fewer bins than the standard SI. With background turbulence, we find that the most favorable conditions for DSI growth are weak turbulence, characterized by $\alpha \lesssim 10^{-6}$ with intermediate-sized grains that settle from one gas scale-height. These conditions could arise during a sudden decrease in disc turbulence following an accretion outburst. Ignoring background turbulence, we performed a parameter survey of local 2D DSI simulations. Particle clumping was either weak or occurred slower than particles settle. Clumping was reduced by a factor of two in a comparison 3D simulation. Overall, our results strongly disfavor the hypothesis that the DSI significantly promotes planetesimal formation. Non-linear simulations of the DSI with different numerical methods could support or challenge these findings.

We investigate the time evolution of dense cores identified in molecular cloud simulations using dendrograms, which are a common tool to identify hierarchical structure in simulations and observations of star formation. We develop an algorithm to link dendrogram structures through time using the three-dimensional density field from magnetohydrodynamical simulations, thus creating histories for all dense cores in the domain. We find that the population-wide distributions of core properties are relatively invariant in time, and quantities like the core mass function match with observations. Despite this consistency, an individual core may undergo large (>40%), stochastic variations due to the redefinition of the dendrogram structure between timesteps. This variation occurs independent of environment and stellar content. We identify a population of short-lived (<200 kyr) overdensities masquerading as dense cores that may comprise ~20% of any time snapshot. Finally, we note the importance of considering the full history of cores when interpreting the origin of the initial mass function; we find that, especially for systems containing multiple stars, the core mass defined by a dendrogram leaf in a snapshot is typically less than the final system stellar mass. This work reinforces that there is no time-stable density contour that defines a star-forming core. The dendrogram itself can induce significant structure variation between timesteps due to small changes in the density field. Thus, one must use caution when comparing dendrograms of regions with different ages or environment properties because differences in dendrogram structure may not come solely from the physical evolution of dense cores.

Blue straggler stars (BSS) have been identified in star clusters and in field populations in our own Milky Way galaxy and in its satellite galaxies. They manifest as stars bluer and more luminous than the dominant old population, and usually have a spatial distribution that follows the old population. Their progenitors are likely to have been close binaries. We investigate trends of the BSS population in dwarf spheroidal galaxies (dSph) and in the bulge of the Milky Way and find an anti-correlation between the relative frequency of BSS and the metallicity of the parent population. The rate of occurrence of BSS in the metal-poor dwarf galaxies is approximately twice that found in the solar-metallicity bulge population. This trend of decreasing relative population of BSS with increasing metallicity mirrors that found for the close-binary fraction in the field population of the Milky Way. We argue that the dominant mode of BSS formation in low-density environments is likely to be mass transfer in close-binary systems. It then follows that the similarity between the trends for BSS in the dSph and field stars in our Galaxy supports the proposal that the small-scale fragmentation during star formation is driven by the same dominant physical process, despite the diversity in environments, plausibly gravitational instability of proto-stellar discs.

We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase.

The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this white paper, we argue that in order to reach the promised critical breakthroughs in astrophysics over the next decade and well beyond, the national agencies must take a serious approach to investment in theoretical astrophysics research. We discuss the role of theory in shaping our understanding of the universe, and then we provide a multi-level strategy, from the grassroots to the national, to address the current underinvestment in theory relative to observational work.

We carry out 2D viscous hydrodynamics simulations of circumbinary disk (CBD) accretion using {\footnotesize AREPO}. We resolve the accretion flow from a large-scale CBD down to the streamers and disks around individual binary components. Extending our recent studies \citep{mun19}, we consider circular binaries with various mass ratios ($0.1\leq q_{\rm{b}}\leq1$) and study accretion from ``infinite'', steady-supply disks and from finite-sized, viscously spreading tori. For ``infinite'' disks, a global steady state can be reached, and the accretion variability has a dominant frequency ${\sim}0.2\Omega_{\rm{b}}$ for $q_{\rm{b}}>0.5$ and $\Omega_{\rm{b}}$ for $q_{\rm{b}}<0.5$, ($\Omega_{\rm{b}}$ is the binary angular frequency). We find that the accretion ``eigenvalue'' $l_0$ -- the net angular momentum transfer from the disk to the binary per unit accreted mass -- is always positive and falls in the range ($0.65$-$0.85)a_{\rm b}^2\Omega_{\rm{b}}$ (with $a_{\rm{b}}$ the binary separation), depending weakly on the mass ratio and viscosity. This leads to binary expansion when $q_{\rm{b}}\gtrsim0.3$. Accretion from a finite torus can be separated into two phases: an initial transient phase, corresponding to the filling of the binary cavity, followed by a viscous pseudo-stationary phase, during which the torus viscously spreads and accretes onto the binary. In the viscous phase, the net torque on the binary per unit accreted mass is close to $l_0$, the value derived for ``infinite'' disks. We conclude that similar-mass binaries accreting from CBDs gain angular momentum and expand over long time scales. This result significantly impacts the coalescence of supermassive binary black holes and newly formed binary stars. We offer a word of caution against conclusions drawn from simulations of transient accretion onto empty circumbinary cavities.

Transitional discs have central regions characterised by significant depletion of both dust and gas compared to younger, optically-thick discs. However, gas and dust are not depleted by equal amounts: gas surface densities are typically reduced by factors of $\sim 100$, but small dust grains are sometimes depleted by far larger factors, to the point of being undetectable. While this extreme dust depletion is often attributed to planet formation, in this paper we show that another physical mechanism is possible: expulsion of grains from the disc by radiation pressure. We explore this mechanism using 2D simulations of dust dynamics, simultaneously solving the equation of radiative transfer with the evolution equations for dust diffusion and advection under the combined effects of stellar radiation and hydrodynamic interaction with a turbulent, accreting background gas disc. We show that, in transition discs that are depleted in both gas and dust fraction by factors of $\sim 100-1000$ compared to minimum mass Solar nebular values, radiative clearing of any remaining $\sim 0.5$ $\mu$m and larger grains is both rapid and inevitable. The process is size-dependent, with smaller grains removed fastest and larger ones persisting for longer times. Our proposed mechanism thus naturally explains the extreme depletion of small grains commonly-found in transition discs. We further suggest that the dependence of this mechanism on grain size and optical properties may explain some of the unusual grain properties recently discovered in a number of transition discs. The simulation code we develop is freely available.

Close binaries suppress the formation of circumstellar (S-type) planets and therefore significantly bias the inferred planet occurrence rates and statistical trends. After compiling various radial velocity and high-resolution imaging surveys, we determine that binaries with a < 1 au fully suppress S-type planets, binaries with a = 10 au host close planets at 15% the occurrence rate of single stars, and wide binaries with a > 200 au have a negligible effect on planet formation. We show that F = 43% +/- 6% of solar-type primaries do not host close planets due to suppression by close stellar companions. By removing spectroscopic binaries from their samples, radial velocity surveys for giant planets boost their detection rates by a factor of 1/(1-F) = 1.8 +/- 0.2 compared to transiting surveys. This selection bias fully accounts for the discrepancy in hot Jupiter occurrence rates inferred from these two detection methods. Correcting for both planet suppression by close binaries and transit dilution by wide binaries, the occurrence rate of small planets orbiting single G-dwarfs is 2.1 +/- 0.3 times larger than the rate inferred from all G-dwarfs in the Kepler survey. About half (but not all) of the observed increase in small, short-period planets toward low-mass hosts can be explained by the corresponding decrease in the binary fraction. Finally, we demonstrate that the apparent enhancement of wide stellar companions to hot Jupiter hosts is due to multiple selection effects. Although very close binaries with secondary masses M_2 > 10 M_J preferentially have wide tertiary companions, genuine hot Jupiters with M_p = 0.2-4 M_J that formed via core accretion do not exhibit a statistically significant excess of wide stellar companions.

Stars rarely form in isolation. Nearly half of the stars in the Milky Way have a companion, and this fraction increases in star-forming regions. However, why some dense cores and filaments form bound pairs while others form single stars remains unclear. We present a set of three-dimensional, gravo-magnetohydrodynamic simulations of turbulent star-forming clouds, aimed at understanding the formation and evolution of multiple-star systems formed through large scale (>~$10^3$ AU) turbulent fragmentation. We investigate three global magnetic field strengths, with global mass-to-flux ratios of $\mu_\phi$=2, 8, and 32. The initial separations of protostars in multiples depends on the global magnetic field strength, with stronger magnetic fields (e.g., $\mu_\phi$=2) suppressing fragmentation on smaller scales. The overall multiplicity fraction (MF) is between 0.4-0.6 for our strong and intermediate magnetic field strengths, which is in agreement with observations. The weak field case has a lower fraction. The MF is relatively constant throughout the simulations, even though stellar densities increase as collapse continues. While the MF rarely exceeds 60% in all three simulations, over 80% of all protostars are part of a binary system at some point. We additionally find that the distribution of binary spin mis-alignment angles is consistent with a randomized distribution. In all three simulations, several binaries originate with wide separations and dynamically evolve to <~ $10^2$ AU separations. We show that a simple model of mass accretion and dynamical friction with the gas can explain this orbital evolution.

We present ALMA (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observations at similar wavelengths. We identify six compact continuum sources within OMC2-FIR4, four in OMC2-FIR3, and one additional source just outside OMC2-FIR4. This continuum emission is tracing the inner envelope and/or disk emission on less than 100~AU scales. HOPS-108 is the only protostar in OMC2-FIR4 that exhibits emission from high-excitation transitions of complex organic molecules (e.g., methanol and other lines) coincident with the continuum emission. HOPS-370 in OMC2-FIR3 with L~$\sim$~360~\lsun, also exhibits emission from high-excitation methanol and other lines. The methanol emission toward these two protostars is indicative of temperatures high enough to thermally evaporate methanol from icy dust grains; overall these protostars have characteristics similar to hot corinos. We do not identify a clear outflow from HOPS-108 in \twco, but find evidence of interaction between the outflow/jet from HOPS-370 and the OMC2-FIR4 region. The multitude of observational constraints indicate that HOPS-108 is likely a low to intermediate-mass protostar in its main mass accretion phase and it is the most luminous protostar in OMC2-FIR4. The high resolution data presented here are essential for disentangling the embedded protostars from their surrounding dusty environments and characterizing them.

A major outstanding question regarding the formation of planetary systems is whether wide-orbit giant planets form differently than close-in giant planets. We aim to establish constraints on two key parameters that are relevant for understanding the formation of wide-orbit planets: 1) the relative mass function and 2) the fraction of systems hosting multiple companions. In this study, we focus on systems with directly imaged substellar companions, and the detection limits on lower-mass bodies within these systems. First, we uniformly derive the mass probability distributions of known companions. We then combine the information contained within the detections and detection limits into a survival analysis statistical framework to estimate the underlying mass function of the parent distribution. Finally, we calculate the probability that each system may host multiple substellar companions. We find that 1) the companion mass distribution is rising steeply toward smaller masses, with a functional form of $N\propto M^{-1.3\pm0.3}$, and consequently, 2) many of these systems likely host additional undetected sub-stellar companions. Combined, these results strongly support the notion that wide-orbit giant planets are formed predominantly via core accretion, similar to the better studied close-in giant planets. Finally, given the steep rise in the relative mass function with decreasing mass, these results suggest that future deep observations should unveil a greater number of directly imaged planets.

The majority of the discovered transiting circumbinary planets are located very near the innermost stable orbits permitted, raising questions about the origins of planets in such perturbed environments. Most favored formation scenarios invoke formation at larger distances and subsequent migration to their current locations. Disk-driven planet migration in multi-planet systems is likely to trap planets in mean motion resonances and drive planets inward into regions of larger dynamical perturbations from the binary. We demonstrate how planet-planet resonances can interact with the binary through secular forcing and mean-motion resonances, driving chaos in the system. We show how this chaos will shape the architecture of circumbinary systems, with specific applications to Kepler 47 and the Pluto-Charon system, limiting maximum possible stable eccentricities and indicating what resonances are likely to exist. We are also able to constrain the minimum migration rates of resonant circumbinary planets.

Significant advances have been made over the past decade in the characterization of multiple protostar systems, enabled by the Karl G. Jansky Very Large Array (VLA), high-resolution infrared observations with the Hubble Space Telescope, and ground-based facilities. To further understand the mechanism(s) of multiple star formation, a combination of statistics, high-angular resolution radio/millimeter continuum imaging, characterization of kinematic structure, magnetic fields via polarimetry, and comparison with numerical simulations are needed. Thus, understanding the origin of stellar multiplicity in different regimes of companion separation will soon be within reach. However, to overcome challenges that studies in this field are now confronted with, a range of new capabilities are required: a new millimeter/centimeter wave facility with 10 mas resolution at {\lambda}=1 cm, space-based near to far-infrared observatories, continued development of low to high resolution spectroscopy on 3m to 10m class telescopes, and an ELT-class telescope with near to mid-infrared imaging/spectroscopic capability.

Multiplicity is a fundamental property that is set early during stellar lifetimes, and it is a stringent probe of the physics of star formation. The distribution of close companions around young stars is still poorly constrained by observations. We present an analysis of stellar multiplicity derived from APOGEE-2 spectra obtained in targeted observations of nearby star-forming regions. This is the largest homogeneously observed sample of high-resolution spectra of young stars. We developed an autonomous method to identify double lined spectroscopic binaries (SB2s). Out of 5007 sources spanning the mass range of $\sim$0.05--1.5 \msun, we find 399 binaries, including both RV variables and SB2s. The mass ratio distribution of SB2s is consistent with a uniform for $q<0.95$ with an excess of twins with $q>0.95$. The period distribution is consistent with what has been observed in close binaries ($<10$ AU) in the evolved populations. Three systems are found to have $q\sim$0.1, with a companion located within the brown dwarf desert. There are not any strong trends in the multiplicity fraction (MF) as a function of cluster age from 1 to 100 Myr. There is a weak dependence on stellar density, with companions being most numerous at $\Sigma_*\sim30$ stars/pc$^{-2}$, and decreasing in more diffuse regions. Finally, disk-bearing sources are deficient in SB2s (but not RV variables) by a factor of $\sim$2; this deficit is recovered by the systems without disks. This may indicate a quick dispersal of disk material in short-period equal mass systems that is less effective in binaries with lower $q$.

Almost all massive stars have bound stellar companions, existing in binaries or higher-order multiples. While binarity is theorized to be an essential feature of how massive stars form, essentially all information about such properties is derived from observations of already formed stars, whose orbital properties may have evolved since birth. Little is known about binarity during formation stages. Here we report high angular resolution observations of 1.3 mm continuum and H30alpha recombination line emission, which reveal a massive protobinary with apparent separation of 180 au at the center of the massive star-forming region IRAS07299-1651. From the line-of-sight velocity difference of 9.5 km/s of the two protostars, the binary is estimated to have a minimum total mass of 18 solar masses, consistent with several other metrics, and maximum period of 570 years, assuming a circular orbit. The H30alpha line from the primary protostar shows kinematics consistent with rotation along a ring of radius of 12 au. The observations indicate that disk fragmentation at several hundred au may have formed the binary, and much smaller disks are feeding the individual protostars.

We report molecular line observations of the massive protostellar source G339.88-1.26 with the Atacama Large Millimeter/Submillimeter Array. The observations reveal a highly collimated SiO jet extending from the 1.3 mm continuum source, which connects to a slightly wider but still highly collimated CO outflow. Rotational features perpendicular to the outflow axis are detected in various molecular emissions, including SiO, SO2, H2S, CH3OH, and H2CO emissions. Based on their spatial distributions and kinematics, we find that they trace different parts of the envelope-disk system. The SiO emission traces the disk and inner envelope in addition to the jet. The CH3OH and H2CO emissions mostly trace the infalling-rotating envelope, and are enhanced around the transition region between envelope and disk, i.e., the centrifugal barrier. The SO2 and H2S emissions are enhanced around the centrifugal barrier, and also trace the outer part of the disk. Envelope kinematics are consistent with rotating-infalling motion, while those of the disk are consistent with Keplerian rotation. The radius and velocity of the centrifugal barrier are estimated to be about 530 au and 6 km/s, leading to a central mass of about 11 solar masses, consistent with estimates based on spectral energy distribution fitting. These results indicate that an ordered transition from an infalling-rotating envelope to a Keplerian disk through a centrifugal barrier, accompanied by changes of types of molecular line emissions, is a valid description of this massive protostellar source. This implies that at least some massive stars form in a similar way as low-mass stars via Core Accretion.

We use a cross-calibration of Hipparcos and Gaia DR2 astrometry for $\beta$ Pic to measure the mass of the giant planet $\beta$ Pic b $(13\pm3$ $M_{\rm Jup})$ in a comprehensive joint orbit analysis that includes published relative astrometry and radial velocities. Our mass uncertainty is somewhat higher than previous work because our astrometry from the Hipparcos-Gaia Catalog of Accelerations accounts for the error inflation and systematic terms that are required to bring the two data sets onto a common astrometric reference frame, and because we fit freely for the host-star mass $(1.84\pm0.05$ $M_{\odot})$. This first model-independent mass for a directly imaged planet is inconsistent with cold-start models given the age of the $\beta$ Pic moving group $(22\pm6$ Myr) but consistent with hot- and warm-start models, concordant with past work. We find a higher eccentricity $(0.24\pm0.06)$ for $\beta$ Pic b compared to previous orbital fits. If confirmed by future observations, this eccentricity may help explain inner edge, scale height, and brightness asymmetry of $\beta$ Pic's disk. It could also potentially signal that $\beta$ Pic b has migrated inward to its current location, acquiring its eccentricity from interaction with the 3:1 outer Lindblad resonance in the disk.

We present the results of the largest $L^{\prime}$ ($3.8~\mu$m) direct imaging survey for exoplanets to date, the Large Binocular Telescope Interferometer (LBTI) Exozodi Exoplanet Common Hunt (LEECH). We observed 98 stars with spectral types from B to M. Cool planets emit a larger share of their flux in $L^{\prime}$ compared to shorter wavelengths, affording LEECH an advantage in detecting low-mass, old, and cold-start giant planets. We emphasize proximity over youth in our target selection, probing physical separations smaller than other direct imaging surveys. For FGK stars, LEECH outperforms many previous studies, placing tighter constraints on the hot-start planet occurrence frequency interior to $\sim20$ au. For less luminous, cold-start planets, LEECH provides the best constraints on giant-planet frequency interior to $\sim20$ au around FGK stars. Direct imaging survey results depend sensitively on both the choice of evolutionary model (e.g., hot- or cold-start) and assumptions (explicit or implicit) about the shape of the underlying planet distribution, in particular its radial extent. Artificially low limits on the planet occurrence frequency can be derived when the shape of the planet distribution is assumed to extend to very large separations, well beyond typical protoplanetary dust-disk radii ($\lesssim50$ au), and when hot-start models are used exclusively. We place a conservative upper limit on the planet occurrence frequency using cold-start models and planetary population distributions that do not extend beyond typical protoplanetary dust-disk radii. We find that $\lesssim90\%$ of FGK systems can host a 7 to 10 $M_{\mathrm{Jup}}$ planet from 5 to 50 au. This limit leaves open the possibility that planets in this range are common.

Bathed in intense ionizing radiation, close-in gaseous planets undergo hydrodynamic atmospheric escape, which ejects the upper extent of their atmospheres into the interplanetary medium. Ultraviolet detections of escaping gas around transiting planets corroborate such a framework. Exposed to the stellar environment, the outflow is shaped by its interaction with the stellar wind and by the planet's orbit. We model these effects using Athena to perform 3-D radiative-hydrodynamic simulations of tidally-locked hydrogen atmospheres receiving large amounts of ionizing extreme-ultraviolet flux in various stellar environments for the low-magnetic-field case. Through a step-by-step exploration of orbital and stellar wind effects on the planetary outflow, we find three structurally distinct stellar wind regimes: weak, intermediate, and strong. We perform synthetic Lyman-$\alpha$ observations and find unique observational signatures for each regime. A weak stellar wind$\textrm{---}$which cannot confine the planetary outflow, leading to a torus of material around the star$\textrm{---}$has a pre-transit, red-shifted dayside arm and a slightly redward-skewed spectrum during transit. The intermediate regime truncates the dayside outflow at large distances from the planet and causes periodic disruptions of the outflow, producing observational signatures that mimic a double transit. The first of these dips is blue-shifted and precedes the optical transit. Finally, strong stellar winds completely confine the outflow into a cometary tail and accelerate the outflow outwards, producing large blue-shifted signals post-transit. Across all three regimes, large signals occur far outside of transit, offering motivation to continue ultraviolet observations outside of direct transit.

Gap-opening planets can generate dust-trapping vortices that may explain some of the latest discoveries of high-contrast crescent-shaped dust asymmetries in transition discs. While planet-induced vortices were previously thought to have concentrated shapes, recent computational work has shown that these features naturally become much more elongated in the gas when simulations account for the relatively long timescale over which planets accrete their mass. In this work, we conduct two-fluid hydrodynamical simulations of vortices induced by slowly-growing Jupiter-mass planets in discs with very low viscosity ($\alpha = 3 \times 10^{-5}$). We simulate the dust dynamics for four particle sizes spanning 0.3 mm to 1 cm in order to produce synthetic ALMA images. In our simulations, we find that an elongated vortex still traps dust, but not directly at its center. With a flatter pressure bump and disruptions from the planet's overlapping spiral density waves, the dust instead circulates around the vortex. This motion (1) typically carries the peak off-center, (2) spreads the dust out over a wider azimuthal extent $\geq 180^{\circ}$, (3) skews the azimuthal profile towards the front of the vortex, and (4) can also create double peaks in newly-formed vortices. In particular, we expect that the most defining observational signature, a peak offset of more than $30^{\circ}$, should be detectable $>30\%$ of the time in observations with a beam diameter of at most the planet's separation from its star.

While many observed debris disks are thought to have gaps suggestive of the presence of planets, direct imaging surveys do not find many high mass planets in these systems. We investigate if divergent migration is a viable mechanism for forming gaps in young debris disks with planets of low enough mass to currently elude detection. We perform numerical integrations of planet pairs embedded in planetesimal disks to assess the conditions for which divergent, planetesimal-driven migration occurs and gaps form within the disk. Gap widths and the migration rate of planets within a pair depend on both disk mass and the degree to which the planets share disk material. We find that planet pairs with planets more massive than Neptune can produce gaps with widths similar to their orbit distance within 10 Myr at orbit separations probed by direct imaging campaigns. Pairs of migrating super-Earths likely cannot form observable gaps on the same time and distance scales, however. Inferring the responsible planet masses from these gaps while neglecting migration could overestimate the mass of planets by more than an order of magnitude.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of multiple protostar systems in the Perseus molecular cloud previously detected by the Karl G. Jansky Very Large Array (VLA). We observed 17 close ($<$600~AU separation) multiple systems at 1.3~mm in continuum and five molecular lines (i.e., \twco, \cateo, \thco, H$_2$CO, SO) to characterize the circum-multiple environments in which these systems are forming. We detect at least one component in the continuum for the 17 multiple systems. In three systems, one companion is not detected, and for two systems the companions are unresolved at our observed resolution. We also detect circum-multiple dust emission toward 8 out of 9 Class 0 multiples. Circum-multiple dust emission is not detected toward any of the 8 Class I multiples. Twelve systems are detected in the dense gas tracers toward their disks/inner envelopes. For these 12 systems, we use the dense gas observations to characterize their formation mechanism. The velocity gradients in the circum-multiple gas are clearly orthogonal to the outflow directions in 8 out of the 12 systems, consistent with disk fragmentation. Moreover, only two systems with separations $<$200~AU are \textit{inconsistent} with disk fragmentation, in addition to the two widest systems ($>$500~AU). Our results suggest that disk fragmentation via gravitational instability is an important formation mechanism for close multiple systems, but further statistics are needed to better determine the relative fraction formed via this method.

We present the full disk-fit results VANDAM survey of all Class 0 and I protostars in the Perseus molecular cloud. We have 18 new protostellar disk candidates around Class 0 and I sources, which are well described by a simple, parametrized disk model fit to the 8 mm VLA dust-continuum observations. 33% of Class 0 protostars and just 11% of Class I protostars have candidate disks, while 78% of Class 0 and I protostars do not have signs of disks within our 12 AU disk diameter resolution limit, indicating that at 8 mm most disks in the Class 0 and I phases are <10 AU in radius. These small radii may be a result of surface brightness sensitivity limits. Modeled 8 mm radii are similar to the radii of known Class 0 disks with detected Keplerian rotation. Since our 8 mm data trace a population of larger dust grains which radially drift towards the protostar and are lower limits on true disk sizes, large disks at early times do not seem to be particularly rare. We find statistical evidence that Class 0 and I disks are likely drawn from the same distribution, meaning disk properties may be defined early in the Class 0 phase and do not undergo large changes through the Class I phase. By combining our candidate disk properties with previous polarization observations, we find a qualitative indication that misalignment between inferred envelope-scale magnetic fields and outflows may indicate disks on smaller scales in Class 0 sources.

There is now strong evidence that the close binary fraction (P < 10$^4$ days; a < 10 AU) of solar-type stars ($M_1$ = 0.6-1.5M$_{\odot}$) decreases significantly with metallicity. Although early surveys showed that the observed spectroscopic binary (SB) fractions in the galactic disk and halo are similar (e.g., Carney-Latham sample), these studies did not correct for incompleteness. In this study, we examine five different surveys and thoroughly account for their underlying selection biases to measure the intrinsic occurrence rate of close solar-type binaries. We re-analyze: (1) a volume-limited sample of solar-type stars, (2) an SB survey of high-proper-motion stars, (3) various SB samples of metal-poor giants, (4) the APOGEE survey of radial velocity (RV) variables, and (5) Kepler eclipsing binaries (EBs). The observed APOGEE RV variability fraction and Kepler EB fraction both decrease by a factor of $\approx$4 across $-$1.0 < [Fe/H] < 0.5 at the 22$\sigma$ and 9$\sigma$ confidence levels, respectively. After correcting for incompleteness, all five samples exhibit a quantitatively consistent anti-correlation between the intrinsic close binary fraction (a < 10 AU) and metallicity: $F_{\rm close}$ = 53%$\pm$12%, 40%$\pm$6%, 24%$\pm$4%, and 10%$\pm$3% at [Fe/H] = $-$3.0, $-$1.0, $-$0.2 (mean field metallicity), and +0.5, respectively. We present fragmentation models that explain why the close binary fraction of solar-type stars strongly decreases with metallicity while the wide binary fraction, close binary fraction of OB stars, and initial mass function are all constant across $-$1.5 < [Fe/H] < 0.5. The majority of solar-type stars with [Fe/H] < $-$1.0 will interact with a stellar companion, which has profound implications for binary evolution in old and metal-poor environments such as the galactic halo, bulge, thick disk, globular clusters, dwarf galaxies, and high-redshift universe.

We have discovered a wide planetary-mass companion to the $\beta$ Pic moving group member 2MASSJ02495639-0557352 (M6 VL-G) using CFHT/WIRCam astrometry from the Hawaii Infrared Parallax Program. In addition, Keck laser guide star adaptive optics aperture-masking interferometry shows that the host is itself a tight binary. Altogether, 2MASSJ0249-0557ABc is a bound triple system with an $11.6^{+1.0}_{-1.3}$ $M_{\rm Jup}$ object separated by $1950\pm200$ AU (40") from a relatively close ($2.17\pm0.22$ AU, 0.04") pair of $48^{+12}_{-13}$ $M_{\rm Jup}$ and $44^{+11}_{-14}$ $M_{\rm Jup}$ objects. 2MASSJ0249-0557AB is one of the few ultracool binaries to be discovered in a young moving group and the first confirmed in the $\beta$ Pic moving group ($22\pm6$ Myr). The mass, absolute magnitudes, and spectral type of 2MASSJ0249-0557 c (L2 VL-G) are remarkably similar to those of the planet $\beta$ Pic b (L2, $13.0^{+0.4}_{-0.3}$ $M_{\rm Jup}$). We also find that the free-floating object 2MASSJ2208+2921 (L3 VL-G) is another possible $\beta$ Pic moving group member with colors and absolute magnitudes similar to $\beta$ Pic b and 2MASSJ0249-0557 c. $\beta$ Pic b is the first directly imaged planet to have a "twin," namely an object of comparable properties in the same stellar association. Such directly imaged objects provide a unique opportunity to measure atmospheric composition, variability, and rotation across different pathways of assembling planetary-mass objects from the same natal material.

Emission from protostars at centimeter radio wavelengths has been shown to trace the free-free emission arising from ionizing shocks as a result of jets and outflows driven by protostars. Therefore, measuring properties of protostars at radio frequencies can provide valuable insights into the nature of their outflows and jets. We present a C-band (4.1 cm and 6.4 cm) survey of all known protostars (Class 0 and Class I) in Perseus as part of the VLA Nascent Disk and Multiplicity (VANDAM) Survey. We examine the known correlations between radio flux density and protostellar parameters such as bolometric luminosity and outflow force, for our sample. We also investigate the relationship between radio flux density and far-infrared line luminosities from Herschel. We show that free-free emission originates most likely from J-type shocks; however, the large scatter indicates that those two types of emission probe different time and spatial scales. Using C-band fluxes, we removed an estimation of free-free contamination from the corresponding Ka-band (9 mm) flux densities that primarily probe dust emission from embedded disks. We find that the compact ($<$~1") dust emission is lower for Class I sources (median dust mass 96 M$_{\oplus}$) relative to Class 0 (248 M$_{\oplus}$), but several times higher than in Class II (5-15 M$_{\oplus}$). If this compact dust emission is tracing primarily the embedded disk, as is likely for many sources, this result provides evidence for decreasing disk masses with protostellar evolution, with sufficient mass for forming giant planet cores primarily at early times.

Gl 758 B is a late-T dwarf orbiting a metal-rich Sun-like star at a projected separation of $\rho$ $\approx$ 1.6" (25 AU). We present four epochs of astrometry of this system with NIRC2 at Keck Observatory spanning 2010 to 2017 together with 630 radial velocities (RVs) of the host star acquired over the past two decades from McDonald Observatory, Keck Observatory, and the Automated Planet Finder at Lick Observatory. The RVs reveal that Gl 758 is accelerating with an evolving rate that varies between 2-5 m s$^{-1}$ yr$^{-1}$, consistent with the expected influence of the imaged companion Gl 758 B. A joint fit of the RVs and astrometry yields a dynamical mass of 42$^{+19}_{-7}$ M$_\mathrm{Jup}$ for the companion with a robust lower limit of 30.5 M$_\mathrm{Jup}$ at the 4-$\sigma$ level. Gl 758 B is on an eccentric orbit ($e$ = 0.26-0.67 at 95% confidence) with a semimajor axis of $a$ = $21.1_{-1.3}^{+2.7}$ AU and an orbital period of $P$ = $96_{-9}^{+21}$ yr, which takes it within $\approx$9 AU from its host star at periastron passage. Substellar evolutionary models generally underpredict the mass of Gl 758 B for nominal ages of 1-6 Gyr that have previously been adopted for the host star. This discrepancy can be reconciled if the system is older---which is consistent with activity indicators and recent isochrone fitting of the host star---or alternatively if the models are systematically overluminous by $\approx$0.1-0.2 dex. Gl 758 B is currently the lowest-mass directly imaged companion inducing a measured acceleration on its host star. In the future, bridging RVs and high-contrast imaging with the next generation of extremely large telescopes and space-based facilities will open the door to the first dynamical mass measurements of imaged exoplanets.

Solar-type binaries with short orbital periods ($P_{\rm close}$ $\equiv$ 1 - 10 days; $a$ $\lesssim$ 0.1 AU) cannot form directly via fragmentation of molecular clouds or protostellar disks, yet their component masses are highly correlated, suggesting interaction during the pre-main-sequence (pre-MS) phase. Moreover, the close binary fraction of pre-MS stars is consistent with that of their MS counterparts in the field ($F_{\rm close}$ = 2.1%). Thus we can infer that some migration mechanism operates during the early pre-MS phase ($\tau$ $\lesssim$ 5 Myr) that reshapes the primordial separation distribution. We test the feasibility of this hypothesis by carrying out a population synthesis calculation which accounts for two formation channels: Kozai-Lidov (KL) oscillations and dynamical instability in triple systems. Our models incorporate (1) more realistic initial conditions compared to previous studies, (2) octupole-level effects in the secular evolution, (3) tidal energy dissipation via weak-friction equilibrium tides at small eccentricities and via non-radial dynamical oscillations at large eccentricities, and (4) the larger tidal radius of a pre-MS primary. Given a 15% triple star fraction, we simulate a close binary fraction from KL oscillations alone of $F_{\rm close}$ $\approx$ 0.4% after $\tau$ = 5 Myr, which increases to $F_{\rm close}$ $\approx$ 0.8% by $\tau$ = 5 Gyr. Dynamical ejections and disruptions of unstable coplanar triples in the disk produce solitary binaries with slightly longer periods $P$ $\approx$ 10 - 100 days. The remaining $\approx$60% of close binaries with outer tertiaries, particularly those in compact coplanar configurations with log $P_{\rm out}$ (days) $\approx$ 2 - 5 ($a_{\rm out}$ $<$ 50 AU), can be explained only with substantial extra energy dissipation due to interactions with primordial gas.

We present a high angular resolution ($\sim 0.2^{\prime\prime}$), high sensitivity ($\sigma \sim 0.2$ mJy) survey of the 870 $\mu$m continuum emission from the circumstellar material around 49 pre-main sequence stars in the $\rho$ Ophiuchus molecular cloud. Because most millimeter instruments have resided in the northern hemisphere, this represents the largest high-resolution, millimeter-wave survey of the circumstellar disk content of this cloud. Our survey of 49 systems comprises 63 stars; we detect disks associated with 29 single sources, 11 binaries, 3 triple systems and 4 transition disks. We present flux and radius distributions for these systems; in particular, this is the first presentation of a reasonably complete probability distribution of disk radii at millimeter-wavelengths. We also compare the flux distribution of these protoplanetary disks with that of the disk population of the Taurus-Auriga molecular cloud. We find that disks in binaries are both significantly smaller and have much less flux than their counterparts around isolated stars. We compute truncation calculations on our binary sources and find that these disks are too small to have been affected by tidal truncation and posit some explanations for this. Lastly, our survey found 3 candidate gapped disks, one of which is a newly identified transition disk with no signature of a dip in infrared excess in extant observations.

Massive circumstellar disks are prone to gravitational instabilities, which trigger the formation of spiral arms that can fragment into bound clumps under the right conditions. Two dimensional simulations of self-gravitating disks are useful starting points for studying fragmentation, allowing for high-resolution simulations of thin disks. However, convergence issues can arise in 2D from various sources. One of these sources is the 2D approximation of self-gravity, which exaggerates the effect of self-gravity on small scales when the potential is not smoothed to account for the assumed vertical extent of the disk. This effect is enhanced by increased resolution, resulting in fragmentation at longer cooling timescales $\beta$. If true, it suggests that the 3D simulations of disk fragmentation may not have the same convergence problem and could be used to examine the nature of fragmentation without smoothing self-gravity on scales similar to the disk scale height. To that end, we have carried out local 3D self-gravitating disk simulations with simple $\beta$ cooling with fixed background irradiation to determine if 3D is necessary to properly describe disk fragmentation. Above a resolution of $\sim 40$ grid cells per scale height, we find that our simulations converge with respect to the cooling timescale. This result converges in agreement with analytic expectations which place a fragmentation boundary at $\beta_\mathrm{crit} = 3$.

The presence of a giant planet in a low-viscosity disc can create a gap edge in the disc's radial density profile sharp enough to excite the Rossby Wave Instability. This instability may evolve into dust-trapping vortices that might explain the "banana-shaped" features in recently observed asymmetric transition discs with inner cavities. Previous hydrodynamical simulations of planet-induced vortices have neglected the timescale of hundreds to thousands of orbits to grow a massive planet to Jupiter-size. In this work, we study the effect of a giant planet's runaway growth timescale on the lifetime and characteristics of the resulting vortex. For two different planet masses (1 and 5 Jupiter masses) and two different disc viscosities ($\alpha$=3$\times 10^{-4}$ and 3$\times10^{-5}$), we compare the vortices induced by planets with several different growth timescales between 10 and 4000 planet orbits. In general, we find that slowly-growing planets create significantly weaker vortices with lifetimes and surface densities reduced by more than $50\%$. For the higher disc viscosity, the longest growth timescales in our study inhibit vortex formation altogether. Additionally, slowly-growing planets produce vortices that are up to twice as elongated, with azimuthal extents well above $180^{\circ}$ in some cases. These unique, elongated vortices likely create a distinct signature in the dust observations that differentiates them from the more concentrated vortices that correspond to planets with faster growth timescales. Lastly, we find that the low viscosities necessary for vortex formation likely prevent planets from growing quickly enough to trigger the instability in self-consistent models.

The Pluto-Charon system has come into sharper focus following the fly by of New Horizons. We use N-body simulations to probe the unique dynamical history of this binary dwarf planet system. We follow the evolution of the debris disc that might have formed during the Charon-forming giant impact. First, we note that in-situ formation of the four circumbinary moons is extremely difficult if Charon undergoes eccentric tidal evolution. We track collisions of disc debris with Charon, estimating that hundreds to hundreds of thousands of visible craters might arise from 0.3--5 km radius bodies. New Horizons data suggesting a dearth of these small craters may place constraints on the disc properties. While tidal heating will erase some of the cratering history, both tidal and radiogenic heating may also make it possible to differentiate disc debris craters from Kuiper belt object craters. We also track the debris ejected from the Pluto-Charon system into the Solar System; while most of this debris is ultimately lost from the Solar System, a few tens of 10--30 km radius bodies could survive as a Pluto-Charon collisional family. Most are plutinos in the 3:2 resonance with Neptune, while a small number populate nearby resonances. We show that migration of the giant planets early in the Solar System's history would not destroy this collisional family. Finally, we suggest that identification of such a family would likely need to be based on composition as they show minimal clustering in relevant orbital parameters.

We present ALMA observations of the GQ Lup system, a young Sun-like star with a substellar mass companion in a wide-separation orbit. These observations of 870 $\mu$m continuum and CO J=3-2 line emission with beam size $\sim0.3''$ ($\sim45$ AU) resolve the disk of dust and gas surrounding the primary star, GQ Lup A, and provide deep limits on any circumplanetary disk surrounding the companion, GQ Lup b. The circumprimary dust disk is compact with a FWHM of $59\pm12$ AU, while the gas has a larger extent with a characteristic radius of $46.5\pm1.8$ AU. By forward-modeling the velocity field of the circumprimary disk based on the CO emission, we constrain the mass of GQ Lup A to be $M_* = (1.03\pm0.05)*(d/156\text{ pc})$ $M_\odot$, where $d$ is a known distance, and determine that we view the disk at an inclination angle of $60.5^\circ\pm0.5^\circ$ and a position angle of $346^\circ \pm1^\circ$. The $3\sigma$ upper limit on the 870 $\mu$m flux density of any circumplanetary disk associated with GQ Lup b of $<0.15$ mJy implies an upper limit on the dust disk mass of $<0.04$ $M_\oplus$ for standard assumptions about optically thin emission. We discuss proposed mechanisms for the formation of wide-separation substellar companions given the non-detection of circumplanetary disks around GQ Lup b and other similar systems.

Binary and multiple star systems are a frequent outcome of the star formation process, and as a result, almost half of all sun-like stars have at least one companion star. Theoretical studies indicate that there are two main pathways that can operate concurrently to form binary/multiple star systems: large scale fragmentation of turbulent gas cores and filaments or smaller scale fragmentation of a massive protostellar disk due to gravitational instability. Observational evidence for turbulent fragmentation on scales of $>$1000~AU has recently emerged. Previous evidence for disk fragmentation was limited to inferences based on the separations of more-evolved pre-main sequence and protostellar multiple systems. The triple protostar system L1448 IRS3B is an ideal candidate to search for evidence of disk fragmentation. L1448 IRS3B is in an early phase of the star formation process, likely less than 150,000 years in age, and all protostars in the system are separated by $<$200~AU. Here we report observations of dust and molecular gas emission that reveal a disk with spiral structure surrounding the three protostars. Two protostars near the center of the disk are separated by 61 AU, and a tertiary protostar is coincident with a spiral arm in the outer disk at a 183 AU separation. The inferred mass of the central pair of protostellar objects is $\sim$1 M$_{sun}$, while the disk surrounding the three protostars has a total mass of $\sim$0.30 M$_{\sun}$. The tertiary protostar itself has a minimum mass of $\sim$0.085 M$_{sun}$. We demonstrate that the disk around L1448 IRS3B appears susceptible to disk fragmentation at radii between 150~AU and 320~AU, overlapping with the location of the tertiary protostar. This is consistent with models for a protostellar disk that has recently undergone gravitational instability, spawning one or two companion stars.

We are developing a stable and precise spectrograph for the Large Binocular Telescope (LBT) named "iLocater." The instrument comprises three principal components: a cross-dispersed echelle spectrograph that operates in the YJ-bands (0.97-1.30 microns), a fiber-injection acquisition camera system, and a wavelength calibration unit. iLocater will deliver high spectral resolution (R~150,000-240,000) measurements that permit novel studies of stellar and substellar objects in the solar neighborhood including extrasolar planets. Unlike previous planet-finding instruments, which are seeing-limited, iLocater operates at the diffraction limit and uses single mode fibers to eliminate the effects of modal noise entirely. By receiving starlight from two 8.4m diameter telescopes that each use "extreme" adaptive optics (AO), iLocater shows promise to overcome the limitations that prevent existing instruments from generating sub-meter-per-second radial velocity (RV) precision. Although optimized for the characterization of low-mass planets using the Doppler technique, iLocater will also advance areas of research that involve crowded fields, line-blanketing, and weak absorption lines.

Transition disks, protoplanetary disks with inner clearings, are promising objects in which to directly image forming planets. The high contrast imaging technique of non-redundant masking is well posed to detect planetary mass companions at several to tens of AU in nearby transition disks. We present non-redundant masking observations of the T Cha and LkCa 15 transition disks, both of which host posited sub-stellar mass companions. However, due to a loss of information intrinsic to the technique, observations of extended sources (e.g. scattered light from disks) can be misinterpreted as moving companions. We discuss tests to distinguish between these two scenarios, with applications to the T Cha and LkCa 15 observations. We argue that a static, forward-scattering disk can explain the T Cha data, while LkCa 15 is best explained by multiple orbiting companions.

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

Direct imaging allows for the detection and characterization of exoplanets via their thermal emission. We report the discovery via imaging of a young Jovian planet in a triple star system and characterize its atmospheric properties through near-infrared spectroscopy. The semi-major axis of the planet is closer relative to that of its hierarchical triple star system than for any known exoplanet within a stellar binary or triple, making HD 131399 dynamically unlike any other known system. The location of HD 131399Ab on a wide orbit in a triple system demonstrates that massive planets may be found on long and possibly unstable orbits in multi-star systems. HD 131399Ab is one of the lowest mass (4+/-1 MJup) and coldest (850+/-50 K) exoplanets to have been directly imaged.

There is substantial evidence for disk formation taking place during the early stages of star formation and for most stars being born in multiple systems; however, protostellar multiplicity and disk searches have been hampered by low resolution, sample bias, and variable sensitivity. We have conducted an unbiased, high-sensitivity Karl G. Jansky Very Large Array (VLA) survey toward all known protostars (n = 94) in the Perseus molecular cloud (d~230 pc), with a resolution of ~15 AU (0.06") at 8 mm. We have detected candidate protostellar disks toward 17 sources (with 12 of those in the Class 0 stage) and we have found substructure on < 50AU scales for three Class 0 disk candidates, possibly evidence for disk fragmentation. We have discovered 16 new multiple systems (or new components) in this survey; the new systems have separations < 500 AU and 3 by < 30 AU. We also found a bi-modal distribution of separations, with peaks at ~75 AU and ~3000 AU, suggestive of formation through two distinct mechanisms: disk and turbulent fragmentation. The results from this survey demonstrate the necessity and utility of uniform, unbiased surveys of protostellar systems at millimeter and centimeter wavelengths.

Transiting circumbinary planets discovered by Kepler provide unique insight into binary star and planet formation. Several features of this new found population, for example the apparent pile-up of planets near the innermost stable orbit, may distinguish between formation theories. In this work, we determine how planet-planet scattering shapes planetary systems around binaries as compared to single stars. In particular, we look for signatures that arise due to differences in dynamical evolution in binary systems. We carry out a parameter study of N-body scattering simulations for four distinct planet populations around both binary and single stars. While binarity has little influence on the final system multiplicity or orbital distribution, the presence of a binary dramatically effects the means by which planets are lost from the system. Most circumbinary planets are lost due to ejections rather than planet-planet or planet-star collisions. The most massive planet in the system tends to control the evolution. Systems similar to the only observed multi-planet circumbinary system, Kepler-47, can arise from much more tightly packed, unstable systems. Only extreme initial conditions introduce differences in the final planet populations. Thus, we suggest that any intrinsic differences in the populations are imprinted by formation.

Low mass, self-gravitating accretion disks admit quasi-steady,`gravito-turbulent' states in which cooling balances turbulent viscous heating. However, numerical simulations show that gravito-turbulence cannot be sustained beyond dynamical timescales when the cooling rate or corresponding turbulent viscosity is too large. The result is disk fragmentation. We motivate and quantify an interpretation of disk fragmentation as the inability to maintain gravito-turbulence due to formal secondary instabilities driven by: 1) cooling, which reduces pressure support; and/or 2) viscosity, which reduces rotational support. We analyze the axisymmetric gravitational stability of viscous, non-adiabatic accretion disks with internal heating, external irradiation, and cooling in the shearing box approximation. We consider parameterized cooling functions in 2D and 3D disks, as well as radiative diffusion in 3D. We show that generally there is no critical cooling rate/viscosity below which the disk is formally stable, although interesting limits appear for unstable modes with lengthscales on the order of the disk thickness. We apply this new linear theory to protoplanetary disks subject to gravito-turbulence modeled as an effective viscosity, and cooling regulated by dust opacity. We find that viscosity renders the disk beyond $\sim 60$AU dynamically unstable on radial lengthscales a few times the local disk thickness. This is coincident with the empirical condition for disk fragmentation based on a maximum sustainable stress. We suggest turbulent stresses can play an active role in realistic disk fragmentation by removing rotational stabilization against self-gravity, and that the observed transition in behavior from gravito-turbulent to fragmenting may reflect instability of the gravito-turbulent state itself.

In the coming years, high contrast imaging surveys are expected to reveal the characteristics of the population of wide-orbit, massive, exoplanets. To date, a handful of wide planetary mass companions are known, but only one such multi-planet system has been discovered: HR8799. For low mass planetary systems, multi-planet interactions play an important role in setting system architecture. In this paper, we explore the stability of these high mass, multi-planet systems. While empirical relationships exist that predict how system stability scales with planet spacing at low masses, we show that extrapolating to super-Jupiter masses can lead to up to an order of magnitude overestimate of stability for massive, tightly packed systems. We show that at both low and high planet masses, overlapping mean motion resonances trigger chaotic orbital evolution, which leads to system instability. We attribute some of the difference in behavior as a function of mass to the increasing importance of second order resonances at high planet-star mass ratios. We use our tailored high mass planet results to estimate the maximum number of planets that might reside in double component debris disk systems, whose gaps may indicate the presence of massive bodies.

[Abridged] Star and planet formation are the complex outcomes of gravitational collapse and angular momentum transport mediated by protostellar and protoplanetary disks. In this review we focus on the role of gravitational instability in this process. We begin with a brief overview of the observational evidence for massive disks that might be subject to gravitational instability, and then highlight the diverse ways in which the instability manifests itself in protostellar and protoplanetary disks: the generation of spiral arms, small scale turbulence-like density fluctuations, and fragmentation of the disk itself. We present the analytic theory that describes the linear growth phase of the instability, supplemented with a survey of numerical simulations that aim to capture the non-linear evolution. We emphasize the role of thermodynamics and large scale infall in controlling the outcome of the instability. Despite apparent controversies in the literature, we show a remarkable level of agreement between analytic predictions and numerical results. We highlight open questions related to (1) the development of a turbulent cascade in thin disks, and (2) the role of mode-mode coupling in setting the maximum angular momentum transport rate in thick disks.

We present a multiplicity study of all known protostars (94) in the Perseus molecular cloud from a Karl G. Jansky Very Large Array (VLA) survey at Ka-band (8 mm and 1 cm) and C-band (4 cm and 6.6 cm). The observed sample has a bolometric luminosity range between 0.1 L$_{\odot}$ and $\sim$33 L$_{\odot}$, with a median of 0.7 L$_{\odot}$. This multiplicity study is based on the Ka-band data, having a best resolution of $\sim$0.065" (15 AU) and separations out to $\sim$43" (10000 AU) can be probed. The overall multiplicity fraction (MF) is found to be of 0.40$\pm$0.06 and the companion star fraction (CSF) is 0.71$\pm$0.06. The MF and CSF of the Class 0 protostars are 0.57$\pm$0.09 and 1.2$\pm$0.2, and the MF and CSF of Class I protostars are both 0.23$\pm$0.08. The distribution of companion separations appears bi-modal, with a peak at $\sim$75 AU and another peak at $\sim$3000 AU. Turbulent fragmentation is likely the dominant mechanism on $>$1000 AU scales and disk fragmentation is likely to be the dominant mechanism on $<$200 AU scales. Toward three Class 0 sources we find companions separated by $<$30 AU. These systems have the smallest separations of currently known Class 0 protostellar binary systems. Moreover, these close systems are embedded within larger (50 AU to 400 AU) structures and may be candidates for ongoing disk fragmentation.

Recently, Gauza et al. (2015) reported the discovery of a companion to the late M-dwarf, VHS J125601.92-125723.9 (VHS 1256-1257). The companion's absolute photometry suggests its mass and atmosphere are similar to the HR 8799 planets. However, as a wide companion to a late-type star, it is more accessible to spectroscopic characterization. We discovered that the primary of this system is an equal-magnitude binary. For an age $\sim300$ Myr the A and B components each have a mass of $64.6^{+0.8}_{-2.0}~M_{\mathrm{Jup}}$, and the b component has a mass of $11.2^{+9.7}_{-1.8}$, making VHS 1256-1257 only the third brown dwarf triple system. There exists some tension between the spectrophotometric distance of $17.2\pm2.6$ pc and the parallax distance of $12.7\pm1.0$ pc. At 12.7 pc VHS1256-1257 A and B would be the faintest known M7.5 objects, and are even faint outliers among M8 types. If the larger spectrophotmetric distance is more accurate than the parallax, then the mass of each component increases. In particular, the mass of the b component increases well above the deuterium burning limit to $\sim35~M_{\mathrm{Jup}}$ and the mass of each binary component increases to $73^{+20}_{-17}~M_{\mathrm{Jup}}$. At 17.1 pc, the UVW kinematics of the system are consistent with membership in the AB~Dor moving group. The architecture of the system resembles a hierarchical stellar multiple suggesting it formed via an extension of the star-formation process to low masses. Continued astrometric monitoring will resolve this distance uncertainty and will provide dynamical masses for a new benchmark system.

We present the first dust emission results toward a sample of seven protostellar disk candidates around Class 0 and I sources in the Perseus molecular cloud from the VLA Nascent Disk and Multiplicity (VANDAM) survey with ~0.05'' or 12 AU resolution. To examine the surface brightness profiles of these sources, we fit the Ka-band 8 mm dust-continuum data in the u,v-plane to a simple, parametrized model based on the Shakura-Sunyaev disk model. The candidate disks are well-fit by a model with a disk-shaped profile and have masses consistent with known Class 0 and I disks. The inner-disk surface densities of the VANDAM candidate disks have shallower density profiles compared to disks around more evolved Class II systems. The best-fit model radii of the seven early-result candidate disks are R_c > 10 AU; at 8 mm, the radii reflect lower limits on the disk size since dust continuum emission is tied to grain size and large grains radially drift inwards. These relatively large disks, if confirmed kinematically, are inconsistent with theoretical models where the disk size is limited by strong magnetic braking to < 10 AU at early times.

We present the first results from our Keck program investigating the orbital architectures of planet-hosting multiple star systems. Kepler-444 is a metal-poor triple star system that hosts five sub-Earth-sized planets orbiting the primary star (Kepler-444A), as well as a spatially unresolved pair of M dwarfs (Kepler-444BC) at a projected distance of 1.8" (66 AU). We combine our Keck/NIRC2 adaptive optics astrometry with multi-epoch Keck/HIRES RVs of all three stars to determine a precise orbit for the BC pair around A, given their empirically constrained masses. We measure minimal astrometric motion ($1.0\pm0.6$ mas yr$^{-1}$, or $0.17\pm0.10$ km s$^{-1}$), but our RVs reveal significant orbital velocity ($1.7\pm0.2$ km s$^{-1}$) and acceleration ($7.8\pm0.5$ m s$^{-1}$ yr$^{-1}$). We determine a highly eccentric stellar orbit ($e=0.864\pm0.023$) that brings the tight M dwarf pair within $5.0^{+0.9}_{-1.0}$ AU of the planetary system. We validate that the system is dynamically stable in its present configuration via n-body simulations. We find that the A$-$BC orbit and planetary orbits are likely aligned (98%) given that they both have edge-on orbits and misalignment induces precession of the planets out of transit. We conclude that the stars were likely on their current orbits during the epoch of planet formation, truncating the protoplanetary disk at $\approx$2 AU. This truncated disk would have been severely depleted of solid material from which to form the total $\approx$1.5 $M_{\rm Earth}$ of planets. We thereby strongly constrain the efficiency of the conversion of dust into planets and suggest that the Kepler-444 system is consistent with models that explain the formation of more typical close-in Kepler planets in normal, not truncated, disks.

Exoplanet detections have revolutionized astronomy, offering new insights into solar system architecture and planet demographics. While nearly 1900 exoplanets have now been discovered and confirmed, none are still in the process of formation. Transition discs, protoplanetary disks with inner clearings best explained by the influence of accreting planets, are natural laboratories for the study of planet formation. Some transition discs show evidence for the presence of young planets in the form of disc asymmetries or infrared sources detected within their clearings, as in the case of LkCa 15. Attempts to observe directly signatures of accretion onto protoplanets have hitherto proven unsuccessful. Here we report adaptive optics observations of LkCa 15 that probe within the disc clearing. With accurate source positions over multiple epochs spanning 2009 - 2015, we infer the presence of multiple companions on Keplerian orbits. We directly detect H{\alpha} emission from the innermost companion, LkCa 15 b, evincing hot (~10,000 K) gas falling deep into the potential well of an accreting protoplanet.

We present continuum and molecular line observations at 230 GHz and 345 GHz from the Sub-millimeter Array (SMA) toward three protostars in the Perseus L1448N region. The data are from the large project "Mass Assembly of Stellar Systems and their Evolution with the SMA" (MASSES). Three dust continuum sources, Source B, Source NW, and Source A, are detected at both frequencies. These sources have corresponding emission peaks in C18O (J=2-1), 13CO (J=2-1), and HCO+ (J=4-3), and have offsets with N2D+ (J=3-2) peaks. High angular resolution data from a complimentary continuum survey with the Karl G. Jansky Very Large Array show that Source B is associated with three 8 mm continuum objects, Source NW with two, and Source A remains single. These results suggest that multiplicity in L1448N exists at different spatial scales from a few thousand AU to < 100 AU. Velocity gradients in each source obtained from two-dimensional fits to the SMA C18O emission are found to be perpendicular to within 20 degrees of the outflow directions as revealed by 12CO (J=2-1). We have observed that Sources B and NW with multiplicity have higher densities than Source A without multiplicity. This suggests that thermal Jeans fragmentation can be relevant in the fragmentation process. However, we have not observed a difference in the ratio between rotational and gravitational energy between sources with and without multiplicity. We also have not observed a trend between non-thermal velocity dispersions and the level of fragmentation. Our study has provided the first direct and comprehensive comparison between multiplicity and core properties in low-mass protostars, although based on small number statistics.

Magnetic fields can regulate disk formation, accretion and jet launching. Until recently, it has been difficult to obtain high resolution observations of the magnetic fields of the youngest protostars in the critical region near the protostar. The VANDAM survey is observing all known protostars in the Perseus Molecular Cloud. Here we present the polarization data of IRAS 4A. We find that with ~ 0.2'' (50 AU) resolution at {\lambda} = 8.1 and 10.3 mm, the inferred magnetic field is consistent with a circular morphology, in marked contrast with the hourglass morphology seen on larger scales. This morphology is consistent with frozen-in field lines that were dragged in by rotating material entering the infall region. The field morphology is reminiscent of rotating circumstellar material near the protostar. This is the first polarization detection of a protostar at these wavelengths. We conclude from our observations that the dust emission is optically thin with {\beta} ~ 1.3, suggesting that mm/cm-sized grains have grown and survived in the short lifetime of the protostar.

Ionizing stellar photons heat the upper regions of planetary atmospheres, driving atmospheric mass loss. Gas escaping from several hot, hydrogen-rich planets has been detected using UV and X-ray transmission spectroscopy. Because these planets are tidally locked, and thus asymmetrically irradiated, escaping gas is unlikely to be spherically symmetric. In this paper, we focus on the effects of asymmetric heating on local outflow structure. We use the Athena code for hydrodynamics to produce 3D simulations of hot Jupiter mass loss that jointly model wind launching and stellar heating via photoionization. Our fiducial planet is an inflated, hot Jupiter with radius $R_p=2.14 R_{\rm Jup}$ and mass $M_p = 0.53 M_{\rm Jup}$. We irradiate the initially neutral, atomic hydrogen atmosphere with 13.6 eV photons and compute the outflow's ionization structure. There are clear asymmetries in the atmospheric outflow, including a neutral shadow on the planet's nightside. Given an incident ionizing UV flux comparable to that of the Sun, we find a steady-state mass loss rate of ~$2\times10^{10}$ g s$^{-1}$. The total mass loss rate and the outflow substructure along the substellar ray show good agreement with earlier 1D models, for two different fluxes. Our 3D data cube can be used to generate the outflow's extinction spectrum during transit. As a proof of concept, we find absorption of stellar Lyman-alpha at Doppler-shifted velocities of up to $\pm 50$ km s$^{-1}$. Our work provides a starting point for further 3D models that can be used to predict observable signatures of hot Jupiter mass loss.

General-purpose, moving-mesh schemes for hydrodynamics have opened the possibility of combining the accuracy of grid-based numerical methods with the flexibility and automatic resolution adaptivity of particle-based methods. Due to their supersonic nature, Keplerian accretion discs are in principle a very attractive system for applying such freely moving mesh techniques. However, the high degree of symmetry of simple accretion disc models can be difficult to capture accurately by these methods, due to the generation of geometric grid noise and associated numerical diffusion, which is absent in polar grids. To explore these and other issues, in this work we study the idealized problem of two-dimensional planet-disc interaction with the moving-mesh code AREPO. We explore the hydrodynamic evolution of discs with planets through a series of numerical experiments that vary the planet mass, the disc viscosity and the mesh resolution, and compare the resulting surface density, vortensity field and tidal torque with results from the literature. We find that the performance of the moving-mesh code in this problem is in accordance with published results, showing good consistency with grid codes written in polar coordinates. We also conclude that grid noise and mesh distortions do not introduce excessive numerical diffusion. Finally, we show how the moving-mesh approach can naturally increase resolution in regions of high densityaround planets and planetary wakes, while retaining the background flow at low resolution. This provides an alternative to the difficult task of implementing adaptive mesh refinement in conventional polar-coordinate codes.

The formation and early evolution of circumstellar discs often occurs within dense, newborn stellar clusters. For the first time, we apply the moving-mesh code AREPO, to circumstellar discs in 3-D, focusing on disc-disc interactions that result from stellar fly-bys. Although a small fraction of stars are expected to undergo close approaches, the outcomes of the most violent encounters might leave an imprint on the discs and host stars that will influence both their orbits and their ability to form planets. We first construct well-behaved 3-D models of self-gravitating discs, and then create a suite of numerical experiments of parabolic encounters, exploring the effects of pericenter separation r_p, disc orientation and disc-star mass ratio (M_d/M_*) on the orbital evolution of the host stars. Close encounters (2r_p<~ disc radius) can truncate discs on very short time scales. If discs are massive, close encounters facilitate enough orbital angular momentum extraction to induce stellar capture. We find that for realistic primordial disc masses M_d<~0.1M_*, non-colliding encounters induce minor orbital changes, which is consistent with analytic calculations of encounters in the linear regime. The same disc masses produce entirely different results for grazing/colliding encounters. In the latter case, rapidly cooling discs lose orbital energy by radiating away the energy excess of the shock-heated gas, thus causing capture of the host stars into a bound orbit. In rare cases, a tight binary with a circumbinary disc forms as a result of this encounter.

We are conducting a Jansky VLA Ka-band (8 mm and 1 cm) and C-band (4 cm and 6.4 cm) survey of all known protostars in the Perseus Molecular Cloud, providing resolution down to $\sim$0.06'' and $\sim$0.35" in Ka-band and C-band, respectively. Here we present first results from this survey that enable us to examine the source NGC 1333 IRAS2A in unprecedented detail and resolve it into a proto-binary system separated by 0.621"$\pm$0.006" ($\sim$143 AU) at 8 mm, 1 cm, and 4 cm. These 2 sources (IRAS2A VLA1 and VLA2) are likely driving the two orthogonal outflows known to originate from IRAS2A. The brighter source IRAS2A VLA1 is extended perpendicular to its outflow in the VLA data, with a deconvolved size of 0.055" ($\sim$13 AU), possibly tracing a protostellar disk. The recently reported candidate companions (IRAS2A MM2 and MM3) are not detected in either our VLA data, CARMA 1.3 mm data, or SMA 850 $\mu$m data. SMA CO ($J=3\rightarrow2$), CARMA CO ($J=2\rightarrow1$), and lower resolution CARMA CO ($J=1\rightarrow0$) observations are used to examine the outflow origins and the nature of the candidate companions to IRAS2A VLA1. The CO ($J=3\rightarrow2$) and ($J=2\rightarrow1$) data show that IRAS2A MM2 is coincident with a bright CO emission spot in the east-west outflow, and IRAS2A MM3 is within the north-south outflow. In contrast, IRAS2A VLA2 lies at the east-west outflow symmetry point. We propose that IRAS2A VLA2 is the driving source of the East-West outflow and a true companion to IRAS2A VLA1, whereas IRAS2A MM2 and MM3 may not be protostellar.

Using recently established empirical calibrations for the impact of chromospheric activity on the radii, effective temperatures, and estimated masses of active low-mass stars and brown dwarfs, we reassess the shape of the initial mass function (IMF) across the stellar/substellar boundary in the Upper Sco star-forming region (age 5-10 Myr). We adjust the observed effective temperatures to warmer values using the observed strength of the chromospheric H$\alpha$ emission, and redetermine the estimated masses of objects using pre--main-sequence evolutionary tracks in the H-R diagram. The effect of the activity-adjusted temperatures is to shift the objects to higher masses by 3-100%. While the slope of the resulting IMF at substellar masses is not strongly changed, the peak of the IMF does shift from ~0.06 to ~0.11 Msun. Moreover, for objects with masses <~0.2 Msun, the ratio of brown dwarfs to stars changes from ~80% to ~33%. These results suggest that activity corrections are essential for studies of the substellar mass function, if the masses are estimated from spectral types or from effective temperatures.

Multiplicity is one of the most fundamental observable properties of massive O-type stars and offers a promising way to discriminate between massive star formation theories. Nevertheless, companions at separations between 1 and 100 mas remain mostly unknown due to intrinsic observational limitations. [...] The Southern MAssive Stars at High angular resolution survey (SMASH+) was designed to fill this gap by providing the first systematic interferometric survey of Galactic massive stars. We observed 117 O-type stars with VLTI/PIONIER and 162 O-type stars with NACO/SAM, respectively probing the separation ranges 1-45 and 30-250mas and brightness contrasts of Delta H < 4 and Delta H < 5. Taking advantage of NACO's field-of-view, we further uniformly searched for visual companions in an 8''-radius down to Delta H = 8. This paper describes the observations and data analysis, reports the discovery of almost 200 new companions in the separation range from 1mas to 8'' and presents the catalog of detections, including the first resolved measurements of over a dozen known long-period spectroscopic binaries. Excluding known runaway stars for which no companions are detected, 96 objects in our main sample (DEC < 0 deg; H<7.5) were observed both with PIONIER and NACO/SAM. The fraction of these stars with at least one resolved companion within 200mas is 0.53. Accounting for known but unresolved spectroscopic or eclipsing companions, the multiplicity fraction at separation < 8'' increases to f_m = 0.91 +/- 0.03. The fraction of luminosity class V stars that have a bound companion reaches 100% at 30mas while their average number of physically connected companions within 8'' is f_c = 2.2 +/- 0.3. This demonstrates that massive stars form nearly exclusively in multiple systems. Additionally, the nine non-thermal (NT) radio emitters observed by SMASH+ are all resolved [...]

The recent discovery of planets orbiting main sequence binaries will provide crucial constraints for theories of binary and planet formation. The formation pathway for these planets is complicated by uncertainties in the formation mechanism of the host stars. In this paper, we compare the dynamical states of single and binary star planetary systems. Specifically, we pose two questions: (1) What does it mean for a circumbinary system to be dynamically packed? (2) How many systems are required to differentiate between a population of packed or sparse planets? We determine when circumbinary systems become dynamically unstable as a function of the separation between the host-stars and the inner planet, and the first and second planets. We show that these represent unique stability constraints compared to single-star systems. We find that although the existing Kepler data is insufficient to distinguish between a population of packed or sparse circumbinary systems, a more thorough study of circumbinary TTVs combined with an order of magnitude increase in the number of systems may prove conclusive. Future space missions such as TESS provide the best opportunity for increasing the sample size.

We show that the Kepler spacecraft in two-reaction wheel mode of operation is very well suited for the study of eclipsing binary star systems. Continued observations of the Kepler field will provide the most enduring and long-term valuable science. It will enable the discovery and characterization of eclipsing binaries with periods greater than 1 year - these are the most important, yet least understood binaries for habitable-zone planet background considerations. The continued mission will also enable the investigation of hierarchical multiple systems (discovered through eclipse timing variations), and provide drastically improved orbital parameters for circumbinary planetary systems.

The majority of solar-type stars reside in multiple systems, especially binaries. They form in dense cores of molecular clouds that are observed to be significantly magnetized. Our previous study shows that magnetic braking can tighten the binary separation during the protostellar mass accretion phase by removing the angular momentum of the accreting material. Recent numerical calculations of single star formation have shown that misalignment between the magnetic field and rotation axis may weaken both magnetic braking and the associated magnetically driven outflows. These two effects allow for disk formation even in strongly magnetized cores. Here we investigate the effects of magnetic field misalignment on the properties of protobinaries. Somewhat surprisingly, the misaligned magnetic field is more efficient at tightening the binary orbit compared to the aligned field. The main reason is that the misalignment weakens the magnetically-driven outflow, which allows more material to accrete onto the binary. Even though the specific angular momentum of this inner material is higher than in the aligned case, it is insufficient to compensate for the additional mass. A corollary of this result is that a weaker field is required to achieve the same degree of inward migration when the field is tilted relative to the rotation axis. Large field misalignment also helps to produce rotationally-supported circumbinary disks even for relatively strong magnetic fields, by weakening the magnetically-dominated structure close to the binary. Our result may provide an explanation for the circumbinary disks detected in recent SMA and ALMA observations.

Studies of solar-type binaries have found coplanarity between the equatorial and orbital planes of systems with $<$40 AU separation. By comparison, the alignment of the equatorial and orbital axes in the substellar regime, and the associated implications for formation theory, are relatively poorly constrained. Here we present the discovery of the rotation period of 3.32 $\pm$ 0.15 hours from 2MASS J0746+20A - the primary component of a tight (2.7 AU) ultracool dwarf binary system (L0+L1.5). The newly discovered period, together with the established period via radio observations of the other component, and the well constrained orbital parameters and rotational velocity measurements, allow us to infer alignment of the equatorial planes of both components with the orbital plane of the system to within 10 degrees. This result suggests that solar-type binary formation mechanisms may extend down into the brown dwarf mass range, and we consider a number of formation theories that may be applicable in this case. This is the first such observational result in the very low mass binary regime. In addition, the detected period of 3.32 $\pm$ 0.15 hours implies that the reported radio period of 2.07 $\pm$ 0.002 hours is associated with the secondary star, not the primary, as was previously claimed. This in turn refutes the claimed radius of 0.78 $\pm$ 0.1 $R_{J}$ for 2MASS J0746+20A, which we demonstrate to be 0.99 $\pm$ 0.03 $R_{J}$.

Physical collisions and close approaches between stars play an important role in the formation of exotic stellar systems. Standard theories suggest that collisions are rare, occurring only via random encounters between stars in dense clusters. We present a different formation pathway, the triple evolution dynamical instability (TEDI), in which mass loss in an evolving triple star system causes orbital instability. The subsequent chaotic orbital evolution of the stars triggers close encounters, collisions, exchanges between the stellar components, and the dynamical formation of eccentric compact binaries (including Sirius like binaries). We demonstrate that the rate of stellar collisions due to the TEDI is approximately 10^{-4} yr^{-1} per Milky-Way Galaxy, which is nearly 30 times higher than the total collision rate due to random encounters in the Galactic globular clusters. Moreover, we find that the dominant type of stellar collisions is qualitatively different; most collisions involve asymptotic giant branch stars, rather than main sequence, or slightly evolved stars, which dominate collisions in globular clusters. The TEDI mechanism should lead us to revise our understanding of collisions and the formation of compact, eccentric binaries in the field.

We present empirical relations for determining the amount by which the effective temperatures and radii---and therefore the estimated masses---of low-mass stars and brown dwarfs are altered due to chromospheric activity. Accurate estimates of stellar radii are especially important in the context of searches for transiting exoplanets, which rely upon the assumed stellar radius/density to infer the planet radius/density. Our relations are based on a large set of well studied low-mass stars in the field and on a set of benchmark low-mass eclipsing binaries. The relations link the amount by which an active object's temperature is suppressed, and its radius inflated, to the strength of its Halpha emission. These relations are found to approximately preserve bolometric luminosity. We apply these relations to the peculiar brown-dwarf eclipsing binary 2M0535-05, in which the active, higher-mass brown dwarf has a cooler temperature than its inactive, lower-mass companion. The relations correctly reproduce the observed temperatures and radii of 2M0535-05 after accounting for the Halpha emission; 2M0535-05 would be in precise agreement with theoretical isochrones were it inactive. The relations that we present are applicable to brown dwarfs and low-mass stars with masses below 0.8 Msun and for which the activity, as measured by Halpha, is in the range -4.6 < log Lha/Lbol < -3.3. We expect these relations to be most useful for correcting radius and mass estimates of low-mass stars and brown dwarfs over their active lifetimes (few Gyr). We also discuss the implications of this work for determinations of young cluster IMFs.

We present empirical relations for determining the amount by which the effective temperatures and radii -- and therefore the estimated masses -- of low-mass stars and brown dwarfs are altered due to chromospheric activity. We base our relations on a large set of low-mass stars in the field with Halpha activity measurements, and on a set of low-mass eclipsing binaries with X-ray activity measurements from which we indirectly infer the Halpha activity. Both samples yield consistent relations linking the amount by which an active object's temperature is suppressed, and its radius inflated, to the strength of its Halpha emission. These relations are found to approximately preserve bolometric luminosity. We apply these relations to the peculiar brown-dwarf eclipsing binary 2M0535-05, in which the active, higher-mass brown dwarf has a cooler temperature than its inactive, lower-mass companion. The relations correctly reproduce the observed temperatures and radii of 2M0535-05 after accounting for the Halpha emission; 2M0535-05 would be in precise agreement with theoretical isochrones were it inactive. The relations that we present are applicable to brown dwarfs and low-mass stars with masses below 0.8 Msun and for which the activity, as measured by the fractional Halpha luminosity, is in the range -4.6 < log LHa/Lbol < -3.3. We expect these relations to be most useful for correcting radius and mass estimates of low-mass stars and brown dwarfs over their active lifetimes (few Gyr) and when the ages or distances (and therefore luminosities) are unknown. We also discuss the implications of this work for improved determinations of young cluster initial mass functions.

The Pluto system provides a unique local laboratory for the study of binaries with multiple low mass companions. In this paper, we study the orbital stability of P4, the most recently discovered moon in the Pluto system. This newfound companion orbits near the plane of the Pluto-Charon binary, roughly halfway between the two minor moons Nix and Hydra. We use a suite of few body integrations to constrain the masses of Nix and Hydra, and the orbital parameters of P4. For the system to remain stable over the age of the Solar System, the masses of Nix and Hydra likely do not exceed 5e16 kg and 9e16 kg, respectively. These upper limits assume a fixed mass ratio between Nix and Hydra at the value implied by their median optical brightness. Our study finds that stability is more sensitive to their total mass and that a downward revision of Charon's eccentricity (from our adopted value of 0.0035) is unlikely to significantly affect our conclusions. Our upper limits are an order of magnitude below existing astrometric limits on the masses of Nix and Hydra. For a density at least that of ice, the albedos of Nix and Hydra would exceed 0.3. This constraint implies they are icy, as predicted by giant impact models. Even with these low masses, P4 only remains stable if its eccentricity e < 0.02. The 5:1 commensurability with Charon is particularly unstable, Combining stability constraints with the observed mean motion places the preferred orbit for P4 just exterior to the 5:1 resonance. These predictions will be tested when the New Horizons satellite visits Pluto. Based on the results for the Pluto-Charon system, we expect that circumbinary, multi-planet systems will be more widely spaced than their singleton counterparts. Further, circumbinary exoplanets close to the three-body stability boundary, such as those found by Kepler, are less likely to have other companions nearby.

Approximately 30% of all massive stars in the Galaxy are runaways with velocities exceeding 30 km/s. Their high speeds allow them to travel ~0.1-1 kpc away from their birth place before they explode at the end of their several Myr lifetimes. At high redshift, when galaxies were much smaller than in the local universe, runaways could venture far from the dense inner regions of their host galaxies. From these large radii, and therefore low column densities, much of their ionizing radiation is able to escape into the intergalactic medium. Runaways may therefore significantly enhance the overall escape fraction of ionizing radiation, fesc, from small galaxies at high redshift. We present simple models of the high-redshift runaway population and its impact on fesc as a function of halo mass, size, and redshift. We find that the inclusion of runaways enhances fesc by factors of ~1.1-8, depending on halo mass, galaxy geometry, and the mechanism of runaway production, implying that runaways may contribute 50-90% of the total ionizing radiation escaping from high-redshift galaxies. Runaways may therefore play an important role in reionizing the universe.

Many planets are observed in stellar binary systems, and their frequency may be comparable to that of planetary systems around single stars. Binary stellar evolution in such systems influences the dynamical evolution of the resident planets. Here we study the evolution of a single planet orbiting one star in an evolving binary system. We find that stellar evolution can trigger dynamical instabilities that drive planets into chaotic orbits. This instability leads to planet-star collisions, exchange of the planet between the binary stars ("star-hoppers"), and ejection of the planet from the system. The means by which planets can be recaptured is similar to the pull-down capture mechanism for irregular solar system satellites. Because planets often suffer close encounters with the primary on the asymptotic giant branch, captures during a collision with the stellar envelope are also possible. Such capture could populate the habitable zone around white dwarfs.

Binary stars produce an array of dramatic astrophysical phenomena. They allow us to probe stellar structure, nuclear physics, and gravitational wave physics. They also produce the powerful supernovae that allow us to measure the scale of the universe. Despite their importance and ubiquity, many questions remain unanswered as to how the star formation process produces the wide array of stellar multiples that we observe. A complete model for binary formation encompasses three main components. We must know: (1) the primordial population of systems, (2) the influence of the dynamical processes that reshape this distribution as stars form and natal clusters disperse, and (3) the role of binary stellar evolution. In this article I review the most prominent theories for binary formation: turbulent core fragmentation, disk fragmentation, and competitive accretion. I argue that turbulent core fragmentation at all masses, with disk fragmentation added in at the upper end of the mass spectrum, might explain the trend towards increasing multiplicity for higher mass stars. In addition, I provide a brief overview of the observational statistics and of some of the important processes that modify the primordial distribution of stellar orbits.

Accretion disks that become gravitationally unstable can fragment into stellar or sub-stellar companions. The formation and survival of these fragments depends on the precarious balance between self-gravity, internal pressure, tidal shearing, and rotation. Disk fragmentation depends on two key factors (1) whether the disk can get to the fragmentation boundary of Q=1, and (2) whether fragments can survive for many orbital periods. Previous work suggests that to reach Q=1, and have fragments survive, a disk must cool on an orbital timescale. Here we show that disks heated primarily by external irradiation always satisfy the standard cooling time criterion. Thus even though irradiation heats disks, and makes them more stable in general, once they reach the fragmentation boundary, they fragment more easily. We derive a new cooling criterion that determines fragment survival, and calculate a pressure modified Hill radius, which sets the maximum size of pressure-supported objects in a Keplerian disk. We conclude that fragmentation in protostellar disks might occur at slightly smaller radii than previously thought, and recommend tests for future simulations that will better predict the outcome of fragmentation in real disks.

Young protostars embedded in circumstellar discs accrete from an angular momentum-rich mass reservoir. Without some braking mechanism, all stars should be spinning at or near break-up velocity. In this paper, we perform simulations of the self-gravitational collapse of an isothermal cloud using the ORION adaptive mesh refinement code and investigate the role that gravitational torques might play in the spin-down of the dense central object. While magnetic effects likely dominate for low mass stars, high mass and Population III stars might be less well magnetised. We find that gravitational torques alone prevent the central object from spinning up to more than half of its breakup velocity, because higher rotation rates lead to bar-like deformations that enable efficient angular momentum transfer to the surrounding medium. We also find that the long-term spin evolution of the central object is dictated by the properties of the surrounding disc. In particular, spiral modes with azimuthal wavenumber $m=2$ couple more effectively to its spin than the lopsided $m=1$ mode, which was found to inhibit spin evolution. We suggest that even in the absence of magnetic fields, gravitational torques may provide an upper limit on stellar spin, and that moderately massive circumstellar discs can cause long-term spin down.

We characterize the infall rate onto protostellar systems forming in self-gravitating radiation-hydrodynamic simulations. Using two dimensionless parameters to determine disks' susceptability to gravitational fragmentation, we infer limits on protostellar system multiplicity and the mechanism of binary formation. We show that these parameters give robust predictions even in the case of marginally resolved protostellar disks. We find that protostellar systems with radiation feedback predominately form binaries via turbulent fragmentation, not disk instability, and we predict turbulent fragmentation is the dominant channel for binary formation for low-mass stars. We clearly demonstrate that systems forming in simulations including radiative feedback have fundamentally different parameters than those in purely hydrodynamic simulations.

Recent direct imaging discoveries suggest a new class of massive, distant planets around A stars. These widely separated giants have been interpreted as signs of planet formation driven by gravitational instability, but the viability of this mechanism is not clear cut. In this paper, we first discuss the local requirements for fragmentation and the initial fragment mass scales. We then consider whether the fragment's subsequent growth can be terminated within the planetary mass regime. Finally, we place disks in the larger context of star formation and disk evolution models. We find that in order for gravitational instability to produce planets, disks must be atypically cold in order to reduce the initial fragment mass. In addition, fragmentation must occur during a narrow window of disk evolution, after infall has mostly ceased, but while the disk is still sufficiently massive to undergo gravitational instability. Under more typical conditions, disk-born objects will likely grow well above the deuterium burning planetary mass limit. We conclude that if planets are formed by gravitational instability, they must be the low mass tail of the distribution of disk-born companions. To validate this theory, on-going direct imaging surveys must find a greater abundance of brown dwarf and M-star companions to A-stars. Their absence would suggest planet formation by a different mechanism such as core accretion, which is consistent with the debris disks detected in these systems.

We study rapidly accreting, gravitationally unstable disks with a series of global, three dimensional, numerical experiments using the code ORION. In this paper we conduct a numerical parameter study focused on protostellar disks, and show that one can predict disk behavior and the multiplicity of the accreting star system as a function of two dimensionless parameters which compare the disk's accretion rate to its sound speed and orbital period. Although gravitational instabilities become strong, we find that fragmentation into binary or multiple systems occurs only when material falls in several times more rapidly than the canonical isothermal limit. The disk-to-star accretion rate is proportional to the infall rate, and governed by gravitational torques generated by low-m spiral modes. We also confirm the existence of a maximum stable disk mass: disks that exceed ~50% of the total system mass are subject to fragmentation and the subsequent formation of binary companions.

Most analytic work to date on protostellar disks has focused on those in isolation from their environments. However, observations are now beginning to probe the earliest, most embedded phases of star formation, during which disks are rapidly accreting from their parent cores and cannot be modeled in isolation. We present a simple, one-zone model of protostellar accretion disks with high mass infall rates. Our model combines a self-consistent calculation of disk temperatures with an approximate treatment of angular momentum transport via two mechanisms. We use this model to survey the properties of protostellar disks across a wide range of stellar masses and evolutionary times, and make predictions for disks' masses, sizes, spiral structure, and fragmentation that will be directly testable by future large-scale surveys of deeply embedded disks. We define a dimensionless accretion-rotation parameter which, in conjunction with the disk's temperature, controls the disk evolution. We track the dominant mode of angular momentum transport, and demonstrate that for stars with final masses greater than roughly one solar mass, gravitational instabilities are the most important mechanism as most of the mass accumulates. We predict that binary formation through disk fission, fragmentation of the disk into small objects, and spiral arm strength all increase in importance to higher stellar masses.

Recent advances in our understanding of massive star formation have made clear the important role of protostellar disks in mediating accretion. Here we describe a simple, semi-analytic model for young, deeply embedded, massive accretion disks. Our approach enables us to sample a wide parameter space of stellar mass and environmental variables, providing a means to make predictions for a variety of sources that next generation telescopes like ALMA and the EVLA will observe. Moreover we include, at least approximately, multiple mechanisms for angular momentum transport, a comprehensive model for disk heating and cooling, and a realistic estimate for the angular momentum in the gas reservoir. We make predictions for the typical sizes, masses, and temperatures of the disks, and describe the role of gravitational instabilities in determining the binarity fraction and upper mass cut-off.

We examine whether massive-star accretion disks are likely to fragment due to self-gravity. Rapid accretion and high angular momentum push these disks toward fragmentation, whereas viscous heating and the high protostellar luminosity stabilize them. We find that for a broad range of protostar masses and for reasonable accretion times, massive disks larger than ~150 AU are prone to fragmentation. We develop an analytical estimate for the angular momentum of accreted material, extending the analysis of Matzner and Levin (2005) to account for strongly turbulent initial conditions. In a core-collapse model, we predict that disks are marginally prone to fragmentation around stars of about four to 15 solar masses -- even if we adopt conservative estimates of the disks' radii and tendency to fragment. More massive stars are progressively more likely to fragment, and there is a sharp drop in the stability of disk accretion at the very high accretion rates expected above 110 solar masses. Fragmentation may starve accretion in massive stars, especially above this limit, and is likely to create swarms of small, coplanar companions.

The major obstacle to the direct detection of companions to nearby stars is the overwhelming brightness of the host star. Current instruments employing the combination of adaptive optics (AO) and coronagraphy can typically detect objects within 2'' of the star that are 10^{4-5} times fainter. Correlated speckle noise is one of the biggest obstacles limiting such high-contrast imaging. We have obtained a series of 284 8 s, AO-corrected, coronagraphically occulted H-band images of the star Vega at the 3.63 m AEOS telescope located on Haleakala, Hawaii. This dataset is unique for studying the temporal behavior of speckle noise, and represents the first time such a study on highly corrected _coronagraphic_ AO images has been carried out in a quantitative way. We find the speckle pattern to be highly stable in both position and time in our data. This is due to the fact that the AO system corrects disturbances to the stellar wave front at the level where the instrumental wave front errors dominate. Because of this, we find that our detection limit is not significantly improved simply with increased exposure time alone. However, we are able to improve our dynamic range by 1.5-2 magnitudes through subtraction of static/quasi-static speckles in two rotating frames: the telescope pupil frame and the deformable mirror frame. Furthermore, from our data, we are able to constrain the mass of any purported companion to Vega to be less than 45 M_J at 8 AU and less than 30 M_J at 16 AU, radii not previously probed at these sensitivities.