-
Introducing two improved methods for approximating radiative cooling in hydrodynamical simulations of accretion discs
Authors:
Alison K. Young,
Maggie Celeste,
Richard A. Booth,
Ken Rice,
Adam Koval,
Ethan Carter,
Dimitris Stamatellos
Abstract:
The evolution of many astrophysical systems depends strongly on the balance between heating and cooling, in particular star formation in giant molecular clouds and the evolution of young protostellar systems. Protostellar discs are susceptible to the gravitational instability, which can play a key role in their evolution and in planet formation. The strength of the instability depends on the rate…
▽ More
The evolution of many astrophysical systems depends strongly on the balance between heating and cooling, in particular star formation in giant molecular clouds and the evolution of young protostellar systems. Protostellar discs are susceptible to the gravitational instability, which can play a key role in their evolution and in planet formation. The strength of the instability depends on the rate at which the system loses thermal energy. To study the evolution of these systems, we require radiative cooling approximations because full radiative transfer is generally too expensive to be coupled to hydrodynamical models. Here we present two new approximate methods for computing radiative cooling that make use of the polytropic cooling approximation. This approach invokes the assumption that each parcel of gas is located within a spherical pseudo-cloud which can then be used to approximate the optical depth. The first method combines the methods introduced by Stamatellos et al. and Lombardi et al. to overcome the limitations of each method at low and high optical depths respectively. The second, the "Modified Lombardi" method, is specifically tailored for self-gravitating discs. This modifies the scale height estimate from the method of Lombardi et al. using the analytical scale height for a self-gravitating disc. We show that the Modified Lombardi method provides an excellent approximation for the column density in a fragmenting disc, a regime in which the existing methods fail to recover the clumps and spiral structures. We therefore recommend this improved radiative cooling method for more realistic simulations of self-gravitating discs.
△ Less
Submitted 9 May, 2024;
originally announced May 2024.
-
Insights into the first and second hydrostatic core stages from numerical simulations
Authors:
Alison K. Young
Abstract:
The theory of how low mass stars form from the collapse of a dense molecular cloud core has been well-established for decades. Thanks to significant progress in computing and numerical modelling, more physical models have been developed and a wider parameter space explored to understand the early stages of star formation more fully. In this review, I describe the expected physical properties of th…
▽ More
The theory of how low mass stars form from the collapse of a dense molecular cloud core has been well-established for decades. Thanks to significant progress in computing and numerical modelling, more physical models have been developed and a wider parameter space explored to understand the early stages of star formation more fully. In this review, I describe the expected physical properties of the first and second core stages and how the inclusion of different physics affects those predicted characteristics. I provide an overview of chemical models and synthetic observations, looking towards the positive identification of the first core in nature, which remains elusive. However, there are a few likely candidate first cores, which are listed, and I briefly discuss the recent progress in characterising the youngest protostellar sources. Chemistry will be instrumental in the firm identification of the first core so we require robust theoretical predictions of the chemical evolution of protostellar cores, especially of the first and second core outflows. Looking ahead, simulations can shed light on how the protostellar collapse phase shapes the evolution of the protostellar disc. Simulations of dust evolution during protostellar core collapse show there is significant enhancement in grain size and abundance towards the centre of the core. Chemical models show that the warm, dense conditions of the first core drive chemical evolution. There is a wide scope for further study of the role that the first and second core stages play in determining the structure and composition of the protostellar disc and envelope and, of course, the eventual influence on the formation of planets.
△ Less
Submitted 5 December, 2023;
originally announced December 2023.
-
On the conditions for warping and breaking protoplanetary discs
Authors:
Alison K. Young,
Struan Stevenson,
C. J. Nixon,
Ken Rice
Abstract:
Recent observations demonstrate that misalignments and other out-of-plane structures are common in protoplanetary discs. Many of these have been linked to a central host binary with an orbit that is inclined with respect to the disc. We present simulations of misaligned circumbinary discs with a range of parameters to gain a better understanding of the link between those parameters and the disc mo…
▽ More
Recent observations demonstrate that misalignments and other out-of-plane structures are common in protoplanetary discs. Many of these have been linked to a central host binary with an orbit that is inclined with respect to the disc. We present simulations of misaligned circumbinary discs with a range of parameters to gain a better understanding of the link between those parameters and the disc morphology in the wave-like regime of warp propagation that is appropriate to protoplanetary discs. The simulations confirm that disc tearing is possible in protoplanetary discs as long as the mass ratio, $μ$, and disc-binary inclination angle, $i$, are not too small. For the simulations presented here this corresponds to $μ> 0.1$ and $i \gtrsim 40^\circ$. For highly eccentric binaries, tearing can occur for discs with smaller misalignment. Existing theoretical predictions provide an estimate of the radial extent of the disc in which we can expect breaking to occur. However, there does not seem to be a simple relationship between the disc properties and the radius within the circumbinary disc at which the breaks appear, and furthermore the radius at which the disc breaks can change as a function of time in each case. We discuss the implications of our results for interpreting observations and suggest some considerations for modelling misaligned discs in the future.
△ Less
Submitted 11 August, 2023; v1 submitted 20 June, 2023;
originally announced June 2023.
-
Characteristics of small protoplanetary disc warps in kinematic observations
Authors:
Alison K. Young,
Richard Alexander,
Giovanni Rosotti,
Christophe Pinte
Abstract:
Many circumstellar discs appear to have misaligned central regions that give rise to shadows seen in scattered light observations. Small warps ($<20^\circ$ misalignment) are probably more common but are also more difficult to detect than the large misalignments studied previously. We present the characteristics of CO emission that may be used to identify a small disc warp, found from synthetic…
▽ More
Many circumstellar discs appear to have misaligned central regions that give rise to shadows seen in scattered light observations. Small warps ($<20^\circ$ misalignment) are probably more common but are also more difficult to detect than the large misalignments studied previously. We present the characteristics of CO emission that may be used to identify a small disc warp, found from synthetic $^{13}$CO maps of a model misaligned circumbinary disc. The spectra are not symmetrical, so fitting a Keplerian model is not appropriate and can hide a warp or lead to spurious features such as spirals appearing in the residuals. We quantify the observed warp structure by fitting sinusoids to concentric annuli of the disc. From this we can trace the radial variation of the peak velocity and of the azimuth of the peak velocity, i.e., the twist. At near face-on inclinations, these radial profiles reveal the warp structure. The twist remains detectable at moderate inclinations (${i_{\rm outer~disc}\lesssim 35^{\circ}}$) in the absence of radial flows but the measured inclination must be accurate to $\lesssim 5^{\circ}$ to allow detection of the radial variation. The observed twist does not provide a direct measure of the warp structure because of its dependence on optical depth. The warp causes broad asymmetries in the channel maps that span several channels and that are distinct from localised features caused by embedded planets and gravitational instability. We suspect that kinematic evidence of warps may have been missed and we suggest a few examples where the data may be revisited.
△ Less
Submitted 23 March, 2022;
originally announced March 2022.
-
Chemical signatures of a warped protoplanetary disc
Authors:
Alison K. Young,
Richard Alexander,
Catherine Walsh,
Rebecca Nealon,
Alice Booth,
Christophe Pinte
Abstract:
Circumstellar discs may become warped or broken into distinct planes if there is a stellar or planetary companion with an orbit that is misaligned with respect to the disc. There is mounting observational evidence for protoplanetary discs with misaligned inner discs and warps that may be caused by such interactions with a previously undetected companion, giving us a tantalising indication of possi…
▽ More
Circumstellar discs may become warped or broken into distinct planes if there is a stellar or planetary companion with an orbit that is misaligned with respect to the disc. There is mounting observational evidence for protoplanetary discs with misaligned inner discs and warps that may be caused by such interactions with a previously undetected companion, giving us a tantalising indication of possible planets forming there. Hydrodynamical and radiative transfer models indicate that the temperature varies azimuthally in warped discs due to the variable angle at which the disc surface faces the star and this impacts the disc chemistry. We perform chemical modelling based on a hydrodynamical model of a protoplanetary disc with an embedded planet orbiting at a 12$^{\circ}$ inclination to the disc. Even for this small misalignment, abundances of species including CO and HCO$^+$ vary azimuthally and this results in detectable azimuthal variations in submillimetre line emission. Azimuthal variations in line emission may therefore indicate the presence of an unseen embedded companion. Nonaxisymmetric chemical abundances should be considered when interpreting molecular line maps of warped or shadowed protoplanetary discs.
△ Less
Submitted 4 June, 2021;
originally announced June 2021.
-
A triple star system with a misaligned and warped circumstellar disk shaped by disk tearing
Authors:
Stefan Kraus,
Alexander Kreplin,
Alison K. Young,
Matthew R. Bate,
John D. Monnier,
Tim J. Harries,
Henning Avenhaus,
Jacques Kluska,
Anna S. E. Laws,
Evan A. Rich,
Matthew Willson,
Alicia N. Aarnio,
Fred C. Adams,
Sean M. Andrews,
Narsireddy Anugu,
Jaehan Bae,
Theo ten Brummelaar,
Nuria Calvet,
Michel Curé,
Claire L. Davies,
Jacob Ennis,
Catherine Espaillat,
Tyler Gardner,
Lee Hartmann,
Sasha Hinkley
, et al. (7 additional authors not shown)
Abstract:
Young stars are surrounded by a circumstellar disk of gas and dust, within which planet formation can occur. Gravitational forces in multiple star systems can disrupt the disk. Theoretical models predict that if the disk is misaligned with the orbital plane of the stars, the disk should warp and break into precessing rings, a phenomenon known as disk tearing. We present observations of the triple…
▽ More
Young stars are surrounded by a circumstellar disk of gas and dust, within which planet formation can occur. Gravitational forces in multiple star systems can disrupt the disk. Theoretical models predict that if the disk is misaligned with the orbital plane of the stars, the disk should warp and break into precessing rings, a phenomenon known as disk tearing. We present observations of the triple star system GWOrionis, finding evidence for disk tearing. Our images show an eccentric ring that is misaligned with the orbital planes and the outer disk. The ring casts shadows on a strongly warped intermediate region of the disk. If planets can form within the warped disk, disk tearing could provide a mechanism for forming wide-separation planets on oblique orbits.
△ Less
Submitted 4 September, 2020; v1 submitted 2 April, 2020;
originally announced April 2020.
-
Synthetic molecular line observations of the first hydrostatic core from chemical calculations
Authors:
Alison K. Young,
Matthew R. Bate,
Tim J. Harries,
David M. Acreman
Abstract:
The first stable object to develop in the low-mass star formation process has long been predicted to be the first hydrostatic core (FHSC). Despite much effort, it has still yet to be definitively observed in nature. More specific observational signatures are required to enable observers to distinguish the FHSC from young, faint, but more evolved protostars. Here we present synthetic spectral line…
▽ More
The first stable object to develop in the low-mass star formation process has long been predicted to be the first hydrostatic core (FHSC). Despite much effort, it has still yet to be definitively observed in nature. More specific observational signatures are required to enable observers to distinguish the FHSC from young, faint, but more evolved protostars. Here we present synthetic spectral line observations for CO, SO, CS and HCO$^+$ that were calculated from radiation (magneto)hydrodynamical models, chemical modelling and Monte Carlo radiative transfer. HCO$^+$ $(1-0)$ and SO $(8_7 - 7_6)$ spectra of the FHSC show variations for observations at a low inclination which may allow a candidate FHSC to be distinguished from a more evolved object. We find that the FHSC outflow is unlikely to be detectable with ALMA, which would discount the observed sources with slow outflows that are currently identified as candidate FHSCs. We compare the results of simulated ALMA observations with observed candidate FHSCs and recommend Oph A SM1N and N6-mm as the most promising candidates to follow up.
△ Less
Submitted 4 June, 2019;
originally announced June 2019.
-
Dust-trapping vortices and a potentially planet-triggered spiral wake in the pre-transitional disk of V1247 Orionis
Authors:
Stefan Kraus,
Alexander Kreplin,
Misato Fukagawa,
Takayuki Muto,
Michael L. Sitko,
Alison K. Young,
Matthew R. Bate,
Timothy Harries,
John D. Monnier,
Matthew Willson,
John Wisniewski
Abstract:
The radial drift problem constitutes one of the most fundamental problems in planet formation theory, as it predicts particles to drift into the star before they are able to grow to planetesimal size. Dust-trapping vortices have been proposed as a possible solution to this problem, as they might be able to trap particles over millions of years, allowing them to grow beyond the radial drift barrier…
▽ More
The radial drift problem constitutes one of the most fundamental problems in planet formation theory, as it predicts particles to drift into the star before they are able to grow to planetesimal size. Dust-trapping vortices have been proposed as a possible solution to this problem, as they might be able to trap particles over millions of years, allowing them to grow beyond the radial drift barrier. Here, we present ALMA 0.04"-resolution imaging of the pre-transitional disk of V1247 Orionis that reveals an asymmetric ring as well as a sharply-confined crescent structure, resembling morphologies seen in theoretical models of vortex formation. The asymmetric ring (at 0.17"=54 au separation from the star) and the crescent (at 0.38"=120 au) seem smoothly connected through a one-armed spiral arm structure that has been found previously in scattered light. We propose a physical scenario with a planet orbiting at $\sim0.3$"$\approx$100 au, where the one-armed spiral arm detected in polarised light traces the accretion stream feeding the protoplanet. The dynamical influence of the planet clears the gap between the ring and the crescent and triggers two vortices that trap mm-sized particles, namely the crescent and the bright asymmetry seen in the ring. We conducted dedicated hydrodynamics simulations of a disk with an embedded planet, which results in similar spiral-arm morphologies as seen in our scattered light images. At the position of the spiral wake and the crescent we also observe $^{12}$CO (3-2) and H$^{12}$CO$^{+}$ (4-3) excess line emission, likely tracing the increased scale-height in these disk regions.
△ Less
Submitted 13 October, 2017;
originally announced October 2017.
-
What can the SEDs of first hydrostatic core candidates reveal about their nature?
Authors:
Alison K. Young,
Matthew R. Bate,
Chris F. Mowat,
Jennifer Hatchell,
Tim J. Harries
Abstract:
The first hydrostatic core (FHSC) is the first stable object to form in simulations of star formation. This stage has yet to be observed definitively, although several candidate FHSCs have been reported. We have produced synthetic spectral energy distributions (SEDs) from 3D hydrodynamical simulations of pre-stellar cores undergoing gravitational collapse for a variety of initial conditions. Varia…
▽ More
The first hydrostatic core (FHSC) is the first stable object to form in simulations of star formation. This stage has yet to be observed definitively, although several candidate FHSCs have been reported. We have produced synthetic spectral energy distributions (SEDs) from 3D hydrodynamical simulations of pre-stellar cores undergoing gravitational collapse for a variety of initial conditions. Variations in the initial rotation rate, radius and mass lead to differences in the location of the SED peak and far-infrared flux. Secondly, we attempt to fit the SEDs of five FHSC candidates from the literature and five newly identified FHSC candidates located in the Serpens South molecular cloud with simulated SEDs. The most promising FHSC candidates are fitted by a limited number of model SEDs with consistent properties, which suggests the SED can be useful for placing constraints on the age and rotation rate of the source. The sources we consider most likely to be in FHSC phase are B1-bN, CB17-MMS, Aqu-MM1 and Serpens South candidate K242. We were unable to fit SerpS-MM22, Per-Bolo 58 and Chamaeleon-MMS1 with reasonable parameters, which indicates that they are likely to be more evolved.
△ Less
Submitted 12 October, 2017;
originally announced October 2017.
-
The plasmoid instability during asymmetric inflow magnetic reconnection
Authors:
Nicholas A. Murphy,
Aleida K. Young,
Chengcai Shen,
Jun Lin,
Lei Ni
Abstract:
Theoretical studies of the plasmoid instability generally assume that the reconnecting magnetic fields are symmetric. We relax this assumption by performing two-dimensional resistive magnetohydrodynamic simulations of the plasmoid instability during asymmetric inflow magnetic reconnection. Magnetic asymmetry modifies the onset, scaling, and dynamics of this instability. Magnetic islands develop pr…
▽ More
Theoretical studies of the plasmoid instability generally assume that the reconnecting magnetic fields are symmetric. We relax this assumption by performing two-dimensional resistive magnetohydrodynamic simulations of the plasmoid instability during asymmetric inflow magnetic reconnection. Magnetic asymmetry modifies the onset, scaling, and dynamics of this instability. Magnetic islands develop preferentially into the weak magnetic field upstream region. Outflow jets from individual X-points impact plasmoids obliquely rather than directly as in the symmetric case. Consequently, deposition of momentum by the outflow jets into the plasmoids is less efficient, the plasmoids develop net vorticity, and shear flow slows down secondary merging between islands. Secondary merging events have asymmetry along both the inflow and outflow directions. Downstream plasma is more turbulent in cases with magnetic asymmetry because islands are able to roll around each other after exiting the current sheet. As in the symmetric case, plasmoid formation facilitates faster reconnection for at least small and moderate magnetic asymmetries. However, when the upstream magnetic field strengths differ by a factor of four, the reconnection rate plateaus at a lower value than expected from scaling the symmetric results. We perform a parameter study to investigate the onset of the plasmoid instability as a function of magnetic asymmetry and domain size. There exist domain sizes for which symmetric simulations are stable but asymmetric simulations are unstable, suggesting that moderate magnetic asymmetry is somewhat destabilizing. We discuss the implications for plasmoid and flux rope formation in solar eruptions, laboratory reconnection experiments, and space plasmas. The differences between symmetric and asymmetric simulations provide some hints regarding the nature of the three-dimensional plasmoid instability.
△ Less
Submitted 15 May, 2013;
originally announced May 2013.