The spectrum of stars in the submillimeter to centimeter wavelength range remains poorly constrained due to a lack of data for most spectral types. An accurate characterization of stellar emission in this regime is needed to test stellar atmosphere models, and is also essential for revealing emission associated with unresolved circumstellar debris. We present ALMA observations of the three nearby, main-sequence, debris-poor, F-type stars $\gamma$ Lep, $\gamma$ Vir A, and $\gamma$ Vir B at 0.87 and 1.29 millimeters. We use these data to constrain semi-empirical atmospheric models. We discuss the atmospheric structure of these stars, explore potential short term variability, and the potential impact on debris disk studies. These results are part of an ongoing campaign to obtain long wavelength observations of debris-poor stars, entitled Measuring the Emission of Stellar Atmospheres at Submillimeter/millimeter wavelengths (MESAS).

A planet's orbital orientation relative to an observer's line of sight determines the chord length for a transiting planet, i.e., the projected distance a transiting planet travels across the stellar disc. For a given circular orbit, the chord length determines the transit duration. Changes in the orbital inclination, the direction of the ascending node, or both, can alter this chord length and thus result in transit duration variations (TDVs). Variation of the full orbital inclination vector can even lead to de-transiting or newly transiting planets for a system. We use Laplace-Lagrange secular theory to estimate the fastest nodal eigenfrequencies for over 100 short-period planetary systems. The highest eigenfrequency is an indicator of which systems should show the strongest TDVs. We further explore five cases (TRAPPIST-1, Kepler-11, K2-138, Kepler-445, and Kepler-334) using direct N-body simulations to characterize possible TDVs and to explore whether de-transiting planets could be possible for these systems. A range of initial conditions are explored, with each realization being consistent with the observed transits. We find that tens of percent of multiplanet systems have fast enough eigenfrequencies to expect large TDVs on decade timescales. Among the directly integrated cases, we find that de-transiting planets could occur on decade timescales and TDVs of 10 minutes per decade should be common.

[Abridged] Asteroid mining is not necessarily a distant prospect. Hayabusa2 and OSIRIS-REx have recently rendezvoused with near-Earth asteroids and will return samples to Earth. While there is significant science motivation for these missions, there are also resource interests. Space agencies and commercial entities are particularly interested in ices and water-bearing minerals that could be used to produce rocket fuel in space. The internationally coordinated roadmaps of major space agencies depend on utilizing the natural resources of such celestial bodies. Several companies have already created plans for intercepting and extracting water and minerals from near-Earth objects, as even a small asteroid could have high economic worth. However, the low surface gravity of asteroids could make the release of mining waste and the subsequent formation of debris streams a consequence of asteroid mining. Strategies to contain material during extraction could still eventually require the purposeful jettison of waste to avoid managing unwanted mass. Using simulations, we explore the formation of mining debris streams by integrating particles released from four select asteroids. Radiation effects are included, and a range of debris sizes are explored. The simulation results are used to investigate the timescales for debris stream formation, the sizes of the streams, and the meteoroid fluxes compared with sporadic meteoroids. We find that for prodigious mining activities resulting in the loss of a few percent of the asteroid's mass or more, it is possible to produce streams that exceed the sporadic flux during stream crossing for some meteoroid sizes. The result of these simulations are intended to highlight potential unintended consequences that could result from NewSpace activity, which could help to inform efforts to develop international space resource guidelines.

This white paper submitted for the 2020 Canadian Long-Range Planning process (LRP2020) presents the prospects for Canada and the Square Kilometre Array (SKA) from 2020-2030, focussing on the first phase of the project (SKA1) scheduled to begin construction early in the next decade. SKA1 will make transformational advances in our understanding of the Universe across a wide range of fields, and Canadians are poised to play leadership roles in several. Canadian key SKA technologies will ensure a good return on capital investment in addition to strong scientific returns, positioning Canadian astronomy for future opportunities well beyond 2030. We therefore advocate for Canada's continued scientific and technological engagement in the SKA from 2020-2030 through participation in the construction and operations phases of SKA1.

There is a vibrant and effective planetary science community in Canada. We do research in the areas of meteoritics, asteroid and trans-Neptunian object orbits and compositions, and space weather, and are involved in space probe missions to study planetary surfaces and interiors. For Canadian planetary scientists to deliver the highest scientific impact possible, we have several recommendations. Our top recommendation is to join LSST and gain access to the full data releases by hosting a data centre, which could be done by adding to the CADC, which is already highly involved in hosting planetary data and supporting computational modelling for orbital studies. We also support MSE, which can provide spectroscopy and thus compositional information for thousands of small bodies. We support a Canadian-led microsatellite, POEP, which will provide small body sizes by measuring occultations. We support the idea of piggybacking space weather instruments on other astronomical space probes to provide data for the space weather community. Many Canadian planetary scientists are involved in space probe missions, but through haphazard and temporary arrangements like co-appointments at US institutions, so we would like the community to support Canadian researchers to participate in these large, international missions.

In this white paper (WP), we highlight several examples of small and moderate aperture telescopes that are being used for education and/or research. We further discuss potential costs for establishing new, small observatories, as well as joining existing international consortia. The WP includes a brief overview of select observing sites, with a discussion on how small telescopes at exceptional observing locations can be competitive, under certain circumstances, with larger and more expensive facilities located at poorer sites. In addition to research, these facilities enable many different types of educational experiences for wide range of people, from high school students to undergraduates to graduate students to postdocs. Canada should remain committed to partnering with large, international observatories such as CFHT, Gemini, and TMT, but it should also negotiate international agreements and commit funding to expand the use of small and moderate research observatories at domestic and international sites through coordination with the NRC, the Tri-Council, and the Canadian Foundation for Innovation. Both capital and operational costs (with site rental costs allowed) need to be included in support possibilities. CASCA should establish and maintain a small to moderate telescope expression of interest database that would help to facilitate Canadian institutions in organizing consortia, particularly for smaller institutions. The astronomical community should work with the NRC to make existing facilities more accessible to the astronomical community for research. This could involve, for example, automating the Plaskett and/or providing travel funds for supporting classical observing modes. Finally, a small to moderate aperture facility in the Arctic would be a world-class observatory and should be advanced over the next decade.

[Highly abridged, from executive summary] As much as NewSpace presents opportunities, there are significant challenges that must be overcome, requiring engagement with policy makers to influence domestic and international space governance. Failure to do so could result in a range of long-lasting negative outcomes for science and space stewardship. How will the Canadian astronomical community engage with NewSpace? What are the implications for NewSpace on the astro-environment, including Earth orbits, lunar and cis-lunar orbits, and surfaces of celestial bodies? This white paper analyzes the rapid changes in space use and what those changes could mean for Canadian astronomers. Our recommendations are as follows: Greater cooperation between the astronomical and the Space Situational Awareness communities is needed. Build closer ties between the astronomical community and Global Affairs Canada (GAC). Establish a committee for evaluating the astro-environmental impacts of human space use, including on and around the Moon and other bodies. CASCA and the Tri-Council should coordinate to identify programs that would enable Canadian astronomers to participate in pay-for-use services at appropriate funding levels. CASCA should continue to foster a relationship with CSA, but also build close ties to the private space industry. Canadian-led deep space missions are within Canada's capabilities, and should be pursued.

