We discuss major challenges in modeling giant impacts between planetary bodies, focusing on the equations of state (EOS). During the giant impact stage of planet formation, rocky planets are melted and partially vaporized. However, most EOS models fail to reproduce experimental constraints on the thermodynamic properties of the major minerals over the required phase space. Here, we present an updated version of the widely-used ANEOS model that includes a user-defined heat capacity limit in the thermal free energy term. Our revised model for forsterite (Mg$_2$SiO$_4$), a common proxy for the mantles of rocky planets, provides a better fit to material data over most of the phase space of giant impacts. We discuss the limitations of this model and the Tillotson equation of state, a commonly used alternative model.

The radii and orbital periods of 4000+ confirmed/candidate exoplanets have been precisely measured by the Kepler mission. The radii show a bimodal distribution, with two peaks corresponding to smaller planets (likely rocky) and larger intermediate-size planets, respectively. While only the masses of the planets orbiting the brightest stars can be determined by ground-based spectroscopic observations, these observations allow calculation of their average densities placing constraints on the bulk compositions and internal structures. Yet an important question about the composition of planets ranging from 2 to 4 Earth radii still remains. They may either have a rocky core enveloped in a H2-He gaseous envelope (gas dwarfs) or contain a significant amount of multi-component, H2O-dominated ices/fluids (water worlds). Planets in the mass range of 10-15 Earth masses, if half-ice and half-rock by mass, have radii of 2.5 Earth radii, which exactly match the second peak of the exoplanet radius bimodal distribution. Any planet in the 2-4 Earth radii range requires a gas envelope of at most a few mass percentage points, regardless of the core composition. To resolve the ambiguity of internal compositions, we use a growth model and conduct Monte Carlo simulations to demonstrate that many intermediate-size planets are water worlds.

Context: The protostellar envelopes, outflow and large-scale chemistry of Class~0 and Class~I objects have been well-studied, but while previous works have hinted at or found a few Keplerian disks at the Class~0 stage, it remains to be seen if their presence in this early stage is the norm. Likewise, while complex organics have been detected toward some Class~0 objects, their distribution is unknown as they could reside in the hottest parts of the envelope, in the emerging disk itself or in other components of the protostellar system, such as shocked regions related to outflows. Aims: In this work, we aim to address two related issues regarding protostars: when rotationally supported disks form around deeply embedded protostars and where complex organic molecules reside in such objects. Methods: We observed the deeply embedded protostar, L483, using Atacama Large Millimeter/submillimeter Array (ALMA) Band~7 data from Cycles~1 and 3 with a high angular resolution down to $\sim$~0.1$^{\prime\prime}$ (20~au) scales. Results: We find that the kinematics of CS~$J=7$--$6$ and H$^{13}$CN~$J=4$--$3$ are best fitted by the velocity profile from infall under conservation of angular momentum and not by a Keplerian profile. The spatial extents of the observed complex organics are consistent with an estimated ice sublimation radius of the envelope at $\sim$~50~au, suggesting that the complex organics exist in the hot corino of L483. Conclusions: We find that L483 does not harbor a Keplerian disk down to at least $15$~au in radius. Instead, the innermost regions of L483 are undergoing a rotating collapse. This result highlights that some Class~0 objects contain only very small disks, or none at all, with the complex organic chemistry taking place on scales inside the hot corino of the envelope, in a region larger than the emerging disk.

[Abridged] The majority of stars form in binary or higher order systems. The Class 0 protostellar system IRAS16293-2422 contains two protostars, 'A' and 'B', separated by ~600 au and embedded in a single, 10^4 au scale envelope. Their relative evolutionary stages have been debated. We aim to study the relation and interplay between the two protostars A and B at spatial scales of 60 to ~1000 au. We selected molecular gas line transitions of CO, H2CO, HCN, CS, SiO, and CCH from the ALMA-PILS spectral imaging survey (329-363 GHz) and used them as tracers of kinematics, density, and temperature in the IRAS16293-2422 system. The angular resolution of the PILS data set allows us to study these quantities at a resolution of 0.5 arcsec (60 au [..]). Line-of-sight velocity maps of both optically thick and optically thin molecular lines reveal: (i) new manifestations of previously known outflows emanating from protostar A; (ii) a kinematically quiescent bridge of dust and gas spanning between the two protostars, with an inferred density between 4 10^4 and 3 10^7 cm^-3; and (iii) a separate, straight filament seemingly connected to protostar B seen only in CCH, with a flat kinematic signature. Signs of various outflows, all emanating from source A, are evidence of high-density and warmer gas; none of them coincide spatially and kinematically with the bridge. We hypothesize that the bridge arc is a remnant of filamentary substructure in the protostellar envelope material from which protostellar sources A and B have formed. One particular morphological structure appears to be due to outflowing gas impacting the quiescent bridge material. The continuing lack of clear outflow signatures unambiguously associated to protostar B and the vertically extended shape derived for its disk-like structure lead us to conclude that source B may be in an earlier evolutionary stage than source A.

