Variable Accretion Research Articles

Direct Collapse Accretion Disks within Dark Matter Halos: Saturation of the Magnetorotational Instability and the Field Expulsion

We have used high-resolution zoom-in simulations of direct collapse to supermassive black hole (SMBH) seeds within dark mater halos in the presence of magnetic fields generated during the collapse, down to 10−5 pc or 2 au. We confirm an efficient amplification of magnetic field during collapse, the formation of a geometrically thick self-gravitating accretion disk inside 0.1 pc, and damping of fragmentation in the disk by the field. This disk differs profoundly from SMBH accretion disks. We find the following: (1) The accretion disk is subject to the magnetorotational instability, which further amplifies the field to near equipartition. No artificial seeding of the disk field has been used. (2) The equipartition toroidal field changes its polarity in the midplane. (3) The nonlinear Parker instability develops, accompanied by the vertical buckling of the field lines, which injects material above the disk, leading to an increase in the disk scale height. (4) With the Coriolis force producing a coherent helicity above the disk, the vertical poloidal field has been generated and amplified. (5) We estimate that the associated outflow will be most probably squashed by accretion. The resulting configuration consists of a magnetized disk with β ≳ 0.1 and its magnetosphere with β ≪ 1, where β = P th/P B is the ratio of thermal to magnetic energy density. (6) The disk is highly variable, due to feeding by variable accretion flow, and strong vortical motions are present. (7) Finally, the negative gradient of the total vertical stress drives an equatorial outflow sandwiched by an inward accretion flow.

ExoplaNeT accRetion mOnitoring sPectroscopic surveY (ENTROPY)

Context. Accretion among planetary mass companions is a poorly understood phenomenon, due to the lack of both observational and theoretical studies. The detection of emission lines from accreting gas giants facilitates detailed investigations into this process. Aims. This work presents a detailed analysis of Balmer lines from one of the few known young, planetary-mass objects with observed emission, the isolated L2γ dwarf 2MASS J11151597+1937266 with a mass between 7 and 21 MJup and an age of 5–45 Myr, located at 45 ± 2 pc. Methods. We obtained the first high-resolution (R ~ 50 000) spectrum of the target with VLT/UVES, an echelle spectrograph operating in the near-ultraviolet to visible wavelengths (3200–6800 Å). Results. We report several resolved hydrogen (H I; H3–H6) and helium (He I; λ5875.6) emission lines in the spectrum. Based on the asymmetric line profiles of Hα and Hβ, the 10% width of Hα (199 ± 1 km s−1), tentative He I λ6678 emission, and indications of a disk from mid-infrared excess, we confirm ongoing accretion at this object. Using the Gaia update of the parallax, we revise its temperature to 1816 ± 63 K and radius to 1.5 ± 0.1 RJup. Analysis of observed H I profiles using a 1D planet-surface shock model implies a pre-shock gas velocity, v0 = 120−40+ 80 km s−1, and a pre-shock density, log(n0/cm−3) = 14−5+ 0. The pre-shock velocity points to a mass, Mp = 6−4+ 8 MJup, for the target. Combining H I line luminosities (Lline) and planetary Lline−Lacc (accretion luminosity) scaling relations, we derived a mass accretion rate, Ṁacc = 1.4−0.9+ 2.8 × 10−8 MJup yr−1. Conclusions. The line-emitting area predicted from the planet-surface shock model is very small (~0.03%), and points to a shock at the base of a magnetospherically induced funnel. The Hα profile exhibits a much stronger flux than predicted by the model that best fits the rest of the H I profiles, indicating that another mechanism than shock emission contributes to the Hα emission. Comparison of line fluxes and Ṁacc from archival moderate-resolution SDSS spectra indicate variable accretion at 2MASS J11151597+1937266.

The effects of stellar and AGN feedback on the cosmic star formation history in the simba simulations