Successful exoplanet surveys in the last decade have revealed that planets are ubiquitous throughout the Milky Way, and show a large diversity in mass, location and composition. At the same time, new facilities such as the Atacama Large Millimeter/submillimeter Array (ALMA) and optical/infrared facilities including Gemini/GPI have provided us with sharper images than ever before of protoplanetary disks around young stars, the birth cradles of planets. The high spatial resolution has revealed astonishing structures in disks, such as rings, gaps, asymmetries and spiral arms, and the enormous jump in sensitivity has provided the tools for both large, statistically relevant surveys and deep, sensitive molecular line studies. These observations have revolutionized our view of planet formation, disk formation and disk evolution, bringing model simulations and observations closer to the same level of detail, with many contributions from Canadian researchers on theoretical, observational and technological sides. The new results have inevitably led to a range of new questions, which require next generation instruments such as the Next Generation Very Large Array (ngVLA) and large scale optical infrared facilities. In this white paper we will discuss the current transformation in our understanding of planet formation and the next steps and challenges in connecting theory with exoplanet demographics and protoplanetary disk observations for Canadian research.

TRAPPIST-1 is an ultracool dwarf (UCD) with a system of 7 terrestrial planets, at least three of which orbit in the habitable zone. The radio emission of such low-mass stars is poorly understood; few UCDs have been detected at radio frequencies at all, and the likelihood of detection is only loosely correlated with stellar properties. Relative to other low-mass stars, UCDs with slow rotation such as TRAPPIST-1 tend to be radio dim, whereas rapidly rotating UCDs tend to have strong radio emission - although this is not always the case. We present radio observations of TRAPPIST-1 using ALMA at 97.5 GHz and the VLA at 44 GHz. TRAPPIST-1 was not detected at either frequency and we place $3 \sigma$ upper flux limits of 10.6 and 16.2 $\mu$Jy, respectively. We use our results to constrain the magnetic properties and possible outgoing high energy particle radiation from the star. The presence of radio emission from UCDs is indicative of a stellar environment that could pose a threat to life on surrounding planets. Gyrosynchrotron emission, discernible at frequencies between 20 and 100 GHz, is one of the only processes that can be used to infer the presence of high energy particles released during magnetic reconnection events. Since M dwarfs are frequent hosts of terrestrial planets, characterizing their stellar emission is a crucial part of assessing habitability. Exposure to outgoing high energy particle radiation - traceable by radio flux - can erode planetary atmospheres. While our results do not imply that the TRAPPIST-1 planets are suitable for life, we find no evidence that they are overtly unsuitable due to proton fluxes.

Modeling the submillimeter to centimeter emission of stars is challenging due to a lack of sensitive observations at these long wavelengths. We launched an ongoing campaign to obtain new observations entitled Measuring the Emission of Stellar Atmospheres at Submillimeter/millimeter wavelengths (MESAS). Here we present ALMA, GBT, and VLA observations of Sirius A, the closest main-sequence A-type star, that span from 1.4 to 9.0 millimeters. These observations complement our previous millimeter data on Sirius A and are entirely consistent with the PHOENIX stellar atmosphere models constructed to explain them. We note that accurate models of long wavelength emission from stars are essential not only to understand fundamental stellar processes, but also to determine the presence of dusty debris in spatially unresolved observations of circumstellar disks.

Our understanding of stellar atmospheres and our ability to infer architectures of extrasolar planetary systems rely on understanding the emission of stars at submillimeter to centimeter wavelengths. In this chapter we describe how unconstrained stellar emission can interfere with the accurate characterization of circumstellar debris. The ngVLA is the only facility with the sensitivity that allows for the observations of a broad range of stellar spectral types in a feasible amount of time. The observations will enable the building and testing of accurate models of stellar emission, which in turn are required for evaluating both the occurrence and abundance of debris over the proposed wavelength range of the ngVLA

At the time of formation, protoplanetary discs likely contain a comparable mass to their host protostars. As a result, gravitational instabilities (GIs) are expected to play a pivotal role in the early phases of disc evolution. However, as these young objects are heavily embedded, confirmation of GIs has remained elusive. Therefore, we use the radiative transfer code LIME to produce line images of a $0.17\,\mathrm{M}_{\odot}$ self-gravitating protosolar-like disc. We note the limitations of using LIME and explore methods to improve upon the default gridding. We synthesise noiseless observations to determine the sensitivities required to detect the spiral flux, and find that the line flux distribution does not necessarily correlate to the abundance density distribution; hence performing radiative transfer calculations is imperative. Moreover, the spiral features are seen in absorption, due to the GI-heated midplane and high extinction, which could be indicative of GI activity. If a small beamsize and appropriate molecular line are used then spatially resolving spirals in a protosolar-like disc should be possible with ALMA for an on-source time of 30 hr. Spectrally resolving non-axisymmetric structure takes only a tenth as long for a reasonable noise level, but attributing this structure to GI-induced activity would be tentative. Finally, we find that identifying finger-like features in PV diagrams of nearly edge-on discs, which are a direct indicator of spirals, is feasible with an on-source time of 19 hr, and hence likely offers the most promising means of confirming GI-driven spiral structure in young, embedded protoplanetary discs.

In the early stages of planet formation, small dust grains grow to become mm sized particles in debris disks around stars. These disks can in principle be characterized by their emission at submillimeter and millimeter wavelengths. Determining both the occurrence and abundance of debris in unresolved circumstellar disks of A-type main-sequence stars requires that the stellar photospheric emission be accurately modeled. To better constrain the photospheric emission for such systems, we present observations of Sirius A, an A-type star with no known debris, from the JCMT, SMA, and VLA at 0.45, 0.85, 0.88, 1.3, 6.7, and 9.0 mm. We use these observations to inform a PHOENIX model of Sirius A's atmosphere. We find the model provides a good match to these data and can be used as a template for the submm/mm emission of other early A-type stars where unresolved debris may be present. The observations are part of an ongoing observational campaign entitled Measuring the Emission of Stellar Atmospheres at Submm/mm wavelengths (MESAS)

We present archival ALMA observations of the HD 141569 circumstellar disk at 345, 230, and 100 GHz. These data detect extended millimeter emission that is exterior to the inner disk. We find through simultaneous visibility modeling of all three data sets that the system's morphology is described well by a two-component disk model. The inner disk ranges from approximately 16 to 45 au with a spectral index of 1.81 (q = 2.95) and the outer disk ranges from 95 to 300 au with a spectral index of 2.28 (q = 3.21). Azimuthally averaged radial emission profiles derived from the continuum images at each frequency show potential emission that is consistent with the visibility modeling. The analysis presented here shows that at ~5 Myr HD 141569's grain size distribution is steeper, and therefore more evolved, in the outer disk than in the inner disk.

Recent laboratory efforts (Fu et al., 2014) have constrained the remanent magnetizations of chondrules and the magnetic field strengths at which the chondrules were exposed to as they cooled below their Curie points. An outstanding question is whether the inferred paleofields represent the background magnetic field of the solar nebula or were unique to the chondrule-forming environment. We investigate the amplification of the magnetic field above background values for two proposed chondrule formation mechanisms, large-scale nebular shocks and planetary bow shocks. Behind large-scale shocks, the magnetic field parallel to the shock front is amplified by factors $\sim 10-30$, regardless of the magnetic diffusivity. Therefore, chondrules melted in these shocks probably recorded an amplified magnetic field. Behind planetary bow shocks, the field amplification is sensitive to the magnetic diffusivity. We compute the gas properties behind a bow shock around a 3000 km-radius planetary embryo, with and without atmospheres, using hydrodynamics models. We calculate the ionization state of the hot, shocked gas, including thermionic emission from dust, and thermal ionization of gas-phase potassium atoms, and the magnetic diffusivity due to Ohmic dissipation and ambipolar diffusion. We find that the diffusivity is sufficiently large that magnetic fields have already relaxed to background values in the shock downstream where chondrules acquire magnetizations, and that these locations are sufficiently far from the planetary embryos that chondrules should not have recorded a significant putative dynamo field generated on these bodies. We conclude that, if melted in planetary bow shocks, chondrules probably recorded the background nebular field.

