{"id":"1cdf181b-7a2d-40c4-b153-e3d48a8096ea","arxiv_id":"2411.12133","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In 89 archival spectra of accreting brown dwarfs and very low-mass stars, 55 are classified as flow-dominated and 15 as shock-dominated using Balmer line ratios, with no clear boundary in physical parameters.","lead":"This paper compares two competing models of how gas falls onto brown dwarfs and very low-mass stars, using hydrogen line ratios from 76 archived objects. Most sources look flow-dominated, but about 15 of 89 observations favor a shock origin, and the choice can change the estimated mass accretion rate by severalfold for the lowest-mass accretors.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The ad hoc 10% error floor in §4.1 directly sets the classification counts; removing it could shift the 15/55 census and the shock-sample accretion-rate comparison.","rationale":"I read the paper as a genuinely useful empirical application of two published emission models to archival Balmer data. The models are prior work, the data are external, and the qualitative conclusion that both shock and flow mechanisms are present is plausible. The reader's CONDITIONAL verdict identified several addressable weaknesses. I agree with the reader that the grid coarseness of the shock model (§3, Figure 1) deserves scrutiny, but I think the single most load-bearing issue is the 10% error floor in §4.1, because it is an explicit, unvalidated parameter injected directly into the classification that produces the central census. The floor is applied to a large subset of the points, including some with extremely small formal errors, and the paper's stated purpose for it is to change classifications ('prevent... neither'). Without a demonstration that the floor does not alter the counts, the 15/55 numbers are not yet trustworthy. The concrete test above is cheap and decisive. I therefore do not change the reader's verdict: the paper should be CONDITIONAL until the census is shown to be stable under reasonable error prescriptions.","tokens_in":31875,"tokens_out":5494,"duration_ms":57303,"concrete_test":"Re-run the §4.1 classification of all 89 points under three error prescriptions: (1) reported errors only, excluding points with no reported errors; (2) a 5% floor instead of 10%; (3) the published 10% floor. Compare the shock/flow/both/neither counts and, crucially, the membership of the 15-point shock sample. If the shock count or the identities of the shock points change materially (more than ~3 points), the census and the §5.2 accretion-rate comparison are not robust to the ad hoc error floor.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The classification in §4.1 is made by asking whether the 1σ error ellipse around each of the 89 data points overlaps the shock locus, the flow locus, both, or neither. For every point without reported errors, and for every point whose reported errors are smaller than 10%, the authors replace the errors with a 10% floor 'to prevent classifying points near the model loci as neither.' This is an arbitrary intervention that inflates the error ellipses of many points with high-SNR measurements (e.g., J08440915–7833457, with formal ratios of 1.32 ± 0.00 and 2.48 ± 0.01, becomes ±0.13 and ±0.25). Larger error ellipses are more likely to overlap one or both model loci, so the floor systematically reduces the 'neither' category and can move points among 'shock,' 'flow,' and 'both.' Because the paper's strongest claim is the census (15 shock vs. 55 flow) and the subsequent accretion-rate comparison for the 15 shock points, an unvalidated floor can directly change the headline numbers. No sensitivity test (e.g., 5% floor, or using only reported errors) is provided, so the central quantitative result is not robust to this choice.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper applies two published emission-line models (the accretion shock model of Aoyama et al. 2018 and the accretion flow model of Kwan & Fischer 2011) to archival Balmer line ratios (Hβ/Hγ versus Hβ/H8) for 89 data points from 76 accreting brown dwarfs and very low-mass stars. It classifies 15 points as shock-dominated, 55 as flow-dominated, 13 as consistent with both, and 6 as consistent with neither, based on 1σ error-ellipse overlap with the model loci. For the 15 shock-dominated points, the shock model yields mass accretion rates up to an order of magnitude higher than the extrapolated stellar scaling at the low-mass end. The paper also reports new Subaru/HDS spectroscopy of three targets, examines the absence of a clear boundary in physical parameters between the two categories, and documents epoch-to-epoch variability of the line-ratio classification for five objects.","tokens_in":32100,"tokens_out":5968,"duration_ms":58741,"significance":"The central claim, if robust, would provide a practical line-ratio diagnostic for identifying the dominant accretion-emission mechanism in substellar objects, which is directly relevant to interpreting accretion-rate estimates for brown dwarfs and giant planets. The use of a relatively large archival sample and the explicit comparison of two physical models is valuable, and the inclusion of new high-resolution spectroscopy adds useful data. The paper is transparent about the ad hoc 10% error floor and the extrapolation of the chromospheric-activity relation, but the lack of sensitivity analysis makes the headline classification counts and the subsequent accretion-rate comparison insufficiently robust. The strength of the paper is its clear formulation of the diagnostic; the main weakness is the unvalidated choices in the classification procedure, which are load-bearing for the quantitative claims.","major_comments":[{"comment":"The 10% error floor assigned to data points without reported errors or with errors below 10% is an arbitrary intervention that directly inflates the error ellipses used for classification. For example, the 2021-04-27 epoch of J08440915−7833457 in Table 2 has reported errors of 0.00 and 0.01 in Hβ/Hγ and Hβ/H8, respectively, which become ±0.13 and ±0.25 under the floor. Because the overlap classification depends on the ellipse size, the floor systematically reduces the 'neither' category and can shift points between shock, flow, and both categories. The paper provides no sensitivity test (e.g., 5% floor, no floor, or a floor based on independently estimated flux-calibration uncertainties). Since the central census (15 shock vs 55 flow) and the accretion-rate comparison for the 15 shock points are built on this classification, the headline quantitative result is not robust to this choice. A sensitivity analysis should be provided, and the floor should be justified with actual measurement-uncertainty information.","section":"§4.1"},{"comment":"The caption admits that the shock-model grid is coarser than the flow-model grid and asserts without demonstration that 'this does not impact our overall results.' If the shock locus is missing regions of physically allowed parameter space because of the coarse grid, data points in those regions could be misclassified as flow, both, or neither. The classification counts and the subsequent accretion-rate comparison for shock-dominated points are therefore potentially sensitive to the grid resolution. The authors should either run a finer shock grid or demonstrate quantitatively that the existing grid covers the relevant parameter space at sufficient resolution, e.g., by showing that no gap in the shock locus is large enough to affect the classification of any of the 89 data points.","section":"§3, Figure 1 caption"},{"comment":"The chromospheric-activity exclusion of 7 data points relies on the Manara et al. (2013, 2017) relation, which is validated only for 3.35 ≲ log(Teff/K) ≲ 3.65. The coolest object, 2M1115 (J11151597+1937266, log Teff ≈ 3.23), is classified as shock-dominated in Table 2, and this classification depends on an extrapolation of the chromospheric-activity locus down to log(Teff/K) = 3.2. If the extrapolation is incorrect, 2M1115 could be chromospherically dominated and should be excluded, reducing the shock sample from 15 to 14 points and changing the accretion-rate comparison in §5.2. The paper should test the sensitivity of the classification to this extrapolation, for example by excluding 2M1115 from the shock sample and re-evaluating the main results, or by using an alternative chromospheric-activity criterion that does not require extrapolation.","section":"§5.1, Appendix B"}],"minor_comments":[{"comment":"The abstract states that the shock model gives accretion rates 'up to several times higher' than the stellar scaling, but Figure 5 shows ratios reaching approximately 10 (one order of magnitude) for the lowest accretion rates; consider rephrasing to 'up to an order of magnitude higher' to match the figure.","section":"Abstract and §5.2"},{"comment":"The classification uses 1σ error ellipses, but the choice of 1σ is not discussed; a brief justification or a check of how the counts change with, e.g., 2σ ellipses would strengthen the interpretation.","section":"§4.1"},{"comment":"In the table caption, the column labeled 'Em?' is not defined; it would be clearer to spell out 'Emission mechanism' in the caption.","section":"Table 2"},{"comment":"The axis labels in panel (b) appear garbled (the text '8 654321 654321' near the axes); please check the figure and ensure the tick labels are rendered correctly.","section":"Figure 1"},{"comment":"The sentence 'Strictly speaking, the chromospheric activity relation in Manara et al. (2013, 2017) cannot be applied to objects with shock-dominated emission' is followed by an argument that the identification remains the same; this argument is not fully convincing and should be supported by a more explicit derivation or a test using both conversions.","section":"§5.1"}],"recommendation":"major_revision","confidential_remarks":"The paper applies existing models developed in part by one of the authors (Aoyama et al. 2018, 2021), which is appropriate for this kind of archival survey, and there is no indication of circular reasoning. The main concern is methodological robustness: the classification census is sensitive to the ad hoc error floor and the coarse shock grid, and the inclusion of 2M1115 relies on an extrapolated chromospheric relation. These issues are fixable with additional analysis, so a major revision is appropriate. The paper is within the scope of the journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is the first systematic application of the shock-vs-flow line-ratio diagnostic to a large sample of accreting brown dwarfs and very low-mass stars. It deserves a real referee, but the quantitative census is softer than the abstract implies.\n\nWhat is new: they compile 96 archival data points for 76 objects, use Hβ/Hγ versus Hβ/H8 to separate the two model loci (H8 is a good addition), and classify 89 points. The result that both shock- and flow-dominated emission occur in this mass regime is credible and interesting; the same conclusion was suggested for individual objects before, but not on a sample this size. They also add Subaru/HDS line profiles for three objects and document epoch-to-epoch category switches. That variability discussion is honest.\n\nThe main soft spot is the 10% error floor in §4.1. For any point without errors or with sub-10% errors, they inflate the 1σ ellipse to 10%. Many high-SNR X-shooter measurements have formal errors around 1%, so the floor is doing real work. Larger ellipses are more likely to overlap a model locus, so the floor systematically shrinks the 'neither' class and can shuffle points among shock/flow/both. Since the headline is 15 shock versus 55 flow, and the accretion-rate comparison is done only for those 15, this is load-bearing. They do not provide a sensitivity test (e.g., 5% floor, or using only reported errors). That should be straightforward to run and needs to be in the paper.\n\nAlso minor: the shock grid is coarser than the flow grid, and the one-sentence claim that this does not matter (Figure 1 caption) is an assertion, not a demonstration. A denser shock grid would make the test fairer. And the final leap—'detection rate is independent of emission mechanism'—goes beyond what the data show; it is a speculation that should be flagged as such.\n\nThe chromospheric-activity selection and the Lacc–Lline conversions are handled with reasonable care, and they explicitly acknowledge the circularity between Manara et al. and Alcalá et al. The single-object Teff extrapolation is minor.\n\nBottom line: cite this for the H8 diagnostic and the qualitative result that both mechanisms exist. Do not quote the 15/55 numbers until they have shown robustness to the error floor. Send it to review; a good referee will ask for the sensitivity test and a denser shock grid.","headline":"A useful first census of Balmer-line emission mechanisms in brown dwarfs, but the headline 15/55 split rests on an unvalidated 10% error floor.","tokens_in":32664,"tokens_out":2601,"would_cite":true,"duration_ms":27663,"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":"For 89 archival data points on accreting brown dwarfs and very low-mass stars, Balmer line ratios identify 15 as shock-dominated and 55 as flow-dominated emission.","keywords":["Accretion","Brown dwarfs","Exoplanet formation","Free floating planets","H I line emission","L dwarfs","M stars","Pre-main sequence stars"],"falsifier":"Compute a radiative-transfer model that includes both the infalling accretion flow and the shock at its base as a single structure; if that hybrid model fills the empty region between the two loci in the $\\mathrm{H}\\beta/\\mathrm{H}\\gamma$ versus $\\mathrm{H}\\beta/\\mathrm{H}8$ diagram, then the binary 15-of-89 versus 55-of-89 classification loses its foundation. Alternatively, observe the five objects that switched categories with simultaneous Hβ and Paschen lines to see whether the flip is reproduced in an independent hydrogen series.","tokens_in":31636,"feed_emoji":"🔭","tokens_out":14930,"duration_ms":135690,"temperature":0.7,"pith_summary":"The paper asks a practical question: when the hydrogen emission lines of a brown dwarf or very low-mass star are used to measure how fast it is swallowing gas from its disk, which emission model should be trusted? The authors apply the two competing models—emission from an accretion shock and emission from the accretion flow—to archival simultaneous Balmer line ratios (Hβ/Hγ versus Hβ/H8) for 96 data points from 76 objects with masses of roughly 0.02 to 0.1 solar masses. They find 15 points that fall only on the shock locus, 55 on the flow locus, and the remaining points are ambiguous or chromospheric. For the shock-dominated points, the shock model gives mass accretion rates up to several times higher than the usual stellar scaling, and the gap widens at lower accretion rates, where planetary-mass objects live. If the paper is right, the choice of emission model is not a detail: it changes the inferred growth rate and the physical picture of how substellar objects assemble.","feed_headline":"15 of 89 brown dwarfs accrete by shock, 55 by inflow","feed_subtitle":"The choice of emission model can raise inferred growth rates severalfold, more toward planetary masses.","key_machinery":"The load-bearing object is the $\\mathrm{H}\\beta/\\mathrm{H}\\gamma$ versus $\\mathrm{H}\\beta/\\mathrm{H}8$ line-ratio diagram. $\\mathrm{H}\\beta/\\mathrm{H}\\gamma$ mainly tracks the temperature of the emitting gas, while the inclusion of H8, the $n=8$ Balmer line at 388.9 nm, separates the two models: in the accretion flow model the higher levels are populated partly by recombination, so $\\mathrm{H}\\beta/\\mathrm{H}8$ is lower than in the shock model. The shock locus comes from the post-shock cooling gas modeled with one-dimensional thermo-hydrodynamics, chemistry, and radiative transfer, and the flow locus comes from a constant-temperature slab with Sobolev escape probabilities. The diagram is what lets the authors turn archival multi-line spectra into a per-object mechanism assignment; without the H8 ratio the two loci overlap too much in the lower-left corner.","core_discovery":"The central claim is that archival Balmer line ratios can classify the dominant accretion-emission mechanism of individual brown dwarfs and very low-mass stars: 15 of 89 usable data points are best described by the accretion shock model, while 55 are best described by the accretion flow model. The diagnostic is the ratio pair $\\mathrm{H}\\beta/\\mathrm{H}\\gamma$ versus $\\mathrm{H}\\beta/\\mathrm{H}8$, which separates the two model loci because the higher hydrogen level H8 is populated by recombination as well as by collisional excitation, lowering $\\mathrm{H}\\beta/\\mathrm{H}8$ in the flow model. For the 15 shock-classified points, converting the Hβ luminosity through the shock model yields mass accretion rates up to several times those obtained by extrapolating stellar scaling relations, and the discrepancy grows as the accretion rate falls; if the empirical $\\dot{M}\\propto M^2$ trend holds into the planetary regime, the difference can reach one to two orders of magnitude. The paper also reports that the ratio classification flips between shock and flow at different epochs for five objects, and that high-resolution line profiles show central absorption in a flow-dominated object and a single peak in a shock-dominated object.","pith_inferences":["Extension: applying the same two-locus test to directly imaged planets using Paschen or Brackett line ratios could test whether the shock fraction grows at planetary masses; the paper only extrapolates the mass-accretion-rate trend and does not claim the fraction itself would grow.","Extension: the five objects that flip categories between epochs hint that the accretion geometry or inflow temperature changes on observable timescales, so a dedicated high-cadence monitoring campaign of Balmer ratios could reveal what drives the switch.","Extension: a radiative-transfer model that includes both the infalling flow and the shock at its base, treated as one structure rather than two alternatives, would show whether intermediate configurations fill the gap between the loci; if they do, the 15/55 split would need reinterpretation."],"forward_implications":["A single multi-line spectrum covering Hβ, Hγ, and H8 can place a brown dwarf or very low-mass star on the shock or flow locus, giving a per-object mechanism assignment rather than a one-size-fits-all model.","For shock-dominated objects, accretion rates estimated from stellar scaling relations should be revised upward; in this sample the correction is up to several times, and the paper's extrapolation to planetary masses widens it to one to two orders of magnitude.","The flow model still describes the majority of the sample, so stellar-origin accretion flow treatment remains relevant for most substellar objects.","Five objects switch between the shock and flow categories at different epochs, so a single-epoch classification is not a permanent label for the object.","Line profiles are a secondary diagnostic: central absorption points to flow emission, while a single central peak can be either shock or high-temperature flow, making ratios the more reliable discriminator."],"supporting_citations":[{"why":"Supplies the accretion shock model whose post-shock cooling calculations produce the shock locus in the Balmer line-ratio diagram.","marker":"Aoyama et al. (2018)"},{"why":"Supplies the accretion flow slab model with Sobolev escape probabilities that defines the flow locus.","marker":"Kwan & Fischer (2011)"},{"why":"Provides the shock-model conversion from H-beta luminosity to accretion luminosity used to compute shock-based mass accretion rates.","marker":"Aoyama et al. (2021)"},{"why":"Provides the stellar scaling conversion from H-beta luminosity to accretion luminosity used for the comparison rates and for the chromospheric-activity diagnostic.","marker":"Alcalá et al. (2017)"},{"why":"Showed on a planetary-mass object that multiple hydrogen-line ratios can distinguish shock from flow emission, the method this paper generalizes.","marker":"Betti et al. (2022)"},{"why":"T Tauri accretion flow simulations whose lower H-beta/H-gamma ratios provide the comparison sample for the statistical test of shock versus flow fractions.","marker":"Muzerolle et al. (1998a)"},{"why":"One of the archival sources of the measured Balmer line fluxes and the Lacc-line calibrations feeding the sample.","marker":"Herczeg & Hillenbrand (2008)"}],"fun_headline_variants":["Shock model lifts brown dwarf accretion rates severalfold","15 brown dwarfs accrete via shock, flow dominates in 55","Accretion shock vs flow: ratio test classifies 70 dwarfs","Balmer lines reveal shock accretion in accreting brown dwarfs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole classification rests on the assumption that shock and flow emission are the only two ways the gas can make the observed Balmer lines, so any point that fits neither locus is written off as chromospheric activity or noise, and that the shock model's coarser grid does not hide additional valid configurations.","fun_headline_variants_meta":{"raw":{"variants":["Shock model lifts brown dwarf accretion rates severalfold","15 brown dwarfs accrete via shock, flow dominates in 55","Accretion shock vs flow: ratio test classifies 70 dwarfs","Balmer lines reveal shock accretion in accreting brown dwarfs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000269,"raw_usage":{"total_tokens":1674,"prompt_tokens":1047,"completion_tokens":627,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":663,"completion_tokens_details":{"reasoning_tokens":563}},"tokens_in":663,"tokens_out":627,"duration_ms":6846,"temperature":1.0,"reasoning_tokens":563,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T17:52:05.010304+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute a radiative-transfer model that includes both the infalling accretion flow and the shock at its base as a single structure; if that hybrid model fills the empty region between the two loci in the $\\mathrm{H}\\beta/\\mathrm{H}\\gamma$ versus $\\mathrm{H}\\beta/\\mathrm{H}8$ diagram, then the binary 15-of-89 versus 55-of-89 classification loses its foundation. Alternatively, observe the five objects that switched categories with simultaneous Hβ and Paschen lines to see whether the flip is reproduced in an independent hydrogen series.","supporting_citations":[],"review_version":1}