Applying the survival function analysis to the planet radius distribution of the Kepler confirmed/candidate planets, we have identified two natural divisions of planet radius at 4 Earth radii and 10 Earth radii. These divisions place constraints on planet formation and interior structure model. The division at 4 Earth radii separates small exoplanets from large exoplanets above. When combined with the recently-discovered radius gap at 2 Earth radii, it supports the treatment of planets 2-4 Earth radii as a separate group, likely water worlds. For planets around solar-type FGK main-sequence stars, we argue that 2 Earth radii is the separation between water-poor and water-rich planets, and 4 Earth radii is the separation between gas-poor and gas-rich planets. We confirm that the slope of survival function in between 4 and 10 Earth radii to be shallower compared to either ends, indicating a relative paucity of planets in between 4-10 Earth radii, namely, the sub-Saturnian desert there. We name them transitional planets, as they form a bridge between the gas-poor small planets and gas giants. Accordingly, we propose the following classification scheme: (<2 Earth radii) rocky planets, (2-4 Earth radii) water worlds, (4-10 Earth radii) transitional planets, and (>10 Earth radii) gas giants.

The size distribution of exoplanets is a bimodal division into two groups: Rocky planet (<2 Earth radii) and water-rich planet (>2 Earth radii) with or without gaseous envelope.

Applying the survival function analysis to the planet radius distribution of the Kepler exoplanet candidates, we have identified two natural divisions of planet radius at 4 Earth radii and 10 Earth radii. These divisions place constraints on planet formation and interior structure model. The division at 4 Earth radii separates small exoplanets from large exoplanets above. When combined with the recently-discovered radius gap at 2 Earth radii, it supports the treatment of planets 2-4 Earth radii as a separate group, likely water worlds. Thus, for planets around solar-type FGK main-sequence stars, we argue that 2 Earth radii is the separation between water-poor and water-rich planets, and 4 Earth radii is the separation between gas-poor and gas-rich planets. We confirm that the slope of survival function in between 4 and 10 Earth radii to be shallower compared to either ends, indicating a relative paucity of planets in between 4-10 Earth radii, namely, the sub-Saturnian desert there. We name them transitional planets, as they form a bridge between the gas-poor small planets and gas giants. Accordingly, we propose the following classification scheme: (<2 Earth radii) rocky planets, (2-4 Earth radii) water worlds, (4-10 Earth radii) transitional planets, and (>10 Earth radii) gas giants.

Applying survival function analysis to the planet orbital period (P) and semi-major axis (a) distribution from the Kepler sample, we find that all exoplanets are uniformly distributed in (ln a) or (ln P), with an average inner cut-off of 0.05 AU to the host star. More specifically, this inner cut-off is 0.04 AU for rocky worlds (1-2 Earth radii) and 0.08 AU for water worlds (2-4 Earth radii). Moreover, the transitional planets (4-10 Earth radii) and gas giants (>10 Earth radii) have a change of slope of survival function at 0.4 AU from -1 to -1/2, suggesting a different statistical distribution uniform in \Sqrt[a] inside 0.4 AU, compared to small exoplanets (<4 Earth radii). This difference in distribution is likely caused by the difference in planet migration mechanism, and susceptibility to host stellar irradiation, for gas-poor (<4 Earth radii) versus gas-rich (>4 Earth radii) planets. Armed with this knowledge and combined with the survival function analysis of planet size distribution, we can make precise estimates of planet occurrence rate and predict the TESS mission yield.

