An official website of the United States government
ABSTRACT We study the growth of stellar discs of Milky Way-sized galaxies using a suite of cosmological simulations. We calculate the half-mass axis lengths and axis ratios of stellar populations split by age in galaxies with stellar mass $$M_{*}=10^7\!-\!10^{10}\, \mathrm{M}_{\odot }$$ at redshifts z > 1.5. We find that in our simulations stars always form in relatively thin discs, and at ages below 100 Myr are contained within half-mass height z1/2 ∼ 0.1 kpc and short-to-long axial ratio z1/2/x1/2 ∼ 0.15. Disc thickness increases with the age of stellar population, reaching median z1/2 ∼ 0.8 kpc and z1/2/x1/2 ∼ 0.6 for stars older than 500 Myr. We trace the same group of stars over the simulation snapshots and show explicitly that their intrinsic shape grows more spheroidal over time. We identify a new mechanism that contributes to the observed disc thickness: rapid changes in the orientation of the galactic plane mix the configuration of young stars. The frequently mentioned ‘upside-down’ formation scenario of galactic discs, which posits that young stars form in already thick discs at high redshift, may be missing this additional mechanism of quick disc inflation. The actual formation of stars within a fairly thin plane is consistent with the correspondingly flat configuration of dense molecular gas that fuels star formation.
more »
« less
3D MHD simulations of accretion on to stars with tilted magnetic and rotational axes
ABSTRACT We present results of global 3D magnetohydrodynamic simulations of accretion on to magnetized stars where both the magnetic and rotational axes of the star are tilted about the rotational axis of the disc. We observed that initially the inner parts of the disc are warped, tilted, and precess due to the magnetic interaction between the magnetosphere and the disc. Later, larger tilted discs form with the size increasing with the magnetic moment of the star. The normal vector to the discs are tilted at different angles, from ∼5°–10° up to ∼30°–40°. Small tilts may result from the winding of the magnetic field lines about the rotational axis of the star and the action of the magnetic force which tends to align the disc. Another possible explanation is the magnetic Bardeen–Petterson effect in which the disc settles in the equatorial plane of the star due to precessional and viscous torques in the disc. Tilted discs slowly precess with the time-scale of the order of ∼50 Keplerian periods at the reference radius (∼3 stellar radii). Our results can be applied to different types of stars where signs of tilted discs and/or slow precession have been observed.
more »
« less
- Award ID(s):
- 2009820
- PAR ID:
- 10346542
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 506
- Issue:
- 1
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 372 to 384
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)ABSTRACT Luminous active galactic nuclei and X-ray binaries often contain geometrically thin, radiatively cooled accretion discs. According to theory, these are – in many cases – initially highly misaligned with the black hole equator. In this work, we present the first general relativistic magnetohydrodynamic simulations of very thin (h/r ∼ 0.015–0.05) accretion discs around rapidly spinning (a ∼ 0.9) black holes and tilted by 45°–65°. We show that the inner regions of the discs with h/r ≲ 0.03 align with the black hole equator, though out to smaller radii than predicted by analytic work. The inner aligned and outer misaligned disc regions are separated by a sharp break in tilt angle accompanied by a sharp drop in density. We find that frame dragging by the spinning black hole overpowers the disc viscosity, which is self-consistently produced by magnetized turbulence, tearing the disc apart and forming a rapidly precessing inner sub-disc surrounded by a slowly precessing outer sub-disc. We find that the system produces a pair of relativistic jets for all initial tilt values. At small distances, the black hole launched jets precess rapidly together with the inner sub-disc, whereas at large distances they partially align with the outer sub-disc and precess more slowly. If the tearing radius can be modeled accurately in future work, emission model independent measurements of black hole spin based on precession-driven quasi-periodic oscillations may become possible.more » « less
-
ABSTRACT We investigate the origin of photometric variability in the classical T Tauri star TW Hya by comparing light curves obtained by Transiting Exoplanet Survey Satellite (TESS) and ground-based telescopes with light curves created using three-dimensional (3D) magnetohydrodynamic (MHD) simulations. TW Hya is modelled as a rotating star with a dipole magnetic moment, which is slightly tilted about the rotational axis. We observed that for various model parameters, matter accretes in the unstable regime and produces multiple hotspots on the star’s surface, which leads to stochastic-looking light curves similar to the observed ones. Wavelet and Fourier spectra of observed and modelled light curves show multiple quasi-periodic oscillations (QPOs) with quasi-periods from less than 0.1 to 9 d. Models show that variation in the strength and tilt of the dipole magnetosphere leads to different periodograms, where the period of the star may dominate or be hidden. The amplitude of QPOs associated with the stellar period can be smaller than that of other QPOs if the tilt of the dipole magnetosphere is small and when the unstable regime is stronger. In models with small magnetospheres, the short-period QPOs associated with rotation of the inner disc dominate and can be mistaken for a stellar period. We show that longer period (5–9 d) QPOs can be caused by waves forming beyond the corotation radius.more » « less
-
ABSTRACT Stars form from the gravitational collapse of turbulent, magnetized molecular cloud cores. Our non-ideal MHD simulations reveal that the intrinsically anisotropic magnetic resistance to gravity during the core collapse naturally generates dense gravomagneto sheetlets within inner protostellar envelopes – disrupted versions of classical sheet-like pseudo-discs. They are embedded in a magnetically dominant background, where less dense materials flow along the local magnetic field lines and accumulate in the dense sheetlets. The sheetlets, which feed the disc predominantly through its upper and lower surfaces, are the primary channels for mass and angular momentum transfer from the envelope to the disc. The protostellar disc inherits a small fraction (up to 10 per cent) of the magnetic flux from the envelope, resulting in a disc-averaged net vertical field strength of 1–10 mG and a somewhat stronger toroidal field, potentially detectable through ALMA Zeeman observations. The inherited magnetic field from the envelope plays a dominant role in disc angular momentum evolution, enabling the formation of gravitationally stable discs in cases where the disc field is relatively well-coupled to the gas. Its influence remains significant even in marginally gravitationally unstable discs formed in the more magnetically diffusive cases, removing angular momentum at a rate comparable to or greater than that caused by spiral arms. The magnetically driven disc evolution is consistent with the apparent scarcity of prominent spirals capable of driving rapid accretion in deeply embedded protostellar discs. The dense gravomagneto sheetlets observed in our simulations may correspond to the ‘accretion streamers’ increasingly detected around protostars.more » « less
-
ABSTRACT We describe the first grid-based simulations of the polar alignment of a circumbinary disc. We simulate the evolution of an inclined disc around an eccentric binary using the grid-based code athena++ . The use of a grid-based numerical code allows us to explore lower disc viscosities than have been examined in previous studies. We find that the disc aligns to a polar orientation when the α viscosity is high, while discs with lower viscosity nodally precess with little alignment over 1000 binary orbital periods. The time-scales for polar alignment and disc precession are compared as a function of disc viscosity, and are found to be in agreement with previous studies. At very low disc viscosities (e.g. α = 10−5), anticyclonic vortices are observed along the inner edge of the disc. These vortices can persist for thousands of binary orbits, creating azimuthally localized overdensities and multiple pairs of spiral arms. The vortex is formed at ∼3–4 times the binary semimajor axis, close to the inner edge of the disc, and orbits at roughly the local Keplerian speed. The presence of a vortex in the disc may play an important role in the evolution of circumbinary systems, such as driving episodic accretion and accelerating the formation of polar circumbinary planets.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1093/mnras/stab1724
