Theoretical models of the ionization state in protoplanetary disks suggest the existence of large areas with low ionization and weak coupling between the gas and magnetic fields. In this regime hydrodynamical instabilities may become important. In this work we investigate the gas and dust structure and dynamics for a typical T Tauri system under the influence of the vertical shear instability (VSI). We use global 3D radiation hydrodynamics simulations covering all $360^\circ$ of azimuth with embedded particles of 0.1 and 1mm size, evolved for 400 orbits. Stellar irradiation heating is included with opacities for 0.1- to 10-$\mu$m-sized dust. Saturated VSI turbulence produces a stress-to-pressure ratio of $\alpha \simeq 10^{-4}$. The value of $\alpha$ is lowest within 30~au of the star, where thermal relaxation is slower relative to the orbital period and approaches the rate below which VSI is cut off. The rise in $\alpha$ from 20 to 30~au causes a dip in the surface density near 35~au, leading to Rossby wave instability and the generation of a stationary, long-lived vortex spanning about 4~au in radius and 40~au in azimuth. Our results confirm previous findings that mm size grains are strongly vertically mixed by the VSI. The scale height aspect ratio for 1mm grains is determined to be 0.037, much higher than the value $H/r=0.007$ obtained from millimeter-wave observations of the HL~Tau system. The measured aspect ratio is better fit by non-ideal MHD models. In our VSI turbulence model, the mm grains drift radially inwards and many are trapped and concentrated inside the vortex. The turbulence induces a velocity dispersion of $\sim 12$~m/s for the mm grains, indicating that grain-grain collisions could lead to fragmentation.

We present a numerical code for radiation hydrodynamics designed as a module for the freely available PLUTO code. We adopt a gray approximation and include radiative transfer following a two-moment approach by imposing the M1 closure to the radiation fields. This closure allows for a description of radiative transport in both the diffusion and free-streaming limits, and is able to describe highly anisotropic radiation transport as can be expected in the vicinity of an accreting planet in a protoplanetary disk. To reduce the computational cost caused by the timescale disparity between radiation and matter fields, we integrate their evolution equations separately in an operator-split way, using substepping to evolve the radiation equations. We further increase the code's efficiency by adopting the reduced speed of light approximation (RSLA). Our integration scheme for the evolution equations of radiation fields relies on implicit-explicit schemes, in which radiation-matter interaction terms are integrated implicitly while fluxes are integrated via Godunov-type solvers. The module is suitable for general astrophysical computations in 1, 2, and 3 dimensions in Cartesian, spherical and cylindrical coordinates, and can be implemented on rotating frames. We demonstrate the algorithm performance on different numerical benchmarks, paying particular attention to the applicability of the RSLA for computations of physical processes in protoplanetary disks. We show 2D simulations of vertical convection in disks and 3D simulations of gas accretion by planetary cores, which are the first of their kind to be solved with a two-moment approach.

Stars form as an end product of the gravitational collapse of cold, dense gas in magnetized molecular clouds. This multi-scale scenario occurs via the formation of two quasi-hydrostatic cores and involves complex physical processes, which require a robust, self-consistent numerical treatment. The aim of this study is to understand the formation and evolution of the second Larson core and the dependence of its properties on the initial cloud core mass. We used the PLUTO code to perform high resolution, 1D and 2D RHD collapse simulations. We include self-gravity and use a grey FLD approximation for the radiative transfer. Additionally, we use for the gas EOS density- and temperature-dependent thermodynamic quantities to account for the effects such as dissociation, ionisation, and molecular vibrations and rotations. Properties of the second core are investigated using 1D studies spanning a wide range of initial cloud core masses from 0.5 to 100 $M_{\odot}$. Furthermore, we expand to 2D collapse simulations for a few cases of 1, 5, 10, and 20 $M_{\odot}$. We follow the evolution of the second core for $\geq$ 100 years after its formation, for each of these non-rotating cases. Our results indicate a dependence of several second core properties on the initial cloud core mass. For the first time, due to an unprecedented resolution, our 2D non-rotating collapse studies indicate that convection is generated in the outer layers of the second core, which is formed due to the gravitational collapse of a 1 $M_{\odot}$ cloud core. Additionally, we find large-scale oscillations of the second accretion shock front triggered by the standing accretion shock instability, which has not been seen before in early evolutionary stages of stars. We predict that the physics within the second core would not be significantly influenced by the effects of magnetic fields or an initial cloud rotation.

The presence of CO gas around 10-50 Myr old A stars with debris discs has sparked debate on whether the gas is primordial or secondary. Since secondary gas released from planetesimals is poor in H$_2$, it was thought that CO would quickly photodissociate never reaching the high levels observed around the majority of A stars with bright debris discs. Kral et al. 2019 showed that neutral carbon produced by CO photodissociation can effectively shield CO and potentially explain the high CO masses around 9 A stars with bright debris discs. Here we present a new model that simulates the gas viscous evolution, accounting for carbon shielding and how the gas release rate decreases with time as the planetesimal disc loses mass. We find that the present gas mass in a system is highly dependant on its evolutionary path. Since gas is lost on long timescales, it can retain a memory of the initial disc mass. Moreover, we find that gas levels can be out of equilibrium and quickly evolving from a shielded onto an unshielded state. With this model, we build the first population synthesis of gas around A stars, which we use to constrain the disc viscosity. We find a good match with a high viscosity ($\alpha\sim0.1$), indicating that gas is lost on timescales $\sim1-10$ Myr. Moreover, our model also shows that high CO masses are not expected around FGK stars since their planetesimal discs are born with lower masses, explaining why shielded discs are only found around A stars. Finally, we hypothesise that the observed carbon cavities could be due to radiation pressure or accreting planets.

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 planet formation requires to discern how dust grows in protoplanetary disks. An important parameter to measure in disks is the maximum dust grain size present. This is usually estimated through measurements of the dust opacity at different millimeter wavelengths assuming optically thin emission and dust opacity dominated by absorption. However, ALMA observations have shown that these assumptions might not be correct in the case of protoplanetary disks, leading to overestimation of particle sizes and to underestimation of the disk's mass. Here, we present an analysis of high quality ALMA and VLA images of the HL Tau protoplanetary disk, covering a wide range of wavelengths, from 0.8 mm to 1 cm, and with a physical resolution of $\sim$7.35 au. We describe a procedure to analyze a set of millimeter images without any assumption about the optical depth of the emission, and including the effects of absorption and scattering in the dust opacity. This procedure allows us to obtain the dust temperature, the dust surface density and the maximum particle size at each radius. In the HL Tau disk, we found that particles have already grown up to a few millimeters in size. We detect differences in the dust properties between dark and bright rings, with dark rings containing low dust density and small dust particles. Different features in the HL Tau disk seem to have different origins. Planet-disk interactions can explain substructure at the external half of the disk, but the internal rings seem to be associated to the presence of snow lines of several molecules.

