We present Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm observations of four young, eruptive star-disk systems at 0\farcs4 resolution: two FUors (V582 Aur and V900 Mon), one EXor (UZ Tau E) and one source with an ambiguous FU/EXor classification (GM Cha). The disks around GM Cha, V900 Mon and UZ Tau E are resolved. These observations increase the sample of FU/EXors observed at sub-arcsecond resolution by 15%. The disk sizes and masses of FU/EXors objects observed by ALMA so far suggest that FUor disks are more massive than Class 0/I disks in Orion and Class II disks in Lupus of similar size. EXor disks in contrast do not seem to be distinguishable from these two populations. We reach similar conclusions when comparing the FU/EXor sample to the Class I and Class II disks in Ophiuchus. FUor disks around binaries are host to more compact disks than those in single-star systems, similar to non-eruptive young disks. We detect a wide-angle outflow around GM Cha in $^{12}$CO emission, wider than typical Class I objects and more similar to those found around some FUor objects. We use radiative transfer models to fit the continuum and line data of the well-studied disk around UZ Tau E. The line data is well described by a keplerian disk, with no evidence of outflow activity (similar to other EXors). The detection of wide-angle outflows in FUors and not in EXors support to the current picture in which FUors are more likely to represent an accretion burst in the protostellar phase (Class I), while EXors are smaller accretion events in the protoplanetary (Class II) phase.

We present Adaptive Optics (AO) near infrared (NIR) observations using VLT/NACO and Keck/NIRC2 of ODISEA targets. ODISEA is an ALMA survey of the entire population of circumstellar discs in the Ophiuchus molecular cloud. From the whole sample of ODISEA we select all the discs that are not already observed in the NIR with AO and that are observable with NACO or NIRC2. The NIR-ODISEA survey consists of 147 stars observed in NIR AO imaging for the first time, as well as revisiting almost all the binary systems of Ophiuchus present in the literature (20 out of 21). In total, we detect 20 new binary systems and one triple system. For each of them we calculate the projected separation and position angle of the companion, as well as their NIR and millimeter flux ratios. From the NIR contrast we derived the masses of the secondaries, finding that 9 of them are in the sub-stellar regime (30-50 \MJup). Discs in multiple systems reach a maximum total dust mass of $\sim$ 50 M$_{\oplus}$, while discs in single stars can reach a dust mass of 200 M$_{\oplus}$. Discs with masses above 10 M$_{\oplus}$ are found only around binaries with projected separations larger than $\sim$ 110 au. The maximum disc size is also larger around single star than binaries. However, since most discs in Ophiuchus are very small and low-mass, the effect of visual binaries is relatively weak in the general disc population.

Using the Atacama Large Millimeter/submillimeter Array (ALMA), we observed the young Herbig star HD 100546, host to a prominent disk with a deep, wide gap in the dust. The high-resolution 1.3 mm continuum observation reveals fine radial and azimuthal substructures in the form of a complex maze of ridges and trenches sculpting a dust ring. The $^{12}$CO(2-1) channel maps are modulated by wiggles or kinks that deviate from Keplerian kinematics particularly over the continuum ring, where deviations span 90$^\circ$ in azimuth, covering 5 km s$^{-1}$. The most pronounced wiggle resembles the imprint of an embedded massive planet of at least 5 M$_{\rm Jup}$ predicted from previous hydrodynamical simulations (Perez, Casassus, & Benitez-Llambay 2018). Such planet is expected to open a deep gap in both gas and dust density fields within a few orbital timescales, yet the kinematic wiggles lie near ridges in the continuum. The lesser strength of the wiggles in the $^{13}$CO and C$^{18}$O isotopologues show that the kinematic signature weakens at lower disk heights, and suggests qualitatively that it is due to vertical flows in the disk surface. Within the gap, the velocity field transitions from Keplerian to strongly non-Keplerian via a twist in position angle, suggesting the presence of another perturber and/or an inner warp. We also present VLT/SPHERE sparse aperture masking data which recovers scattered light emission from the gap's edges but shows no evidence for signal within the gap, discarding a stellar binary origin for its opening.

The Art, Astronomy, Technology and Society (AATS) project is an artscience collaboration which dwells in the inbetweenness layer of scientific understanding of the origins of planets and the embodied, intuited ways of knowing. One of the oldest and at the same time one of the newest concerns of humanity is "How did the Earth and the planets come to be?". Planet formation, an intricate and potentially chaotic process, is also very efficient. Every star harbors at least one planet, as evidenced by the high frequency of exoplanet detections. Planet formation must therefore be a frequent process. Yet, learning about the origins of planets has been difficult. Here I share my experience in the AATS artscience project and discuss the nature of artscience collaborations.

