Disk Evolution Research Articles

Evolution of the Magnetic Field in High- and Low-β Disks with Initially Toroidal Fields

Abstract We present the results from a pair of high-resolution, long-timescale (∼105 GM/c 3), global, three-dimensional magnetohydrodynamical accretion disk simulations with differing initial magnetic plasma β in order to study the effects of the initial toroidal field strength on the production of a large-scale poloidal field. We initialize our disks in approximate equilibrium with purely toroidal magnetic fields of strength β 0 = 5 and 200. We also perform a limited resolution study. We find that simulations of differing field strengths diverge early in their evolution and remain distinct over the time studied, indicating that the initial magnetic conditions leave a persistent imprint in our simulations. Neither simulation enters the magnetically arrested disk regime. Both simulations are able to produce poloidal fields from initially toroidal fields, with the β 0 = 5 simulation evolving clear signs of a large-scale poloidal field. We make a cautionary note that computational artifacts in the form of large-scale vortices may be introduced in the combination of initially weak field and disk-internal mesh refinement boundaries, as evidenced by the production of an m = 1 mode overdensity in the weak field simulation. Our results demonstrate that the initial toroidal field strength plays a vital role in the simulated disk evolution for the models studied.

Evolution of eccentric stellar discs around supermassive black holes: the complex disc disruption dynamics and the milliparsec stars

ABSTRACT We study the 10 Myr evolution of parsec-scale stellar discs with initial masses of Mdisc = 1.0–$7.5 \times 10^4\, \mathrm{M}_\odot$ and eccentricities einit = 0.1–0.9 around supermassive black holes (SMBHs). Our disc models are embedded in a spherical background potential and have top-heavy single and binary star initial mass functions (IMF) with slopes of 0.25–1.7. The systems are evolved with the N-body code BIFROST, including post-Newtonian (PN) equations of motion and simplified stellar evolution. All discs are unstable and evolve on Myr time-scales towards similar eccentricity distributions peaking at e⋆ ∼ 0.3–0.4. Models with high einit also develop a very eccentric (e⋆ ≳ 0.9) stellar population. For higher disc masses Mdisc ≳ 3 × 104 M⊙, the disc disruption dynamics is more complex than the standard secular eccentric disc instability with opposite precession directions at different disc radii – a precession direction instability. We present an analytical model describing this behaviour. A milliparsec population of N ∼ 10–100 stars forms around the SMBH in all models. For low einit, stars migrate inward while for einit ≳ 0.6 stars are captured by the Hills mechanism. Without PN, after 6 Myr, the captured stars have a sub-thermal eccentricity distribution. We show that including PN effects prevents this thermalization by suppressing resonant relaxation effects and cannot be ignored. The number of tidally disrupted stars is similar or larger than the number of milliparsec stars. None of the simulated models can simultaneously reproduce the kinematic and stellar population properties of the Milky Way centre clockwise disc and the S-cluster.

Role of redox conditions and thermal metamorphism in the preservation of Cr isotopic anomalies in components of non-carbonaceous chondrites

Mass-independent isotopic anomalies in meteorites are probes to the dynamics and evolution of the protoplanetary disc and the reservoirs from which planetary bodies accreted. Variable Cr nucleosynthetic compositions between meteorite groups are observed but the cause is not well understood. Investigations on presolar carriers in unequilibrated chondrites are thus required to establish these relationships. For that purpose, we analysed the Cr isotope compositions of leachates (stepwise dissolution solutions) of one ordinary chondrite and two enstatite chondrites of EL and EH subgroups previously measured for Nd and Sm isotopes. The leachates of the ordinary chondrite display large variations of their Cr isotope compositions, whereas the leachates of the enstatite chondrites show very limited variations of their nucleosynthetic compositions, confirming other studies. The Cr isotope compositions of leachates of unequilibrated chondrites are significantly modified by thermal processing and aqueous alteration under various conditions. The leachates of the enstatite chondrites are relatively homogeneous for Cr isotope composition, suggesting that presolar carriers of Cr are not stable under very reduced conditions and are likely oxidised phases. After careful examination of the role of parent body processing, we conclude that there is one predominant presolar carrier enriched in 54Cr that produced most of the nucleosynthetic variations observed in leachates. This carrier is likely depleted in the carbonaceous reservoir relative to the non-carbonaceous reservoir. The drastically different isotope compositions of Cr and heavy refractory lithophile elements such as Nd of leachates of chondrites reflect relative preservation biases of their respective presolar carriers, namely a still poorly characterised presolar oxide and silicon carbide (SiC), respectively.

