Outflows driven by large-scale magnetic fields likely play an important role in the evolution and dispersal of protoplanetary disks, and in setting the conditions for planet formation. We extend our 2-D axisymmetric non-ideal MHD model of these outflows by incorporating radiative transfer and simplified thermochemistry, with the twin aims of exploring how heating influences wind launching, and illustrating how such models can be tested through observations of diagnostic spectral lines. Our model disks launch magnetocentrifugal outflows primarily through magnetic tension forces, so the mass-loss rate increases only moderately when thermochemical effects are switched on. For typical field strengths, thermochemical and irradiation heating are more important than magnetic dissipation. We furthermore find that the entrained vertical magnetic flux diffuses out of the disk on secular timescales as a result of non-ideal MHD. Through post-processing line radiative transfer, we demonstrate that spectral line intensities and moment-1 maps of atomic oxygen, the HCN molecule, and other species show potentially observable differences between a model with a magnetically driven outflow and one with a weaker, photoevaporative outflow. In particular, the line shapes and velocity asymmetries in the moment-1 maps could enable the identification of outflows emanating from the disk surface.

Context. The relationship between outflow launching and formation of accretion disks around young stellar objects is still not entirely understood, which is why spectrally and spatially resolved observations are needed. Recently, the Atacama Large Millimetre/sub-millimetre Array (ALMA) has carried out long-baseline observations towards a handful of sources, revealing connections between outflows and the inner regions of disks. Aims. Here we aim to determine the small-scale kinematic and morphological properties of the outflow from the isolated protostar B335 for which no Keplerian disk has, so far, been observed on scales down to 10 au. Methods. We use ALMA in its longest-baseline configuration to observe emission from CO isotopologs, SiO, SO$_2$ and CH$_3$OH. The proximity of B335 provides a resolution of ~3 au (0.03''). We also combine our long-baseline data with archival data to produce a high-fidelity image covering scales up to 700 au (7''). Results. $^{12}$CO has a X-shaped morphology with arms ~50 au in width that we associate with the walls of an outflow cavity, similar to what is observed on larger scales. Long-baseline continuum emission is confined to <7 au of the protostar, while short-baseline continuum emission follows the $^{12}$CO outflow and cavity walls. Methanol is detected within ~30 au of the protostar. SiO is also detected in the vicinity of the protostar, but extended along the outflow. Conclusions. The $^{12}$CO outflow shows no clear signs of rotation at distances $\gtrsim$30 au from the protostar. SiO traces the protostellar jet on small scales, but without obvious rotation. CH$_3$OH and SO$_2$ trace a region <16 au in diameter, centred on the continuum peak, which is clearly rotating. Using episodic, high-velocity, $^{12}$CO features, we estimate the launching radius of the outflow to be <0.1 au and dynamical timescales on the order of a few years.

We present 2.5-D global, ideal MHD simulations of magnetically and rotationally driven protostellar jets from Keplerian accretion discs, wherein only the initial magnetic field strength at the inner radius of the disc, $B_{\rm i}$, is varied. Using the AMR-MHD code AZEUS, we self-consistently follow the jet evolution into the observational regime ($>10^3\,\mathrm{AU}$) with a spatial dynamic range of $\sim6.5\times10^5$. The simulations reveal a three-component outflow: 1) A hot, dense, super-fast and highly magnetised 'jet core'; 2) a cold, rarefied, trans-fast and highly magnetised 'sheath' surrounding the jet core and extending to a tangential discontinuity; and 3) a warm, dense, trans-slow and weakly magnetised shocked ambient medium entrained by the advancing bow shock. The simulations reveal power-law relationships between $B_{\rm i}$ and the jet advance speed, $v_{\rm jet}$, the average jet rotation speed, $\langle v_\varphi\rangle$, as well as fluxes of mass, momentum, and kinetic energy. Quantities that do not depend on $B_{\rm i}$ include the plasma-$\beta$ of the transported material which, in all cases, seems to asymptote to order unity. Jets are launched by a combination of the 'magnetic tower' and 'bead-on-a-wire' mechanisms, with the former accounting for most of the jet acceleration---even for strong fields---and continuing well beyond the fast magnetosonic point. At no time does the leading bow shock leave the domain and, as such, these simulations generate large-scale jets that reproduce many of the observed properties of protostellar jets including their characteristic speeds and transported fluxes.

From molecular clouds to protoplanetary disks, non-ideal magnetic effects are important in many astrophysical environments. Indeed, in star and disk formation processes, it has become clear that these effects are critical to the evolution of the system. The efficacy of non-ideal effects are, however, determined by the complex interplay between magnetic fields, ionising radiation, cosmic rays, microphysics, and chemistry. In order to understand these key microphysical parameters, we present a one-dimensional non-ideal magnetohydrodynamics code and apply it to a model of a time-dependent, oblique, magnetic shock wave. By varying the microphysical ingredients of the model, we find that cosmic rays and dust play a major role, and that, despite the uncertainties, the inclusion of microphysics is essential to obtain a realistic outcome in magnetic astrophysical simulations.

