Disk Evolution Research Articles

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.

Impact of neutron star spin on Poynting–Robertson drag during a type I X-ray burst

ABSTRACT External irradiation of a neutron star (NS) accretion disc induces Poynting–Robertson (PR) drag, removing angular momentum and increasing the mass accretion rate. Recent simulations show PR drag significantly enhancing the mass accretion rate during Type I X-ray bursts, which could explain X-ray spectral features such as an increase in the persistent emission and a soft excess. However, prograde spin of the NS is expected to weaken PR drag, challenging its importance during bursts. Here, we study the effect of spin on PR drag during X-ray bursts. We run four simulations, with two assuming a non-spinning NS and two using a spin parameter of a* = 0.2, corresponding to a rotation frequency of 500 Hz. For each scenario, we simulate the disc evolution subject to an X-ray burst and compare it to the evolution found with no burst. PR drag drains the inner disc region during a burst, moving the inner disc radius outwards by ≈1.6 km in the a* = 0 and by ≈2.2 km in the a* = 0.2 simulation. The burst enhances the mass accretion rate across the innermost stable circular orbit ≈7.9 times when the NS is not spinning and ≈11.2 times when it is spinning. The explanation for this seemingly contradictory result is that the disc is closer to the NS when a* = 0.2, and the resulting stronger irradiating flux offsets the weakening effect of spin on the PR drag. Hence, PR drag remains a viable explanation for the increased persistent emission and soft excess observed during X-ray bursts in spinning NS systems.

Linear and non-linear eccentric mode evolution in unstratified MHD discs

ABSTRACT In this paper, we develop a framework for studying unstratified, magnetized eccentric discs and compute uniformly precessing eccentric modes in a cylindrical annulus which provide convenient initial conditions for numerical simulations. The presence of a magnetic field in an eccentric disc can be described by an effective gas with a modified equation of state. At magnetic field strengths relevant to the magneto-rotational instability the magnetic field has negligible influence on the evolution of the eccentric disc, however, the eccentric disc can significantly enhance the magnetic field strength over that in the a circular disc. We verify the suitability of these eccentric disc solutions by carrying out 2D simulations in ramses. Our simulated modes (in 2D) follow a similar evolution to the purely hydrodynamical modes, matching theoretical expectations, provided they are adequately resolved. Such solutions will provide equilibrium states for studies of the eccentric magneto-rotational instability and magnetized parametric instability in unstratified discs and are useful for exploring the response of disc turbulence on top of a fluid flow varying on the orbital time-scale.

Metal enrichment due to embedded stars in AGN discs

ABSTRACT We separately assess elemental abundances in active galactic nuclei's (AGNs) broad and narrow emission line regions (BLR and NLR), based on a critical assessment of published results together with new photoionization models. We find (1) He/H enhancements in some AGN, exceeding what can be explained by normal chemical evolution and confirm, (2) super-solar α abundance, though to a lesser degree than previously reported. We also reaffirm, (3) an N/O ratio consistent with secondary production, (4) solar or slightly sub-solar Fe abundance, and (5) red-shift independent metallicity, in contrast with galactic chemical evolution. We interpret (6) the larger metallicity in the BLR than NRL in terms of an in situ stellar evolution and pollution in AGN discs (SEPAD) model. We attribute (a) the redshift independence to the heavy element pollutants being disposed into the disc and accreted onto the central supermassive black hole (SMBH), (b) the limited He excess to the accretion–wind metabolism of a top-heavy population of evolving massive main sequence stars, (c) the super-solar CNO enrichment to the nuclear synthesis during their post-main-sequence evolution, (d) the large N/O to the byproduct of multiple stellar generations, and (e) the Mg, Si, and Fe to the ejecta of type II supernovae in the disc. These results provide supporting evidence for (f) ongoing self-regulated star formation, (g) adequate stellar luminosity to maintain marginal gravitational stability, (h) prolific production of seeds, and (i) dense coexistence of subsequently grown residual black hole populations in AGN discs.

Close-in ice lines and the super-stellar C/O ratio in discs around very low-mass stars