The giant impact hypothesis remains the leading theory for lunar origin. However, current models struggle to explain the Moon's composition and isotopic similarity with Earth. Here we present a new lunar origin model. High-energy, high-angular momentum giant impacts can create a post-impact structure that exceeds the corotation limit (CoRoL), which defines the hottest thermal state and angular momentum possible for a corotating body. In a typical super-CoRoL body, traditional definitions of mantle, atmosphere and disk are not appropriate, and the body forms a new type of planetary structure, named a synestia. Using simulations of cooling synestias combined with dynamic, thermodynamic and geochemical calculations, we show that satellite formation from a synestia can produce the main features of our Moon. We find that cooling drives mixing of the structure, and condensation generates moonlets that orbit within the synestia, surrounded by tens of bars of bulk silicate Earth (BSE) vapor. The moonlets and growing moon are heated by the vapor until the first major element (Si) begins to vaporize and buffer the temperature. Moonlets equilibrate with BSE vapor at the temperature of silicate vaporization and the pressure of the structure, establishing the lunar isotopic composition and pattern of moderately volatile elements. Eventually, the cooling synestia recedes within the lunar orbit, terminating the main stage of lunar accretion. Our model shifts the paradigm for lunar origin from specifying a certain impact scenario to achieving a Moon-forming synestia. Giant impacts that produce potential Moon-forming synestias were common at the end of terrestrial planet formation.

Context. The Class 0 protostellar binary IRAS 16293-2422 is an interesting target for (sub)millimeter observations due to, both, the rich chemistry toward the two main components of the binary and its complex morphology. Its proximity to Earth allows the study of its physical and chemical structure on solar system scales using high angular resolution observations. Such data reveal a complex morphology that cannot be accounted for in traditional, spherical 1D models of the envelope. Aims. The purpose of this paper is to study the environment of the two components of the binary through 3D radiative transfer modeling and to compare with data from the Atacama Large Millimeter/submillimeter Array. Such comparisons can be used to constrain the protoplanetary disk structures, the luminosities of the two components of the binary and the chemistry of simple species. Methods. We present 13CO, C17O and C18O J=3-2 observations from the ALMA Protostellar Interferometric Line Survey (PILS), together with a qualitative study of the dust and gas density distribution of IRAS 16293-2422. A 3D dust and gas model including disks and a dust filament between the two protostars is constructed which qualitatively reproduces the dust continuum and gas line emission. Results and conclusions. Radiative transfer modeling of source A and B, with the density solution of an infalling, rotating collapse or a protoplanetary disk model, can match the constraints for the disk-like emission around source A and B from the observed dust continuum and CO isotopologue gas emission. If a protoplanetary disk model is used around source B, it has to have an unusually high scale-height in order to reach the dust continuum peak emission value, while fulfilling the other observational constraints. Our 3D model requires source A to be much more luminous than source B; LA ~ 18 $L_\odot$ and LB ~ 3 $L_\odot$.

Exoplanets smaller than Neptune are numerous, but the nature of the planet populations in the 1-4 Earth radii range remains a mystery. The complete Kepler sample of Q1-Q17 exoplanet candidates shows a radius gap at ~ 2 Earth radii, as reported by us in January 2017 in LPSC conference abstract #1576 (Zeng et al. 2017). A careful analysis of Kepler host stars spectroscopy by the CKS survey allowed Fulton et al. (2017) in March 2017 to unambiguously show this radius gap. The cause of this gap is still under discussion (Ginzburg et al. 2017; Lehmer & Catling 2017; Owen & Wu 2017). Here we add to our original analysis the dependence of the radius gap on host star type.

One of the open questions in astrochemistry is how complex organic and prebiotic molecules are formed. Aims. Our aim is to start the process of compiling an inventory of oxygen-bearing complex organic molecules toward the solar-type Class 0 protostellar binary IRAS16293-2422 from an unbiased spectral survey with ALMA (PILS). Here we focus on the new detections of ethylene oxide (c-C$_2$H$_4$O), acetone (CH$_3$COCH$_3$), and propanal (C$_2$H$_5$CHO). Methods. With ALMA, we surveyed the spectral range from 329 to 363 GHz at 0.5$"$ (60 AU diameter) resolution. Using a simple model for the molecular emission in LTE, the excitation temperatures and column densities of each species were constrained. Results. We successfully detect propanal (44 lines), ethylene oxide (20 lines) and acetone (186 lines) toward one component of the protostellar binary, IRAS16293B. The high resolution maps demonstrate that the emission for all investigated species originates from the compact central region close to the protostar. This, along with a derived common excitation temperature of $\sim$ 125 K, is consistent with a coexistence of these molecules in the same gas. Conclusions. The observations mark the first detections of acetone, propanal and ethylene oxide toward a low-mass protostar. The relative abundance ratios of the two sets of isomers (CH$_3$COCH$_3$/C$_2$H$_5$CHO $\sim$ 8 and CH$_3$CHO/c-C$_2$H$_4$O $\sim$ 12) are comparable to previous observations toward high-mass protostars. The majority of observed abundance ratios from these results as well as those measured toward high-mass protostars are up to an order of magnitude above the predictions from chemical models. This may reflect either missing reactions or uncertain rates in the chemical networks. The physical conditions, such as temperatures or densities, used in the models, may not be applicable to solar-type protostars either.

