{"id":"375abaca-c4bf-4e81-bff4-094b8f0e904d","arxiv_id":"2507.19698","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A 60-day GBT campaign on the magnetar Swift J1818.0-1607 shows a stable single-peaked pulse profile, slower spin-down, and revised characteristic age of about 2500 years and surface field of about 10^14 G.","lead":"We measured the radio pulses of the magnetar Swift J1818.0-1607 over 60 days with the Green Bank Telescope and found a stable single-peaked profile that differs from the double-peaked shape seen four months earlier. The new timing solution yields a characteristic age near 2500 years and a surface magnetic field near 10^14 G, roughly 2.5 times older and 3 times weaker than earlier estimates, illustrating how short-term spin-down changes distort age and field estimates.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Timing-solution robustness is the load-bearing risk: the phase-connected fit formally includes a large positive F2 term yet drops the first epoch for unquantified timing noise, so the inferred nu-dot, age, and B-field may not be unique.","rationale":"I agree with the reader that this is a solid observational paper and that the cross-frequency pulse-profile comparison weakens the mode-switching claim. However, I identify the timing solution itself as the more load-bearing assumption: the age, B-field, and the spin-down comparison in Figure 6 all rest on the phase connection and on the reality of the positive F2 term. The paper explicitly states that the first epoch could not be included because of timing noise, but does not show that the remaining ToAs yield a unique phase connection, and Table A.3 shows spin-frequency measurements that scatter by far more than their quoted uncertainties (e.g., MJD 59568.709 vs 59568.737 differ by 1.5e-5 Hz, ~20 sigma), indicating unmodeled red noise. A short-span fit with an unexplained, formally highly significant second derivative is exactly the situation where a cycle-count ambiguity or a couple of outlier points can change F1 by order one. The concrete refits I propose would settle this directly. I therefore keep CONDITIONAL, matching the reader's verdict, but for a different primary reason.","tokens_in":18385,"tokens_out":1675,"duration_ms":20981,"concrete_test":"Refit the Table A.1 ToAs with TEMPO2 using a brute-force search over trial phase wraps (e.g., the NX01 method) and repeat the fit (a) including the MJD 59520.76 ToAs, (b) after removing the MJD 59572.63 ToAs, and (c) with F2 forced to zero. If the recovered F1 moves outside the quoted 1-sigma uncertainty, or if the best solutions differ by an integer cycle count, then the phase-connected solution is not unique and the derived age and B-field should be presented only as epoch-dependent estimates.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central result is the phase-connected timing solution in Table 3, from which nu-dot = -4.4855e-12 Hz/s yields tau_c ~ 2500 yr and B_field ~ 1e14 G. The fit includes a positive second derivative F2 = 1.2286e-19 Hz/s^2 quoted with uncertainty 8.5e-20, i.e. formally ~14.5 sigma, which is not physically interpreted and dominates the spin evolution over the 42-day span. Yet the first observing epoch, MJD 59520.76, is excluded from the fit, with only a terse statement that timing noise prevented a phase-connected solution, even though Table A.1 lists five ToAs at that epoch and Table A.3 gives a spin-frequency there. No timing noise amplitude is quantified, no fit including that epoch is shown, and no cycle-count or leave-one-out check is presented. If the phase connection is off by an integer wrap, or the F2 term is driven by a couple of outliers such as the MJD 59572 ToAs (uncertainties 768 and 683 microseconds, the largest in the table), then F1, and hence tau_c and B, change by factors, not percentages. The abstract's headline age and field therefore hinge on a single short-span fit whose uniqueness is not demonstrated. The reader's mode-switching concern is real, but the timing solution is more load-bearing because every derived quantity depends on it.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports GBT 2.0 GHz observations of the radio-loud magnetar Swift J1818.0-1607 taken in November-December 2021, roughly 20 months after its 2020 outburst. Over the 60-day dense set, the integrated pulse profile is stable, single-peaked with a small precursor, in contrast to the double-peaked profile seen at 1.4 GHz in the preceding campaign of Rajwade et al. (2022). The spin-down rate is slower than at the end of that campaign, and the flux density is low and stable. From a phase-connected timing solution the authors derive F1 = -4.4855e-12 Hz/s, a characteristic age of about 2500 years, and a surface dipole field of about 1e14 G, roughly 2.5 times older and three times weaker than the most recent previous estimates. The paper interprets the profile, spin-down, and flux changes as evidence for a likely mode-switching event during the ~100-day gap between the two campaigns.","tokens_in":18773,"tokens_out":6684,"duration_ms":74900,"significance":"If the timing solution is robust, this is a valuable post-outburst measurement of a young radio-loud magnetar, showing that short-term spin-down variability can change the inferred characteristic age and magnetic field by large factors. The paper publishes the full set of times of arrival in Table A.1, uses standard and well-documented software (PRESTO, TEMPO2, DSPSR/PSRCHIVE), and transparently acknowledges that the characteristic age assumes a constant braking index and can be affected by spin-down variability. These are real strengths that make the central measurement reproducible. The main concerns are that the phase-connected solution depends on a large, uninterpreted second frequency derivative and on the exclusion of the first epoch, and that the mode-switching claim compares profiles at different observing frequencies without quantifying frequency-dependent morphology.","major_comments":[{"comment":"The load-bearing timing solution excludes the first epoch at MJD 59520.76 even though Table A.1 lists five ToAs for that epoch, and the only justification in §3.1 is 'significant timing noise,' which is neither quantified nor demonstrated. The fit includes F2 = 1.2286e-19 Hz/s^2 with a quoted 1-sigma uncertainty of 8.5e-20, i.e., formally about 14.5 sigma, yet no physical interpretation, cycle-count check, or leave-one-out test is presented. Over the 42-day span of the fit, the F2 term contributes nearly one full rotation, so a small number of outliers could change F1, and therefore tau_c and B, by factors rather than percentages. In particular, the MJD 59572.63 ToAs have uncertainties of 768 and 683 microseconds, the largest in Table A.1, and the MJD 59568.73670 spin-frequency entry in Table A.3 differs from the MJD 59568.70900 entry by about 1.5e-5 Hz, far outside either error bar. The authors should show a fit that includes the first epoch, quantify the timing noise amplitude, and demonstrate that F1 and F2 are stable under removal of individual epochs or ToAs.","section":"§3.1, Table 3; Table A.1; Table A.3"},{"comment":"The claimed mode switch between the Rajwade et al. (2022) campaign and this campaign rests on comparing a double-peaked 1.4 GHz profile with a single-peaked 2.0 GHz profile, but the paper does not model or quantify how the pulse profile of Swift J1818.0-1607 depends on observing frequency. Since the same physical emission state can appear different at different frequencies, as the paper itself notes for flux density in §4.1, the profile difference alone does not uniquely imply a mode switch. The authors should either quantify frequency-dependent profile morphology using contemporaneous multi-frequency data (e.g., Huang et al. 2021, Bansal et al. 2023) or soften the mode-switching claim to an explicitly conditional interpretation.","section":"§4.1, Figures 3 and 4"},{"comment":"The inferred characteristic age and magnetic field use only F1 from Table 3, while F2 is formally highly significant. The paper should state how tau_c and B would change if F2 were absorbed into the error budget or omitted. It should also avoid the statement that the magnetar has 'effectively aged by about 2000 years over the past two years,' which conflates a change in the measured spin-down rate with a physical age increase; the authors do later acknowledge the caveat, but Table 5 and Figure 7 treat the successive P-Pdot points as directly comparable measurements.","section":"§4.4, Table 5"}],"minor_comments":[{"comment":"The abstract contains the typo 'with with' in the second sentence, and the text has several other typos ('obsevations', 'surpising', 'Febraury', 'monotic', 'adminstrated', 'Throghout') that should be corrected in the proof stage.","section":"Abstract"},{"comment":"Table 2 lists the first epoch at MJD 59520.76, but Table 3 states the timing-solution date range begins at MJD 59536.713; the text in §3.1 should explicitly note this discrepancy and clarify how the first epoch was used for the spin-frequency comparison in Figure 6.","section":"Table 3 and Table 2"},{"comment":"The caption of Figure 7 says the inset shows four data points, while Table 5 lists five campaign estimates; the caption should clarify whether the initial Champion et al. (2020a) point is omitted from the inset.","section":"Figure 7 caption"},{"comment":"The fixed RA, Dec, and DM values are taken from the literature without a sensitivity test; a sentence reporting how much the timing solution changes when these are varied within their quoted uncertainties would be helpful.","section":"§3.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is well within the scope of ApJ and the central measurement is reproducible because the full ToA set is published. The main risk is that the phase-connected solution, and hence the headline age and magnetic field, may not be unique over this short span; robustness tests of the type requested in the major comments should be feasible without new observations. No concerns about novelty or citation practice."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a straightforward observational paper, and most of it is solid. The genuinely new piece is a phase-connected timing solution for Swift J1818.0-1607 covering MJD 59536.7-59578.8 at 2 GHz, with the full ToA list published in the appendix. That is a concrete, useful data point, and the paper makes a fair point that short-term spin-down variability can change the inferred characteristic age by a factor of 2.5 within two years.\n\nThe analysis uses standard, well-documented tools (TEMPO2, PRESTO, PSRCHIVE), and the authors are honest that the derived tau_c and B-field reflect the current spin-down rate rather than intrinsic changes. Placing the result on the P-Pdot diagram and comparing with earlier campaigns is well handled.\n\nThe main soft spot is the timing solution itself. The fit includes a large positive F2 = 1.23e-19 Hz/s^2, formally ~14.5 sigma, which is not physically interpreted. More concerning, the first epoch, MJD 59520.76, is excluded from the fit with only a terse 'significant timing noise' explanation. A quick extrapolation of the Table 3 solution back to that epoch gives a ~4-sigma difference from the measured spin frequency in Table A.3. That does not necessarily break the solution, but it means the derived F1, and hence tau_c and B, are sensitive to the chosen span and possibly to red noise absorbed by F2. The referee should ask the authors to show a fit including that epoch, quantify the timing noise, and run leave-one-out checks on F1 and F2. This is a robustness issue, not a fatal flaw, and the published ToAs make it checkable.\n\nThe mode-switching claim is the other weak spot. The profile comparison is between a double-peaked profile at 1.4 GHz from Rajwade et al. and a single-peaked profile at 2.0 GHz here. Frequency-dependent profile morphology could produce exactly this difference. The authors are aware of this and lean on the accompanying changes in spin-down rate and flux density, which makes the claim plausible but not conclusive. A near-simultaneous multi-frequency observation would have settled it.\n\nMinor typos and some overlong discussion do not affect the substance. This paper is for the magnetar timing community and deserves a serious referee, but the referee should push on the timing robustness before accepting the headline numbers. I would not desk-reject it.","headline":"Useful and reproducible timing update for Swift J1818.0-1607, but the headline age/field values rest on a short span that excludes one epoch, and the mode-switch claim crosses frequencies.","tokens_in":19241,"tokens_out":3812,"would_cite":true,"duration_ms":45312,"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":"A phase-connected timing solution from a 60-day Green Bank Telescope campaign places the magnetar Swift J1818.0-1607 at a characteristic age of about 2,500 years, roughly 2.5 times older than the previous estimate.","keywords":["magnetars","Swift J1818.0-1607","radio pulsars","timing solution","spin-down","mode switching","pulse profile","neutron stars"],"falsifier":"Observe Swift J1818.0-1607 simultaneously at 1.4 GHz and 2.0 GHz during a single epoch: if the 1.4 GHz profile is double-peaked while the 2.0 GHz profile is single-peaked, the mode-switch interpretation would be weakened in favor of frequency-dependent profile morphology.","tokens_in":18165,"feed_emoji":"🕰️","tokens_out":6170,"duration_ms":60285,"temperature":0.7,"pith_summary":"The paper uses Green Bank Telescope observations at 2.0 GHz, taken about 20 months after Swift J1818.0-1607's 2020 outburst, to build a phase-connected timing solution spanning roughly 60 days. From that solution it infers a characteristic age of about 2,500 years and a surface dipole magnetic field of about $1\\times 10^{14}$ G, both substantially revised from earlier estimates. The paper also argues that the magnetar switched emission modes in the ~100-day gap before these observations, because the pulse profile changed from double-peaked to a stable single narrow peak with a precursor, the spin-down rate slowed, and the flux stayed low. The broader point is that for a young magnetar with highly variable spin-down, inferred age and field strength depend strongly on when the measurement is made rather than on the magnetar's intrinsic properties.","feed_headline":"Radio campaign revises young magnetar's age to 2,500 years","feed_subtitle":"Two months of 2 GHz observations find a stable single-peaked profile and a slower spin-down, pointing to a mode switch.","key_machinery":"The load-bearing object is the phase-connected timing solution built from arrival times extracted with PRESTO and fitted with TEMPO2 using the phase-connection technique of Freire & Ridolfi (2018). The inference chain then runs through two standard identities: the characteristic age $\\tau_c = P/(2\\dot{P})$ assuming a braking index of $n = 3$, and the spin-down-inferred surface dipole field $B = 3.2\\times 10^{19}\\sqrt{P\\dot{P}}$ G. These identities turn the measured spin frequency and its derivative into age and field estimates, which is why short-term changes in the spin-down rate directly change the inferred properties. The mode-switching argument is carried by comparing the pulse profile, spin-down rate, and flux density across epochs.","core_discovery":"Over MJD 59536.7–59578.8 the authors derive a phase-connected timing solution for Swift J1818.0-1607 with spin frequency $F_0 = 0.7326046915(5)$ Hz, spin-down $F_1 = -4.4855(11)\\times 10^{-12}$ Hz s$^{-1}$, and a small second derivative, holding the dispersion measure fixed at $710\\pm1$ pc cm$^{-3}$. Using the standard magnetar formulas $\\tau_c = P/(2\\dot{P})$ and $B = 3.2\\times 10^{19}\\sqrt{P\\dot{P}}$ G, they infer $\\tau_c \\sim 2500$ years and $B \\sim 1\\times 10^{14}$ G, about 2.5 times older and nearly three times weaker than the most recent published values. Throughout the campaign the integrated pulse profile remained stable: a single narrow peak with $W_{50}\\sim 18$–23 ms, a small precursor component, no postcursor, a flat radio spectrum ($\\alpha \\gtrsim -1$), and flux densities of roughly 0.2–0.3 mJy. Comparing these with the preceding campaign's double-peaked profile, faster spin-down, and declining flux, the paper concludes that a mode-switching event likely occurred during the ~100-day gap (MJD 59426–59520) and that the source had returned to a state resembling Mode 3 rather than Mode 4.","pith_inferences":["If characteristic ages of young magnetars fluctuate this much with spin-down state, population-level age and magnetic-field estimates drawn from P–$\\dot{P}$ diagrams for recently outbursting magnetars carry systematic uncertainties that independent age indicators, such as supernova remnant expansion or kinematic measurements, would be needed to break.","The stability of the single-peaked profile across the full 60 days, in contrast to the minute-timescale mode switching seen shortly after the outburst, suggests the magnetosphere may settle into a more stable configuration during reduced activity; a longer multi-frequency campaign could test whether that stability persists.","A natural extension is to measure the braking index once the post-outburst relaxation finishes; if it is less than 3, as seen in some magnetars, the true age would be even larger than the characteristic age.","The persistently flat spectrum at $\\alpha \\gtrsim -1$ supports the idea that Swift J1818.0-1607's radio emission is becoming more rotation-powered-pulsar-like over time, potentially making it a bridge between the magnetar and radio pulsar populations."],"forward_implications":["The magnetar's inferred characteristic age and surface dipole field are epoch-dependent: the same object has been assigned $\\tau_c$ from roughly 265 to 2500 years and $B$ from about $3.4\\times 10^{14}$ G to $1\\times 10^{14}$ G depending on when the spin-down was measured.","A mode-switching episode occurred between the last published campaign and this one, during the ~100-day gap, leaving the magnetar in a state resembling the single-peaked, slower-spinning Mode 3.","The spin-down rate measured here is an order of magnitude slower than during Mode 4, reinforcing that post-outburst magnetar spin-down is dominated by short-term fluctuations rather than steady electromagnetic braking.","Continued monitoring will determine whether Swift J1818.0-1607 follows the fading and reactivation cycles seen in PSR J1622-4950 and XTE J1810-197 or becomes more sporadic like 1E 1547.0-5408."],"supporting_citations":[{"why":"Supplies the prior timing solution, the Mode 3 and Mode 4 definitions, and the 1.4 GHz double-peaked profile that this campaign is directly compared against.","marker":"Rajwade et al. (2022)"},{"why":"Provides the first spin-down measurements and the initial characteristic age and magnetic field estimates that this work revises.","marker":"Champion et al. (2020b)"},{"why":"Contributes earlier multi-frequency spin-down and spectral-index measurements used to track the source's evolution and to compare flux densities.","marker":"Huang et al. (2021)"},{"why":"Provides the phase-connection technique used to construct the timing solution from the extracted arrival times.","marker":"Freire & Ridolfi (2018)"},{"why":"Documents the two emission modes below 4 GHz seen during the 2020 outburst, establishing the magnetar's mode-switching behavior.","marker":"Lower et al. (2021)"},{"why":"Sets the magnetar population context and the standard formulas for characteristic age and spin-down-inferred magnetic field.","marker":"Kaspi & Beloborodov (2017)"}],"fun_headline_variants":["Magnetar's age jumps to 2,500 years in new timing analysis","Swift J1818: mode switch and a more mature magnetar","Radio campaign reveals magnetar's true age and a mode shift","Magnetar's spin-down change points to mode switching","Young magnetar looks older after 60 days of radio watching"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The mode-switching claim assumes that the difference between the earlier double-peaked 1.4 GHz profile and this campaign's single-peaked 2.0 GHz profile reflects a genuine change in emission state, rather than the same state looking different at a different observing frequency.","fun_headline_variants_meta":{"raw":{"variants":["Magnetar's age jumps to 2,500 years in new timing analysis","Swift J1818: mode switch and a more mature magnetar","Radio campaign reveals magnetar's true age and a mode shift","Magnetar's spin-down change points to mode switching","Young magnetar looks older after 60 days of radio watching"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000397,"raw_usage":{"total_tokens":2169,"prompt_tokens":1129,"completion_tokens":1040,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":745,"completion_tokens_details":{"reasoning_tokens":950}},"tokens_in":745,"tokens_out":1040,"duration_ms":10139,"temperature":1.0,"reasoning_tokens":950,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T14:07:44.526006+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe Swift J1818.0-1607 simultaneously at 1.4 GHz and 2.0 GHz during a single epoch: if the 1.4 GHz profile is double-peaked while the 2.0 GHz profile is single-peaked, the mode-switch interpretation would be weakened in favor of frequency-dependent profile morphology.","supporting_citations":[{"cited_title":"2021, MNRAS, 505, 1311, doi: 10.1093/mnras/stab1362 —","cited_arxiv_id":null,"evidence_quote":"Contributes earlier multi-frequency spin-down and spectral-index measurements used to track the source's evolution and to compare flux densities."}],"review_version":1}