Revealing a multi-zone circumnuclear dust response in UGC 11487 through phase-resolved WISE W1–W2 diagnostics

We analysed epoch-averaged WISE/NEOWISE $W1$ and $W2$ photometry of UGC 11487 from 2014 to 2024. To constrain the pre-flare MIR level, we also queried the AllWISE multi-epoch photometry table (allwise_p3as_mep), which provides two 2010 $W1/W2$ epoch-averaged measurements. Together with the first pre-rise NEOWISE epochs, these data define the quiescent MIR baseline adopted for the baseline-subtracted excess light curves.

Single-exposure NEOWISE measurements were retained when both $W1$ and $W2$ had valid photometry, $S/N(W1)>5$, $S/N(W2)>3$, $\texttt{cc\_flags}=00$, and $\texttt{qual\_frame}>0$. The cleaned measurements were grouped according to the natural NEOWISE cadence. The retained NEOWISE measurements have stable PSF-fit diagnostics, $\texttt{w1fitr}=\texttt{w2fitr}=7.5\arcsec$, and no saturation flags, $\texttt{w1sat}=\texttt{w2sat}=0$.

The fiducial pre-flare baseline was defined as the median flux density of the two 2010 AllWISE epochs and the first two pre-rise NEOWISE epochs, yielding $F_{W1,{\rm base}}=8.453\penalty 10000\ {\rm mJy}$ and $F_{W2,{\rm base}}=5.428\penalty 10000\ {\rm mJy}$. The adopted baseline combines the earliest AllWISE and pre-rise NEOWISE measurements to reduce sensitivity to possible low-level variability and to provide a robust estimate of the quiescent MIR emission. The 2010 points were used only to define the baseline and were excluded from the flare fluence and energy integration.

WISE Vega magnitudes were converted to flux densities as

using $F_{0,W1}=309.54\penalty 10000\ {\rm Jy}$ and $F_{0,W2}=171.787\penalty 10000\ {\rm Jy}$. The baseline-subtracted excess flux density was defined as

All subsequent colour, timing, fluence, and colour-blackbody diagnostics were based on these baseline-subtracted excess light curves.

For model fitting and Monte Carlo resampling, we adopted an effective flux-density uncertainty

where $\sigma_{F_{\nu}}$ is the formal WISE flux-density uncertainty and $\sigma_{\rm floor}=1.10\penalty 10000\ {\rm mJy}$. This value accounts for residual systematic effects and prevents unrealistically small formal errors from dominating the fits. The uncertainty of the adopted pre-flare baseline was included separately in the Monte Carlo procedure by bootstrap resampling of the baseline epochs and perturbing their fluxes within their formal photometric uncertainties.

Additional robustness tests are provided in Appendix E.

To characterise the internal phase-dependent MIR evolution, we used empirical $W1-W2$ diagnostics based on colour evolution, colour–magnitude behaviour, flux–flux relations, and model-independent timing measures. This set of diagnostics traces both the spectral and temporal evolution of the MIR flare and quantifies the degree of phase-dependent complexity in the WISE response.

Observed-frame epochs were converted to rest-frame time relative to the first NEOWISE epoch,

The WISE colour index $W1-W2$ and the colour–magnitude relation were used to follow the MIR spectral evolution. The colour–magnitude relation was fitted using the York regression (York et al., 2004), accounting for correlated uncertainties between $W1-W2$ and $W1$. Flux–flux fits were performed with orthogonal distance regression for both absolute and baseline-subtracted flux densities.

The rise and decline occupy different loci in the flux–flux plane. For this reason, we fitted these phases separately to quantify potential hysteresis and deviations from a single-component evolution. We applied both a broad rise–decline split and a stricter division into the clean rise, turnover, and post-peak decline phases (Fig. 4).

In the absence of a well-sampled optical or UV driver, classical reverberation mapping cannot be applied. We therefore used internal MIR timing diagnostics based solely on the relative evolution of $W1$ and $W2$. The $W2$–$W1$ peak-time offset was defined as

and the positive-fluence centroid offset was defined as

where

The corresponding positive MIR fluence in each WISE band was

The integrals were evaluated numerically from the epoch-averaged light curves using trapezoidal integration. The 2010 AllWISE points were used only to define and resample the pre-flare baseline and were excluded from the flare fluence and energy integrations.

The positive-fluence centroid was adopted as the preferred timing measure because it incorporates the full positive-excess light curve. This reduces the dependence on a single sampled maximum and makes the estimate more robust against the sparse WISE cadence.

The corresponding timing offsets were expressed as $c\Delta t_{\rm rest}$ to place the internal $W1$–$W2$ timing difference on an approximate light-travel-distance scale. This quantity represents a characteristic separation between the $W1$- and $W2$-weighted dust responses. It is distinct from an absolute dust radius or an optical–MIR reverberation radius. The sensitivity of the adopted timing diagnostics to the baseline definition is tested in Appendix E.

We also used a shift-and-broaden $W1$-to-$W2$ template test to characterise the shape difference between the two MIR bands. In this test, the $W2$ excess light curve was represented as an amplitude-scaled, time-shifted, and temporally smoothed version of the $W1$ excess light curve,

