M. Guolo, A. Mummery, A. Ingram, M. Nicholl, S. Gezari, E. Nathan

We present the implementation of a fully time-dependent relativistic disk model - based on the light curve fitting package FitTeD - into the X-ray spectral fitting environment, pyXspec. This implementation enables simultaneous fitting of multi-epoch and multi-wavelength spectral data, where the only free parameters are those describing the black hole and the initial conditions, while the subsequent evolution is fully governed by the dynamical equations of an evolving accretion flow. We apply the model to fit seven epochs of X-ray spectra and two epochs of UV spectra of the 'long-lived' tidal disruption event (TDE) and quasi-periodic eruption (QPE) source GSN 069, spanning from its discovery in 2010 through late-2019. Our results show that such 'long-lived', X-ray-bright TDEs - of which GSN 069 is a prime, but not unique, example-can be naturally explained within the same framework as events with shorter-lived X-ray emission, such as ASASSN-14li and AT2019dsg. Their key distinction lies in the 'viscous' timescale parameter - tied to the disk's stress and angular momentum transport efficiency - which should be treated as a free parameter when modeling the disk evolution of transient events. We examine the implications for QPE models by tracking the time evolution of disk properties such as mass surface density and accretion rate. We show that existing QPE models cannot reproduce the observed connection between the presence (e.g., 2018) or absence (e.g., 2014) of eruptions and disk properties. In the context of orbiter-disk collision models, the change in mass surface density appears insufficient to explain the needed variation in the eruption's temperature. The absence of eruptions in GSN 069 in 2014 remains a major challenge for QPE models.