ABSTRACT Using several variants of the cosmological simba simulations, we investigate the impact of different feedback prescriptions on the cosmic star formation history. Adopting a global-to-local approach, we link signatures seen in global observables, such as the star formation rate density (SFRD) and the galaxy stellar mass function (GSMF), to feedback effects in individual galaxies. We find a consistent picture: stellar feedback mainly suppresses star formation below halo masses of $M_{\rm H} = 10^{12} \rm \, {\rm M}_{\odot }$ and before $z = 2$, whereas AGN feedback quenches the more massive systems after $z = 2$. Among simba’s AGN feedback modes, AGN jets are the dominant quenching mechanism and set the shape of the SFRD and the GSMF at late times. AGN-powered winds only suppress the star formation rate in intermediate-mass galaxies ($M_{\rm \star } = 10^{9.5 - 10} \rm \, {\rm M}_{\odot }$), without affecting the overall stellar mass-assembly significantly. At late times, the AGN X-ray feedback mode mainly quenches residual star formation in massive galaxies. Our analysis reveals that this mode is also necessary to produce the first fully quenched galaxies before $z=2$, where the jets alone are inefficient. These initially highly star-forming galaxies contain relatively large black holes, likely strengthening the X-ray-powered heating and ejection of gas from the dense, central region of galaxies. Such extra heating source quenches the local star formation and produces a more variable accretion rate. More generally, this effect also causes the break down of correlations between the specific star formation rate, the accretion rate and the black hole mass.

Highly variable magmatic accretion at the ultraslow-spreading Gakkel Ridge

Crustal accretion at mid-ocean ridges governs the creation and evolution of the oceanic lithosphere. Generally accepted models1–4 of passive mantle upwelling and melting predict notably decreased crustal thickness at a spreading rate of less than 20 mm year−1. We conducted the first, to our knowledge, high-resolution ocean-bottom seismometer (OBS) experiment at the Gakkel Ridge in the Arctic Ocean and imaged the crustal structure of the slowest-spreading ridge on the Earth. Unexpectedly, we find that crustal thickness ranges between 3.3 km and 8.9 km along the ridge axis and it increased from about 4.5 km to about 7.5 km over the past 5 Myr in an across-axis profile. The highly variable crustal thickness and relatively large average value does not align with the prediction of passive mantle upwelling models. Instead, it can be explained by a model of buoyant active mantle flow driven by thermal and compositional density changes owing to melt extraction. The influence of active versus passive upwelling is predicted to increase with decreasing spreading rate. The process of active mantle upwelling is anticipated to be primarily influenced by mantle temperature and composition. This implies that the observed variability in crustal accretion, which includes notably varied crustal thickness, is probably an inherent characteristic of ultraslow-spreading ridges.

What drives the variability in AGN? Explaining the UV-Xray disconnect through propagating fluctuations

ABSTRACT Intensive broad-band reverberation mapping campaigns have shown that AGN variability is significantly more complex than expected from disc reverberation of the variable X-ray illumination. The UV/optical variability is highly correlated and lagged, with longer lags at longer wavelengths as predicted, but the observed time-scales are longer than expected. Worse, the UV/optical light curves are not well correlated with the X-rays, which should drive them. Instead, we consider an intrinsically variable accretion disc, where slow mass accretion rate fluctuations are generated in the optical-UV disc, propagating down to modulate intrinsically faster X-ray variability from the central regions. We match our model to Fairall 9, a well-studied AGN with L ∼ 0.1LEdd, where the spectrum is dominated by the UV/EUV. Our model produces light curves where the X-rays and UV have very different fast variability, yet are well correlated on longer time-scales, as observed. It predicts that the intrinsic variability has optical/UV leading the X-rays, but including reverberation of the variable EUV from an inner wind produces a lagged bound-free continuum that matches the observed UV-optical lags. We conclude that optical/UV AGN variability is likely driven by intrinsic fluctuations within the disc, not X-ray reprocessing: the observed longer than expected lags are produced by reverberation of the EUV illuminating a wind, not by X-ray illumination of the disc: the increasing lag with increasing wavelength is produced by the increased contribution of the (constant lag) bound-free continuum to the spectrum, rather than indicating intrinsically larger reverberation distances for longer wavelengths.

The evolution and impact of ∼3000 M⊙ stars in the early Universe