This work aims at exploring the scaling relations among rocky exoplanets. With the assumption of internal gravity increasing linearly in the core, and staying constant in the mantle, and tested against numerical simulations, a simple model is constructed, applicable to rocky exoplanets of core mass fraction (CMF) $\in 0.1\sim0.4$ and mass $\in 0.1 \sim 10 M_{\oplus}$. Various scaling relations are derived: (1) core radius fraction $\rm CRF \approx \sqrt{CMF}$, (2) Typical interior pressure $P_{\text{typical}} \sim g_{\rm{s}}^2$ (surface gravity squared), (3) core formation energy $E_{\rm diff} \sim \frac{1}{10} E_{\rm grav}$ (the total gravitational energy), (4) effective heat capacity of the mantle $C_{\rm p}\approx \left( \frac{M_{\rm{p}}}{M_{\oplus}} \right) \cdot 7 \cdot 10^{27}$ J K$^{-1}$, and (5) the moment of inertia $I\approx \frac{1}{3} \cdot M_{\rm{p}} \cdot R_{\rm{p}}^2$. These scaling relations, though approximate, are handy for quick use owing to their simplicity and lucidity, and provide insights into the interior structures of rocky exoplanets. As examples, this model is applied to several planets including Earth, GJ 1132b, Kepler-93b, and Kepler-20b, and made comparison with the numerical method.

How protostars accrete mass is one of the fundamental problems of star formation. High column densities and complex kinematical structures make direct observations challenging and they only provide a snapshot. Chemical tracers provide an interesting alternative to characterise the infall histories of protostars. Previous observations of H13CO+ towards the low-mass protostar IRAS15398-3359 showed a depression in the abundance. This is a sign of destruction of HCO+ by an enhanced presence of gaseous water in an extended region, possibly related to a recent burst in the accretion. Direct observations of water vapour can determine the exact extent of the emission and confirm the hypothesis that HCO+ is indeed a good tracer of the water snow-line. IRAS15398 was observed using ALMA at 0.5" resolution. Maps of HDO(101-000) and H218O(414-321) were taken simultaneously with observations of the CS(8-7) and N2H+(5-4) lines and continuum at 0.65 and 0.75 mm. The maps were interpreted using dust radiative transfer calculations of the protostellar infalling envelope with an outflow cavity. HDO is clearly detected and extended over the scales of the H13CO+ depression, although it is displaced by ~500 AU in the direction of the outflow. H218O is tentatively detected towards the red-shifted outflow lobe, but otherwise it is absent from the mapped region, which suggests that temperatures are low. Based on the temperature structure obtained from dust radiative transfer models, we conclude that the water was most likely released from the grains in an extended hour-glass configuration during a recent accretion burst. HDO is only detected in the region closest to the protostar, at distances of up to 500 AU. These signatures can only be explained if the luminosity has recently been increased by orders of magnitudes. Additionally, the densities in the outflow cones must be sufficiently low.

In the past few years, the number of confirmed planets has grown above 2000. It is clear that they represent a diversity of structures not seen in our own solar system. In addition to very detailed interior modeling, it is valuable to have a simple analytical framework for describing planetary structures. Variational principle is a fundamental principle in physics, entailing that a physical system follows the trajectory which minimizes its action. It is alternative to the differential equation formulation of a physical system. Applying this principle to planetary interior can beautifully summarize the set of differential equations into one, which provides us some insight into the problem. From it, a universal mass-radius relation, an estimate of error propagation from equation of state to mass-radius relation, and a form of virial theorem applicable to planetary interiors are derived.