where $A$ is an amplitude scale factor, $\tau$ is the compact $W1$-to-$W2$ shift, $G_{s}$ is a Gaussian smoothing kernel with smoothing timescale $s$, and $\ast$ denotes convolution. The fitted shift, $\tau$, describes the compact $W1$-to-$W2$ displacement after allowing for additional smoothing of the $W2$ response, while the smoothing timescale captures the broader temporal width of $W2$. This diagnostic separates the compact $W1$-to-$W2$ displacement from the additional temporal broadening of the $W2$ response. The fit was performed on a grid of $\tau$ and $s$, with the optimal amplitude obtained by weighted least squares for each grid point.

Uncertainties for the empirical diagnostics in this subsection were estimated through Monte Carlo resampling of the WISE fluxes, including bootstrap resampling of the pre-flare baseline.

Together, the colour, flux–flux, timing, and shift-and-broaden diagnostics probe different aspects of the MIR evolution. They are complemented by effective colour-blackbody quantities, which provide an empirical description of the temperature and emitting-scale evolution of the MIR flare.

We estimated effective two-band colour-blackbody quantities from the baseline-subtracted $W1$ and $W2$ excess fluxes, $\Delta F_{W}=F_{W}-F_{W,{\rm base}}$. These quantities provide empirical diagnostics of the temperature, emitting scale, and luminosity evolution of the MIR flare. They should be interpreted as descriptive thermal parameters of the WISE excess emission.

For epochs with $\Delta F_{W1}>0$ and $\Delta F_{W2}>0$, the effective colour temperature was obtained from the ratio of the $W1$ and $W2$ excess flux densities,

where $\nu_{\rm rest}=(1+z)\nu_{\rm obs}$. After $T_{\rm eff}$ was determined, the effective blackbody radius was estimated from the $W1$ excess flux density,

The approximate bolometric IR luminosity was then calculated as

and the integrated effective IR energy as

with the integration performed over rest-frame time. Uncertainties were evaluated through Monte Carlo resampling of both the WISE fluxes and the pre-flare baseline.

An epoch was classified as reliable if it was not used in the baseline definition, both $W1$ and $W2$ excess fluxes satisfied $\Delta F/\sigma_{\rm eff}\geq 2$, and the resulting Monte Carlo distributions were sufficiently constrained. Specifically, we required fractional uncertainties below 0.35 for $T_{\rm eff}$ and below 0.75 for both $R_{\rm BB}$ and $L_{\rm IR}$. Finite colour-blackbody solutions that did not satisfy these criteria were retained separately as poorly constrained estimates.

The derived blackbody radii represent effective emitting scales. They depend on the two-band colour temperature, the adopted baseline, and the assumption of a blackbody-like MIR excess. For this reason, they provide empirical constraints on the thermal evolution of the flare, while the underlying dust geometry remains degenerate.

The effective IR energy was also used for an energetic consistency check. Combining the effective IR energy with the published kinetic-energy scale of the radio-emitting outflow, we estimated the accreted-mass scale as

where $E_{\rm K,radio}$ is the published kinetic-energy scale of the radio-emitting outflow, $f_{\rm cov}$ is the effective dust reprocessing factor, and $\eta_{\rm eff}$ is an illustrative radiative efficiency. The corresponding disrupted-stellar mass scale was estimated as

where $f_{\rm acc}$ is the fraction of the disrupted stellar mass that is ultimately accreted. These quantities provide order-of-magnitude energetic consistency estimates and depend on the adopted values of $f_{\rm cov}$, $\eta_{\rm eff}$, and $f_{\rm acc}$.

Together, these effective quantities provide an independent empirical view of the thermal evolution of the MIR-emitting material and serve as a complementary test of the structured dust-response scenario explored below.

We employed two complementary levels of modelling. First, phenomenological multi-component models were used to quantify the level of smooth temporal complexity required by the WISE light curves. Second, physically motivated hydro-clumpy dust models driven by hydrodynamical TDE fallback histories were used as physical consistency tests for the empirical diagnostics.

The phenomenological decomposition describes the morphology of the MIR response in a compact parametric form. It provides a way to compare one-, two-, and three-component smooth representations of the $W1$ and $W2$ excess light curves.

We fitted the $W1$ and $W2$ excess light curves simultaneously with one-, two-, and three-component smooth phenomenological models. Each component was represented by a smooth, asymmetric temporal profile, with the amplitudes of $W1$ and $W2$ linked by the temperature-dependent Planck ratio. The component definitions are given in Appendix A.

The models were compared using the Akaike information criterion (AIC) and the Bayesian information criterion (BIC) (Akaike, 1974; Schwarz, 1978). These criteria assess whether additional smooth components improve the description of the two-band WISE data after penalising the increased number of free parameters. The decomposition is used to describe the MIR light-curve morphology and to provide a compact way to quantify the phase-dependent $W1-W2$ response.

As a physical consistency test, we also constructed a clumpy dust-echo forward model driven by hydrodynamical TDE fallback curves. The goal of these calculations is to examine whether the empirical WISE diagnostics are compatible with a physically plausible TDE-driven dust response.

Because only two WISE bands are available, substantial degeneracies between the input luminosity history, accretion smoothing, and dust geometry are expected. The input mass-supply histories were taken from the STARS fallback-rate library (Law-Smith et al., 2020a, b), broadened by an effective accretion/circularisation smoothing timescale, and used to illuminate a clumpy, radially stratified dust distribution. The model follows the standard dust-echo interpretation of TDE-powered infrared flares (e.g. Jiang et al., 2016; Lu et al., 2016; van Velzen et al., 2016, 2021) and uses modified-blackbody dust emission (e.g. Draine, 2003).