We present evolutionary models of massive, accreting population III stars with constant and variable accretion rates until the end of silicon burning, with final masses of ~ 1000–3000 M⊙. In all our models, after the core-hydrogen-burning phase, the star expands towards the red side of the Hertzsprung-Russell diagram is where it spends the rest of its evolution. During core helium burning, the models exhibit an outer convective envelope as well as many large intermediate convective zones. These intermediate zones allow for strong internal mixing to occur which enriches the surface in helium. The effect of increasing metallicity at a constant accretion rate of 10−3 M⊙yr−1 shows an increase in the lifetime, final mass and distribution of helium in the envelope. Our fiducial model with mass of 3000 M⊙ has a final surface helium abundance of 0.74 and 9% of its total mass or 50% of the core mass, has a value of Γ1 < 4/3 at the end of core silicon burning. If the collapse of the core is accompanied by the ejection of the envelope above the carbon-oxygen core, this could have a significant impact on the chemical evolution of the surroundings and subsequent stellar generations. The model has a final log(N/O) ≈ 0.45, above the lower limit in the recently detected high-redshift galaxy GN-z11. We discuss the impact of a single 3000 M⊙ star on chemical, mechanical and radiative feedback, and present directions for future work.

Morphologic Adjustment of a River Reach with Groynes to Channel Bypassing

AbstractThis article is focused on the investigation of the spatio-temporal variability of the Danube River reach’s vertical accretion thickness due to the response of the Danube River reach to bypassing. Five groyne-induced benches (GIBs) of the bypassed channel developed after water diversion in 1992 was studied by changes in topography for three-time spans (for the original gravel surface, for the surface before the 2013 flood and for the surface after the 2013 flood). The allostratigraphic approach was applied to 548 drilling probes at all GIBs and toptop, supra-platform, tail and backchannel geomorphic units have been identified at each GIB. The main to side-channel system connectivity increase sedimentation rates and the accretion was controlled by large flood events. The 100-year flood in 2013 contributed to the total volume by almost 26%. During study period 1992–2017, totally 1,146,589 m3 was accreted on five GIBs, of this 209,752 m3 during flood event in 2013 and 267,700 m3 in post flood period 2014–2017. The top geomorphic unit exhibits the highest median values of vertical accretion and for all GIBs accretion thickness are not related to the height above the mean channel water level. The thickness of accretion varied, likely because the variability of the vegetation cover caused variable hydraulic conditions and accretion rate span from 3.8 cm.year−1 to 5.3 cm.year−1. The investigation of the sediment thickness over long time spans 24 years and a large flood event) allowed us to conclude that the thickness is spatially variable for individual GIBs; however, its trend over time remains constants depending on the intake of sediments during large floodsd events. This article also provides a methodological template for the detailed investigation of river channel adjustment due to bypassing.

Shedding Light on the Ejection History of Molecular Outflows: Multiple Velocity Modes and Precession