The origin of the elevated C/O ratios in discs around late M dwarfs compared to discs around solar-type stars is not well understood. Here we endeavour to reproduce the observed differences in the disc C/O ratios as a function of stellar mass using a viscosity-driven disc evolution model and study the corresponding atmospheric composition of planets that grow inside the water-ice line in these discs. We carried out simulations using a coupled disc evolution and planet formation code that includes pebble drift and evaporation. We used a chemical partitioning model for the dust composition in the disc midplane. Inside the water-ice line, the disc’s C/O ratio initially decreases to sub-stellar due to the inward drift and evaporation of water-ice-rich pebbles before increasing again to super-stellar values due to the inward diffusion of carbon-rich vapour. We show that this process is more efficient for very low-mass stars compared to solar-type stars due to the closer-in ice lines and shorter disc viscous timescales. In high-viscosity discs, the transition from sub-stellar to super-stellar takes place faster due to the fast inward advection of carbon-rich gas. Our results suggest that planets accreting their atmospheres early (when the disc C/O is still sub-stellar) will have low atmospheric C/O ratios, while planets that accrete their atmospheres late (when the disc C/O has become super-stellar) can obtain high C/O ratios. Our model predictions are consistent with observations, under the assumption that all stars have the same metallicity and chemical composition, and that the vertical mixing timescales in the inner disc are much shorter than the radial advection timescales. This further strengthens the case for considering stellar abundances alongside disc evolution in future studies that aim to link planet (atmospheric) composition to disc composition.

Evolution of fallback discs around millisecond magnetars: effect of supercritical accretion on GRB afterglows

ABSTRACT Some models of gamma-ray bursts (GRBs) invoke nascent millisecond magnetars as the central engine and address the X-ray afterglows with the interaction of magnetar magnetospheres with fallback discs. We study the evolution of fallback discs interacting with the millisecond magnetars. Initially, the accretion rate in the fallback disc is very high, well above the rate required for the Eddington limit. The inner parts of such a disc, even if it is cooling by the neutrino emission, get spherical due to the radiation pressure, which regulates the mass-accretion rate within the spherization radius. Such a disc can not penetrate the light cylinder radius for the typical magnetic fields, and the initial spin frequencies invoked for the magnetars. As a result of the autoregulation of the accretion flow, the fallback disc can not interact directly with the magnetar’s magnetosphere within the first few days. This has implications for the fallback disc models of GRB afterglows since the accretion and propeller luminosities, in the presence of radiation pressure, are too low to address the typical luminosities of X-ray afterglows.

Do all gaps in protoplanetary discs host planets?

Following the assumption that the disc substructures observed in protoplanetary discs originate from the interaction between the disc and the forming planets embedded therein, we aim to test if these putative planets could represent the progenitors of the currently observed giant exoplanets. We performed N-body simulations initially assuming three, four, five, or seven planets. Our model includes pebble and gas accretion, migration, damping of eccentricities and inclinations, disc-planet interaction, and disc evolution. We located the planets in the positions where the gaps in protoplanetary discs have been observed and we evolved the systems for 100 Myr including a few million years of gas disc evolution, while also testing three values of α viscosity. For planetary systems with initially three and four planets, we find that most of the growing planets lie beyond the radial-velocity (RV) detection limit of 5AU and only a small fraction of them migrate into the inner region. We also find that these systems have final eccentricities that are too low to be in agreement with the observed giant planet population. Systems initially consisting of five or seven planets become unstable after ≈40 Kyr of integration time. This clearly shows that not every gap can host a planet. The general outcome of our simulations – eccentricities that are too low – is independent of the disc’s viscosity and surface density. Further observations could either confirm the existence of an undetected population of wide-orbit giants or exclude the presence of such an undetected population to constrain how many planets hide in gaps even further.

High-resolution [O i] Line Spectral Mapping of TW Hya Consistent with X-Ray-driven Photoevaporation

Theoretical models indicate that photoevaporative and magnetothermal winds play a crucial role in the evolution and dispersal of protoplanetary disks and affect the formation of planetary systems. However, it is still unclear what wind-driving mechanism is dominant or if both are at work, perhaps at different stages of disk evolution. Recent spatially resolved observations by Fang et al. of the [O i] 6300 Å spectral line, a common disk wind tracer in TW Hya, revealed that about 80% of the emission is confined to the inner few astronomical units of the disk. In this work, we show that state-of-the-art X-ray-driven photoevaporation models can reproduce the compact emission and the line profile of the [O i] 6300 Å line. Furthermore, we show that the models also simultaneously reproduce the observed line luminosities and detailed spectral profiles of both the [O i] 6300 Å and the [Ne ii] 12.8 μm lines. While MHD wind models can also reproduce the compact radial emission of the [O i] 6300 Å line, they fail to match the observed spectral profile of the [O i] 6300 Å line and underestimate the luminosity of the [Ne ii] 12.8 μm line by a factor of 3. We conclude that, while we cannot exclude the presence of an MHD wind component, the bulk of the wind structure of TW Hya is predominantly shaped by a photoevaporative flow.

An inflationary disk phase to explain extended protoplanetary dust disks