We explore the effects of an undetected outer giant planet on the dynamics, observability, and stability of Systems with Tightly-packed Inner Planets (STIPs). We use direct numerical simulations along with secular theory and synthetic secular frequency spectra to analyze how analogues of Kepler-11 and Kepler-90 behave in the presence of a nearly co-planar, Jupiter-like outer perturber with semi-major axes between 1 and 5.2 au. Most locations of the outer perturber do not affect the evolution of the inner planetary systems, apart from altering precession frequencies. However, there are locations at which an outer planet causes system instability due to, in part, secular eccentricity resonances. In Kepler-90, there is a range of orbital distances for which the outer perturber drives planets b and c, through secular interactions, onto orbits with inclinations that are $\sim16^\circ$ away from the rest of the planets. Kepler-90 is stable in this configuration. Such secular resonances can thus affect the observed multiplicity of transiting systems. We also compare the synthetic apsidal and nodal precession frequencies with the secular theory and find some misalignment between principal frequencies, indicative of strong interactions between the planets (consistent with the system showing TTVs). First-order libration angles are calculated to identify MMRs in the systems, for which two near-MMRs are shown in Kepler-90, with a 5:4 between b and c, as well as a 3:2 between g and h.

We present VLA 9 mm (33 GHz) observations of the HD 141569 system from semester 16A. The observations achieve a resolution of 0.25 arcsec ($\sim28$ au) and a sensitivity of $4.7~\mu \rm Jy~beam^{-1}$. We find (1) a $52\pm 5~\mu$Jy point source at the location of HD 141569A that shows potential variability, (2) the detected flux is contained within the SED-inferred central clearing of the disc meaning the spectral index of the dust disc is steeper than previously inferred, and (3) the M dwarf companions are also detected and variable. Previous lower-resolution VLA observations (semester 14A) found a higher flux density, interpreted as solely dust emission. When combined with ALMA observations, the VLA 14A observations suggested the spectral index and grain size distribution of HD 141569's disc was shallow and an outlier among debris systems. Using archival ALMA observations of HD 141569 at 0.87 mm and 2.9 mm we find a dust spectral index of $\alpha_{\rm mm} = 1.81\pm 0.20$. The VLA 16A flux corresponds to a brightness temperature of $\sim5\times10^{6}$ K, suggesting strong non-disc emission is affecting the inferred grain properties. The VLA 16A flux density of the M2V companion HD 141569B is $149\pm9~\mu$Jy, corresponding to a brightness temperature of $\sim2\times10^{8}$ K and suggesting significant stellar variability when compared to the VLA14A observations, which are smaller by a factor of $\sim6$.

Recent observations of HL Tau revealed remarkably detailed structure within the system's circumstellar disc. A range of hypotheses have been proposed to explain the morphology, including, e.g., planet-disc interactions, condensation fronts, and secular gravitational instabilities. While embedded planets seem to be able to explain some of the major structure in the disc through interactions with gas and dust, the substructure, such as low-contrast rings and bands, are not so easily reproduced. Here, we show that dynamical interactions between three planets (only two of which are modelled) and an initial population of large planetesimals can potentially explain both the major and minor banded features within the system. In this context, the small grains, which are coupled to the gas and reveal the disc morphology, are produced by the collisional evolution of the newly-formed planetesimals, which are ubiquitous in the system and are decoupled from the gas.

In this paper, we model the chemical evolution of a 0.25 M$_{\odot}$ protoplanetary disc surrounding a 1 M$_{\odot}$ star that undergoes fragmentation due to self-gravity. We use Smoothed Particle Hydrodynamics including a radiative transfer scheme, along with time-dependent chemical evolution code to follow the composition of the disc and resulting fragments over approximately 4000 yrs. Initially, four quasi-stable fragments are formed, of which two are eventually disrupted by tidal torques in the disc. From the results of our chemical modelling, we identify species that are abundant in the fragments (e.g. H$_{\rm 2}$O, H$_{\rm 2}$S, HNO, N$_{\rm 2}$, NH$_{\rm 3}$, OCS, SO), species that are abundant in the spiral shocks within the disc (e.g. CO, CH$_{\rm 4}$, CN, CS, H$_{\rm 2}$CO), and species which are abundant in the circumfragmentary material (e.g. HCO$^{\rm +}$). Our models suggest that in some fragments it is plausible for grains to sediment to the core before releasing their volatiles into the planetary envelope, leading to changes in, e.g., the C/O ratio of the gas and ice components. We would therefore predict that the atmospheric composition of planets generated by gravitational instability should not necessarily follow the bulk chemical composition of the local disc material.

The growth and migration of planetesimals in a young protoplanetary disc are fundamental to planet formation. In all models of early growth, there are several processes that can inhibit grains from reaching larger sizes. Nevertheless, observations suggest that growth of planetesimals must be rapid. If a small number of 100 km sized planetesimals do manage to form in the disc, then gas drag effects could enable them to efficiently accrete small solids from beyond their gravitationally focused cross-section. This gas drag-enhanced accretion can allow planetesimals to grow at rapid rates, in principle. We present self-consistent hydrodynamics simulations with direct particle integration and gas drag coupling to estimate the rate of planetesimal growth due to pebble accretion. Wind tunnel simulations are used to explore a range of particle sizes and disc conditions. We also explore analytic estimates of planetesimal growth and numerically integrate planetesimal drift due to the accretion of small solids. Our results show that, for almost every case that we consider, there is a clearly preferred particle size for accretion that depends on the properties of the accreting planetesimal and the local disc conditions. For solids much smaller than the preferred particle size, accretion rates are significantly reduced as the particles are entrained in the gas and flow around the planetesimal. Solids much larger than the preferred size accrete at rates consistent with gravitational focusing. Our analytic estimates for pebble accretion highlight the timescales that are needed for the growth of large objects under different disc conditions and initial planetesimal sizes.

Gravitational instabilities (GIs) are most likely a fundamental process during the early stages of protoplanetary disc formation. Recently, there have been detections of spiral features in young, embedded objects that appear consistent with GI-driven structure. It is crucial to perform hydrodynamic and radiative transfer simulations of gravitationally unstable discs in order to assess the validity of GIs in such objects, and constrain optimal targets for future observations. We utilise the radiative transfer code LIME to produce continuum emission maps of a $0.17\,\mathrm{M}_{\odot}$ self-gravitating protosolar-like disc. We note the limitations of using LIME as is and explore methods to improve upon the default gridding. We use CASA to produce synthetic observations of 270 continuum emission maps generated across different frequencies, inclinations and dust opacities. We find that the spiral structure of our protosolar-like disc model is distinguishable across the majority of our parameter space after 1 hour of observation, and is especially prominent at 230$\,$GHz due to the favourable combination of angular resolution and sensitivity. Disc mass derived from the observations is sensitive to the assumed dust opacities and temperatures, and therefore can be underestimated by a factor of at least 30 at 850$\,$GHz and 2.5 at 90$\,$GHz. As a result, this effect could retrospectively validate GIs in discs previously thought not massive enough to be gravitationally unstable, which could have a significant impact on the understanding of the formation and evolution of protoplanetary discs.