The implementation of a new particle module describing the physics of dust grains coupled to the gas via drag forces is the subject of this work. The proposed particle-gas hybrid scheme has been designed to work in Cartesian as well as in cylindrical and spherical geometries. The numerical method relies on a Godunov-type second-order scheme for the fluid and an exponential midpoint rule for dust particles which overcomes the stiffness introduced by the linear coupling term. Besides being time-reversible and globally second-order accurate in time, the exponential integrator provides energy errors which are always bounded and it remains stable in the limit of arbitrarily small particle stopping times yielding the correct asymptotic solution. Such properties make this method preferable to the more widely used semi-implicit or fully implicit schemes at a very modest increase in computational cost. Coupling between particles and grid quantities is achieved through particle deposition and field-weighting techniques borrowed from Particle-In-Cell simulation methods. In this respect, we derive new weight factors in curvilinear coordinates that are more accurate than traditional volume- or area-weighting. A comprehensive suite of numerical benchmarks is presented to assess the accuracy and robustness of the algorithm in Cartesian, cylindrical and spherical coordinates. Particular attention is devoted to the streaming instability which is analyzed in both local and global disk models. The module is part of the PLUTO code for astrophysical gas-dynamics and it is mainly intended for the numerical modeling of protoplanetary disks in which solid and gas interact via aerodynamic drag.

We present a detailed multi-wavelength characterization of the multi-ring disk of HD 169142. We report new ALMA observations at 3 mm and analyze them together with archival 0.89 and 1.3 mm data. Our observations resolve three out of the four rings in the disk previously seen in high-resolution ALMA data. A simple parametric model is used to estimate the radial profile of the dust optical depth, temperature, density, and particle size distribution. We find that the multiple ring features of the disk are produced by annular accumulations of large particles, probably associated with gas pressure bumps. Our model indicates that the maximum dust grain size in the rings is $\sim1$ cm, with slightly flatter power-law size distributions than the ISM-like size distribution ($p\sim3.5$) found in the gaps. In particular, the inner ring ($\sim26$ au) is associated with a strong and narrow buildup of dust particles that could harbor the necessary conditions to trigger the streaming instability. According to our analysis, the snowlines of the most important volatiles do not coincide with the observed substructures. We explore different ring formation mechanisms and find that planet-disk interactions are the most likely scenario to explain the main features of HD 169142. Overall, our multi-wavelength analysis provides some of the first unambiguous evidence of the presence of radial dust traps in the rings of HD 169142. A similar analysis in a larger sample of disks could provide key insights on the impact that disk substructures have on the dust evolution and planet formation processes.

ALMA observations revealed recently polarised radiation of several protoplanetary disks in the (sub-)millimetre wavelength range. Besides self-scattering of large particles, thermal emission by elongated grains is a potential source for the detected polarisation signal. We calculate the wavelength dependent absorption and intrinsic polarisation of spheroidally shaped, micrometre and sub-millimetre sized dust grains using the discrete dipole approximation. In particular, we analyse the impact of dust grain porosity which appears to be present in disks when small grains coagulate to form larger aggregates. For the first time our results show that (a) the intrinsic polarisation decreases for increasing grain porosity and (b) the polarisation orientation flips by 90 degree for certain ratios of wavelength to grain size. We present a new method to constrain grain porosity and the grain size in protoplanetary disks using multi-wavelength polarisation observations in the far-infrared to millimetre wavelengths. Finally, we find that moderate grain porosities ($\mathcal{P}\lesssim0.7$) potentially explain the observed polarisation fraction in the system HD 142527 while highly porous grains ($\mathcal{P}>0.7$) fail unless the grain's axis ratio is extraordinarily large.

PDS 70b is the most robustly detected young planet imaged in the gap of a transition disk so far, found at a separation of ~195 mas (~22 au) from its host star and at a position angle of about 155 deg. This system is therefore a unique laboratory to characterize the properties of young planetary systems at the stage of their formation. We aim to trace direct and indirect imprints of PDS 70b on the gas and dust emission of the circumstellar disk in order to study the properties of this ~5 Myr young planetary system. We obtained ALMA band 7 observations of PDS 70 in dust continuum and $^{12}CO$ (3-2) and combined them with archival data resulting in an unprecedented angular resolution of about 70 mas (~8 au). We derive an upper limit on circumplanetary material at the location of PDS 70b of ~0.01 $M_{\oplus}$ and find a highly structured circumstellar disk in both dust and gas. The outer dust ring peaks at ~0.65" (74 au) and reveals a possible second unresolved peak at ~0.53" (60 au). The integrated intensity of CO also shows evidence of a depletion of emission at ~0.2" (23 au) with a width of ~0.1" (11 au). The gas kinematics show evidence of a deviation from Keplerian rotation inside $\lesssim$0.8" (91 au). This implies a pressure gradient that can account for the location of the dust ring well beyond the location of PDS 70b. Further in, we detect an inner disk which appears to be connected to the outer disk by a possible bridge feature in the North-West region in both gas and dust. We compare the observations to hydrodynamical simulations including a planet with different masses covering the estimated mass range previously derived from NIR photometry (~5-9 M$_{Jup}$). We find that even a planet with a mass of 10 M$_{Jup}$ may not be sufficient to explain the extent of the wide gap and an additional low-mass companion may be needed to account for the observed disk morphology.

We study dust concentration in axisymmetric gas rings in protoplanetary disks. Given the gas surface density, we derived an analytical total dust surface density by taking into account the differential concentration of all the grain sizes. This model allows us to predict the local dust-to-gas mass ratio and the slope of the particle size distribution, as a function of radius. We test this analytical model comparing it with a 3D magneto-hydrodynamical simulation of dust evolution in an accretion disk. The model is also applied to the disk around HD 169142. By fitting the disk continuum observations simultaneously at $\lambda = 0.87$, 1.3, 3.0 mm, we obtain a global dust-to-gas mass ratio $\epsilon_{\rm global} = 1.05 \times 10^{-2}$ and a viscosity coefficient $\alpha = 1.35 \times 10^{-2}$. This model can be easily implemented in numerical simulations of accretion disks.

