{"id":"7c9a65e3-8000-4caf-a992-35cb2f5bc05c","arxiv_id":"2412.17046","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":11,"one_line_summary":"SDSS J000118.70+003314.0, a z=1.4 quasar, shows a nine-year brightening and fading pattern that is modeled as a 1.9 solar-mass star tidally disrupted by a 6.5x10^7 solar-mass black hole.","lead":"Astronomers report a possible tidal disruption event: a star shredded by a supermassive black hole in a distant quasar at redshift 1.4, traced through nine years of SDSS light curves. If real, it is one of the most distant such events found inside a broad-line quasar and supports the idea that some quasar flickering is caused by shredded stars.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 0.009% DRW probability ignores the 7253-quasar search and is computed with CAR parameters fitted to the same light curve, overstating TDE confidence.","rationale":"The reader's verdict is CONDITIONAL, and the weakest assumption identified is the per-object DRW probability being used without a trials correction and with parameters fitted to the same light curve. I agree that this is the most load-bearing weakness. The paper's central quantitative support for the TDE interpretation is the 0.009% probability from the CAR-process simulation in Section 5.3. However, the simulation does not account for the fact that SDSS J0001 was selected because its light curve showed a rise-and-decline pattern among 7253 quasars visually inspected, nor for the fact that the CAR parameters are estimated from the full light curve including the purported flare. A proper null would either fix CAR parameters from a pre-flare segment or use ensemble quasar variability parameters from MacLeod et al. (2010), and would propagate the search trials. The reduced chi-square of 4.5 also means the TDE model is not a particularly good fit, so the acceptance criterion is lenient. These issues do not disprove the TDE candidate, but they strip the 0.009% number of its stated statistical force. The candidate remains plausible and worth follow-up, especially since it is one of very few high-redshift TDE candidates in a broad-line AGN and the light-curve shape, color evolution, and order-of-magnitude BH mass are broadly consistent with TDE expectations. The recommendation is therefore to retain the CONDITIONAL verdict: the paper should report a trials-corrected or clearly per-object probability, re-estimate the CAR null without the flare segment, and ideally release the light curves and fitting code for independent verification.","tokens_in":28150,"tokens_out":3258,"duration_ms":30148,"concrete_test":"Recompute the false-positive probability with a validated two-stage procedure. First, fit the CAR(1) model to the g-band light curve of SDSS J0001 using only the pre-flare baseline (e.g., MJD-50500 < 2500), separate from the decline phase, and generate 10^5 mock light curves with these uncontaminated parameters; count how many yield a TDE-model fit with χ²/dof < 4.5, giving an unbiased per-object p. Then repeat the same mock-generation and TDE-fit pipeline on the other 19 visually selected Stripe82 candidates (or a random sample from the 7253 quasars) to estimate the experiment-wide false-positive rate in the parent sample. If the expected number of false positives in 7253 trials is ≳0.1, or if the per-object p computed with pre-flare CAR parameters exceeds 0.1%, the claimed 0.009% significance is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing claim is the quoted probability p≈0.009% that SDSS J0001's long-term variability is intrinsic AGN variability (Abstract; Section 5.3). This p-value is not a valid experiment-wide significance because it is computed for one object selected from a visual search of 7253 z>1 Stripe82 quasars (Section 2). With 9/10^5 mock light curves passing the TDE-model criterion for this object, the per-object false-positive rate is ~9×10^-5; multiplying by 7253 independent trials gives ~0.65 expected false positives. Thus the reported value cannot exclude the hypothesis that the flare is a chance DRW fluctuation in the parent sample. The issue is compounded because the CAR(1) parameters τ≈470±260 d and σ are fit with JAVELIN to the same g-band light curve that contains the candidate flare, and the mock light curves are then generated from those data-derived parameters (Section 5.3). If the flare is real, the fitted CAR process is contaminated by the very signal being tested; if the flare is intrinsic DRW noise, fitting CAR to a light curve that by selection contains an extreme fluctuation biases the mock ensemble toward reproducing that fluctuation. Either way the simulation is not an independent null. A third, secondary, point is that the 'TDE-like' criterion is the reduced chi-square of the best TDE fit, χ²/dof≈4.5, which is a poor fit by conventional standards, so the 9 mock acceptances are relative to a very loose threshold. The combination of an uncorrected search trials factor and a data-fitted null means the 0.009% probability substantially overstates confidence in the TDE interpretation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"SDSS J000118.70+003314.0 (SDSS J0001) is a z=1.404 quasar with a broad Mg II line. It was one of 20 flare-shaped candidates found by visual inspection of the light curves of 7253 z>1 Stripe82 quasars, and is analyzed here because it has the highest redshift among them. The paper fits its nine-year ugriz photometry (about 60 epochs per band from Stripe82 plus PHOTOOBJALL) with the publicly available MOSFIT/TDEFIT tidal disruption event model and reports a best fit with a main-sequence star of about 1.9 M_sun and a black hole of 6.5x10^7 M_sun, at reduced chi^2/dof of about 4.5. The authors then test the null hypothesis of intrinsic AGN variability by fitting a damped random walk (CAR) model with JAVELIN to the g-band light curve (tau about 470 days), generating 10^5 mock light curves, and counting how many the TDE model fits as well as the data; 9 of 10^5 pass, giving the quoted 0.009% probability. Alternative explanations (dust extinction, microlensing, accretion) are discussed and argued to be disfavored. The paper concludes that SDSS J0001 is a high-redshift TDE candidate in an AGN, supporting the view that some AGN variability is produced by tidal disruption events.","tokens_in":28551,"tokens_out":24984,"duration_ms":209672,"significance":"If confirmed, this would be one of the very few TDE candidates hosted by a broad-line quasar at z>1 and would add empirical support to the idea that central stellar disruptions contribute to AGN variability. The manuscript has genuine strengths: the TDE modeling is not ad hoc, since it uses published MOSFIT/TDEFIT hydrodynamical templates and standard scaling relations; the flare characterization is multi-pronged (Gaussian-rise/exponential-decay fit, Weibull shape parameters, g-r color evolution); the DRW/CAR test is a serious attempt to quantify the intrinsic-variability null rather than asserting it; and each alternative scenario receives a concrete physical argument. The analysis is transparent and reproducible in structure: the light-curve sources, the SQL query, the fitting codes (kmpfit, emcee, JAVELIN), and the template interpolations (Appendix B) are all described. These strengths make the candidate worth keeping in the literature, provided the statistical significance issues described in the major comments are resolved.","major_comments":[{"comment":"The quoted probability p about 0.009% is a per-object false-positive rate and does not account for the way SDSS J0001 was selected. Section 2 reports a visual inspection of the light curves of 7253 z>1 Stripe82 quasars, from which 20 flare-shaped candidates were chosen, with SDSS J0001 analyzed as the highest-redshift example. Section 5.3 then simulates 10^5 mock DRW light curves for this single object and finds 9 that pass the TDE-fit criterion. If the per-object rate is about 9x10^-5, the parent sample of 7253 quasars yields about 0.65 expected false positives, and the authors' own alternative calculation with tau uniformly drawn from [50,5000] days (0.144%) yields about 10 expected false positives. The observation is therefore fully consistent with the null of intrinsic DRW variability once the search is taken into account. I recommend that the authors either analyze all 20 selected candidates with the same machinery, apply an explicit trial correction, or clearly state that 0.009% is a per-object conditional rate and refrain from presenting it in the Abstract as the probability that the event results from intrinsic variability of quasars.","section":"Section 5.3; Abstract"},{"comment":"The DRW/CAR null simulations are calibrated on the very light curve under test. The parameters tau about 470+260-178 days and sigma are fitted with JAVELIN to the same g-band light curve that is hypothesized to contain the TDE (Section 5.3, Figure 9), and the mock light curves are then generated from these data-derived parameters using the variance (0.07 mag^2) of that same light curve. If the flare is real, the fitted CAR process is contaminated by the signal; if the flare is an extreme DRW fluctuation, the object was selected from 7253 quasars precisely because of that fluctuation. In neither case is the simulation an independent null, and the resulting 0.009% is a conditional probability given the fitted parameters, not the probability that the event is intrinsic variability as the Abstract states. Depending on the direction of the bias, this could either over- or under-state the significance. I recommend re-running the simulation with population-level priors on (tau, sigma), or with parameters fitted only to the pre- and post-flare portions of the light curve, and reporting how the resulting probability changes.","section":"Section 5.3; Eq. (10)"},{"comment":"The fit quality does not support the Abstract's statement that the light curves can be described by the conventional TDE model. Section 4 reports reduced chi^2/dof of about 4.5; with roughly 300 photometric points (Section 2 states about 60 per band in five bands) and 13 free parameters (Table 1), this corresponds to chi^2 of about 1300 for about 287 degrees of freedom and is a formally poor fit. The MCMC parameter uncertainties in Table 1 (e.g., log(M_star/M_sun)=0.28 with very small asymmetric errors) are therefore underestimated, because no rescaling for chi^2/dof>1 is applied. Please quantify the likely sources of excess scatter (underestimated photometric errors, correlated systematics, model deficiencies), rescale the uncertainties accordingly, and justify the chi^2/dof<4.5 acceptance threshold used for the mock light curves in Section 5.3 in light of this assessment.","section":"Section 4; Table 1"}],"minor_comments":[{"comment":"The claim that SDSS J0001 is the highest-redshift optical TDE candidate known in AGN is internally inconsistent with the Introduction, which lists SDSS J120414.37+351800.5 at z about 2.359 as a TDE candidate in a N-rich quasar; please qualify the statement (e.g., among candidates with well-sampled multi-band light curves) or correct it.","section":"Section 2"},{"comment":"The spectrum was taken near the flare peak, so the continuum luminosity entering the Shen et al. (2011) virial Mg II estimator is likely flare-boosted; part of the factor-7.5 discrepancy between the virial and TDE black-hole masses may therefore arise from the virial estimator itself, a possibility not discussed alongside the proposed non-virial dynamics of TDE debris.","section":"Section 5.2"},{"comment":"The dust-extinction test in the left panel of Figure 11 derives E(B-V) from binned i-band mean magnitudes under the assumption that all i-band variability is extinction; the test would be biased if the i-band contains intrinsic variability, and the binning and interpolation choices are not described in detail.","section":"Section 5.4"},{"comment":"No fit statistics are reported for the 1S2L microlensing model; a chi^2 value or a residual plot is needed to support the claim that microlensing is disfavored by the achromaticity argument.","section":"Section 5.4; Figure 11"},{"comment":"Equation (10) is not written as a proper stochastic differential equation: the drift term is missing the dt factor and the noise term notation is ambiguous; please correct the expression so that the CAR process is unambiguously defined.","section":"Eq. (10)"},{"comment":"The comparison with the Weibull parameters of simulated lensing flares in Graham et al. (2017) is purely qualitative; a figure or table showing the location of SDSS J0001 relative to the simulated 1S1L and lensing distributions would make the argument against lensing quantitative.","section":"Section 2"}],"recommendation":"major_revision","confidential_remarks":"The paper's central quantitative claim, the 0.009% probability of intrinsic variability, is a per-object simulation result that ignores the 7253-object visual search and is calibrated on CAR parameters derived from the same light curve. These are load-bearing issues for the headline claim, but they are fixable: the authors could compute experiment-wide rates over the 20 selected candidates, apply a trial correction, or reframe the claim as a per-object conditional rate. The TDE fit itself is not circular, uses published templates, and is transparently documented, so I regard the manuscript as a candidate for major revision rather than rejection. A second issue to watch in revision is the fit quality (chi^2/dof about 4.5) and the correspondingly underestimated parameter uncertainties. If the authors present the result as a candidate with appropriately caveated statistics, the paper would be a useful contribution to the literature on TDEs in AGNs."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nQuick take: this is a candidate, not a discovery. The paper reports the highest-redshift optical TDE candidate in a broad-line quasar (z=1.404), and the headline evidence is a 0.009% DRW/CAR probability that the flare is intrinsic AGN variability. The object is new, and the nine-year Stripe82/PHOTOOBJALL light curves do show a rise-to-peak then decline that the MOSFIT TDE model can roughly reproduce. The authors also give more than a passing thought to dust and microlensing alternatives, which is better than many candidate papers. The 7.5x discrepancy between the TDE-derived BH mass and the virial Mg II mass is interesting, though the non-virial debris explanation is speculative.\n\nThe soft spot is exactly where the stress-test lands. Section 5.3 computes that 0.009% by fitting JAVELIN CAR parameters to the same g-band light curve that supposedly contains the TDE, then simulating 1e5 light curves from those contaminated parameters. Only 9 of the mocks pass the TDE-fit threshold (chi2/dof<4.5), giving a per-object false positive rate of 9e-5. But SDSS J0001 was visually selected from 7253 z>1 quasars, with 20 flare-shaped candidates flagged. Per-object rate times 7253 trials gives ~0.65 expected false positives, so the quoted probability is not the experiment-wide significance. It is also not a clean per-object probability because the null model is trained on the signal. To make matters worse, the reduced chi2 of 4.5 for the real light curve is poor—the 'TDE-like' threshold is loose. The alternative calculation with tau drawn uniformly from 50-5000 days gives 0.144%, which is still not trial-corrected. So the statistical case is overstated in the abstract, even though the candidate might be real.\n\nWho should read this? Anyone tracking TDE candidates in AGN or the high-z TDE population. It is a useful data point to know, but not something I would build on until the statistics are redone. The paper is readable and honest about alternatives, but the data and code are not released beyond the public databases, which limits reproducibility.\n\nRecommendation: send it to a serious referee, but with the clear expectation of major revision on the statistics. The object deserves follow-up; the 0.009% claim does not deserve to stand as is.","headline":"A plausible but unproven high-z TDE candidate in a quasar; the headline 0.009% probability is not meaningful as stated.","tokens_in":29127,"tokens_out":4711,"would_cite":false,"duration_ms":44835,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The nine-year $ugriz$ light curves of quasar SDSS J0001 can be modeled as a main-sequence star of about 1.9 solar masses being tidally disrupted by a black hole of about $6.5\\times10^7$ solar masses, making it the highest-redshift optical…","keywords":["tidal disruption event","broad-line quasar","high redshift","quasar variability","damped random walk","Stripe82 light curves","black hole mass","SDSS J000118.70+003314.0"],"falsifier":"A decisive observation would be continued photometric and spectroscopic monitoring of SDSS J0001 over the next several years. A tidal disruption event should fade smoothly toward the pre-flare level and not re-flare on timescales of years, while the damped random-walk model predicts continuing stochastic variability; the first clear re-brightening or erratic upturn would contradict the TDE interpretation. In addition, a spectrum taken after the flare has faded could check whether the Mg II line width and shift follow the ordinary quasar virial relation, which the debris scenario predicts they should not.","tokens_in":27942,"feed_emoji":"🔭","tokens_out":10422,"duration_ms":88231,"temperature":0.7,"pith_summary":"The paper reports a high-redshift tidal disruption event (TDE) candidate in the broad-line quasar SDSS J000118.70+003314.0 at redshift $z=1.404$. It argues that the nine-year SDSS $ugriz$ light curves, which show a rise to peak followed by a smooth decline, can be described by the conventional TDE model: a main-sequence star of about $1.9\\,M_\\odot$ disrupted by a black hole of about $6.5\\times10^7\\,M_\\odot$. It further estimates, through damped-random-walk simulations of quasar variability, that the chance the flare is intrinsic quasar churn is only about 0.009 per cent. If correct, SDSS J0001 would be the highest-redshift optical TDE candidate found in a broad-line quasar, supporting the idea that TDEs occur in active galaxies and that some AGN variability is powered by shredded stars.","feed_headline":"Quasar's nine-year flare matches a star being shredded","feed_subtitle":"At z=1.4 it would be the most distant such flare in a broad-line quasar; odds it is ordinary quasar flicker: 0.009%.","key_machinery":"The argument runs on the standard tidal-disruption fallback model in which a star is shredded at the tidal radius, the debris returns to the black hole on a spread of binding energies, and the resulting accretion rate is convolved with a viscous delay to produce a blackbody-photosphere light curve. The paper uses hydrodynamically computed fallback-rate templates, scales them to arbitrary black-hole and stellar masses through the known scaling relations, and adopts a main-sequence mass-radius relation, then fits the model to all five SDSS bands with MCMC. A second machine, the damped random walk (CAR) stochastic process, generates $10^5$ mock light curves with the quasar's measured variance and correlation time; counting how many of those mocks can also be fitted by the TDE model yields the claimed 0.009 per cent probability that the flare is intrinsic AGN variability.","core_discovery":"On the paper's own terms, the long-term variability of SDSS J0001 is a tidal disruption event rather than a coincidence of AGN activity. The observed $ugriz$ light curves are well fitted by the standard viscous-delayed fallback model, yielding a disrupted main-sequence star of mass $1.905^{+0.023}_{-0.009}\\,M_\\odot$ and a central black hole of mass $6.5^{+3.5}_{-2.6}\\times10^7\\,M_\\odot$, with about $0.78\\,M_\\odot$ accreted and $1.12\\,M_\\odot$ ejected. The flare's peak luminosity and photosphere temperature are moderate compared with known optical TDEs, and its unusually long observed timescale follows from a large impact parameter, a relatively massive black hole, and redshift time dilation. The virial black-hole mass from the broad Mg II line is about 7.5 times larger than the TDE fit, which the paper attributes to non-virial motions of TDE debris contributing to the broad-line emission. Finally, mock light curves generated with the damped random walk (CAR) process and fitted with the TDE model give a probability of about 0.009 per cent that the observed flare is intrinsic quasar variability.","pith_inferences":["Applying the paper's TDE-template fit to the other 19 visually selected flare candidates from the same Stripe82 parent sample would reveal whether the highest-redshift case is a one-off or the bright end of a population.","If the TDE interpretation is right, the same archival method should find more such flares in high-redshift quasars, and the fraction of AGN variability attributable to TDEs could be measured from the parent sample.","The non-virial broad-line interpretation could be tested by high-resolution follow-up spectroscopy: TDE debris should produce velocity-offset or time-variable line profiles tied to the continuum decay, whereas a normal broad-line region would not."],"forward_implications":["SDSS J0001 becomes the highest-redshift optical TDE candidate found in a broad-line quasar at $z=1.404$, extending TDE searches beyond quiescent galaxies.","The $7.5\\times$ gap between the TDE-fitted black-hole mass and the Mg II virial mass implies that broad emission lines in TDE-hosting AGN may contain non-virial debris components, so virial masses in such objects can be biased.","The long observed timescale of the flare follows from the combination of a high impact parameter, a large black-hole mass, and time dilation, so similar long flares at high redshift need not require exotic physics.","The low mock-light-curve probability supports the interpretation that some apparent quasar flares are central TDEs rather than stochastic accretion variability."],"supporting_citations":[{"why":"Supplies the conventional TDE model, including the viscous-delayed fallback accretion scaling relations and blackbody photosphere treatment used for all light-curve fits.","marker":"Mockler et al. (2019)"},{"why":"Provides the fallback-rate templates and the impact-parameter dependence that set the flare's rise and decay shape.","marker":"Guillochon & Ramirez-Ruiz (2013)"},{"why":"Provides hydrodynamical simulations of debris stream evolution that underpin the viscous-delayed accretion template.","marker":"Guillochon et al. (2014)"},{"why":"Defines the CAR/DRW stochastic process used to model quasar variability and to generate mock light curves.","marker":"Kelly et al. (2009)"},{"why":"Supplies the parent sample of 7253 high-redshift Stripe82 quasars and the likelihood criterion for testing whether a light curve is long enough to measure the DRW timescale.","marker":"MacLeod et al. (2010)"},{"why":"Provides the Gaussian-rise/exponential-decay phenomenological fit and the reference TDE sample for comparing timescales, luminosity, temperature, and radius.","marker":"van Velzen et al. (2021)"},{"why":"Supplies the virial black-hole mass from the broad Mg II line that is compared with the TDE-fitted black-hole mass.","marker":"Shen et al. (2011)"},{"why":"Supplies the main-sequence mass-radius relation used to scale the disrupted star's radius in the TDE fallback model.","marker":"Tout et al. (1996)"}],"fun_headline_variants":["Star-shredding event found in quasar at redshift 1.4","Quasar's nine-year flare is a TDE, not AGN noise","Record TDE in broad-line quasar: odds of flicker 0.009%","Quasar flare best explained by star being shredded at z=1.4"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the 0.009 per cent probability computed from mock light curves for this single object, without accounting for the 7253 visually inspected light curves or for fitting the CAR parameters to the same data, is a valid measure of the chance that the flare is intrinsic quasar variability rather than a TDE.","fun_headline_variants_meta":{"raw":{"variants":["Star-shredding event found in quasar at redshift 1.4","Quasar's nine-year flare is a TDE, not AGN noise","Record TDE in broad-line quasar: odds of flicker 0.009%","Quasar flare best explained by star being shredded at z=1.4"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00029,"raw_usage":{"total_tokens":1790,"prompt_tokens":1131,"completion_tokens":659,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":747,"completion_tokens_details":{"reasoning_tokens":573}},"tokens_in":747,"tokens_out":659,"duration_ms":6835,"temperature":1.0,"reasoning_tokens":573,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T05:51:12.263206+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive observation would be continued photometric and spectroscopic monitoring of SDSS J0001 over the next several years. A tidal disruption event should fade smoothly toward the pre-flare level and not re-flare on timescales of years, while the damped random-walk model predicts continuing stochastic variability; the first clear re-brightening or erratic upturn would contradict the TDE interpretation. In addition, a spectrum taken after the flare has faded could check whether the Mg II line width and shift follow the ordinary quasar virial relation, which the debris scenario predicts they should not.","supporting_citations":[],"review_version":1}