We investigate the effects of radiative energy transfer on a series of nested-grid, high-resolution hydrodynamic simulations of gas and particle dynamics in the vicinity of an Earth-mass planetary embryo. We include heating due to the accretion of solids and the subsequent convective motions. Using a constant embryo surface temperature, we show that radiative energy transport results in a tendency to reduce the entropy in the primordial atmosphere, but this tendency is alleviated by an increase in the strength of convective energy transport, triggered by a correspondingly increased super-adiabatic temperature gradient. As a consequence, the amplitude of the convective motions increase by roughly an order of magnitude in the vicinity of the embryo. In the cases investigated here, where the optical depth towards the disk surface is larger than unity, the reduction of the temperature in the outer parts of the Hill sphere relative to cases without radiative energy transport is only $\sim$100K, while the mass density increase is on the order of a factor of two in the inner parts of the Hill sphere. Our results demonstrate that, unless unrealistically low dust opacities are assumed, radiative cooling in the context of primordial rocky planet atmospheres can only become important after the disk surface density has dropped significantly below minimum-mass-solar-nebula values.

We present nested-grid, high-resolution hydrodynamic simulations of gas and particle dynamics in the vicinity of Mars- to Earth-mass planetary embryos. The simulations extend from the surface of the embryos to a few vertical disk scale heights, with \rev{a spatial} dynamic range \rev{of} $\sim\! 1.4\times 10^5$. Our results confirm that "pebble"-sized particles are readily accreted, with accretion rates continuing to increase up to metre-size "boulders" for a 10\% MMSN surface density model. The gas mass flux in and out of the Hill sphere is consistent with the Hill rate, $\Sigma\Omega R_\mathrm{H}^2 = 4\, 10^{-3}$ M$_\oplus$ yr$^{-1}$. While smaller size particles mainly track the gas, a net accretion rate of $\approx 2\,10^{-5}$ M$_\oplus$ yr$^{-1}$ is reached for 0.3--1 cm particles, even though a significant fraction leaves the Hill sphere again. Effectively all pebble-sized particles that cross the Bondi sphere are accreted. The resolution of these simulations is sufficient to resolve accretion-driven convection. Convection driven by a nominal accretion rate of $10^{-6}$ M$_\oplus$ yr$^{-1}$ does not significantly alter the pebble accretion rate. We find that, due to cancellation effects, accretion rates of pebble-sized particles are nearly independent of disk surface density. As a result, we can estimate accurate growth times for specified particle sizes. For 0.3--1 cm size particles, the growth time from a small seed is $\sim$0.15 million years for an Earth mass planet at 1 AU and $\sim$0.1 million years for a Mars mass planet at 1.5 AU.

We describe a simple and effective algorithm for solving Poisson's equation in the context of self-gravity within the DISPATCH astrophysical fluid framework. The algorithm leverages the fact that DISPATCH stores multiple time slices and uses asynchronous time-stepping to produce a scheme that does not require any explicit global communication or sub-cycling, only the normal, local communication between patches and the iterative solution to Poisson's equation. We demonstrate that the implementation is suitable for both collections of patches of a single resolution and for hierarchies of adaptively resolved patches. Benchmarks are presented that demonstrate the accuracy, effectiveness and efficiency of the scheme.

The growth of dust grains in protoplanetary disks is a necessary first step towards planet formation. This growth has been inferred via observations of thermal dust emission towards mature protoplanetary systems (age >2 million years) with masses that are, on average, similar to Neptune3. In contrast, the majority of confirmed exoplanets are heavier than Neptune. Given that young protoplanetary disks are more massive than their mature counterparts, this suggests that planet formation starts early, but evidence for grain growth that is spatially and temporally coincident with a massive reservoir in young disks remains scarce. Here, we report observations on a lack of emission of carbon monoxide isotopologues within the inner ~15 au of a very young (age ~100,000 years) disk around the Solar-type protostar TMC1A. By using the absence of spatially resolved molecular line emission to infer the gas and dust content of the disk, we conclude that shielding by millimeter-size grains is responsible for the lack of emission. This suggests that grain growth and millimeter-size dust grains can be spatially and temporally coincident with a mass reservoir sufficient for giant planet formation. Hence, planet formation starts during the earliest, embedded phases in the life of young stars.

We introduce a high-performance simulation framework that permits the semi-independent, task-based solution of sets of partial differential equations, typically manifesting as updates to a collection of `patches' in space-time. A hybrid MPI/OpenMP execution model is adopted, where work tasks are controlled by a rank-local `dispatcher' which selects, from a set of tasks generally much larger than the number of physical cores (or hardware threads), tasks that are ready for updating. The definition of a task can vary, for example, with some solving the equations of ideal magnetohydrodynamics (MHD), others non-ideal MHD, radiative transfer, or particle motion, and yet others applying particle-in-cell (PIC) methods. Tasks do not have to be grid-based, while tasks that are, may use either Cartesian or orthogonal curvilinear meshes. Patches may be stationary or moving. Mesh refinement can be static or dynamic. A feature of decisive importance for the overall performance of the framework is that time steps are determined and applied locally; this allows potentially large reductions in the total number of updates required in cases when the signal speed varies greatly across the computational domain, and therefore a corresponding reduction in computing time. Another feature is a load balancing algorithm that operates `locally' and aims to simultaneously minimise load and communication imbalance. The framework generally relies on already existing solvers, whose performance is augmented when run under the framework, due to more efficient cache usage, vectorisation, local time-stepping, plus near-linear and, in principle, unlimited OpenMP and MPI scaling.

