scholarly journals Ring formation and dust dynamics in wind-driven protoplanetary discs: global simulations

2020 ◽  
Vol 639 ◽  
pp. A95 ◽  
Author(s):  
A. Riols ◽  
G. Lesur ◽  
F. Menard

Large-scale vertical magnetic fields are believed to play a key role in the evolution of protoplanetary discs. Associated with non-ideal effects, such as ambipolar diffusion, they are known to launch a wind that could drive accretion in the outer part of the disc (R >  1 AU). They also potentially lead to self-organisation of the disc into large-scale axisymmetric structures, similar to the rings recently imaged by sub-millimetre or near-infrared instruments (ALMA and SPHERE). The aim of this paper is to investigate the mechanism behind the formation of these gaseous rings, but also to understand the dust dynamics and its emission in discs threaded by a large-scale magnetic field. To this end, we performed global magneto-hydrodynamics (MHD) axisymmetric simulations with ambipolar diffusion using a modified version of the PLUTO code. We explored different magnetisations with the midplane β parameter ranging from 105 to 103 and included dust grains -treated in the fluid approximation- ranging from 100 μm to 1 cm in size. We first show that the gaseous rings (associated with zonal flows) are tightly linked to the existence of MHD winds. Secondly, we find that millimetre-size dust is highly sedimented, with a typical scale height of 1 AU at R = 100 AU for β = 104, compatible with recent ALMA observations. We also show that these grains concentrate into pressure maxima associated with zonal flows, leading to the formation of dusty rings. Using the radiative transfer code MCFOST, we computed the dust emission and make predictions on the ring-gap contrast and the spectral index that one might observe with interferometers like ALMA.

2018 ◽  
Vol 617 ◽  
pp. A117 ◽  
Author(s):  
A. Riols ◽  
G. Lesur

Context. Magnetohydrodynamic (MHD) turbulence plays a crucial role in the dust dynamics of protoplanetary discs. It affects planet formation, vertical settling, and is one possible origin of the large scale axisymmetric structures, such as rings, recently imaged by ALMA and SPHERE. Among the variety of MHD processes in discs, the magnetorotational instability (MRI) has raised particular interest since it provides a source of turbulence and potentially organizes the flow into large scale structures. However, the weak ionization of discs prevents the MRI from being excited beyond 1 AU. Moreover, the low velocity dispersion observed in CO and strong sedimentation of millimetre dust measured in T-Tauri discs are in contradiction with predictions based on ideal MRI turbulence. Aims. In this paper, we study the effects of non-ideal MHD and magnetized winds on the dynamics and sedimentation of dust grains. We consider a weakly ionized plasma subject to ambipolar diffusion characterizing the disc outer regions (≫1 AU). Methods. To compute the dust and gas motions, we performed numerical MHD simulations in the stratified shearing box, using a modified version of the PLUTO code. We explored different grain sizes from micrometre to few centimetres and different disc vertical magnetizations with plasma beta ranging from 103 to 105. Results. Our simulations show that the mm-cm dust is contained vertically in a very thin layer, with typical heightscale ≲0.4 AU at R = 30 AU, compatible with recent ALMA observations. Horizontally, the grains are trapped within the pressure maxima (or zonal flows) induced by ambipolar diffusion, leading to the formation of dust rings. For micrometre grains and strong magnetization, we find that the dust layer has a size comparable to the disc heightscale H. In this regime, dust settling cannot be explained by a simple 1D diffusion theory but results from a large scale 2D circulation induced by both MHD winds and zonal flows. Conclusions. Our results suggest that non-ideal MHD effects and MHD winds associated with zonal flows play a major role in shaping the radial and vertical distribution of dust in protoplanetary discs. Leading to effective accretion efficiency α ≃ 10−3–10−1, non-ideal MHD models are also a promising avenue to reconcile the low turbulent activity measured in discs with their relatively high accretion rates.


2019 ◽  
Vol 622 ◽  
pp. A132 ◽  
Author(s):  
A. V. Mosenkov ◽  
M. Baes ◽  
S. Bianchi ◽  
V. Casasola ◽  
L. P. Cassarà ◽  
...  

