Since the discovery of active galactic nuclei (AGN), their light curves have proved to be intrinsically variable. In the optical/UV band, this variability is consistent with correlated or red noise and is particularly well described by the damped random walk (DRW) model. In this work, we evaluate the feasibility of a new method for identifying spatially unresolved couples of active galactic nuclei through a fully Bayesian time-domain analysis of the observed light curves. More specifically, we check whether observed light curves are better described by a single DRW, which we interpret as emitted by a single massive black hole, or a pair of independent damped random walks, generated by a pair of massive black holes. We test the method on mock light curves associated with single massive black holes and pairs generated with different cadences and lengths of observational campaigns. We constrained the occurrence of false positives, that is, the percentage of single MBH light curves that show substantial evidence in favour of the unresolved MBH pair scenario, finding a fraction of 0.2% and 0.59% in the even and uneven sampling scenarios. We discuss how well the method recovers the model parameters, showing that about 51% and 7% of the simulated light curves have all the recovered parameters within 20% of their true values in our best scenario of evenly sampled light curves for the single massive black hole and massive black hole pair scenarios, respectively. We finally study the region of the parameter space in which the detection of a massive black hole pair is possible, finding that such objects can be correctly identified if the timescales of the process describing the noise are very different, with a ratio smaller than similar to 0.2, and the variability amplitudes are similar, with their ratio bigger than similar to 0.2. When limiting to such a region of the parameter space, the fraction of pairs with all the recovered parameters within 20% of the injected values increases up to about 14% and 8% for evenly and unevenly sampled light curves, respectively.

Gas within the influence sphere of accreting massive black holes is responsible for the emission of the broad lines observed in optical-UV spectra of unobscured active galactic nuclei. Since the region contributing the most to the broad emission lines (i.e. the broad line region) depends on the active galactic nucleus luminosity, the study of broad line reverberation to a varying continuum can map the morphology and kinematics of gas at sub-pc scales. In this study, we modify a preexisting model for disc-like broad line regions, including non-axisymmetric structures, by adopting an emissivity profile that mimics the observed luminosity-radius relation. This makes our implementation particularly well suited for the analysis of multi-epoch spectroscopic campaigns. After validating the model, we use it to check if strongly non-axisymmetric, single broad line regions could mimic the short time-scale evolution expected from massive black hole binaries. We explore different orientations and anisotropy degrees of the broad line region, as well as different light curve patterns of the continuum to which the broad line region responds. Our analysis confirms that recently proposed algorithms designed to search for massive black hole binaries in large multi-epoch spectroscopic data are not contaminated by false positives ascribed to anisotropic broad line regions around single massive black holes.

Variability studies of the broad emission lines of active galactic nuclei (AGN) and quasars show stochastic radial velocity variations (i.e., fluctuations in the centroid of the line), “jitter,” on timescales of weeks to months. This jitter may be intrinsic as the broad-line region (BLR) reverberates from the AGN continuum. There are also coordinated variations in the width of the broad emission lines and the luminosity of the central source (“breathing” or “antibreathing”) that remain unexplained. These can be used as a tool for testing models of the BLR. We have constructed a pipeline to simulate a disk-like BLR geometry that reverberates in response to various chosen continuum light curves and produce synthetic emission line profiles. These profiles can then be characterized by measured shape parameters (centroid velocity shift, velocity dispersion, and Pearson skewness coefficient) and compared to observed time series of those same parameters. We have found that through our pipeline, we can recreate the velocity jitter at similar variations found in observations. The computational tools presented in this paper will also be applicable to case studies of quasars observed under the Sloan Digital Sky Survey V Black Hole Mapper reverberation mapping program. This paper is the first in a series of papers—in this paper, we present the model and pipeline, and in future papers, we will present applications.

