Small planets are common around late-M dwarfs and can be detected through highly precise photometry by the transit method. Planets orbiting nearby stars are particularly important as they are often the best-suited for future follow-up studies. We present observations of three nearby M-dwarfs referred to as EIC-1, EIC-2, and EIC-3, and use them to search for transits and set limits on the presence of planets. On most nights our observations are sensitive to Earth-sized transiting planets, and photometric precision is similar to or better than TESS for faint late-M dwarfs of the same magnitude (I=15 mag). We present our photometry and transit search pipeline, which utilizes simple median detrending in combination with transit least squares based transit detection (Hippke & Heller 2019).For these targets, and transiting planets between one and two Earth radii, we achieve an average transit detection probability of 60% between periods of 0.5 and 2 days, 30% between 2 and 5 days,and 10% between 5 and 10 days. These sensitivities are conservative compared to visual searches.

The formation of planetesimals is an unsolved problem in planet formation theory. A prominent scenario for overcoming dust growth barriers in dead zones is the gravitational collapse of locally over-dense regions, shown to robustly produce $\sim$100 km sized objects. Still, the conditions under which planetesimal formation occurs remain unclear. For collapse to proceed, the self-gravity of an overdensity must overcome stellar tidal disruption on large scales and turbulent diffusion on small scales. Here, we relate the scales of streaming and Kelvin-Helmholtz instability, which both regulate particle densities on the scales of gravitational collapse, directly to planetesimal formation. We support our analytic findings by performing 3D hydrodynamical simulations of streaming and Kelvin-Helmholtz instability and planetesimal formation. We find that the vertical extent of the particle mid-plane layer and the radial width of streaming instability filaments are set by the same characteristic length scale, thus governing the strength of turbulent diffusion on the scales of planetesimal formation. We present and successfully test a collapse criterion: $0.1 Q \beta \epsilon^{-1}Z^{-1} \lesssim 1$ and show that even for Solar metallicities, planetesimals can form in dead zones of sufficiently massive disks. For a given gas Toomre-parameter $Q$, pressure gradient $\beta$, metallicity $Z$ and local particle enhancement $\epsilon$, the collapse criterion also provides a range of unstable scales, instituting a promising path for studying initial planetesimal mass distributions. Streaming instability is not required for planetesimal collapse, but by increasing $\epsilon$, can evolve a system to instability.

Due to the gas rich environments of early circumstellar disks, the gravitational collapse of cool, dense regions of the disk form fragments largely composed of gas. During formation, disk fragments may attain increased metallicities as they interact with the surrounding disk material, whether through particle migration to pressure maxima or through mutual gravitational interaction. In this paper, we investigate the ability of fragments to collect and retain a significant solid component through gas-particle interactions in high-resolution 3D self-gravitating shearing box simulations. The formation of axissymmetric perturbations associated with gravitational instabilities allows particles of intermediate sizes to concentrate through aerodynamic drag forces. By the onset of fragmentation, the mass of local particle concentrations within the fragment are comparable to that of the gas component and the sebsequent gravitational collapse results in the formation of a solid core. We find that these cores can be up to several tens of Earth masses, depending on grain size, before the fragment center reaches temperatures which would sublimate solids. The solid fraction and total mass of the fragment also depend on the metallicity of the young parent protoplanetary disk, with higher initial metallicities resulting in larger fragments and larger solid cores. Additionally, the extended atmospheres of these soon-to-be gas giants or brown dwarfs are occasionally enriched above the initial metallicity, provided no solid core forms in the center and are otherwise lacking in heavier elements when a core does form.

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$.

Protoplanetary disks fragment due to gravitational instability when there is enough mass for self-gravitation, described by the Toomre parameter, and when heat can be lost at a rate comparable to the local dynamical timescale, described by t_c=beta Omega^-1. Simulations of self-gravitating disks show that the cooling parameter has a rough critical value at beta_crit=3. When below beta_crit, gas overdensities will contract under their own gravity and fragment into bound objects while otherwise maintaining a steady state of gravitoturbulence. However, previous studies of the critical cooling parameter have found dependence on simulation resolution, indicating that the simulation of self-gravitating protoplanetary disks is not so straightforward. In particular, the simplicity of the cooling timescale t_c prevents fragments from being disrupted by pressure support as temperatures rise. We alter the cooling law so that the cooling timescale is dependent on local surface density fluctuations, a means of incorporating optical depth effects into the local cooling of an object. For lower resolution simulations, this results in a lower critical cooling parameter and a disk more stable to gravitational stresses suggesting the formation of large gas giants planets in large, cool disks is generally suppressed by more realistic cooling. At our highest resolution however, the model becomes unstable to fragmentation for cooling timescales up to beta = 10.