Protostellar discs fed by dense collapsing gravomagneto sheetlets

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.

Insights into the first and second hydrostatic core stages from numerical simulations

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.

Testing external photoevaporation in the σ-Orionis cluster with spectroscopy and disk mass measurements (Corrigendum)

{The evolution of protoplanetary disks is regulated by an interplay of several processes, either internal to the system or related to the environment. As most of the stars and planets, including our own Solar System, have formed in massive stellar clusters that contain OB-type stars, studying the effects of UV radiation on disk evolution is of paramount importance.} % aims heading (mandatory) {Here we test the impact of external photoevaporation on the evolution of disks in the mid-age ($sim$3--5 Myr) $sigma$-Orionis cluster by conducting the first combined large-scale UV to IR spectroscopic and mm-continuum survey of this region.} % methods heading (mandatory) {We study a sample of 50 targets located at increasing distances from the central, massive OB system $sigma$-Ori. We combine new spectra obtained with VLT/X-Shooter, used to measure mass accretion rates and stellar masses, with new and previously published ALMA measurements of disk dust and gas fluxes and masses.} % results heading (mandatory) {We confirm the previously found decrease of $M_{rm dust}$ in the inner $sim$0.5 pc of the cluster, and discover an additional dependence with the stellar host mass. The disks around the more massive stars ($ge$ 0.4msun) and located in the inner part ($<$ 0.5 pc) of the cluster have $M_{rm dust}$ about an order of magnitude lower than the more distant ones, in contrast to the rather flat $M_{rm dust}$ distribution as a function of the projected separation from $sigma$-Ori, for the disks around the lowest-mass stars in the sample. About half of the sample is located in the region of the $dot{M}_{rm acc}$ vs $M_{rm disk}$ expected by models of external photoevaporation, i.e. showing shorter disk lifetimes than expected for their ages. The shorter disk lifetimes is observed for all targets with projected separation from $sigma$-Ori $<$ 0.5 pc, proving that the presence of a massive stellar system affects disk evolution. } % conclusions heading (optional), leave it empty if necessary {External photoevaporation is a viable mechanism to explain the observed shorter disk lifetimes and lower $M_{rm dust}$ in the inner $sim$0.5 pc of the $sigma$-Orionis cluster, where the effects of this process are more pronounced. Follow-up observations of the low stellar mass targets are crucial to constrain disk dispersion time scales in the cluster and to confirm the dependence of the external photoevaporation process with stellar host mass. This work confirms that the effects of external photoevaporation are significant down to at least impinging radiation as low as $sim 10^{4}$ G$_0$. }

A Primordial Origin for the Gas-rich Debris Disks around Intermediate-mass Stars

While most debris disks consist of dust with little or no gas, a fraction have significant amounts of gas detected via emission lines of CO, ionized carbon, and/or atomic oxygen. Almost all such gaseous debris disks known are around A-type stars with ages up to 50 Myr. We show, using semianalytic disk evolution modeling, that this can be understood if the gaseous debris disks are remnant protoplanetary disks that have become depleted of small grains compared to the interstellar medium. Photoelectric heating by the A stars’ far-UV (FUV) radiation is then inefficient, while the stars’ extreme-UV (EUV) and X-ray emissions are weak owing to a lack of surface convective zones capable of driving magnetic activity. In this picture, it is relatively difficult for stars outside the range of spectral types from A through early B to have such long-lived gas disks. Less-massive stars have stronger magnetic activity in the chromosphere, transition region, and corona with resulting EUV and X-ray emission, while more-massive stars have photospheres hot enough to produce strong EUV radiation. In both cases, primordial disk gas is likely to photoevaporate well before 50 Myr. These results come from 0D disk evolution models where we incorporate internal accretion stresses, MHD winds, and photoevaporation by EUV and X-ray photons with luminosities that are functions of the stellar mass and age. A key issue this work leaves open is how some disks become depleted in small dust so that FUV photoevaporation slows. Candidates include the grains’ growth, settling, radial drift, radiation force, and incorporation into planetary systems.

