Efficient Magnetic Braking Research Articles

The effects of surface fossil magnetic fields on massive star evolution: V. Models at low metallicity

Abstract At metallicities lower than that of the Small Magellanic Cloud, it remains essentially unexplored how fossil magnetic fields, forming large-scale magnetospheres, could affect the evolution of massive stars, thereby impacting the fundamental building blocks of the early Universe. We extend our stellar evolution model grid with representative calculations of main-sequence, single-star models with initial masses of 20 and 60 M⊙, including appropriate changes for low-metallicity environments (Z = 10−3–10−6). We scrutinise the magnetic, rotational, and chemical properties of the models. When lowering the metallicity, the rotational velocities can become higher and the tendency towards quasi-chemically homogeneous evolution increases. While magnetic fields aim to prevent the development of this evolutionary channel, the weakening stellar winds lead to less efficient magnetic braking in our models. Since the stellar radius is almost constant during a blueward evolution caused by efficient chemical mixing, the surface magnetic field strength remains unchanged in some models. We find core masses at the terminal-age main sequence between 22 and 52 M⊙ for initially 60 M⊙ models. This large difference is due to the vastly different chemical and rotational evolution. We conclude that in order to explain chemical species and, in particular, high nitrogen abundances in the early Universe, the adopted stellar models need to be under scrutiny. The assumptions regarding wind physics, chemical mixing, and magnetic fields will strongly impact the model predictions.

Constraints on Magnetic Braking from the G8 Dwarf Stars 61 UMa and τ Cet

During the first half of their main-sequence lifetimes, stars rapidly lose angular momentum to their magnetized winds, a process known as magnetic braking. Recent observations suggest a substantial decrease in the magnetic braking efficiency when stars reach a critical value of the Rossby number, the stellar rotation period normalized by the convective overturn timescale. Cooler stars have deeper convection zones with longer overturn times, reaching this critical Rossby number at slower rotation rates. The nature and timing of the transition to weakened magnetic braking have previously been constrained by several solar analogs and two slightly hotter stars. In this Letter, we derive the first direct constraints from stars cooler than the Sun. We present new spectropolarimetry of the old G8 dwarf τ Cet from the Large Binocular Telescope, and we reanalyze a published Zeeman Doppler image of the younger G8 star 61 UMa, yielding the large-scale magnetic field strengths and morphologies. We estimate mass-loss rates using archival X-ray observations and inferences from Lyα measurements, and we adopt other stellar properties from asteroseismology and spectral energy distribution fitting. The resulting calculations of the wind braking torque demonstrate that the rate of angular momentum loss drops by a factor of 300 between the ages of these two stars (1.4–9 Gyr), well above theoretical expectations. We summarize the available data to help constrain the value of the critical Rossby number, and we identify a new signature of the long-period detection edge in recent measurements from the Kepler mission.

Modeling disks and magnetic outflows around a forming massive star

Context. Similar to their lower mass siblings, massive protostars can be expected to: (a) be surrounded by circumstellar disks, and (b) launch magnetically driven jets and outflows. The disk formation and global evolution is thereby controlled by advection of angular momentum from large scales, the efficiency of magnetic braking and the resistivity of the medium, and the internal thermal and magnetic pressures of the disk. Aims. We determine the dominant physical mechanisms that shape the appearance of these circumstellar disks, their sizes, and aspect ratios. Methods. We performed a series of 30 simulations of a massive star forming from the gravitational collapse of a molecular cloud threaded by an initially uniform magnetic field, starting from different values for the mass of the cloud, its initial density and rotation profiles, its rotational energy content, the magnetic field strength, and the resistivity of the material. The gas and dust was modeled with the methods of resistive magnetohydrodynamics, also considering radiation transport of thermal emission and self-gravity. We checked for the impact of spatial resolution in a dedicated convergence study. Results. After the initial infall phase dominated by the gravitational collapse, an accretion disk was formed, shortly followed by the launching of magnetically driven outflows. Two layers can be distinguished in the accretion disk: a thin layer, vertically supported by thermal pressure, and a thick layer, vertically supported by magnetic pressure. Both regions exhibit Keplerian-like rotation and grow outward over time. We observed the effects of magnetic braking in the inner ~50 au of the disk at late times in our fiducial case. The parameter study reveals that the size of the disk is mostly determined by the density and rotation profiles of the initial mass reservoir and not by the magnetic field strength. We find that the disk size and protostellar mass gain scale with the initial mass of the cloud. Magnetic pressure can slightly increase the size of the accretion disk, while magnetic braking is more relevant in the innermost parts of the disk as opposed to the outer disk. From the parameter study, we infer that multiple initial conditions for the onset of gravitational collapse are able to produce a given disk size and protostellar mass.

