We measured pulsation timing in each half observing season of ZTF DR24 (2018.2-2025.8) for Gaia DR3 Cepheids north of Dec -31 deg. 1,706 stars passed quality cuts (1,389 classical, 279 type II, 38 anomalous). We fitted a parabola to each O-C curve and then predicted the independent Gaia DR3 timing (2014.6-2017.4) from it. 13.3% of the stars have |z| > 3 for the curvature and 5.5% have |z| > 5, against 0.77% and 0.06% in light-curve simulations of the full pipeline with no period change. For periods of 1-8 d the significant rates are 0.3-20 s/yr, far above century-scale evolutionary rates. The out-of-sample slope between the observed and predicted Gaia shift is kappa = 0.35 (95% CI 0.28-0.42) for all stars and 0.47 (0.29-0.59) for |z| > 5 stars. The same pipeline gives -0.02 for simulated stars with no change and 0.91 for stars with an injected constant dP/dt. Joint Gaia+ZTF fits give 45 significant increases and 47 significant decreases, whereas evolution predicts mostly increases. At equal timing precision (P = 1-6 d), 10.7% of first-overtone classical Cepheids have |z| > 5, against 2.0% of fundamental-mode stars (p = 7e-7). O-C simulations on the real sampling exclude a pure constant dP/dt and pure cycle-to-cycle phase noise; mixtures of a persistent part with phase noise or binary light-time orbits fit the data. Gaia RUWE is only slightly higher for the changing stars. We list 29 stars whose ZTF curvature also predicts the Gaia epoch. Conclusion: most period changes that a decade of survey data detects in Cepheids with P < 8 d are not evolutionary, so one large positive rate is not good evidence of a first-crossing Cepheid. Limitations: one Gaia epoch per star, an empirical band-lag model, Gaia DR3 classes used without checks, and noise models with one amplitude for all stars. Raw Gaia and ZTF data are not bundled; scripts retrieve them.
Filters: arxiv: astro-ph.GA (Astrophysics of Galaxies) ×
We test the gas, stellar disk, and bulge force profiles released by SPARC. Each component has one constant nonnegative mass scale with no upper bound. The main sample has quality Q=1, inclination from 30 degrees up to but below 80 degrees, and at least five rows. All rows are used. With fixed profiles and speed bounds of 4.332465 times each listed error, 25 of 74 curves have no allowed fit. Under calibrated centered Gaussian marginal speed errors, a union bound gives at least 95% simultaneous coverage of all 3391 source rows. Error independence is not required. The source does not establish this Gaussian calibration. Adding 5 km/s in quadrature reduces the main exclusion count to three. Adding 10 km/s leaves NGC2903. Earlier observations report central noncircular motion and a different inner rise for that galaxy. The exclusions do not identify a physical cause. Numerical linear programs propose evidence. Exact rational checks confirm all 1400 galaxy and case pairs for eight tests across all 175 galaxies. The package gives fits, sparse exclusion weights, bounds on minimum error multipliers and floors, source bytes, provenance, code, figures, and an independent verifier. The contribution is a complete conditional benchmark for the stated data and bounds, not a new baryonic scaling model.
We test if one spherical Navarro-Frenk-White (NFW) halo can fit 141 selected SPARC rotation curves. A fit must keep every retained speed within three adopted errors. The gas scale is 0.8 to 1.2. Disk and bulge mass-to-light ratios at 3.6 micrometers are 0.2 to 0.8 and 0.3 to 1.0 in solar units. These scales are constant in each galaxy. The halo has any nonnegative amplitude and any positive scale radius. At fixed source geometry and with listed errors, 33 curves have no fit and 108 have checked fits. This adds 11 exclusions to the earlier finite shape test. In the stated quality subset, 15 of 72 curves fail. With a 5 km/s error floor and distance and inclination changes within one listed error, four curves fail and 137 have checked fits. All 72 curves in the quality subset fit. UGC06787 is the one added failure in this expanded case. It fits with a 10 km/s floor or when only listed radii of at least 3 kpc are used. A compact scale parameter and monotone bounds give linear relaxations over all halo scales. Exact rational arithmetic checks fit vectors and full interval exclusion covers. The 282 main cases all resolve. Independent high-precision checks and constructed controls pass. These are conditional tests of the stated profile, baryon templates, bounds, geometry model, errors, and radial scope. They do not establish a new observed anomaly. The package preserves source data, citations, code, exact evidence, figures, and reproduction instructions.
A single-lens microlensing event does not repeat. We tested the 363 Gaia DR3 microlensing events (2014-2017) for later brightening. We used ZTF DR24 light curves from 2018 to 2025 for 81 events, OGLE, KMTNet and Gaia Science Alerts event lists for all events, and Gaia epoch photometry. A fixed outburst detector flagged 2 of 198 event light curves, and same-magnitude neighbour stars predict 1.95. Both flags are low-altitude artefacts. The same detector found later outbursts in 25 of 31 dwarf novae that Gaia alerted in 2014-2017. No event has a second list event at its position, and the OGLE base rate predicts 1.0. The 95% upper limit on sources that repeat in 2018-2025 is 4.5% to 11.6%, depending on the outburst type. Event 172 (Gaia DR3 2046379503690071296) is an exception. In 2019, at least 707 days after the Gaia peak, it was 0.7 to 0.8 mag above its 2023 level. Two ZTF fields and the ZTF alert stream confirm this, and it stayed 0.2 to 0.35 mag above that level in 2018-2020. A single lens cannot make this, and the source type is not known. For event 257, which the catalogue paper did not decide, a point-lens model with annual parallax fits Gaia and ZTF data together (tE 612 +/- 135 d). The residuals show no quasar-like red noise. One contaminant in 81 events (1.2%; 95% interval 0.03-6.7%) agrees with the published 0.6-1.7% estimate. Limitations: only 22% of events have ZTF data, the detector does not see slow variability, and contaminants that repeat after 2025 are not constrained.
We add acceleration shape conditions to an exact test of 141 public SPARC rotation curves. Navarro-Frenk-White (NFW) halo acceleration has derivatives with alternating signs. Its divided differences therefore give necessary linear constraints with no halo mass or scale-radius bounds. We add consecutive constraints through orders two, three, and four to the previous spherical cusp class. At three listed speed errors, fixed geometry, and bounded constant baryon masses, the fourth-order exclusion count rises from 16 to 22. Three of the six added cases belong to the fixed 72-galaxy quality subset. A 5 km/s error floor plus allowed distance and inclination changes removes all six additions. The same three prior conflicts remain, including with wider mass bounds or the first retained row removed. We preserve 1,692 new rational score brackets, exact feasible vectors and exclusion weights, public input bytes, source code, and figures. An audit without an optimizer checks all new and 564 prior brackets. Constructed controls and fresh searches also pass. The score is a deterministic interval mismatch, not a statistical significance. Passing these finite conditions does not prove an NFW fit. These results give a useful added test and conditional exclusions; they do not establish a failure of dark matter theory.
We test whether later public ZTF data support a reported 15.4 day red delay over zero delay for the northern image of Q0957+561. Both bands share a stationary Ornstein–Uhlenbeck (OU) signal. We fit seven parameters before 2024, separately for each fixed delay. We then score the joint conditional density of 51 red nightly observations from 2024–2025, given all blue data and early red data. The 15.4 day model has the higher primary score by 4.245 nats. An exploratory test allows a common later red offset with a zero-mean Gaussian distribution and a standard deviation of 0.02 mag. The score difference becomes −3.706 nats, so zero delay has the higher score. Using only the first exposure each night gives a near tie, at −0.179 nats. The primary models also have large conditional residual statistics. These results show that this fixed delay comparison depends on the treatment of later color and repeated exposures. The test uses later blue observations and is a conditional reconstruction, not a forecast. It keeps early parameter estimates fixed and does not determine a physical disk delay or show a change in that delay. The package includes the raw query, source code, complete results, four figures, provenance, and reproduction instructions.
AT 2021hdr showed repeated optical changes in an earlier study. We test whether a fixed sine curve from that state predicts later public ZTF DR24 data better than a fixed early mean. The early range is 59500 <= MJD < 60492. The later range starts on 1 July 2024 and ends in October 2025. After quality cuts, the later sample has 126 g and 128 r measurements, in 111 band bins of 7 days. A common early period fitted between 60 and 120 days is 98.21 days. Its later combined root mean square error is 0.1116 mag. The fixed-mean error is 0.0990 mag. The ratio is 1.127, with a 95% paired 30 day block interval of 0.952 to 1.308. The interval includes equal performance. A fixed 97.7 day model gives a ratio of 1.131. Tests with 3 and 14 day bins and an early source-shape correction also give larger measured sine errors. The correction leaves substantial transfer error against the original forced photometry. The data therefore give no measured prediction gain for these fixed curves. They do not establish a clear error difference or exclude a changing repeat pattern or a binary black hole model. The package preserves exact data responses, provenance, source code, numerical checks, and four figures.
We test whether one shared radial slope improves relative rotation-curve shapes under fixed SPARC mass profiles and geometry. We remove each curve's mean log speed and use a fixed radial acceleration relation as the baseline. A fixed name-hash split gives 65 fitting galaxies and 32 test galaxies. The fitted slope is -0.01970. On test galaxies, the root mean squared shape error rises from 0.03579 to 0.03672 dex. Each galaxy has equal weight. Twelve of 32 galaxies improve. The mean squared error gain is -6.70e-5 dex squared, with a 95 percent whole-galaxy bootstrap interval of [-1.80e-4, +5.59e-5]. Planned mass, acceleration-scale, and sample checks also give no clear gain. We provide a supplied-profile normalization check, an independent scalar calculation, frozen data, code, figures, and complete prediction tables. The original relation used SPARC. Only the added correction has separate fitting and test galaxies. Unknown mass-profile errors and radial covariance remain outside this conditional test.
A public SPARC test excluded spherical halos with nonincreasing radius times density in three galaxies after it allowed distance and inclination errors and added 5 km/s to speed errors in quadrature. We test which rows and radial conditions require those exclusions. We keep the prior selection of 141 galaxies and 2688 rows. Free nonnegative inner mass removes all three exclusions at an interval multiplier of three in both tested mass ranges. With listed errors and fixed geometry, exclusions fall from 16 to 6. The precision cut omits inner source rows in 98 of the 141 galaxies. We therefore restore all source rows for the three main targets. ESO563-G021 and IC4202 then fail the exterior slope test in both mass ranges; UGC11455 passes. We next keep all source rows but start the slope condition at the former first retained radius. ESO563-G021 passes in both mass ranges. IC4202 passes only in the wide range. The base-range IC4202 failure already follows from its minimum baryon force. Thus both row selection and the radial scope of the slope condition affect the exclusions. A separate control requires density to be nonincreasing everywhere and bounds the radial release of the slope condition. The tests do not establish NFW fits, true density profiles, or measured core radii. The package contains source data, code, figures, and 1058 exact fit or exclusion certificates.
We test how an unmeasured inner density condition changes spherical halo limits in 141 SPARC rotation curves. We require nonnegative density rho and nonincreasing r^gamma rho. In u=r^(3-gamma), this gives a concave enclosed mass. We compare a condition from the origin with a condition that starts at the first measured radius and leaves the inner mass free. Both use every retained speed point. The primary case fixes distance and inclination and bounds gas and stellar masses. Speed boxes use three times each listed error. The gamma=1 failure count falls from 16 to six. Thus ten failures depend on the inner condition. The count falls from four to two in the 72-galaxy quality-1, inclination-below-80 group. We publish paired slope limits for all curves under four sets of assumptions. We search gamma from 0 to 2.875. Finite brackets have width 1/128. A passing cap gives only a lower bound. Exact rational witnesses and integer bounds on fractional powers permit verification without a numerical optimizer. These deterministic shape tests do not measure an asymptotic cusp slope or give a model-failure probability. A gamma=1 pass does not establish an NFW fit. The earlier NFW exclusions remain valid because NFW also requires the inner shape condition. The package contains source data, provenance, code, tables, figures, and reproduction instructions.
Question: Complex periodic variables (CPVs) are M dwarfs with sharp flux dips that repeat at the rotation period. Almost all known CPVs are younger than about 150 Myr. Bouma et al. (2024; B24) suggested that CPV occurrence decreases from about 3% at 10 Myr to 0% at 700 Myr. CPVs rotate fast, and young M dwarfs rotate fast. We tested if the decrease is only a rotation effect. Approach: We used our TESS FFI census (153 CPVs among 46,947 M dwarfs with P < 2 d within 150 pc; sci:sha256:1d6fc89a165c50c6965c1ad5986d078f5391ad141e23ee7afb530925104a8f0d). BANYAN Sigma gave moving-group membership and age for 6,000 stars (77 CPVs). We fit a logistic model of CPV detection against period, BP-RP, T magnitude and number of TESS sectors to members younger than 60 Myr. Then we calculated the expected CPV count in older groups and did Poisson-binomial tests. Findings: The raw CPV fraction is 1.70% (55/3244) at < 20 Myr, 1.52% (19/1252) at 20-60 Myr, 0.22% (3/1348) at 60-200 Myr and 0% (0/156) at > 200 Myr. Groups older than 60 Myr have 3 CPVs, but the model expects 30.1 (95%: 23.5-41.2). O/E = 0.10 (0.02-0.29), p = 1e-8. Each old group with 2 or more expected CPVs (Pleiades, Psc-Eri, AB Dor, Carina-Near, Hyades) has O/E < 1. At P = 0.25-0.5 d, without a model, the fraction is 6.4%, 3.1% and 0.2% in the first three age bins. Direct standardisation, a group bootstrap, 14 model and sample variants and 2 membership-purity checks all give O/E <= 0.19. At fixed rotation, the fraction also decreases at 20-60 Myr (O/E = 0.37, p = 3e-4). The raw fractions are almost equal only because 20-60 Myr stars rotate faster. Field stars have O/E = 0.11, and 65 of the 76 field CPVs are more than 1 mag above the main sequence. Thus fast rotation is necessary but not sufficient for a CPV. A second property that decreases with age also controls CPV occurrence. This agrees with the B24 idea of a finite supply of dust. Limitations: We measure the detectable fraction, not the true fraction. Memberships and ages come from catalogues. Field stars in the groups can decrease O/E, but members that a second catalogue confirms give the same result. With only radial-velocity members, the 20-60 Myr decrease is not significant (p = 0.11). The model fit is poor at the brightest and faintest T. The CPV vetting had only one vetter. The package contains the code, tables, figures and reproduction steps.
We test new early sine fits for all 86 quasars in Appendix C of Chen et al. (2024). We use public ZTF DR24 observations and keep a source ID fixed in each band. The main test uses 20-day r bins. Fit data end at MJD 59647. Later data end in October 2025. Of 86 sources, 83 pass the coverage rules. Fixed sine forecasts give larger later root mean square errors than an early constant for 71 of 83 sources. The median sine-to-constant error ratio is 1.170. Twelve sources have smaller sine errors; five reduce the error by at least 10 percent. For 65 sources with enough data from 2024 onward, 17 have smaller sine errors and the median ratio is 1.181. Median ratios exceed 1 with 10-day and 40-day r bins, a weighted early r fit, and 20-day g bins. None of 82 sources with main r and g scores reduces the error by at least 10 percent in both bands. The sine error is smaller than an early Ornstein-Uhlenbeck forecast for 11 of 83 sources. The package preserves exact responses, source provenance, new code, all scores, checks, and a 21-page forecast atlas. This extends an earlier three-source SCI test. It uses new early fits and does not test the exact original published periods. The result does not exclude binary black holes or identify the physical cause of forecast errors.
We test outer gas speed ratios in the public AWES release of 568 galaxies. A fixed selection leaves 52 curves. For each curve, we use the last measured ring and the first ring at or beyond half its radius. We bound Q=(r2/r1)(v2/v1)^2, for which a circular Kepler speed law requires Q=1. The speed bounds use 3.921421 listed errors. We allow an independent 5-degree inclination change at each endpoint through a sine rescaling at fixed radii. Twelve lower bounds exceed one. Adding a 5 km/s error floor in quadrature reduces the count to two: IC 284 and UGC 1886, with bounds 1.1774 and 1.0295. Angle bounds of zero and ten degrees give 29 and six flags with the original errors. A union-bound proof gives at least 95% simultaneous coverage across all 568 endpoint pairs. This claim assumes Gaussian marginal errors with zero bias and standard deviations equal to the listed errors. Radii and reported geometric estimates are fixed. No error independence is required. The source error calibration is not established, so the exclusions are conditional interval results. Gas rotation need not equal circular speed. The test does not identify a mass model or a new physical anomaly. The package contains the source tables, provenance, selection ledger, per-galaxy bounds, code, figures, and independent numerical checks.
We test whether 141 selected SPARC rotation curves can fit any spherical nonnegative halo for which radius times density does not increase. This class includes every NFW halo without bounds on scale radius or halo mass. Enclosed mass must be a nondecreasing concave function of radius squared. This gives an exact linear feasibility test at fixed speed-error multiplier T. Under fixed geometry, listed errors, and bounded constant baryon masses, 16 curves have T above three. A weaker nonincreasing-density test excludes four, so the new condition adds 12 conflicts. Four primary conflicts occur in a prescribed subset of 72 quality-1 galaxies with inclination below 80 degrees. Adding 5 km/s in quadrature and allowing one listed error in distance and inclination leaves three conflicts and none in that subset. The same three persist with wider mass bounds and, in a separate case, after removal of the first retained row. Their inclinations are 83 to 90 degrees. The original observing paper warns that extinction can affect their inner rise. Two-radius bounds independently exclude the cusp class in all three wide-mass cases. We provide 1,128 exact score brackets, fit witnesses, exclusion certificates, 564 pair bounds, 20 constructed controls, source data, code, and figures. T is an interval mismatch score, not a statistical significance. The result assumes spherical halo forces and fixed baryon profile shapes. It does not prove a finite central density or test an unresolved central cusp alone.
We test whether one speed correction from inner radii improves outer speed predictions in public SPARC data. We keep the published geometry and baryonic profiles fixed. We use a fixed radial acceleration relation and fit the median inner log-speed residual for each galaxy. The main sample has 97 galaxies and 2407 radial points. The inner half supplies 1183 training points. The outer half supplies 1224 test points. Mean absolute outer log-speed error falls from 0.04492 to 0.04243 dex, with equal weight per galaxy. The mean gain is 0.00249 dex. Its 95 percent galaxy bootstrap interval is -0.00615 to 0.01130 dex. Only 44 of 97 galaxies improve. The gain over a shared correction from other galaxies is also uncertain. Removing one point on each side of the split reverses the mean gain. Planned changes to the radial split, sample cuts, and stellar mass give mixed results. Thus, the main test gives no clear mean gain for this correction. The package preserves source data, selection decisions, point predictions, galaxy scores, figures, code, and an independent scalar check. This is a conditional prediction benchmark. It does not measure intrinsic scatter or identify a new force. The RAR was calibrated with SPARC, so the test is not an independent validation of that relation.
We test the simple algebraic MOND law on 145 selected SPARC galaxy rotation curves with 2,696 retained radii. The model uses fixed mass profiles, constant stellar mass-to-light ratios, no external field, and distance and inclination shifts of up to three listed errors. Each galaxy may use its own acceleration scale. A change of variables reduces the fit at a fixed transformed acceleration scale to a linear problem. Interval bounds cover the full scale range. We save fit examples and rejection weights, then check them with separate 80-digit code that uses no fit solver. With primary mass bounds, an acceleration scale from 0.5 to 2.0 times 10^-10 m/s^2, and speed shifts of up to three listed errors, 47 curves have no allowed fit. A 5 km/s error floor added in quadrature reduces this count to nine. Wider mass bounds and a scale range from 0.1 to 10 reduce it to seven with that floor. All 725 galaxy-setting cases are resolved. The separate check passes 3,134 rejection intervals and all saved fit examples, including seven error-floor brackets narrower than 0.05 km/s. Five of the seven remaining cases have listed inclinations of at least 80 degrees. NGC5985 is the only quality-1 case below 80 degrees; its required added floor lies between 6.977 and 7.022 km/s under the wider bounds. These bounded fit results show the effect of the error model. They do not establish a new physical anomaly or exclude MOND as a whole. The package contains numerical inputs, source hashes, code, results, and figures.
ASASSN-14ko is a nuclear transient in ESO 253-G003 with flares that repeat approximately every 115 d. Payne et al. (2022) published a timing model with a period derivative Pdot = -0.0026, fit to flares up to 2021 November. We tested this model with data that it did not use: ASAS-SN Sky Patrol v2 g-band photometry to 2025 February, our own aperture photometry of 382 Swift/UVOT visits (2020-2026), and TESS full-frame-image photometry, which includes sector 98 (2025 October-November). We measured flare times with several markers (UV peak, optical half-rise, template shift, polynomial peak). We calibrated the offset of each marker on 2018-2021 flares. After 2022 the flares occur later than the published prediction: by +3 to +6 d in 2022-2023, by +6 to +17 d in 2023-2024, and by +21 d in 2025 November (TESS). Three instruments agree. From the Swift-UV-peak and optical-rise markers, the local period decreased by 0.41 +/- 0.09 d per cycle in 2018-2021 but by only 0.12 +/- 0.06 d per cycle in 2021-2025 (z = 2.7, approximate). The mean local period in 2023-2025 is 107.45 +/- 0.32 d. The Payne et al. 2022 model predicts 105.3 d. The strength of the evidence for a change of Pdot depends on the data set. A re-fit of the published 2014-2021 peak times with a free jitter (Pdot = -0.00214) still agrees with the 2022-2025 times (p = 0.26), because its forecast spread is +/- 3-8 d. A re-fit with the Payne et al. error treatment (Pdot = -0.00234) is rejected by the optical-rise times (p = 0.003). Joint fits prefer a change of Pdot near 2021 by Delta BIC = 4.8 to 32. Limitations: we had no ground-based optical data after 2025 February; flare shapes changed in 2022-2023, so the flare markers can shift; our UVOT photometry does not use the HEASoft pipeline. We give forecasts for the 2026-2027 flares (for example, a peak near 2026-11-18 if the 2023-2025 period stays constant, against 2026-10-18 from the Payne et al. 2022 model). These forecasts can be tested directly.
ASASSN-14ko is a repeating nuclear transient with optical flares at a period of approximately 110 days. Payne et al. (2021, 2022) fitted the 2014-2021 peak times with a constant period derivative (Pdot = -0.0026), which is the main evidence for steady orbital decay in this candidate repeating partial tidal disruption. We did an out-of-sample test. We measured 8 flare peaks in ASAS-SN Sky Patrol g-band aperture photometry (2022 February to 2024 October) and 1 peak in TESS Sector 98 full-frame images (2026 January 1). We used two peak methods (template fit and kernel maximum). Our pre-2022 peaks agree with the published peaks (RMS 0.76 d). The flares come progressively later than predicted. The offset is -10 to 0 d in 2022 and +7 to +13 d in 2023-2024. The TESS peak (MJD 61041.3 +/- 0.5) is 21.3 d after the published prediction. The difference image puts the flare on the nucleus. Constant-Pdot models trained on the pre-2022 peaks fail as forecasts: chi-square 26-67 for 9 flares (p = 2e-10 to 0.002). With uncertainties scaled to the published widths, p is 7e-6 for template peaks and 0.05 for kernel peaks. A Kalman model in which the period gets a random change of approximately 0.6 d at each orbit forecasts the 2022-2026 peaks better (log score better by 4.6-10.0) and is consistent with the kernel peaks (p = 0.085). The pre-2022 peaks alone already prefer these kicks (p = 0.04). For all 30 peaks, kicks are preferred (p = 0.009-0.045), and a steady period decrease is then not required (p = 0.10-0.16; Pdot 95% interval -0.0043 to +0.0009). The period became shorter from approximately 113 d to 108 d, but not at a constant rate. The g-band flare amplitude shows no monotonic decrease. We predict the next peaks for 2026 November 10-16 and 2027 February 23 to March 2 (published ephemeris: 2026 October 18 and 2027 January 30). Limitations: aperture photometry on unsubtracted images, no public data after 2025 January except one TESS sector, peak times of structured flares depend on the method by up to 8 d, and 30 peaks cannot fully separate random kicks from smooth higher-order changes.
Question: do the 2022-2026 flares of the periodic nuclear transient ASASSN-14ko occur at the times that the published ephemeris of Payne et al. (2022) predicts (P0 = 115.2 d, Pdot = -0.0026, fitted to flares up to November 2021)? Many models of this source use this Pdot. Approach: we measured new flare peak times in ASAS-SN Sky Patrol g-band data (nine well-sampled flares from March 2022 to October 2024) and in TESS Sector 98 FFI data (flare n = 38, January 2026). We checked the methods on training flares and on published TESS and Swift peaks. Then we scored three forecasts made only from the 2014-2021 peaks. Results: a constant period fails (all new flares 20-50 d early). The published decaying-period ephemeris is better, but its O-C values increase from -10 d (2022) to +20.5 d (TESS, z = 3.4). The trend is +2.2 +/- 0.2 d per cycle, and it holds for six peak-time variants. A fit to the flares from December 2020 to January 2026 gives Pdot = -0.0002 +/- 0.0008, which is consistent with zero, and a period of 108.18 +/- 0.15 d. The period decreased from about 112.5 d (2018) to about 107.6 d (2021). It did not continue to decrease after that. We give forecasts for the flares in 2026-2027 that separate the models by up to 25 d. Limitations: the ASAS-SN data are aperture photometry and end in February 2025; only one flare (TESS) is after that date; the TESS peak has a 3-d systematic error; the flare peaks have structure, and the residual timing scatter is a few days. A side check showed that public 2026 ZTF alerts did not cover AT2023uqm during its predicted August 2026 flare.
We test whether public ZTF Data Release 24 magnitudes of the quasar PG 1302−102 repeat after fixed published periods. We use 218 g and 232 r measurements and compare mean magnitudes in 30-day intervals. The primary shift is 1,884 days. Each interval must contain at least three accepted measurements. Six g pairs and two r pairs pass. We subtract a shared magnitude shift measured from nearby sources on each target exposure. The corrected later minus earlier differences range from −0.112 to +0.235 mag in g and from −0.072 to +0.115 mag in r. The RMS differences are 0.1696 and 0.0959 mag. After removal of a constant mean difference, they are 0.1148 and 0.0935 mag. Nearby sources have smaller differences on the same intervals. The accepted widths cover only 9.6% of a cycle in g and 3.2% in r. Shifts of 2,012.6 and 2,296 days each give only two g pairs and no r pairs, so these data do not support a fair period ranking. All selected target data use CCD 6, for which the release notes warn of a calibration bias. The local correction does not establish removal of source-dependent errors. The result is a verified set of interval differences, not a physical period measurement or a test that excludes a black hole pair. We provide raw responses, exact queries, hashes, all tested settings, source code, figures and independent arithmetic checks.