Photoevaporation obfuscates the distinction between wind and viscous angular momentum transport in protoplanetary discs

ABSTRACT How protoplanetary discs evolve remains an unanswered question. Competing theories of viscosity and magnetohydrodynamic disc winds have been put forward as the drivers of angular momentum transport in protoplanetary discs. These two models predict distinct differences in the disc mass, radius, and accretion rates over time, that could be used to distinguish them. However that expectation is built on models that do not include another important process – photoevaporation, both internally by the host star and externally by neighbouring stars. In this work we produce numerical models of protoplanetary discs including viscosity, magnetohydrodynamic disc winds, and internal and external photoevaporation. We find that even weak levels of external photoevaporation can significantly affect the evolution of protoplanetary discs, influencing the observable features such as disc radii, that might otherwise distinguish between viscous and wind driven discs. Including internal photoevaporation further suppresses differences in evolution between viscous and wind driven discs. This makes it much more difficult than previously anticipated, to use observations of nearby star forming regions to determine whether discs are viscous or wind driven. Interestingly we find that evolved protoplanetary discs in intermediate FUV environments may be the best cases for differentiating whether they evolve through viscosity or magnetohydrodynamic disc winds. Ultimately this work demonstrates the importance of understanding what are the key evolutionary processes and including as many of those as possible when exploring the evolution of protoplanetary discs.

The effect of Jupiter on the CAI storage problem

Context. Meteorites preserve an imprint of conditions in the early Solar System. By studying the distribution of calcium-aluminium-rich inclusions (CAIs) that are embedded within meteorites, we can learn about the dynamical history of the protoplanetary disk from which our Solar System formed. A long-standing problem concerning CAIs is the CAI storage problem. CAIs are thought to have formed at high temperatures near the Sun, but they are primarily found in carbonaceous chondrites, which formed much further out, beyond the orbit of Jupiter. Additionally, radial drift of CAI particles should have removed them from the solar protoplanetary disk several million years before the parent bodies of meteorites in which they are encountered would have accreted. Aims. We revisit a previously suggested solution to the CAI storage problem by Desch, Kalyaan, and Alexander which proposed that CAIs were mixed radially outward through the disk and subsequently got trapped in a pressure maximum created by Jupiter’s growing core opening a planet gap. Our aim is to investigate whether their solution still works when we take into account the infall phase during which the disk builds up from the collapse of a molecular cloud core. Methods. We build a 1D numerical code in Python using the DISKLAB package to simulate the evolution of the solar protoplanetary disk, starting with a collapsing molecular cloud. CAIs are created in the model by thermal processing of solar nebula composition dust, and subsequently transported through the disk by turbulent diffusion, radial drift and advection by the gas. Results. We find that outward transport of CAIs during the infall phase is very efficient, possibly mixing them all the way into the far outer disk. Subsequent inward radial drift collects CAIs in the pressure maximum beyond Jupiter’s orbit while draining the inner disk, roughly reproducing parts of the result by Desch et al. By introducing CAI formation so early, abundances out to 100 AU remain significant, possibly not consistent with some meteoritic data. It is possible to create a disk that does not expand as far out and also does not push CAIs as far out by using a very slowly rotating cloud.

Millimeter emission in photoevaporating disks is determined by early substructures