Variable accretion has been well studied in the evolved stages of low-mass star formation. However, the accretion history in the initial phases of star formation is still a seldom studied topic. The outflows and jets emerging from protostellar objects could shed some light on their accretion history. We consider the recently studied case of W43-MM1, a protocluster containing 46 outflows driven by 27 protostellar cores. The outflow kinematics of the individual cores and associated knots in W43-MM1 indicate episodic protostellar ejection. We take the observed parameters of an individual core system (core #8) and perform 3D hydrodynamic simulations of such a system, including episodic changes in the velocity of the outflow. The simulations consist of a collimated jet emerging from a core, taking into account one- and two-velocity modes in the variation of the ejection velocity of the jet. In addition, we investigated the effect of including the precession of the jet in the one- and two-velocity-mode models. From the simulations, we constructed position–velocity diagrams and compared them with the observations. We find that including a second mode in the ejection velocity, as well as the precession, are required to explain the positions of the outflow knots and other position–velocity features observed in core #8 in W43-MM1.

The birth of backbarrier marshes in Culatra Island (Ria Formosa, South Portugal)

The rapid elongation of Culatra Island, a sandy barrier in the Ria Formosa chain (S. Portugal), since the mid-1940s led to the formation of three new embayments in its backbarrier that were gradually colonised by halophytic vegetation. This provided a rare opportunity to collect information and data on the very early stages of backbarrier marsh plant establishment and evolution. Sediment (surface and subsurface) sampling in two of the recently formed bays, combined with information extracted from vertical aerial photographs, allowed us to assess modern sedimentation characteristics and vertical accretion rates since the shift from a bare sandflat to a vegetated marsh platform. Present-day topography appears largely inherited by overwash or/and inlet-related tidal deposits that provided the necessary sediment pulse for the formation of an intertidal sandy substrate, suitable for colonisation. The variability in accretion rates, noted even within the same embayment, as well as the differences in accretion balance with similarly young backbarrier marshes, highlight the importance of local conditions (sediment import, distance to creeks and marsh edge, storm frequency and intensity) to marsh build-up, even during the very early stages. Variable accretion rates were also identified over intertidal seagrass patches, indicating similar influences. Organic deposition rates were very low in all vegetated intertidal habitats, indicating the dominance of mineral deposition to the vertical growth. A lag, ranging from roughly 10–30 years, was observed between the formation of the intertidal sandy platform and plant establishment in all embayments. The different timescales in the observed lag are likely linked to differences in hydrodynamic conditions, promoted by the embayment morphology (opening width). The lowest lag was observed in protected embayments, which could reflect a ‘typical’ delay for plant establishment in the system, while the highest lag was associated with higher energy backbarrier environments.

Shock-driven periodic variability in a low-mass-ratio supermassive black hole binary

ABSTRACT We investigate the time-varying electromagnetic emission of a low-mass-ratio supermassive black hole binary (SMBHB) embedded in a circumprimary disc, with a particular interest in variability of shocks driven by the binary. We perform a 2D, locally isothermal hydrodynamics simulation of an SMBHB with mass ratio q = 0.01 and separation a = 100 Rg, using a physically self-consistent steady disc model. We estimate the electromagnetic variability from the system by monitoring accretion on to the secondary and using an artificial viscosity scheme to capture shocks and monitor the energy dissipated. The SMBHB produces a wide, eccentric gap in the disc, previously only observed for larger mass ratios, which we attribute to our disc model being much thinner (H/R ≈ 0.01 near the secondary) than is typical of previous works. The eccentric gap drives periodic accretion on to the secondary SMBH on a time-scale matching the orbital period of the binary, $t_{\rm {bin}}\approx 0.1\,\,\rm {yr}$, implying that the variable accretion regime of the SMBHB parameter space extends to lower mass ratios than previously established. Shocks driven by the binary are periodic, with a period matching the orbital period, and the shocks are correlated with the accretion rate, with peaks in the shock luminosity lagging peaks in the accretion rate by 0.43 tbin. We propose that the correlation of these quantities represents a useful identifier of SMBHB candidates, via observations of correlated variability in X-ray and ultraviolet monitoring of active galactic nuclei, rather than single-waveband periodicity alone.

Critical accretion rates for rapidly growing massive Population III stars

Efforts to understand the origin and growth of massive black holes observed in the early Universe have spurred significant interest in the evolution and fate of rapidly accreting primordial (metal-free) stars. Here, we investigate the evolution of such Population III (PopIII) stars under variable accretion rates, focusing on the thermal response and stellar structure, the impact of the luminosity wave encountered early in the pre-main sequence (pre-MS) phase, and the influence of accretion on their subsequent evolution. We employ the Geneva stellar evolution code and simulate ten models with varying accretion histories, covering a final mass range from 491 M⊙ to 6127 M⊙. Our findings indicate that the critical accretion rate delineating the red and blue supergiant regimes during the pre-MS evolution is approximately 2.5 × 10−2 M⊙ yr−1. Once core hydrogen burning commences, the value of this critical accretion rate drops to 7.0 × 10−3 M⊙ yr−1. Moreover, we also confirm that the Kelvin–Helmholtz timescale in the outer surface layers is the more relevant timescale for determining the transition between red and blue phases. Regarding the luminosity wave, we find that it affects only the early pre-MS phase of evolution and does not directly influence the transition between red and blue phases, which primarily depends on the accretion rate. Finally, we demonstrate that variable accretion rates significantly impact the lifetimes, surface enrichment, and final mass of the PopIII stars, as well as the time they spend in the red phase. Our study provides a comprehensive understanding of the intricate evolutionary patterns of PopIII stars subjected to variable accretion rates.

A new 2D stochastic methodology for simulating variable accretion discs: propagating fluctuations and epicyclic motion

ABSTRACT Accretion occurs across a large range of scales and physical regimes. Despite this diversity in the physics, the observed properties show remarkable similarity. The theory of propagating fluctuations, in which broad-band variability within an accretion disc travel inwards and combine, has long been used to explain these phenomena. Recent numerical work has expanded on the extensive analytical literature but has been restricted to using the 1D diffusion equation for modelling the disc behaviour. In this work we present a novel numerical approach for 2D (vertically integrated), stochastically driven α-disc simulations, generalizing existing 1D models. We find that the theory of propagating fluctuations translates well to 2D. However, the presence of epicyclic motion in 2D (which cannot be captured within the diffusion equation) is shown to have an important impact on local disc dynamics. Additionally, there are suggestions that for sufficiently thin discs the log-normality of the light curves changes. As in previous work, we find that the break frequency in the luminosity power spectrum is strongly dependent on the driving time-scale of the stochastic perturbations within the disc, providing a possible observational signature for probing the magnetorotational instability dynamo. We also find that thinner discs are significantly less variable than thicker ones, providing a compelling explanation for the greater variability seen in the hard state versus the soft state of X-ray binaries. Finally, we consider the wide-ranging applications of our numerical model for use in other simulations.

Photometric and spectroscopic monitoring of YSOs in nearby star-forming regions – I. Eruptive YSOs

ABSTRACTMid-infrared (mid-IR) variability in young stellar objects (YSOs) is driven by several physical mechanisms, which produce a variety of amplitudes and light-curve shapes. One of these mechanisms, variable disc accretion, is predicted by models of episodic accretion to drive secular variability, including in the mid-IR. Because the largest accretion bursts are rare, adding new objects to the YSO eruptive variable class aids our understanding of the episodic accretion phenomenon and its possible impact on stellar and planetary formation. A previous analysis of 6.5 yr of NeoWISE light curves (3–5 μm) of ∼7000 nearby YSOs found an increase in the fraction of variability and variability amplitude for objects at younger stages of evolution. To help interpret these light curves, we have obtained low- and high-resolution near-IR spectra of 78 objects from this sample of YSOs. In this work, we present the analysis of nine nearby YSOs (d<1 kpc) that show the characteristics of known classes of eruptive variable YSOs. We find one FU Orionis (FUor)-like source, one EX Lupi type object, and six YSOs with mixed characteristics or V1647 Ori like objects. The varied characteristics observed in our sample are consistent with recent discoveries of eruptive YSOs. We discuss how a wide range in YSO outburst parameters (central mass, maximum accretion rate during outburst, evolutionary stage, and/or instability leading to the outburst) may play a significant role in the observed spectrophotometric properties of YSO outbursts.

NEOWISE Variability of the New Gaia Herbig Ae/Be and Classical Be Stars—Variability Sample and Statistical Analysis

Light variability is an essential attribute of young stellar objects (YSOs) at the pre-main-sequence (PMS) stage, which can be treated as an indicator to trace star formation. We present a statistical study on the infrared NEOWISE variability of new and large-sized PMS stars as well as Herbig Ae/Be (HAeBe) and classical Be (CBe) stars identified from Gaia Data Release 2. A total of 3792 (including 714 HAeBe), 417, and 473 NEOWISE variable sources were detected from the catalog of 8470 PMS (including 1361 HAeBe), 693 CBe, and 1309 uncertain (possibly belonging to either HAeBe or CBe) stars, respectively. We found that 80% of the PMS variables are dominated by irregular infrared variability, whereas variables in the CBe and uncertain categories exhibit comparable contributions from secular (linear, curved, and periodic) and stochastic (burst, drop, and irregular) variables. This, in combination with infrared W1–W2 colors, W4 luminosities, and the spectral energy distribution slope at near- and mid-infrared wavelengths, suggests different evolutionary stages of the PMS (including HAeBe), uncertain, and CBe stars, from less to more evolved. A statistical analysis of variable source properties (e.g., period and fractional amplitude) further reveals a potential sequence possibly associated with disk evolution for different types of light curves (from linear, to curved, periodic, stochastic, and nonvariable stars). Possible mechanisms for the variable sources are suggested to be related to either extinction via circumstellar matter, or contribution of multiple systems surrounding the YSOs, as opposed to variable accretion (e.g., episodic accretion or accretion burst).

Line-driven winds from variable accretion discs

ABSTRACTWe use numerical hydrodynamics simulations to study line-driven winds launched from an accreting α-disc. Building on previous work where the driving radiation field is static, we compute a time-dependent radiation flux from the local, variable accretion rate of the disc. We find that prior to the establishment of a steady state in the disc, variations of $\sim 15{{\ \rm per\ cent}}$ in disc luminosity correlate with variations of ∼ 2–3 in the mass flux of the wind. After a steady state is reached, when luminosity variations drop to $\sim 3{{\ \rm per\ cent}}$, these correlations vanish as the variability in the mass flux is dominated by the intrinsic variability of the winds. This is especially evident in lower luminosity runs where intrinsic variability is higher due to a greater prevalence of failed winds. The changing mass flux occurs primarily due to the formation of clumps and voids near the disc atmosphere that propagate out into the low velocity part of the flow, a process that can be influenced by local variations in disc intensity. By computing the normalized standard deviation of the mass outflow, we show that the impact of luminosity variations on mass outflow is more visible at higher luminosity. However, the absolute change in mass outflow due to luminosity increases is larger for lower luminosity models due to the luminosity-mass flux scaling relation becoming steeper. We further discuss implications for CVs and AGN and observational prospects.

Directly Tracing Cool Filamentary Accretion over >100 kpc into the Interstellar Medium of a Quasar Host at z = 1

We report the discovery of giant (50−100 kpc) [O ii] emitting nebulae with MUSE in the field of TXS 0206−048, a luminous quasar at z = 1.13. “Down-the-barrel” UV spectra of the quasar show absorption at velocities coincident with those of the extended nebulae, enabling new insights into inflows and outflows around the quasar host. One nebula exhibits a filamentary morphology extending over 120 kpc from the halo toward the quasar and intersecting with another nebula surrounding the quasar host with a radius of 50 kpc. This is the longest cool filament observed to date and arises at higher redshift and in a less massive system than those in cool-core clusters. The filamentary nebula has line-of-sight velocities >300 km s−1 from nearby galaxies but matches that of the nebula surrounding the quasar host where they intersect, consistent with accretion of cool intergalactic or circumgalactic medium or cooling hot halo gas. The kinematics of the nebulae surrounding the quasar host are unusual and complex, with redshifted and blueshifted spiral-like structures. The emission velocities at 5−10 kpc from the quasar match those of inflowing absorbing gas observed in UV spectra of the quasar. Together, the extended nebulae and associated redshifted absorption represent a compelling case of cool, filamentary gas accretion from halo scales into the extended interstellar medium and toward the nucleus of a massive quasar host. The inflow rate implied by the combined emission and absorption constraints is well below levels required to sustain the quasar’s radiative luminosity, suggesting anisotropic or variable accretion.

Radiative feedback on supermassive star formation: the massive end of the Population III initial mass function

ABSTRACT Supermassive stars with masses of M* ≃ 104–105 M⊙ are invoked as possible seeds of high-redshift supermassive black holes, but it remains under debate whether their protostar indeed acquires sufficient mass via gas accretion overcoming radiative feedback. We investigate protostellar growth in dynamically heated atomic cooling haloes (ACHs) found in recent cosmological simulations, performing three-dimensional radiation hydrodynamical simulations that consider stellar evolution under variable mass accretion. We find that one of the ACHs feeds the central protostar at rates exceeding a critical value, above which the star evolves in a cool bloating phase and hardly produces ionizing photons. Consequently, the stellar mass reaches M* ≳ 104 M⊙ unimpeded by radiative feedback. In the other ACH, where the mass supply rate is lower, the star evolves almost as a hot main-sequence star, emitting intense ionizing radiation. Then, the stellar mass growth is terminated around 500 M⊙ by photoevaporation of the circumstellar disc. Our simulations provide a formula of the final stellar mass determined either by stellar feedback or their lifetime as a function of the mass supply rate from the parent cloud. Combining the results with the statistical properties of star-forming clouds in high-redshift quasar progenitor haloes, we construct a top-heavy mass distribution of primordial stars over M* ≃ 100–105 M⊙, approximately following a power-law spectrum of ${\propto} M_\ast ^{-1.3}$. Their black hole remnants would be further fed via the dense debris disc, powering ‘milliquasars’ with a bolometric luminosity of Lbol ≳ 1043 erg s−1.

Infrared variability of young solar analogues in the Lagoon Nebula

ABSTRACT T Tauri stars are low-mass pre-main sequence stars that are intrinsically variable. Due to the intense magnetic fields they possess, they develop dark spots on their surface that, because of rotation, introduce a periodic variation of brightness. In addition, the presence of surrounding discs could generate flux variations by variable extinction or accretion. Both can lead to a brightness decrease or increase, respectively. Here, we have compiled a catalogue of light curves for 379 T Tauri stars in the Lagoon Nebula (M8) region, using VVVX survey data in the Ks-band. All these stars were already classified as pre-MS stars based on other indicators. The data presented here are spread over a period of about eight years, which gives us a unique follow-up time for these sources at this wavelength. The light curves were classified according to their degree of periodicity and asymmetry, to constrain the physical processes responsible for their variation. Periods were compared with the ones found in literature, on a much shorter baseline. This allowed us to prove that for 126 stars, the magnetically active regions remain stable for several years. Besides, our near-IR data were compared with the optical Kepler/K2 light curves, when available, giving us a better understanding of the mechanisms responsible for the brightness variations observed and how they manifest at different bands. We found that the periodicity in both bands is in fairly good agreement, but the asymmetry will depend on the amplitude of the bursts or dips events and the observation cadence.

Simulations of black hole fueling in isolated and merging galaxies with an explicit, multiphase ISM

ABSTRACT We study gas inflows on to supermassive black holes using hydrodynamics simulations of isolated galaxies and idealized galaxy mergers with an explicit, multiphase interstellar medium (ISM). Our simulations use the recently developed ISM and stellar evolution model called Stars and MUltiphase Gas in GaLaxiEs (SMUGGLE). We implement a novel super-Lagrangian refinement scheme that increases the gas mass resolution in the immediate neighbourhood of the black holes (BHs) to accurately resolve gas accretion. We do not include black hole feedback in our simulations. We find that the complex and turbulent nature of the SMUGGLE ISM leads to highly variable BH accretion. BH growth in SMUGGLE converges at gas mass resolutions ≲3 × 103 M⊙. We show that the low resolution simulations combined with the super-Lagrangian refinement scheme are able to produce central gas dynamics and BH accretion rates very similar to that of the uniform high resolution simulations. We further explore BH fueling by simulating galaxy mergers. The interaction between the galaxies causes an inflow of gas towards the galactic centres and results in elevated and bursty star formation. The peak gas densities near the BHs increase by orders of magnitude resulting in enhanced accretion. Our results support the idea that galaxy mergers can trigger AGN activity, although the instantaneous accretion rate depends strongly on the local ISM. We also show that the level of merger-induced enhancement of BH fueling predicted by the SMUGGLE model is much smaller compared to the predictions by simulations using an effective equation of state model of the ISM.

A study of natural frequencies in a dynamic corona – disk system

Context. Black-hole X-ray binaries (BHXRBs) in the hard and hard-intermediate spectral (and temporal) states exhibit in their power spectra characteristic frequencies called type-C quasi-periodic oscillations (QPOs). Various models that can explain them with various degrees of success have been proposed, but a definitive answer is still missing. Aims. The hot Comptonizing corona interacting with the cold accretion disk, both of which are central in understanding BHXRBs, is essentially a dynamical system. Our aim is to investigate if the radiative coupling between the two components can produce QPOs. Methods. We write and solve the time-dependent equations that describe energy conservation in the system corona – accretion disk. We examine both constant and variable mass accretion rates. By necessity, in this first investigation we use a simple model, but it contains all the essential ingredients. Results. For a constant mass accretion rate and certain justifiable conditions, the dynamic corona – disk system exhibits oscillations, which die out after a few cycles. The characteristic frequencies of these oscillations are similar to the ones observed in the power spectra of BHXRBs. For most parameters, the natural frequencies persist even in the case of variable accretion rates. Conclusions. We argue that type-C QPOs in BHXRBs could, in principle, arise from the interaction of the hot Comptonizing corona with the much colder accretion disk. If this picture is correct, it has immediate implications for other systems that contain the above constituents, such as active galactic nuclei.