The SR24 multi-star system hosts both circumprimary and circumsecondary disks, which are strongly misaligned from each other. The circumsecondary disk is circumbinary in nature. Interestingly, both disks are interacting, and they possibly rotate in opposite directions. To investigate the nature of this unique twin disk system, we present 0.''1 resolution near-infrared polarized intensity images of the circumstellar structures around SR24, obtained with HiCIAO mounted on the Subaru 8.2 m telescope. Both the circumprimary disk and the circumsecondary disk are resolved and have elongated features. While the position angle of the major axis and radius of the NIR polarization disk around SR24S are 55$^{\circ}$ and 137 au, respectively, those around SR24N are 110$^{\circ}$ and 34 au, respectively. With regard to overall morphology, the circumprimary disk around SR24S shows strong asymmetry, whereas the circumsecondary disk around SR24N shows relatively strong symmetry. Our NIR observations confirm the previous claim that the circumprimary and circumsecondary disks are misaligned from each other. Both the circumprimary and circumsecondary disks show similar structures in $^{12}$CO observations in terms of its size and elongation direction. This consistency is because both NIR and $^{12}$CO are tracing surface layers of the flared disks. As the radius of the polarization disk around SR24N is roughly consistent with the size of the outer Roche lobe, it is natural to interpret the polarization disk around SR24N as a circumbinary disk surrounding the SR24Nb-Nc system.

FU Orionis objects are low-mass pre-main sequence stars characterized by dramatic outbursts of several magnitudes in brightness. These outbursts are linked to episodic accretion events in which stars gain a significant portion of their mass. The physical processes behind these accretion events are not yet well understood. The archetypical FU Ori system, FU Orionis, is composed of two young stars with detected gas and dust emission. The continuum emitting regions have not been resolved until now. Here, we present 1.3 mm observations of the FU Ori binary system with ALMA. The disks are resolved at 40 mas resolution. Radiative transfer modeling shows that the emission from FU Ori north (primary) is consistent with a dust disk with a characteristic radius of $\sim$11 au. The ratio between major and minor axes shows that the inclination of the disk is $\sim$37 deg. FU Ori south is consistent with a dust disk of similar inclination and size. Assuming the binary orbit shares the same inclination angle as the disks, the deprojected distance between north and south components is 0.6'', i.e. $\sim$250 au. Maps of $^{12}$CO emission show a complex kinematic environment with signatures disk rotation at the location of the northern component, and also (to a lesser extent) for FU Ori south. The revised disk geometry allows us to update FU Ori accretion models (Zhu et al.), yielding a stellar mass and mass accretion rate of FU Ori north of 0.6 M$_{\odot}$ and 3.8$\times10^{-5}$ M$_{\odot}$ yr$^{-1}$, respectively.

The migration of planetary cores embedded in a protoplanetary disk is an important mechanism within planet-formation theory, relevant for the architecture of planetary systems. Consequently, planet migration is actively discussed, yet often results of independent theoretical or numerical studies are unconstrained due to the lack of observational diagnostics designed in light of planet migration. In this work we follow the idea of inferring the migration behavior of embedded planets by means of the characteristic radial structures that they imprint in the disk's dust density distribution. We run hydrodynamical multifluid simulations of gas and several dust species in a locally isothermal $\alpha$-disk in the low-viscosity regime ($\alpha=10^{-5}$) and investigate the obtained dust structures. In this framework, a planet of roughly Neptune mass can create three (or more) rings in which dust accumulates. We find that the relative spacing of these rings depends on the planet's migration speed and direction. By performing subsequent radiative transfer calculations and image synthesis we show that - always under the condition of a near-inviscid disk - different migration scenarios are, in principle, distinguishable by long-baseline, state-of-the-art ALMA observations.

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 Solar System gas giants are each surrounded by many moons, with at least 50 prograde satellites thought to have formed from circumplanetary material. Just like the Sun is not the only star surrounded by planets, extrasolar gas giants are likely surrounded by satellite systems. Here, we report on ALMA observations of four <40 Myr old stars with directly-imaged companions: PZ Tel, AB Pic, 51 Eri, and $\kappa$ And. Continuum emission at 1.3 mm is undetected for any of the systems. Since these are directly-imaged companions, there is knowledge of their temperatures, masses and locations. These allow for upper limits on the amount of circumplanetary dust to be derived from detailed radiative transfer models. These protolunar disc models consider two disc sizes: 0.4 and 0.04 times the exoplanet's Hill radius. The former is representative of hydrodynamic simulations of circumplanetary discs while the latter a case with significant radial drift of solids. The more compact case is also motivated by the semi-major axis of Callisto, enclosing Jupiter's Galilean satellites. All upper limits fall below the expected amount of dust required to explain regular satellite systems ($10^{-4}$ times the mass of their central planet). Upper limits are compared with viscous evolution and debris disc models. Our analysis suggests that the non detections can be interpreted as evidence of dust growth beyond metre sizes to form moonetesimals in timescales <10 Myr. This sample increases by 50% the number of ALMA non-detections of young companions available in the literature.

