GW Ori is a hierarchical triple system which has a rare circumtriple disk. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and 12CO J=2-1 molecular gas emission of the disk. For the first time, we identify three dust rings in the disk at ~46, 188, and 338 AU, with estimated dust mass of ~70-250 Earth masses, respectively. To our knowledge, the outer ring in GW Ori is the largest dust ring ever found in protoplanetary disks. We use visibility modelling of dust continuum to show that the disk has misaligned parts and the innermost dust ring is eccentric. The disk misalignment is also suggested by the CO kinematics modelling. We interpret these substructures as evidence of ongoing dynamical interactions between the triple stars and the circumtriple disk.

We observe the dust continuum at 225 GHz and CO isotopologue (12CO, 13CO, and C18O) J=2-1 emission lines toward the CR Cha protoplanetary disk using the Atacama Large Millimeter/Submillimeter Array (ALMA). The dust continuum image shows a dust gap-ring structure in the outer region of the dust disk. A faint dust ring is also detected around 120 au beyond the dust gap. The CO isotopologue lines indicate that the gas disk is more extended than the dust disk. The peak brightness temperature of the 13CO line shows a small bump around 130 au while 12CO and C18O lines do not show. We investigate two possible mechanisms for reproducing the observed dust gap-ring structure and a gas temperature bump. First, the observed gap structure can be opened by a Jupiter mass planet using the relation between the planet mass and the gap depth and width. Meanwhile, the radiative transfer calculations based on the observed dust surface density profile show that the observed dust ring could be formed by dust accumulation at the gas temperature bump, that is, the gas pressure bump produced beyond the outer edge of the dust disk.

Major components of chondrites are chondrules and matrix. Measurements of the volatile abundance in Semarkona chondrules suggest that chondrules formed in a dense clump that had a higher solid density than the gas density in the solar nebula. We investigate collisions between chondrules and matrix in the surface region of dense clumps using fluffy aggregate growth models. Our simulations show that the collisional growth of aggregates composed of chondrules and matrix takes place in the clumps well before they experience gravitational collapse. The internal structure of chondrite parent bodies (CPBs) can be thereby determined by aggregate growth. We find that the aggregate growth generates two scales within CPBs. The first scale is involved with the small scale distribution of chondrules and determined by the early growth stage, where chondrules accrete aggregates composed of matrix grains. This accretion can reproduce the thickness of the matrix layer around chondrules found in chondrites. The other scale is related to the large scale distribution of chondrules. Its properties (e.g., the abundance of chondrules and the overall size) depend on the gas motion within the clump, which is parameterized in this work. Our work thus suggests that the internal structure of CPBs may provide important clues about their formation conditions and mechanisms.

The increasing number of newly detected exoplanets at short orbital periods raises questions about their formation and migration histories. A particular puzzle that requires explanation arises from one of the key results of the Kepler mission, namely the increase in the planetary occurrence rate with orbital period up to 10 days for F, G, K and M stars. We investigate the conditions for planet formation and migration near the dust sublimation front in protostellar disks around young Sun-like stars. For this analysis we use iterative 2D radiation hydrostatic disk models which include irradiation by the star, and dust sublimation and deposition depending on the local temperature and vapor pressure. We perform a parameter study by varying the magnetized turbulence onset temperature, the accretion stress, the dust mass fraction, and the mass accretion rate. Our models feature a gas-only inner disk, a silicate sublimation front and dust rim starting at around 0.08 au, an ionization transition zone with a corresponding density jump, and a pressure maximum which acts as a pebble trap at around 0.12 au. Migration torque maps show Earth- and super-Earth-mass planets halt in our model disks at orbital periods ranging from 10 to 22 days. Such periods are in good agreement with both the inferred location of the innermost planets in multiplanetary systems, and the break in planet occurrence rates from the Kepler sample at 10 days. In particular, models with small grains depleted produce a trap located at a 10-day orbital period, while models with a higher abundance of small grains present a trap at around a 17-day orbital period. The snow line lies at 1.6 au, near where the occurrence rate of the giant planets peaks. We conclude that the dust sublimation zone is crucial for forming close-in planets, especially when considering tightly packed super-Earth systems.

Understanding chondrule formation provides invaluable clues about the origin of the solar system. Recent studies suggest that planetesimal collisions and the resulting impact melts are promising for forming chondrules. Given that the dynamics of planetesimals is a key in impact-based chondrule formation scenarios, we here perform direct $N$-body simulations to examine how the presence of Jupiter affects the properties of chondrule-forming collisions. Our results show that the absence/presence of Jupiter considerably changes the properties of high velocity collisions whose impact velocities are higher than 2.5 km s$^{-1}$; high velocity collisions occur due to impacts between protoplanets and planetesimals for the case without Jupiter; for the case with Jupiter, eccentricities of planetesimals are pumped up by the secular and resonant perturbations from Jupiter. We also categorize the resulting planetesimal collisions and find that most of high velocity collisions are classified as grazing ones for both cases. To examine the effect of Jupiter on chondrule formation directly, we adopt the impact jetting scenario and compute the resulting abundance of chondrules. Our results show that for the case without Jupiter, chondrule formation proceeds in the inside-out manner, following the growth of protoplanets. If Jupiter is present, the location and timing of chondrule formation are determined by Jupiter's eccentricity, which is treated as a free parameter in our simulations. Thus, the existence of Jupiter is the key parameter for specifying when and where chondrule formation occurs for impact-based scenarios.

The Wide Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) will be the first high-performance stellar coronagraph using active wavefront control for deep starlight suppression in space, providing unprecedented levels of contrast, spatial resolution, and sensitivity for astronomical observations in the optical. One science case enabled by the CGI will be taking images and(R~50)spectra of faint interplanetary dust structures present in the habitable zone of nearby sunlike stars (~10 pc) and within the snow-line of more distant ones(~20pc), down to dust density levels commensurate with that of the solar system zodiacal cloud. Reaching contrast levels below~10-7 for the first time, CGI will cross an important threshold in debris disks physics, accessing disks with low enough optical depths that their structure is dominated by transport phenomena than collisions. Hence, CGI results will be crucial for determining how exozodiacal dust grains are produced and transported in low-density disks around mature stars. Additionally, CGI will be able to measure the brightness level and constrain the degree of asymmetry of exozodiacal clouds around individual nearby sunlike stars in the optical, at the ~10x solar zodiacal emission level. This information will be extremely valuable for optimizing the observational strategy of possible future exo-Earth direct imaging missions, especially those planning to operate at optical wavelengths, such as Habitable Exoplanet Observatory (HabEx) and the Large Ultraviolet/Optical/Infrared Surveyor (LUVOIR).

We report new Atacama Large Millimeter/submillimeter Array Band 3 (86-100 GHz; $\sim$80 mas angular resolution) and Band 4 (146-160 GHz; $\sim$50 mas angular resolution) observations of the dust continuum emission towards the archetypal and ongoing accretion burst young stellar object FU Ori, which simultaneously covered its companion, FU Ori S. In addition, we present near-infrared (2-2.45 $\mu$m) observations of FU Ori taken with the General Relativity Analysis via VLT InTerferometrY (GRAVITY; $\sim$1 mas angular resolution) instrument on the Very Large Telescope Interferometer (VLTI). We find that the emission in both FU Ori and FU Ori S at (sub)millimeter and near infrared bands is dominated by structures inward of $\sim$10 au radii. We detected closure phases close to zero from FU Ori with VLTI/GRAVITY, which indicate the source is approximately centrally symmetric and therefore is likely viewed nearly face-on. Our simple model to fit the GRAVITY data shows that the inner 0.4 au radii of the FU Ori disk has a triangular spectral shape at 2-2.45 $\mu$m, which is consistent with the H$_{2}$O and CO absorption features in a $\dot{M}\sim$10$^{-4}$ $M_{\odot}\,yr^{-1}$, viscously heated accretion disk. At larger ($\sim$0.4-10 au) radii, our analysis shows that viscous heating may also explain the observed (sub)millimeter and centimeter spectral energy distribution when we assume a constant, $\sim$10$^{-4}$ $M_{\odot}\,yr^{-1}$ mass inflow rate in this region. This explains how the inner 0.4 au disk is replenished with mass at a modest rate, such that it neither depletes nor accumulates significant masses over its short dynamic timescale. Finally, we tentatively detect evidence of vertical dust settling in the inner 10 au of the FU Ori disk, but confirmation requires more complete spectral sampling in the centimeter bands.

The origin of close-in Jovian planets is still elusive. We examine the in-situ gas accretion scenario as a formation mechanism of these planets. We reconstruct natal disk properties from the occurrence rate distribution of close-in giant planets, under the assumption that the occurrence rate may reflect the gas accretion efficiency onto cores of these planets. We find that the resulting gas surface density profile becomes an increasing function of the distance from the central star with some structure at $r \simeq 0.1$ au. This profile is quite different from the standard minimum-mass solar nebula model, while our profile leads to better reproduction of the population of observed close-in super-Earths based on previous studies. We compute the resulting magnetic field profiles and find that our profiles can be fitted by stellar dipole fields ($\propto r^{-3}$) in the vicinity of the central star and large-scale fields ($\propto r^{-2}$) at the inner disk regions, either if the isothermal assumption breaks down or if nonideal MHD effects become important. For both cases, the transition between these two profiles occurs at $r \simeq 0.1$ au, which corresponds to the period valley of giant exoplanets. Our work provides an opportunity to test the in-situ gas accretion scenario against disk quantities, which may constrain the gas distribution of the minimum-mass {\it extra}solar nebula.

Identification of the main planet formation site is fundamental to understanding how planets form and migrate to the current locations. We consider the heavy-element content trend of observed exoplanets derived from improved measurements of mass and radius, and explore how this trend can be used as a tracer of their formation sites. Using gas accretion recipes obtained from detailed hydrodynamical simulations, we confirm that the disk-limited gas accretion regime is most important for reproducing the heavy-element content trend. Given that such a regime is specified by two characteristic masses of planets, we compute these masses as a function of the distance ($r$) from the central star, and then examine how the regime appears in the mass-semimajor axis diagram. Our results show that a plausible solid accretion region emerges at $r \simeq 0.6$ au and expands with increasing $r$, using the conventional disk model. Given that exoplanets that possess the heavy-element content trend distribute currently near their central stars, our results imply the importance of planetary migration that would occur after solid accretion onto planets might be nearly completed at $r \geq 0.6$ au. Self-consistent simulations would be needed to verify the predictions herein.