We present ALMA Band 6 observations (1.3 mm/233 GHz) of Fomalhaut and its debris disc. The observations achieve a sensitivity of 17 $\mu$Jy and a resolution of 0.28 arcsec (2.1 au at a distance of 7.66 pc), which are the highest resolution observations to date of the millimetre grains in Fomalhaut's main debris ring. The ring is tightly constrained to $139^{+2}_{-3}$ au with a FWHM of $13\pm3$ au, following a Gaussian profile. The millimetre spectral index is constrained to $\alpha_{mm} = -2.62\pm0.12$. We explore fitting debris disc models in the image plane, as well as fitting models using visibility data directly. The results are compared and the potential advantages/disadvantages of each approach are discussed. The detected central emission is indistinguishable from a point source, with a most probable flux of $0.90\pm 0.12$ mJy (including calibration uncertainties). This implies that any inner debris structure, as was inferred from far-Infrared observations, must contribute little to the total central emission. Moreover, the stellar flux is less than 70\% of that predicted by extrapolating a black body from the constrained stellar photosphere temperature. This result emphasizes that unresolved inner debris components cannot be fully characterized until the behaviour of the host star's intrinsic stellar emission at millimetre wavelengths is properly understood.

We present ALMA band 7 (345 GHz) continuum and $^{12}$CO(J = 3-2) observations of the circumstellar disk surrounding HD141569. At an age of about 5 Myr, the disk has a complex morphology that may be best interpreted as a nascent debris system with gas. Our $870\rm~\mu m$ ALMA continuum observations resolve a dust disk out to approximately $ 56 ~\rm au$ from the star (assuming a distance of 116 pc) with $0."38$ resolution and $0.07 ~ \rm mJy~beam^{-1}$ sensitivity. We measure a continuum flux density for this inner material of $3.8 \pm 0.4 ~ \rm mJy$ (including calibration uncertainties). The $^{12}$CO(3-2) gas is resolved kinematically and spatially from about 30 to 210 au. The integrated $^{12}$CO(3-2) line flux density is $15.7 \pm 1.6~\rm Jy~km~s^{-1}$. We estimate the mass of the millimeter debris and $^{12}$CO(3-2) gas to be $\gtrsim0.04~\rm M_{\oplus}$ and $\sim2\times 10^{-3}~\rm M_{\oplus}$, respectively. If the millimeter grains are part of a collisional cascade, then we infer that the inner disk ($<50$ au) has $\sim 160~\rm M_{\oplus}$ contained within objects less than 50 km in radius, depending on the planetesimal size distribution and density assumptions. MCMC modeling of the system reveals a disk morphology with an inclination of $53.4^{\circ}$ centered around a $\rm M=2.39~ M_{\odot}$ host star ($\rm Msin(i)=1.92~ M_{\odot}$). We discuss whether the gas in HD141569's disk may be second generation. If it is, the system can be used to study the clearing stages of planet formation.

We use radiation hydrodynamics with direct particle integration to explore the feasibility of chondrule formation in planetary embryo bow shocks. The calculations presented here are used to explore the consequences of a Mars-size planetary embryo traveling on a moderately excited orbit through the dusty, early environment of the solar system. The embryo's eccentric orbit produces a range of supersonic relative velocities between the embryo and the circularly orbiting gas and dust, prompting the formation of bow shocks. Temporary atmospheres around these embryos, which can be created via volatile outgassing and gas capture from the surrounding nebula, can non-trivially affect thermal profiles of solids entering the shock. We explore the thermal environment of solids that traverse the bow shock at different impact radii, the effects that planetoid atmospheres have on shock morphologies, and the stripping efficiency of planetoidal atmospheres in the presence of high relative winds. Simulations are run using adiabatic and radiative conditions, with multiple treatments for the local opacities. Shock speeds of 5, 6, and 7 km/s are explored. We find that a high-mass atmosphere and inefficient radiative conditions can produce peak temperatures and cooling rates that are consistent with the constraints set by chondrule furnace studies. For most conditions, the derived cooling rates are potentially too high to be consistent with chondrule formation.

The disk around HD 141569 is one of a handful of systems whose weak infrared emission is consistent with a debris disk, but still has a significant reservoir of gas. Here we report spatially resolved mm observations of the CO(3-2) and CO(1-0) emission as seen with the SMA and CARMA. We find that the excitation temperature for CO is lower than expected from cospatial blackbody grains, similar to previous observations of analogous systems, and derive a gas mass that lies between that of gas-rich primordial disks and gas-poor debris disks. The data also indicate a large inner hole in the CO gas distribution and an outer radius that lies interior to the outer scattered light rings. This spatial distribution, with the dust rings just outside the gaseous disk, is consistent with the expected interactions between gas and dust in an optically thin disk. This indicates that gas can have a significant effect on the location of the dust within debris disks.

We propose that two of the most surprising results so far among exoplanet discoveries are related: the existences of both hot Jupiters and the high frequency of multi-planet systems with periods $P\lesssim200$~days. In this paradigm, the vast majority of stars rapidly form along with multiple close-in planets in the mass range of Mars to super-Earths/mini-Neptunes. Such systems of tightly-packed inner planets (STIPs) are metastable, with the time scale of the dynamical instability having a major influence on final planet types. In most cases, the planets consolidate into a system of fewer, more massive planets, but long after the circumstellar gas disk has dissipated. This can yield planets with masses above the traditional critical core of $\sim$10 $M_\oplus$, yielding short-period giants that lack abundant gas. A rich variety of physical states are also possible given the range of collisional outcomes and formation time of the close-in planets. However, when dynamical consolidation occurs before gas dispersal, a critical core can form that then grows via gas capture into a short-period gas giant. In this picture the majority of Hot and Warm Jupiters formed locally, rather than migrating down from larger distances.

Recent high-resolution, near-infrared images of protoplanetary disks have shown that these disks often present spiral features. Spiral arms are among the structures predicted decades ago by numerical simulations of disk-planet interaction and thus it is tempting to suspect that planetary perturbers are responsible for the observed signatures. However, such interpretation is not free of problems. The spirals are found to have large pitch angles, and in at least one case (HD 100546) the spiral feature appears effectively unpolarized, which implies thermal emission of the order of 1000K (465$\pm$40K at closer inspection). We have recently shown in two-dimensional models that shock dissipation in the supersonic wake of high-mass planets can lead to significant heating if the disk is sufficiently adiabatic. In this paper we extend this analysis to three dimensions in thermodynamically evolving disks. We use the Pencil Code in spherical coordinates for our models, with a prescription for thermal cooling based on the optical depth of the local vertical gas column. We use a 5$M_J$ planet, and show that shocks in the region around the planet where the Lindblad resonances occur heat the gas to substantially higher temperatures than the ambient disk gas at that radius. The gas is accelerated vertically away from the midplane by the shocks to form shock bores, and the gas falling back toward the midplane breaks up into a turbulent surf near the Lindblad resonances. This turbulence, although localized, has high $\alpha$ values, reaching 0.05 in the inner Lindblad resonance, and 0.1 in the outer one. We also find evidence that the disk regions heated up by the planetary shocks eventually becomes superadiabatic, generating convection far from the planet's orbit.

