Knowledge of the midplane temperature of protoplanetary disks is one of the key ingredients in theories of dust growth and planet formation. However, direct measurement of this quantity is complicated, and often depends on the fitting of complex models to the data. In this paper we demonstrate a method to directly measure the midplane gas temperature from an optically thick molecular line, if the disk is moderately inclined. The only model assumption that enters is that the line is very optically thick, also in the midplane region where we wish to measure the temperature. Freeze-out of the molecule onto dust grains could thwart this. However, in regions that are expected to be warm enough to avoid freeze-out, this method should work. We apply the method to the CO 2-1 line channel maps of the disk around HD 163296. We find that the midplane temperature between 100 and 400 au drops only mildly from 25 K down to 18 K. While we see no direct evidence of the midplane being optically thin due to strong CO depletion by freeze-out, we cannot rule it out either. The fact that the inferred temperatures are close to the expected CO freeze-out temperature could be an indication of this. Incidently, for the disk around HD 163296 we also find dynamic evidence for a rather abrupt outer edge of the disk, suggestive of outside-in photoevaporation or truncation by an unseen companion.

The grand question of star and planet formation is the distribution of magnetic flux in the protoplanetary disks. To answer it, a detailed self-consistent chemical evolution is needed to describe the magnetic dissipation in the collapsing core accurately. We use a chemo-dynamical version of RAMSES to follow the evolution of collapsing dense cores for a range in the dust size assumptions. The number density of dust and it's mean size are affecting the efficiency of the charge capturing and thus the chemistry. The chemical abundances for the range of dust sizes are produced by RAMSES and serve as an input to calculations of Ohmic, ambipolar and Hall diffusivity terms. We find that Ohmic resistivity only plays a role at the late stage of the collapse where gas density exceeds few times $10^{13}\rm cm^{-3}$. Ambipolar diffusion is a dominant magnetic diffusivity term in cases where mean dust size is a typical ISM value or larger. We show that the assumption of a fixed 'dominant ion' mass can change the ambipolar diffusion up to one order of magnitude. 'Negative' Hall effect is dominant during the collapse in case of mean dust size of 0.02 $\mu$m and smaller, the effect which we connect to the dominance of negatively charged grains. We find that the Hall effect's sign reversal is depending on relative contribution of negatively charged dust to the hall conductivity. The dust grain mean size appears to be the parameter dividing the collapsing clouds in Hall-dominated and ambipolar-dominated clouds, and thus affecting the size of the new-born disks. We propose to link the dust properties and the occurance and size of disk structures in Class 0 YSO's. The proper accounting for dust grain growth in the radiative magneto-hydrodynamical collapse models appear to be be as important as coupling the dynamics of the collapse with the chemistry.

Combining the co-evolving chemistry, hydrodynamics and radiative transfer is an important step for star formation studies. It allows both a better link to observations and a self-consistent monitoring of the magnetic dissipation in the collapsing core. Our aim is to follow a chemo-dynamical evolution of collapsing dense cores with a reduced gas-grain chemical network. We present the results of radiative hydrodynamic (RHD) simulations of 1 M$_\odot$ isolated dense core collapse. The physical setup includes RHD and dynamical evolution of a chemical network. To perform those simulations, we merged the multi-dimensional adaptive-mesh-refinement code RAMSES and the thermo-chemistry Paris-Durham shock code. We simulate the formation of the first hydro-static core (FHSC) and the co-evolution of 56 species describing mainly H-C-O chemistry. Accurate benchmarking is performed, testing the reduced chemical network against a well-establiched complex network. We show that by using a compact set of reactions, one can match closely the CO abundances with results of a much more complex network. Our main results are: (a) We find that gas-grain chemistry post-processing can lead to one order of magnitude lower CO gas-phase abundances compared to the dynamical chemistry, with strongest effect during the isothermal phase of collapse. (b) The free-fall time has little effect on the chemical abundances for our choice of the parameters. (c) Dynamical chemical evolution is required to describe the CO gas phase abundance as well as the CO ice formation for the mean grain size larger then 1$\mu{}$m. (d) Furthermore, dust mean size and size distribution have a strong effect on chemical abundances and hence on the ionization degree and magnetic dissipation. We conclude that dust grain growth in the collapse simulations can be as important as coupling the collapse with chemistry.

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 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.

The edges of magnetically-dead zones in protostellar disks have been proposed as locations where density bumps may arise, trapping planetesimals and helping form planets. Magneto-rotational turbulence in magnetically-active zones provides both accretion of gas on the star and transport of mass to the dead zone. We investigate the location of the magnetically-active regions in a protostellar disk around a solar-type star, varying the disk temperature, surface density profile, and dust-to-gas ratio. We also consider stellar masses between 0.4 and 2 $M_\odot$, with corresponding adjustments in the disk mass and temperature. The dead zone's size and shape are found using the Elsasser number criterion with conductivities including the contributions from ions, electrons, and charged fractal dust aggregates. The charged species' abundances are found using the approach proposed by S. Okuzumi. The dead zone is in most cases defined by the ambipolar diffusion. In our maps, the dead zone takes a variety of shapes, including a fish-tail pointing away from the star and islands located on and off the midplane. The corresponding accretion rates vary with radius, indicating locations where the surface density will increase over time, and others where it will decrease. We show that density bumps do not readily grow near the dead zone's outer edge, independently of the disk parameters and the dust properties. Instead, the accretion rate peaks at the radius where the gas-phase metals freeze out. This could lead to clearing a valley in the surface density, and to a trap for pebbles located just outside the metal freeze-out line.

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 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 have analyzed two epochs of HST/WFPC2 observations of the young Galactic starburst cluster in NGC 3603 with the aim to study its internal dynamics and stellar population. Relative proper motions measured over 10.15 yrs of more than 800 stars enable us to distinguish cluster members from field stars. The best-fitting isochrone yields Av=4.6-4.7 mag, a distance of 6.6-6.9 kpc, and an age of 1 Myr for NGC 3603 Young Cluster (NYC). We identify pre-main-sequence/main-sequence transition stars located in the short-lived radiative-convective gap, which in the NYC occurs in the mass range 3.5-3.8 Msun. We also identify a sparse population of stars with an age of 4 Myr, which appear to be the lower mass counterparts to previously discovered blue supergiants located in the giant HII region NGC 3603. For the first time, we are able to measure the internal velocity dispersion of a starburst cluster from 234 stars with I < 18.5 mag to {\sigma}_pm1D=141+/-27 {\mu}as/yr (4.5+/-0.8 km/s at a distance of 6.75 kpc). As stars with masses between 1.7 and 9 Msun all exhibit the same velocity dispersion, the cluster stars have not yet reached equipartition of kinetic energy (i.e., the cluster is not in virial equilibrium). The results highlight the power of combining high-precision astrometry and photometry, and emphasize the role of NYC as a benchmark object for testing stellar evolution models and dynamical models for young clusters and as a template for extragalactic starburst clusters.

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.