In recent years, the disk populations in a number of young star-forming regions have been surveyed with ALMA. Understanding the disk properties and their correlation with those of the central star is critical to understand planet formation. In particular, a decrease of the average measured disk dust mass with the age of the region has been observed. We conducted high-sensitivity continuum ALMA observations of 43 Class II young stellar objects in CrA at 1.3 mm (230 GHz). The typical spatial resolution is 0.3". The continuum fluxes are used to estimate the dust masses of the disks, and a survival analysis is performed to estimate the average dust mass. We also obtained new VLT/X-Shooter spectra for 12 of the objects in our sample. 24 disks are detected, and stringent limits have been put on the average dust mass of the non-detections. Accounting for the upper limits, the average disk mass in CrA is $6\pm3\,\rm M_\oplus$, significantly lower than that of disks in other young (1-3 Myr) star forming regions (e.g. Lupus) and appears consistent with the 5-10 Myr old Upper Sco. The position of the stars in our sample on the HR diagram, however, seems to confirm that that CrA has age similar to Lupus. Neither external photoevaporation nor a lower than usual stellar mass distribution can explain the low disk masses. On the other hand, a low-mass disk population could be explained if the disks are small, which could happen if the parent cloud has a low temperature or intrinsic angular momentum, or if the the angular momentum of the cloud is removed by some physical mechanism such as magnetic braking. In order to fully explain and understand the dust mass distribution of protoplanetary disks and their evolution, it may also be necessary to take into consideration the initial conditions of star and disk formation process, which may vary from region to region, and affect planet formation.

Planned astronomical observatories of the 2020s will be capable of obtaining reflected light photometry and spectroscopy of cool extrasolar giant planets. Here we explain that such data are valuable both for understanding the origin and evolution of giant planets as a whole and for preparing for the interpretation of similar datasets from potentially habitable extrasolar terrestrial planets in the decades to follow.

As we begin to discover rocky planets in the habitable zone of nearby stars with missions like TESS and CHEOPS, we will need quick advancements on instrumentation and observational techniques that will enable us to answer key science questions, such as What are the atmospheric characteristics of habitable zone rocky planets? How common are Earth-like biosignatures in rocky planets?} How similar or dissimilar are those planets to Earth? Over the next decade we expect to have discovered several Earth-analog candidates, but we will not have the tools to study the atmospheres of all of them in detail. Ground-based ELTs can identify biosignatures in the spectra of these candidate exo-Earths and understand how the planets' atmospheres compare to the Earth at different epochs. Transit spectroscopy, high-resolution spectroscopy, and reflected-light direct imaging on ELTs can identify multiple biosignatures for habitable zone, rocky planets around M stars at optical to near-infrared wavelengths. Thermal infrared direct imaging can detect habitable zone, rocky planets around AFGK stars, identifying ozone and motivating reflected-light follow-up observations with NASA missions like HabEx/LUVOIR. Therefore, we recommend that the Astro2020 Decadal Survey Committee support: (1) the search for Earth-like biosignatures on rocky planets around nearby stars as a key science case; (2) the construction over the next decade of ground-based Extremely Large Telecopes (ELTs), which will provide the large aperture and spatial resolution necessary to start revealing the atmospheres of Earth-analogues around nearby stars; (3) the development of instrumentation that optimizes the detection of biosignatures; and (4) the generation of accurate line lists for potential biosignature gases, which are needed as model templates to detect those molecules.

Surveys reveal that terrestrial- to Neptune-sized planets (1 $< R <$ 4 R$_{\rm{Earth}}$) are the most common type of planets in our galaxy. Detecting and characterizing such small planets around nearby stars holds the key to understanding the diversity of exoplanets and will ultimately address the ubiquitousness of life in the universe. The following fundamental questions will drive research in the next decade and beyond: (1) how common are terrestrial to Neptune-sized planets within a few AU of their host star, as a function of stellar mass? (2) How does planet composition depend on planet mass, orbital radius, and host star properties? (3) What are the energy budgets, atmospheric dynamics, and climates of the nearest worlds? Addressing these questions requires: a) diffraction-limited spatial resolution; b) stability and achievable contrast delivered by adaptive optics; and c) the light-gathering power of extremely large telescopes (ELTs), as well as multi-wavelength observations and all-sky coverage enabled by a comprehensive US ELT Program. Here we provide an overview of the challenge, and promise of success, in detecting and comprehensively characterizing small worlds around the very nearest stars to the Sun with ELTs. This white paper extends and complements the material presented in the findings and recommendations published in the National Academy reports on Exoplanet Science Strategy and Astrobiology Strategy for the Search for Life in the Universe.

We highlight how guaranteed time observations (GTOs) and early release science (ERS) will advance understanding of exoplanet atmospheres and provide a glimpse into what transiting exoplanet science will be done with JWST during its first year of operations. These observations of 27 transiting planets will deliver significant insights into the compositions, chemistry, clouds, and thermal profiles of warm-to-hot gas-dominated planets well beyond what we have learned from HST, Spitzer, and other observatories to date. These data and insights will in turn inform our understanding of planet formation, atmospheric transport and climate, and relationships between various properties. Some insight will likely be gained into rocky planet atmospheres as well. JWST will be the most important mission for characterizing exoplanet atmospheres in the 2020s, and this should be considered in assessing exoplanet science for the 2020s and 2030s and future facilities.

Transmission spectra probe the atmospheres of transiting exoplanets, but these observations are also subject to signals introduced by magnetic active regions on host stars. Here we outline scientific opportunities in the next decade for providing useful constraints on stellar photospheres and inform interpretations of transmission spectra of the smallest ($R<4\,R_{\odot}$) exoplanets. We identify and discuss four primary opportunities: (1) refining stellar magnetic active region properties through exoplanet crossing events; (2) spectral decomposition of active exoplanet host stars; (3) joint retrievals of stellar photospheric and planetary atmospheric properties with studies of transmission spectra; and (4) continued visual transmission spectroscopy studies to complement longer-wavelength studies from $\textit{JWST}$. We make five recommendations to the Astro2020 Decadal Survey Committee: (1) identify the transit light source (TLS) effect as a challenge to precise exoplanet transmission spectroscopy and an opportunity ripe for scientific advancement in the coming decade; (2) include characterization of host star photospheric heterogeneity as part of a comprehensive research strategy for studying transiting exoplanets; (3) support the construction of ground-based extremely large telescopes (ELTs); (4) support multi-disciplinary research teams that bring together the heliophysics, stellar physics, and exoplanet communities to further exploit transiting exoplanets as spatial probes of stellar photospheres; and (5) support visual transmission spectroscopy efforts as complements to longer-wavelength observational campaigns with $\textit{JWST}$.

We have now accumulated a wealth of observations of the planet-formation environment and of mature planetary systems. These data allow us to test and refine theories of gas-giant planet formation by placing constraints on the conditions and timescale of this process. Yet a number of fundamental questions remain unanswered about how protoplanets accumulate material, their photospheric properties and compositions, and how they interact with protoplanetary disks. While we have begun to detect protoplanet candidates during the last several years, we are presently only sensitive to the widest separation, highest mass / accretion rate cases. Current observing facilities lack the angular resolution and inner working angle to probe the few-AU orbital separations where giant planet formation is thought to be most efficient. They also lack the contrast to detect accretion rates that would form lower mass gas giants and ice giants. Instruments and telescopes coming online over the next decade will provide high contrast in the inner giant-planet-forming regions around young stars, allowing us to build a protoplanet census and to characterize planet formation in detail for the first time.

The processes that transform gas and dust in circumstellar disks into diverse exoplanets remain poorly understood. One key pathway is to study exoplanets as they form in their young ($\sim$few~Myr) natal disks. Extremely Large Telescopes (ELTs) such as GMT, TMT, or ELT, can be used to establish the initial chemical conditions, locations, and timescales of planet formation, via (1)~measuring the physical and chemical conditions in protoplanetary disks using infrared spectroscopy and (2)~studying planet-disk interactions using imaging and spectro-astrometry. Our current knowledge is based on a limited sample of targets, representing the brightest, most extreme cases, and thus almost certainly represents an incomplete understanding. ELTs will play a transformational role in this arena, thanks to the high spatial and spectral resolution data they will deliver. We recommend a key science program to conduct a volume-limited survey of high-resolution spectroscopy and high-contrast imaging of the nearest protoplanetary disks that would result in an unbiased, holistic picture of planet formation as it occurs.

Provided that sufficient resources are deployed, we can look forward to an extraordinary future in which we will characterize potentially habitable planets. Until now, we have had to base interpretations of observations on habitability hypotheses that have remained untested. To test these theories observationally, we propose a statistical comparative planetology approach to questions of planetary habitability. The key objective of this approach will be to make quick and cheap measurements of critical planetary characteristics on a large sample of exoplanets, exploiting statistical marginalization to answer broad habitability questions. This relaxes the requirement of obtaining multiple types of data for a given planet, as it allows us to test a given hypothesis from only one type of measurement using the power of an ensemble. This approach contrasts with a "systems science" approach, where a few planets would be extensively studied with many types of measurements. A systems science approach is associated with a number of difficulties which may limit overall scientific return, including: the limited spectral coverage and noise of instruments, the diversity of exoplanets, and the extensive list of potential false negatives and false positives. A statistical approach could also be complementary to a systems science framework by providing context to interpret extensive measurements on planets of particular interest. We strongly recommend future missions with a focus on exoplanet characterization, and with the capability to study large numbers of planets in a homogenous way, rather than exclusively small, intense studies directed at a small sample of planets.

Ground-based telescopes coupled with adaptive optics (AO) have been playing a leading role in exoplanet direct imaging science and technological development for the past two decades and will continue to have an indispensable role for the next decade and beyond. Over the next decade, extreme AO systems on 8-10m telescopes will 1) mitigate risk for WFIRST-CGI by identifying numerous planets the mission can spectrally characterize, 2) validate performance requirements and motivate improvements to atmosphere models needed to unambiguously characterize solar system-analogues from space, and 3) mature novel technological innovations useful for space. Extremely Large Telescopes can deliver the first thermal infrared (10 $\mu m$) images of rocky planets around Sun-like stars and identify biomarkers. These data provide a future NASA direct imaging flagship mission (i.e. HabEx, LUVOIR) with numerous exo-Earth candidates and critical ancillary information to help clarify whether these planets are habitable.