Context. Photoevaporation and dust-trapping are individually considered to be important mechanisms in the evolution and morphology of protoplanetary disks. However, it is not yet clear what kind of observational features are expected when both processes operate simultaneously. Aims. We studied how the presence (or absence) of early substructures, such as the gaps caused by planets, affects the evolution of the dust distribution and flux in the millimeter continuum of disks that are undergoing photoevaporative dispersal. We also tested if the predicted properties resemble those observed in the population of transition disks. Methods. We used the numerical code Dustpy to simulate disk evolution considering gas accretion, dust growth, dust-trapping at substructures, and mass loss due to X-ray and EUV (XEUV) photoevaporation and dust entrainment. Then, we compared how the dust mass and millimeter flux evolve for different disk models. Results. We find that, during photoevaporative dispersal, disks with primordial substructures retain more dust and are brighter in the millimeter continuum than disks without early substructures, regardless of the photoevaporative cavity size. Once the photoevaporative cavity opens, the estimated fluxes for the disk models that are initially structured are comparable to those found in the bright transition disk population (Fmm > 30 mJy), while the disk models that are initially smooth have fluxes comparable to the transition disks from the faint population (Fmm < 30 mJy), suggesting a link between each model and population. Conclusions. Our models indicate that the efficiency of the dust trapping determines the millimeter flux of the disk, while the gas loss due to photoevaporation controls the formation and expansion of a cavity, decoupling the mechanisms responsible for each feature. In consequence, even a planet with a mass comparable to Saturn could trap enough dust to reproduce the millimeter emission of a bright transition disk, while its cavity size is independently driven by photoevaporative dispersal.

Magnetic support, wind-driven accretion, coronal heating, and fast outflows in a thin magnetically arrested disc

ABSTRACT Accretion discs properties should deviate from standard theory when magnetic pressure exceeds the thermal pressure. To quantify these deviations, we present a systematic study of the dynamical properties of magnetically arrested discs (MADs), the most magnetized type of accretion disc. Using an artificial cooling function to regulate the gas temperature, we study MADs of three different thermal thicknesses, hth/r = 0.3, 0.1, and 0.03. We find that the radial structure of the disc is never mostly supported by the magnetic field. In fact, thin MADs are very near Keplerian. However, as discs gets colder, they become more magnetized and the largest deviations from standard theory appear in our thinnest disc with hth/r = 0.03. In this case, the disc is much more extended vertically and much less dense than in standard theory because of vertical support from the turbulent magnetic pressure and wind-driven angular momentum transport that enhances the inflow speed. The thin disc also dissipates a lot of thermal energy outside of z/r = ±0.03 and a significant fraction of this dissipation happens in mildly relativistic winds. The enhanced dissipation in low-density regions could possibly feed coronae in X-ray binaries (XRBs) and active galactic nuclei (AGNs). Wind-driven accretion will also impact the dynamical evolution of accretion discs and could provide a mechanism to explain the rapid evolution of changing-look AGN and the secular evolution of XRBs. Finally, our MAD winds have terminal velocities and mass-loss rates in good agreement with the properties of ultrafast outflows observed in AGN.

Photoevaporation and Close Encounters: How the Environment around Cygnus OB2 Affects the Evolution of Protoplanetary Disks

In our Galaxy, star formation occurs in a variety of environments, with a large fraction of stars formed in clusters hosting massive stars. OB stars have an important feedback on the evolution of protoplanetary disks orbiting around nearby young stars and likely on the process of planet formation occurring in them. The nearby massive association Cyg OB2 is an outstanding laboratory to study this feedback. It is the closest massive association to our Sun and hosts hundreds of massive stars and thousands of low-mass members, both with and without disks. In this paper, we analyze the spatial variation of the disk fraction (i.e., the fraction of cluster members bearing a disk) in Cyg OB2 and study its correlation with the local values of far-ultraviolet (FUV) and extreme-ultraviolet (EUV) radiation fields and the local stellar surface density. We present definitive evidence that disks are more rapidly dissipated in the regions of the association characterized by intense local UV fields and large stellar density. In particular, the FUV radiation dominates disk dissipation timescales in the proximity (i.e., within 0.5 pc) of the O stars. In the rest of the association, EUV photons potentially induce a significant mass loss from the irradiated disks across the entire association, but the efficiency of this process is reduced at increasing distances from the massive stars owing to absorption by the intervening intracluster material. We find that disk dissipation due to close stellar encounters is negligible in Cyg OB2 and likely to have affected 1% or fewer of the stellar population. Disk dissipation is instead dominated by photoevaporation. We also compare our results to what has been found in other young clusters with different massive populations, concluding that massive associations like Cyg OB2 are potentially hostile to protoplanetary disks but that the environments where disks can safely evolve in planetary systems are likely quite common in our Galaxy.