Magnetised winds in transition discs

Context. Protoplanetary discs are cold, dense, and weakly ionised environments that witness planetary formation. Among these discs, transition discs (TDs) are characterised by a wide cavity (up to tens of au) in the dust and gas distribution. Despite this lack of material, a considerable fraction of TDs are still strongly accreting onto their central star, possibly indicating that a mechanism is driving fast accretion in TD cavities. Aims. The presence of radially extended ‘dead zones’ in protoplanetary discs has recently revived interest in magnetised disc winds (MDWs), where accretion is driven by a large magnetic field extracting angular momentum from the disc. We propose that TDs could be subject to similar disc winds, and that these could naturally explain the fast-accreting and long-lived cavities inferred in TDs. Methods. We present the results of the first 2.5D global numerical simulations of TDs harbouring MDWs using the PLUTO code. We imposed a cavity in the gas distribution with various density contrasts, and considered a power-law distribution for the large-scale magnetic field strength. We assume the disc is weakly ionised and is therefore subject to ambipolar diffusion, as expected in this range of densities and temperatures. Results. We find that our simulated TDs always reach a steady state with an inner cavity and an outer ‘standard’ disc. These models also maintain an approximately constant accretion rate through the entire structure, reaching 10−7 M⊚ yr−1 for typical surface density values. The MDW launched from the cavity is more magnetised and has a significantly larger lever arm (up to 10) than the MDW launched from the outer disc. The material in the cavity is accreted at sonic velocities, and the cavity itself is rotating at 70% of the Keplerian velocity due to the efficient magnetic braking imposed by the MDW. Overall, our cavity matches the dynamical properties of an inner jet emitting disc (JED) and of magnetically arrested discs (MADs) in black-hole physics. Finally, we observe that the cavity is subject to recurring accretion bursts that may be driven by a magnetic Rayleigh-Taylor instability of the cavity edge. Conclusions. Some strongly accreting TDs could be the result of magnetised wind sculpting protoplanetary discs. Kinematic diagnostics of the disc or the wind (orbital velocity, wind speeds, accretion velocities) could disentangle classic photo-evaporation from MDW models.

Impact of magnetic braking on high-mass close binary formation

ABSTRACT Combining numerical simulations and analytical modelling, we investigate whether close binary systems form by the effect of magnetic braking. Using magnetohydrodynamics simulations, we calculate the cloud evolution with a sink, for which we do not resolve the binary system or binary orbital motion to realize long-term time integration. Then, we analytically estimate the binary separation using the accreted mass and angular momentum obtained from the simulation. In unmagnetized clouds, wide binary systems with separations of >100 au form, in which the binary separation continues to increase during the main accretion phase. In contrast, close binary systems with separations of <100 au can form in magnetized clouds. Since the efficiency of magnetic braking strongly depends on both the strength and configuration of the magnetic field, they also affect the formation conditions of a close binary. In addition, the protostellar outflow has a negative impact on close binary formation, especially when the rotation axis of the prestellar cloud is aligned with the global magnetic field. The outflow interrupts the accretion of gas with small angular momentum, which is expelled from the cloud, while gas with large angular momentum preferentially falls from the side of the outflow on to the binary system and widens the binary separation. This study shows that a cloud with a magnetic field that is not parallel to the rotation axis is a favourable environment for the formation of close binary systems.

The effects of surface fossil magnetic fields on massive star evolution – III. The case of τ Sco