The hydro-clumpy models were selected using a combined objective function that includes both the $W1+W2$ light-curve residuals and the model-independent WISE diagnostics. The objective function also incorporates quantities such as the $W2/W1$ fluence ratio and centroid timing, which helps connect the forward models to the empirical MIR diagnostics. The detailed objective-function definition and implementation details are given in Appendix B. Additional toy $W1$-to-$W2$ response tests and the Monte Carlo parameters of the preferred two-component phenomenological model are given in Appendix D.

These hydro-clumpy calculations provide physical consistency tests for a structured multi-zone response. The resulting models represent plausible dust-response scenarios that reproduce the observed $W1-W2$ evolution within the adopted assumptions.

Thus, the phenomenological and hydro-clumpy approaches serve different purposes. The phenomenological models quantify the temporal complexity required by the WISE data, while the hydro-clumpy models assess whether that complexity can be reproduced within physically plausible TDE-driven dust scenarios.

We tested the stability of the $W2-W1$ centroid scale against three baseline definitions: pre-rise NEOWISE only, the fiducial AllWISE+NEOWISE baseline, and AllWISE 2010 only. We also tested the dependence of the smooth component comparison on the adopted effective error floor.

The possible late-time $W2$ residual around epochs 17–19, marked in the numbered WISE light curve in Fig. 1, was assessed by comparing those epochs with a smooth post-peak decline fitted after excluding them. This test evaluates the significance of possible deviations from a smooth decline and treats the late-time $W2$ excess as a residual feature whose interpretation remains model-dependent. Details of these tests are given in Appendix E.

To characterise the host-galaxy environment and assess the possible contribution of a persistent AGN component, we constructed an aperture-matched UV-to-radio spectral energy distribution. This analysis provides global host-galaxy properties and constrains the time-averaged non-transient emission.

We used archival data for UGC 11487 from GALEX, Pan-STARRS1, 2MASS, WISE, IRAS, and NVSS. For the photometric measurements, we matched full-galaxy apertures where angular resolution allowed, with separate low-resolution apertures for WISE W3–W4, IRAS, and NVSS. The resulting broadband fluxes were modelled with CIGALE to estimate the time-averaged stellar mass, star formation rate, attenuation, and dust luminosity. Particular attention was paid to the AGN fraction, since the presence of a dominant persistent AGN would complicate the interpretation of the MIR flare. The detailed aperture definitions, masking procedure, photometric conversions, treatment of IRAS upper limits, and validation images are given in Appendix F.

Thus, the host-galaxy SED modelling provides the environmental context for the MIR flare and constrains the contribution of any persistent AGN component, complementing the phase-resolved diagnostics presented above.

The epoch-averaged $W1$ and $W2$ light curves of UGC 11487 show a broad, asymmetric MIR flare lasting nearly a decade (Fig. 1). The peak-to-peak magnitude amplitude measured from the epoch-averaged light curves is larger in $W2$, $A_{W2}\simeq 2.27$ mag, than in $W1$, $A_{W1}\simeq 1.39$ mag. The $W2$ excess is also more persistent: integrating the baseline-subtracted positive excess over the flare gives a $W2$-to-$W1$ positive-fluence ratio of $2.23\pm 0.12$ (Table 2). These results show chromatic MIR variability with a substantial contribution from cooler dust traced most strongly by $W2$.

The peak of the excess emission occurs first in $W1$ and later in $W2$: the $W1$ excess reaches its maximum at epoch 6, whereas the $W2$ excess peaks at epoch 7 (Fig. 1). The $W2$ light curve also shows a slower decline and remains above the baseline after the $W1$ excess has substantially faded. This morphology already suggests that the MIR evolution is phase dependent and motivates the internal $W1$–$W2$ timing analysis presented below.

The WISE $W1-W2$ versus $W1$ colour–magnitude relation for UGC 11487 shows a redder-when-brighter (RWB) trend (Fig. 3). The York fit gives a slope of $-0.42\pm 0.16$, with a Pearson correlation coefficient of approximately $-0.63$. This trend is consistent with a stronger relative contribution from cooler, $W2$-emitting dust during the brighter MIR phases.

The phase dependence is more apparent in the WISE flux–flux trajectory (Fig. 4). During the rise, the source follows a lower branch in the $F_{W2}$–$F_{W1}$ plane. During the decline, it returns along an upper branch, with larger $F_{W2}$ at a similar $F_{W1}$. This hysteresis-like trajectory demonstrates that the MIR spectral evolution depends on phase and suggests a delayed and more persistent $W2$-emitting component.

The phase-resolved flux–flux fits quantify this behaviour. For the broad rise–decline division, the rising branch, defined by epochs 0–6, has $m=1.51\pm 0.09$, whereas the declining branch, defined by epochs 7–21, has $m=2.01\pm 0.18$. A stricter split excluding the turnover region, epochs 5–8, gives $m_{\rm rise}=1.43\pm 0.12$ for epochs 0–4 and $m_{\rm decline}=2.34\pm 0.06$ for epochs 9–21. The declining branch is therefore substantially more $W2$-enhanced than the rising branch, quantitatively confirming the phase dependence of the MIR evolution. Together, these empirical diagnostics show that the MIR evolution of UGC 11487 is phase dependent and departs from a self-similar single-component trajectory, motivating the timing analysis presented in Sect. 4.2.