Rings and radial gaps are ubiquitous in protoplanetary disks, yet their possible connection to planet formation is currently subject to intense debates. In principle, giant planet formation leads to wide gaps which separate the gas and dust mass reservoir in the outer disk, while lower mass planets lead to shallow gaps which are manifested mainly on the dust component. We used the Atacama Large Millimeter/submillimeter Array (ALMA) to observe the star HD169142, host to a prominent disk with deep wide gaps that sever the disk into inner and outer regions. The new ALMA high resolution images allow for the outer ring to be resolved as three narrow rings. The HD169142 disk thus hosts both the wide gaps trait of transition disks and a narrow ring system similar to those observed in the TW Hya and HL Tau systems. The mass reservoir beyond a deep gap can thus host ring systems. The observed rings are narrow in radial extent (width/radius of 1.5/57.3, 1.8/64.2 and 3.4/76.0, in au) and have asymmetric mutual separations: the first and middle ring are separated by 7 au while the middle and outermost ring are distanced by ~12 au. Using hydrodynamical modeling we found that a simple explanation, involving a single migrating low mass planet (10 M$_\oplus$), entirely accounts for such an apparently complex phenomenon. Inward migration of the planet naturally explains the ring's asymmetric mutual separation. The isolation of HD169142's outer rings thus allows a proof of concept to interpret the detailed architecture of the outer region of protoplanetary disks with low mass planet formation of mini-Neptune's size, i.e. as in the protosolar nebula.

Resolved ALMA and VLA observations indicate the existence of two dust traps in the protoplanetary disc MWC 758. By means of 2D gas+dust hydrodynamical simulations post-processed with 3D dust radiative transfer calculations, we show that the spirals in scattered light, the eccentric, asymmetric ring and the crescent-shaped structure in the (sub)millimetre can all be caused by two giant planets: a 1.5-Jupiter mass planet at 35 au (inside the spirals) and a 5-Jupiter mass planet at 140 au (outside the spirals). The outer planet forms a dust-trapping vortex at the inner edge of its gap (at ~85 au), and the continuum emission of this dust trap reproduces the ALMA and VLA observations well. The outer planet triggers several spiral arms which are similar to those observed in polarised scattered light. The inner planet also forms a vortex at the outer edge of its gap (at ~50 au), but it decays faster than the vortex induced by the outer planet, as a result of the disc's turbulent viscosity. The vortex decay can explain the eccentric inner ring seen with ALMA as well as the low signal and larger azimuthal spread of this dust trap in VLA observations. Finding the thermal and kinematic signatures of both giant planets could verify the proposed scenario.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

The effect of magnetic fields in the Equations of State (EoS) of compact objects is the splitting of the pressure in two components, one parallel and the other perpendicular to the magnetic field. This anisotropy suggests the necessity of using structure equations considering the axial symmetry of the magnetized system. In this work, we consider an axially symmetric metric in spherical coordinates, the $\gamma$-metric, and construct a system of equations to describe the structure of spheroidal compact objects. In addition, we connect the geometrical parameter $\gamma$ linked to the spheroid's radii, with the source of the anisotropy. So, the model relates the shape of the compact object to the physics that determines the properties of the composing matter. To illustrate how our structure equations work, we obtain the mass-radii solutions for magnetized White Dwarfs. Our results show that the main effect of the magnetic field anisotropy in White Dwarfs structure is to cause a deformation of these objects. Since this effect is only relevant at low densities, it does not affect the maximum values of magnetized White Dwarf's masses, which remain under Chandrasekhar limit.

We introduce the Ophiuchus DIsc Survey Employing ALMA (ODISEA), a project aiming to study the entire population of Spitzer-selected protoplanetary discs in the Ophiuchus Molecular Cloud (~300 objects) from both millimeter continuum and CO isotopologues data. Here we present 1.3 mm/230 GHz continuum images of 147 targets at 0.2" (28 au) resolution and a typical rms of 0.15 mJy. We detect a total of 133 discs, including the individual components of 11 binary systems and 1 triple system. Fifty-three of these discs are spatially resolved. We find clear substructures (inner cavities, rings, gaps, and/or spiral arms) in 8 of the sources and hints of such structures in another 4 discs. We construct the disc luminosity function for our targets and perform comparisons to other regions. A simple conversion between flux and dust mass (adopting standard assumptions) indicates that all discs detected at 1.3 mm are massive enough to form one or more rocky planets. In contrast, only ~50 discs (~1/3 of the sample) have enough mass in the form of dust to form the canonical 10 M_Earth core needed to trigger runaway gas accretion and the formation of gas giant planets, although the total mass of solids already incorporated into bodies larger than cm scales is mostly unconstrained. The distribution in continuum disc sizes in our sample is heavily weighted towards compact discs: most detected discs have radii < 15 au, while only 23 discs (~15% of the targets) have radii > 30 au.