Massive black hole binaries are expected to be observable as periodic active galactic nuclei in time-domain photometric surveys. Periodicities may originate from different physical processes, including the intermittent gas feeding of the black holes caused by the time-varying non-axisymmetric binary potential, the Doppler boosting of the flux emitted by individual accretion discs bound to the orbiting black holes, and the gravitational lensing of the accretion disc of one black hole by another such disc. Only the Doppler boost scenario applied to circular binaries with non-modulated accretion predicts a sinusoidal light curve, while in the general case, binary signals are expected to show more complex periodic patterns. Current searches for massive black hole binaries rely on techniques tailored to quasi-sinusoidal light curves, but fail to identify the more complex predicted periodicities. We present an alternative method that leverages Gaussian processes, making use of a generic periodic kernel that is flexible enough to identify arbitrary periodicities in unevenly sampled light curves with realistic quasar noise. We demonstrate that it outperforms previously proposed strategies in identifying general periodicities by analysing mock light curves with different baselines. Specifically, our analysis can detect non-sinusoidal periodicities (e.g. sawtooth-shaped or symmetric flares) and retrieves a higher fraction of true periodicities when compared to a periodogram analysis or a Gaussian process analysis with less flexible periodic kernels. Furthermore, by comparing the retrieved fraction of periodicities between mock Palomar transient factory light curves and mock Legacy Survey of Space and Time light curves, we found that our analysis is most sensitive to the number of observed cycles. The application of this analysis has the potential to greatly increase the scientific return of current and upcoming large time-domain photometric surveys.

The interaction between massive black hole (MBH) binaries and circumbinary discs (CBDs) affects the spin orientations of merging binaries - a key observable for the future LISA mission. While gas accretion is known to align BH spins with the orbital angular momentum via the Bardeen-Petterson effect, the impact of AGN feedback on this process has remained largely unexplored. We present hydrodynamical simulations of an equal-mass MBH binary embedded in a CBD, using the GIZMO code with a subgrid model that self-consistently evolves BH mass, spin, and anisotropic AGN feedback. We explore different BH spin magnitudes and orientations to identify the configurations producing the strongest feedback effect on the disc structure and binary evolution. AGN winds substantially modify the CBD and minidisc structure, enlarging the central cavity and, after 20-25 binary orbits, destroying the minidiscs in all feedback simulations. Feedback-driven cavity excavation suppresses the gravitational torques driving orbital evolution, causing the binary to stall. AGN feedback dramatically affects MBH accretion, reducing the Eddington ratio by one to two orders of magnitude before ultimately suppressing it entirely. This strongly inhibits Bardeen-Petterson spin alignment, delaying it well beyond the feedback-free timescale. Extending one simulation to 70 binary orbits reveals a feedback-regulated duty cycle of alternating active and quiescent phases, with the binary spending most of its time in low-density, quiescent conditions that further suppress alignment. These results indicate that AGN feedback is an important, previously overlooked channel for preserving spin misalignment in MBH binaries prior to coalescence, with direct implications for interpreting LISA gravitational-wave observations.

The Laser Interferometer Space Antenna is expected to observe gravitational waves from massive black hole binaries across cosmic time. Many are anticipated to be detectable hours to weeks before coalescence. We present a fast algorithm for the premerger detection and preliminary characterization of such binaries. The method performs a search for excess power with a chirping time-frequency morphology in short-time Fourier transform spectrograms. By tiling the time-frequency plane with slices defined by the quadrupole frequency evolution, we define a signal significance relative to a fitted background distribution of instrumental noise and Galactic foreground. Individual search triggers are followed by a coherence tracker that groups triggers consistent with the same physical signal over time. Doing so, our analysis provides progressively refined estimates of the chirp mass and coalescence time. We validate our algorithm on the SangriaHM LISA Data Challenge dataset, successfully detecting all 15 injected massive black hole binaries; 14 of them hours to weeks before merger, while one is only detected after the binary coalescence. The algorithm yields chirp mass relative errors below 3% for high-SNR sources and coalescence time uncertainties of up to a few hours. With a computational cost of less than a second to process a 10-day data segment on a single core, our approach is suitable for generating real-time alerts, triggering protected observational periods, and providing informative priors for Bayesian parameter estimation.

Since their discovery, AGN light curves are known to be intrinsically variable. In the optical/UV band, this variability is consistent with correlated or red noise and is particularly well described by the damped random walk (DRW) model. In this work, we evaluate the feasibility of a new method for identifying spatially unresolved couples of AGN through a fully Bayesian time-domain analysis of the observed light curves (LCs). More specifically, we check whether observed LCs are better described by a single DRW, which we interpret as emitted by a single massive black hole (MBH), or a pair of independent DRWs, generated by a pair of MBHs. We test the method on mock LCs associated with a single MBH and pairs generated with different cadences and lengths of observational campaigns. We constrained the occurrence of false positives, that is, the percentage of single MBH LCs that show substantial evidence in favour of the unresolved MBH pair scenario, finding a fraction of 0.2