The $W2$ emission is delayed relative to $W1$, as independently indicated by both peak-based and positive-fluence centroid diagnostics. The $W2$ excess peak occurs $191.8$ d after the $W1$ excess peak, corresponding to $c\Delta t_{\rm peak}=0.161$ pc (Table 3). This peak-based estimate is sensitive to the discrete NEOWISE sampling and is therefore used only as a supporting timing diagnostic.

The positive-fluence centroid provides the preferred timing measure because it uses the full positive-excess light curve and is less dependent on a single sampled maximum. For the combined AllWISE 2010 plus pre-rise NEOWISE baseline, the Monte Carlo resampling with bootstrap resampling of the baseline candidates gives $\Delta t_{\rm cent}=224\pm 45$ d, corresponding to $c\Delta t_{\rm cent}=0.188\pm 0.038$ pc (Table 3). We adopt this value as the characteristic $W2$–$W1$ MIR timing scale. The difference between the Monte Carlo estimate and the corresponding fixed-baseline value reflects the bootstrap resampling of the pre-flare baseline candidates included in the uncertainty analysis.

The $W2$–$W1$ centroid offset persists across alternative baseline choices. In fixed-baseline recalculations, using only the first two pre-rise NEOWISE epochs gives $\Delta t_{\rm cent}\simeq 207$ d, corresponding to $c\Delta t_{\rm cent}\simeq 0.173$ pc; using the combined AllWISE 2010 plus NEOWISE epochs 0–1 baseline gives $\Delta t_{\rm cent}\simeq 234$ d, corresponding to $c\Delta t_{\rm cent}\simeq 0.196$ pc; and using only the two 2010 AllWISE epochs gives $\Delta t_{\rm cent}\simeq 257$ d, corresponding to $c\Delta t_{\rm cent}\simeq 0.215$ pc (Appendix E). Although the exact value depends slightly on the adopted baseline, all tested definitions consistently indicate a delayed $W2$ response on sub-parsec scales.

The $W2$ excess is also broader in time than $W1$. The positive-excess temporal widths are approximately $697$ d for $W1$ and $781$ d for $W2$, with a width ratio of $\simeq 1.12$. Together with the larger $W2$ fluence, this demonstrates that the $W2$ response is both broader and more energetic than that of $W1$.

As an independent model-dependent characterisation of the light-curve morphology, we fitted a shift-and-broaden $W1$-to-$W2$ template model, in which the $W2$ excess is represented as a shifted and temporally smoothed version of the $W1$ excess. This gives a compact, effective shift of $\tau\simeq 110\pm 20$ d, corresponding to $c\tau\simeq 0.092\pm 0.017$ pc. The existence of two characteristic scales reflects that the $W2$ evolution combines a compact delayed component with substantial temporal broadening. The centroid scale traces the broader $W2$-weighted dust response, whereas $\tau$ characterises the more compact shift component of the $W1$-to-$W2$ evolution. Taken together, these timing diagnostics show that the $W2$ evolution is broader than a simple shifted replica of $W1$, indicating an extended MIR response.

The effective $W1$–$W2$ colour-blackbody analysis provides an independent empirical view of the thermal evolution of the MIR flare. Figure 5 shows the evolution of the effective colour temperature, blackbody radius, and approximate bolometric IR luminosity derived from the baseline-subtracted WISE excess fluxes.

Finite $T_{\rm eff}$ solutions are obtained for 20 non-baseline epochs. Of these, 12 epochs satisfy the adopted reliability cuts on the excess-flux signal-to-noise ratio and on the fractional uncertainties in $T_{\rm eff}$, $R_{\rm BB}$, and $L_{\rm IR}$. The remaining eight finite estimates are retained as poorly constrained, mainly because the $W1$ excess approaches the baseline at late times while $W2$ remains detectable. The first two NEOWISE epochs are used as part of the baseline definition and are therefore not assigned reliable colour-blackbody quantities.

The effective colour temperature decreases during the event. During the rise, $T_{\rm eff}$ is of order $900$ K. Near the MIR maximum it is typically $\sim 700$–$850$ K, and during the decline it decreases to $\sim 500$–$600$ K. Over the same period, the effective blackbody radius increases from $\sim 0.03$–$0.05$ pc on the rise to $\sim 0.1$–$0.15$ pc near and after the MIR maximum (Fig. 5). This evolution suggests a transition from a hotter compact MIR component to a cooler and more persistent $W2$-dominated response. At late times, the uncertainties in $R_{\rm BB}$ and $L_{\rm IR}$ increase because the $W1$ excess approaches the baseline while $W2$ remains detectable. These late-time estimates are therefore retained as poorly constrained effective emitting-radius estimates and are not used to infer precise dust radii.

The approximate bolometric IR luminosity reaches $\sim 3\times 10^{43}$ erg s^-1. Integrating the effective $L_{\rm IR}(t)$ over rest-frame time gives $E_{\rm IR}\simeq 2.9\times 10^{51}$ erg for the fiducial reliable-epoch estimate. Including all finite non-baseline estimates gives $E_{\rm IR}\simeq 3.8\times 10^{51}$ erg. We therefore regard the two-band effective IR energy as being of order a few $10^{51}$ erg.

The combined decrease in effective temperature, increase in emitting scale, and persistence of the $W2$-dominated emission independently support a thermally stratified and temporally extended MIR response.

The smooth joint $W1+W2$ phenomenological modelling shows that the observed WISE morphology cannot be adequately reproduced by a single smooth component. With the effective-error treatment, the one-component model gives $\chi^{2}_{\rm eff}=321.70$ and is strongly disfavoured relative to the two-component model, with $\Delta{\rm BIC}=281.7$ (Table 5). The smooth two-component model provides the preferred compact phenomenological description, with $\chi^{2}_{\rm eff}=21.06$, ${\rm BIC}_{\rm eff}=58.91$, and $\Delta{\rm BIC}=0$.