We investigate the rest-frame UV morphologies of a large sample of Lyman-a emitters (LAEs) from z~2 to z~6, selected in a uniform way with 16 different narrow- and medium-bands over the full COSMOS field. We use 3045 LAEs with HST coverage in a stacking analysis and find that they have M_UV~-20, below M*_UV at these redshifts. We also focus our analysis on a subsample of 780 individual galaxies with i_AB<25 for which GALFIT converges for 429 of them. The individual median size (re~1 kpc), ellipticities (slightly elongated with (b/a)~0.45), S\'ersic index (disk-like with n<2) and light concentration (comparable to that of disk or irregular galaxies, with C~2.7) of LAEs show mild evolution from z~2 to z~6. LAEs with the highest rest-frame equivalent widths (EW) are the smallest/most compact (re~0.8 kpc, compared to re~1.5 kpc for the lower EW LAEs). When stacking our samples in bins of fixed Lya luminosity and Lya EW we find evidence for redshift evolution in n and C, but not in galaxy sizes. The evolution seems to be stronger for LAEs with 25<EW<100 {\AA}. When compared to other SFGs, LAEs are found to be smaller at all redshifts. The difference between the two populations changes with redshift, from a factor of ~1 at z>5 to SFGs being a factor of ~2-4 larger than LAEs for z<2. This means that at the highest redshifts, where typical sizes approach those of LAEs, the fraction of galaxies showing Lya in emission (and with a high Lya escape fraction) should be much higher, consistent with observations.

CO observations allow estimations of the gas content of molecular clouds, which trace the reservoir of cold gas fuelling star formation, as well as to determine extinction via H$_2$ column density, N(H$_2$). Here, we studied millimetric and optical properties at 26 supernovae (SNe) locations of different types in a sample of 23 nearby galaxies by combining molecular $^{12}$C$^{16}$O (J = 1 $\rightarrow$ 0) resolved maps from the EDGE survey and optical Integral Field Spectroscopy from the CALIFA survey. We found an even clearer separation between type II and type Ibc SNe in terms of molecular gas than what we found in the optical using H$\alpha$ emission as a proxy for current SF rate, which reinforces the fact that SNe Ibc are more associated with SF-environments. While A$_V$ at SN locations is similar for SNe II and SNe Ibc, and higher compared to SNe Ia, N(H$_2$) is significantly higher for SNe Ibc than for SNe II and SNe Ia. When compared to alternative extinction estimations directly from SN photometry and spectroscopy, we find that our SNe Ibc have also redder color excess but showed standard Na I D absorption pseudo-equivalent widths ($\sim$1 \AA). In some cases we find no extinction when estimated from the environment, but high amounts of extinction when measured from SN observations, which suggests that circumstellar material or dust sublimation may be playing a role. This work serves as a benchmark for future studies combining last generation millimeter and optical IFS instruments to reveal the local environmental properties of extragalactic SNe.

We carried out a 12CO(3-2) survey of 52 southern stars with a wide range of IR excesses (LIR/L*) using the single dish telescopes APEX and ASTE. The main aims were (1) to characterize the evolution of molecular gas in circumstellar disks using LIR/L* values as a proxy of disk dust evolution, and (2) to identify new gas-rich disk systems suitable for detailed study with ALMA. About 60% of the sample (31 systems) have LIR/L* > 0.01 typical of T-Tauri or Herbig AeBe stars, and the rest (21 systems) have LIR/L* < 0.01 typical of debris disks. We detect CO(3-2) emission from 20 systems, and 18 (90%) of these have LIR/L* > 0.01. However, the spectra of only four of the newly detected systems appear free of contamination from background or foreground emission from molecular clouds. These include the early-type stars HD 104237 (A4/5V, 116 pc) and HD 98922 (A2 III, 507 pc, as determined in this work), where our observations reveal the presence of CO-rich circumstellar disks for the first time. Of the other detected sources, many could harbor gaseous circumstellar disks, but our data are inconclusive. For these two newly discovered gas-rich disks, we present radiative transfer models that simultaneously reproduce their spectral energy distributions and the 12CO(3-2) line profiles. For both of these systems, the data are fit well by geometrically flat disks, placing them in the small class of non-flaring disks with significant molecular gas reservoirs.