ABSTRACT τ Sco, a well-studied magnetic B-type star in the Upper Sco association, has a number of surprising characteristics. It rotates very slowly and shows nitrogen excess. Its surface magnetic field is much more complex than a purely dipolar configuration which is unusual for a magnetic massive star. We employ the cmfgen radiative transfer code to determine the fundamental parameters and surface CNO and helium abundances. Then, we employ mesa and genec stellar evolution models accounting for the effects of surface magnetic fields. To reconcile τ Sco’s properties with single-star models, an increase is necessary in the efficiency of rotational mixing by a factor of 3–10 and in the efficiency of magnetic braking by a factor of 10. The spin-down could be explained by assuming a magnetic field decay scenario. However, the simultaneous chemical enrichment challenges the single-star scenario. Previous works indeed suggested a stellar merger origin for τ Sco. However, the merger scenario also faces similar challenges as our magnetic single-star models to explain τ Sco’s simultaneous slow rotation and nitrogen excess. In conclusion, the single-star channel seems less likely and versatile to explain these discrepancies, while the merger scenario and other potential binary-evolution channels still require further assessment as to whether they may self-consistently explain the observables of τ Sco.

The JCMT BISTRO Survey: Alignment between Outflows and Magnetic Fields in Dense Cores/Clumps

Abstract We compare the directions of molecular outflows of 62 low-mass Class 0 and I protostars in nearby (<450 pc) star-forming regions with the mean orientations of the magnetic fields on 0.05–0.5 pc scales in the dense cores/clumps where they are embedded. The magnetic field orientations were measured using the JCMT POL-2 data taken by the BISTRO-1 survey and from the archive. The outflow directions were observed with interferometers in the literature. The observed distribution of the angles between the outflows and the magnetic fields peaks between 15° and 35°. After considering projection effects, our results could suggest that the outflows tend to be misaligned with the magnetic fields by 50° ± 15° in three-dimensional space and are less likely (but not ruled out) randomly oriented with respect to the magnetic fields. There is no correlation between the misalignment and the bolometric temperatures in our sample. In several sources, the small-scale (1000–3000 au) magnetic field is more misaligned with the outflow than the large-scale magnetic field, suggesting that the small-scale magnetic field has been twisted by the dynamics. In comparison with turbulent MHD simulations of core formation, our observational results are more consistent with models in which the energy densities in the magnetic field and the turbulence of the gas are comparable. Our results also suggest that the misalignment alone cannot sufficiently reduce the efficiency of magnetic braking to enable formation of the observed number of large Keplerian disks with sizes larger than 30–50 au.

LBT/PEPSI Spectropolarimetry of a Magnetic Morphology Shift in Old Solar-type Stars*

Abstract Solar-type stars are born with relatively rapid rotation and strong magnetic fields. Through a process known as magnetic braking, the rotation slows over time as stellar winds gradually remove angular momentum from the system. The rate of angular momentum loss depends sensitively on the magnetic morphology, with the dipole field exerting the largest torque on the star. Recent observations suggest that the efficiency of magnetic braking may decrease dramatically in stars near the middle of their main-sequence lifetimes. One hypothesis to explain this reduction in efficiency is a shift in magnetic morphology from predominantly larger to smaller spatial scales. We aim to test this hypothesis with spectropolarimetric measurements of two stars that sample chromospheric activity levels on opposite sides of the proposed magnetic transition. As predicted, the more active star (HD 100180) exhibits a significant circular polarization signature due to a nonaxisymmetric large-scale magnetic field, while the less active star (HD 143761) shows no significant signal. We identify analogs of the two stars among a sample of well-characterized Kepler targets, and we predict that the asteroseismic age of HD 143761 from future Transiting Exoplanet Survey Satellite observations will substantially exceed the age expected from gyrochronology. We conclude that a shift in magnetic morphology likely contributes to the loss of magnetic braking in middle-aged stars, which appears to coincide with the shutdown of their global dynamos.

Kinematics of neutral and ionzied gas in the candidate protostar with efficient magnetic braking B335

