We present a new method of modelling time-series data based on the running optimal average (ROA). By identifying the effective number of parameters for the ROA model, in terms of the shape and width of its window function and the times and accuracies of the data, we enable a Bayesian analysis, optimising the ROA width, along with other model parameters, by minimising the Bayesian Information Criterion (BIC) and sampling joint posterior parameter distributions using MCMC methods. For analysis of quasar lightcurves, our implementation of ROA modelling can measure time delays among lightcurves at different wavelengths or from different images of a lensed quasar and, in future work, be used to inter-calibrate lightcurve data from different telescopes and estimate the shape and thus the power-density spectrum of the lightcurve. Our noise model implements a robust treatment of outliers and error-bar adjustments to account for additional variance or poorly-quantified uncertainties. Tests with simulated data validate the parameter uncertainty estimates. We compare ROA delay measurements with results from cross-correlation and from JAVELIN, which models lightcurves with a prior on the power-density spectrum. We analyse published COSMOGRAIL lightcurves of multi-lensed quasar lightcurves and present the resulting measurements of the inter-image time delays and detection of microlensing effects.

Context. A complex environment exists in the inner few astronomical units of planet-forming disks. High-angular-resolution observations play a key role in our understanding of the disk structure and the dynamical processes at work. Aims. In this study we aim to characterize the mid-infrared brightness distribution of the inner disk of the young intermediate-mass star HD 163296 from early VLTI/MATISSE observations taken in the L- and N-bands. We put special emphasis on the detection of potential disk asymmetries. Methods. We use simple geometric models to fit the interferometric visibilities and closure phases. Our models include a smoothed ring, a flat disk with an inner cavity, and a 2D Gaussian. The models can account for disk inclination and for azimuthal asymmetries as well. We also perform numerical hydrodynamical simulations of the inner edge of the disk. Results. Our modeling reveals a significant brightness asymmetry in the L-band disk emission. The brightness maximum of the asymmetry is located at the NW part of the disk image, nearly at the position angle of the semimajor axis. The surface brightness ratio in the azimuthal variation is 3.5 ± 0.2. Comparing our result on the location of the asymmetry with other interferometric measurements, we confirm that the morphology of the r < 0.3 au disk region is time-variable. We propose that this asymmetric structure, located in or near the inner rim of the dusty disk, orbits the star. To find the physical origin of the asymmetry, we tested a hypothesis where a vortex is created by Rossby wave instability, and we find that a unique large-scale vortex may be compatible with our data. The half-light radius of the L-band-emitting region is 0.33 ±0.01 au, the inclination is 52°−7°+5°, and the position angle is 143° ± 3°. Our models predict that a non-negligible fraction of the L-band disk emission originates inside the dust sublimation radius for μm-sized grains. Refractory grains or large (≳10 μm-sized) grains could be the origin of this emission. N-band observations may also support a lack of small silicate grains in the innermost disk (r ≲ 0.6 au), in agreement with our findings from L-band data.

Polycyclic aromatic hydrocarbons (PAHs) are carbon-based molecules that are ubiquitous in a variety of astrophysical objects and environments. In this work we use JWST/MIRI MRS spectroscopy of three Seyferts to compare their nuclear PAH emission with that of star-forming (SF) regions. This study represents the first of its kind to use sub-arcsecond angular resolution data of local luminous Seyferts (Lbol >  1044.46 erg s−1) with a wide wavelength coverage (4.9–28.1 μm). We present an analysis of their nuclear PAH properties by comparing the observed ratios with PAH diagnostic model grids derived from theoretical spectra. Our results show that a suite of PAH features is present in the innermost parts of luminous Seyfert galaxies (∼0.45″ at 12 μm; in the inner ∼142–245 pc). We find that the nuclear regions of active galactic nuclei (AGN) lie at different positions of the PAH diagnostic diagrams, whereas the SF regions are concentrated around the average values of SF galaxies. In particular, we find that the nuclear PAH emission mainly originates in neutral PAHs. In contrast, PAH emission originating in the SF regions favours ionised PAH grains. The observed PAH ratios in the nuclear region of the AGN-dominated galaxy NGC 6552 indicate the presence of larger PAH molecules compared with those of the SF regions. Therefore, our results provide evidence that the AGN have a significant impact on the ionisation state (and probably the size) of the PAH grains on scales of ∼142–245 pc.

Compact Obscured Nuclei (CONs) potentially hide extreme supermassive black hole (SMBH) growth behind large column densities of gas/dust. We present a spectroscopic analysis of the heavily obscured nucleus and the surrounding environment of the eastern region of the nearby (z = 0.02007) interacting galaxy VV 114 with the JWST Mid-InfraRed Instrument (MIRI). We model the spectrum from 4.9 - 28 μm to extract Polycyclic Aromatic Hydrocarbon (PAH) emission and the underlying obscured continuum. We find that the NE nucleus (A) is highly obscured where the low PAH equivalent width (EW) ratio, EW(12.7)/EW(11.3), reveals a dust enshrouded continuum source. This is confirmed by decomposing the continuum into nuclear and star-forming where the nuclear component is found to be typical of CONs. The 11.3/6.2 PAH flux ratio is consistent with originating in star-forming regions rather than typical AGN. The second nucleus (B) is much less obscured, with PAH flux ratios also typical of star-forming regions. We do not detect any high ionisation lines such as [Ne V] or [Ne VI] which suggests that if an AGN is present it must be highly obscured. Additionally, we detect a shock front south of the secondary nucleus (B) in the [Fe II] (5.34 μm) line and in warm molecular hydrogen. The 6.2 PAH emission does not spatially coincide with the low-J transitions of H2 but rather appears strong at the shock front which may suggest destruction of the ionised PAHs in the post-shock gas behind the shock front.

We measure the black hole mass and investigate the accretion flow around the local (z = 0.0502) quasar PG 1119+120. Spectroscopic monitoring with Calar Alto provides Hβ lags and linewidths from which we estimate a black hole mass of log (M•/M⊙) = 7.0, uncertain by ∼0.4 dex. High cadence photometric monitoring over two years with the Las Cumbres Observatory provides lightcurves in 7 optical bands suitable for intensive continuum reverberation mapping. We identify variability on two timescales. Slower variations on a 100-day timescale exhibit excess flux and increased lag in the u′ band and are thus attributable to diffuse bound-free continuum emission from the broad line region. Faster variations that we attribute to accretion disc reprocessing lack a u′-band excess and have flux and delay spectra consistent with either τ∝λ4/3, as expected for a temperature structure of T(R)∝R−3/4 for a thin accretion disc, or τ∝λ2 expected for a slim disc. Decomposing the flux into variable (disc) and constant (host galaxy) components, we find the disc SED to be flatter than expected with $f_{\nu } \sim \rm {const}$. Modelling the SED predicts an Eddington ratio of λEdd &gt; 1, where the flat spectrum can be reproduced by a slim disc with little dust extinction or a thin disc which requires more dust extinction. While this accretion is super-Eddington, the geometry is still unclear, however a slim disc is expected due to the high radiation pressure at these accretion rates, and is entirely consistent with our observations.