No Arabic abstract
Planets form in protoplanetary discs. Their masses, distribution, and orbits sensitively depend on the structure of the protoplanetary discs. However, what sets the initial structure of the discs in terms of mass, radius and accretion rate is still unknown. We perform non-ideal MHD numerical simulations using the adaptive mesh refinement code Ramses, of a collapsing, one solar mass, molecular core to study the disc formation and early, up to 100 kyr, evolution, paying great attention to the impact of numerical resolution and accretion scheme. We found that while the mass of the central object is almost independent of the numerical parameters such as the resolution and the accretion scheme onto the sink particle, the disc mass, and to a lower extent its size, heavily depend on the accretion scheme, which we found, is itself resolution dependent. This implies that the accretion onto the star and through the disc are largely decoupled. For a relatively large domain of initial conditions (except at low magnetisation), we found that the properties of the disc do not change too significantly. In particular both the level of initial rotation and turbulence do not influence the disc properties provide the core is sufficiently magnetized. After a short relaxation phase, the disc settles in a stationary state. It then slowly grows in size but not in mass. The disc itself is weakly magnetized but its immediate surrounding is on the contrary highly magnetized. Our results show that the disc properties directly depend on the inner boundary condition, i.e. the accretion scheme onto the central object, suggesting that the disc mass is eventually controlled by the small scale accretion process, possibly the star-disc interaction. Because of ambipolar diffusion and its significant resistivity, the disc diversity remains limited and except for low magnetisation, their properties are (abridged).
The formation of protoplanetary discs during the collapse of molecular dense cores is significantly influenced by angular momentum transport, notably by the magnetic torque. In turn, the evolution of the magnetic field is determined by dynamical processes and non-ideal MHD effects such as ambipolar diffusion. Considering simple relations between various timescales characteristic of the magnetized collapse, we derive an expression for the early disc radius, $ r simeq 18 , {rm AU} , left({eta_{rm AD} / 0.1 , {rm s}} right)^{2/9} left({B_z / 0.1, {rm G}} right) ^{-4/9} left({M / 0.1 msol} right) ^{1/3},$ where $M$ is the total disc plus protostar mass, $eta_mathrm{AD}$ is the ambipolar diffusion coefficient and $B_z$ is the magnetic field in the inner part of the core. This is about significantly smaller than the discs that would form if angular momentum was conserved. The analytical predictions are confronted against a large sample of 3D, non-ideal MHD collapse calculations covering variations of a factor 100 in core mass, a factor 10 in the level of turbulence, a factor 5 in rotation, and magnetic mass-to-flux over critical mass-to-flux ratios 2 and 5. The disc radius estimates are found to agree with the numerical simulations within less than a factor 2. A striking prediction of our analysis is the weak dependence of circumstellar disc radii upon the various relevant quantities, suggesting weak variations among class-0 disc sizes. In some cases, we note the onset of large spiral arms beyond this radius.
I use volume- and mass-limited subsamples and recently published data from the Spitzer Survey of Stellar Structure in Galaxies (S4G) to investigate how the size of bars depends on galaxy properties. The known correlation between bar semi-major-axis $a$ and galaxy stellar mass (or luminosity) is actually *bimodal*: for $log M_{star} < 10.1$, bar size is almost independent of stellar mass ($a propto M_{star}^{0.1}$), while it is a strong function for higher masses ($a propto M_{star}^{0.6}$). Bar size is a slightly stronger function of galaxy half-light radius $r_{e}$ and (especially) exponential disc scale length $h$ ($a propto h^{0.8}$). Correlations between stellar mass and galaxy size can explain the bar-size--$M_{star}$ correlation -- but only for galaxies with $log M_{star} < 10.1$; at higher masses, there is an extra dependence of bar size on $M_{star}$ itself. Despite theoretical arguments that the presence of gas can affect bar growth, there is no evidence for any residual dependence of bar size on (present-day) gas mass fraction. The traditional dependence of bar size on Hubble type (longer bars in early-type discs) can be explained as a side-effect of stellar-mass--Hubble-type correlations. Finally, I show that galaxy size ($r_{e}$ or $h$) can be modeled as a function of stellar mass and both bar presence and bar size: barred galaxies tend to be more extended than unbarred galaxies of the same mass, with larger bars correlated with larger sizes.
Many stars form in regions of enhanced stellar density, wherein the influence of stellar neighbours can have a strong influence on a protoplanetary disc (PPD) population. In particular, far ultraviolet (FUV) flux from massive stars drives thermal winds from the outer edge of PPDs, accelerating disc destruction. In this work, we present a novel technique for constraining the dynamical history of a star forming environment using PPD properties in a strongly FUV irradiated environment. Applying recent models for FUV induced mass loss rates to the PPD population of Cygnus OB2, we constrain how long ago primordial gas was expelled from the region; $ 0.5$ Myr ago if the Shakura & Sunyaev $alpha$-viscosity parameter is $alpha = 10^{-2}$ (corresponding to a viscous timescale of $tau_mathrm{visc} approx 0.5$ Myr for a disc of scale radius $40$ au around a $1, M_odot$ star). This value of $alpha$ is effectively an upper limit, since it assumes efficient extinction of FUV photons throughout the embedded phase. With this gas expulsion timescale we are able to produce a full dynamical model that fits kinematic and morphological data as well as disc fractions. We suggest Cygnus OB2 was originally composed of distinct massive clumps or filaments, each with a stellar mass $sim 10^4 , M_odot$. Finally we predict that in regions of efficient FUV induced mass loss, disc mass $M_mathrm{disc}$ as a function of stellar host mass $m_mathrm{star}$ follows a power law with $M_mathrm{disc} propto m_mathrm{star}^beta$, where $beta gtrsim 2.7$ (depending on disc initial conditions and FUV exposure). This is steeper than observed correlations in regions of moderate FUV flux ($1 < beta <1.9$), and offers a promising diagnostic to establish the influence of external photoevaporation in a given region.
The inner boundary of protoplanetary discs is structured by the dramatic opacity changes at the transition from the dust-containing to a dust-free zone. This paper explores the variety and limits of inner rim structures in passively heated dusty discs. For this study, we implemented detailed sublimation physics in a fast Monte Carlo radiative transfer code. We show that the inner rim in dusty discs is not an infinitely sharp wall but a diffuse region which may be narrow or wide. Furthermore, high surface densities and large silicate grains as well as iron and corundum grains decrease the rim radius, from a 2.2AU radius for small silicates around a 47 Solar luminosity Herbig Ae star typically to 0.4AU and as close as 0.2AU. A passive disc with grain growth and a diverse dust composition must thus have a small inner rim radius. Finally, an analytical expression is presented for the rim location as a function of dust, disc and stellar properties.
The metallicity and its relationship with other galactic properties is a fundamental probe of the evolution of galaxies. In this work, we select about 750,000 star-forming spatial pixels from 1122 blue galaxies in the MaNGA survey to investigate the global stellar mass - local stellar mass surface density - gas-phase metallicity ($M_*$ - $Sigma_*$ - $Z$ ) relation. At a fixed $M_*$, the metallicity increases steeply with increasing $Sigma_*$. Similarly, at a fixed $Sigma_*$, the metallicity increases strongly with increasing $M_*$ at low mass end, while this trend becomes less obvious at high mass end. We find the metallicity to be more strongly correlated to $Sigma_*$ than to $M_*$. Furthermore, we construct a tight (0.07 dex scatter) $M_*$ - $Sigma_*$ - $Z$ relation, which reduces the scatter in the $Sigma_*$ - $Z$ relation by about 30$%$ for galaxies with $7.8 < {rm log}(M_*/M_odot) < 11.0$, while the reduction of scatter is much weaker for high-mass galaxies. This result suggests that, especially for low-mass galaxies, the $M_*$ - $Sigma_*$ - $Z$ relation is largely more fundamental than the $M_*$ - $Z$ and $Sigma_*$ - $Z$ relations, meaning that both $M_*$ and $Sigma_*$ play important roles in shaping the local metallicity. We also find that the local metallicity is probably independent on the local star formation rate surface density at a fixed $M_*$ and $Sigma_*$. Our results are consistent with the scenario that the local metallicities in galaxies are shaped by the combination of the local stars formed in the history and the metal loss caused by galactic winds.