AbstractAmbipolar diffusion can cause a velocity drift between ions and neutrals. This is one of the non-ideal MHD effects proposed to enable the formation of large Keplerian disks with sizes of tens of au (Zhao et al. 2018). To observationally study ambipolar diffusion in collapsing protostellar envelopes, we analyzed the ALMA H13CO+ (3–2) and C18O (2–1) data of the protostar B335, which is a candidate source with efficient magnetic braking (Yen et al. 2015). We constructed kinematical models to fit the velocity structures observed in H13CO+ and C18O. With our kinematical models, the infalling velocities in H13CO+ and C18O are both measured to be 0.85 ± 0.2 km s−1 at a radius of 100 au, suggesting that the velocity drift between the ionized and neutral gas is at most 0.3 km s−1 at a radius of 100 au in B335. The Hall parameter for H13CO+ is estimated to be ≫1 on a 100 au scale in B335, so that H13CO+ is expected to be attached to the magnetic field. Our non-detection or upper limit of the velocity drift between the ionized and neutral gas could suggest that the magnetic field remains rather well coupled to the bulk neutral material on a 100 au scale in B335, and that any significant field-matter decoupling, if present, likely occurs only on a smaller scale, leading to an accumulation of magnetic flux and thus efficient magnetic braking in the inner envelope in B335.

Magnetically regulated collapse in the B335 protostar ?

AbstractThe role of the magnetic field during protostellar collapse is still poorly constrained from an observational point of view, and only few constraints exist that shed light on the magnetic braking efficiency during the main accretion phase. I presented our ALMA polarimetric observations of the thermal dust continuum emission at 1.3 mm, towards the B335 Class 0 protostar (Maury et al. 2018a). Linearly polarized dust emission is detected at all scales probed by our observations (50 to 1000 au). The magnetic field structure has a very ordered topology in the inner envelope, with a transition from a large-scale poloidal magnetic field, in the outflow direction, to strongly pinched in the equatorial direction. We compared our data to a family of magnetized protostellar collapse models. We show that only models with an initial core mass-to-flux ratio μ∼5-6 are able to reproduce the observed properties of B335, especially the upper-limits on its disk size, its large-scale envelope rotation β and the pronounced magnetic field lines pinching observed in our ALMA data. In these MHD models, the magnetic field is dynamically relevant to regulate the typical outcome of protostellar collapse, suggesting a magnetically-regulated disk formation scenarios is at work in B335.

Constraint on ion–neutral drift velocity in the Class 0 protostar B335 from ALMA observations

Aims. Ambipolar diffusion can cause a velocity drift between ions and neutrals. This is one of the non-ideal magnetohydrodynamics (MHD) effects proposed to enable the formation of large-scale Keplerian disks with sizes of tens of au. To observationally study ambipolar diffusion in collapsing protostellar envelopes, we compare here gas kinematics traced by ionized and neutral molecular lines and discuss the implication on ambipolar diffusion. Methods. We analyzed the data of the H13CO+ (3–2) and C18O (2–1) emission in the Class 0 protostar B335 obtained with our ALMA observations. We constructed kinematical models to fit the velocity structures observed in the H13CO+ and C18O emission and to measure the infalling velocities of the ionized and neutral gas on a 100 au scale in B335. Results. A central compact (~1′′–2′′) component that is elongated perpendicular to the outflow direction and exhibits a clear velocity gradient along the outflow direction is observed in both lines and most likely traces the infalling flattened envelope. With our kinematical models, the infalling velocities in the H13CO+ and C18O emission are both measured to be 0.85 ± 0.2 km s−1 at a radius of 100 au, suggesting that the velocity drift between the ionized and neutral gas is at most 0.3 km s−1 at a radius of 100 au in B335. Conclusions. The Hall parameter for H13CO+ is estimated to be ≫1 on a 100 au scale in B335, so that H13CO+ is expected to be attached to the magnetic field. Our non-detection or upper limit of the velocity drift between the ionized and neutral gas could suggest that the magnetic field remains rather well coupled to the bulk neutral material on a 100 au scale in this source, and that any significant field-matter decoupling, if present, likely occurs only on a smaller scale, leading to an accumulation of magnetic flux and thus efficient magnetic braking in the inner envelope. This result is consistent with the expectation from the MHD simulations with a typical ambipolar diffusivity and those without ambipolar diffusion. On the other hand, the high ambipolar drift velocity of 0.5–1.0 km s−1 on a 100 au scale predicted in the MHD simulations with an enhanced ambipolar diffusivity by removing small dust grains, where the minimum grain size is 0.1 μm, is not detected in our observations. However, because of our limited angular resolution, we cannot rule out a significant ambipolar drift only in the midplane of the infalling envelope. Future observations with higher angular resolutions (~0. ′′1) are needed to examine this possibility and ambipolar diffusion on a smaller scale.

