Can a time-evolving, asymmetric broad line region mimic a massive black hole binary?

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.