To date, most simulations of the chemistry in protoplanetary discs have used 1+1D or 2D axisymmetric $\alpha$-disc models to determine chemical compositions within young systems. This assumption is inappropriate for non-axisymmetric, gravitationally unstable discs, which may be a significant stage in early protoplanetary disc evolution. Using 3D radiative hydrodynamics, we have modelled the physical and chemical evolution of a 0.17 M$_{\odot}$ self-gravitating disc over a period of 2000 yr. The 0.8 M$_{\odot}$ central protostar is likely to evolve into a solar-like star, and hence this Class 0 or early Class I young stellar object may be analogous to our early Solar System. Shocks driven by gravitational instabilities enhance the desorption rates, which dominate the changes in gas-phase fractional abundances for most species. We find that at the end of the simulation, a number of species distinctly trace the spiral structure of our relatively low-mass disc, particularly CN. We compare our simulation to that of a more massive disc, and conclude that mass differences between gravitationally unstable discs may not have a strong impact on the chemical composition. We find that over the duration of our simulation, successive shock heating has a permanent effect on the abundances of HNO, CN and NH$_3$, which may have significant implications for both simulations and observations. We also find that HCO$^+$ may be a useful tracer of disc mass. We conclude that gravitational instabilities induced in lower mass discs can significantly, and permanently, affect the chemical evolution, and that observations with high-resolution instruments such as ALMA offer a promising means of characterising gravitational instabilities in protosolar discs.

Recent observations of gaps and non-axisymmetric features in the dust distributions of transition disks have been interpreted as evidence of embedded massive protoplanets. However, comparing the predictions of planet-disk interaction models to the observed features has shown far from perfect agreement. This may be due to the strong approximations used for the predictions. For example, spiral arm fitting typically uses results that are based on low-mass planets in an isothermal gas. In this work, we describe two-dimensional, global, hydrodynamical simulations of disks with embedded protoplanets, with and without the assumption of local isothermality, for a range of planet-to-star mass ratios 1-10 M_jup for a 1 M_sun star. We use the Pencil Code in polar coordinates for our models. We find that the inner and outer spiral wakes of massive protoplanets (M>5 M_jup) produce significant shock heating that can trigger buoyant instabilities. These drive sustained turbulence throughout the disk when they occur. The strength of this effect depends strongly on the mass of the planet and the thermal relaxation timescale; for a 10 M_jup planet embedded in a thin, purely adiabatic disk, the spirals, gaps, and vortices typically associated with planet-disk interactions are disrupted. We find that the effect is only weakly dependent on the initial radial temperature profile. The spirals that form in disks heated by the effects we have described may fit the spiral structures observed in transition disks better than the spirals predicted by linear isothermal theory.

We present a solution to the long outstanding meter barrier problem in planet formation theory. As solids spiral inward due to aerodynamic drag, they will enter disk regions that are characterized by high temperatures, densities, and pressures. High partial pressures of rock vapor can suppress solid evaporation, and promote collisions between partially molten solids, allowing rapid growth. This process should be ubiquitous in planet-forming disks, which may be evidenced by the abundant class of Systems with Tightly-packed Inner Planets (STIPs) discovered by the NASA Kepler mission.

Giant planet formation in the core accretion (CA) paradigm is predicated by the formation of a core, assembled by the coagulation of grains and later by planetesimals within a protoplanetary disc. In contrast, in the disc instability paradigm, giant planet formation is believed to be independent of core formation: massive self-gravitating gas fragments cool radiatively and collapse as a whole. We show that giant planet formation in the disc instability model may be also enhanced by core formation for reasons physically very similar to the CA paradigm. In the model explored here, efficient grain sedimentation within an initial fragment (rather than the disc) leads to the formation of a core composed of heavy elements. We find that massive atmospheres form around cores and undergo collapse as a critical core mass is exceeded, analogous to CA theory. The critical mass of the core to initiate such a collapse depends on the fragment mass and metallicity, as well as core luminosity, but ranges from less than 1 to as much as $\sim80$ Earth masses. We therefore suggest that there are two channels for the collapse of a gaseous fragment to planetary scales within the disc instability model: (i) H$_2$ dissociative collapse of the entire gaseous clump, and (ii) core-assisted gas capture, as presented here. We suggest that the first of these two is favoured in metal-poor environments and for fragments more massive than $\sim 5-10$ Jupiter masses, whereas the second is favored in metal-rich environments and fragments of lower mass. [Abridged]

The large number of detected giant exoplanets offers the opportunity to improve our understanding of the formation mechanism, evolution, and interior structure of gas giant planets. The two main models for giant planet formation are core accretion and disk instability. There are substantial differences between these formation models, including formation timescale, favorable formation location, ideal disk properties for planetary formation, early evolution, planetary composition, etc. First, we summarize the two models including their substantial differences, advantages, and disadvantages, and suggest how theoretical models should be connected to available (and future) data. We next summarize current knowledge of the internal structures of solar- and extrasolar- giant planets. Finally, we suggest the next steps to be taken in giant planet exploration.

A fundamental, unsolved problem in Solar System formation is explaining the melting and crystallization of chondrules found in chondritic meteorites. Theoretical models of chondrule melting in nebular shocks has been shown to be consistent with many aspects of thermal histories inferred for chondrules from laboratory experiments; but, the mechanism driving these shocks is unknown. Planetesimals and planetary embryos on eccentric orbits can produce bow shocks as they move supersonically through the disk gas, and are one possible source of chondrule-melting shocks. We investigate chondrule formation in bow shocks around planetoids through 3D radiation hydrodynamics simulations. A new radiation transport algorithm that combines elements of flux-limited diffusion and Monte Carlo methods is used to capture the complexity of radiative transport around bow shocks. An equation of state that includes the rotational, vibrational, and dissociation modes of H$_2$ is also used. Solids are followed directly in the simulations and their thermal histories are recorded. Adiabatic expansion creates rapid cooling of the gas, and tail shocks behind the embryo can cause secondary heating events. Radiative transport is efficient, and bow shocks around planetoids can have luminosities $\sim$few$\times10^{-8}$ L$_{\odot}$. While barred and radial chondrule textures could be produced in the radiative shocks explored here, porphyritic chondrules may only be possible in the adiabatic limit. We present a series of predicted cooling curves that merit investigation in laboratory experiments to determine whether the solids produced by bow shocks are represented in the meteoritic record by chondrules or other solids.

The NASA Kepler mission has revealed an abundant class of Systems with Tightly-packed Inner Planets (STIPs). The current paradigm for planet formation suggests that small planetesimals will quickly spiral into the host star due to aerodynamic drag, preventing rocky planet formation. In contrast, we find that aerodynamic drift, when acting on an ensemble of solids, can concentrate mass at short orbital periods in gaseous disks. Sublimation fronts may further aid this process. Kepler data suggest that the innermost known planets are found near the silicate sublimation zone. STIP planets should have a wide range of volatile fractions due to aerodynamic drift and H2 dissociation-driven gas accretion. We further propose that the low mass of Mars is evidence that the Solar System was once a proto-STIP.

The formation and earliest stages of protoplanetary discs remain poorly constrained by observations. ALMA will soon revolutionise this field. Therefore, it is important to provide predictions which will be valuable for the interpretation of future high sensitivity and high angular resolution observations. Here we present simulated ALMA observations based on radiative transfer modelling of a relatively massive (0.39 M_solar) self-gravitating disc embedded in a 10 M_solar dense core, with structure similar to the pre-stellar core L1544. We focus on simple species and conclude that C17O 3-2, HCO+ 3-2, OCS 26-25 and H2CO 404-303 lines can be used to probe the disc structure and kinematics at all scales.