Decoupling of magnetic fields in collapsing protostellar envelopes and disc formation and fragmentation

Efficient magnetic braking is a formidable obstacle to the formation of rotationally supported disks (RSDs) around protostars in magnetized dense cores. We have previously shown, through 2D (axisymmetric) non-ideal MHD simulations, that removing very small grains (VSGs: ~10 \AA$~$to few 100 \AA) can greatly enhance ambipolar diffusion and enable the formation of RSDs. Here we extend the simulations of disk formation enabled by VSG removal to 3D. We find that the key to this scenario of disk formation is that the drift velocity of the magnetic field almost cancels out the infall velocity of the neutrals in the $10^2$-$10^3$ AU-scale "pseudo-disk" where the field lines are most severely pinched and most of protostellar envelope mass infall occurs. As a result, the bulk neutral envelope matter can collapse without dragging much magnetic flux into the disk-forming region, which lowers the magnetic braking efficiency. We find that the initial disks enabled by VSG removal tend to be Toomre-unstable, which leads to the formation of prominent spiral structures that function as centrifugal barriers. The piling-up of infall material near the centrifugal barrier often produces dense fragments of tens of Jupiter masses, especially in cores that are not too strongly magnetized. Some fragments accrete onto the central stellar object, producing bursts in mass accretion rate. Others are longer lived, although whether they can survive long-term to produce multiple systems remains to be ascertained. Our results highlight the importance of dust grain evolution in determining the formation and properties of protostellar disks and potentially multiple systems.

ANGULAR MOMENTUM LOSS IN THE ENVELOPE–DISK TRANSITION REGION OF THE HH 111 PROTOSTELLAR SYSTEM: EVIDENCE FOR MAGNETIC BRAKING?

ABSTRACT HH 111 is a Class I protostellar system at a distance of ∼400 pc, with the central source VLA 1 associated with a rotating disk deeply embedded in a flattened envelope. Here we present the observations of this system at ∼0.″6 (240 au) resolution in C18O (J = 2 — 1) and a 230 GHz continuum obtained with the Atacama Large Millimeter/Submillimeter Array, and in SO ( ) obtained with the Submillimeter Array. The observations show for the first time how a Keplerian rotating disk can be formed inside a flattened envelope. The flattened envelope is detected in C18O, extending out to ≳2400 au from the VLA 1 source. It has a differential rotation, with the outer part (≳2000 au) better described by a rotation that has constant specific angular momentum, and the innermost part (≲160 au) by a Keplerian rotation. The rotationally supported disk is therefore relatively compact in this system, which is consistent with the dust continuum observations. Most interestingly, if the flow is in steady state, there is a substantial drop in specific angular momentum in the envelope–disk transition region from 2000 to 160 au, by a factor of ∼3. Such a decrease is not expected outside a disk formed from simple hydrodynamic core collapse, but can happen naturally if the core is significantly magnetized, because magnetic fields can be trapped in the transition region outside the disk by the ram pressure of the protostellar accretion flow, which can lead to efficient magnetic braking. In addition, SO shock emission is detected around the outer radius of the disk and could trace an accretion shock around the disk.

STELLAR EVIDENCE THAT THE SOLAR DYNAMO MAY BE IN TRANSITION