The three-component model reduces $\chi^{2}_{\rm eff}$ to 14.80 and follows the marginal late $W2$ shoulder more closely. However, after penalising the additional free parameters, it is disfavoured relative to the two-component model by $\Delta{\rm BIC}=12.7$. We therefore adopt the smooth two-component model as the preferred phenomenological description of the WISE morphology, while treating the three-component case as a diagnostic of possible late-time residual structure. The Monte Carlo parameter estimates for the preferred smooth two-component model are listed in Table 10.

In this phenomenological description, the first component is hotter and contributes more strongly to the early $W1$ emission, while the second component is cooler and longer-lived, contributing mainly to the $W2$-bright decline. This decomposition captures the main phase-dependent structure of the MIR flare and is consistent with a radially or thermally stratified dust response.

We also applied the clumpy dust-echo model driven by hydrodynamical TDE fallback curves described in Sect. 3.4 and Appendix B. The aim was to test whether the empirical MIR diagnostics of UGC 11487 can be reproduced by a smoothed TDE-like luminosity input reprocessed by a structured dust distribution. The adopted luminosity normalisation gives an effective bolometric input energy of $E_{\rm input}=2.44\times 10^{52}$ erg. For a fiducial radiative efficiency of $\eta_{\rm rad}=0.1$, this corresponds to an accreted-mass scale of order $0.1\,M_{\odot}$, consistent with the MIR energy-budget estimate discussed in Sect. 4.3.

The STARS fallback curves were not selected a priori to match the WISE light curves. We explored hydrodynamical fallback cases with $\beta=2.0$ and a grid of effective accretion-smoothing parameters. The representative solutions in Table 6 were selected using the combined objective function described in Appendix B. These representative solutions demonstrate that fallback-driven luminosity inputs can reproduce the empirical WISE diagnostics after accretion smoothing and dust reprocessing.

Within the explored parameter grid, representative three-zone solutions reproduce the combined empirical diagnostics more consistently than the corresponding two-zone cases. Representative three-zone solutions reproduce the main empirical diagnostics: the $W2$–$W1$ centroid offset, the $W2/W1$ fluence ratio, the $W2/W1$ temporal-width ratio, and the $W2$-enhanced declining branch. The detailed fallback curve, accretion-smoothing timescale, and dust-zone parameters remain degenerate, as expected for sparsely sampled two-band MIR data.

Case A gives the lowest objective value and reproduces the $W2$–$W1$ centroid offset and fluence ratio well, but it requires the highest peak luminosity, $L_{\rm peak}/L_{\rm Edd}\simeq 0.84$, and underestimates the global flux–flux slope relative to the empirical value. Case B is adopted as the fiducial compromise model. It has a lower peak luminosity, $L_{\rm peak}/L_{\rm Edd}\simeq 0.56$, while preserving the main timing, fluence, and width diagnostics: $\Delta t_{\rm cent}^{\rm model}(W2-W1)=216.4$ d, $\Phi_{W2}/\Phi_{W1}=2.218$, and $\sigma_{t,W2}/\sigma_{t,W1}=1.106$. Case C provides a still lower-Eddington alternative, with $L_{\rm peak}/L_{\rm Edd}\simeq 0.42$ and a somewhat larger global flux–flux slope, but with a higher objective value and a colder outer response component.

The fiducial Case B model is illustrated in Figs. 6 and 7. The schematic clumpy dust distribution and corresponding delay-response functions are shown in Fig. 6. The same model is compared with the WISE excess light curves in Fig. 7. It reproduces the delayed and broader $W2$ response, the larger $W2$ fluence, and the $W2$-enhanced decline. The qualitative conclusion is stable across the selected hydrodynamical drivers: the WISE data are consistently reproduced by a structured multi-zone response, with an early component contributing to $W1$ and a broader, cooler component producing the delayed and more persistent $W2$ emission.

We therefore use these hydro-clumpy solutions as representative physical consistency models. Within the constraints of the available two-band WISE data, they represent plausible structured dust-response scenarios for the circumnuclear environment. The different roles of the phenomenological decomposition and hydro-clumpy modelling are discussed in Appendix B.

We tested whether the apparent late-time $W2$ residual around epochs 17–19 in Fig. 1 represents a significant secondary excess. We fitted the surrounding post-peak $W2$ decline with a smooth exponential-plus-constant profile after excluding epochs 17–19, and then measured the residuals of the excluded epochs relative to this smooth decline. The details of this local residual test are given in Appendix E.

The $W2$ residual over epochs 17–19 is positive at the $\sim 2.1\sigma$ level, while the corresponding $W1$ residual is not significant, at $\sim 0.9\sigma$. The largest positive $W2$ residuals occur at epochs 17 and 18, whereas epoch 19 is consistent with the smooth decline within the uncertainties. We therefore identify this feature as a marginal late-time $W2$ shoulder within the long MIR decline. Its $W2$-dominated character is consistent with a cooler and more persistent MIR component, although the statistical significance remains modest.