Context. The mechanisms governing the opening of cavities in transition disks are not fully understood. Several processes have been proposed but their occurrence rate is still unknown. Aims. We present spatially resolved observations of two transition disks and aim at constraining their vertical and radial structure using multiwavelength observations. Methods. We have obtained near-IR scattered light observations with VLT/SPHERE of the transition disks J1608 and J1852. We complement our datasets with ALMA observations and with unresolved photometric observations covering a wide range of wavelengths. We performed radiative transfer modeling to analyze the morphology of the disks and compare the results with a sample of 20 other transition disks observed with both SPHERE and ALMA. Results. The scattered light image of J1608 reveals a very inclined disk, with two bright lobes and a large cavity. J1852 shows an inner ring extending beyond the coronagraphic radius up to 15au, a gap and a second ring at 42au. Our radiative transfer model of J1608 indicates that the millimeter-sized grains are less extended vertically and radially than the micron-sized grains, indicating advanced settling and radial drift. We find good agreement with the observations of J1852 with a similar model, but due to the low inclination of the system, the model remains partly degenerate. The analysis of 22 transition disks shows that, in general, the cavities observed in scattered light are smaller than the ones detected at millimeter wavelengths. Conclusions. The analysis of a sample of transition disks indicates that the small grains can flow inward of the region where millimeter grains are trapped. While 15 out of the 22 cavities in our sample could be explained by a planet of less than 13 Jupiter masses, the others either require the presence of a more massive companion or of several low-mass planets.

We present simulations of the capabilities of the Next Generation Very Large Array to image at high angular resolution substructures in the dust emission of protoplanetary disks. The main goal of this study is to investigate the kinds of substructures that are expected by state-of-the-art 3D simulations of disks and that an instrument like the ngVLA, with its current design, can detect. The disk simulations adopted in this investigation consist of global 3D radiation-hydrodynamics models with embedded particles, the latter representing dust grains. Our work shows that the ngVLA can detect and spatially resolve, down to sub-astronomical unit scales in disks in nearby star forming regions, the dust continuum emission at 3mm from azimuthal asymmetric structures, as well as from weak rings and gaps produced in these models as a consequence of the vertical shear instability (VSI). This hydrodynamical instability has been proposed to generate turbulence in regions of weak coupling between the disk gas and magnetic field, as well as to form vortices which may be preferred locations of planetesimal formation.

We perform simulations of the dust and gas disk evolution to investigate the observational features of a dust-pileup at the dead-zone inner edge. We show that the total mass of accumulated dust particles is sensitive to the turbulence strength in the dead zone, $\alpha_{\rm dead}$, because of the combined effect of turbulence-induced particle fragmentation (which suppresses particle radial drift) and turbulent diffusion. For a typical critical fragmentation velocity of silicate dust particles of $1~{\rm m~s^{-1}}$, the stress to pressure ratio $\alpha_{\rm dead}$ needs to be lower than $3 \times 10^{-4}$ for dust trapping to operate. The obtained dust distribution is postprocessed using the radiative transfer code RADMC-3D to simulate infrared scattered-light images of the inner part of protoplanetary disks with a dust pileup. We find that a dust pileup at the dead-zone inner edge, if present, casts a shadow extending out to $\sim 10~{\rm au}$. In the shadowed region the temperature significantly drops, which in some cases yields even multiple water snow lines. We also find that even without a dust pileup at the dead-zone inner edge, the disk surface can become thermally unstable, and the excited waves can naturally produce shadows and ring-like structures in observed images. This mechanism might account for the ring-like structures seen in the scattered-light images of some disks, such as the TW Hya disk.

Dwarf novae (DNe) are accreting white dwarfs that show eruptions due to a thermal-viscous instability in the accretion disk. The outburst timescales constrain $\alpha$, the ratio of the viscous stress to the thermal pressure, and so the mechanism of angular momentum transport. The eruptive state has $\alpha\approx0.1$ while the quiescent state has $\alpha\approx0.03$. Turbulent transport due to the magneto-rotational instability (MRI) is generally considered to be the source of angular momentum transport in DNe. Here, we perform 3D local magnetohydrodynamic (MHD) shearing box simulations including vertical stratification, radiative transfer and a net constant vertical magnetic flux to investigate how transport changes between the outburst and quiescent states of DNe. We find that a constant $B_z$ provides a higher $\alpha$ in quiescence than in outburst, in opposition to what is expected. Including resistivity quenches MRI turbulence in quiescence, suppressing transport, unless the magnetic field is high enough, which again leads to $\alpha\approx0.1$. A major difference between simulations with a net poloidal flux and simulations without is that angular momentum transport in the former is shared between turbulent and wind-driven transport. We find that wind-driven transport dominates in quiescence even for low magnetic fields $\sim 1$ G. This can have a major impact on observational signatures since wind-driven transport does not heat the disk. Furthermore, wind transport cannot be reduced to an $\alpha$ prescription. We provide fits for $\alpha$ and the wind torque with $\beta$, ratio of thermal to magnetic pressure. We conclude that the evolution of the thermal-viscous instability, and its consequences on the outburst cycles of CVs, needs to be seriously revised to take into account that most of the accretion energy may be carried away by a wind instead of being locally dissipated.

We present a survey of optical [OI] emission at 6300 \AA\ toward 65 TTauri stars at the spectral resolution of $\sim7$ km/s. Past work identified a highly blueshifted velocity component (HVC) tracing microjets, and a less blueshifted low-velocity-component (LVC) attributed to winds. We focus here on the LVC kinematics to investigate links between winds, jets, accretion, and disk dispersal. We track the behavior of four types of LVC components: a broad and narrow component ("BC" and "NC") in LVCs that are decomposed into two Gaussians, which typically have an HVC, and the single-Gaussian LVC profiles separated into those that have an HVC ("SCJ") and those that do not ("SC"). LVC centroid velocities and line widths correlate with HVC equivalent width and accretion luminosity, suggesting that LVC/winds and HVC/jets are kinematically linked and connected to accretion. The deprojected HVC velocity correlates with accretion luminosity, showing that faster jets come with higher accretion. BC and NC kinematics correlate and their blueshifts are maximum at $\sim 35^{\circ}$, suggesting a conical wind geometry with this semi-opening angle. Only SCs include $n_{13-31}$ up to $\sim3$ and their properties correlate with this infrared index, showing that [OI] emission recedes to larger radii as the inner dust is depleted, tracing less dense/hot gas and a decrease in wind velocity. All together, these findings support a scenario where optically thick, accreting inner disks launch radially-extended MHD disk winds that feed jets, and where inner disk winds recede to larger radii and jets disappear in concert with dust depletion.

