{"id":"e2f7eef8-4f3a-4341-9d8c-6eb7ea381c79","arxiv_id":"1908.08583","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Using near-infrared spectra and Zeeman-broadening models, the authors derive a 1.36 kG surface magnetic field for the class I protostar V347 Aur, the first such measurement for this object, and a 2.5 kG field for BP Tau.","lead":"This paper measures the magnetic fields of two young stars using near-infrared spectra from the iSHELL instrument on NASA's IRTF telescope. It reports the first magnetic field measurement for the class I protostar V347 Aur and confirms a previously measured field for BP Tau.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Solar-calibrated line parameters may absorb model defects and transfer poorly to the cool, low-gravity V347 Aur, so the 1.36 kG detection needs an independent test and its error bars must include the 0.31 kG systematic floor.","rationale":"The reader's weakest assumption correctly identifies the solar-calibrated line-list modifications (Appendix B) as the most fragile link, and the paper's own cautionary sentence confirms the risk. I agree with the reader that this is an addressable concern rather than a fatal flaw: the BP Tau recovery of 2.5 kG against literature values and the worse nonmagnetic fits provide independent support for the method's qualitative magnetic sensitivity. However, the transferability gap is real and specific: no standard star in the sample has both cool temperature and low gravity, so the VdW and log(gf) adjustments tuned on the Sun (5778 K, log g = 4.43) are being applied unvalidated at Teff = 3233 K, log g = 3.25, where pressure broadening and molecular opacities differ substantially. The mechanism by which this could bias B is concrete: several heavily modified lines (Ti I 2.2239, 2.2317; Ca I 2.2614, 2.2657) are precisely the magnetically sensitive lines that carry the Zeeman signal, and their VdW constants were changed by up to 1.2 dex. The second issue, omission of the 0.31 kG systematic floor from the B error bars, is an internal-consistency problem: the floor is quoted in the same paper (Section 4.5) and should be propagated into the headline measurement, not merely used as a pass-fail threshold. With the floor added, 1.36 ± 0.31 kG is still a detection, but a far more modest one. Since the reader already conditioned acceptance on these points and my analysis does not reveal a deeper flaw, the verdict remains conditional. The concrete test I propose (refitting with the unmodified line list and with magnetically insensitive lines) would decisively separate a genuine Zeeman signal from a line-list-induced pseudo-broadening and is simple enough to run with the same MCMC machinery.","tokens_in":26909,"tokens_out":2388,"duration_ms":27618,"concrete_test":"Re-derive Teff, log g, and <B> for V347 Aur with three modifications in sequence: (1) use the unmodified VALD3 line list (keeping the fit's free microturbulence and v sin i), (2) restrict the comparison to the CO (2-0) and other magnetically insensitive lines to fix Teff/log g independently, and (3) add the 0.31 kG null-field floor in quadrature to the B uncertainty. If the unmodified list shifts <B> toward the 0.31 kG floor, or if the insensitive-line fit gives Teff/log g that differ by more than the quoted errors, the solar-calibrated parameter modifications are likely absorbing model defects and the V347 Aur detection is not yet validated. A stronger alternative is to fit the same spectra with an independent magnetic radiative transfer code and independently measured line parameters (e.g., COSSAM or Zeeman2 with APOGEE line lists).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central new result is the first magnetic field measurement of V347 Aur (1.36 kG), inferred from Zeeman broadening in K-band lines fitted with MoogStokes. That inference rests on 26 solar-calibrated log(gf) and van der Waals modifications (Appendix B, Table 8). As the authors themselves warn, these adjustments 'could be hiding defects associated with the stellar atmospheric models, the radiative transfer code, or even the measured instrument spectral profile.' The calibration star is the Sun (Teff ~ 5778 K, log g = 4.43); V347 Aur is cool (3233 K) and low-gravity (log g = 3.25). The modified lines include exactly the magnetically sensitive Ti I and Ca I lines that carry the Zeeman signal, and several VdW constants were changed by 0.3–1.2 dex. In a cool, low-gravity atmosphere, van der Waals broadening is much stronger than in the Sun, so an absorbed VdW error would create a Teff/log g/broadening bias that mimics or cancels real Zeeman width. The paper does not validate the transfer: the standard-star sample covers the Teff range but not the combination of cool Teff with log g ~ 3.2, and the null-field floor of 0.31 kG was measured on giants and the Sun, not at YSO parameters. Separately, the quoted error bars for V347 Aur's B (±0.05–0.06 kG) come only from MCMC posteriors; the 0.31 kG systematic floor is stated as a detection threshold but never added to the uncertainty. If that floor applies at V347 Aur's parameters, the true uncertainty is ~0.31 kG, reducing the detection significance from ~20 sigma to ~4 sigma. These two issues jointly weaken the claim that 1.36 kG is robustly above the detection limit.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a method for deriving stellar parameters and surface magnetic fields of young stars by fitting iSHELL K-band spectra with the MoogStokes Zeeman-broadening code. The method is validated on nine main- and post-main-sequence stars plus the Sun, giving a temperature scatter of 91 K and a gravity scatter of 0.14 dex, and a null-field detection floor of 0.31 kG. Applied to BP Tau, the method recovers B = 2.5 kG, consistent with previous literature. For the class I source V347 Aur, the paper reports the first magnetic field measurement (B = 1.36 kG), along with Teff = 3233 K and log g = 3.25, and shows in a nonmagnetic control test that excluding B degrades the fit to the magnetically sensitive lines while increasing v sin(i) and microturbulence. The paper also discusses the influence of starspots on masses and ages derived from pre-main-sequence evolutionary tracks.","tokens_in":27227,"tokens_out":6417,"duration_ms":62050,"significance":"If the calibration concerns are resolved, this is a valuable contribution: it extends Zeeman-broadening measurements to the K band with broad wavelength coverage, provides an externally consistent BP Tau benchmark, and delivers a new class I magnetic field measurement for V347 Aur. The empirical instrumental profile, the tabulated line-list modifications, the standard-star validation, and the nonmagnetic control test are explicit and reproducible steps that give the method a solid empirical base. The main risk is systematic rather than statistical: the solar-calibrated line parameters and the 0.31 kG null-field floor must be shown to transfer to V347 Aur's cool, low-gravity parameters before the reported 1.36 kG can be taken at face value.","major_comments":[{"comment":"The quoted uncertainties for V347 Aur's magnetic field (B = 1.36+0.06-0.05 kG in Table 5) are posterior percentiles only and do not include the 0.31 kG null-field floor established in Section 4.5. Because that floor is the strongest field recovered on stars where a null field was expected, it should be added in quadrature as a systematic uncertainty, or the authors should demonstrate that it does not apply at V347 Aur's Teff = 3233 K and log g = 3.25. Even after adding the floor the detection remains significant, but the reported precision would be materially different and the paper should present the combined error budget.","section":"Section 4.5 and Table 5"},{"comment":"The V347 Aur detection relies on the Zeeman broadening of Ti I and Ca I lines whose log(gf) and van der Waals constants were adjusted on the Sun (Table 8). Some VdW changes are large, up to about 1.2 dex (e.g., Ti I 22627.394 Å), and because van der Waals broadening is stronger in a cool, low-gravity atmosphere, an error absorbed into these constants at solar conditions could mimic or cancel the Zeeman width at V347 Aur's parameters. The authors themselves caution in Appendix B that the adjustments 'could be hiding defects associated with the stellar atmospheric models, the radiative transfer code, or even the measured instrument spectral profile.' I ask for a concrete external check: fit a standard star near V347 Aur's Teff and log g that was not used in the solar calibration, or explicitly propagate a plausible range of VdW errors into the B measurement, to quantify the resulting systematic uncertainty.","section":"Appendix B, Table 8"},{"comment":"The statement that V347 Aur's parameters 'are contained within the range of stellar parameters we investigated in Section 4' is inaccurate for temperature: the coolest standard star in Table 4 is GJ 436 at 3401 K, while V347 Aur is 3233 K, and no standard star combines Teff below about 3400 K with log g around 3.2. The null-field floor of 0.31 kG in Section 4.5 was measured on warmer giants and on the Sun. The transfer of both the adjusted line list and the detection limit to V347 Aur is therefore an extrapolation and should be identified as such, or supported by additional anchor points.","section":"Section 5.2 with Section 4"}],"minor_comments":[{"comment":"The text says 'we adopted, for simplicity, a solar composition in all our models,' but Section 3.3 and Table 4 treat [M/H] as a free parameter; please clarify that 'solar composition' refers to element abundance ratios rather than a fixed metallicity.","section":"Section 3.1"},{"comment":"There is a missing space before 'λ' in the sentence below Equation (1), and the wavelength unit is given as microns there while other parts of the paper use Ångströms; please harmonize the notation.","section":"Equation (1) and surrounding text"},{"comment":"Please state the units and sign convention for the 'Waals' column and provide the origin of the default values, so that the modifications are fully reproducible.","section":"Table 8"},{"comment":"The sentence 'if the covering factor ... iw known' contains a typo; it should read 'is known.'","section":"Section 6.1"},{"comment":"The note says that actual errors are significantly larger than the tabulated formal uncertainties, but the paper does not state the final combined uncertainty for each standard star; consider adding a column or equation showing how sigma_Teff = 91 K, sigma_logg = 0.14, and the 0.31 kG magnetic floor combine with the MCMC errors.","section":"Table 4 note"}],"recommendation":"major_revision","confidential_remarks":"This is a technically solid paper that fits the journal's scope. The central new result (V347 Aur's field) is plausible, but the error budget and the transfer of the solar-calibrated line list to cool, low-gravity parameters need to be addressed in revision; these issues are fixable within the manuscript's scope, so major revision is appropriate rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the genuinely new result is the first field measurement of V347 Aur, 1.36 kG, for a class I protostar. I think the detection is real, but the quoted error bars are smaller than the actual systematic uncertainty, and the solar-calibrated line list is the main thing to worry about.\n\nThe paper does real work. It measures iSHELL's spectral profile, recalibrates 26 K-band lines against solar spectra, validates on nine standard stars with 91 K and 0.14 dex scatter, recovers BP Tau's known 2.5 kG field, and runs a nonmagnetic control that forces v sin i and microturbulence to absorb some of the Zeeman broadening. That is a careful, honest setup. The authors flag the line-list caveat themselves in Appendix B in plain terms, which I credit. The citation pattern looks fine and engages the relevant YSO Zeeman literature.\n\nThe two soft spots are real and should be fixable. First, the 0.31 kG null-field floor is stated as a detection limit but never added to the B uncertainties. For V347 Aur that means the true uncertainty is closer to 0.3 kG than 0.05 kG, so the detection is roughly 4 sigma, not 20 sigma. Still detected, but the abstract's error bars overstate precision. Second, the line-list adjustments are calibrated on the Sun and include large van der Waals changes on exactly the Ti I and Ca I lines that carry the Zeeman signal. The standard-star sample does not cover the cool, low-gravity combination of V347 Aur (Teff ~ 3230 K, log g ~ 3.25), where van der Waals broadening is much stronger, so there is no direct check that the recalibration transfers. The nonmagnetic control helps, but it uses the same line list. This is a concern, not a fatal flaw; 1.36 kG is plausible and still above the floor even with the floor added in quadrature.\n\nMinor point: the starspot interpretation for BP Tau's cool NIR temperature is reasonable but not uniquely established, and the authors say as much. The evolutionary-model discussion is illustrative and appropriately hedged.\n\nThis paper is for observers working on YSO magnetic fields and near-infrared parameter pipelines. It deserves a serious referee, though conditional acceptance is the right outcome: the B uncertainty should include the 0.31 kG systematic floor, and the line-list transfer should either be justified with a cool, low-gravity benchmark or discussed more cautiously.","headline":"The first magnetic field measurement of the class I protostar V347 Aur is a real result, but the quoted error bars understate the systematic floor and the solar-calibrated line list needs an independent check at cool, low-gravity parameters.","tokens_in":27851,"tokens_out":2750,"would_cite":true,"duration_ms":28840,"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":"Fitting high-resolution near-infrared spectra with a magnetized radiative-transfer model recovers stellar parameters to within 91 K and measures the first magnetic field of the protostar V347 Aur at 1.36 kG.","keywords":["young stellar objects","stellar magnetic fields","Zeeman effect","infrared spectroscopy","pre-main-sequence evolution","starspots","T Tauri stars","protostars"],"falsifier":"Fit the V347 Aur spectrum with a different atmospheric model grid or with line parameters that were not solar-calibrated; if the recovered field drops to the 0.31 kG detection limit, the 1.36 kG detection was an artifact of the calibration assumptions. Alternatively, obtain an independent spectropolarimetric or optical Zeeman measurement of V347 Aur; disagreement beyond the quoted uncertainties would falsify the transferability of the method.","tokens_in":26634,"feed_emoji":"🧲","tokens_out":6486,"duration_ms":61479,"temperature":0.7,"pith_summary":"The paper establishes that high-resolution near-infrared K-band spectra can be modeled with a magnetized radiative-transfer code to recover the fundamental parameters of young stars in one self-consistent fit. The method is validated on nine main-sequence and giant stars plus the Sun, giving temperature uncertainties of 91 K and gravity uncertainties of 0.14 dex, with a magnetic detection limit of 0.31 kG. Applied to the well-studied class II star BP Tau, it recovers a surface field of 2.5 kG in agreement with earlier Zeeman studies, while its lower temperature points to starspots dominating the infrared. Applied to the class I protostar V347 Aur, it yields the first measurement of its magnetic field, 1.36 kG, with log g = 3.25. If correct, these results extend magnetic-field measurements to the earliest observable protostellar stage and show that starspot contamination must be handled when deriving masses and ages.","feed_headline":"First magnetic field measured for protostar V347 Aur: 1.36 kG","feed_subtitle":"K-band Zeeman fitting also reproduces BP Tau's known 2.5 kG field and points to starspots in cooler infrared temperatures.","key_machinery":"The Zeeman effect is the load-bearing mechanism: in a magnetic field, a spectral line's components separate by an amount proportional to lambda squared times the Lande factor times the field strength, so line widths carry the magnetic signal. Six selected K-band windows include lines with high magnetic sensitivity, such as Ti lines, and lines with low sensitivity, such as CO rovibrational transitions, and all six are fitted simultaneously. The fitting chain combines a plane-parallel LTE magnetic synthesis code, model atmospheres, an empirically measured instrument profile, and a Markov chain Monte Carlo parameter search to recover temperature, gravity, field strength, rotation, microturbulence, veiling, and CO abundance together.","core_discovery":"The paper claims that the broad wavelength coverage of a high-resolution K-band spectrograph, combined with a magnetic radiative-transfer code that includes Zeeman broadening, is sufficient to determine the fundamental parameters of young stars in one consistent fit. Validated against nine main- and post-main-sequence stars and the Sun, the method returns effective temperatures with a scatter of 91 K and surface gravities with a scatter of 0.14 dex, and it sets a detection limit of 0.31 kG. Applied to the class II star BP Tau, it recovers a surface magnetic field of 2.5 kG, consistent with earlier optical Zeeman studies, but a temperature about 400 K cooler, which the authors attribute to starspot emission dominating the K band. Applied to the class I protostar V347 Aur, it yields the first measurement of that object's field, 1.36 kG, together with log g = 3.25, and it argues that dropping the magnetic term from the fit inflates rotation and turbulence estimates while degrading the fit.","pith_inferences":["If starspots systematically depress infrared temperatures, then masses derived from infrared-only fits to heavily spotted young stars are biased low unless spot filling factors are modeled explicitly.","The method's detection limit suggests that a survey of class I protostars using the same six spectral windows could map the distribution of fields at the earliest ages and test whether V347 Aur's relatively weak 1.36 kG field is typical.","A direct extension would be to fit the same spectra with two-temperature spotted models; the difference between the single-temperature and two-temperature recovered fields would indicate how much of the Zeeman broadening is masked by spot contrasts."],"forward_implications":["Magnetic field measurements can be extended to class I protostars, the youngest pre-main-sequence phase, with a uniform parameter fit rather than assumed temperatures and gravities.","Because Zeeman broadening scales as wavelength squared, the near-infrared approach can detect weaker fields than optical studies, making it a route to a larger protostellar magnetic-field sample.","The 400 K difference between optical and infrared temperatures of BP Tau means single-temperature masses and ages of spotted young stars are uncertain by up to roughly a factor of two, depending on the evolutionary model.","Nonmagnetic fits to strongly magnetic young stars overestimate projected rotation and microturbulence, so magnetic terms cannot be ignored when measuring rotation in such objects."],"supporting_citations":[{"why":"Supplies the magnetic radiative transfer code that introduces Zeeman broadening into the spectral synthesis.","marker":"Deen 2013"},{"why":"Provides the MARCS model atmospheres used as input for all synthetic spectra.","marker":"Gustafsson et al. 2008"},{"why":"Supplies the VALD3 atomic line parameters that the paper calibrates against solar observations.","marker":"Ryabchikova et al. 2015"},{"why":"Provides the prior BP Tau magnetic field measurement that the new method must reproduce.","marker":"Johns-Krull et al. 1999"},{"why":"Supplies the second benchmark BP Tau field measurement and the earlier Zeeman modeling approach.","marker":"Johns-Krull 2007"},{"why":"Provides the prior spectral classification and low-gravity indication for V347 Aur, the main new target.","marker":"Connelley & Greene 2010"},{"why":"Describes the iSHELL spectrograph whose large K-band coverage enables the multi-line fitting.","marker":"Rayner et al. 2016"},{"why":"Supplies the MCMC sampler used to estimate stellar parameters and their uncertainties.","marker":"Foreman-Mackey et al. 2013"},{"why":"Provides the nonmagnetic pre-main-sequence evolutionary models used for mass and age comparisons.","marker":"Baraffe et al. 2015"},{"why":"Provides magnetic evolutionary models that show how convection suppression changes inferred masses and ages.","marker":"Feiden 2016"}],"fun_headline_variants":["First magnetic field for protostar V347 Aur: 1.36 kG","K-band Zeeman fit measures protostar magnetic field for first time","Starspots explain cooler infrared temperatures in young stars","Zeeman fitting in K-band yields stellar parameters and magnetic fields","New method measures stellar magnetic fields from K-band spectra"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The solar-calibrated changes to 26 atomic line parameters are physically real and carry over to the cooler, lower-gravity photosphere of V347 Aur, rather than absorbing errors in the atmospheric models, radiative transfer, or instrument profile.","fun_headline_variants_meta":{"raw":{"variants":["First magnetic field for protostar V347 Aur: 1.36 kG","K-band Zeeman fit measures protostar magnetic field for first time","Starspots explain cooler infrared temperatures in young stars","Zeeman fitting in K-band yields stellar parameters and magnetic fields","New method measures stellar magnetic fields from K-band spectra"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001139,"raw_usage":{"total_tokens":4804,"prompt_tokens":1095,"completion_tokens":3709,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":711,"completion_tokens_details":{"reasoning_tokens":3621}},"tokens_in":711,"tokens_out":3709,"duration_ms":23350,"temperature":1.0,"reasoning_tokens":3621,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:35:22.237976+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fit the V347 Aur spectrum with a different atmospheric model grid or with line parameters that were not solar-calibrated; if the recovered field drops to the 0.31 kG detection limit, the 1.36 kG detection was an artifact of the calibration assumptions. Alternatively, obtain an independent spectropolarimetric or optical Zeeman measurement of V347 Aur; disagreement beyond the quoted uncertainties would falsify the transferability of the method.","supporting_citations":[{"cited_title":"P.\\ 2013, , 146, 51","cited_arxiv_id":null,"evidence_quote":"Supplies the magnetic radiative transfer code that introduces Zeeman broadening into the spectral synthesis."},{"cited_title":"M., Valenti, J","cited_arxiv_id":null,"evidence_quote":"Provides the prior BP Tau magnetic field measurement that the new method must reproduce."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Describes the iSHELL spectrograph whose large K-band coverage enables the multi-line fitting."},{"cited_title":"A.\\ 2016, , 593, A99 [F16]","cited_arxiv_id":null,"evidence_quote":"Provides magnetic evolutionary models that show how convection suppression changes inferred masses and ages."}],"review_version":1}