We also compared the MIR-derived scales with the radio properties reported by Golay et al. (2025), which independently probe the compact nuclear environment of WTP14adeqka. The VLBA-resolved radio size is $0.111$–$0.127$ pc, while the VLA SED radio radii are $0.057$–$0.062$ pc. The preferred model-independent $W2$–$W1$ centroid scale, $c\Delta t_{\rm cent}=0.188\pm 0.038$ pc, traces the broader $W2$-weighted dust response, whereas the shift-and-broaden $W1$-to-$W2$ analysis gives a shorter compact-shift scale, $c\tau=0.092\pm 0.017$ pc. The latter lies above the VLA SED radio radii and close to the VLBA-resolved size.

This scale ordering supports a structured sub-parsec nuclear environment in which WISE traces thermal dust reprocessing over a radially or thermally extended dust distribution, while the radio emission traces synchrotron-emitting outflow or shock material on comparable compact nuclear scales. For a launch epoch near the onset of the positive MIR excess, the velocity required to reach the compact MIR shift scale by the $W2$-shoulder epoch is $\beta_{\rm req}\sim 0.040$, comparable in magnitude to the published radio expansion velocities for WTP14adeqka (Table 13). The relative placement of the characteristic MIR and radio scales is shown in Fig. 8. Their possible physical connection is discussed in Sect. 5.5.

The broadband UV–radio SED of UGC 11487 was modelled with CIGALE using the aperture-matched photometry described in Sect. 3.6 and in detail in Appendix F. The initial grid parameters and the CIGALE modules used for SED modelling are described in Appendix G. The fit reproduces the observed host-galaxy SED with a reduced $\chi^{2}_{\nu}\simeq 2.27$ (Fig. 9; Table 7). The Bayesian estimates indicate a moderately star-forming host galaxy, with $M_{\star}=(2.01\pm 0.15)\times 10^{10}\,M_{\odot}$, ${\rm SFR}=2.92\pm 0.15\,M_{\odot}\,{\rm yr}^{-1}$, and ${\rm sSFR}\simeq 1.45\times 10^{-10}\,{\rm yr}^{-1}$.

The fitted SED is consistent with moderate stellar-continuum attenuation and stronger attenuation towards the nebular-emitting regions: $E(B-V)_{\star}=0.176\pm 0.001$, $E(B-V)_{\rm lines}=0.400\pm 0.003$, and $E(B-V)_{\star}/E(B-V)_{\rm lines}\simeq 0.44$. The corresponding FUV attenuation is $A_{\rm FUV}=2.17\pm 0.02$ mag. The infrared emission is dominated by star-formation-heated dust, with $L_{\rm dust}=(1.05\pm 0.05)\times 10^{37}\,{\rm W}$, $\log(L_{\rm dust}/L_{\odot})\simeq 10.44$, and $\alpha_{\rm dust}=1.67\pm 0.08$, placing UGC 11487 below the classical LIRG regime. The low-significance IRAS 12, 25, and 60 $\mu{\rm m}$ measurements were treated as upper limits, and the 100 $\mu{\rm m}$ point as a marginal far-infrared constraint, as described in Appendix F.

Overall, the host-galaxy properties are consistent with a moderately star-forming isolated barred disc in which the MIR flare appears as a transient phenomenon superposed on an otherwise normal galaxy.

UGC 11487/WTP14adeqka is a long-lived MIR-selected nuclear flare whose evolution is dominated by dust-reprocessed emission. The decade-scale WISE light curve shows a larger peak-to-peak amplitude and a larger time-integrated excess in $W2$ than in $W1$, together with an RWB trend and phase-dependent flux–flux behaviour (Table 2). These properties are consistent with a transient heating event followed by a structured dust response, as expected when a substantial fraction of the primary UV/optical radiation from a nuclear transient is absorbed and re-emitted by circumnuclear dust (Jiang et al., 2016, 2021; Masterson et al., 2024; Necker et al., 2025).

The MIR evolution is consistent with a dust-obscured nuclear transient and, in particular, with a TDE-like accretion episode reprocessed by circumnuclear dust. A definitive classification as a classical TDE is limited by the absence of a directly and densely sampled primary optical/UV flare. Alternative nuclear-accretion scenarios, including an accretion-state change or a changing-look AGN episode, therefore remain viable. We consequently interpret UGC 11487 as an obscured TDE candidate or TDE-like nuclear flare.

The host-galaxy SED modelling presented in Sect. 4.6 provides an important context for interpreting the MIR flare. The time-averaged UV–far-infrared emission of UGC 11487 is consistent with a moderately star-forming barred disc whose broadband energy budget is dominated by stellar light and star-formation-heated dust. In the adopted CIGALE model, the persistent broadband AGN contribution is negligible, with a best-fit value consistent with zero.

This result separates the time-averaged host-galaxy emission from the transient MIR component. The decade-long WISE flare, strong $W2$ excess, colour hysteresis, and delayed $W2$ response cannot be explained as static features of the host-galaxy SED and require an additional time-dependent nuclear component. UGC 11487 therefore appears to be a moderately star-forming isolated barred disc that underwent a temporary nuclear dust-reprocessing episode, compatible with an obscured TDE-like or other short-lived nuclear accretion event.

TDEs in the Milky Way would be associated with Sgr A*, although no Galactic TDE has been directly observed to date. Evidence for episodic past activity of Sgr A* has been obtained from X-ray observations of molecular clouds in the Galactic Centre region. In particular, studies by Sunyaev and collaborators demonstrated that several dense clouds, including Sgr B2, act as natural X-ray reflection nebulae, reprocessing radiation emitted during past high-luminosity flares of the Galactic SMBH (Revnivtsev et al., 2004; Churazov et al., 2017; Khabibullin et al., 2022; Marin et al., 2023). Variability, spectral properties, and polarisation measurements of the reflected emission indicate that Sgr A* experienced significantly higher activity in the recent past than in its present quiescent state.