Magnetohydrodynamic (MHD) and photoevaporative winds are thought to play an important role in the evolution and dispersal of planet-forming disks. We report the first high-resolution ($\Delta v\sim$6\kms) analysis of [S II] $\lambda$4068, [O I] $\lambda$5577, and [O I] $\lambda$6300 lines from a sample of 48 T Tauri stars. Following Simon et al. (2016), we decompose them into three kinematic components: a high-velocity component (HVC) associated with jets, and a low-velocity narrow (LVC-NC) and broad (LVC-BC) components. We confirm previous findings that many LVCs are blueshifted by more than 1.5 kms$^{-1}$ thus most likely trace a slow disk wind. We further show that the profiles of individual components are similar in the three lines. We find that most LVC-BC and NC line ratios are explained by thermally excited gas with temperatures between 5,000$-$10,000 K and electron densities $\sim10^{7}-10^{8}$ cm$^{-3}$. The HVC ratios are better reproduced by shock models with a pre-shock H number density of $\sim10^{6}-10^{7}$ cm$^{-3}$. Using these physical properties, we estimate $\dot{M}_{\rm wind}/\dot{M}_{\rm acc}$ for the LVC and $\dot{M}_{\rm jet}/\dot{M}_{\rm acc}$ for the HVC. In agreement with previous work, the mass carried out in jets is modest compared to the accretion rate. With the likely assumption that the NC wind height is larger than the BC, the LVC-BC $\dot{M}_{\rm wind}/\dot{M}_{\rm acc}$ is found to be higher than the LVC-NC. These results suggest that most of the mass loss occurs close to the central star, within a few au, through an MHD driven wind. Depending on the wind height, MHD winds might play a major role in the evolution of the disk mass.

Observations of multiple rotational transitions from a single molecule allow for unparalleled constraints on the physical conditions of the emitting region. We present an analysis of CS in TW~Hya using the $J=7-6$, $5-4$ and $3-2$ transitions imaged at $\sim 0.5^{\prime\prime}$ spatial resolution, resulting in a temperature and column density profile of the CS emission region extending out to 230~au, far beyond previous measurements. In addition, the 15~kHz resolution of the observations and the ability to directly estimate the temperature of the CS emitting gas, allow for one of the most sensitive searches for turbulent broadening in a disk to date. Limits of $v_{\rm turb} \lesssim 0.1 c_s$ can be placed across the entire radius of the disk. We are able to place strict limits of the local H$_2$ density due to the collisional excitations of the observed transitions. From these we find that a minimum disk mass of $3 \times 10^{-4}~M_{\rm sun}$ is required to be consistent with the CS excitation conditions and can uniquely constrain the gas surface density profile in the outer disk.

Young circumstellar disks are of prime interest to understand the physical and chemical conditions under which planet formation takes place. Only very few detections of planet candidates within these disks exist, and most of them are currently suspected to be disk features. In this context, the transition disk around the young star PDS 70 is of particular interest, due to its large gap identified in previous observations, indicative of ongoing planet formation. We aim to search for the presence of planets and search for disk structures indicative for disk-planet interactions and other evolutionary processes. We analyse new and archival near-infrared (NIR) images of the transition disk PDS 70 obtained with the VLT/SPHERE, VLT/NaCo and Gemini/NICI instruments in polarimetric differential imaging (PDI) and angular differential imaging (ADI) modes. We detect a point source within the gap of the disk at about 195 mas (about 22 au) projected separation. The detection is confirmed at five different epochs, in three filter bands and using different instruments. The astrometry results in an object of bound nature, with high significance. The comparison of the measured magnitudes and colours to evolutionary tracks suggests that the detection is a companion of planetary mass. We confirm the detection of a large gap of about 54 au in size within the disk in our scattered light images, and detect a signal from an inner disk component. We find that its spatial extent is very likely smaller than about 17 au in radius. The images of the outer disk show evidence of a complex azimuthal brightness distribution which may in part be explained by Rayleigh scattering from very small grains. Future observations of this system at different wavelengths and continuing astrometry will allow us to test theoretical predictions regarding planet-disk interactions, planetary atmospheres and evolutionary models.

In the present study we aim to investigate the circumstellar environment of the spectroscopic binary T Tauri star CS Cha. From unresolved mid- to far-infrared photometry it is predicted that CS Cha hosts a disk with a large cavity. In addition, SED modeling suggests significant dust settling, pointing towards an evolved disk that may show signs of ongoing or completed planet formation. We observed CS Cha with the high contrast imager VLT/SPHERE in polarimetric differential imaging mode to resolve the circumbinary disk in near infrared scattered light. These observations were followed-up by VLT/NACO L-band observations and complemented by archival VLT/NACO K-band and HST/WFPC2 I-band data. We resolve the compact circumbinary disk around CS Cha for the first time in scattered light. We find a smooth, low inclination disk with an outer radius of $\sim$55 au (at 165 pc). We do not detect the inner cavity but find an upper limit for the cavity size of $\sim$15 au. Furthermore, we find a faint co-moving companion with a projected separation of 210 au from the central binary outside of the circumbinary disk. The companion is detected in polarized light and shows an extreme degree of polarization (13.7$\pm$0.4 \% in J-band). The companion's J- and H-band magnitudes are compatible with masses of a few M$_\mathrm{Jup}$. However, K-, L- and I-band data draw this conclusion into question. We explore with radiative transfer modeling whether an unresolved circum-companion disk can be responsible for the high polarization and complex photometry. We find that the set of observations is best explained by a heavily extincted low mass ($\sim 20 \mathrm{M}_\mathrm{Jup}$) brown dwarf or high mass planet with an unresolved disk and dust envelope.

The discovery of thousands of exoplanets over the last couple of decades has shown that the birth of planets is a very efficient process in nature. Theories invoke a multitude of mechanisms to describe the assembly of planets in the disks around pre-main-sequence stars, but observational constraints have been sparse on account of insufficient sensitivity and resolution. Understanding how planets form and interact with their parental disk is crucial also to illuminate the main characteristics of a large portion of the full population of planets that is inaccessible to current and near-future observations. This White Paper describes some of the main issues for our current understanding of the formation and evolution of planets, and the critical contribution expected in this field by the Next Generation Very Large Array.

The eruptive cycles of dwarf novae (DN) are thought to be due to a thermal-viscous instability in the accretion disk surrounding the white dwarf (WD). This model has long been known to imply a stress to pressure ratio \alpha ~0.1 in outburst compared to \alpha ~ 0.01 in quiescence. Such an enhancement in $\alpha$ has recently been observed in simulations of turbulent transport driven by the magneto-rotational instability (MRI) when convection is present, without requiring a net magnetic flux. We independently recover this result by carrying out PLUTO MHD simulations of vertically stratified, radiative, shearing boxes with the thermodynamics and opacities appropriate to DN. The results are robust against the choice of vertical boundary conditions. The thermal equilibrium solutions found by the simulations trace the well-known S-curve in the density-temperature plane. We confirm that the high values of \alpha ~ 0.1 occur near the tip of the hot branch of the S-curve, where convection is active. However, we also present thermally-stable simulations at lower temperatures that have standard values of \alpha ~ 0.03 despite the presence of vigorous convection. We find no simple relationship between \alpha and the strength of the convection, as measured by the ratio of convective to radiative flux. The cold branch is only very weakly ionized so, in the second part of this work, we studied the impact of non-ideal MHD effects on transport. We include resistivity in the simulations and find that the MRI-driven transport is quenched (\alpha ~ 0) below the critical density at which the magnetic Reynolds number R_m \leq 10^4. This is problematic as X-ray emission observed in quiescent systems requires ongoing accretion onto the WD. We verify that these X-rays cannot self-sustain MRI-driven turbulence by photo-ionizing the disk and discuss possible solutions to the issue of accretion in quiescence.