Twenty-five Years of Accretion onto the Classical T Tauri Star TW Hya

Accretion plays a central role in the physics that governs the evolution and dispersal of protoplanetary disks. The primary goal of this paper is to analyze the stability over time of the mass accretion rate onto TW Hya, the nearest accreting solar-mass young star. We measure veiling across the optical spectrum in 1169 archival high-resolution spectra of TW Hya, obtained from 1998–2022. The veiling is then converted to accretion rate using 26 flux-calibrated spectra that cover the Balmer jump. The accretion rate measured from the excess continuum has an average of 2.51 × 10−9 M ⊙ yr−1 and a Gaussian distribution with an FWHM of 0.22 dex. This accretion rate may be underestimated by a factor of up to 1.5 because of uncertainty in the bolometric correction and another factor of 1.7 because of excluding the fraction of accretion energy that escapes in lines, especially Lyα. The accretion luminosities are well correlated with He line luminosities but poorly correlated with Hα and Hβ luminosity. The accretion rate is always flickering over hours but on longer timescales has been stable over 25 years. This level of variability is consistent with previous measurements for most, but not all, accreting young stars.

Three-dimensional magnetic field imaging of protoplanetary disks using Zeeman broadening and linear polarization observations

Context. Magnetic fields are predicted to have a crucial impact on the structure, evolution, and chemistry of protoplanetary disks. However, a direct detection of the magnetic field toward these objects has yet to be achieved. Aims. In order to characterize the magnetic fields of protoplanetary disks, we investigate the impact of the Zeeman effect on the (polarized) radiative transfer of emission from paramagnetic molecules excited in protoplanetary disks. Methods. While the effects of the Zeeman effect are commonly studied in the circular polarization of spectral lines, we also performed a comprehensive modeling of the Zeeman-induced broadening of spectral lines and their linear polarization. We developed simplified radiative transfer models adapted to protoplanetary disks, which we compare to full three-dimensional polarized radiative transfer simulations. Results. We find that the radiative transfer of circular polarization is heavily affected by the expected polarity change of the magnetic field between opposite sides of the disk. In contrast, Zeeman broadening and linear polarization are relatively unaffected by this sign change due to their quadratic dependence on the magnetic field. We can match our simplified radiative transfer models to full polarization modeling with high fidelity, which in turn allows us to prescribe straightforward methods to extract magnetic field information from Zeeman broadening and linear polarization observations. Conclusions. We find that Zeeman broadening and linear polarization observations are highly advantageous methods to characterize protoplanetary disk magnetic fields as they are both sensitive probes of the magnetic field and are marginally affected by any sign change of the disk magnetic field. Applying our results to existing circular polarization observations of protoplanetary disk specral lines suggests that the current upper limits on the toroidal magnetic field strengths have to be raised.

Long-term Evolution of Massive-star Post-common-envelope Circumbinary Disks and the Environments of Fast Luminous Transients