In general, the expected TDE rate is low, of the order of $\Gamma_{\rm TDE}\sim 10^{-5}-10^{-4}\,{\rm yr}^{-1}$, corresponding to approximately one event per $10^{4}$–$10^{5}$ yr (Komossa, 2015; Alexander, 2005). Nevertheless, episodic TDE-driven energy injection has been discussed as a possible contributor to the past activity of Sgr A* and to large-scale Galactic-centre structures such as the Fermi and eROSITA bubbles (Sazonov and Sunyaev, 2012; Bland-Hawthorn et al., 2013; Su et al., 2010; Predehl et al., 2020).

In this context, TDEs may represent short-lived episodes in the evolution of otherwise normal spiral galaxies. Systems such as UGC 11487 are therefore of particular interest, since they demonstrate that quiescent Milky-Way-like galaxies can occasionally undergo energetic accretion events associated with stellar disruption and temporarily exhibit enhanced nuclear activity. Such episodes may provide a glimpse of processes that could have occurred in the recent history of the Milky Way itself. Consequently, TDE candidates among isolated barred galaxies offer a useful bridge between studies of transient accretion phenomena and evolutionary diagnostics of MWNeTs (Vavilova et al., 2026)

Several independent diagnostics consistently indicate that a single self-similar MIR response is inadequate for UGC 11487. The $W2$ emission peaks later than $W1$, declines more slowly, and has a larger positive fluence. The preferred model-independent timing diagnostic is the positive-fluence centroid offset, $\Delta t_{\rm cent}=224\pm 45$ d, corresponding to $c\Delta t_{\rm cent}=0.188\pm 0.038$ pc. This scale traces the broader $W2$-weighted dust response. The shift-and-broaden $W1$-to-$W2$ analysis gives a shorter compact-shift scale, $\tau\simeq 110\pm 20$ d, corresponding to $c\tau\simeq 0.092\pm 0.017$ pc, showing that the $W2$ evolution combines a delayed response with substantial temporal broadening.

The colour–magnitude and flux–flux behaviour support the same picture. The RWB colour–magnitude trend is accompanied by a $W2$-enhanced declining branch in the flux–flux plane. This hysteresis-like behaviour suggests that cooler dust continues to emit after the hotter $W1$-dominated component has faded. The effective colour-blackbody diagnostics show a coherent cooling trend, from $\sim 900$ K during the rise to $\sim 500$–$600$ K during the decline, while the effective emitting radius increases to sub-parsec scales. These two-band quantities trace the temperature, emitting scale, and luminosity evolution of the MIR-emitting dust.

The smooth phenomenological modelling also favours a structured response. A single smooth $W1+W2$ component is strongly disfavoured, while the preferred two-component phenomenological description provides the most compact representation once the adopted effective-error treatment and BIC penalty are included. The late $W2$ shoulder remains a marginal residual within the long MIR decline, with insufficient significance to define a separate flare. The robust observational result is therefore a delayed, broader, and more persistent $W2$ response, consistent with a multi-zone circumnuclear dust structure.

The hydrodynamically driven clumpy dust-echo models provide a physical consistency test linking these empirical WISE diagnostics to a TDE-like accretion scenario. In this framework, the main MIR properties of UGC 11487 can be reproduced when a hydrodynamical fallback input is smoothed by an effective accretion/circularisation timescale and reprocessed by a multi-zone circumnuclear dust distribution. The observable MIR emission is expected to differ from the prompt fallback rate because the event can be reshaped by circularisation, shocks in incompletely circularised debris, outflows, reprocessing layers, viewing geometry, light-travel-time delays, dust temperature stratification, and the radial distribution of the reprocessing material (Roth et al., 2020b). With only $W1$ and $W2$ photometry, the detailed fallback history, smoothing timescale, and dust geometry remain degenerate.

The representative hydro-clumpy models require an accreted mass of order $0.1\,M_{\odot}$, consistent with the MIR energy-budget estimate discussed in Sect. 4.3 and with partial accretion of bound debris in a TDE-like event. Across the tested models, $W2$ remains delayed, broader, and more persistent than $W1$, and the declining branch remains $W2$-enhanced. The hydro-clumpy modelling therefore supports the interpretation of UGC 11487 as a TDE-like accretion transient reprocessed by structured circumnuclear dust.

The strength of this interpretation comes from the combination of diagnostics: timing, colour evolution, flux–flux hysteresis, effective cooling, phenomenological model comparison, and hydro-clumpy consistency tests all point to a structured multi-zone circumnuclear response.

Independent radio observations of WTP14adeqka provide an external probe of the compact nuclear environment (Golay et al., 2025). Together with the WISE data, they allow a comparison between thermal dust reprocessing and synchrotron-emitting plasma associated with jets, outflows, unbound debris, or shocks in TDEs (Alexander et al., 2020; Roth et al., 2020b). The MIR and radio observations therefore provide complementary diagnostics of the circumnuclear response to the same nuclear transient.

Delayed radio emission is increasingly recognised as an important stage of TDE evolution. Radio emission emerging hundreds to thousands of days after discovery has been reported in both optically and X-ray-selected TDEs, favouring delayed outflows in many systems (Cendes et al., 2024; Goodwin et al., 2025). Long-lived radio emission may also contain multiple temporal components and be affected by host contamination (Christy et al., 2026), while years-delayed radio flares in slowly evolving TDE candidates provide a relevant comparison for WTP14adeqka (Zhang et al., 2026).