Planets are born in protostellar disks, which are now observed with enough resolution to address questions about internal gas flows. Candidates for driving the flows include magnetic forces, but ionization state estimates suggest much of the gas mass decouples from magnetic fields. Thus, hydrodynamical instabilities could play a major role. We investigate disk dynamics under conditions typical for a T Tauri system, using global 3D radiation hydrodynamics simulations with embedded particles and a resolution of 70 cells per scale height. Stellar irradiation heating is included with realistic dust opacities. The disk starts in joint radiative balance and hydrostatic equilibrium. The vertical shear instability (VSI) develops into turbulence that persists up to at least 1600 inner orbits (143 outer orbits). Turbulent speeds are a few percent of the local sound speed at the midplane, increasing to 20%, or 100 m/s, in the corona. These are consistent with recent upper limits on turbulent speeds from optically thin and thick molecular line observations of TW Hya and HD 163296. The predominantly vertical motions induced by the VSI efficiently lift particles upwards. Grains 0.1 and 1 mm in size achieve scale heights greater than expected in isotropic turbulence. We conclude that while kinematic constraints from molecular line emission do not directly discriminate between magnetic and nonmagnetic disk models, the small dust scale heights measured in HL Tau and HD 163296 favor turbulent magnetic models, which reach lower ratios of the vertical kinetic energy density to the accretion stress.