Context. Understanding planetesimal formation is an essential first step towards understanding planet formation. The distribution of these first solid bodies drives the locations where planetary embryos, which eventually form fully-fledged planets, grow. Aims. We seek to understand the parameter space of possible protoplanetary disk formation and evolution models of our Solar System. A good protoplanetary disk scenario for the Solar System must meet at least the following three criteria: (1) It must produce an extended gas and dust disk (e.g. 45 au for the dust); (2) within the disk, the local dust-to-gas ratio in at least two distinct locations must sufficiently increase to explain the early formation of the parent bodies of non-carbonaceous and carbonaceous iron meteorites; and (3) dust particles, which have condensed at high temperatures (i.e. calcium–aluminium-rich inclusions), must be transported to the outer disk. Though current protoplanetary disk models are able to satisfy one or two of these criteria, none have been successful in recreating all three. We aim to find scenarios that satisfy all three. Methods. In this study we used a 1D disk model that tracks the evolution of the gas and dust disks. Planetesimals are formed within the disk at locations where the streaming instability can be triggered. We explored a large parameter space to study the effect of the disk viscosity, the timescale of infall of material into the disk, the distance within which material is deposited into the disk, and the fragmentation threshold of dust particles. Results. We find that scenarios with a large initial disk viscosity (α > 0.05), a relatively short infall timescale (Tinfall < 100–200kyr), and a small centrifugal radius (RC ~ 0.4 au; i.e. the distance within which material falls into the disk) result in disks that satisfy all three criteria needed to represent the protoplanetary disk of the Solar System. The large initial viscosity and short infall timescale result in a rapid initial expansion of the disk, which we dub the ‘inflationary phase’ of the disk. Furthermore, a temperature-dependent fragmentation threshold, which accounts for cold icy particles breaking more easily, results in larger and more massive disks. This, in turn, results in more ‘icy’ than ‘rocky’ planetesimals. Such scenarios are also better in line with our Solar System, which has small terrestrial planets and massive giant planet cores. Finally, we find that scenarios with large RC cannot transport calcium–aluminium-rich inclusions to the outer disk and do not produce planetesimals at two locations within the disk.

The Effect of Dust Evolution and Traps on Inner Disk Water Enrichment

Abstract Substructures in protoplanetary disks can act as dust traps that shape the radial distribution of pebbles. By blocking the passage of pebbles, the presence of gaps in disks may have a profound effect on pebble delivery into the inner disk, crucial for the formation of inner planets via pebble accretion. This process can also affect the delivery of volatiles (such as H2O) and their abundance within the water snow line region (within a few au). In this study, we aim to understand what effect the presence of gaps in the outer gas disk may have on water vapor enrichment in the inner disk. Building on previous work, we employ a volatile-inclusive disk evolution model that considers an evolving ice-bearing drifting dust population, sensitive to dust traps, which loses its icy content to sublimation upon reaching the snow line. We find that the vapor abundance in the inner disk is strongly affected by the fragmentation velocity (v f) and turbulence, which control how intense vapor enrichment from pebble delivery is, if present, and how long it may last. Generally, for disks with low to moderate turbulence (α ≤ 1 × 10−3) and a range of v f, radial locations and gap depths (especially those of the innermost gaps) can significantly alter enrichment. Shallow inner gaps may continuously leak material from beyond it, despite the presence of additional deep outer gaps. We finally find that for realistic v f (≤10 m s−1), the presence of gaps is more important than planetesimal formation beyond the snow line in regulating pebble and volatile delivery into the inner disk.

How Large Is a Disk—What Do Protoplanetary Disk Gas Sizes Really Mean?

It remains unclear what mechanism is driving the evolution of protoplanetary disks. Direct detection of the main candidates, either turbulence driven by magnetorotational instabilities or magnetohydrodynamical disk winds, has proven difficult, leaving the time evolution of the disk size as one of the most promising observables able to differentiate between these two mechanisms. But to do so successfully, we need to understand what the observed gas disk size actually traces. We studied the relation between R CO,90%, the radius that encloses 90% of the 12CO flux, and R c , the radius that encodes the physical disk size, in order to provide simple prescriptions for conversions between these two sizes. For an extensive grid of thermochemical models, we calculate R CO,90% from synthetic observations and relate properties measured at this radius, such as the gas column density, to bulk disk properties, such as R c and the disk mass M disk. We found an empirical correlation between the gas column density at R CO,90% and disk mass: . Using this correlation we derive an analytical prescription of R CO,90% that only depends on R c and M disk. We derive R c for disks in Lupus, Upper Sco, Taurus, and the DSHARP sample, finding that disks in the older Upper Sco region are significantly smaller (〈R c 〉 = 4.8 au) than disks in the younger Lupus and Taurus regions (〈R c 〉 = 19.8 and 20.9 au, respectively). This temporal decrease in R c goes against predictions of both viscous and wind-driven evolution, but could be a sign of significant external photoevaporation truncating disks in Upper Sco.