We present a new, Bayesian analysis of the highest-resolution optical spectrum of the supermassive black hole (SMBH) binary candidate PG 1302-102, obtained with ESPRESSO@VLT (R \simeq 138, 000). Our methodology, based on robust Bayesian model selection, reveals the presence of multiple narrow emission lines at the expected redshift of the source and confirms (for H\{beta}) and detects (for H{\gamma}) the presence of redshifted broad components. Additionally, we have discovered a very broad and, if it is associated with the H\{beta}, very redshifted component at {\lambda} \simeq 5000{\AA}. We evaluate two scenarios for explaining the observed broad emission line (BEL) in PG 1302-102. In the case in which the redshifted BEL asymmetry arises from the orbital motion of a putative binary, our measurements coupled with simple estimates of the broad-line region (BLR) sizes suggest that the individual black hole BLRs are either settled in a single BLR or in the process of merging and, therefore truncated and highly disturbed. Alternatively, in the scenario of a single SMBH, we explain the distorted emission of the BELs with a nonsymmetric distribution of the BLR clouds; namely, a thin disk with a spiral perturbation. This BLR configuration is statistically preferred over any empirical multi-Gaussian fit and simultaneously explains the asymmetric emission of the H{\beta} and H{\gamma} close to the bulk of the line and any additional excess (or the lack of it, in the case of the H{\gamma}) at much longer wavelengths. The physical origins of the perturbation are unclear, and a connection with the possible presence of a black hole binary cannot be ruled out. Given the growing evidence from theoretical and observational works demonstrating the common presence of disturbed BLRs in active galactic nuclei, we argue that an origin related to self-gravitating instabilities may be more plausible.

Supermassive black hole binaries (SBHBs) are a natural outcome of galaxy mergers, and they are expected to be among the loudest gravitational-wave sources at low frequencies. The source SDSS J2320+0024 was recently proposed as a promising SBHB candidate due to a possible periodicity in its light curve and variability in the MgII emission line. In this work, we reanalysed the optical (gr, and i bands) light curves of J2320+0024 within the framework of Bayesian model selection. When periodicity was searched for together with red noise, analysis of the g-band light curve reveals a peak in the posterior of the period at similar to 290 days. The posterior profile is too broad to yield a preference for periodic models over models that include only red noise. Furthermore, the same peak is not present in the analysis of the r-band and i-band light curve. A periodic model without red noise identified a different (similar to 1100 days) periodicity, but this model is statistically significantly disfavoured relative to the other models tested. In summary, we find no significant evidence in favour of a true periodic signal over red-noise variability. Our analysis questions the robustness of the previously proposed periodicity and emphasises the importance of rigorous statistical treatment. While our findings challenge the binary interpretation for J2320+0024, they do not rule it out. A statistically robust joint analysis of the photometric light curves and evolving broad-line profiles would shed further light on the true nature of this object.

Massive black hole binaries are expected to be observable as periodic AGN in time-domain photometric surveys. Periodicities may originate from different physical processes, including the intermittent gas feeding of the black holes caused by the time-varying non-axisymmetric binary potential, the Doppler boosting of the flux emitted by individual accretion discs bound to the orbiting BHs, and the gravitational lensing of the accretion disc of one black hole due to the presence of the other. Only the Doppler boost scenario applied to circular binaries with non-modulated accretion predicts a sinusoidal light curve, while in the general case, binary signals are expected to show more complex periodic patterns. Current searches for massive black hole binaries rely on techniques tailored to quasi-sinusoidal light curves, but fail to identify the more complex periodicities predicted. We present an alternative method that leverages Gaussian processes, making use of a generic periodic kernel flexible enough for the identification of arbitrary periodicities in unevenly sampled light curves with realistic quasar noise. We demonstrate that it outperforms previously proposed strategies in identifying general periodicities by analysing mock light curves with different baselines. Specifically, we find that our analysis can detect non-sinusoidal periodicities (e.g., sawtooth-shaped) and retrieves a higher fraction of true periodicities when compared to periodogram analysis or Gaussian processes analysis with less flexible periodic kernels. Furthermore, by comparing the retrieved fraction of periodicities between mock PTF light curves and mock LSST light curves, we find that our analysis is most sensitive to the number of observed cycles. The application of this analysis has the potential to greatly increase the scientific return of current and upcoming large time-domain photometric surveys.