We recommend an intensive effort to survey and understand the obliquity distribution of small close-in extrasolar planets over the coming decade. The orbital obliquities of exoplanets--i.e., the relative orientation between the planetary orbit and the stellar rotation--is a key tracer of how planets form and migrate. While the orbital obliquities of smaller planets are poorly explored today, a new generation of facilities coming online over the next decade will make such observations possible en masse. Transit spectroscopic observations with the extremely large telescopes will enable us to measure the orbital obliquities of planets as small as $\sim2R_{\oplus}$ around a wide variety of stars, opening a window into the orbital properties of the most common types of planets. This effort will directly contribute to understanding the formation and evolution of planetary systems, a key objective of the National Academy of Sciences' Exoplanet Science Strategies report.

The balloon-borne very long baseline interferometry (VLBI) experiment is a technical feasibility study for performing radio interferometry in the stratosphere. The flight model has been developed. A balloon-borne VLBI station will be launched to establish interferometric fringes with ground-based VLBI stations distributed over the Japanese islands at an observing frequency of approximately 20 GHz as the first step. This paper describes the system design and development of a series of observing instruments and bus systems. In addition to the advantages of avoiding the atmospheric effects of absorption and fluctuation in high frequency radio observation, the mobility of a station can improve the sampling coverage (`uv-coverage') by increasing the number of baselines by the number of ground-based counterparts for each observation day. This benefit cannot be obtained with conventional arrays that solely comprise ground-based stations. The balloon-borne VLBI can contribute to a future progress of research fields such as black holes by direct imaging.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the dust continuum emission at 1.3 mm and 12CO J=2-1 line emission of the transitional disk around DM Tau. DM Tau's disk is thought to possess a dust-free inner cavity inside a few au, from the absence of near-infrared excess on its spectral energy distribution (SED). Previous submillimeter observations were, however, unable to detect the cavity; instead, a dust ring ~20 au in radius was seen. The excellent angular resolution achieved in the new ALMA observations, 43 x 31 mas, allows discovery of a 4 au radius inner dust ring, confirming previous SED modeling results. This inner ring is symmetric in continuum emission, but asymmetric in 12CO emission. The known (outer) dust ring at ~20 au is recovered and shows azimuthal asymmetry with a strong-weak side contrast of ~1.3. The gap between these two rings is depleted by a factor of ~40 in dust emission relative to the outer ring. An extended outer dust disk is revealed, separated from the outer ring by another gap. The location of the inner ring is comparable to that of the main asteroid belt in the solar system. As a disk with a "proto-asteroid belt," the DM Tau system offers valuable clues to disk evolution and planet formation in the terrestrial planet forming region.

In spite of making a small contribution to total protoplanetary disk mass, dust affects the disk temperature by controlling absorption of starlight. As grains grow from their initial ISM-like size distribution, settling depletes the disk's upper layers of dust and decreases the optical depth, cooling the interior. Here we investigate the effect of collisional growth of dust grains and their dynamics on the thermal and optical profile of the disk, and explore the possibility that cooling induced by grain growth and settling could lead to gravitational instability. We develop a Monte Carlo dust collision model with a weighting technique and allow particles to collisionally evolve through sticking and fragmentation, along with vertical settling and turbulent mixing. We explore two disk models, the MMEN (minimum-mass extrasolar nebula), and a "heavy" disk with higher surface density than the MMEN, and perform simulations for both constant and spatially variable turbulence efficiency profile $\alpha(R,z)$. We then calculate mean wavelength-dependent opacities for the evolving disks and perform radiative transfer to calculate the temperature profile $T(R,z)$. Finally, we calculate the Toomre Q parameter, a measure of the disk's stability against self-gravity, for each disk model after it reaches a steady state dust-size distribution. We find that even weak turbulence can keep sub-micron sized particles stirred in the disk's upper layer, affecting its optical and thermal profiles, and the growth of large particles in the midplane can make a massive disk optically thick at millimeter wavelengths, making it difficult to calculate the surface density of dust available for planet formation in the inner disk. Also, for an initially massive disk, grain settling and growth can produce a drop in the Toomre Q parameter, driving the disk to $Q < 1.4$ and possibly triggering spiral instabilities.

We explore the origin of the trend of heavy elements in observed massive exoplanets. Coupling of better measurements of mass ($M_p$) and radius of exoplanets with planet structure models enables estimating the total heavy element mass ($M_Z$) in these planets. The corresponding relation is characterized by a power-law profile, $M_Z \propto M_p^{3/5}$. We develop a simplified, but physically motivated analysis to investigate how the power-law profile can be produced under the current picture of planet formation. Making use of the existing semi-analytical formulae of accretion rates of pebbles and planetesimals, our analysis shows that the relation can be reproduced well if it traces the final stage of planet formation. In the stage, planets accrete solids from gapped planetesimal disks and gas accretion is limited by disk evolution. We also find that dust accretion accompanying with gas accretion does not contribute to $M_Z$ for planets with $M_p < 10^3 M_{\oplus}$. Our findings are broadly consistent with that of previous studies, yet we explicitly demonstrate how planetesimal dynamics is crucial for better understanding the relation. While our approach is simple, we can also reproduce the trend of a correlation between planet metallicity and $M_p$ that is obtained by detailed population synthesis calculations, when the same assumption is adopted. Our analysis suggests that pebble accretion would not play a direct role at the final stage of planet formation, whereas radial drift of pebbles might be important indirectly for metal enrichment of planets. Detailed numerical simulations and more observational data are required for confirming our analysis.

Chondrites are one of the most primitive objects in the solar system, and keep the record of the degree of thermal metamorphism experienced in their parent bodies. This thermal history can be classified by the petrologic type. We investigate the thermal evolution of planetesimals to account for the current abundances (known as the fall statistics) of petrologic types 3 - 6 ordinary chondrites. We carry out a number of numerical calculations in which formation times and sizes of planetesimals are taken as parameters. We find that planetesimals that form within 2.0 Myr after the formation of Ca-Al-rich inclusions (CAIs) can contain all petrologic types of ordinary chondrites. Our results also indicate that plausible scenarios of planetesimal formation, which are consistent with the fall statistics, are that planetesimals with radii larger than 60 km start to form around 2.0 Myr after CAIs and/or that ones with radii less than 50 km should be formed within 1.5 Myr after CAIs. Thus, thermal modelling of planetesimals is important for revealing the occurrence and amount of metamorphosed chondrites, and for providing invaluable insights into planetesimal formation.

The Greenland Telescope project has recently participated in an experiment to image the supermassive black hole shadow at the center of M87 using Very Long Baseline Interferometry technique in April of 2018. The antenna consists of the 12-m ALMA North American prototype antenna that was modified to support two auxiliary side containers and to withstand an extremely cold environment. The telescope is currently at Thule Air Base in Greenland with the long-term goal to move the telescope over the Greenland ice sheet to Summit Station. The GLT currently has a single cryostat which houses three dual polarization receivers that cover 84-96 GHz, 213-243 GHz and 271-377 GHz bands. A hydrogen maser frequency source in conjunction with high frequency synthesizers are used to generate the local oscillator references for the receivers. The intermediate frequency outputs of each receiver cover 4-8 GHz and are heterodyned to baseband for digitization within a set of ROACH-2 units then formatted for recording onto Mark-6 data recorders. A separate set of ROACH-2 units operating in parallel provides the function of auto-correlation for real-time spectral analysis. Due to the stringent instrumental stability requirements for interferometry a diagnostic test system was incorporated into the design. Tying all of the above equipment together is the fiber optic system designed to operate in a low temperature environment and scalable to accommodate a larger distance between the control module and telescope for Summit Station. A report on the progress of the above electronics instrumentation system will be provided.

Spatially resolved structures in protoplanetary disks hint at unseen planets. Previous imaging observations of the transitional disk around MWC 758 revealed an inner cavity, a ring-like outer disk, emission clumps, and spiral arms, all possibly generated by companions. We present ALMA dust continuum observations of MWC 758 at 0.87 millimeter (mm) wavelength with 43$\times$39 mas angular resolution (6.9$\times$6.2 AU) and 20 $\mu$Jy beam$^{-1}$ rms. The central sub-mm emission cavity is revealed to be eccentric; once deprojected, its outer edge can be well-fitted by an ellipse with an eccentricity of 0.1 and one focus on the star. The broad ring-like outer disk is resolved into three narrow rings with two gaps in between. The outer two rings tentatively show the same eccentricity and orientation as the innermost ring bounding the inner cavity. The two previously known dust emission clumps are resolved in both the radial and azimuthal directions, with radial widths equal to $\sim$4$\times$ the local scale height. Only one of the two spiral arms previously imaged in near-infrared (NIR) scattered light is revealed in ALMA dust emission, at a slightly larger stellocentric distance owing to projection effects. We also submit evidence of disk truncation at $\sim$100 AU based on comparing NIR imaging observations with models. The spirals, the north clump, and the truncated disk edge are all broadly consistent with the presence of one companion exterior to the spirals at roughly 100 AU.

We performed the new JVLA full polarization observations at 40-48 GHz (6.3-7.5 mm) towards the nearby ($d$ $=$147$\pm$3.4 pc) Class 0 YSO IRAS 16293-2422, and compare with the previous SMA observations reported by Rao et al. (2009; 2014). We observed the quasar J1407+2827 which is weakly polarized and can be used as a leakage term calibrator for $<$9 GHz observations, to gauge the potential residual polarization leakage after calibration. We did not detect Stokes Q, U, and V intensities from the observations of J1407+2827, and constrain (3-$\sigma$) the residual polarization leakage after calibration to be $\lesssim$0.3\%. We detect linear polarization from one of the two binary components of our target source, IRAS\,16293-2422\,B. The derived polarization position angles from our observations are in excellent agreement with those detected from the previous observations of the SMA, implying that on the spatial scale we are probing ($\sim$50-1000 au), the physical mechanisms for polarizing the continuum emission do not vary significantly over the wavelength range of $\sim$0.88-7.5 mm. We hypothesize that the observed polarization position angles trace the magnetic field which converges from large scale to an approximately face-on rotating accretion flow. In this scenario, magnetic field is predominantly poloidal on $>$100 au scales, and becomes toroidal on smaller scales. However, this interpretation remains uncertain due to the high dust optical depths at the central region of IRAS\,16293-2422\,B and the uncertain temperature profile. We suggest that dust polarization at wavelengths comparable or longer than 7\,mm may still trace interstellar magnetic field. Future sensitive observations of dust polarization in the fully optically thin regime will have paramount importance for unambiguously resolving the magnetic field configuration.

The next generation ground-based extremely large telescopes (ELTs) present incredible opportunities to discover and characterize diverse planetary systems, even potentially habitable worlds. Adaptive-optics assisted thermal-IR (3-14 micron) imaging is a powerful tool to study exoplanets with extant 6-12 meter telescopes. ELTs have the spatial resolution and sensitivity that offer an unparalleled expansion of the available discovery space. AO-assisted thermal-IR instruments on ELTs will be superior to JWST for high contrast imaging in the thermal-IR, and complementary to high contrast observations at shorter wavelengths, in space or with second-generation extreme AO instruments. With appropriate investments in instrumentation and pre-cursor observations, thermal-IR equipped ELTs could image the first terrestrial and super-earth planets around nearby stars, opening the door to characterization of potentially habitable planets from the ground and space.

We present ALMA 0.87 mm continuum, HCO+ J=4--3 emission line, and CO J=3--2 emission line data of the disk of material around the young, Sun-like star PDS 70. These data reveal the existence of a possible two component transitional disk system with a radial dust gap of 0."2 +/- 0."05, an azimuthal gap in the HCO+ J=4--3 moment zero map, as well as two bridge-like features in the gas data. Interestingly these features in the gas disk have no analogue in the dust disk making them of particular interest. We modeled the dust disk using the Monte Carlo radiative transfer code HOCHUNK3D (Whitney et al. 2013) using a two disk components. We find that there is a radial gap that extends from 15-60 au in all grain sizes which differs from previous work.

Dozens of habitable zone, approximately earth-sized exoplanets are known today. An emerging frontier of exoplanet studies is identifying which of these habitable zone, small planets are actually habitable (have all necessary conditions for life) and, of those, which are earth-like. Many parameters and processes influence habitability, ranging from the orbit through detailed composition including volatiles and organics, to the presence of geological activity and plate tectonics. While some properties will soon be directly observable, others cannot be probed by remote sensing for the foreseeable future. Thus, statistical understanding of planetary systems' formation and evolution is a key supplement to the direct measurements of planet properties. Probabilistically assessing parameters we cannot directly measure is essential to reliably assess habitability, to prioritizing habitable-zone planets for follow-up, and for interpreting possible biosignatures.

This whitepaper discusses the diversity of exoplanets that could be detected by future observations, so that comparative exoplanetology can be performed in the upcoming era of large space-based flagship missions. The primary focus will be on characterizing Earth-like worlds around Sun-like stars. However, we will also be able to characterize companion planets in the system simultaneously. This will not only provide a contextual picture with regards to our Solar system, but also presents a unique opportunity to observe size dependent planetary atmospheres at different orbital distances. We propose a preliminary scheme based on chemical behavior of gases and condensates in a planet's atmosphere that classifies them with respect to planetary radius and incident stellar flux.

A collision between two molecular clouds is one possible candidate for high-mass star formation. The HII region RCW 36, located in the Vela molecular ridge, contains a young star cluster with two O-type stars. We present new CO observations of RCW 36 with NANTEN2, Mopra, and ASTE using $^{12}$CO($J$ = 1-0, 2-1, 3-2) and $^{13}$CO($J$ = 2-1) line emissions. We have discovered two molecular clouds lying at the velocities $V_\mathrm{LSR} \sim$5.5 and 9 km s$^{-1}$. Both clouds are likely to be physically associated with the star cluster, as verified by the good spatial correspondence among the two clouds, infrared filaments, and the star cluster. We also found a high intensity ratio of $\sim$0.6-1.2 for CO $J$ = 3-2 / 1-0 toward both clouds, indicating that the gas temperature has been increased due to heating by the O-type stars. We propose that the O-type stars in RCW 36 were formed by a collision between the two clouds, with a relative velocity separation of 5 km s$^{-1}$. The complementary spatial distributions and the velocity separation of the two clouds are in good agreement with observational signatures expected for O-type star formation triggered by a cloud-cloud collision. We also found a displacement between the complementary spatial distributions of the two clouds, which we estimate to be 0.3 pc assuming the collision angle to be 45$^{\circ}$ relative to the line-of-sight. We estimate the collision timescale to be $\sim$10$^5$ yr. It is probable that the cluster age by Ellerbroek et al. (2013b) is dominated by the low-mass members which were not formed under the triggering by cloud-cloud collision, and that the O-type stars in the center of the cluster are explained by the collisional triggering independently from the low-mass star formation.

We made CO ($J$ = 1--0, 2--1, and 3--2) observations toward an H{\sc ii} region RCW~32 in the Vela Molecular Ridge. The CO gas distribution associated with the H{\sc ii} region was revealed for the first time at a high resolution of 22 arcsec. The results revealed three distinct velocity components which show correspondence with the optical dark lanes and/or H$\alpha$ distribution. Two of the components show complementary spatial distribution which suggests collisional interaction between them at a relative velocity of $\sim$4 km~s$^{-1}$. Based on these results, we present a hypothesis that cloud-cloud collision determined the cloud distribution and triggered formation of the exciting star ionizing RCW~32. The collision time scale is estimated from the cloud size and the velocity separation to be $\sim$2 Myrs and the collision terminated $\sim$1 Myr ago, which is consistent with an age of the exciting star and the associated cluster. By combing the previous works on the H{\sc ii} regions in the Vela Molecular Ridge, we argue that the majority, at least four, of the H{\sc ii} regions in the Ridge were formed by triggering of cloud-cloud collision.

We report a possibility that the high-mass star located in the HII region RCW 34 was formed by a triggering induced by a collision of molecular clouds. Molecular gas distributions of the $^{12}$CO and $^{13}$CO $J=$2-1, and $^{12}$CO $J=$3-2 lines toward RCW 34 were measured by using the NANTEN2 and ASTE telescopes. We found two clouds with the velocity ranges of 0-10 km s$^{-1}$ and 10-14 km s$^{-1}$. Whereas the former cloud as massive as ~2.7 x 10$^{4}$ Msun has a morphology similar to the ring-like structure observed in the infrared wavelengths, the latter cloud with the mass of ~10$^{3}$ Msun, which has not been recognized by previous observations, distributes just likely to cover the bubble enclosed by the other cloud. The high-mass star with the spectral types of O8.5V is located near the boundary of the two clouds. The line intensity ratio of $^{12}$CO $J=$3-2 / $J=$2-1 yields high values (~1.5) in the neighborhood of the high-mass star, suggesting that these clouds are associated with the massive star. We also confirmed that the obtained position-velocity diagram shows a similar distribution with that derived by a numerical simulation of the supersonic collision of two clouds. Using the relative velocity between the two clouds (~5 km s$^{-1}$), the collisional time scale is estimated to be $\sim$0.2 Myr with the assumption of the distance of 2.5 kpc. These results suggest that the high-mass star in RCW 34 was formed rapidly within a time scale of ~0.2 Myr via a triggering of cloud-cloud collision.

The recent rapid progress in observations of circumstellar disks and extrasolar planets has reinforced the importance of understanding an intimate coupling between star and planet formation. Under such a circumstance, it may be invaluable to attempt to specify when and how planet formation begins in star-forming regions and to identify what physical processes/quantities are the most significant to make a link between star and planet formation. To this end, we have recently developed a couple of projects. These include an observational project about dust growth in Class 0 YSOs and a theoretical modeling project of the HL Tauri disk. For the first project, we utilize the archive data of radio interferometric observations, and examine whether dust growth, a first step of planet formation, occurs in Class 0 YSOs. We find that while our observational results can be reproduced by the presence of large ($\sim$ mm) dust grains for some of YSOs under the single-component modified blackbody formalism, an interpretation of no dust growth would be possible when a more detailed model is used. For the second project, we consider an origin of the disk configuration around HL Tauri, focusing on magnetic fields. We find that magnetically induced disk winds may play an important role in the HL Tauri disk. The combination of these attempts may enable us to move towards a comprehensive understanding of how star and planet formation are intimately coupled with each other.

The mechanism of angular momentum transport in protoplanetary disks is fundamental to understand the distributions of gas and dust in the disks. The unprecedented, high spatial resolution ALMA observations taken toward HL Tau and subsequent radiative transfer modeling reveal that a high degree of dust settling is currently achieved at the outer part of the HL Tau disk. Previous observations however suggest a high disk accretion rate onto the central star. This configuration is not necessarily intuitive in the framework of the conventional viscous disk model, since efficient accretion generally requires a high level of turbulence, which can suppress dust settling considerably. We develop a simplified, semi-analytical disk model to examine under what condition these two properties can be realized in a single model. Recent, non-ideal MHD simulations are utilized to realistically model the angular momentum transport both radially via MHD turbulence and vertically via magnetically induced disk winds. We find that the HL Tau disk configuration can be reproduced well when disk winds are properly taken into account. While the resulting disk properties are likely consistent with other observational results, such an ideal situation can be established only if the plasma $\beta$ at the disk midplane is $\beta_0 \simeq 2 \times 10^4$ under the assumption of steady accretion. Equivalently, the vertical magnetic flux at 100 au is about 0.2 mG. More detailed modeling is needed to fully identify the origin of the disk accretion and quantitatively examine plausible mechanisms behind the observed gap structures in the HL Tau disk.

M16, the Eagle Nebula, is an outstanding HII region where extensive high-mass star formation is taking place in the Sagittarius Arm, and hosts the remarkable "pillars" observed with HST. We made new CO observations of the region in the 12CO J=1--0 and J=2--1 transitions with NANTEN2. These observations revealed for the first time that a giant molecular cloud of $\sim 1.3 \times 10^5$ \Msun \ is associated with M16, which is elongated vertically to the Galactic plane over 35 pc at a distance of 1.8 kpc. We found a cavity of the molecular gas of $\sim 10$ pc diameter toward the heart of M16 at \lbeq (16.95\degree, 0.85\degree), where more than 10 O-type stars and $\sim 400$ stars are associated, in addition to a close-by molecular cavity toward a Spitzer bubble N19 at \lbeq (17.06\degree, 1.0\degree). We found three velocity components which show spatially complementary distribution in the entire M16 giant molecular cloud (GMC) including NGC6611 and N19, suggesting collisional interaction between them. Based on the above results we frame a hypothesis that collision between the red-shifted and blue-shifted components at a relative of $\sim 10$ \kms \ triggered formation of the O-type stars in the M16 GMC in the last 1-2 Myr. The collision is two fold in the sense that one of the collisional interactions is major toward the M16 cluster and the other toward N19 with a RCW120 type, the former triggered most of the O star formation with almost full ionization of the parent gas, and the latter an O star formation in N19.

We are motivated by the recent measurements of dust opacity indices beta around young stellar objects (YSOs), which suggest that efficient grain growth may have occurred earlier than the Class I stage. The present work makes use of abundant archival interferometric observations at submillimeter,millimeter, and centimeter wavelength bands to examine grain growth signatures in the dense inner regions (<1000 AU) of nine Class 0 YSOs. A systematic data analysis is performed to derive dust temperatures, optical depths, and dust opacity indices based on single-component modified black body fittings to the spectral energy distributions (SEDs). The fitted dust opacity indices (beta) are in a wide range of 0.3 to 2.0 when single-component SED fitting is adopted. Four out of the nine observed sources show beta lower than 1.7, the typical value of the interstellar dust. Low dust opacity index (or spectral index) values may be explained by the effect of dust grain growth, which makes beta<1.7. Alternatively, the very small observed values of beta may be interpreted by the presence of deeply embedded hot inner disks, which only significantly contribute to the observed fluxes at long wavelength bands. This possibility can be tested by the higher angular resolution imaging observations of ALMA, or more detailed sampling of SEDs in the millimeter and centimeter bands. The beta values of the remaining five sources are close to or consistent with 1.7, indicating that grain growth would start to significantly reduce the values of beta no earlier than the late-Class 0 stage for these YSOs.

We present GPI polarized intensity imagery of HD 100453 in Y-, J-, and K1 bands which reveals an inner gap ($9 - 18$ au), an outer disk ($18-39$ au) with two prominent spiral arms, and two azimuthally-localized dark features also present in SPHERE total intensity images (Wagner 2015). SED fitting further suggests the radial gap extends to $1$ au. The narrow, wedge-like shape of the dark features appears similar to predictions of shadows cast by a inner disk which is misaligned with respect to the outer disk. Using the Monte Carlo radiative transfer code HOCHUNCK3D (Whitney 2013), we construct a model of the disk which allows us to determine its physical properties in more detail. From the angular separation of the features we measure the difference in inclination between the disks 45$^{\circ}$, and their major axes, PA = 140$^{\circ}$ east of north for the outer disk and 100$^{\circ}$for the inner disk. We find an outer disk inclination of $25 \pm 10^{\circ}$ from face-on in broad agreement with the Wagner 2015 measurement of 34$^{\circ}$. SPHERE data in J- and H-bands indicate a reddish disk which points to HD 100453 evolving into a young debris disk.

Chondrules are primitive materials in the Solar System. They are formed in the first about 3 Myr of the Solar System's history. This timescale is longer than that of Mars formation, and it is conceivable that protoplanets, planetesimals and chondrules might have existed simultaneously in the solar nebula. Due to protoplanets perturbation on the planetesimal dynamics and chondrule accretion on them, all the formed chondrules are unlikely to be accreted by planetesimals. We investigate the amount of chondrules accreted by planetesimals in such a condition. We assume that a protoplanet is in oligarchic growth, and we perform analytical calculations of chondrule accretion both by a protoplanet and by planetesimals. Through the oligarchic growth stage, planetesimals accrete about half of the formed chondrules. The smallest planetesimals get the largest amount of the chondrules, compared with the amount accreted by more massive planetesimals. We perform a parameter study and find that this fraction is not largely changed for a wide range of parameter sets.

We resolved FU Ori at 29-37 GHz using the JVLA with $\sim$0$''$.07 resolution, and performed the complementary JVLA 8-10 GHz observations, the SMA 224 GHz and 272 GHz observations, and compared with archival ALMA 346 GHz observations to obtain the SEDs. Our 8-10 GHz observations do not find evidence for the presence of thermal radio jets, and constrain the radio jet/wind flux to at least 90 times lower than the expected value from the previously reported bolometric luminosity-radio luminosity correlation. The emission at $>$29 GHz may be dominated by the two spatially unresolved sources, which are located immediately around FU Ori and its companion FU Ori S, respectively. Their deconvolved radii at 33 GHz are only a few au. The 8-346 GHz SEDs of FU Ori and FU Ori S cannot be fit by constant spectral indices (over frequency). The more sophisticated models for SEDs suggest that the $>$29 GHz emission is contributed by a combination of free-free emission from ionized gas, and thermal emission from optically thick and optically thin dust components. We hypothesize that dust in the innermost parts of the disks ($\lesssim$0.1 au) has been sublimated, and thus the disks are no more well shielded against the ionizing photons. The estimated overall gas and dust mass based on SED modeling, can be as high as a fraction of a solar mass, which is adequate for developing disk gravitational instability. Our present explanation for the observational data is that the massive inflow of gas and dust due to disk gravitational instability or interaction with a companion/intruder, was piled up at the few au scale due to the development of a deadzone with negligible ionization. The piled up material subsequently triggered the thermal and the MRI instabilities when the ionization fraction in the inner sub-au scale region exceeded a threshold value, leading to the high protostellar accretion rate.

Presolar grains are small particles found in meteorites through their isotopic compositions which are considerably different from those of materials in the Solar System. If some isotopes in presolar grains diffused out beyond their grain sizes when they were embedded in parent bodies of meteorites, their isotopic compositions could be washed out, and hence the grains cannot be identified as presolar grains any more. We explore this possibility for the first time by self-consistently simulating the thermal evolution of planetesimals and the diffusion length of $^{18}$O in presolar silicate grains. Our results show that presolar silicate grains smaller than $\sim$ 0.03 $\mu m$ cannot keep their original isotopic compositions even if the host planetesimals experienced maximum temperature as low as 600 $^{\circ}$C. Since this temperature corresponds to the one experienced by petrologic type 3 chondrites, the isotopic diffusion can constrain the size of presolar silicate grains discovered in such chondrites to be larger than $\sim$ 0.03 $\mu m$. We also find that the diffusion lengths of $^{18}$O reach $\sim$ 0.3-2 $\mu m$ in planetesimals that were heated up to 700-800 $^{\circ}$C. This indicates that, if the original size of presolar grains spans a range from $\sim$ 0.001 $\mu m$ to $\sim$ 0.3 $\mu m$ like that in the interstellar medium, the isotopic records of the presolar grains may be almost completely lost in such highly thermalized parent bodies. We propose that isotopic diffusion could be a key process to control the size distribution and abundance of presolar grains in some types of chondrites.

Chondritic meteorites contain unique spherical materials named chondrules: sub-mm sized silicate grains once melted in a high temperature condition in the solar nebula. We numerically explore one of chondrule forming processes, planetesimal collisions. Previous studies found that impact jetting via protoplanet-planetesimal collisions make chondrules with an amount of 1 % of impactors' mass, when impact velocity exceeds 2.5 km s$^{-1}$. Based on the mineralogical data of chondrules, undifferentiated planetesimals would be more suitable for chondrule-forming collisions than potentially differentiated protoplanets. We examine planetesimal-planetesimal collisions using a shock physics code and find two things: one is that planetesimal-planetesimal collisions produce the nearly same amount of chondrules as protoplanet-planetesimal collisions ($\sim$ 1 %). The other is that the amount of produced chondrules becomes larger as the impact velocity increases when two planetesimals collide with each other. We also find that progenitors of chondrules can be ejected from deeper regions of large targets (planetesimals or protoplanets) than small impactors (planetesimals). The composition of targets is therefore important to fully account for the mineralogical data of currently sampled chondrules.

We explore whether close-in super-Earths were formed as rocky bodies that failed to grow fast enough to become the cores of gas giants before the natal protostellar disk dispersed. We model the failed cores' inward orbital migration in the low-mass or type I regime, to stopping points at distances where the tidal interaction with the protostellar disk applies zero net torque. The three kinds of migration traps considered are those due to the dead zone's outer edge, the ice line, and the transition from accretion to starlight as the disk's main heat source. As the disk disperses, the traps move toward final positions near or just outside 1~au. Planets at this location exceeding about 3~M$_\oplus$ open a gap, decouple from their host trap, and migrate inward in the high-mass or type II regime to reach the vicinity of the star. We synthesize the population of planets formed in this scenario, finding that some fraction of the observed super-Earths can be failed cores. Most super-Earths formed this way have more than 4~M$_\oplus$, so their orbits when the disk disperses are governed by type II migration. These planets have solid cores surrounded by gaseous envelopes. Their subsequent photoevaporative mass loss is most effective for masses originally below about 6 M$_\oplus$. The failed core scenario suggests a division of the observed super-Earth mass-radius diagram into five zones according to the inferred formation history.

Dust trapping accelerates the coagulation of dust particles, and thus it represents an initial step toward the formation of planetesimals. We report $H$-band (1.6 um) linear polarimetric observations and 0.87 mm interferometric continuum observations toward a transitional disk around LkHa 330. As results, a pair of spiral arms were detected in the $H$-band emission and an asymmetric (potentially arm-like) structure was detected in the 0.87 mm continuum emission. We discuss the origin of the spiral arm and the asymmetric structure, and suggest that a massive unseen planet is the most plausible explanation. The possibility of dust trapping and grain growth causing the asymmetric structure was also investigated through the opacity index (beta) by plotting the observed SED slope between 0.87 mm from our SMA observation and 1.3 mm from literature. The results imply that grains are indistinguishable from ISM-like dust in the east side ($beta = 2.0 pm 0.5$), but much smaller in the west side $beta = 0.7^{+0.5}_{-0.4}$, indicating differential dust size distribution between the two sides of the disk. Combining the results of near-infrared and submillimeter observations, we conjecture that the spiral arms exist at the upper surface and an asymmetric structure resides in the disk interior. Future observations at centimeter wavelengths and differential polarization imaging in other bands (Y to K) with extreme AO imagers are required to understand how large dust grains form and to further explore the dust distribution in the disk.

Chondritic meteorites provide valuable opportunities to investigate the origins of the solar system. We explore impact jetting as a mechanism of chondrule formation and subsequent pebble accretion as a mechanism of accreting chondrules onto parent bodies of chondrites, and investigate how these two processes can account for the currently available meteoritic data. We find that when the solar nebula is $\le 5$ times more massive than the minimum-mass solar nebula at $a \simeq 2-3$ AU and parent bodies of chondrites are $\le 10^{24}$ g ($\le$ 500 km in radius) in the solar nebula, impact jetting and subsequent pebble accretion can reproduce a number of properties of the meteoritic data. The properties include the present asteroid belt mass, the formation timescale of chondrules, and the magnetic field strength of the nebula derived from chondrules in Semarkona. Since this scenario requires a first generation of planetesimals that trigger impact jetting and serve as parent bodies to accrete chondrules, the upper limit of parent bodies' masses leads to the following implications: primordial asteroids that were originally $\ge 10^{24}$ g in mass were unlikely to contain chondrules, while less massive primordial asteroids likely had a chondrule-rich surface layer. The scenario developed from impact jetting and pebble accretion can therefore provide new insights into the origins of the solar system.

We report new Karl G. Jansky Very Large Array (JVLA), 0$"$.5 angular resolution observations of linearly polarized continuum emission at 6.9 mm, towards the Class 0 young stellar object (YSO) NGC1333 IRAS4A. This target source is a collapsing dense molecular core, which was resolved at short wavelengths to have hourglass shaped B-field configuration. We compare these 6.9 mm observations with previous polarization Submillimeter Array (SMA) observations at 0.88 mm, which have comparable angular resolution ($\sim$0$"$7). We found that at the same resolution, the observed polarization position angles at 6.9 mm are slightly deviated from those observed at 0.88 mm. Due to the lower optical depth of the emission at 6.9 mm, and the potential effect of dust grain growth, the new JVLA observations are likely probing B-field alignments in regions interior to those sampled by the previous polarization observations at higher frequencies. Our understanding can be improved by more sensitive observations, and observations for the more extended spatial scales.

We report Submillimeter Array (SMA) 1.3 mm high angular resolution observations towards the four EXor type outbursting young stellar objects (YSOs) VY Tau, V1118 Ori, V1143 Ori, and NY Ori. The data mostly show low dust masses $M_{dust}$ in the associated circumstellar disks. Among the sources, NY Ori possesses a relatively massive disk with $M_{dust} \sim 9 \times 10^{-4}$ $M_{\odot}$. V1118 Ori has a marginal detection equivalent to $M_{dust} \sim 6 \times 10^{-5}$ $M_{\odot}$. V1143 Ori has a non-detection also equivalent to $M_{dust} < 6 \times 10^{-5}$ $M_{\odot}$. For the nearest source VY Tau, we get a surprising non-detection which provides a stringent upper limit $M_{dust} < 6 \times 10^{-6}$ $M_{\odot}$. We interpret our findings as suggesting that the gas and dust reservoirs that feed the short duration, repetitive optical outbursts seen in some EXors may be limited to the small scale, innermost region of their circumstellar disks. This hot dust may have escaped our detection limits. Follow-up, more sensitive millimeter observations are needed to improve our understanding of the triggering mechanisms of EXor type outbursts.

We re-process the Atacama Large Millimeter/Submillimeter Array (ALMA) long-baseline science verification data taken toward HL Tauri. As shown by the previous work, we confirm that the high spatial resolution (~ 0."019, corresponding to ~ 2.7 AU) dust continuum images at \lambda = 0.87, 1.3, and 2.9 mm exhibit a multiple ring-like gap structure in the circumstellar disk. Assuming that the observed gaps are opened up by currently forming, unseen bodies, we estimate the mass of such hypothetical bodies based on following two approaches; the Hill radius analysis and a more elaborated approach developed from the angular momentum transfer analysis in gas disks. For the former, the measured gap widths are used for calibrating the mass of the bodies, while for the latter, the measured gap depths are utilized. We show that their masses are likely comparable to or less than the mass of Jovian planets, and then discuss an origin of the observed gap structure. By evaluating Toomre's gravitational instability (GI) condition and cooling effect, we find that the GI might be a possible mechanism to form the bodies in the outer region of the disk. As the disk might be gravitationally unstable only in the outer region of the disk, inward planetary migration would be needed to construct the current architecture of the hypothetical bodies. We estimate the gap-opening mass and show that type II migration might be able to play such a role. Combining GIs with inward migration, we conjecture that all of the observed gaps may be a consequence of bodies that might have originally formed at the outer part of the disk, and have subsequently migrated to the current locations. While ALMA's unprecedented high spatial resolution observations can revolutionize our picture of planet formation, more dedicated observational and theoretical studies are needed in order to fully understand the HL Tau images.

Chondrules are one of the most primitive elements that can serve as a fundamental clue as to the origin of our Solar system. We investigate a formation scenario of chondrules that involves planetesimal collisions and the resultant impact jetting. Planetesimal collisions are the main agent to regulate planetary accretion that corresponds to the formation of terrestrial planets and cores of gas giants. The key component of this scenario is that ejected materials can melt when the impact velocity between colliding planetesimals exceeds about 2.5 km s$^{-1}$. The previous simulations show that the process is efficient enough to reproduce the primordial abundance of chondrules. We examine this scenario carefully by performing semi-analytical calculations that are developed based on the results of direct $N$-body simulations. As found by the previous work, we confirm that planetesimal collisions that occur during planetary accretion can play an important role in forming chondrules. This arises because protoplanet-planetesimal collisions can achieve the impact velocity of about 2.5 km s$^{-1}$ or higher, as protoplanets approach the isolation mass ($M_{p,iso}$). Assuming that the ejected mass is a fraction ($F_{ch}$) of colliding planetesimals' mass, we show that the resultant abundance of chondrules is formulated well by $F_{ch}M_{p,iso}$, as long as the formation of protoplanets is completed within a given disk lifetime. We perform a parameter study and examine how the abundance of chondrules and their formation timing change. We find that the impact jetting scenario generally works reasonably well for a certain range of parameters, while more dedicated work would be needed to include other physical processes that are neglected in this work and to examine their effects on chondrule formation.

Layered accretion is one of the inevitable ingredients in protoplanetary disks when disk turbulence is excited by magnetorotational instabilities (MRIs). In the accretion, disk surfaces where MRIs fully operate have a high value of disk accretion rate ($\dot{M}$), while the disk midplane where MRIs are generally quenched ends up with a low value of $\dot{M}$. Significant progress on understanding MRIs has recently been made by a number of dedicated MHD simulations, which requires improvement of the classical treatment of $\alpha$ in 1D disk models. To this end, we obtain a new expression of $\alpha$ by utilizing an empirical formula that is derived from recent MHD simulations of stratified disks with Ohmic diffusion. It is interesting that this new formulation can be regarded as a general extension of the classical $\alpha$. Armed with the new $\alpha$, we perform a linear stability analysis of protoplanetary disks that undergo layered accretion, and find that a viscous instability can occur around the outer edge of dead zones. Disks become stable in using the classical $\alpha$. We identify that the difference arises from $\Sigma-$dependence of $\dot{M}$; whereas $\Sigma$ is uniquely determined for a given value of $\dot{M}$ in the classical approach, the new approach leads to $\dot{M}$ that is a multi-valued function of $\Sigma$. We confirm our finding both by exploring a parameter space as well as by performing the 1D, viscous evolution of disks. We finally discuss other non-ideal MHD effects that are not included in our analysis, but may affect our results.

We developed a dual-linear-polarization HEMT (High Electron Mobility Transistor) amplifier receiver system of the 45-GHz band (hereafter Z45), and installed it in the Nobeyama 45-m radio telescope. The receiver system is designed to conduct polarization observations by taking the cross correlation of two linearly-polarized components, from which we process full-Stokes spectroscopy. We aim to measure the magnetic field strength through the Zeeman effect of the emission line of CCS ($J_N=4_3-3_2$) toward pre-protostellar cores. A linear-polarization receiver system has a smaller contribution of instrumental polarization components to the Stokes $V$ spectra than that of the circular polarization system, so that it is easier to obtain the Stokes $V$ spectra. The receiver has an RF frequency of 42 $-$ 46 GHz and an intermediate frequency (IF) band of 4$-$8 GHz. The typical noise temperature is about 50 K, and the system noise temperature ranges from 100 K to 150K over the frequency of 42 $-$ 46 GHz. The receiver system is connected to two spectrometers, SAM45 and PolariS. SAM45 is a highly flexible FX-type digital spectrometer with a finest frequency resolution of 3.81 kHz. PolariS is a newly-developed digital spectrometer with a finest frequency resolution of 60 Hz, having a capability to process the full-Stokes spectroscopy. The Half Power Beam Width (HPBW) of the beam was measured to be 37$"$ at 43 GHz. The main beam efficiency of the Gaussian main beam was derived to be 0.72 at 43 GHz. The SiO maser observations show that the beam pattern is reasonably round at about 10 \% of the peak intensity and the side-lobe level was less than 3 \% of the peak intensity. Finally, we present some examples of astronomical observations using Z45.

A few particles of presolar Al2O3 grains with sizes above 0.5 mum are believed to have been produced in the ejecta of core-collapse supernovae (SNe). In order to clarify the formation condition of such large Al2O3 grains, we investigate the condensation of Al2O3 grains for wide ranges of the gas density and cooling rate. We first show that the average radius and condensation efficiency of newly formed Al2O3 grains are successfully described by a non-dimensional quantity "Lambda_on" defined as the ratio of the timescale with which the supersaturation ratio increases to the collision timescale of reactant gas species at dust formation. Then, we find that the formation of submicron-sized Al2O3 grains requires at least ten times higher gas densities than those presented by one-dimensional SN models. This indicates that presolar Al2O3 grains identified as a SN origin might be formed in dense gas clumps, allowing us to propose that the measured sizes of presolar grains can be a powerful tool to constrain the physical conditions in which they formed. We also briefly discuss the survival of newly formed Al2O3 grains against the destruction in the shocked gas within the SN remnants.

Photoevaporation due to high-energy stellar photons is thought to be one of the main drivers of protoplanetary disk dispersal. The fully or partially ionized disk surface is expected to produce free-free continuum emission at centimeter (cm) wavelengths that can be routinely detected with interferometers such as the upgraded Very Large Array (VLA). We use deep (rms noise down to 8 $\mu$Jy beam$^{-1}$ in the field of view center) 3.5 cm maps of the nearby (130 pc) Corona Australis (CrA) star formation (SF) region to constrain disk photoevaporation models. We find that the radio emission from disk sources in CrA is surprisingly faint. Only 3 out of 10 sources within the field of view are detected, with flux densities of order $10^2$ $\mu$Jy. However, a significant fraction of their emission is non-thermal. Typical upper limits for non-detections are $3\sigma\sim 60~\mu$Jy beam$^{-1}$. Assuming analytic expressions for the free-free emission from extreme-UV (EUV) irradiation, we derive stringent upper limits to the ionizing photon luminosity impinging on the disk surface $\Phi_\mathrm{EUV}<1-4\times10^{41}$ s$^{-1}$. These limits constrain $\Phi_\mathrm{EUV}$ to the low end of the values needed by EUV-driven photoevaporation models to clear protoplanetary disks in the observed few Myr timescale. Therefore, at least in CrA, EUV-driven photoevaporation is not likely to be the main agent of disk dispersal. We also compare the observed X-ray luminosities $L_X$ of disk sources with models in which photoevaporation is driven by such photons. Although predictions are less specific than for the EUV case, most of the observed fluxes (upper limits) are roughly consistent with the (scaled) predictions. Deeper observations, as well as predictions spanning a wider parameter space, are needed to properly test X-ray driven photoevaporation.

We present a new method of analysis for determining the surface geometry of five protoplanetary disks observed with near-infrared imaging polarimetry using Subaru-HiCIAO. Using as inputs the observed distribution of polarized intensity (PI), disk inclination, assumed properties for dust scattering, and other reasonable approximations, we calculate a differential equation to derive the surface geometry. This equation is numerically integrated along the distance from the star at a given position angle. We show that, using these approximations, the local maxima in the PI distribution of spiral arms (SAO 206462, MWC 758) and rings (2MASS J16042165-2130284, PDS 70) are associated with local concave-up structures on the disk surface. We also show that the observed presence of an inner gap in scattered light still allows the possibility of a disk surface that is parallel to the light path from the star, or a disk that is shadowed by structures in the inner radii. Our analysis for rings does not show the presence of a vertical inner wall as often assumed in studies of disks with an inner gap. Finally, we summarize the implications of spiral and ring structures as potential signatures of ongoing planet formation.

Massive exoplanets are observed preferentially around high metallicity ([Fe/H]) stars while low-mass exoplanets do not show such an effect. This so-called planet-metallicity correlation generally favors the idea that most observed gas giants at $r<10$ AU are formed via a core accretion process. We investigate the origin of this phenomenon using a semi-analystical model, wherein the standard core accretion takes place at planet traps in protostellar disks where rapid type I migrators are halted. We focus on the three major exoplanetary populations - hot-Jupiters, exo-Jupiters located at $r \simeq 1$ AU, and the low-mass planets. We show using a statistical approach that the planet-metallicity correlations are well reproduced in these models. We find that there are specific transition metallicities with values [Fe/H]$=-0.2$ to $-0.4$, below which the low-mass population dominates, and above which the Jovian populations take over. The exo-Jupiters significantly exceed the hot-Jupiter population at all observed metallicities. The low-mass planets formed via the core accretion are insensitive to metallicity, which may account for a large fraction of the observed super-Earths and hot-Neptunes. Finally, a controlling factor in building massive planets is the critical mass of planetary cores ($M_{c,crit}$) that regulates the onset of runaway gas accretion. Assuming the current data is roughly complete at [Fe/H]$>-0.6$, our models predict that the most likely value of the "mean" critical core mass of Jovian planets is $\braket{M_{c,crit}} \simeq 5 M_{\oplus}$ rather than $10 M_{\oplus}$. This implies that grain opacities in accreting envelopes should play an important role in lowering $M_{c,crit}$.

The Serpens South infrared dark cloud consists of several filamentary ridges, some of which fragment into dense clumps. On the basis of CCS ($J_N=4_3-3_2$), HC$_3$N ($J=5-4$), N$_2$H$^+$ ($J=1-0$), and SiO ($J=2-1, v=0$) observations, we investigated the kinematics and chemical evolution of these filamentary ridges. We find that CCS is extremely abundant along the main filament in the protocluster clump. We emphasize that Serpens South is the first cluster-forming region where extremely-strong CCS emission is detected. The CCS-to-N$_2$H$^+$ abundance ratio is estimated to be about 0.5 toward the protocluster clump, whereas it is about 3 in the other parts of the main filament. We identify six dense ridges with different $V_{\rm LSR}$. These ridges appear to converge toward the protocluster clump, suggesting that the collisions of these ridges may have triggered cluster formation. The collisions presumably happened within a few $\times \ 10^5$ yr because CCS is abundant only in such a short time. The short lifetime agrees with the fact that the number fraction of Class I objects, whose typical lifetime is $0.4 \times \ 10^5$ yr, is extremely high as about 70 percent in the protocluster clump. In the northern part, two ridges appear to have partially collided, forming a V-shape clump. In addition, we detected strong bipolar SiO emission that is due to the molecular outflow blowing out of the protostellar clump, as well as extended weak SiO emission that may originate from the filament collisions.

The ubiquity of planets poses an interesting question: when first planets are formed in galaxies. We investigate this problem by adopting a theoretical model developed for understanding the statistical properties of exoplanets. Our model is constructed as the combination of planet traps with the standard core accretion scenario in which the efficiency of forming planetary cores directly relates to the dust density in disks or the metallicity ([Fe/H]). We statistically compute planet formation frequencies (PFFs) as well as the orbital radius ($<R_{rapid}>$) within which gas accretion becomes efficient enough to form Jovian planets. The three characteristic exoplanetary populations are considered: hot Jupiters, exo-Jupiters densely populated around 1 AU, and low-mass planets such as super-Earths. We explore the behavior of the PFFs as well as $<R_{rapid}>$ for the three different populations as a function of metallicity ($-2 \leq$[Fe/H]$\leq -0.6$). We show that the total PFFs increase steadily with metallicity, which is the direct outcome of the core accretion picture. For the entire range of the metallicity considered here, the population of the low-mass planets dominates over the Jovian planets. The Jovian planets contribute to the PFFs above [Fe/H]-1. We find that the hot Jupiters form at lower metallcities than the exo-Jupiters. This arises from the radially inward transport of planetary cores by their host traps, which is more effective for lower metallicity disks due to the slower growth of the cores. The PFFs for the exo-Jupiters exceed those for the hot Jupiters around [Fe/H]-0.7. Finally, we show that the critical metallicity for forming Jovian planets is [Fe/H]-1.2, which is evaluated by comparing the values of $<R_{rapid}>$ between the hot Jupiters and the low-mass planets. The comparison intrinsically links to the different gas accretion efficiency between them.

Dust in protoplanetary disks is recognized as the building blocks of planets. In the core accretion scenario, the abundance of dust in disks (or metallicity) is crucial for forming cores of gas giants. We present our recent progress on the relationship between the metallicity and planet formation, wherein planet formation frequencies (PFFs) as well as the critical mass of planetary cores ($M_{c,crit}$) that can initiate gas accretion are statistically examined. We focus on three different planetary populations that are prominent for observed exoplanets in the mass-semimajor axis diagram: hot Jupiters, exo-Jupiters that are densely populated around 1 AU, and low-mass planets in tight orbits. We show that the resultant PFFs for both Jovian planets are correlated positively with the metallicity whereas low-mass planets form efficiently for a wide range of metallicities. This is consistent with the observed Planet-Metallicity correlation. Examining the statistically averaged value of $M_{c,crit}$ (defined as $<M_{c,crit}>$), we find that the correlation originates from the behavior of $<M_{c,crit}>$ that increases steadily with metallicity for two kinds of the Jovian planets while the low-mass planets obtain a rather constant value for $<M_{c,crit}>$. Such a difference in $<M_{c,crit}>$ can define transition metallicities (TMs) above which the Jovian planets gain a larger value of $<M_{c,crit}>$ than the low-mass planets, and hence gas giant formation takes place more efficiently. We find that TMs are sensitive to the important parameter that involves $M_{c,crit}$. We show, by comparing with the observations, that a most likely value of $M_{c,crit}$ is $\simeq 5M_{\oplus}$, which is smaller than the conventional value in the literature ($\simeq 10M_{\oplus}$). Our results suggest that opacities in planetary atmospheres play an important role for lowering $M_{c,crit}$.

We examine the effect of the dust size distribution on Kelvin-Helmholtz instabilities in the protoplanetary disk with dust sedimentation. With newly taking into account the dust size distribution, the growth rate of the Kelvin-Helmholtz instability is calculated using the linear stability analysis with the dust density distribution consistent with sedimentation. Dust abundance required for gravitational instabilities before the Kelvin-Helmholtz instability is derived from the linear stability analysis, and it is found that the required dust abundance significantly coincides with that estimated from the Richardson number. It is also found that when the dust size distribution is taken into account, the critical Richardson number for the onset of the Kelvin-Helmholtz instability tends to increase with dust abundance. This result is different from that in the case without the dust size distribution.

The rapid growth in the number of known exoplanets has revealed the existence of several distinct planetary populations in the observed mass-period diagram. Two of the most surprising are, (1) the concentration of gas giants around 1AU and (2) the accumulation of a large number of low-mass planets with tight orbits, also known as super-Earths and hot Neptunes. We have recently shown that protoplanetary disks have multiple planet traps that are characterized by orbital radii in the disks and halt rapid type I planetary migration. By coupling planet traps with the standard core accretion scenario, we showed that one can account for the positions of planets in the mass-period diagram. In this paper, we demonstrate quantitatively that most gas giants formed at planet traps tend to end up around 1 AU with most of these being contributed by dead zones and ice lines. In addition, we show that a large fraction of super-Earths and hot Neptunes are formed as "failed" cores of gas giants - this population being constituted by comparable contributions from dead zone and heat transition traps. Our results are based on the evolution of forming planets in an ensemble of disks where we vary only the lifetimes of disks as well as their mass accretion rates onto the host star. We show that a statistical treatment of the evolution of a large population of planetary cores initially caught in planet traps accounts for the existence of three distinct exoplantary populations - the hot Jupiters, the more massive planets at roughly orbital radii around 1 AU orbital, and the short period SuperEarths and hot Neptunes. There are very few evolutionary tracks that feed into the large orbital radii characteristic of the imaged Jovian planet, which agrees with recent surveys. Finally, we find that low-mass planets in tight orbits become the dominant planetary population for low mass stars ($M_* \le 0.7 M_{\odot}$).

Recent high-resolution observations show that protoplanetary disks have various kinds of structural properties or inhomogeneities. These are the consequence of a mixture of a number of physical and chemical processes taking place in the disks. Here, we discuss the results of our comprehensive investigations on how disk inhomogeneities affect planetary migration. We demonstrate that disk inhomogeneities give rise to planet traps - specific sites in protoplanetary disks at which rapid type I migration is halted. We show that up to three types of traps (heat transitions, ice lines and dead zones) can exist in a single disk, and that they move differently as the disk accretion rate decreases with time. We also demonstrate that the position of planet traps strongly depends on stellar masses and disk accretion rates. This indicates that host stars establish preferred (initial) scales of their planetary systems. Finally, we discuss the possible observational signatures of disk inhomogeneities.

Planetary migration is one of the most serious problems to systematically understand the observations of exoplanets. We clarify that the theoretically predicted type II migration is too fast, as well as type I migration, by developing detailed analytical arguments in which the timescale of type II migration is compared with the disk lifetime. In the disk-dominated regime, the type II migration timescale is characterized by a {\it local} viscous diffusion timescale, while the disk lifetime characterized by a {\it global} diffusion timescale that is much longer than the local one. Even in the planet-dominated regime where the inertia of the planet mass reduces the migration speed, the timescale is still shorter than the disk lifetime except in the final disk evolution stage where the total disk mass decays below the planet mass. This suggests that most giant planets plunge into the central stars within the disk lifetime, and it contradicts the exoplanet observations that gas giants are piled up at $r \ga 1$AU. We examine additional processes that may arise in protoplanetary disks: dead zones, photoevaporation of gas, and gas flow across a gap formed by a type II migrator. Although they make the type II migration timescale closer to the disk lifetime, we show that none of them can act as an effective barrier for rapid type II migration with the current knowledge of these processes. We point out that gas flow across a gap and the fraction of the flow accreted onto the planets are uncertain and they may have a potential to solve the problem. Much more detailed investigation for each process may be needed to explain the observed distribution of gas giants in extrasolar planetary systems.

In this paper, we reexamine the possibility of shear-driven turbulence during sedimentation of dust grains in the protoplanetary disk. The shear-driven turbulence is expected to occur before the onset of the gravitational instability for MMSN model. While according to previous studies without taking account of growth of dust grains, with the larger abundance of dust grains, the gravitational instability is indicated to occur before shear-driven turbulence. In this paper, the case with dust growth is considered, and it is found that the Kelvin-Helmholtz instability tends to occur before the gravitational instability even in the case with large abundance of dust grains. This is different from previous results without the dust growth.

The large number of observed exoplanets ($\gtrsim $ 700) provides important constraints on their origin as deduced from the mass-period diagram of planets. The most surprising features in the diagram are 1) the (apparent) pile up of gas giants at a period of $\sim 500$ days ($\sim1$ AU) and 2) the so-called mass-period relation which indicates that planetary mass is an increasing function of orbital period. We construct the evolutionary tracks of growing planets at planet traps in evolving protoplanetary disks and show that they provide a good physical understanding of how these observational properties arise. The fundamental feature of our model is that inhomogeneities in protoplanetary disks give rise to multiple (up to 3) trapping sites for rapid (type I) planetary migration of planetary cores. The viscous evolution of disks results in the slow radial movement of the traps and their cores from large to small orbital periods. In our model, the slow inward motion of planet traps is coupled with the standard core accretion scenario for planetary growth. As planets grow, type II migration takes over. Planet growth and radial movement are ultimately stalled by the dispersal of gas disks via photoevaporation. Our model makes a number of important predictions: that distinct sub-populations of planets that reflect the properties of planet traps where they have grown result in the mass-period relation; that the presence of these sub-populations naturally explains a pile-up of planets at $\sim 1$ AU; and that evolutionary tracks from the ice line do put planets at short periods and fill an earlier claimed "planet desert" - sparse population of planets in the mass-semi-major axis diagram.

Planetary migration provides a theoretical basis for the observed diversity of exoplanetary systems. We demonstrate that dust settling - an inescapable feature of disk evolution - gives even more rapid type I migration by up to a factor of about 2 than occurs in disks with fully mixed dust. On the other hand, type II migration becomes slower by a factor of 2 due to dust settling. This even more problematic type I migration can be resolved by the presence of a dead zone; the inner, high density region of a disk which features a low level of turbulence. We show that enhanced dust settling in the dead zone leaves a dusty wall at its outer edge. Back-heating of the dead zone by this wall produces a positive radial gradient for the disk temperature, which acts as a barrier for type I migration.

The structure of planetary systems around their host stars depends on their initial formation conditions. Massive planets will likely be formed as a consequence of rapid migration of planetesimals and low mass cores into specific trapping sites in protoplanetary discs. We present analytical modeling of inhomogeneities in protoplanetary discs around a variety of young stars, - from Herbig Ae/Be to classical T Tauri and down to M stars, - and show how they give rise to planet traps. The positions of these traps define the initial orbital distribution of multiple protoplanets. We investigate both corotation and Lindblad torques, and show that a new trap arises from the (entropy-related) corotation torque. This arises at that disc radius where disc heating changes from viscous to stellar irradiation dominated processes. We demonstrate that up to three traps (heat transitions, ice lines and dead zones) can exist in a single disc, and that they move differently as the disc accretion rate $\dot{M}$ decreases with time. The interaction between the giant planets which grow in such traps may be a crucial ingredient for establishing planetary systems. We also demonstrate that the position of planet traps strongly depends on stellar masses and disc accretion rates. This indicates that host stars establish preferred scales of planetary systems formed around them. We discuss the potential of planet traps induced by ice lines of various molecules such as water and CO, and estimate the maximum and minimum mass of planets which undergo type I migration. We finally apply our analyses to accounting for the initial conditions proposed in the Nice model for the origin of our Solar system.

Planetary migration is essential to explain the observed mass-period relation for exoplanets. Without some stopping mechanism, the tidal, resonant interaction between planets and their gaseous disc generally causes the planets to migrate inward so efficiently that they plunge into the host star within the gaseous disc lifetime ($\sim $ 1-3 Myrs). We investigate planetary migration by analytically calculating the migration rate and time within self-consistently computed, radiatively heated discs around M stars in which the effects of dust settling are included. We show that dust settling lowers the disc temperature and raises the gas density in the mid-plane. This inescapable evolution of disc structure speeds up type I planetary migration for lower mass bodies by up to a factor of about 2. We also examine the effects of dust settling on the gap-opening mass and type II migration, and find that the gap-opening mass is reduced by a factor of 2 and type II migration becomes slower by a factor of 2. While dust settling can somewhat alleviate the problem of planetary migration for more massive planets, the more rapid migration of low mass planets and planetary cores requires a robust slowing mechanism.

IceCube Collaboration Contributions to the 2009 International Cosmic Ray Conference

Over 5,000 PMTs are being deployed at the South Pole to compose the IceCube neutrino observatory. Many are placed deep in the ice to detect Cherenkov light emitted by the products of high-energy neutrino interactions, and others are frozen into tanks on the surface to detect particles from atmospheric cosmic ray showers. IceCube is using the 10-inch diameter R7081-02 made by Hamamatsu Photonics. This paper describes the laboratory characterization and calibration of these PMTs before deployment. PMTs were illuminated with pulses ranging from single photons to saturation level. Parameterizations are given for the single photoelectron charge spectrum and the saturation behavior. Time resolution, late pulses and afterpulses are characterized. Because the PMTs are relatively large, the cathode sensitivity uniformity was measured. The absolute photon detection efficiency was calibrated using Rayleigh-scattered photons from a nitrogen laser. Measured characteristics are discussed in the context of their relevance to IceCube event reconstruction and simulation efforts.

Planetary migration in standard models of gaseous protoplanetary disks is known to be very rapid ($\sim 10^5$ years) jeopardizing the existence of planetary systems. We present a new mechanism for significantly slowing rapid planetary migration, discovered by means of radiative transfer calculations of the thermal structure of protoplanetary disks irradiated by their central stars. Rapid dust settling in a disk's dead zone - a region with very little turbulence - leaves a dusty wall at its outer edge. We show that the back-heating of the dead zone by this irradiated wall produces a positive gradient of the disk temperature which acts as a thermal barrier to planetary migration which persists for the disk lifetime. Although we analyze in detail the migration of a Super-Earth in a low mass disk around an M star, our findings can apply to wide variety of young planetary systems. We compare our findings with other potentially important stopping mechanisms and show that there are large parameter spaces for which dead zones are likely to play the most important role for reproducing the observed mass-period relation in longer planetary periods.

The irradiation of protoplanetary discs by central stars is the main heating mechanism for discs, resulting in their flared geometric structure. In a series of papers, we investigate the deep links between 2D self-consistent disc structure and planetary migration in irradiated discs, focusing particularly on those around M stars. In this first paper, we analyse the thermal structure of discs that are irradiated by an M star by solving the radiative transfer equation by means of a Monte Carlo code. Our simulations of irradiated hydrostatic discs are realistic and self-consistent in that they include dust settling with multiple grain sizes (N=15), the gravitational force of an embedded planet on the disc, and the presence of a dead zone (a region with very low levels of turbulence) within it. We show that dust settling drives the temperature of the mid-plane from an $r^{-3/5}$ distribution (well mixed dust models) toward an $r^{-3/4}$. The dead zone, meanwhile, leaves a dusty wall at its outer edge because dust settling in this region is enhanced compared to the active turbulent disc at larger disc radii. The disc heating produced by this irradiated wall provides a positive gradient region of the temperature in the dead zone in front of the wall. This is crucially important for slowing planetary migration because Lindblad torques are inversely proportional to the disc temperature. Furthermore, we show that low turbulence of the dead zone is self-consistently induced by dust settling, resulting in the Kelvin-Helmholtz instability (KHI). We show that the strength of turbulence arising from the KHI in the dead zone is $\alpha=10^{-5}$.