The statistical properties of protostellar discs and their dependence on metallicity

We present the analysis of the properties of large samples of protostellar discs formed in four radiation hydrodynamical simulations of star cluster formation. The four calculations have metallicities of 0.01, 0.1, 1 and 3 times solar metallicity. The calculations treat dust and gas temperatures separately and include a thermochemical model of the diffuse interstellar medium. We find the radii of discs of bound protostellar systems tend to decrease with decreasing metallicity, with the median characteristic radius of discs in the 0.01 and 3 times solar metallicity calculations being ≈20 and ≈65 au, respectively. Disc masses and radii of isolated protostars also tend to decrease with decreasing metallicity. We find that the circumstellar discs and orbits of bound protostellar pairs, and the two spins of the two protostars are all less well aligned with each other with lower metallicity than with higher metallicity. These variations with metallicity are due to increased small scale fragmentation due to lower opacities and greater cooling rates with lower metallicity, which increase the stellar multiplicity and increase dynamical interactions. We compare the disc masses and radii of protostellar systems from the solar metallicity calculation with recent surveys of discs around Class 0 and I objects in the Orion and Perseus star-forming regions. The masses and radii of the simulated discs have similar distributions to the observed Class 0 and I discs.

Do we need non-ideal MHD to model protostellar discs?

We investigate and discuss protostellar discs in terms of where the various non-ideal magnetohydrodynamics (MHD) processes are important. We find that the traditional picture of a magnetised disc (where Ohmic resistivity is dominant near the mid-plane, surrounded by a region dominated by the Hall effect, with the remainder of the disc dominated by ambipolar diffusion) is a great oversimplification. In simple parameterised discs, we find that the Hall effect is typically the dominant term throughout the majority of the disc. More importantly, we find that in much of our parameterised discs, at least two non-ideal processes have coefficients within a factor of 10 of one another, indicating that both are important and that naming a dominant term underplays the importance of the other terms. Discs that were self-consistently formed in our previous studies are also dominated by the Hall effect, and the ratio of ambipolar diffusion and Hall coefficients is typically less than 10, suggesting that both terms are equally important and listing a dominant term is misleading. These conclusions become more robust once the magnetic field geometry is taken into account. In agreement with the literature we review, we conclude that non-ideal MHD processes are important for the formation and evolution of protostellar discs. Ignoring any of the non-ideal processes, especially ambipolar diffusion and the Hall effect, yields an incorrect description of disc evolution.

Inferring (sub)millimetre dust opacities and temperature structure in edge-on protostellar discs from resolved multiwavelength continuum observations: the case of the HH 212 disc

ABSTRACT (Sub)millimetre dust opacities are required for converting the observable dust continuum emission to the mass, but their values have long been uncertain, especially in discs around young stellar objects. We propose a method to constrain the opacity κν in edge-on discs from a characteristic optical depth τ0,ν, the density ρ0, and radius R0 at the disc outer edge through κν = τ0,ν/(ρ0R0), where τ0,ν is inferred from the shape of the observed flux along the major axis, ρ0 from gravitational stability considerations, and R0 from direct imaging. We applied the 1D semi-analytical model to the embedded, Class 0, HH 212 disc, which has high-resolution data in Atacama Large Millimetre/submillimetre Array (ALMA) bands 9, 7, 6, and 3 and Very Large Array Ka band (λ = 0.43, 0.85, 1.3, 2.9, and 9.1 mm). The modelling is extended to 2D through RADMC-3D radiative transfer calculations. We find a dust opacity of κν ≈ 1.9 × 10−2, 1.3 × 10−2, and 4.9 × 10−3 cm2 g−1 of gas and dust for ALMA bands 7, 6, and 3, respectively, with uncertainties dependent on the adopted stellar mass. The inferred opacities lend support to the widely used prescription κλ = 2.3 × 10−2(1.3mm/λ) cm2 g−1 . We inferred a temperature of ∼45 K at the disc outer edge that increases radially inwards. It is well above the sublimation temperatures of ices such as CO and N2, which supports the notion that the disc chemistry cannot be completely inherited from the protostellar envelope.

Donald Lynden-Bell. 5 April 1935— 6 February 2018