The inner regions of the envelopes surrounding young protostars are characterised by a complex chemistry, with prebiotic molecules present on the scales where protoplanetary disks eventually may form. This paper introduces a systematic survey, "Protostellar Interferometric Line Survey (PILS)" of the Class 0 protostellar binary IRAS 16293-2422 using the Atacama Large Millimeter/submillimeter Array (ALMA). The survey covers the full frequency range from 329 to 363 GHz (0.8 mm) with additional targeted observations at 3.0 and 1.3 mm. More than 10,000 features are detected toward one component in the protostellar binary. Glycolaldehyde, its isomers, methyl formate and acetic acid, and its reduced alcohol, ethylene glycol, are clearly detected. For ethylene glycol both lowest state conformers, aGg' and gGg', are detected, the latter for the first time in the ISM. The abundance of glycolaldehyde is comparable to or slightly larger than that of ethylene glycol. In comparison to the Galactic Center, these two species are over-abundant relative to methanol, possibly an indication of formation at low temperatures in CO-rich ices. Both 13C and deuterated isotopologues of glycolaldehyde are detected, also for the first time ever in the ISM. For the deuterated species, a D/H ratio of approximately 5% is found with no differences between the deuteration in the different functional groups of glycolaldehyde. Measurements of the 13C-species lead to a 12C:13C ratio of approximately 30, lower than the typical ISM value. This low ratio may reflect an enhancement of 13CO in the ice due to either ion-molecule reactions in the gas before freeze-out or differences in the temperatures where 12CO and 13CO ices sublimate. The results reinforce the importance of low-temperature grain surface chemistry for the formation of prebiotic molecules seen here in the gas after sublimation of the entire ice mantle.

Formamide (NH$_2$CHO) has previously been detected in several star-forming regions and is thought to be a precursor for different prebiotic molecules. Its formation mechanism is still debated, however. Observations of formamide, related species and their isopotologues may provide useful clues to the chemical pathways leading to their formation. The Protostellar Interferometric Line Survey (PILS) represents an unbiased high angular resolution and sensitivity spectral survey of the low-mass protostellar binary IRAS 16293-2422 with the Atacama Large Millimeter/submillimeter Array (ALMA). We detect for the first time the three singly deuterated forms of NH$_2$CHO (NH$_2$CDO, cis- and trans-NHDCHO) as well as DNCO towards the component B of this binary source. The images reveal that the different isotopologues are all present in the same region. Based on observations of the $^{13}$C isotopologues of formamide and a standard $^{12}$C/$^{13}$C ratio, the deuterium fractionation is found to be similar for the three different forms with a value of about 2%. The DNCO/HNCO ratio is also comparable to the D/H ratio of formamide ($\sim$1%). These results are in agreement with the hypothesis that NH$_2$CHO and HNCO are chemically related through grain surface formation.

Several small dense exoplanets are now known, inviting comparisons to Earth and Venus. Such comparisons require translating their masses and sizes to composition models of evolved multi-layer-interior planets. Such theoretical models rely on our understanding of the Earth's interior, as well as independently derived equations of state (EOS), but have so far not involved direct extrapolations from Earth's seismic model -PREM. In order to facilitate more detailed compositional comparisons between small exoplanets and the Earth, we derive here a semi-empirical mass-radius relation for two-layer rocky planets based on PREM: ${\frac{R}{R_\oplus}} = (1.07-0.21\cdot \text{CMF})\cdot (\frac{M}{M_\oplus})^{1/3.7}$, where CMF stands for Core Mass Fraction. It is applicable to 1$\sim$8 M$_{\oplus}$ and CMF of 0.0$\sim$0.4. Applying this formula to Earth and Venus and several known small exoplanets with radii and masses measured to better than $\sim$30\% precision gives a CMF fit of $0.26\pm0.07$.

We analyze a sample of multiple-exoplanet systems which contain at least 3 transiting planets detected by the Kepler mission ("Kepler multiples"). We use a generalized Titius-Bode relation to predict the periods of 228 additional planets in 151 of these Kepler multiples. These Titius-Bode-based predictions suggest that there are, on average, ~2 planets in the habitable zone of each star. We estimate the inclination of the invariable plane for each system and prioritize our planet predictions by their geometric probability to transit. We highlight a short list of 77 predicted planets in 40 systems with a high geometric probability to transit, resulting in an expected detection rate of ~15%, ~3 times higher than the detection rate of our previous Titius-Bode-based predictions.