If the envelope of a massive star is not entirely removed during common envelope (CE) interaction with an orbiting compact (e.g., black hole (BH) or neutron star (NS)) companion, the residual bound material eventually cools, forming a centrifugally supported disk around the binary containing the stripped He core. We present a time-dependent height-integrated model for the long-term evolution of post-CE circumbinary disks (CBDs), accounting for mass and angular momentum exchange with the binary, irradiation heating by the He core, and photoevaporation wind mass loss. A large fraction of the CBD’s mass is accreted prior to its outwards viscous spreading and wind dispersal on a timescale of ∼104–105 yr, driving significant orbital migration, even for disks containing ∼10% of the original envelope mass. Insofar that the CBD lifetime is comparable to the thermal (and, potentially, nuclear) timescale of the He core, over which a second mass-transfer episode onto the companion can occur, the presence of the CBD could impact the stability of this key phase. Disruption of the core by the BH/NS would result in a jetted energetic explosion into the dense gaseous CBD (≲1015 cm) and its wind (≳1016 cm), consistent with the environments of luminous fast blue optical transients like AT2018cow. Evolved He cores that undergo core collapse still embedded in their CBD could generate Type Ibn/Icn supernovae. Thousands of dusty wind-shrouded massive-star CBDs may be detectable as extragalactic luminous infrared sources with the Roman Space Telescope; synchrotron radio nebulae powered by the CBD-fed BH/NS may accompany these systems.

High-resolution Chemical Abundances of the Nyx Stream

Nyx is a nearby, prograde, and high-eccentricity stellar stream physically contained in the thick disk, but its origin is unknown. Nyx could be the remnant of a disrupted dwarf galaxy, in which case the associated dark matter substructure could affect terrestrial dark matter direct-detection experiments. Alternatively, Nyx could be a signature of the Milky Way’s disk formation and evolution. To determine the origin of Nyx, we obtained high-resolution spectroscopy of 34 Nyx stars using Keck/HIRES and Magellan/MIKE. A differential chemical abundance analysis shows that most Nyx stars reside in a metal-rich ([Fe/H] > −1) high-α component that is chemically indistinguishable from the thick disk. This rules out the originally suggested scenario that Nyx is the remnant of a single massive dwarf galaxy merger. However, we also identify 5 substantially more metal-poor stars ([Fe/H] ∼ −2.0) whose chemical abundances are similar to those of the metal-weak thick disk. It remains unclear how stars that are chemically identical to the thick disk can be on such prograde, high-eccentricity orbits. We suggest two most likely scenarios: that Nyx is the result of an early minor dwarf galaxy merger, or that it is a record of the early spin-up of the Milky Way disk—although neither perfectly reproduces the chemodynamic observations. The most likely formation scenarios suggest that future spectroscopic surveys should find Nyx-like structures outside of the solar neighborhood.

Radio-continuum decrements associated to shadowing from the central warp in transition disc DoAr 44

ABSTRACT Warps have often been used to explain disc properties, but well-characterized examples are important due to their role in disc evolution. Scattered light images of discs with central gaps have revealed sharp warps, such that the outer rings are shadowed by tilted inner discs. The near-infrared intensity drops along the ring around T-Tauri star DoAr 44 have been interpreted in terms of a central warp. We report new ALMA observations of DoAr 44 in the continuum at 230 and 350 GHz (at ∼10 au), along with a new epoch of Spectro Polarimetric High-contrast Exoplanet REsearch (SPHERE)/Infrared Dual-band Imager and Spectrograph differential polarized imaging taken during excellent weather conditions. The ALMA observations resolve the ring and confirm the decrements proposed from deconvolution of coarse 336 GHz data. The scattered light image constrains the dips, which correspond to a misaligned inner disc with a relative inclination ξ = 21.4 $^{+6.7}_{-8.3}$ deg. The SPHERE intensity profile shows a morphological change compared to a previous epoch that may be interpreted as a variable orientation of the inner disc, from ξ ∼ 30 deg to ξ ∼ 20 deg. The intensity dips probably correspond to temperature decrements, as their mm-spectral index, $\alpha ^{230 \textrm {GHz}}_{350 \textrm {GHz}} \sim$ 2.0 ± 0.1, is indicative of optically thick emission. The azimuth of the two temperature decrements are leading clockwise relative to the infrared-dips, by η = 14.95 deg and η = 7.92 deg. For a retrograde disc, such shifts are expected from a thermal lag and imply gas surface densities of Σg = 117 ± 10 g cm−2 and Σg = 48 ± 10 g cm−2. A lopsided disc, with contrast ratio fr = 2.4 ± 0.5, is also consistent with the large continuum crescent.