Most radiative transfer models assume that dust in spiral galaxies is distributed exponentially. In this paper our goal is to verify this assumption by analysing the two-dimensional large-scale distribution of dust in galaxies from the DustPedia sample. For this purpose, we have made use of Herschel imaging in five bands, from 100 to 500 μm, in which the cold dust constituent is primarily traced and makes up the bulk of the dust mass in spiral galaxies. For a subsample of 320 disc galaxies, we successfully performed a simultaneous fitting with a single Sérsic model of the Herschel images in all five bands using the multi-band modelling code GALFITM. We report that the Sérsic index n, which characterises the shape of the Sérsic profile, lies systematically below 1 in all Herschel bands and is almost constant with wavelength. The average value at 250 μm is 0.67 ± 0.37 (187 galaxies are fitted with n250 ≤ 0.75, 87 galaxies have 0.75 < n250 ≤ 1.25, and 46 – with n250 >  1.25). Most observed profiles exhibit a depletion in the inner region (at r <  0.3−0.4 of the optical radius r25) and are more or less exponential in the outer part. We also find breaks in the dust emission profiles at longer distances (0.5−0.6) r25 which are associated with the breaks in the optical and near-infrared. We assumed that the observed deficit of dust emission in the inner galaxy region is related to the depression in the radial profile of the HI surface density in the same region because the atomic gas reaches high enough surface densities there to be transformed into molecular gas. If a galaxy has a triggered star formation in the inner region (for example, because of a strong bar instability, which transfers the gas inwards to the centre, or a pseudobulge formation), no depletion or even an excess of dust emission in the centre is observed.


Author(s):  
Clément Baruteau ◽  
Gaylor Wafflard-Fernandez ◽  
Romane Le Gal ◽  
Florian Debras ◽  
Andrés Carmona ◽  
...  

Abstract Predicting how a young planet shapes the gas and dust emission of its parent disc is key to constraining the presence of unseen planets in protoplanetary disc observations. We investigate the case of a 2 Jupiter mass planet that becomes eccentric after migrating into a low-density gas cavity in its parent disc. Two-dimensional hydrodynamical simulations are performed and post-processed by three-dimensional radiative transfer calculations. In our disc model, the planet eccentricity reaches ∼0.25, which induces strong asymmetries in the gas density inside the cavity. These asymmetries are enhanced by photodissociation and form large-scale asymmetries in 12CO J=3→2 integrated intensity maps. They are shown to be detectable for an angular resolution and a noise level similar to those achieved in ALMA observations. Furthermore, the planet eccentricity renders the gas inside the cavity eccentric, which manifests as a narrowing, stretching and twisting of iso-velocity contours in velocity maps of 12CO J=3→2. The planet eccentricity does not, however, give rise to detectable signatures in 13CO and C18O J=3→2 inside the cavity because of low column densities. Outside the cavity, the gas maintains near-circular orbits, and the vertically extended optically thick CO emission displays a four-lobed pattern in integrated intensity maps for disc inclinations $\gtrsim$ 30○. The lack of large and small dust inside the cavity in our model further implies that synthetic images of the continuum emission in the sub-millimetre, and of polarized scattered light in the near-infrared, do not show significant differences when the planet is eccentric or still circular inside the cavity.


2012 ◽  
Vol 10 (H16) ◽  
pp. 615-615
Author(s):  
Michael D. Pavel ◽  
Robert C. Marchwinski ◽  
Dan P. Clemens

Marchwinski et al. (2012) mapped the magnetic field strength across the quiescent cloud GRSMC 45.60+0.30 (shown in Figure 1 subtending 40x10 pc at a distance of 1.88 kpc) with the Chandrasekhar-Fermi method CF; Chandrasekhar & Fermi 1953) using near-infrared starlight polarimetry from the Galactic Plane Infrared Polarization Survey (Clemens et al.2012a, b) and gas properties from the Galactic Ring Survey (Jackson et al.2006). The large-scale magnetic field is oriented parallel to the gas-traced ‘spine’ of the cloud. Seven ‘magnetic cores’ with high magnetic field strength were identified and are coincident with peaks in the gas column density. Calculation of the mass-to-flux ratio (Crutcher 1999) shows that these cores are exclusively magnetically subcritical and that magnetostatic pressure can support them against gravitational collapse.