We present a suite of zoom-in cosmological simulations of Milky Way-like galaxies with a prominent disc component and a strong bar in their centre based on a subsample of barred galaxies from the TNG50 magneto-hydrodynamic simulation. We modified the physical models that regulate star formation, namely, the supernova feedback and black hole quasar feedback, in order to examine how they affect disc and bar formation. We find that independent of the feedback prescriptions, all galaxies show a similar morphology, which is dominant in comparison with the bulge mass. The black hole quasar feedback models used in this study do not affect bar formation, although they can affect the bar strength and length. The energy released by the supernovae causes a delay in the time of bar formation, and in models with the strongest feedback, galaxies form stable discs against bar formation. This could be expected since supernova feedback influences disc and bulge assembly and results in discs with a lower mass content, radial velocity dispersion, and a larger size as the supernova feedback strength increases. We studied disc stability using three bar instability criteria proposed in the literature. We find that galaxies with varied supernovae and black hole quasar feedback satisfy these criteria at the moment of bar formation except in extreme cases where the galaxy lacks or has weak supernova feedback. In these models, two of the three criteria fail to forecast the existence (or absence) of a bar, probably because they do not account for the influence of a massive compact bulge. Our findings provide insights into the physical processes behind bar formation and highlight the importance of additional conditions other than a massive compact disc that promote bar formation.

Emission from two massive black holes (MBHs) bound in a close binary is expected to be modulated by different processes, such as the Doppler boost due to the orbital motion, accretion rate variability generated by the interaction with a circumbinary disc, and binary gravitational self-lensing. When the binary is compact enough, the two black holes are thought to be surrounded by a common broad-line region that reprocesses the impinging periodically varying ionising flux, creating broad emission lines with variable line shapes. Therefore, the study of broad emission line variability through multi-epoch spectroscopic campaigns is of paramount importance for the unambiguous identification of a binary. In this work, we study the response of a disc-like broad-line region to the Doppler-boosted ionising flux emitted by sub-milliparsec MBH binaries on a circular orbit and compare it with the response of a broad-line region illuminated by a single MBH with a periodically but isotropically varying intrinsic luminosity. We show that in the binary case, the time lags of the blue and red wings of the broad emission lines, arising from diametrically opposite sides of the circumbinary disc, are out of phase by half of the binary's orbital period, as they each respond to the periodic 'lighthouse' modulation from the binary's continuum emission. This asymmetric time lag represents a new binary signature that cannot be mimicked by a single MBH.

We study the impact on the cosmic microwave background (CMB) landscape of peculiar rotational general relativistic effects. These effects, on galactic scales, do not possess a Newtonian analogue, and therefore could a priori impact CMB analysis. We find that the velocity inferred from the CMB dipole, under the kinematic interpretation, coincides with that measured by a stationary observer within the Milky Way and not with the one measured by the zero angular momentum observer. We show that the galaxy peculiar frame-dragging effects do not impact the standard CMB analysis, as these modify the multipole coefficients only at higher orders with respect to the dominant terms. Moreover, we prove that no general relativistic framework at the galactic scale patched within the standard cosmological model can account for the current tension on the CMB quadrupole amplitude.

We present JWST/NIRSpec integral field spectroscopy in the rest-frame optical bands of the system PJ308-21, a quasar at z = 6.2342 caught as its host galaxy interacts with companion galaxies. We detect the spatially extended emission of several emission lines (H alpha, H beta, [O III], [N II], [S II], and He II), which we used to study the properties of the ionized phase of the interstellar medium: the source and hardness of the photoionizing radiation field, metallicity, dust reddening, electron density and temperature, and star formation. We also marginally detected continuum starlight emission associated with the companion sources. We find that at least two independent satellite galaxies are part of the system. While the quasar host appears highly enriched and obscured, with photoionization conditions typical of an Active Galactic Nucleus, the western companion shows minimal dust extinction, low metallicity (Z similar to 0.4 Z(circle dot)), and star formation driven photoionization. The eastern companion shows higher extinction and metallicity (Z similar to 0.8 Z(circle dot)) compared to the western companion, and it is at least partially photoionized by the nearby quasar. We do not find any indication of AGN in the companion sources. Our study shows that while the quasar host galaxy is already very massive (M-dyn > 10(11) M-circle dot), it is still rapidly building up by accreting two relatively massive (M-star similar to 10(10) M-circle dot) companion sources. This dataset showcases the power of JWST in exposing the buildup of massive galaxies in the first gigayear of the Universe.