Accretion disc evolution in GRO J1655−40 and LMC X-3 with relativistic and non-relativistic disc models

ABSTRACT Black hole X-ray binaries are ideal environments to study the accretion phenomena in strong gravitational potentials. These systems undergo dramatic accretion state transitions and analysis of the X-ray spectra is used to probe the properties of the accretion disc and its evolution. In this work, we present a systematic investigation of ∼1800 spectra obtained by Rossi X-Ray Timing Explorer Proportional Counter Array observations of GRO J1655−40 and LMC X-3 to explore the nature of the accretion disc via non-relativistic and relativistic disc models describing the thermal emission in black hole X-ray binaries. We demonstrate that the non-relativistic modelling throughout an outburst with the phenomenological multicolour disc model DISKBB yields significantly lower and often unphysical inner disc radii and correspondingly higher (∼50–60 per cent) disc temperatures compared to its relativistic counterparts KYNBB and KERRBB. We obtained the dimensionless spin parameters of a* = 0.774 ± 0.069 and a* = 0.752 ± 0.061 for GRO J1655−40 with KERRBB and KYNBB, respectively. We report a spin value of a* = 0.098 ± 0.063 for LMC X-3 using the updated black hole mass of 6.98 M⊙. Both measurements are consistent with the previous studies. Using our results, we highlight the importance of self-consistent modelling of the thermal emission, especially when estimating the spin with the continuum-fitting method which assumes the disc terminates at the innermost stable circular orbit at all times.

Constraining the Evolution of the Unstable Accretion Disk in SMC X-1 with NICER

Neutron star high-mass X-ray binaries with superorbital modulations in luminosity host warped inner accretion disks that occult the neutron star during precession. In SMC X-1, the instability in the warped disk geometry causes superorbital period “excursions”: times of instability when the superorbital period decreases from its typical value of 55 to ∼40 days. Disk instability makes SMC X-1 an ideal system in which to investigate the effects of variable disk geometry on the inner accretion flow. Using the high-resolution spectral and timing capabilities of the Neutron Star Interior Composition Explorer, we examined the high state of four different superorbital cycles of SMC X-1 to search for changes in spectral shape and connections to the unstable disk geometry. We performed pulse phase-averaged and phase-resolved spectroscopy to closely compare the changes in spectral shape and any cycle-to-cycle variations. While some parameters, including the photon index and absorbing column density, show slight variations with superorbital phase, these changes are most evident during the intermediate state of the superorbital cycle. Few spectral changes are observed within the high state of the superorbital cycle, possibly indicating the disk instability does not significantly change SMC X-1's accretion process.

How accretion of planet-forming disks influences stellar abundances

Millimeter-sized dust grains experience radial velocities exceeding the gas velocities by orders of magnitude. The viscous evolution of the accretion disk adds disk material onto the central star’s convective envelope, influencing its elemental abundances, [X/H]. At the same time, the envelope mass shrinks as the stellar age increases, amplifying the rate of abundance change. Therefore, the elemental abundances of the star are sensitive to disk processes that alter the composition and timing of disk accretion. We performed numerical 1D log-radial simulations integrating the disk advection-diffusion equation, while accounting for evaporation and condensation of chemical species at the evaporation fronts. They reveal a peak of refractory abundance within the first 2 Myr of Δ[X/H] ~ 5 × 10−2 if grain growth is significant, but subsequent accretion diminishes previous refractory abundance increases for long-lived disks. Planet formation can reduce the abundance of dust species whose evaporation fronts lie within the planet’s orbit by preventing solids from reaching the inner edge once the planet starts opening a gap exerting a pressure bump exterior to its orbit and consequently blocking inward drifting pebbles. We expect the accretion of the solar protoplanetary disk with Jupiter present to have changed the Sun’s elemental abundances by ~1 × 10−2 throughout its lifetime. These considerations were also applied to the HD106515 wide binary system. We find that measurements of Δ[X/H] are in reasonable agreement with results from simulations where the observed giant planet around HD106515 A is included and if HD106515B’s disk formed planetesimals more efficiently. Simulations where the planet formed inside the water ice line are more favorable to agree with observations. Even though the general changes in the stellar abundances due to disk accretion are small, they are detectable at current sensitivities, indicating that the methods presented here can be used to constrain the planet formation pathway.