Using radiation hydrodynamics simulations, we explore the evolution of circumplanetary disks around wide-orbit proto-gas giants. At large distances from the star (~100 AU), gravitational instability followed by disk fragmentation can form low-mass substellar companions (massive gas giants and/or brown dwarfs) that are likely to host large disks. We examine the initial evolution of these subdisks and their role in regulating the growth of their substellar companions, as well as explore consequences of their interactions with circumstellar material. We find that subdisks that form in the context of GIs evolve quickly from a very massive state. Long-term accretion rates from the subdisk onto the proto-gas giant reach ~0.3 Jupiter masses per kyr. We also find consistency with previous simulations, demonstrating that subdisks are truncated at ~1/3 of the companion's Hill radius and are thick, with (h/r) of great than or equal to 0.2. The thickness of subdisks draws to question the use of thin-disk approximations for understanding the behavior of subdisks, and the morphology of subdisks has implications for the formation and extent of satellite systems. These subdisks create heating events in otherwise cold regions of the circumstellar disk, and serve as planet formation beacons that can be detected by instruments such as ALMA.

The strength and direction of migration of low mass embedded planets depends on the disk's thermodynamic state, where the internal dissipation is balanced by radiative transport, and the migration can be directed outwards, a process which extends the lifetime of growing embryos. Very important parameters determining the structure of disks, and hence the direction of migration, are the viscosity and the adiabatic index. In this paper we investigate the influence of different viscosity prescriptions (alpha-type and constant) and adiabatic indices on disk structures and how this affects the migration rate of planets embedded in such disks. We perform 3D numerical simulations of accretion disks with embedded planets. We use the explicit/implicit hydrodynamical code NIRVANA that includes full tensor viscosity and radiation transport in the flux-limited diffusion approximation, as well as a proper equation of state for molecular hydrogen. The migration of embedded 20Earthmass planets is studied. Low-viscosity disks have cooler temperatures and the migration rates of embedded planets tend toward the isothermal limit. In these disks, planets migrate inwards even in the fully radiative case. The effect of outward migration can only be sustained if the viscosity in the disk is large. Overall, the differences between the treatments for the equation of state seem to play a more important role in disks with higher viscosity. A change in the adiabatic index and in the viscosity changes the zero-torque radius that separates inward from outward migration. For larger viscosities, temperatures in the disk become higher and the zero-torque radius moves to larger radii, allowing outward migration of a 20 Earth-mass planet to persist over an extended radial range. In combination with large disk masses, this may allow for an extended period of the outward migration of growing protoplanetary cores.

We present Swarm-NG, a C++ library for the efficient direct integration of many n-body systems using highly-parallel Graphics Processing Unit (GPU), such as NVIDIA's Tesla T10 and M2070 GPUs. While previous studies have demonstrated the benefit of GPUs for n-body simulations with thousands to millions of bodies, Swarm-NG focuses on many few-body systems, e.g., thousands of systems with 3...15 bodies each, as is typical for the study of planetary systems. Swarm-NG parallelizes the simulation, including both the numerical integration of the equations of motion and the evaluation of forces using NVIDIA's "Compute Unified Device Architecture" (CUDA) on the GPU. Swarm-NG includes optimized implementations of 4th order time-symmetrized Hermite integration and mixed variable symplectic integration, as well as several sample codes for other algorithms to illustrate how non-CUDA-savvy users may themselves introduce customized integrators into the Swarm-NG framework. To optimize performance, we analyze the effect of GPU-specific parameters on performance under double precision. Applications of Swarm-NG include studying the late stages of planet formation, testing the stability of planetary systems and evaluating the goodness-of-fit between many planetary system models and observations of extrasolar planet host stars (e.g., radial velocity, astrometry, transit timing). While Swarm-NG focuses on the parallel integration of many planetary systems,the underlying integrators could be applied to a wide variety of problems that require repeatedly integrating a set of ordinary differential equations many times using different initial conditions and/or parameter values.

We present 3D smoothed particle hydrodynamics simulations of the collapse of clumps formed through gravitational instability in the outer part of a protoplanetary disc. The initial conditions are taken directly from a global disc simulation, and a realistic equation of state is used to follow the clumps as they contract over several orders of magnitude in density, approaching the molecular hydrogen dissociation stage. The effects of clump rotation, asymmetries, and radiative cooling are studied. Rotation provides support against fast collapse, but non-axisymmetric modes develop and efficiently transport angular momentum outward, forming a circumplanetary disc. This transport helps the clump reach the dynamical collapse phase, resulting from molecular hydrogen dissociation, on a thousand-year timescale, which is smaller than timescales predicted by some previous spherical 1D collapse models. Extrapolation to the threshold of the runaway hydrogen dissociation indicates that the collapse timescales can be shorter than inward migration timescales, suggesting that clumps could survive tidal disruption and deliver a proto-gas giant to distances of even a few AU from the central star.

The chance that a planetary system will interact with another member of its host star's nascent cluster would be greatly increased if gas giant planets form in situ on wide orbits. In this paper, we explore the outcomes of planet-planet scattering for a distribution of multiplanet systems that all have one of the planets on an initial orbit of 100 AU. The scattering experiments are run with and without stellar flybys. We convolve the outcomes with distributions for protoplanetary disk and stellar cluster sizes to generalize the results where possible. We find that the frequencies of large mutual inclinations and high eccentricities are sensitive to the number of planets in a system, but not strongly to stellar flybys. However, flybys do play a role in changing the low and moderate portions of the mutual inclination distributions, and erase dynamically cold initial conditions on average. Wide-orbit planets can be mixed throughout the planetary system, and in some cases, can potentially become hot Jupiters, which we demonstrate using scattering experiments that include a tidal damping model. If planets form on wide orbits in situ, then there will be discernible differences in the proper motion distributions of a sample of wide-orbit planets compared with a pure scattering formation mechanism. Stellar flybys can enhance the frequency of ejections in planetary systems, but auto-ionization is likely to remain the dominant source of free-floating planets.

Recent isotopic studies of Martian meteorites by Dauphas & Pourmond (2011) have established that large (~ 3000 km radius) planetary embryos existed in the solar nebula at the same time that chondrules - millimeter-sized igneous inclusions found in meteorites - were forming. We model the formation of chondrules by passage through bow shocks around such a planetary embryo on an eccentric orbit. We numerically model the hydrodynamics of the flow, and find that such large bodies retain an atmosphere, with Kelvin-Helmholtz instabilities allowing mixing of this atmosphere with the gas and particles flowing past the embryo. We calculate the trajectories of chondrules flowing past the body, and find that they are not accreted by the protoplanet, but may instead flow through volatiles outgassed from the planet's magma ocean. In contrast, chondrules are accreted onto smaller planetesimals. We calculate the thermal histories of chondrules passing through the bow shock. We find that peak temperatures and cooling rates are consistent with the formation of the dominant, porphyritic texture of most chondrules, assuming a modest enhancement above the likely solar nebula average value of chondrule densities (by a factor of 10), attributable to settling of chondrule precursors to the midplane of the disk or turbulent concentration. We calculate the rate at which a planetary embryo's eccentricity is damped and conclude that a single planetary embryo scattered into an eccentric orbit can, over ~ 10e5 years, produce ~ 10e24 g of chondrules. In principle, a small number (1-10) of eccentric planetary embryos can melt the observed mass of chondrules in a manner consistent with all known constraints.