We present JWST/NIRSpec integral field data of the quasar PJ308-21 at $z=6.2342$. As shown by previous ALMA and HST imaging, the quasar has two companion sources, interacting with the quasar host galaxy. The high-resolution G395H/290LP NIRSpec spectrum covers the $2.87-5.27\ \rm \mu m$ wavelength range and shows the rest-frame optical emission of the quasar with exquisite quality ($S/N\sim 100-400$ per spectral element). Based on the H$\beta$ line from the broad line region, we obtain an estimate of the black hole mass $M_{\rm BH,H\beta}\sim 2.7\times 10^{9}\ \rm M_{\odot}$. This value is within a factor $\lesssim 1.5$ of the H$\alpha$-based black hole mass from the same spectrum ($M_{\rm BH, H\alpha}\sim 1.93\times 10^{9}\ \rm M_{\odot}$) and is consistent with a previous estimate relying on the MgII $\lambda 2799$ ($M_{\rm BH, MgII}\sim 2.65\times 10^{9}\ \rm M_{\odot}$). All these $M_{\rm BH}$ are within the $\sim 0.5$ dex intrinsic scatter of the adopted mass calibrations. The high Eddington ratio of PJ308-21 $\lambda_{\rm Edd,H\beta}\sim 0.67$ ($\lambda_{\rm Edd,H\alpha}\sim 0.96$) is in line with the overall quasar population at $z \gtrsim 6$. The relative strengths of the [OIII], FeII and H$\beta$ lines are consistent with the empirical "Eigenvector 1" correlations as observed for low redshift quasars. We find evidence for blueshifted [OIII] $\lambda 5007$ emission with a velocity offset $\Delta v_{\rm [OIII]}=-1922\pm 39$ km s$^{-1}$ from the systemic velocity and a $\rm FWHM([OIII])=2776^{+75}_{-74}$ km s$^{-1}$. This may be the signature of an outflow from the nuclear region, despite the true values of $\Delta v_{\rm [OIII]}$ and $\rm FWHM([OIII])$ are likely more uncertain due to the blending with H$\beta$ and FeII lines. Our study demonstrates the unique capabilities of NIRSpec in capturing quasar spectra at cosmic dawn and studying their properties in unprecedented detail.

The compact object in quasar 3C 186 is one of the most promising recoiling black hole candidates, exhibiting both an astrometric displacement between the quasar and the host galaxy as well as a spectroscopic shift between broad and narrow lines. 3C 186 also presents a radio jet that, when projected onto the plane of the sky, appears to be perpendicular to the quasar-galaxy displacement. Assuming a gravitational-wave kick is indeed responsible for the properties of 3C 186 and using state-of-the-art relativistic modeling, we show that current observations allow for exquisite modeling of the recoiling black hole. Most notably, we find that the kick velocity and the black hole spin are almost collinear with the line of sight and the two former vectors appear perpendicular to each other only because of a strong projection effect. The targeted configuration requires substantial fine-tuning: while there is a region in the black hole binary parameter space that is compatible with 3C 186, the observed system appears to be a rare occurrence. Using archival radio observations, we explored different strategies that could potentially confirm or rule out our interpretation. In particular, we developed two observational tests that rely on the brightness ratio between the approaching and receding jet as well as the asymmetry of the jet lobes. While the available radio data provide loose constraints, deeper observations have the unique potential of unveiling the nature of 3C 186.

