{"id":"0d4aa084-7c1d-4495-bf2b-b3e4a7bacb08","arxiv_id":"2501.04866","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Flare occultations observed with optical photometry can constrain the latitude of M dwarf flares, with an estimated 3 to 22 detectable events expected in existing TESS data.","lead":"This paper develops a method to locate stellar flares using the brief dimming caused when a transiting planet or eclipsing companion passes in front of an ongoing flare, and applies the method to estimate how many such events should already exist in TESS data. It also presents one unconfirmed candidate event in the eclipsing binary CM Draconis that could hint at an equatorial flare location.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 3–22 detection estimate and the method's applicability rest on an uncalibrated 'detectability' criterion (no noise model) and a single statistically unconfirmed candidate, given the known intrinsic flare morphologies that mimic occultations.","rationale":"The reader identified the uniform flare distribution as the weakest assumption. That assumption is acknowledged in the paper and mainly affects the predicted probabilities, not the geometric inversion itself. A more load-bearing concern is that the entire 'detectability' claim—both the method's applicability and the 3–22 event estimate—has not been validated end-to-end. The simulation's detectability criterion ignores photometric noise, so 'detectable occultations' are really geometric occultations that may be invisible in actual TESS data. Even setting noise aside, the paper's own Sect. 7.1 lists numerous intrinsic flare behaviors that mimic the predicted occultation morphologies, and the only empirical candidate is not statistically significant (4.8% false-positive probability, and not an outlier among 125 flares in Fig. 11). Without an injection-recovery study or a confirmed detection, the central claim that optical occultations can be reliably identified and used for latitude constraints is not yet established. This does not change the reader's CONDITIONAL verdict—the method is plausible and worth developing—but it identifies a different and arguably more fundamental gap than the uniform distribution assumption. The concrete test proposed would directly quantify this gap.","tokens_in":22694,"tokens_out":12201,"duration_ms":117437,"concrete_test":"Perform an injection-recovery test on the TESS 2-minute and 20-second light curves of CM Dra and GJ 3236. Inject synthetic occultation signals generated by the paper's own simulation at random times, with amplitudes drawn from the same 3–20% distribution used in the paper, add the real photometric noise, then run the same detrending and visual/flare-fit detection procedure used in Sect. 5. Measure both the recovery fraction of injected events and the false-positive rate from intrinsic flares. If the recovery fraction is low (e.g., below 50%) or if intrinsic flares produce false positives at a rate comparable to the expected 3–22 event rate, then the estimate and the method's empirical demonstration would require substantial revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's 'detectable occultation' criterion in Sect. 2.3 is purely geometric: a >10% occultation of the flare area, and per Sect. 4.5, at least two cadence points during the occultation. No photometric noise is included in the Monte Carlo simulations. For the fiducial low-amplitude flares (3% of stellar flux), a 10% area occultation produces a dip of only 0.3% of the stellar flux, comparable to or below TESS 2-minute noise for typical M-dwarf targets. The 3–22 event prediction in Sect. 6 therefore counts geometric events, not events that would actually be detected in the data. Furthermore, the only empirical demonstration, the CM Dra candidate in Sect. 5.5, has a false-positive probability of 4.8% when compared to the Howard & MacGregor (2022) flare sample, and the paper itself catalogs multiple intrinsic flare substructures (sympathetic flares, quasi-periodic pulsations, peak-bump flares, rapid-decay flares) in Sect. 7.1 that mimic the predicted occultation morphologies. The candidate is not a statistically significant outlier among the 125 CM Dra flares (Fig. 11). Thus the central claim that the technique 'can be used' and that 3–22 events are 'detectable' currently rests on an uncalibrated detection criterion and a single unconfirmed candidate.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes using occultations of flares by transiting planets or eclipsing binaries in optical photometry to constrain flare latitudes. It models transit geometry and flare areas, uses Davenport/Tovar Mendoza flare templates to generate light-curve morphologies, and runs Monte Carlo simulations to compute occultation probabilities as functions of radius ratio, impact parameter, flare rate, orbital period, and cadence. It applies the framework to TRAPPIST-1, AU Mic, CM Draconis, GJ 3236, and selected TOIs, searches TESS short-cadence data of GJ 3236 and CM Draconis for in-eclipse flare occultations, presents one candidate in CM Draconis, and estimates that 3–22 detectable occultations should exist in TESS primary mission data, mostly in eclipsing binaries.","tokens_in":22885,"tokens_out":7659,"duration_ms":73300,"significance":"The geometric idea is timely and potentially valuable: if validated, it would convert existing TESS photometry into direct constraints on M dwarf flare latitudes, complementing rotation-modulation methods. The simulation is clearly described, the use of established flare templates is appropriate, and the paper is candid about limitations, repeatedly labeling the CM Draconis event a candidate and cataloging false-positive morphologies and spin-orbit caveats in Sections 7.1 and 7.5. The probabilistic framework is internally consistent. However, the headline estimate and the demonstration rest on an uncalibrated geometric detectability criterion, ad hoc flare parameter distributions, and a single unconfirmed candidate; the numerical predictions are not yet supported at the level claimed.","major_comments":[{"comment":"The “detectability” criterion used to produce the headline 3–22 estimate is purely geometric. In §2.3 an occultation is counted when at least 10% of the flare area is covered, and §4.5 adds the requirement of at least two in-occultation cadence points; no photometric noise, systematics, or detection-significance threshold enters the Monte Carlo. For the low-amplitude flares in the simulation (3% of stellar flux, §2.2), a 10% area occultation corresponds to a 0.3% flux dip, comparable to or below the TESS 2-minute photometric precision for typical faint M-dwarf targets. Consequently, the events counted in Table 2 and the “3–22 detectable occultations” statement in §6.2 are geometric occultations, not detectable events. The central quantitative claim needs either an explicit noise model with an SNR threshold or a clear reframing as geometric occultations whose actual detectability remains to be calibrated.","section":"§2.3, §4.5, §6.2"},{"comment":"The flare ensembles are generated from a vaguely described “seeded value from a normal distribution” without quantitative calibration. The amplitude range (3–20%), FWHM thresholds (10–45 min for low-amplitude flares, 2–4 h for high-amplitude flares), and Earth-to-Jupiter flare radii are stated without justification or reference. Since P_occultation and the expected counts in Table 2 depend on the adopted flare area, duration, and amplitude distributions, the 3–22 estimate is sensitive to these unvalidated inputs. The authors should calibrate the distributions to observed flare samples (e.g., Günther et al. 2020b; Howard & MacGregor 2022) or at least provide a sensitivity analysis over plausible parameter ranges.","section":"§2.2"},{"comment":"§3 assumes a spherically uniform flare distribution and §6.1 carries this assumption into the sample predictions. The assumption is stated, but the probabilities for individual systems in Table 1 and the expected-count range in Table 2 are all computed under it. Given that Ilin et al. (2021) infer polar flare latitudes for rapidly rotating M dwarfs, and §7.5 itself questions whether tight binaries have active longitudes, the uniform assumption is not benign. The authors should add sensitivity tests with polar-concentrated, equatorial, and banded latitude distributions and show how P_occultation, Table 1, and Table 2 change; without this, the predicted 3–22 range should be presented as conditional on a uniform latitude prior.","section":"§3, §6.1"},{"comment":"The CM Draconis demonstration does not currently validate the method. As the paper states in §5.5 and §7.3, the candidate has at best a 4.8% tail probability relative to the Howard & MacGregor (2022) decay/FWHM distribution and is not a strong outlier among the 125 flares in Fig. 11. Moreover, §7.1 lists sympathetic flares, quasi-periodic pulsations, peak-bump flares, and rapid-decay flares that produce morphologies matching the predicted occultation signals. The text nevertheless says in §5.5 that the candidate “strengthens the case” and later uses it to infer an equatorial flare location. The demonstration would need a full false-alarm analysis calibrated on the non-occulted flare sample, or the equatorial-flare inference should be explicitly labeled as a speculative illustration rather than an empirical result.","section":"§5.5, §7.3"}],"minor_comments":[{"comment":"Equation (4) defines erfc but the text says “the error function for which we use the Scipy package scipy.special.erf”; this is inconsistent, and the correct function used in Eq. (2) should be stated.","section":"Eq. (4)"},{"comment":"The text says the candidate resembles morphology “Fig. 1c”, but §5.5 and the caption of Fig. 1 describe the candidate as matching morphology (b); the reference should be corrected.","section":"§7.3"},{"comment":"In §6.1 the authors state that for eclipsing binaries “the primary and secondary stars will have the same radii”, but the M4–M6 eclipsing-binary row of Table 2 lists R2 = 0.446 R☉ against R1 = 0.196 R☉; if this is a typo it should be corrected, since it affects the expected count for that bin.","section":"Table 2"},{"comment":"The requirement of at least two data points during the occultation is introduced in §4.5 but is not stated in the detection criteria in §2.3; the criteria should be defined consistently in one place.","section":"§4.5"},{"comment":"The caption calls 4.8% a “false positive chance,” but §7.3 describes a percentile within a comparison flare sample; these are not equivalent, and the caption should be reworded.","section":"Fig. 10 caption"},{"comment":"The paper states that simulation and observational data “may be provided upon reasonable request”; releasing the simulation code would improve reproducibility and allow readers to check the Monte Carlo calibration.","section":"Data Availability"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Key thing to know: the geometric machinery for flare occultations is sound and genuinely new in the optical/exoplanet context, but the headline number of 3-22 detectable events is not calibrated against photometric noise. Read it as a search strategy plus an optimistic upper bound, not as a detection prediction.\n\nWhat's new: it extends Sanz-Forcada's X-ray eclipse-occultation idea to optical TESS photometry and transiting exoplanets, and it quantifies which systems, cadences, and impact parameters are favorable. The Monte Carlo simulation with circle-circle occultation geometry is straightforward and internally consistent. The parameter study is useful. The search of CM Dra and GJ 3236 is honest work: they find one candidate in CM Dra, stack all 125 flares, note that the rapid decay is not a big outlier, report a 4.8% false-positive rate, and explicitly do not claim confirmation. The discussion of degeneracies with sympathetic flares, quasi-periodic pulsations, and peak-bump flares is thorough and shows real care.\n\nSoft spots. The main one is the 'detectable' criterion in Sect. 2.3: a 10% area occultation and at least two cadence points, purely geometric, with no noise model. For a fiducial 3% flare, that is a 0.3% flux dip, right at the TESS 2-minute noise floor for typical M dwarfs. So the 3-22 estimate is really a count of geometric occultations that might be recoverable, not events an actual pipeline would find. An injection-recovery test into real TESS light curves would fix this, and the language should be softened until then.\n\nSecond, the uniform flare-latitude prior is load-bearing for the individual probabilities and for the equatorial reading of the CM Dra candidate. Stated clearly, but the numbers would move if M dwarf flares are pole-concentrated (Ilin et al. 2021). Third, the simulation inputs—flare amplitude range, FWHM thresholds, radius ranges, the 10% threshold, the 4x FWHM window—are hand-set without sensitivity testing. Not fatal, but a robustness check would help. Finally, no code or data is shipped, and the data availability statement is 'on reasonable request'—a reproducibility miss.\n\nWho is it for: anyone working on M dwarf flare latitudes and habitability. The method could turn TESS archives into a new constraint on flare locations, so it deserves referee time. My recommendation: peer review, yes, and I'd push for a major revision that adds a noise model or at least rescales the 'detectable' claim, plus sensitivity tests on the priors. The core geometry holds up; the quantitative overreach is fixable.","headline":"The occultation geometry is sound, the 3-22 'detectable' estimate is not noise-calibrated, and the CM Dra candidate is an honest non-confirmation.","tokens_in":23543,"tokens_out":3917,"would_cite":true,"duration_ms":37621,"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":"This paper argues that when a transiting planet or eclipsing companion passes in front of an ongoing flare, the timing and shape of the resulting light-curve dip can reconstruct the flare's latitude and longitude on an M dwarf, and that 3…","keywords":["flare occultations","M dwarfs","stellar activity","exoplanets","eclipsing binaries","TESS","flare latitudes","transit photometry"],"falsifier":"A systematic search of all TESS 20-second and 120-second data for the 55 M-dwarf eclipsing binaries defined in the paper: if the number of confirmed type-B occultations after a complete search is statistically consistent with zero while the model predicts at least three, the uniform-flare assumption or the simulated flare morphology templates would be falsified.","tokens_in":22380,"feed_emoji":"🔭","tokens_out":8609,"duration_ms":75022,"temperature":0.7,"pith_summary":"The paper proposes a new way to locate flares on M dwarf stars: watch for a transiting planet or eclipsing companion to pass in front of an ongoing flare, blocking part of its light. From the timing and geometry of that occultation—the impact parameter, the sizes of the two bodies, and how deep into the eclipse the dip occurs—the flare's latitude and longitude can be reconstructed. The authors simulate the expected light-curve shapes, compute detection probabilities for known systems, and estimate that 3–22 detectable occultations should already be hiding in TESS primary-mission photometry, most of them in eclipsing binaries. They analyze one candidate in the CM Draconis system, showing that its unusually fast decay is consistent with an occultation but is not yet statistically confirmed. If the technique works, it gives a direct, purely photometric way to map where M dwarfs flare, which matters because equatorial flares are much more threatening to orbiting planets' atmospheres than polar ones.","feed_headline":"Transit shadows can map M dwarf flare latitudes","feed_subtitle":"Model predicts 3–22 flare occultations hiding in TESS data, mostly in eclipsing binaries.","key_machinery":"The core is a geometric Monte Carlo simulation of a transit across a star with circular flare patches. The occulted area is computed with the circle-circle intersection formula (Eq. 1), which takes the distance between the flare center and the transiting body's center with their radii and returns the overlap area; an event is recorded as detectable when the overlap covers at least 10% of the flare's area, occurs after the flare's peak and within four FWHM durations, and leaves at least one data point with a partially visible flare. Flare light curves are generated with the Tovar Mendoza et al. (2022) template, producing three recognizable morphologies: a dip in decay (A), an abrupt cut-off of decay (B), and a rising-phase occultation that mimics a low-amplitude flare (C). The probabilities run 50,000 transits per parameter set, with $P_{\\rm occultation} = N_{\\rm success}/N_{\\rm runs}$, then $P_{\\rm obs} = 1-(1-P_{\\rm occultation})^{N_{\\rm transit}}$ for multi-transit windows. Flare positions are drawn from a uniform distribution across the stellar sphere—an explicitly stated modeling assumption that does not enter the geometric method itself.","core_discovery":"The paper's central claim is that an occultation event—a transiting body passing in front of an ongoing flare—can be identified in optical photometry and converted into a latitude and longitude measurement for the flare. The geometry is fixed by the known transit parameters: the impact parameter gives the latitude band swept out, and the time of the occultation within the transit gives the longitude. Using this, the authors build a simulation that generates realistic flare light curves, compute single-transit and per-sector occultation probabilities for known M dwarf systems, and find that eclipsing binaries dominate the expected yield. The quantitative headline is that 3–22 detectable occultations should already be present in TESS primary-mission photometry. The CM Draconis candidate is presented as a demonstration of the whole chain—geometry, light-curve morphology, and model comparison—but explicitly not as a confirmed detection.","pith_inferences":["If the uniform-flare assumption is wrong, the predicted counts shift: polar-concentrated flares would lower the occultation rate for low-impact-parameter transits, so a carefully defined null result from a TESS search would itself constrain the flare latitude distribution.","The same geometric machinery applies to other high-cadence photometric surveys and to stars of other spectral types, since nothing in the geometry depends on spectral type once radii, orbital parameters, and flare rates are known.","By stacking occultation chords from multiple transits of the same star, one could build a coarse 'flare latitude histogram' that distinguishes equatorial, intermediate, and polar flare populations without high-resolution stellar imaging.","The visual search could be automated by fitting every in-transit flare with an injected-occultation template and comparing Bayesian evidence, which would turn the technique into a systematic population survey rather than a by-eye candidate hunt."],"forward_implications":["Eclipsing binaries are the prime targets: CM Draconis reaches a 99.8% chance of having at least one occultation across its 16 TESS sectors, and GJ 3236 an 87.4% chance across its 4 sectors.","For transiting planets, larger radius ratios and shorter periods give the best odds; Neptune-sized planets with 1-day periods could produce up to 3 detectable occultations in current TESS data, while Earth-sized planets are essentially hopeless in a single sector.","Cadence matters: 20–120 second data can catch occultations, but 10–30 minute cadence reduces the probability to effectively zero, so searches should prioritize short-cadence observations.","The three predicted light-curve morphologies give observers a template: type B (a flare that disappears mid-decay) is the most reliable occultation signature, while type A risks confusion with sympathetic flares and quasi-periodic pulsations, and type C is indistinguishable from a low-amplitude flare.","If the CM Draconis candidate is real, it would be a near-equatorial M dwarf flare, a data point set against the polar latitudes inferred from rotation-modulation studies."],"supporting_citations":[{"why":"First demonstrated flare occultations as a latitude probe at X-ray wavelengths; the present optical method is a direct extension of this idea.","marker":"Sanz-Forcada et al. 2006"},{"why":"Detailed the geometric modeling of X-ray flare occultations that this paper's simulation adapts to optical photometry.","marker":"Sanz-Forcada et al. 2007"},{"why":"Reported polar flare latitudes on rapidly rotating M dwarfs, the main comparison point for the CM Draconis candidate and the motivation for constraining flare latitudes.","marker":"Ilin et al. 2021"},{"why":"Provided the classic white-light flare template whose morphology this paper's simulated occultation light curves build on.","marker":"Davenport et al. 2014"},{"why":"Supplied the updated flare morphology model used to generate simulated occultation shapes and to fit the CM Draconis candidate.","marker":"Tovar Mendoza et al. 2022"},{"why":"Gives the M dwarf flare occurrence rates used to bin the TESS sample and compute expected occultation counts.","marker":"Günther et al. 2020b"},{"why":"Defines the TESS eclipsing binary catalog from which the binary sample for the expected-count estimate is drawn.","marker":"Prsa et al. 2022"},{"why":"Provides the TESS Objects of Interest flare catalog used to select the transiting planet sample.","marker":"Howard 2022"},{"why":"Identified the 125 CM Draconis flares, including the candidate occultation event that this paper analyzes in detail.","marker":"Martin et al. 2023"}],"fun_headline_variants":["Transit shadows map M dwarf flare latitudes","Occultations pinpoint flare latitude on M dwarfs","TESS data likely hides 3–22 flare occultations","Flare occultations: a map for M dwarf activity","Eclipsing binaries reveal more flare occultations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The numerical predictions assume flares are spread uniformly over the star's surface; if M dwarf flares actually prefer high latitudes, the expected 3–22 events and the equatorial reading of the CM Draconis candidate would both shift.","fun_headline_variants_meta":{"raw":{"variants":["Transit shadows map M dwarf flare latitudes","Occultations pinpoint flare latitude on M dwarfs","TESS data likely hides 3–22 flare occultations","Flare occultations: a map for M dwarf activity","Eclipsing binaries reveal more flare occultations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001093,"raw_usage":{"total_tokens":4570,"prompt_tokens":954,"completion_tokens":3616,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":570,"completion_tokens_details":{"reasoning_tokens":3539}},"tokens_in":570,"tokens_out":3616,"duration_ms":21958,"temperature":1.0,"reasoning_tokens":3539,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:22:25.239901+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A systematic search of all TESS 20-second and 120-second data for the 55 M-dwarf eclipsing binaries defined in the paper: if the number of confirmed type-B occultations after a complete search is statistically consistent with zero while the model predicts at least three, the uniform-flare assumption or the simulated flare morphology templates would be falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"First demonstrated flare occultations as a latitude probe at X-ray wavelengths; the present optical method is a direct extension of this idea."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Detailed the geometric modeling of X-ray flare occultations that this paper's simulation adapts to optical photometry."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the TESS eclipsing binary catalog from which the binary sample for the expected-count estimate is drawn."},{"cited_title":"The Benchmark M Dwarf Eclipsing Binary CM Draconis With TESS: Spots, Flares and Ultra-Precise Parameters","cited_arxiv_id":"2301.10858","evidence_quote":"Identified the 125 CM Draconis flares, including the candidate occultation event that this paper analyzes in detail."}],"review_version":1}