The dynamical evolution of planetary systems leaves observable signatures in debris disks. Optical images trace micron-sized grains, which are strongly affected by stellar radiation and need not coincide with their parent body population. Observations of mm-size grains accurately trace parent bodies, but previous images lack the resolution and sensitivity needed to characterize the ring's morphology. Here we present ALMA 350 GHz observations of the Fomalhaut debris ring. These observations demonstrate that the parent body population is 13-19 AU wide with a sharp inner and outer boundary. We discuss three possible origins for the ring, and suggest that debris confined by shepherd planets is the most consistent with the ring's morphology.

Until now, axisymmetric, alpha-disc models have been adopted for calculations of the chemical composition of protoplanetary discs. While this approach is reasonable for many discs, it is not appropriate when self-gravity is important. In this case, spiral waves and shocks cause temperature and density variations that affect the chemistry. We have adopted a dynamical model of a solar-mass star surrounded by a massive (0.39 Msun), self-gravitating disc, similar to those that may be found around Class 0 and early Class I protostars, in a study of disc chemistry. We find that for each of a number of species, e.g. H2O, adsorption and desorption dominate the changes in the gas-phase fractional abundance; because the desorption rates are very sensitive to temperature, maps of the emissions from such species should reveal the locations of shocks of varying strengths. The gas-phase fractional abundances of some other species, e.g. CS, are also affected by gas-phase reactions, particularly in warm shocked regions. We conclude that the dynamics of massive discs have a strong impact on how they appear when imaged in the emission lines of various molecular species.

We present high-resolution 3D smoothed particle hydrodynamics simulations of the formation and evolution of protostellar discs in a turbulent molecular cloud. Using a piecewise polytropic equation of state, we perform two sets of simulations. In both cases we find that isolated systems undergo a fundamentally different evolution than members of binary or multiple systems. When formed, isolated systems must accrete mass and increase their specific angular momentum, leading to the formation of massive, extended discs, which undergo strong gravitational instabilities and are susceptible to disc fragmentation. Fragments with initial masses of 5.5 M_jup, 7.4 M_jup and 12 M_jup are produced in our simulations. In binaries and small clusters, we observe that due to competition for material from the parent core, members do not accrete significant amounts of high specific angular momentum gas relative to isolated systems. We find that discs in multiple systems are strongly self-gravitating but that they are stable against fragmentation due to disc truncation and mass profile steeping by tides, accretion of high specific angular momentum gas by other members, and angular momentum being redirected into members' orbits. In general, we expect disc fragmentation to be less likely in clusters and to be more a feature of isolated systems.

Characterization of migration in gravitationally unstable disks is necessary to understand the fate of protoplanets formed by disk instability. As part of a larger study, we are using a 3D radiative hydrodynamics code to investigate how an embedded gas giant planet interacts with a gas disk that undergoes gravitational instabilities (GIs). This Letter presents results from simulations with a Jupiter-mass planet placed in orbit at 25 AU within a 0.14 $M_{\odot}$ disk. The disk spans 5 to 40 AU around a 1 $M_{\odot}$ star and is initially marginally unstable. In one simulation, the planet is inserted prior to the eruption of GIs; in another, it is inserted only after the disk has settled into a quasi-steady GI-active state, where heating by GIs roughly balances radiative cooling. When the planet is present from the beginning, its own wake stimulates growth of a particular global mode with which it strongly interacts, and the planet plunges inward six AU in about 10$^3$ years. In both cases with embedded planets, there are times when the planet's radial motion is slow and varies in direction. At other times, when the planet appears to be interacting with strong spiral modes, migration both inward and outward can be relatively rapid, covering several AUs over hundreds of years. Migration in both cases appears to stall near the inner Lindblad resonance of a dominant low-order mode. Planet orbit eccentricities fluctuate rapidly between about 0.02 to 0.1 throughout the GI-active phases of the simulations.

Transit surveys combined with Doppler data have revealed a class of gas giant planets that are massive and highly enriched in heavy elements (e.g., HD149026b, GJ436b, and HAT-P-20b). It is tempting to consider these planets as validation of core accretion plus gas capture because it is often assumed that disk instability planets should be of nebular composition. We show in this paper, to the contrary, that gas giants that form by disk instability can have a variety of heavy element compositions, ranging from sub- to super-nebular values. High levels of enrichment can be achieved through one or multiple mechanisms, including enrichment at birth, planetesimal capture, and differentiation plus tidal stripping. As a result, the metallicity of an individual gas giant cannot be used to discriminate between gas giant formation modes.

We investigate the coupling between rock-size solids and gas during the formation of gas giant planets by disk fragmentation in the outer regions of massive disks. In this study, we use three-dimensional radiative hydrodynamics simulations and model solids as a spatial distribution of particles. We assume that half of the total solid fraction is in small grains and half in large solids. The former are perfectly entrained with the gas and set the opacity in the disk, while the latter are allowed to respond to gas drag forces, with the back reaction on the gas taken into account. To explore the maximum effects of gas-solid interactions, we first consider 10cm-size particles. We then compare these results to a simulation with 1 km-size particles, which explores the low-drag regime. We show that (1) disk instability planets have the potential to form large cores due to aerodynamic capturing of rock-size solids in spiral arms before fragmentation; (2) that temporary clumps can concentrate tens of $M_{\oplus}$ of solids in very localized regions before clump disruption; (3) that the formation of permanent clumps, even in the outer disk, is dependent on the grain-size distribution, i.e., the opacity; (4) that nonaxisymmetric structure in the disk can create disk regions that have a solids-to-gas ratio greater than unity; (5) that the solid distribution may affect the fragmentation process; (6) that proto-gas giants and proto-brown dwarfs can start as differentiated objects prior to the H$_2$ collapse phase; (7) that spiral arms in a gravitationally unstable disk are able to stop the inward drift of rock-size solids, even redistributing them to larger radii; and, (8) that large solids can form spiral arms that are offset from the gaseous spiral arms. We conclude that planet embryo formation can be strongly affected by the growth of solids during the earliest stages of disk accretion.

We explore the initial conditions for fragments in the extended regions ($r\gtrsim 50$ AU) of gravitationally unstable disks. We combine analytic estimates for the fragmentation of spiral arms with 3D SPH simulations to show that initial fragment masses are in the gas giant regime. These initial fragments will have substantial angular momentum, and should form disks with radii of a few AU. We show that clumps will survive for multiple orbits before they undergo a second, rapid collapse due to H$_2$ dissociation and that it is possible to destroy bound clumps by transporting them into the inner disk. The consequences of disrupted clumps for planet formation, dust processing, and disk evolution are discussed. We argue that it is possible to produce Earth-mass cores in the outer disk during the earliest phases of disk evolution.

We place constraints on the formation redshifts for blue globular clusters (BGCs), independent of the details of hydrodynamics and population III star formation. The observed radial distribution of BGCs in the Milky Way Galaxy suggests that they formed in biased dark matter halos at high redshift. As a result, simulations of a ~1 Mpc box up to z~10 must resolve BGC formation in LCDM. We find that most halo stars could be produced from destroyed BGCs and other low-mass clusters that formed at high redshift. We present a proof-of-concept simulation that captures the formation of globular-like star clusters.