We compared the MIR-derived timing and response scales with the independent radio constraints reported by Golay et al. (2025), using the radio measurements only as external reference quantities. The compact MIR shift scale from the shift-and-broaden analysis, $c\tau=0.092\pm 0.017$ pc, is comparable to the VLBA characteristic size of approximately $0.11$–$0.13$ pc reported by Golay et al. (2025). The preferred $W2$–$W1$ centroid scale, $c\Delta t_{\rm cent}=0.188\pm 0.038$ pc, is larger and traces the broader $W2$-weighted dust response. Since the MIR and radio quantities are derived from different observables, their agreement should be interpreted as an order-of-magnitude consistency of characteristic sub-parsec scales.

The marginal late-time $W2$ shoulder occurs during the radio-bright phase. Assuming a launch epoch near the onset of the positive MIR excess, the velocity required to reach the compact MIR shift scale by the $W2$-shoulder epoch is $\beta_{\rm req}\sim 0.040$, comparable to the VLBA size-increase velocity reported for WTP14adeqka (Table 13). Although this travel-time estimate depends on the assumed launch epoch, the temporal ordering and characteristic scales are compatible with an interaction between a radio-emitting outflow or shock and compact dusty circumnuclear material.

Several physical scenarios can account for delayed radio activity. In the debris–torus collision model of Lei et al. (2024), unbound TDE debris impacts a dusty torus and produces a delayed radio outburst. Alternatively, Linial et al. (2026) proposed that viscous expansion of the TDE accretion disc can lead to repeated collisions with a pre-existing stellar EMRI, launching a slow and energetic outflow after an intrinsic delay. They identified WTP14adeqka as a plausible example of this mechanism. Within the structured dust environment inferred from the WISE evolution, either promptly produced debris or an intrinsically delayed outflow could contribute to late-time dust heating or reprocessing while propagating through compact circumnuclear material.

The late-time $W2$ shoulder remains only marginally significant, and the MIR and radio measurements trace physically distinct emitting components. The present observations therefore establish temporal and order-of-magnitude scale consistency, whereas the launching mechanism of the radio outflow and its causal connection to the late-time MIR evolution remain model-dependent.

The IR and radio energetics provide an additional consistency check for the TDE-like interpretation. Using the effective IR energy derived from the WISE light curves together with the published kinetic-energy scale of the radio-emitting outflow, we obtain $M_{\rm acc}\simeq 0.04{-}0.22\,M_{\odot}$ for $E_{\rm IR}\simeq(2.9{-}3.8)\times 10^{51}$ erg, $E_{\rm K,radio}\sim 10^{50.7}$ erg, $\eta_{\rm eff}=0.1$, and $f_{\rm cov}=0.1{-}0.5$. For an accreted fraction $f_{\rm acc}=0.1{-}0.5$, this corresponds to an illustrative disrupted-star mass scale of $M_{\star,\rm scale}\simeq 0.07{-}2.2\,M_{\odot}$. These order-of-magnitude estimates indicate that the energetics inferred independently from the infrared flare and the radio-emitting outflow are mutually compatible with a TDE-like event involving a star with a sub-solar to a few-solar-mass scale.

The main limitations follow directly from the available data. WISE provides only two MIR bands, so $T_{\rm eff}$, $R_{\rm BB}$, $L_{\rm IR}$, and $E_{\rm IR}$ are effective two-band colour-blackbody diagnostics of the MIR evolution. The NEOWISE cadence is sparse, making the positive-fluence centroid a more stable timing measure than the discrete $W1$ and $W2$ peak epochs. The pre-flare baseline also introduces a systematic uncertainty, which motivated the baseline-sensitivity tests of the $W2$–$W1$ centroid scale.

More detailed dust-echo and reprocessing modelling can provide stronger geometric constraints when the driving optical/UV flare is well sampled. For example, convex-ring or torus-remnant models applied to AT 2019qiz (Wu et al., 2025), and multi-component disk-reprocessing and MIR dust-echo modelling of AT 2020nov (Earl et al., 2025), demonstrate how broader multiwavelength coverage can separate the primary emission, reprocessing structure, and dust response. In UGC 11487, the primary optical/UV driver was not observed, and the MIR coverage is limited to $W1$ and $W2$. The preferred two-component phenomenological $W1+W2$ description therefore constrains the presence of a structured MIR response, while the detailed dust geometry remains degenerate.

The radio measurements provide an important independent consistency check. They were not incorporated into the fitting procedure and therefore cannot remove the intrinsic degeneracies of the infrared modelling. In addition, the radio sizes reported by Golay et al. (2025) are VLBA-derived characteristic Gaussian scales, whereas the MIR quantities are light-travel response scales inferred from WISE timing. Their comparison should therefore be treated as a characteristic-scale consistency test.

The main observational conclusions are robust within these limitations: $W2$ peaks later than $W1$, declines more slowly, has a larger positive fluence, shows phase-dependent flux–flux behaviour, and traces a cooling effective dust component.

The present data leave the detailed three-dimensional dust structure degenerate. Even with this limitation, phase-resolved $W1$–$W2$ diagnostics recover characteristic sub-parsec response scales and provide valuable constraints on the internal structure of dust reprocessing in obscured accretion-powered transients.