We investigate the relevance of kinematically identified bulges, discs and their role relative to galaxy quenching. We utilize an analysis of the SDSS-MaNGA survey conducted with the GPU-based code BANG which simultaneously models galaxy photometry and kinematics to decompose galaxies into their structural components. Below M~1011 Msun, galaxies exhibit a wide range of dynamical properties, determined by the relative prominence of a dispersion-supported inner region and a rotationally-supported disc. Our analysis reveals a natural separation between these classes, with only a minor fraction of stellar mass retained by structures exhibiting intermediate dynamical support. When examining galaxies in terms of their star formation activity, an apparent decrease in rotational support is observed as they move below the star-forming main sequence. This behaviour is evident with luminosity-weighted tracers of kinematics, while it almost vanishes with mass-weighted tracers. Luminosity-weighted quantities not only capture differences in kinematics but also in the stellar population, potentially leading to biased interpretations of galaxy dynamical properties and quenching. Our findings suggest that quenching does not imply almost any structural transformation in galaxies below M~10^11 Msun. Processes as disc fading more likely account for observed differences in mass-weighted and luminosity-weighted galaxy properties; when the galactic disc ceases star formation, its mass-to-light ratio grows without any significant morphological transformation. The picture is remarkably different above M~10^11 Msun. Regardless of the tracer used, a substantial increase in galaxy dispersion support is observed along with a significant structural change. A different quenching mechanism, most likely associated with mergers, dominates. Notably, this mechanism is confined to a very limited range of high masses.

Observations by the James Webb Space Telescope of the Universe at $z\gtrsim 4$ have shown that massive black holes (MBHs) appear extremely overmassive compared to the local correlation for active galactic nuclei. In some cases, these objects might even reach half the stellar mass inferred for the galaxy. Understanding how such objects formed and grew to this masses has then become a big challenge for theoretical models, with different ideas ranging from heavy seed to super-Eddington accretion phases. Here, we take a different approach, and try to infer how accurate these MBH mass estimates are and whether we really need to revise our physical models. By considering how the emerging spectrum (both the continuum and the broad lines) of an accreting MBH changes close to and above the Eddington limit, we infer a much larger uncertainty in the MBH mass estimates relative to that of local counterparts, up to an order of magnitude, and a potential preference for lower masses and higher accretion rates, which i) move them closer to the local correlations, and ii) might indicate that we are witnessing for the first time a widespread phase of very rapid accretion.

Two-body relaxation may drive stars onto near-radial orbits around a massive black hole, resulting in a tidal disruption event (TDE). In some circumstances, stars are unlikely to undergo a single terminal disruption, but rather to have a sequence of many grazing encounters with the black hole. It has long been unclear what is the physical outcome of this sequence: each of these encounters can only liberate a small amount of stellar mass, but may significantly alter the orbit of the star. We study the phenomenon of repeating partial tidal disruptions (pTDEs) by building a semi-analytical model that accounts for mass loss and tidal excitations. In the empty loss cone regime, where two-body relaxation is weak, we estimate the number of consecutive partial disruption that a star can undergo, on average, before being significantly affected by two-body encounters. We find that in this empty loss cone regime, a star will be destroyed in a sequence of weak pTDEs, possibly explaining the tension between the low observed TDE rate and its higher theoretical estimates.

According to the hierarchical formation paradigm, galaxies form through mergers of smaller entities and massive black holes (MBHs), if present at their centers, migrate to the nucleus of the newly formed galaxy, where they form binary systems. The formation and evolution of MBH binaries, and in particular their coalescence timescale, is highly relevant for current and future facilities aimed at detecting the gravitational wave signal produced by the MBHs close to coalescence. While most of the studies targeting this process are based on hydrodynamic simulations, the high computational cost makes a complete parameter space exploration prohibitive. Semianalytic approaches represent a valid alternative, but they require ad hoc prescriptions for the mass loss of the merging galaxies in minor mergers due to tidal stripping, which is not commonly considered or is at best modelled assuming very idealised geometries. In this work we propose a novel, effective model for the tidal stripping in axisymmetric potentials, to be implemented in semi-analytic models. We validated our semi-analytic approach against N-body simulations considering different galaxy sizes, inclinations, and eccentricities, finding only a moderate dependence on the orbit eccentricity. In particular, we find that, for almost circular orbits, our model mildly overestimates the mass loss, and this is due to the adjustment of the stellar distribution after the mass is removed. Nonetheless, the model exhibits a very good agreement with simulations in all the considered conditions, and thus represents an extremely powerful addition to semi-analytic calculations.