I argue for two modes of gas giant planet formation and discuss the conditions under which each mode operates. Gas giant planets at disk radii $r>100$ AU are likely to form in situ by disk instability, while core accretion plus gas capture remains the dominant formation mechanism for $r<100$ AU. During the mass accretion phase, mass loading can push disks toward fragmentation conditions at large $r$. Massive, extended disks can fragment into clumps of a few to tens of Jupiter masses. This is confirmed by radiation hydrodynamics simulations. The two modes of gas giant formation should lead to a bimodal distribution of gas giant semi-major axes. Because core accretion is expected to be less efficient in low-metallicity systems, the ratio of gas giants at large $r$ to planets at small $r$ should increase with decreasing metallicity.

Using analytic arguments and numerical simulations, we examine whether chondrule formation and the FU Orionis phenomenon can be caused by the burst-like onset of gravitational instabilities (GIs) in dead zones. At least two scenarios for bursting dead zones can work, in principle. If the disk is on the verge of fragmention, GI activation near $r\sim4$ to 5 AU can produce chondrule-forming shocks, at least under extreme conditions. Mass fluxes are also high enough during the onset of GIs to suggest that the outburst is related to an FU Orionis phenomenon. This situation is demonstrated by numerical simulations. In contrast, as supported by analytic arguments, if the burst takes place close to $r\sim1$ AU, then even low pitch angle spiral waves can create chondrule-producing shocks and outbursts. We also study the stability of the massive disks in our simulations against fragmentation and find that although disk evolution is sensitive to changes in opacity, the disks we study do not fragment, even at high resolution and even for extreme assumptions.

It is generally thought that protoplanetary disks embedded in envelopes are more massive and thus more susceptible to gravitational instabilities (GIs) than exposed disks. We present three-dimensional radiative hydrodynamics simulations of protoplanetary disks with the presence of envelope irradiation. For a disk with a radius of 40 AU and a mass of 0.07 Msun around a young star of 0.5 Msun, envelope irradiation tends to weaken and even suppress GIs as the irradiating flux is increased. The global mass transport induced by GIs is dominated by lower-order modes, and irradiation preferentially suppresses higher-order modes. As a result, gravitational torques and mass inflow rates are actually increased by mild irradiation. None of the simulations produce dense clumps or rapid cooling by convection, arguing against direct formation of giant planets by disk instability, at least in irradiated disks. However, dense gas rings and radial mass concentrations are produced, and these might be conducive to accelerated planetary core formation. Preliminary results from a simulation of a massive embedded disk with physical characteristics similar to one of the disks in the embedded source L1551 IRS5 indicate a long radiative cooling time and no fragmentation. The GIs in this disk are dominated by global two and three-armed modes.

In this dissertation, I describe theoretical and numerical studies that address the three-dimensional behavior of spiral shocks in protoplanetary disks and the controversial topic of gas giant formation by disk instability. For this work, I discuss characteristics of gravitational instabilities (GIs) in bursting and asymptotic phase disks; outline a theory for the three-dimensional structure of spiral shocks, called shock bores, for isothermal and adiabatic gases; consider convection as a source of cooling for protoplanetary disks; investigate the effects of opacity on disk cooling; use multiple analyses to test for disk stability against fragmentation; test the sensitivity of GI behavior to radiation boundary conditions; measure shock strengths and frequencies in GI-bursting disks; evaluate temperature fluctuations in unstable disks; and investigate whether spiral shocks can form chondrules when GIs activate. The numerical methods developed for these studies are discussed, including a radiation transport routine that explicitly couples the low and high optical depth regimes and a routine that models ortho and parahydrogen. Finally, I explore the hypothesis that chondrule formation and the FU Ori phenomenon are driven by GI activation in dead zones.

Recent three-dimensional radiative hydrodynamics simulations of protoplanetary disks report disparate disk behaviors, and these differences involve the importance of convection to disk cooling, the dependence of disk cooling on metallicity, and the stability of disks against fragmentation and clump formation. To guarantee trustworthy results, a radiative physics algorithm must demonstrate the capability to handle both the high and low optical depth regimes. We develop a test suite that can be used to demonstrate an algorithm's ability to relax to known analytic flux and temperature distributions, to follow a contracting slab, and to inhibit or permit convection appropriately. We then show that the radiative algorithm employed by Meji\'a (2004) and Boley et al. (2006) and the algorithm employed by Cai et al. (2006) and Cai et al. (2007, in prep.) pass these tests with reasonable accuracy. In addition, we discuss a new algorithm that couples flux-limited diffusion with vertical rays, we apply the test suite, and we discuss the results of evolving the Boley et al. (2006) disk with this new routine. Although the outcome is significantly different in detail with the new algorithm, we obtain the same qualitative answers. Our disk does not cool fast due to convection, and it is stable to fragmentation. We find an effective $\alpha\approx 10^{-2}$. In addition, transport is dominated by low-order modes.

The gas equation of state may be one of the critical factors for the disk instability theory of gas giant planet formation. This letter addresses the treatment of H$_2$ in hydrodynamical simulations of gravitationally unstable disks. In our discussion, we point out possible consequences of erroneous specific internal energy relations, approximate specific internal energy relations with discontinuities, and assumptions of constant $\Gamma_1$. In addition, we consider whether the ortho/para ratio for H$_2$ in protoplanetary disks should be treated dynamically as if the species are in equilibrium. Preliminary simulations indicate that the correct treatment is particularly critical for the study of gravitational instability when $T = 30$-50 K.

This paper presents a fully three-dimensional radiative hydrodymanics simulation with realistic opacities for a gravitationally unstable 0.07 Msun disk around a 0.5 Msun star. We address the following aspects of disk evolution: the strength of gravitational instabilities under realistic cooling, mass transport in the disk that arises from GIs, comparisons between the gravitational and Reynolds stresses measured in the disk and those expected in an alpha-disk, and comparisons between the SED derived for the disk and SEDs derived from observationally determined parameters. The mass transport in this disk is dominated by global modes, and the cooling times are too long to permit fragmentation for all radii. Moreover, our results suggest a plausible explanation for the FU Ori outburst phenomenon.

Observational studies show that the probability of finding gas giant planets around a star increases with the star's metallicity. Our latest simulations of disks undergoing gravitational instabilities (GIs) with realistic radiative cooling indicate that protoplanetary disks with lower metallicity generally cool faster and thus show stronger overall GI-activity. More importantly, the global cooling times in our simulations are too long for disk fragmentation to occur, and the disks do not fragment into dense protoplanetary clumps. Our results suggest that direct gas giant planet formation via disk instabilities is unlikely to be the mechanism that produced most observed planets. Nevertheless, GIs may still play an important role in a hybrid scenario, compatible with the observed metallicity trend, where structure created by GIs accelerates planet formation by core accretion.

In this paper, we describe the nonlinear outcome of spiral shocks in protoplanetary disks. Spiral shocks, for most protoplanetary disk conditions, create a loss of vertical force balance in the post-shock region and result in rapid expansion of the gas perpendicular to the disk midplane. This expansion has characteristics similar to hydraulic jumps, which occur in incompressible fluids. We present a theory to describe the behavior of these hybrids between shocks and hydraulic jumps (shock bores) and then compare the theory to three-dimensional hydrodynamics simulations. We discuss the fully three-dimensional shock structures that shock bores produce and discuss possible consequences for disk mixing, turbulence, and evolution of solids.