Donald Lynden-Bell's many contributions to astrophysics encompass general relativity, galactic dynamics, telescope design and observational astronomy. In the 1960s, his papers on stellar dynamics led to fundamental insights into the equilibria of elliptical galaxies, the growth of spiral patterns in disc galaxies and the stability of differentially rotating, self-gravitating flows. Donald introduced the ideas of ‘violent relaxation’ and ‘the gravothermal catastrophe’ in pioneering work on the thermodynamics of galaxies and negative heat capacities. He shared the inaugural Kavli Prize in Astrophysics in 2008 for his contributions to our understanding of quasars. His prediction that dead quasars or supermassive black holes may reside in the nuclei of nearby galaxies has been confirmed by multiple pieces of independent evidence. His work on accretion discs led to new insights into their workings, as well as the realization that the infrared excess in T Tauri stars was caused by protostellar discs around these young stars. He introduced the influential idea of monolithic collapse of a gas cloud as a formation mechanism for the Milky Way Galaxy. As this gave way to modern ideas of merging and accretion as drivers of galaxy formation, Donald was the first to realize the importance of tidal streams as measures of the past history and present-day gravity field of the Galaxy. Though primarily a theorist, Donald participated in one of the first observational programmes to measure the large-scale streaming of nearby galaxies. This led to the discovery of the ‘Great Attractor’. The depth and versatility of his contributions mark Donald out as one of the most influential and pre-eminent astronomers of his day.

Magnetic reconnection in partially ionized plasmas

Magnetic reconnection has been intensively studied in fully ionized plasmas. However, plasmas are often partially ionized in astrophysical environments. The interactions between the neutral particles and ionized plasmas might strongly affect the reconnection mechanisms. We review magnetic reconnection in partially ionized plasmas in different environments from theoretical, numerical, observational and experimental points of view. We focus on mechanisms which make magnetic reconnection fast enough to compare with observations, especially on the reconnection events in the low solar atmosphere. The heating mechanisms and the related observational evidence of the reconnection process in the partially ionized low solar atmosphere are also discussed. We describe magnetic reconnection in weakly ionized astrophysical environments, including the interstellar medium and protostellar discs. We present recent achievements about fast reconnection in laboratory experiments for partially ionized plasmas.

On the early evolution of massive star clusters: the case of cloud D1 and its embedded cluster in NGC 5253

ABSTRACT We discuss a theoretical model for the early evolution of massive star clusters and confront it with the ALMA, radio, and infrared observations of the young stellar cluster highly obscured by the molecular cloud D1 in the nearby dwarf spheroidal galaxy NGC 5253. We show that a large turbulent pressure in the central zones of D1 cluster may cause individual wind-blown bubbles to reach pressure confinement before encountering their neighbours. In this case, stellar winds energy is added to the hot shocked wind pockets of gas around individual massive stars that leads them to meet and produce a cluster wind in time-scales less than 105 yr. In order to inhibit the possibility of cloud dispersal, or the early negative star formation feedback, one should account for mass loading that may come, for example, from pre-main-sequence (PMS) low-mass stars through photoevaporation of their protostellar discs. Mass loading at a rate in excess of 8 × 10−9 M⊙ yr−1 per each PMS star is required to extend the hidden star cluster phase in this particular cluster. In this regime, the parental cloud remains relatively unperturbed, while pockets of molecular, photoionized and hot gas coexist within the star-forming region. Nevertheless, the most likely scenario for cloud D1 and its embedded cluster is that the hot shocked winds around individual massive stars should merge at an age of a few million of years when the PMS star protostellar discs vanish and mass loading ceases that allows a cluster to form a global wind.

Planet formation around M dwarfs via disc instability

