Some Jupiter-mass exoplanets contain $\sim$$100\, M_\oplus$ of metals, well above the $\sim$$10\, M_\oplus$ typically needed in a solid core to trigger giant planet formation by runaway gas accretion. We demonstrate that such `heavy-metal Jupiters' can result from planetary mergers near $\sim$10 au. Multiple cores accreting gas at runaway rates gravitationally perturb one another onto crossing orbits such that the average merger rate equals the gas accretion rate. Concurrent mergers and gas accretion implies the core mass scales with the total planet mass as $M_{\rm core} \propto M^{1/5}$ - heavier planets harbour heavier cores, in agreement with the observed mass-metallicity relation. While the average gas giant merges about once to double its core, others may merge multiple times, as merger trees grow chaotically. We show that the dispersion of outcomes inherent in mergers can reproduce the large scatter in observed planet metallicities, assuming $3-30\, M_\oplus$ pre-runaway cores. Mergers potentially correlate metallicity, eccentricity, and spin.

Planets open gaps in discs. Gap opening is typically modeled by considering the planetary Lindblad torque which repels disc gas away from the planet's orbit. But gaps also clear because the planet consumes local material. We present a simple, easy-to-use, analytic framework for calculating how gaps deplete and how the disc's structure as a whole changes by the combined action of Lindblad repulsion and planetary consumption. The final mass to which a gap-embedded gas giant grows is derived in tandem. The analytics are tested against 1D numerical experiments and calibrated using published multi-dimensional simulations. In viscous alpha discs, the planet, while clearing a gap, initially accretes practically all of the gas that tries to diffuse past, rapidly achieving super-Jupiter if not brown dwarf status. By contrast, in inviscid discs---that may still accrete onto their central stars by, say, magnetized winds---planets open deep, repulsion-dominated gaps. Then only a small fraction of the disc accretion flow is diverted onto the planet, which grows to a fraction of a Jupiter mass. Transitional disc cavities might be cleared by families of such low-mass objects opening inviscid, repulsion-dominated, overlapping gaps which allow most of the outer disc gas to flow unimpeded onto host stars.

We present $L$'-band imaging of the PDS 70 planetary system with Keck/NIRC2 using the new infrared pyramid wavefront sensor. We detected both PDS 70 b and c in our images, as well as the front rim of the circumstellar disk. After subtracting off a model of the disk, we measured the astrometry and photometry of both planets. Placing priors based on the dynamics of the system, we estimated PDS 70 b to have a semi-major axis of $20^{+3}_{-4}$~au and PDS 70 c to have a semi-major axis of $34^{+12}_{-6}$~au (95\% credible interval). We fit the spectral energy distribution (SED) of both planets. For PDS 70 b, we were able to place better constraints on the red half of its SED than previous studies and inferred the radius of the photosphere to be 2-3~$R_{Jup}$. The SED of PDS 70 c is less well constrained, with a range of total luminosities spanning an order of magnitude. With our inferred radii and luminosities, we used evolutionary models of accreting protoplanets to derive a mass of PDS 70 b between 2 and 4 $M_{\textrm{Jup}}$ and a mean mass accretion rate between $3 \times 10^{-7}$ and $8 \times 10^{-7}~M_{\textrm{Jup}}/\textrm{yr}$. For PDS 70 c, we computed a mass between 1 and 3 $M_{\textrm{Jup}}$ and mean mass accretion rate between $1 \times 10^{-7}$ and $5 \times~10^{-7} M_{\textrm{Jup}}/\textrm{yr}$. The mass accretion rates imply dust accretion timescales short enough to hide strong molecular absorption features in both planets' SEDs.

Understanding the anomalous radii of many transiting hot gas giant planets is a fundamental problem of planetary science. Recent detections of re-inflated warm Jupiters orbiting post-main-sequence stars and the re-inflation of hot Jupiters while their host stars evolve on the main-sequence may help constrain models for the anomalous radii of hot Jupiters. In this work, we present evolution models studying the re-inflation of gas giants to determine how varying the depth and intensity of deposited heating affects both main-sequence re-inflation of hot Jupiters and post-main-sequence re-inflation of warm Jupiters. We find that deeper heating is required to re-inflate hot Jupiters than is needed to suppress their cooling, and that the timescale of re-inflation decreases with increasing heating rate and depth. We find a strong degeneracy between heating rate and depth, with either strong shallow heating or weak deep heating providing an explanation for main-sequence re-inflation of hot Jupiters. This degeneracy between heating rate and depth can be broken in the case of post-main-sequence re-inflation of warm Jupiters, as the inflation must be rapid to occur within post-main-sequence evolution timescales. We also show that the dependence of heating rate on incident stellar flux inferred from the sample of hot Jupiters can explain re-inflation of both warm and hot Jupiters. TESS will obtain a large sample of warm Jupiters orbiting post-main-sequence stars, which will help to constrain the mechanism(s) causing the anomalous radii of gas giant planets.

Black widows are close binary systems in which a millisecond pulsar is orbited by a companion a few per cent the mass of the sun. It has been suggested that the pulsar's rotationally powered $\gamma$-ray luminosity gradually evaporates the companion, eventually leaving behind an isolated millisecond pulsar. The evaporation efficiency is determined by the temperature $T_{\rm ch}\propto F^{2/3}$ to which the outflow is heated by the flux $F$ on a dynamical time-scale. Evaporation is most efficient for companions that fill their Roche lobes. In this case, the outflow is dominated by a cap around the L1 point with an angle $\theta_g\sim (T_{\rm ch}/T_g)^{1/2}$, and the evaporation time is $t_{\rm evap}=0.46(T_{\rm ch}/T_g)^{-2}\textrm{ Gyr}$, where $T_g>T_{\rm ch}$ is the companion's virial temperature. We apply our model to the observed black widow population, which has increased substantially over the last decade, considering each system's orbital period, companion mass, and pulsar spin-down power. While the original (Fruchter et al. 1988) black widow evaporates its companion on a few Gyr time-scale, direct evaporation on its own is too weak to explain the overall population. We propose instead that the evaporative wind couples to the companion's magnetic field, removes angular momentum from the binary, and maintains stable Roche-lobe overflow. While a stronger wind carries more mass, it also reduces the Alfv\'en radius, making this indirect magnetic braking mechanism less dependent on the flux $t_{\rm mag}\propto t_{\rm evap}^{1/3}$. This reduces the scatter in evolution times of observed systems, thus better explaining the combined black widow and isolated millisecond pulsar populations.

During the runaway phase of their formation, gas giants fill their gravitational spheres of influence out to Bondi or Hill radii. When runaway ends, planets shrink several orders of magnitude in radius until they are comparable in size to present-day Jupiter; in 1D models, the contraction occurs on the Kelvin-Helmholtz time-scale $t_{\rm KH}$, which is initially a few thousand years. However, if angular momentum is conserved, contraction cannot complete, as planets are inevitably spun up to their breakup periods $P_{\rm break}$. We consider how a circumplanetary disc (CPD) can de-spin a primordially magnetized gas giant and remove the centrifugal barrier, provided the disc is hot enough to couple to the magnetic field, a condition that is easier to satisfy at later times. By inferring the planet's magnetic field from its convective cooling luminosity, we show that magnetic spin-down times are shorter than contraction times throughout post-runaway contraction: $t_{\rm mag}/t_{\rm KH}\sim(P_{\rm break}/t_{\rm KH})^{1/21}\lesssim 1$. Planets can spin down until they corotate with the CPD's magnetospheric truncation radius, at a period $P_{\rm max}/P_{\rm break} \sim (t_{\rm KH}/P_{\rm break})^{1/7}$. By the time the disc disperses, $P_{\rm max}/P_{\rm break}\sim 20-30$; further contraction at fixed angular momentum can spin planets back up to $\sim$$10 P_{\rm break}$, potentially explaining observed rotation periods of giant planets and brown dwarfs.

Giant planets are thought to form by runaway gas accretion onto solid cores. Growth must eventually stop running away, ostensibly because planets open gaps (annular cavities) in their surrounding discs. Typical models stop runaway by artificially capping the accretion rate and lowering it to zero over an arbitrarily short time-scale. In reality, post-runaway accretion persists as long as the disc remains. During this final and possibly longest phase of formation, when the planet is still emerging from the disc, its mass can more than double, and its radius contracts by orders of magnitude. By drawing from the theory of how gaps clear, we find that post-runaway accretion luminosities diverge depending on disc viscosity: luminosities fall in low-viscosity discs but continue to rise past runaway in high-viscosity discs. This divergence amounts to a factor of $10^2$ by the time the disc disperses. Irrespective of the specifics of how planets interact with discs, the observed luminosity and age of an accreting planet can be used to calculate its instantaneous mass, radius, and accretion rate. We perform this exercise for the planet candidates embedded within the discs orbiting PDS 70, HD 163296, and MWC 758, inferring masses of $1-10\, M_{\rm J}$, accretion rates of $0.1-10\, M_{\rm J}$/Myr, and radii of $1-10\, R_{\rm J}$. Our radii are computed self-consistently from the planet's concurrent contraction and accretion and do not necessarily equal the value of $2R_{\rm J}$ commonly assumed; in particular, the radius depends on the envelope opacity as $R \propto \kappa^{0.5}$.

Gas giants are thought to form by runaway accretion: an instability driven by the self-gravity of growing atmospheres that causes accretion rates to rise super-linearly with planet mass. Why runaway should stop at a Jupiter or any other mass is unknown. We consider the proposal that final masses are controlled by circumstellar disc gaps (cavities) opened by planetary gravitational torques. We develop a fully time-dependent theory of gap formation and couple it self-consistently to planetary growth rates. When gaps first open, planetary torques overwhelm viscous torques, and gas depletes as if it were inviscid. In low-viscosity discs, of the kind motivated by recent observations and theory, gaps stay predominantly in this inviscid phase and planet masses finalize at $M_{\rm final}/M_\star\sim(\Omega t_{\rm disc})^{0.07}(H/a)^{2.73}(G\rho_0/\Omega^2)^{1/3}$, with $M_\star$ the host stellar mass, $\Omega$ the planet's orbital angular velocity, $t_{\rm disc}$ the gas disc's lifetime, $H/a$ its aspect ratio, and $\rho_0$ its unperturbed density. This final mass is independent of the dimensionless viscosity $\alpha$ and applies to large orbital distances, typically beyond $\sim$10 AU, where disc scale heights exceed planet radii. It evaluates to a few Jupiter masses at 10-100 AU, increasing gradually with distance as gaps become harder to open.

We present visible-light and ultraviolet (UV) observations of the supernova PTF12glz. The SN was discovered and monitored in near-UV and R bands as part of a joint GALEX and Palomar Transient Factory campaign. It is among the most energetic Type IIn supernovae observed to date (~10^51erg). If the radiated energy mainly came from the thermalization of the shock kinetic energy, we show that PTF12glz was surrounded by ~1 solar mass of circumstellar material (CSM) prior to its explosive death. PTF12glz shows a puzzling peculiarity: at early times, while the freely expanding ejecta are presumably masked by the optically thick CSM, the radius of the blackbody that best fits the observations grows at ~8000km/s. Such a velocity is characteristic of fast moving ejecta rather than optically thick CSM. This phase of radial expansion takes place before any spectroscopic signature of expanding ejecta appears in the spectrum and while both the spectroscopic data and the bolometric luminosity seem to indicate that the CSM is optically thick. We propose a geometrical solution to this puzzle, involving an aspherical structure of the CSM around PTF12glz. By modeling radiative diffusion through a slab of CSM, we show that an aspherical geometry of the CSM can result in a growing effective radius. This simple model also allows us to recover the decreasing blackbody temperature of PTF12glz. SLAB-Diffusion, the code we wrote to model the radiative diffusion of photons through a slab of CSM and evaluate the observed radius and temperature, is made available on-line.

Planets can open cavities (gaps) in the protoplanetary gaseous discs in which they are born by exerting gravitational torques. Viscosity counters these torques and limits the depletion of the gaps. We present a simple one-dimensional scheme to calculate the gas density profile inside gaps by balancing the gravitational and viscous torques. By generalizing the results of Goodman & Rafikov (2001), our scheme properly accounts for the propagation of angular momentum by density waves. This method allows us to easily study low-viscosity discs, which are challenging for full hydrodynamical simulations. We complement our numerical integration by analytical equations for the gap's steady-state depth and width as a function of the planet's to star's mass ratio $\mu$, the gas disc's aspect ratio $h$, and its Shakura & Sunyaev viscosity parameter $\alpha$. Specifically, we focus on low-mass planets ($\mu<\mu_{\rm th}\equiv h^3$) and identify a new low-viscosity regime, $\alpha<h(\mu/\mu_{\rm th})^5$, in which the classical analytical scaling relations are invalid. Equivalently, this low-viscosity regime applies to every gap that is depleted by more than a factor of $(\mu_{\rm th}/\mu)^3$ relative to the unperturbed density. We show that such gaps are significantly deeper and wider than previously thought, and consequently take a longer time to reach equilibrium.

Recent observations identify a valley in the radius distribution of small exoplanets, with planets in the range $1.5-2.0\,{\rm R}_{\oplus}$ significantly less common than somewhat smaller or larger planets. This valley may suggest a bimodal population of rocky planets that are either engulfed by massive gas envelopes that significantly enlarge their radius, or do not have detectable atmospheres at all. One explanation of such a bimodal distribution is atmospheric erosion by high-energy stellar photons. We investigate an alternative mechanism: the luminosity of the cooling rocky core, which can completely erode light envelopes while preserving heavy ones, produces a deficit of intermediate sized planets. We evolve planetary populations that are derived from observations using a simple analytical prescription, accounting self-consistently for envelope accretion, cooling and mass loss, and demonstrate that core-powered mass loss naturally reproduces the observed radius distribution, regardless of the high-energy incident flux. Observations of planets around different stellar types may distinguish between photoevaporation, which is powered by the high-energy tail of the stellar radiation, and core-powered mass loss, which depends on the bolometric flux through the planet's equilibrium temperature that sets both its cooling and mass-loss rates.

The orbits of many observed hot Jupiters are decaying rapidly due to tidal interaction, eventually reaching the Roche limit. We analytically study the ensuing coupled mass loss and orbital evolution during the Roche-lobe overflow and find two possible scenarios. Planets with light cores $M_c\lesssim 6M_\oplus$ (assuming a nominal tidal dissipation factor $Q\sim 10^6$ for the host star) are transformed into Neptune-mass gas planets, orbiting at a separation (relative to the stellar radius) $a/R_\star\approx 3.5$. Planets with heavier cores $M_c\gtrsim 6M_\oplus$ plunge rapidly until they are destroyed at the stellar surface. Remnant gas-Neptunes, which are stable to photo-evaporation, are absent from the observations, despite their unique transit radius ($5-10R_\oplus$). This result suggests that $M_c\gtrsim 6M_\oplus$, providing a useful constraint on the poorly-known core mass that may distinguish between different formation theories of gas giants. Alternatively, if one assumes a prior of $M_c\approx 6 M_\oplus$ from the core-accretion theory, our results suggest that $Q$ does not lie in the range $10^6\lesssim Q\lesssim 10^7$.

Short-period Earth to Neptune size exoplanets (super-Earths) with voluminous gas envelopes seem to be very common. These gas atmospheres are thought to have originated from the protoplanetary disk in which the planets were embedded during their first few Myr. The accretion rate of gas from the surrounding nebula is determined by the ability of the gas to cool and radiate away its gravitational energy. Here we demonstrate that heat from the tidal interaction between the star and the young (and therefore inflated) planet can inhibit the gas cooling and accretion. Quantitatively, we find that the growth of super-Earth atmospheres halts for planets with periods of about 10 days, provided that their initial eccentricities are of the order of 0.2. Thus, tidal heating provides a robust and simple mechanism that can simultaneously explain why these planets did not become gas giants and account for the deficit of low-density planets closer to the star, where the tides are even stronger. We suggest that tidal heating may be as important as other factors (such as the nebula's lifetime and atmosphere evaporation) in shaping the observed super-Earth population.

Some recently discovered short-period Earth to Neptune sized exoplanets (super Earths) have low observed mean densities which can only be explained by voluminous gaseous atmospheres. Here, we study the conditions allowing the accretion and retention of such atmospheres. We self-consistently couple the nebular gas accretion onto rocky cores and the subsequent evolution of gas envelopes following the dispersal of the protoplanetary disk. Specifically, we address mass-loss due to both photo-evaporation and cooling of the planet. We find that planets shed their outer layers (dozens of percents in mass) following the disk's dispersal (even without photo-evaporation), and their atmospheres shrink in a few Myr to a thickness comparable to the radius of the underlying rocky core. At this stage, atmospheres containing less particles than the core (equivalently, lighter than a few % of the planet's mass) can be blown away by heat coming from the cooling core, while heavier atmospheres cool and contract on a timescale of Gyr at most. By relating the mass-loss timescale to the accretion time, we analytically identify a Goldilocks region in the mass-temperature plane in which low-density super Earths can be found: planets have to be massive and cold enough to accrete and retain their atmospheres, while not too massive or cold, such that they do not enter runaway accretion and become gas giants (Jupiters). We compare our results to the observed super-Earth population and find that low-density planets are indeed concentrated in the theoretically allowed region. Our analytical and intuitive model can be used to investigate possible super-Earth formation scenarios.

Recent analyses of the orbits of some Kuiper Belt objects hypothesize the presence of an undiscovered Neptune-size planet at a very large separation from the Sun. The energy budget of Neptunes on such distant orbits is dominated by the internal heat released by their cooling rather than solar irradiation (making them effectively "isolated"). The blackbody radiation that these planets emit as they cool may provide the means for their detection. Here we use an analytical toy model to study the cooling and radiation of isolated Neptunes. This model can translate a detection (or a null detection) to a constraint on the size and composition of the hypothesized "Planet Nine". Specifically, the thick gas atmosphere of Neptune-like planets serves as an insulating blanket which slows down their cooling. Therefore, a measurement of the blackbody temperature, $T_{\rm eff}\sim 50\textrm{K}$, at which a Neptune emits can be used to estimate the mass of its atmosphere, $M_{\rm atm}$. Explicitly, we find the relation $T_{\rm eff}\propto M_{\rm atm}^{1/12}$. Despite this weak relation, a measurement of the flux at the Wien tail can constrain the atmospheric mass, at least to within a factor of a few, and provide useful limits to possible formation scenarios of these planets. Finally, we constrain the size and composition of Planet Nine by combining our model with the null results of recent all-sky surveys.

Many giant exoplanets in close orbits have observed radii which exceed theoretical predictions. One suggested explanation for this discrepancy is heat deposited deep inside the atmospheres of these "hot Jupiters". Here, we study extended power sources which distribute heat from the photosphere to the deep interior of the planet. Our analytical treatment is a generalization of a previous analysis of localized "point sources". We model the deposition profile as a power law in the optical depth and find that planetary cooling and contraction halt when the internal luminosity (i.e. cooling rate) of the planet drops below the heat deposited in the planet's convective region. A slowdown in the evolutionary cooling prior to equilibrium is possible only for sources which do not extend to the planet's center. We estimate the Ohmic dissipation resulting from the interaction between the atmospheric winds and the planet's magnetic field, and apply our analytical model to Ohmically heated planets. Our model can account for the observed radii of most inflated planets which have equilibrium temperatures $\approx 1500\textrm{K}-2500\textrm{K}$, and are inflated to a radius $\approx 1.6 R_J$. However, some extremely inflated planets remain unexplained by our model. We also argue that Ohmically inflated planets have already reached their equilibrium phase, and no longer contract. Following Wu & Lithwick who argued that Ohmic heating could only suspend and not reverse contraction, we calculate the time it takes Ohmic heating to re-inflate a cold planet to its equilibrium configuration. We find that while it is possible to re-inflate a cold planet, the re-inflation timescales are longer by a factor of $\approx 30$ than the cooling time.

Some extrasolar giant planets in close orbits---"hot Jupiters"---exhibit larger radii than that of a passively cooling planet. The extreme irradiation $L_{\rm eq}$ these hot Jupiters receive from their close in stars creates a thick isothermal layer in their envelopes, which slows down their convective cooling, allowing them to retain their inflated size for longer. This is yet insufficient to explain the observed sizes of the most inflated planets. Some models invoke an additional power source, deposited deep in the planet's envelope. Here we present an analytical model for the cooling of such irradiated, and internally heated gas giants. We show that a power source $L_{\rm dep}$, deposited at an optical depth $\tau_{\rm dep}$, creates an exterior convective region, between optical depths $L_{\rm eq}/L_{\rm dep}$ and $\tau_{\rm dep}$, beyond which a thicker isothermal layer exists, which in extreme cases may extend to the center of the planet. This convective layer, which occurs only for $L_{\rm dep}\tau_{\rm dep}>L_{\rm eq}$, further delays the cooling of the planet. Such a planet is equivalent to a planet irradiated with $L_{\rm eq}\left(1+L_{\rm dep}\tau_{\rm dep}/L_{\rm eq}\right)^\beta$, where $\beta\approx 0.35$ is an effective power-law index describing the radiative energy density as function of the optical depth for a convective planet $U\propto\tau^\beta$. Our simple analytical model reproduces the main trends found in previous numerical works, and provides an intuitive understanding. We derive scaling laws for the cooling rate of the planet, its central temperature, and radius. These scaling laws can be used to estimate the effects of tidal or Ohmic dissipation, wind shocks, or any other mechanism involving energy deposition, on sizes of hot Jupiters.

Recent observations suggest that some supernovae may be the result of an explosion into an optically thick circumstellar material, the product of pre-explosion mass-loss (wind) by the progenitor star. This scenario has been studied previously both analytically and numerically. However, many previous studies base their analysis on the diffusion approximation for radiation transfer, which is inappropriate in the optically thin outer layers of the wind. Here we study the deviations from diffusion, and calculate light curves more accurately using a Monte Carlo approach to photon transfer. We distinguish between "compact" winds, for which the diffusion approximation is appropriate, and "extended" winds, which require a more delicate treatment of the radiation. We show that this effect is more significant than that of the light travel time difference to a distant observer, which has a secondary influence on the light curves of extended-wind systems. We also comment on the applicability of the widely used flux-limited diffusion approximation in this context: we find that it generally does not reproduce the Monte Carlo results. The flux-limited diffusion approximation leads to results which are not only quantitatively, but also qualitatively wrong, in the extended-wind regime.

Observations from the last decade have indicated the existence of a general class of superluminous supernovae (SLSNe), in which the peak luminosity exceeds 10^{44} erg/s. Here we focus on a subclass of these events, where the light curve is also tens of days wide, so the total radiated energy is order 10^{51} erg. If the origin of these SLSNe is a core-collapse-driven explosion of a massive star, then the mechanism which converts the explosion energy into radiation must be very efficient (much more than in typical core collapse SNe, where this efficiency is of order one percent). We examine the scenario where the radiated luminosity is due to efficient conversion of kinetic energy of the ejected stellar envelope into radiation by interaction with an optically thick, pre-existing circumstellar material (CSM), presumably the product of a steady wind from the progenitor. We base the analysis on a simple, numerically solved, hydrodynamic diffusion model, which allows us to identify the qualitative behavior of the observable light curves, and to relate them to the parameters of the wind. We specifically show that a wide and superluminous supernova requires the mass of the relevant wind material to be comparable to that of ejected material from the exploding progenitor. We find the wind parameters which explain the peak luminosity and width of the bolometric light curves of three particular SLSNe, namely, SN 2005ap, SN 2006gy, and SN 2010gx, and show that they are best fitted with a wind that extends to a radius of order 10^{15} cm. These results serve as an additional indication that at least some SLSNe are powered by interaction of the ejected material with a steady wind of similar mass.