Young stars are associated with prominent outflows of molecular gas. The ejection of gas via these outflows is believed to remove angular momentum from the protostellar system, thus permitting young stars to grow by accretion of material from the protostellar disk. The underlying mechanism for outflow ejection is not yet understood, but is believed to be closely linked to the protostellar disk. Assorted scenarios have been proposed to explain protostellar outflows; the main difference between these models is the region where acceleration of material takes place: close to the protostar itself ('X-wind', or stellar wind), in a larger region throughout the protostellar disk (disk wind), or at the interface between. Because of the limits of observational studies, outflow launching regions have so far only been probed by indirect extrapolation. Here we report observations of carbon monoxide toward the outflow associated with the TMC1A protostellar system. These data show that gas is ejected from a region extending up to a radial distance of 25 astronomical units from the central protostar, and that angular momentum is removed from an extended region of the disk. This demonstrates that the outflowing gas is launched by an extended disk wind from a Keplerian disk. Hence, we rule out X-wind and stellar wind launching scenarios as the source of the emission on the scales we observe.

Aims. The importance of radiation to the physical structure of protoplanetary disks cannot be understated. However, protoplanetary disks evolve with time, and so to understand disk evolution and by association, disk structure, one should solve the combined and time-dependent equations of radiation hydrodynamics. Methods. We implement a new implicit radiation solver in the AZEuS adaptive mesh refinement magnetohydrodynamics fluid code. Based on a hybrid approach that combines frequency-dependent ray-tracing for stellar irradiation with non-equilibrium flux limited diffusion, we solve the equations of radiation hydrodynamics while preserving the directionality of the stellar irradiation. The implementation permits simulations in Cartesian, cylindrical, and spherical coordinates, on both uniform and adaptive grids. Results. We present several hydrostatic and hydrodynamic radiation tests which validate our implementation on uniform and adaptive grids as appropriate, including benchmarks specifically designed for protoplanetary disks. Our results demonstrate that the combination of a hybrid radiation algorithm with AZEuS is an effective tool for radiation hydrodynamics studies, and produces results which are competitive with other astrophysical radiation hydrodynamics codes.

Context. The growth process of dust particles in protoplanetary disks can be modeled via numerical dust coagulation codes. In this approach, physical effects that dominate the dust growth process often must be implemented in a parameterized form. Due to a lack of these parameterizations, existing studies of dust coagulation have ignored the effects a hydrodynamical gas flow can have on grain growth, even though it is often argued that the flow could significantly contribute either positively or negatively to the growth process. Aims. We intend to provide a quantification of hydrodynamical effects on the growth of dust particles, such that these effects can be parameterized and implemented in a dust coagulation code. Methods. We numerically integrate the trajectories of small dust particles in the flow of disk gas around a proto-planetesimal, sampling a large parameter space in proto-planetesimal radii, headwind velocities, and dust stopping times. Results. The gas flow deflects most particles away from the proto-planetesimal, such that its effective collisional cross section, and therefore the mass accretion rate, is reduced. The gas flow however also reduces the impact velocity of small dust particles onto a proto-planetesimal. This can be beneficial for its growth, since large impact velocities are known to lead to erosion. We also demonstrate why such a gas flow does not return collisional debris to the surface of a proto-planetesimal. Conclusions. We predict that a laminar hydrodynamical flow around a proto-planetesimal will have a significant effect on its growth. However, we cannot easily predict which result, the reduction of the impact velocity or the sweep-up cross section, will be more important. Therefore, we provide parameterizations ready for implementation into a dust coagulation code.

A new adaptive mesh refinement (AMR) version of the ZEUS-3D astrophysical magnetohydrodynamical (MHD) fluid code, AZEuS, is described. The AMR module in AZEuS has been completely adapted to the staggered mesh that characterises the ZEUS family of codes, on which scalar quantities are zone-centred and vector components are face-centred. In addition, for applications using static grids, it is necessary to use higher-order interpolations for prolongation to minimise the errors caused by waves crossing from a grid of one resolution to another. Finally, solutions to test problems in 1-, 2-, and 3-dimensions in both Cartesian and spherical coordinates are presented.

We present the first 2.5-D MHD simulations of protostellar jets that include both the region in which the jet is launched magnetocentrifugally at scale lengths $<0.1$ AU, and where the propagating jet is observed at scale lengths $>10^3$ AU. These simulations, performed with the new AMR-MHD code AZEuS, reveal interesting relationships between conditions at the disc surface, such as the magnetic field strength, and direct observables such as proper motion, jet rotation, jet radius, and mass flux. By comparing these quantities with observed values, we present direct numerical evidence that the magnetocentrifugal launching mechanism is capable, by itself, of launching realistic protostellar jets.