Recent observations of the protoplanetary disk around the Herbig Be star HD 100546 show two bright features in infrared (H and L' bands) at about 50 AU, with one so far unexplained. We explore the observational signatures of a high mass planet causing shock heating in order to determine if it could be the source of the unexplained infrared feature in HD 100546. More fundamentally, we identify and characterize planetary shocks as an extra, hitherto ignored, source of luminosity in transition disks. The RADMC-3D code is used to perform dust radiative transfer calculations on the hydrodynamical disk models, including volumetric heating. A stronger shock heating rate by a factor 20 would be necessary to qualitatively reproduce the morphology of the second infrared source. Instead, we find that the outer edge of the gap carved by the planet heats up by about 50% relative to the initial reference temperature, which leads to an increase in the scale height. The bulge is illuminated by the central star, producing a lopsided feature in scattered light, as the outer gap edge shows an asymmetry in density and temperature attributable to a secondary spiral arm launched not from the Lindblad resonances but from the 2:1 resonance. We conclude that high-mass planets lead to shocks in disks that may be directly observed, particularly at wavelengths of 10 mm or longer, but that they are more likely to reveal their presence in scattered light by puffing up their outer gap edges and exciting multiple spiral arms.

We conducted a detailed radiative transfer modeling of the dust emission from the circumstellar disk around HL Tau. The goal of our study is to derive the surface density profile of the inner disk and its structure. In addition to the Atacama Large Millimeter/submillimeter Array images at Band 3 (2.9mm), Band 6 (1.3mm), and Band 7 (0.87mm), the most recent Karl G. Jansky Very Large Array (VLA) observations at 7mm were included in the analysis. A simulated annealing algorithm was invoked to search for the optimum model. The radiative transfer analysis demonstrates that most radial components (i.e., >6AU) of the disk become optically thin at a wavelength of 7mm, which allows us to constrain, for the first time, the dust density distribution in the inner region of the disk. We found that a homogeneous grain size distribution is not sufficient to explain the observed images at different wavelengths simultaneously, while models with a shallower grain size distribution in the inner disk work well. We found clear evidence that larger grains are trapped in the first bright ring. Our results imply that dust evolution has already taken place in the disk at a relatively young (i.e., ~1Myr) age. We compared the midplane temperature distribution, optical depth, and properties of various dust rings with those reported previously. Using the Toomre parameter, we briefly discussed the gravitational instability as a potential mechanism for the origin of the dust clump detected in the first bright ring via the VLA observations.

Context. Recent high angular resolution polarimetric continuum observations of circumstellar disks provide new insights into their magnetic field. However, direct constraints are limited to the plane of sky component of the magnetic field. Observations of Zeeman split spectral lines are a potential approach to enhance these insights by providing complementary information. Aims. We investigate which constraints for magnetic fields in circumstellar disks can be obtained from Zeeman observations of the $113~\mathrm{GHz}$ CN lines. Furthermore, we analyze the requirements to perform these observations and their dependence on selected quantities. Methods. We simulate the Zeeman splitting with the radiative transfer (RT) code POLARIS (Reissl et al. 2016) extended by our Zeeman splitting RT extension ZRAD (Brauer et al. 2017), which is based on the line RT code Mol3D (Ober et al. 2015). Results. We find that Zeeman observations of the $113~\mathrm{GHz}$ CN lines provide significant insights into the magnetic field of circumstellar disks. However, with the capabilities of recent and upcoming instrument/observatories, even spatially unresolved observations would be challenging. Nevertheless, these observations are feasible for the most massive disks with a strong magnetic field and high abundance of CN/H. The most restrictive quantity is the magnetic field strength, which should be at least in the order of $\sim1~\mathrm{mG}$. In addition, the inclination of the disk should be around $60\deg$ to preserve the ability to derive the line-of-sight (LOS) magnetic field strength and to obtain a sufficiently high circularly polarized flux.

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.

We analytically derive the expressions for the structure of the inner region of protoplanetary disks based on the results from the recent hydrodynamical simulations. The inner part of a disk can be divided into four regions: dust-free region with gas temperature in the optically thin limit, optically thin dust halo, optically thick condensation front and the classical optically thick region in order from the inside. We derive the dust-to-gas mass ratio profile in the dust halo using the fact that partial dust condensation regulates the temperature to the dust evaporation temperature. Beyond the dust halo, there is an optically thick condensation front where all the available silicate gas condenses out. The curvature of the condensation surface is determined by the condition that the surface temperature must be nearly equal to the characteristic temperature $\sim 1200{\,\rm K}$. We derive the mid-plane temperature in the outer two regions using the two-layer approximation with the additional heating by the condensation front for the outermost region. As a result, the overall temperature profile is step-like with steep gradients at the borders between the outer three regions. The borders might act as planet traps where the inward migration of planets due to gravitational interaction with the gas disk stops. The temperature at the border between the two outermost regions coincides with the temperature needed to activate magnetorotational instability, suggesting that the inner edge of the dead zone must lie at this border. The radius of the dead-zone inner edge predicted from our solution is $\sim$ 2-3 times larger than that expected from the classical optically thick temperature.

Many planets orbit within an AU of their stars, raising questions about their origins. Particularly puzzling are the planets found near the silicate sublimation front. We investigate conditions near the front in the protostellar disk around a young intermediate-mass star, using the first global 3-D radiation non-ideal MHD simulations in this context. We treat the starlight heating; the silicate grains sublimation and deposition at the local, time-varying temperature and density; temperature-dependent Ohmic dissipation; and various initial magnetic fields. The results show magnetorotational turbulence around the sublimation front at 0.5 AU. The disk interior to 0.8 AU is turbulent, with velocities exceeding 10% of the sound speed. Beyond 0.8 AU is the dead zone, cooler than 1000 K and with turbulence orders of magnitude weaker. A local pressure maximum just inside the dead zone concentrates solid particles, favoring their growth. Over many orbits, a vortex develops at the dead zones inner edge, increasing the disks thickness locally by around 10%. We synthetically observe the results using Monte Carlo transfer calculations, finding the sublimation front is near-infrared bright. The models with net vertical magnetic fields develop extended, magnetically-supported atmospheres that reprocess extra starlight, raising the near-infrared flux 20%.The vortex throws a nonaxisymmetric shadow on the outer disk. At wavelengths > 2 micron, the flux varies several percent on monthly timescales. The variations are more regular when the vortex is present. The vortex is directly visible as an arc at ultraviolet through near-infrared wavelengths, given sub-AU spatial resolution.

[Abridged] Regions of low ionisation where the activity of the magneto-rotational instability is suppressed, the so-called dead zones, have been suggested to explain gaps and asymmetries of transition disks. We investigate the gas and dust evolution simultaneously assuming simplified prescriptions for a dead zone and a magnetohydrodynamic (MHD) wind acting on the disk. We explore whether the resulting gas and dust distribution can create signatures similar to those observed in transition disks. For the dust evolution, we included the transport, growth, and fragmentation of dust particles. To compare with observations, we produced synthetic images in scattered optical light and in thermal emission at mm wavelengths. In all models with a dead zone, a bump in the gas surface density is produced that is able to efficiently trap large particles ($\gtrsim 1$ mm) at the outer edge of the dead zone. The gas bump reaches an amplitude of a factor of $\sim5$, which can be enhanced by the presence of an MHD wind that removes mass from the inner disk. While our 1D simulations suggest that such a structure can be present only for $\sim$1 Myr, the structure may be maintained for a longer time when more realistic 2D/3D simulations are performed. In the synthetic images, gap-like low-emission regions are seen at scattered light and in thermal emission at mm wavelengths, as previously predicted in the case of planet-disk interaction. As a conclusion, main signatures of transition disks can be reproduced by assuming a dead zone in the disk, such as gap-like structure in scattered light and millimetre continuum emission, and a lower gas surface density within the dead zone. Previous studies showed that the Rossby wave instability can also develop at the edge of such dead zones, forming vortices and also creating asymmetries.

Magnetic fields significantly influence the evolution of protoplanetary disks and the formation of planets, following the predictions of numerous magnetohydrodynamic (MHD) simulations. However, these predictions are yet observationally unconstrained. To validate the predictions on the influence of magnetic fields on protoplanetary disks, we apply 3D radiative transfer simulations of the polarized emission of aligned aspherical dust grains that directly link 3D global non-ideal MHD simulations to ALMA observations. Our simulations show that it is feasible to observe the predicted toroidal large-scale magnetic field structures, not only in the ideal observations but also with high-angular resolution ALMA observations. Our results show further that high angular resolution observations by ALMA are able to identify vortices embedded in outer magnetized disk regions.

Many stars host planets orbiting within a few astronomical units (AU). The occurrence rate and distributions of masses and orbits vary greatly with the host stars mass. These close planets origins are a mystery that motivates investigating protoplanetary disks central regions. A key factor governing the conditions near the star is the silicate sublimation front, which largely determines where the starlight is absorbed, and which is often called the inner rim. We present the first radiation hydrodynamical modeling of the sublimation front in the disks around the young intermediate-mass stars called Herbig Ae stars. The models are axisymmetric, and include starlight heating, silicate grains sublimating and condensing to equilibrium at the local, time-dependent temperature and density, and accretion stresses parametrizing the results of MHD magneto-rotational turbulence models. The results compare well with radiation hydrostatic solutions, and prove to be dynamically stable. Passing the model disks into Monte Carlo radiative transfer calculations, we show that the models satisfy observational constraints on the inner rims location. A small optically-thin halo of hot dust naturally arises between the inner rim and the star. The inner rim has a substantial radial extent, corresponding to several disk scale heights. While the fronts overall position varies with the stellar luminosity, its radial extent depends on the mass accretion rate. A pressure maximum develops near the location of thermal ionization at temperatures about 1000 K. The pressure maximum is capable of halting solid pebbles radial drift and concentrating them in a zone where temperatures are sufficiently high for annealing to form crystalline silicates.

The first long-baseline ALMA campaign resolved the disk around the young star HL Tau into a number of axisymmetric bright and dark rings. Despite the very young age of HL Tau these structures have been interpreted as signatures for the presence of (proto)planets. The ALMA images triggered numerous theoretical studies based on disk-planet interactions, magnetically driven disk structures, and grain evolution. Of special interest are the inner parts of disks, where terrestrial planets are expected to form. However, the emission from these regions in HL Tau turned out to be optically thick at all ALMA wavelengths, preventing the derivation of surface density profiles and grain size distributions. Here, we present the most sensitive images of HL Tau obtained to date with the Karl G. Jansky Very Large Array at 7.0 mm wavelength with a spatial resolution comparable to the ALMA images. At this long wavelength the dust emission from HL Tau is optically thin, allowing a comprehensive study of the inner disk. We obtain a total disk dust mass of 0.001 - 0.003 Msun, depending on the assumed opacity and disk temperature. Our optically thin data also indicate fast grain growth, fragmentation, and formation of dense clumps in the inner densest parts of the disk. Our results suggest that the HL Tau disk may be actually in a very early stage of planetary formation, with planets not already formed in the gaps but in the process of future formation in the bright rings.

Dust grains with sizes around (sub)mm are expected to couple only weakly to the gas motion in regions beyond 10 au of circumstellar disks. In this work, we investigate the influence of the spatial distribution of such grains on the (sub)mm appearance of magnetized protoplanetary disks. We perform non-ideal global 3D magneto-hydrodynamic stratified disk simulations including particles of different sizes (50 micron to 1 cm), using a Lagrangian particle solver. We calculate the spatial dust temperature distribution, including the dynamically coupled submicron-sized dust grains, and derive ideal continuum re-emission maps of the disk through radiative transfer simulations. Finally, we investigate the feasibility to observe specific structures in the thermal re-emission maps with the Atacama Large Millimeter/submillimeter Array (ALMA). The pressure bump close to the outer edge of the dead-zone leads to particle trapping in ring structures. More specifically, vortices in the disk concentrate the dust and create an inhomogeneous distribution of the solid material in the azimuthal direction. The large-scale disk perturbations are preserved in the (sub)mm re-emission maps. The observable structures are very similar to those expected to result from planet-disk interaction. The larger dust particles increase the brightness contrast between gap and ring structures. We find that rings, gaps and the dust accumulation in the vortex could be traced with ALMA down to a scale of a few astronomical units in circumstellar disks located in nearby star-forming regions. Finally, we present a brief comparison of these structures with those recently found with ALMA in the young circumstellar disks of HL Tau and Oph IRS 48.

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

Recent observations by the Atacama Large Millimeter/submillimeter Array (ALMA) of disks around young stars revealed distinct asymmetries in the dust continuum emission. In this work we want to study axisymmetric and non-axisymmetric structures, evocated by the magneto-rotational instability in the outer regions of protoplanetary disks. We combine the results of state-of-the-art numerical simulations with post-processing radiative transfer (RT) to generate synthetic maps and predictions for ALMA. We performed non-ideal global 3D MHD stratified simulations of the dead-zone outer edge using the FARGO MHD code PLUTO. The stellar and disk parameters are taken from a parameterized disk model applied for fitting high-angular resolution multi-wavelength observations of circumstellar disks. The 2D temperature and density profiles are calculated consistently from a given surface density profile and Monte-Carlo radiative transfer. The 2D Ohmic resistivity profile is calculated using a dust chemistry model. The magnetic field is a vertical net flux field. The resulting dust reemission provides the basis for the simulation of observations with ALMA. The fiducial model develops a large gap followed by a jump in surface density located at the dead-zone outer edge. The jump in density and pressure is strong enough to stop the radial drift of particles. In addition, we observe the generation of vortices by the Rossby wave instability (RWI) at the jumps location close to 60 AU. The vortices are steadily generated and destroyed at a cycle of 40 local orbits. The RT results and simulated ALMA observations predict the feasibility to observe such large scale structures appearing in magnetized disks without having a planet.

Our aim is to study the thermal and dynamical evolution of protoplanetary disks in global simulations, including the physics of radiation transfer and magneto-hydrodynamic (MHD) turbulence caused by the magneto-rotational instability. We develop a radiative transfer method based on the flux-limited diffusion approximation that includes frequency dependent irradiation by the central star. This hybrid scheme is implemented in the PLUTO code. The focus of our implementation is on the performance of the radiative transfer method. Using an optimized Jacobi preconditioned BiCGSTAB solver, the radiative module is three times faster than the MHD step for the disk setup we consider. We obtain weak scaling efficiencies of 70% up to 1024 cores. We present the first global 3D radiation MHD simulations of a stratified protoplanetary disk. The disk model parameters are chosen to approximate those of the system AS 209 in the star-forming region Ophiuchus. Starting the simulation from a disk in radiative and hydrostatic equilibrium, the magnetorotational instability quickly causes MHD turbulence and heating in the disk. For the disk parameters we use, turbulent dissipation heats the disk midplane and raises the temperature by about 15% compared to passive disk models. A roughly flat vertical temperature profile establishes in the disk optically thick region close to the midplane. We reproduce the vertical temperature profile with a viscous disk models for which the stress tensor vertical profile is flat in the bulk of the disk and vanishes in the disk corona. The present paper demonstrates for the first time that global radiation MHD simulations of turbulent protoplanetary disks are feasible with current computational facilities. This opens up the windows to a wide range of studies of the dynamics of protoplanetary disks inner parts, for which there are significant observational constraints.

We investigate the characteristic properties of self-sustained MRI turbulence in low-ionized proto-planetary disks. We study the transition regime between active and dead-zone, performing 3D global non-ideal MHD simulations of stratified disk covering range of magnetic Reynolds number between 2700 < Rm < 6600. We found converged and saturated MRI turbulence for Rm > 5000 with a strength of alpha ~ 0.01. Below Rm < 5000 the MRI starts to decay at the midplane, having Elsasser numbers below one. We find a transition regime between 3300 < Rm < 5000 where the MRI turbulence is still sustained but damped. At around Rm < 3000 the MRI turbulence decays but could reestablished due to the accumulation of toroidal magnetic field or the radial transport of magnetic field from the active region. Below Rm < 3000 the MRI cannot be sustained and is decaying. Here hydro-dynamical motions, like density waves dominate. We observe anti-cyclonic vortices in the transition between dead-zone and active zone.

Explicit numerical computations of super-fast differentially rotating disks are subject to the time-step constraint imposed by the Courant condition. When the bulk orbital velocity largely exceeds any other wave speed the time step is considerably reduced and a large number of steps may be necessary to complete the computation. We present a robust numerical scheme to overcome the Courant limitation by extending the algorithm previously known as FARGO (Fast Advection in Rotating Gaseous Objects) to the equations of magnetohydrodynamics (MHD). The proposed scheme conserves total angular momentum and energy to machine precision and works in Cartesian, cylindrical, or spherical coordinates. The algorithm is implemented in the PLUTO code for astrophysical gasdynamics and is suitable for local or global simulations of accretion or proto-planetary disk models. By decomposing the total velocity into an average azimuthal contribution and a residual term, the algorithm solves the MHD equations through a linear transport step in the orbital direction and a standard nonlinear solver applied to the MHD equations written in terms of the residual velocity. Since the former step is not subject to any stability restriction, the Courant condition is computed only in terms of the residual velocity, leading to substantially larger time steps. The magnetic field is advanced in time using the constrained transport method in order to preserve the divergence-free condition. Conservation of total energy and angular momentum is enforced at the discrete level by properly expressing the source terms in terms of upwind fluxes available during the standard solver. Our results show that applications of the proposed orbital-advection scheme to problems of astrophysical relevance provides, at reduced numerical cost, equally accurate and less dissipative results than standard time-marching schemes.

We investigate the significance of large scale azimuthal, magnetic and velocity modes for the MRI turbulence in accretion disks. We perform 3D global ideal MHD simulations of global stratified proto-planetary disk models. Our domains span azimuthal angles of \pi/4, \pi/2, \pi and 2\pi. We observe up to 100% stronger magnetic fields and stronger turbulence for the restricted azimuthal domain models \pi/2 and \pi/4 compared to the full 2\pi model. We show that for those models, the Maxwell Stress is larger due to strong axisymmetric magnetic fields, generated by the \alpha \Omega dynamo. Large radial extended axisymmetric toroidal fields trigger temporal magnification of accretion stress. All models display a positive dynamo-\alpha in the northern hemisphere (upper disk). The parity is distinct in each model and changes on timescales of 40 local orbits. In model 2\pi, the toroidal field is mostly antisymmetric in respect to the midplane. The eddies of the MRI turbulence are highly anisotropic. The major wavelengths of the turbulent velocity and magnetic fields are between one and two disk scale heights. At the midplane, we find magnetic tilt angles around 8-9 degree increasing up to 12-13 degree in the corona. We conclude that an azimuthal extent of \pi is sufficient to reproduce most turbulent properties in 3D global stratified simulations of magnetised accretion disks.

We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine astrophysical MHD methods actively used to model star formation. The set of nine codes includes: ENZO, FLASH, KT-MHD, LL-MHD, PLUTO, PPML, RAMSES, STAGGER, and ZEUS. We present a comprehensive set of statistical measures designed to quantify the effects of numerical dissipation in these MHD solvers. We compare power spectra for basic fields to determine the effective spectral bandwidth of the methods and rank them based on their relative effective Reynolds numbers. We also compare numerical dissipation for solenoidal and dilatational velocity components to check for possible impacts of the numerics on small-scale density statistics. Finally, we discuss convergence of various characteristics for the turbulence decay test and impacts of various components of numerical schemes on the accuracy of solutions. We show that the best performing codes employ a consistently high order of accuracy for spatial reconstruction of the evolved fields, transverse gradient interpolation, conservation law update step, and Lorentz force computation. The best results are achieved with divergence-free evolution of the magnetic field using the constrained transport method, and using little to no explicit artificial viscosity. Codes which fall short in one or more of these areas are still useful, but they must compensate higher numerical dissipation with higher numerical resolution. This paper is the largest, most comprehensive MHD code comparison on an application-like test problem to date. We hope this work will help developers improve their numerical algorithms while helping users to make informed choices in picking optimal applications for their specific astrophysical problems.

We performed 3D MHD simulations of planet migration in stratified disks using the Godunov code PLUTO, where the disk is turbulent due to the magnetorotational instability. We study the migration for planets with different planet-star mass ratios $q=M_{p}/M_{s}$. In agreement with previous studies, for the low-mass planet cases ($q=5\times10^{-6}$ and $10^{-5}$), migration is dominated by random fluctuations in the torque. For a Jupiter-mass planet $(q=M_{p}/M_{s}=10^{-3}$ for $M_{s}=1M_{\odot})$, we find a reduction of the magnetic stress inside the orbit of the planet and around the gap region. After an initial stage where the torque on the planet is positive, it reverses and we recover migration rates similar to those found in disks where the turbulent viscosity is modelled by an $\alpha$ viscosity. For the intermediate-mass planets ($q=5\times10^{-5}, 10^{-4}$ and $2\times10^{-4}$) we find a new and so far unexpected behavior. In some cases they experience sustained and systematic outwards migration for the entire duration of the simulation. For this case, the horseshoe region is resolved and torques coming from the corotation region can remain unsaturated due to the stresses in the disk. These stresses are generated directly by the magnetic field. The magnitude of the horseshoe drag can overcome the negative Lindblad contribution when the local surface density profile is flat or increasing outwards, which we see in certain locations in our simulations due to the presence of a zonal flow. The intermediate-mass planet is migrating radially outwards in locations where there is a positive gradient of a pressure bump (zonal flow).

We present full 2 Pi global 3-D stratified MHD simulations of accretion disks. We interpret our results in the context of proto-planetary disks. We investigate the turbulence driven by the magneto-rotational instability (MRI) using the PLUTO Godunov code in spherical coordinates with the accurate and robust HLLD Riemann solver. We follow the turbulence for more than 1500 orbits at the innermost radius of the domain to measure the overall strength of turbulent motions and the detailed accretion flow pattern. We find that regions within two scale heights of the midplane have a turbulent Mach number of about 0.1 and a magnetic pressure two to three orders of magnitude less than the gas pressure, while outside three scale heights the magnetic pressure equals or exceeds the gas pressure and the turbulence is transonic, leading to large density fluctuations. The strongest large-scale density disturbances are spiral density waves, and the strongest of these waves has m=5. No clear meridional circulation appears in the calculations because fluctuating radial pressure gradients lead to changes in the orbital frequency, comparable in importance to the stress gradients that drive the meridional flows in viscous models. The net mass flow rate is well-reproduced by a viscous model using the mean stress distribution taken from the MHD calculation. The strength of the mean turbulent magnetic field is inversely proportional to the radius, so the fields are approximately force-free on the largest scales. Consequently the accretion stress falls off as the inverse square of the radius.

We briefly overview our newly developed radiation transport module for MHD simulations and two actual applications. The method combines the advantage of the speed of the Flux-Limited Diffusion approximation and the high accuracy obtained in ray-tracing methods.

The poorly-ionized interior of the protoplanetary disk is the location where dust coagulation processes may be most efficient. However even here, planetesimal formation may be limited by the loss of solid material through radial drift, and by collisional fragmentation of the particles. Our aim is to investigate the possibility that solid particles are trapped at local pressure maxima in the dynamically evolving disk. We perform the first 3-D global non-ideal MHD calculations of the disk treating the turbulence driven by the magneto-rotational instability. The domain contains an inner MRI-active region near the young star and an outer midplane dead zone, with the transition between the two modeled by a sharp increase in the magnetic diffusivity. The azimuthal magnetic fields generated in the active zone oscillate over time, changing sign about every 150 years. We thus observe the radial structure of the `butterfly pattern' seen previously in local shearing-box simulations. The mean magnetic field diffuses from the active zone into the dead zone, where the Reynolds stress nevertheless dominates. The greater total accretion stress in the active zone leads to a net reduction in the surface density, so that after 800 years an approximate steady state is reached in which a local radial maximum in the midplane pressure lies near the transition radius. We also observe the formation of density ridges within the active zone. The dead zone in our models possesses a mean magnetic field, significant Reynolds stresses and a steady local pressure maximum at the inner edge, where the outward migration of planetary embryos and the efficient trapping of solid material are possible.

We employ the PLUTO code for computational astrophysics to assess and compare the validity of different numerical algorithms on simulations of the magneto-rotational instability in 3D accretion disks. In particular we stress on the importance of using a consistent upwind reconstruction of the electro-motive force (EMF) when using the constrained transport (CT) method to avoid the onset of numerical instabilities. We show that the electro-motive force (EMF) reconstruction in the classical constrained transport (CT) method for Godunov schemes drives a numerical instability. The well-studied linear growth of magneto-rotational instability (MRI) is used as a benchmark for an inter-code comparison of PLUTO and ZeusMP. We reproduce the analytical results for linear MRI growth in 3D global MHD simulations and present a robust and accurate Godunov code which can be used for 3D accretion disk simulations in curvilinear coordinate systems.