1990 ◽  
Vol 139 ◽  
pp. 35-47
Author(s):  
G. G. Fazio ◽  
T. M. Dame ◽  
S. Kent

The near-infrared region of the spectrum is a particularly advantageous window for observing the distribution of old, evolved stars in the galactic disk and bulge. These stars are important because they provide an excellent tracer of the overall stellar mass distribution. At shorter wavelengths extinction is a serious problem, and at longer wavelengths the flux is dominated by dust emission. A summary of the large-scale diffuse near-infrared observations of the Galaxy is presented, as is a summary of the results obtained from these data on the structure of the galactic disk and bulge. The importance of combining CO and near-infrared maps of similar resolution to determine a three-dimensional model of galactic extinction is demonstrated. The Spacelab-2 Infrared Telescope (IRT) data are used in conjunction with a proposed galactic model to make preliminary measurements of the global scale parameters of the Galaxy.


2010 ◽  
Vol 6 (S276) ◽  
pp. 89-94
Author(s):  
Anders Johansen ◽  
Hubert Klahr ◽  
Thomas Henning

AbstractWe present high resolution computer simulations of dust dynamics and planetesimal formation in turbulence triggered by the magnetorotational instability. Particles representing approximately meter-sized boulders clump in large scale overpressure regions in the simulation box. These overdensities readily contract due to the combined gravity of the particles to form gravitationally bound clusters with masses ranging from a few to several ten times the mass of the dwarf planet Ceres. Gravitationally bound clumps are observed to collide and merge at both moderate and high resolution. The collisional products form the top end of a distribution of planetesimal masses ranging from less than one Ceres mass to 35 Ceres masses. It remains uncertain whether collisions are driven by dynamical friction or underresolution of clumps.


2020 ◽  
Vol 499 (1) ◽  
pp. 993-1004
Author(s):  
Cameron P M Bell ◽  
Maria-Rosa L Cioni ◽  
A H Wright ◽  
Stefano Rubele ◽  
David L Nidever ◽  
...  

ABSTRACT We present a map of the total intrinsic reddening across ≃34 deg2 of the Small Magellanic Cloud (SMC) derived using optical (ugriz) and near-infrared (IR; YJKs) spectral energy distributions (SEDs) of background galaxies. The reddening map is created using a subsample of 29 274 galaxies with low levels of intrinsic reddening based on the lephare χ2 minimization SED-fitting routine. We find statistically significant enhanced levels of reddening associated with the main body of the SMC compared with regions in the outskirts [ΔE(B − V) ≃ 0.3 mag]. A comparison with literature reddening maps of the SMC shows that, after correcting for differences in the volume of the SMC sampled, there is good agreement between our results and maps created using young stars. In contrast, we find significant discrepancies between our results and maps created using old stars or based on longer wavelength far-IR dust emission that could stem from biased samples in the former and uncertainties in the far-IR emissivity and the optical properties of the dust grains in the latter. This study represents one of the first large-scale categorizations of extragalactic sources behind the SMC and as such we provide the lephare outputs for our full sample of ∼500 000 sources.


2018 ◽  
Vol 614 ◽  
pp. A101 ◽  
Author(s):  
M. S. Väisälä ◽  
F. A. Gent ◽  
M. Juvela ◽  
M. J. Käpylä

Context.Efforts to compare polarization measurements with synthetic observations from magnetohydrodynamics (MHD) models have previously concentrated on the scale of molecular clouds.Aims.We extend the model comparisons to kiloparsec scales, taking into account hot shocked gas generated by supernovae and a non-uniform dynamo-generated magnetic field at both large and small scales down to 4 pc spatial resolution.Methods.We used radiative transfer calculations to model dust emission and polarization on top of MHD simulations. We computed synthetic maps of column densityNH, polarization fractionp, and polarization angle dispersionS, and studied their dependencies on important properties of MHD simulations. These include the large-scale magnetic field and its orientation, the small-scale magnetic field, and supernova-driven shocks.Results.Similar filament-like structures ofSas seen in thePlanckall-sky maps are visible in our synthetic results, although the smallest scale structures are absent from our maps. Supernova-driven shock fronts andSdo not show significant correlation. Instead,Scan clearly be attributed to the distribution of the small-scale magnetic field. We also find that the large-scale magnetic field influences the polarization properties, such that, for a given strength of magnetic fluctuation, a strong plane of the sky mean field weakens the observedS, while strengtheningp. The anticorrelation ofpandS, and decreasingpas a function ofNHare consistent across all synthetic observations. The magnetic fluctuations follow an exponential distribution, rather than Gaussian characteristic of flows with intermittent repetitive shocks.Conclusions.The observed polarization properties and column densities are sensitive to the line-of-sight distance over which the emission is integrated. Studying synthetic maps as the function of maximum integration length will further help with the interpretation of observations. The effects of the large-scale magnetic field orientation on the polarization properties are difficult to be quantified from observations solely, but MHD models might turn out to be useful for separating the effect of the large-scale mean field.