ABSTRACT Precise photometry from the Kepler space telescope allows not only the measurement of rotation in solar-type field stars, but also the determination of reliable masses and ages from asteroseismology. These critical data have recently provided the first opportunity to calibrate rotation–age relations for stars older than the Sun. The evolutionary picture that emerges is surprising: beyond middle-age the efficiency of magnetic braking is dramatically reduced, implying a fundamental change in angular momentum loss beyond a critical Rossby number (Ro ∼ 2). We compile published chromospheric activity measurements for the sample of Kepler asteroseismic targets that were used to establish the new rotation–age relations. We use these data along with a sample of well-characterized solar analogs from the Mount Wilson HK survey to develop a qualitative scenario connecting the evolution of chromospheric activity to a fundamental shift in the character of differential rotation. We conclude that the Sun may be in a transitional evolutionary phase, and that its magnetic cycle might represent a special case of stellar dynamo theory.

THE MAGNETIC FIELD IN THE CLASS 0 PROTOSTELLAR DISK OF L1527

We present subarcsecond (~0.35") resolved observations of the 1.3 mm dust polarization from the edge-on circumstellar disk around the Class 0 protostar L1527. The inferred magnetic field is consistent with a dominantly toroidal morphology; there is no significantly detected vertical poloidal component to which observations of an edge-on disk are most sensitive. This suggests that angular momentum transport in Class 0 protostars (when large amounts of material are fed down to the disk from the envelope and accreted onto the protostar) is driven mainly by magnetorotational instability rather than magnetocentrifugal winds at 50 AU scales. In addition, with the data to date there is an early, tentative trend that R>30 AU disks have so far been found in Class 0 systems with average magnetic fields on the 1000 AU scale strongly misaligned with the rotation axis. The absence of such a disk in the aligned case could be due to efficient magnetic braking that disrupts disk formation. If this is the case, this implies that candidate Class 0 disk systems could be identified by the average magnetic field direction at ~1000 AU spatial scales.

Accretion and magnetic field morphology around Class 0 stage protostellar discs

We analyse simulations of turbulent, magnetised molecular cloud cores focussing on the formation of Class 0 stage protostellar discs and the physical conditions in their surroundings. We show that for a wide range of initial conditions Keplerian discs are formed in the Class 0 stage already. In particular, we show that even subsonic turbulent motions reduce the magnetic braking efficiency sufficiently in order to allow rotationally supported discs to form. We therefore suggest that already during the Class 0 stage the fraction of Keplerian discs is significantly higher than 50%, consistent with recent observational trends but significantly higher than predictions based on simulations with misaligned magnetic fields, demonstrating the importance of turbulent motions for the formation of Keplerian discs. We show that the accretion of mass and angular momentum in the surroundings of protostellar discs occurs in a highly anisotropic manner, by means of a few narrow accretion channels. The magnetic field structure in the vicinity of the discs is highly disordered, revealing field reversals up to distances of 1000 AU. These findings demonstrate that as soon as even mild turbulent motions are included, the classical disc formation scenario of a coherently rotating environment and a well-ordered magnetic field breaks down. Hence, it is highly questionable to assess the magnetic braking efficiency based on non-turbulent collapse simulation. We strongly suggest that, in addition to the global magnetic field properties, the small-scale accretion flow and detailed magnetic field structure have to be considered in order to assess the likelihood of Keplerian discs to be present.

Turbulence-induced disc formation in strongly magnetized cloud cores

We present collapse simulations of strongly magnetised, turbulent molecular cloud cores with masses ranging from 2.6 to 1000 M_sun in order to study the influence of the initial conditions on the turbulence-induced disc formation mechanism proposed recently by Seifried et al. 2012. We find that Keplerian discs are formed in all cases independently of the core mass, the strength of turbulence, or the presence of global rotation. The discs appear within a few kyr after the formation of the protostar, are 50 - 150 AU in size, and have masses between 0.05 and a few 0.1 M_sun. During the formation of the discs the mass-to-flux ratio stays well below the critical value of 10 for Keplerian disc formation. Hence, flux-loss alone cannot explain the formation of Keplerian discs. The formation of rotationally supported discs at such early phases is rather due to the disordered magnetic field structure and due to turbulent motions in the surroundings of the discs, two effects lowering the classical magnetic braking efficiency. Binary systems occurring in the discs are mainly formed via the disc capturing mechanism rather than via disc fragmentation, which is largely suppressed by the presence of magnetic fields.

Disc formation in turbulent massive cores: circumventing the magnetic braking catastrophe

ABSTRACT We present collapse simulations of 100 M⊙ turbulent cloud cores threaded by a strong magnetic field. During the initial collapse phase, filaments are generated which fragment quickly and form several protostars. Around these protostars Keplerian discs with typical sizes of up to 100 au build up in contrast to previous simulations neglecting turbulence. We examine three mechanisms potentially responsible for lowering the magnetic braking efficiency and therefore allowing for the formation of Keplerian discs. Analysing the condensations in which the discs form, we show that the build-up of Keplerian discs is neither caused by magnetic flux loss due to turbulent reconnection nor by the misalignment of the magnetic field and the angular momentum. It is rather a consequence of the turbulent surroundings of the disc which exhibit no coherent rotation structure while strong local shear flows carry large amounts of angular momentum. We suggest that the ‘magnetic braking catastrophe’, i.e. the formation of sub-Keplerian discs only, is an artefact of the idealized non-turbulent initial conditions and that turbulence provides a natural mechanism to circumvent this problem.

Magnetic fields during the early stages of massive star formation - I. Accretion and disc evolution

We present simulations of collapsing 100 M_\sun mass cores in the context of massive star formation. The effect of variable initial rotational and magnetic energies on the formation of massive stars is studied in detail. We focus on accretion rates and on the question under which conditions massive Keplerian disks can form in the very early evolutionary stage of massive protostars. For this purpose, we perform 12 simulations with different initial conditions extending over a wide range in parameter space. The equations of magnetohydrodynamics (MHD) are solved under the assumption of ideal MHD. We find that the formation of Keplerian disks in the very early stages is suppressed for a mass-to-flux ratio normalised to the critical value \mu below 10, in agreement with a series of low-mass star formation simulations. This is caused by very efficient magnetic braking resulting in a nearly instantaneous removal of angular momentum from the disk. For weak magnetic fields, corresponding to \mu > 10, large-scale, centrifugally supported disks build up with radii exceeding 100 AU. A stability analysis reveals that the disks are supported against gravitationally induced perturbations by the magnetic field and tend to form single stars rather than multiple objects. We find protostellar accretion rates of the order of a few 10^-4 M_\sun yr^-1 which, considering the large range covered by the initial conditions, vary only by a factor of ~ 3 between the different simulations. We attribute this fact to two competing effects of magnetic fields. On the one hand, magnetic braking enhances accretion by removing angular momentum from the disk thus lowering the centrifugal support against gravity. On the other hand, the combined effect of magnetic pressure and magnetic tension counteracts gravity by exerting an outward directed force on the gas in the disk thus reducing the accretion onto the protostars.

Disk formation during collapse of magnetized protostellar cores

In the context of star and planet formation, understanding the formation of disks is of fundamental importance. Previous studies found that the magnetic field has a very strong impact on the collapse of a prestellar cloud, particularly in possibly suppressing the formation of a disk even for relatively modest values of the magnetic intensity. Since observations infer that cores have a substantial level of magnetization, this raises the question of how disks form. However, most studies have been restricted to the case in which the initial angle, $\alpha$, between the magnetic field and the rotation axis equals 0$^\circ$. We explore and analyse the influence of non aligned configurations on disk formation. We perform 3D ideal MHD, AMR numerical simulations for various values of $\mu$, the ratio of the mass-to-flux to the critical mass-to-flux, and various values of $\alpha$. We find that disks form more easily as $\alpha$ increases from 0 to 90$^\circ$. We propose that as the magnetized pseudo-disks become thicker with increasing $\alpha$, the magnetic braking efficiency is lowered. We also find that even small values of $\alpha$ ($\simeq$ 10-20$^\circ$) show significant differences with the alligned case. Within the framework of ideal MHD and for our choice of initial conditions, centrifugally supported disks cannot form for values of $\mu$ smaller than $\simeq$3, if the magnetic field and the rotation axis are perpendicular, and smaller than about $\simeq$5-10 when they are perfectly aligned.