Context. Around 30 per cent of the observed exoplanets that orbit M dwarf stars are gas giants that are more massive than Jupiter. These planets are prime candidates for formation by disc instability. Aims. We want to determine the conditions for disc fragmentation around M dwarfs and the properties of the planets that are formed by disc instability. Methods. We performed hydrodynamic simulations of M dwarf protostellar discs in order to determine the minimum disc mass required for gravitational fragmentation to occur. Different stellar masses, disc radii, and metallicities were considered. The mass of each protostellar disc was steadily increased until the disc fragmented and a protoplanet was formed. Results. We find that a disc-to-star mass ratio between ~0.3 and ~0.6 is required for fragmentation to happen. The minimum mass at which a disc fragment increases with the stellar mass and the disc size. Metallicity does not significantly affect the minimum disc fragmentation mass but high metallicity may suppress fragmentation. Protoplanets form quickly (within a few thousand years) at distances around ~50 AU from the host star, and they are initially very hot; their centres have temperatures similar to the ones expected at the accretion shocks around planets formed by core accretion (up to 12 000 K). The final properties of these planets (e.g. mass and orbital radius) are determined through long-term disc-planet or planet–planet interactions. Conclusions. Disc instability is a plausible way to form gas giant planets around M dwarfs provided that discs have at least 30% the mass of their host stars during the initial stages of their formation. Future observations of massive M dwarf discs or planets around very young M dwarfs are required to establish the importance of disc instability for planet formation around low-mass stars.

Magnetorotational instability in diamagnetic, misaligned protostellar discs

ABSTRACT In this study, we addressed the question of how the growth rate of the magnetorotational instability is modified when the radial component of the stellar dipole magnetic field is taken into account in addition to the vertical component. Considering a fiducial radius in the disc where diamagnetic currents are pronounced, we carried out a linear stability analysis to obtain the growth rates of the magnetorotational instability for various parameters such as the ratio of the radial-to-vertical component and the gradient of the magnetic field, the Alfvenic Mach number, and the diamagnetization parameter. Our results show that the interaction between the diamagnetic current and the radial component of the magnetic field increases the growth rate of the magnetorotational instability and generates a force perpendicular to the disc plane that may induce a torque. It is also shown that considering the radial component of the magnetic field and taking into account a radial gradient in the vertical component of the magnetic field causes an increase in the magnitudes of the growth rates of both the axisymmetric (m = 0) and the non-axisymmetric (m = 1) modes.

Spiral-arm instability – III. Fragmentation of primordial protostellar discs

ABSTRACT We study the gravitational instability and fragmentation of primordial protostellar discs by using high-resolution cosmological hydrodynamics simulations. We follow the formation and evolution of spiral arms in protostellar discs, examine the dynamical stability, and identify a physical mechanism of secondary protostar formation. We use linear perturbation theory based on the spiral-arm instability (SAI) analysis in our previous studies. We improve the analysis by incorporating the effects of finite thickness and shearing motion of arms, and derive the physical conditions for SAI in protostellar discs. Our analysis predicts accurately the stability and the onset of arm fragmentation that is determined by the balance between self-gravity and gas pressure plus the Coriolis force. Formation of secondary and multiple protostars in the discs is explained by the SAI, which is driven by self-gravity and thus can operate without rapid gas cooling. We can also predict the typical mass of the fragments, which is found to be in good agreement with the actual masses of secondary protostars formed in the simulation.

There is no magnetic braking catastrophe: low-mass star cluster and protostellar disc formation with non-ideal magnetohydrodynamics

We present results from the first radiation non-ideal magnetohydrodynamics (MHD) simulations of low-mass star cluster formation that resolve the fragmentation process down to the opacity limit. We model 50 M⊙ turbulent clouds initially threaded by a uniform magnetic field with strengths of 3, 5 10, and 20 times the critical mass-to-magnetic flux ratio, and at each strength, we model both an ideal and non-ideal (including Ohmic resistivity, ambipolar diffusion, and the Hall effect) MHD cloud. Turbulence and magnetic fields shape the large-scale structure of the cloud, and similar structures form regardless of whether ideal or non-ideal MHD is employed. At high densities (106 ≲ nH ≲ 1011 cm−3), all models have a similar magnetic field strength versus density relation, suggesting that the field strength in dense cores is independent of the large-scale environment. Albeit with limited statistics, we find no evidence for the dependence of the initial mass function on the initial magnetic field strength, however, the star formation rate decreases for models with increasing initial field strengths; the exception is the strongest field case where collapse occurs primarily along field lines. Protostellar discs with radii ≳ 20 au form in all models, suggesting that disc formation is dependent on the gas turbulence rather than on magnetic field strength. We find no evidence for the magnetic braking catastrophe, and find that magnetic fields do not hinder the formation of protostellar discs.