2021 ◽  
Vol 502 (3) ◽  
pp. 3942-3954
Author(s):  
D Hung ◽  
B C Lemaux ◽  
R R Gal ◽  
A R Tomczak ◽  
L M Lubin ◽  
...  

ABSTRACT We present a new mass function of galaxy clusters and groups using optical/near-infrared (NIR) wavelength spectroscopic and photometric data from the Observations of Redshift Evolution in Large-Scale Environments (ORELSE) survey. At z ∼ 1, cluster mass function studies are rare regardless of wavelength and have never been attempted from an optical/NIR perspective. This work serves as a proof of concept that z ∼ 1 cluster mass functions are achievable without supplemental X-ray or Sunyaev-Zel’dovich data. Measurements of the cluster mass function provide important contraints on cosmological parameters and are complementary to other probes. With ORELSE, a new cluster finding technique based on Voronoi tessellation Monte Carlo (VMC) mapping, and rigorous purity and completeness testing, we have obtained ∼240 galaxy overdensity candidates in the redshift range 0.55 &lt; z &lt; 1.37 at a mass range of 13.6 &lt; log (M/M⊙) &lt; 14.8. This mass range is comparable to existing optical cluster mass function studies for the local universe. Our candidate numbers vary based on the choice of multiple input parameters related to detection and characterization in our cluster finding algorithm, which we incorporated into the mass function analysis through a Monte Carlo scheme. We find cosmological constraints on the matter density, Ωm, and the amplitude of fluctuations, σ8, of $\Omega _{m} = 0.250^{+0.104}_{-0.099}$ and $\sigma _{8} = 1.150^{+0.260}_{-0.163}$. While our Ωm value is close to concordance, our σ8 value is ∼2σ higher because of the inflated observed number densities compared to theoretical mass function models owing to how our survey targeted overdense regions. With Euclid and several other large, unbiased optical surveys on the horizon, VMC mapping will enable optical/NIR cluster cosmology at redshifts much higher than what has been possible before.


2021 ◽  
Vol 503 (1) ◽  
pp. 270-291
Author(s):  
F Navarete ◽  
A Damineli ◽  
J E Steiner ◽  
R D Blum

ABSTRACT W33A is a well-known example of a high-mass young stellar object showing evidence of a circumstellar disc. We revisited the K-band NIFS/Gemini North observations of the W33A protostar using principal components analysis tomography and additional post-processing routines. Our results indicate the presence of a compact rotating disc based on the kinematics of the CO absorption features. The position–velocity diagram shows that the disc exhibits a rotation curve with velocities that rapidly decrease for radii larger than 0.1 arcsec (∼250 au) from the central source, suggesting a structure about four times more compact than previously reported. We derived a dynamical mass of 10.0$^{+4.1}_{-2.2}$ $\rm {M}_\odot$ for the ‘disc + protostar’ system, about ∼33 per cent smaller than previously reported, but still compatible with high-mass protostar status. A relatively compact H2 wind was identified at the base of the large-scale outflow of W33A, with a mean visual extinction of ∼63 mag. By taking advantage of supplementary near-infrared maps, we identified at least two other point-like objects driving extended structures in the vicinity of W33A, suggesting that multiple active protostars are located within the cloud. The closest object (Source B) was also identified in the NIFS field of view as a faint point-like object at a projected distance of ∼7000 au from W33A, powering extended K-band continuum emission detected in the same field. Another source (Source C) is driving a bipolar $\rm {H}_2$ jet aligned perpendicular to the rotation axis of W33A.


Sign in / Sign up

Export Citation Format

Share Document