Azimuthal temperature variations in ISO-Oph 2 from multifrequency ALMA observations

ABSTRACT Environmental effects, such as stellar fly-bys and external irradiation, are thought to affect the evolution of protoplanetary discs in clustered star formation. Previous Atacama Large Millimetre/submillimetre Array (ALMA) images at 225 GHz of the ISO-Oph 2 binary revealed a peculiar morphology in the disc of the primary, perhaps due to a possible fly-by with the secondary. Here, we report on new ALMA continuum observations of this system at 97.5, 145, and 405 GHz, which reveal strong morphological variations. Multifrequency positional alignment allows us to interpret these spectral variations in terms of underlying physical conditions. ISO-Oph 2A is remarkably offset from the centroid of its ring, at all frequencies, and the disc is lopsided, pointing at gravitational interactions. However, the dust temperature also varies in azimuth, with two peaks whose direction connects with HD 147889, the earliest-type star in the Ophiuchus complex, suggesting that it is the dominant heat source. The stellar environment of ISO-Oph 2 appears to drive both its density structure and its thermal balance.

The dynamical evolution of protoplanetary discs and planets in dense star clusters

ABSTRACT Most stars are born in dense stellar environments, where the formation and early evolution of planetary systems may be significantly perturbed by encounters with neighbouring stars. To investigate the fate of circumstellar gas discs and planets around young stars dense stellar environments, we numerically evolve star–disc–planet systems. We use the N-body codes NBODY6++ GPU and SnIPES for the dynamical evolution of the stellar population, and the SPH-based code GaSPH for the dynamical evolution of protoplanetary discs. The secular evolution of a planetary system in a cluster differs from that of a field star. Most stellar encounters are tidal, adiabatic, and nearly-parabolic. The parameters that characterize the impact of an encounter include the orientation of the protoplanetary disc and planet relative to the orbit of the encountering star, and the orbital phase, and the semimajor axis of the planet. We investigate this dependence for close encounters (rp/a ≤ 100, where rp is the periastron distance of the encountering star and a is the semimajor axis of the planet). We also investigate distant perturbers (rp/a ≫ 100), which have a moderate effect on the dynamical evolution of the planet and the protoplanetary disc. We find that the evolution of protoplanetary discs in star clusters differs significantly from that of isolated systems. When interpreting the outcome of the planet formation process, it is thus important to consider their birth environments.

Decoding the age–chemical structure of the Milky Way disc: an application of copulas and elicitable maps

ABSTRACT In the Milky Way, the distribution of stars in the [α/Fe] versus [Fe/H] and [Fe/H] versus age planes holds essential information about the history of star formation, accretion, and dynamical evolution of the Galactic disc. We investigate these planes by applying novel statistical methods called copulas and elicitable maps to the ages and abundances of red giants in the Apache Point Observatory Galactic Evolution Experiment survey. We find that the high- and low-α disc stars have a clean separation in copula space and use this to provide an automated separation of the α sequences using a purely statistical approach. This separation reveals that the high-α disc ends at the same [α/Fe] and age at high [Fe/H] as the low-[Fe/H] start of the low-α disc, thus supporting a sequential formation scenario for the high- and low-α discs. We then combine copulas with elicitable maps to precisely obtain the correlation between stellar age τ and metallicity [Fe/H] conditional on Galactocentric radius R and height z in the range 0 < R < 20 kpc and |z| < 2 kpc. The resulting trends in the age–metallicity correlation with radius, height, and [α/Fe] demonstrate a ≈0 correlation wherever kinematically cold orbits dominate, while the naively expected negative correlation is present where kinematically hot orbits dominate. This is consistent with the effects of spiral-driven radial migration, which must be strong enough to completely flatten the age–metallicity structure of the low-α disc.