{"id":"7416fd3a-5e0f-467f-b3a0-8a3df179a693","arxiv_id":"2607.16691","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In a viscously expanding disk, the relative locations of the reaction line and the stagnation line determine whether chemically processed dust is retained or exported; this can explain cometary crystalline silicates with a compact, moderately massive, weakly turbulent protosolar disk.","lead":"This paper uses simulations to show that in a protoplanetary disk that spreads outward over time, the position of the chemical 'reaction line' relative to the flow-reversing 'stagnation line' controls whether processed dust is exported to the outer disk. A smart generalist might read it because it offers a natural explanation for the puzzling presence of crystalline silicates in comets and suggests the early solar nebula was compact, moderately massive, and weakly turbulent.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The stagnation-line mechanism may be an artifact of splicing a ν∝r self-similar velocity field onto an inner region where ν∝r^{3/5}; the compact-disk inference is not secure until the radial velocity is recomputed self-consistently.","rationale":"The paper's central claim and its Solar System inference both rest on the reaction-line/stagnation-line dichotomy. That dichotomy is implemented through Eq. (9) and Eq. (37), which are properties of the self-similar outer solution, not of the combined model used for temperatures. The mismatch is not a matter of missing microphysics; it is an internal consistency issue in the transport model itself. The reader's weakest assumption (dust as a perfect gas tracer) is also a real limitation, but it concerns a standard simplification and could be tested separately by adding size-dependent drift/pebble transport. The velocity inconsistency is more fundamental because even for perfectly coupled dust the stagnation line used for classification may not be the real one. If the concern survives the proposed check, the paper would need major revision; if it does not, the conditional acceptance stands. I therefore keep the reader's CONDITIONAL verdict (UNCHANGED) while emphasizing the specific check.","tokens_in":15890,"tokens_out":14248,"duration_ms":147568,"concrete_test":"Take one representative compact/massive case (e.g., α=10^{-3}, rd0=2 au, Md0=0.1 or 0.2 M_sun) and replace the prescribed Eq. (9) by the radial velocity obtained from mass conservation with the combined Σ, e.g., solve the 1-D viscous evolution equation using the adopted ν(r) and the inner boundary condition, or compute v_r=-(3/(Σ√r)) ∂(Σν√r)/∂r from the actual combined Σ field. Re-run the Monte Carlo and re-measure the crystalline fraction in 3-100 au, 1-3 Myr. If the zero-advection radius is substantially interior to rd(t)/2 or outward transport drops, the compact-disk conclusion fails; if the crystallinity maps are essentially unchanged, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing concern is that the Monte Carlo transport uses a gas velocity field that is not the one implied by the paper's own disk model. Section 2 splices an inner viscous-heating solution (Eqs. 20-23) onto the outer self-similar solution. In the inner region T_vis,c ∝ r^{-9/10} and H_vis ∝ r^{21/20}, so ν=αHc_s ∝ r^{3/5}, not ∝ r. Yet the radial advection velocity (Eq. 9) and the stagnation line r_stag=rd(t)/2 (Eq. 37) are derived for the ν∝r self-similar solution, and Eqs. (29)-(30) apply this velocity at all radii. For the compact, massive disks that produce the Solar System match (e.g., α=10^{-3}, rd0=2 au, Md0=0.2 M_sun), the viscous-heating region initially extends beyond rd/2≈1 au; at those radii the inner-region solution would give v_r≈-3ν/(2r)<0, whereas Eq. (9) gives outward flow. Thus the predicted outward export of crystalline silicates beyond the 'stagnation line' may be an artifact of inconsistent splicing, and the inferred constraints rd0<5 au, Md0>0.05 M_sun, α≲10^{-3} are not secure until the velocity is recomputed self-consistently.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a time-dependent, viscously evolving protoplanetary disk model that splices an inner viscous-heating solution onto an outer self-similar irradiation-dominated solution, and couples it to Monte Carlo particle-tracking simulations with irreversible reaction kinetics. The central idea is that the spatiotemporal distribution of reacted dust is controlled by the relative positions of a \"reaction line\" (where a reaction completes efficiently) and a \"stagnation line\" (where radial advection changes direction). The model is applied to crystalline silicates in comets, and the authors conclude that the protosolar disk was likely compact (rd0 < 5 au), moderately massive (Md0 > 0.05 M_sun), and not strongly turbulent.","tokens_in":16323,"tokens_out":3379,"duration_ms":37818,"significance":"If the central claim holds, the paper provides a useful time-dependent framework for connecting irreversible dust chemistry to disk evolution, extending the reaction-line idea from steady disks to expanding disks and making testable predictions for spatially resolved observations. The use of experimentally calibrated JMA kinetics and a large parameter grid is a strength. However, the main physical conclusion rests on a velocity field whose applicability to the inner, reaction-relevant region is not established, and on an acknowledged disk-model inconsistency. The quantitative constraints on the protosolar disk are therefore not secure in the present form.","major_comments":[{"comment":"The central mechanism for outward export of reacted dust is the stagnation line r_stag = rd(t)/2, derived from the radial velocity vr = -3ν/(2r)[1 - r/(rd/2)]. This velocity is the self-similar solution for the irradiation-dominated outer disk, where ν∝r^{15/14}≈r. But in the inner viscous-heating region the paper's own model gives T_vis,c ∝ r^{-9/10} and H_vis ∝ r^{21/20}, hence ν=αHc_s ∝ r^{3/5}, not ∝r. Equations (29)–(30) nevertheless use the outer-disk vr at all radii. For the compact, massive disks that produce the claimed match (e.g., α=10^{-3}, rd0=2 au, Md0=0.2 M_sun), the viscous region initially extends well beyond rd/2, so the inner-region velocity would be inward rather than outward. The inferred constraint rd0<5 au, and the entire two-line classification, may be an artifact of this inconsistent splicing. The radial velocity must be recomputed self-consistently for the combi","section":"§2, Eqs. (9), (29)–(30), (37); §3.4"},{"comment":"The manuscript explicitly acknowledges that the chosen splicing T=max(T_vis,T_irr), Σ=min(Σ_vis,Σ_irr) produces a \"mass defect\" that causes an artificial increase in total disk mass and hence unphysical mass evolution. The response that the Monte Carlo simulation is the appropriate reference because it tracks the actual mass accretion rate does not resolve the inconsistency: the simulation uses advection velocities derived from the self-similar solution, while the surface density is modified by Eq. (25). The two are not mutually consistent. This affects both the particle trajectories and the quantitative comparison with cometary crystallinity in Fig. 4.","section":"§2, Eqs. (24)–(25), and the following paragraph"},{"comment":"The extension of the predictive reaction-line formula to the combined disk introduces an ad hoc floor on the accretion rate, Mdot=10^{-9} M_sun/yr, whenever Eq. (36) gives a value near zero or negative. This is explicitly outside the original steady-accretion formulation. Although the authors argue the reaction-line temperature is only weakly sensitive to Mdot, the floor is applied precisely in the compact-massive-disk cases that drive the main Solar System conclusion. A sensitivity test varying this floor (e.g., 10^{-10} or 10^{-8}) is needed to show that the inferred disk-parameter range does not depend on this choice.","section":"§3.1, Eq. (36) and the paragraph following it"},{"comment":"The load-bearing assumption that dust grains are perfect tracers of gas, with D=ν and no relative drift, is stated but not relaxed. If cometary crystalline silicates were transported as grains embedded in larger pebbles, or underwent radial drift, the effective outward transport and the inferred disk parameters (rd0<5 au, Md0>0.05 M_sun, α≲10^{-3}) would not directly apply. This is a scope limitation, but it should be stated more prominently in the conclusions rather than only as a modeling simplification.","section":"§2.1, Eqs. (26)–(31)"}],"minor_comments":[{"comment":"Notation is inconsistent between rd0, r0, and rd(t); e.g., §3.4 uses r0<5 au while the abstract and §2 use rd0. Please unify.","section":"Throughout"},{"comment":"Typographical issues: \"T able 1\" in the caption and \"initally\" in the text. Also \"r0 <5,au\" contains a misplaced comma.","section":"§2, Table 1"},{"comment":"The statement \"not strongly turbulent\" is vague; the grid only tests α=10^{-2} and 10^{-3}. Please state explicitly that the conclusion distinguishes these two values, not a continuous range.","section":"§3.4"},{"comment":"The adopted definition of the comet-forming region (r=3–100 au, t=1–3 Myr) is arbitrary. Some discussion or sensitivity test regarding this choice would strengthen the application.","section":"§3.4, Fig. 4 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is potentially interesting, but the central claim is not yet safe. The velocity-field inconsistency identified in the stress-test note is real and load-bearing: the stagnation-line mechanism is derived for a ν∝r self-similar disk, while the inner region in which reactions occur has ν∝r^{3/5}. The authors should be asked to recompute the radial velocity self-consistently and rerun the analysis, or at minimum to demonstrate that the discrepancy does not alter the two-line classification. The acknowledged mass-defect issue and the ad hoc accretion-rate floor are additional concerns that need explicit sensitivity checks."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe genuinely new stuff: the authors extend the reaction-line formalism from their steady-disk paper to a time-dependent viscously spreading disk, and show with Monte Carlo runs that the final radial distribution of processed dust is shaped by whether the reaction line sits inside or beyond the stagnation line (where advection reverses). That classification is a nice organizing device, and the paper is transparent about its own limits — the mass-defect issue in Section 2, the fixed opacity, the prototype reaction kinetics. Credit where due: this is a step beyond steady-disk treatments and it makes concrete, falsifiable predictions about the spatiotemporal distribution of processed dust.\n\nThe soft spots, in order of severity:\n\n1. The stagnation-line mechanism may be partly an artifact of splicing. The radial velocity, Eq. (9), and r_stag = rd(t)/2 are derived from the ν∝r self-similar solution in the outer irradiated disk. But in the inner viscous-heating region, their own T and H scalings give ν∝r^{3/5}. They apply the same Eq. (9) everywhere, including the compact, massive disk models that later produce the Solar System match. For those models the viscous region initially extends past rd/2, so at those radii the inner-region velocity would be inward, not outward. If that is right, the outward export of crystalline silicates could be an artifact of an inconsistent velocity field, not a physical outcome of disk evolution. That is the main thing I would want fixed before trusting the inferred constraints rd0 < 5 au, Md0 > 0.05 Msun, α ≲ 1e-3.\n\n2. The tracer assumption: dust is treated as perfectly coupled gas (D = ν, no drift, no pebbles). The paper states this plainly, but the application to comets depends on it. If the crystalline grains were transported as pebbles, as Okamoto & Ida found, the disk parameter constraints do not directly transfer. The comparison to the observed 10–60% crystalline fractions is also qualitative — no error bars and no specific comet formation scenario. That is a fair caution, not a fatal one.\n\n3. No code or data are released. With this many Monte Carlo runs, the key figures are hard to reproduce from the text alone.\n\nThe predictive reaction-line formula is borrowed from their own prior work, but that earlier formula is independently calibrated and the present tests are against new simulations, so I do not see a circularity problem.\n\nOverall: the two-line framework is a useful contribution and the paper deserves a serious referee, but the Solar System conclusion is not secure until the velocity field is recomputed self-consistently across the inner/outer splice. I would not build on the inferred disk parameters as-is. I would send it to review with a clear request for that calculation.","headline":"Plausible new mechanism for cometary crystalline silicates, but the disk inferences rest on a velocity-field splice that may produce the effect, and on treating dust as a perfect gas tracer.","tokens_in":16723,"tokens_out":6692,"would_cite":false,"duration_ms":67836,"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 claims that in a viscously evolving protoplanetary disk, the radial distribution of chemically processed dust is controlled by the relative positions of the reaction line and the stagnation line, and that the crystalline silicate","keywords":["protoplanetary disks","dust transport","crystalline silicates","comets","irreversible reactions","reaction line","stagnation line","viscous disk evolution"],"falsifier":"Run the same Monte Carlo model with dust radial drift included (Stokes numbers of order 0.01–1) and ask whether crystalline fractions of 10–60% still reach 3–100 au in a compact, weakly turbulent disk; if they cannot, the gas-tracer premise and the derived protosolar disk parameters collapse. A robust observational determination that the protosolar disk's initial characteristic radius exceeded 5 au would also falsify the compact-disk inference.","tokens_in":15813,"feed_emoji":"☄️","tokens_out":7686,"duration_ms":74954,"temperature":0.7,"pith_summary":"This paper asks where dust that has undergone an irreversible chemical reaction—such as crystallization of amorphous silicates—ends up in a protoplanetary disk that is viscously spreading and cooling. The authors claim that the final radial distribution of this reacted dust is set by the relative positions of two radii: the reaction line, inside which the reaction completes efficiently, and the stagnation line, outside which gas and well-coupled dust move outward. When the reaction line lies near or beyond the stagnation line, reacted dust is efficiently exported to the cold outer disk; when the two lines are widely separated, processed dust stays trapped in the inner disk. Applied to the 10–60% crystalline silicates observed in Solar System comets, the model implies the protosolar disk was compact (initial characteristic radius below 5 au), moderately massive (above 0.05 solar masses), and not strongly turbulent. The significance is that cometary crystallinity can arise naturally from disk evolution alone, without invoking extra transport mechanisms such as sticking of grains onto icy pebbles.","feed_headline":"Comet crystals imply a compact, calm protosolar disk","feed_subtitle":"Model shows crystals reach comet-forming regions only when the reaction line meets the stagnation line.","key_machinery":"The controlling object is the pair of radii the paper calls the reaction line and the stagnation line. The reaction line is the temperature-defined boundary where an irreversible reaction completes efficiently; a predictive formula originally developed for steady accretion disks is extended to viscously expanding disks by evaluating the local accretion rate at the characteristic reaction temperature derived from reaction kinetics. The stagnation line is the radius r_d(t)/2 at which the radial advection velocity changes sign: particles inside drift inward and accrete, particles outside drift outward. The Monte Carlo particle-tracking model moves dust with the gas (advection plus diffusion, wi","core_discovery":"The central claim is that the spatiotemporal distribution of dust that completed an irreversible reaction in a viscously evolving disk is governed by the relative locations of the reaction line and the stagnation line. The reaction line—the temperature where a reaction completes efficiently—moves inward as the disk cools, while the stagnation line, where radial advection reverses direction, moves outward as the disk viscously spreads. If the reaction line lies far inside the stagnation line, reacted dust remains confined to the inner disk; if the reaction line lies near or beyond the stagnation line, reacted dust located near or between the two lines is carried outward efficiently, producing","pith_inferences":["If the compact-initial-disk picture is right, young protoplanetary disks with small characteristic radii should show crystalline silicate emission at large radii within the first few million years—a direct test with spatially resolved mid-infrared observations.","The main caveat is the gas-tracer assumption: including radial drift of larger grains or transport inside pebbles could change the inferred disk parameters, possibly widening the allowed range of disk size and turbulence.","The same reaction-line/stagnation-line criterion could be applied to other irreversible tracers such as isotopic anomalies or organic-matter destruction, potentially turning primitive-meteorite and comet records into a disk-evolution chronometer.","The predicted fossil population effect implies that inner-disk samples found today are survivors of accretion, not records of peak processing; interpretations of thermal histories in meteoritic materials should account for this selection effect."],"forward_implications":["Outward transport of crystalline silicates emerges naturally from the viscous spreading of an initially compact disk, so cometary crystallinity need not require additional transport mechanisms such as grains sticking onto icy pebbles.","The reaction-line formula gives a fast analytic way to predict where any irreversible reaction's products will end up in an evolving disk, from kinetic parameters and disk conditions alone.","In strongly turbulent disks the model predicts a fossil population of reacted dust surviving in the outer disk after the hot inner region cools and accretes onto the star, so processed and unprocessed material become spatially decoupled.","The observed 10–60% crystalline fraction and comet-to-comet olivine-to-pyroxene diversity translate into concrete protosolar disk constraints: initial characteristic radius below 5 au, initial mass above 0.05 solar masses, and weak turbulence.","Gas-phase molecules released by thermal decomposition of refractory organics should show radial distributions correlated with processed dust, offering a testable link to spatially resolved observations."],"fun_headline_variants":["Reaction and stagnation line alignment seeds comet crystals","Crossing disk lines spread reacted dust to outer regions","Comet crystals imply compact, calm protosolar disk","Disk lines' meeting point sends crystals outward to comets","Crystalline comets trace a compact, non-turbulent disk"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"Dust grains are treated as perfect gas tracers—sharing the gas velocity, with diffusivity set equal to viscosity and no radial drift relative to the gas (Section 2.1); if the comet crystals were transported inside larger drifting pebbles, the inferred compact, moderately massive, weakly turbulent protosolar disk would not follow.","fun_headline_variants_meta":{"raw":{"variants":["Reaction and stagnation line alignment seeds comet crystals","Crossing disk lines spread reacted dust to outer regions","Comet crystals imply compact, calm protosolar disk","Disk lines' meeting point sends crystals outward to comets","Crystalline comets trace a compact, non-turbulent disk"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000237,"raw_usage":{"total_tokens":1370,"prompt_tokens":794,"completion_tokens":576,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":538,"completion_tokens_details":{"reasoning_tokens":496}},"tokens_in":538,"tokens_out":576,"duration_ms":6732,"temperature":1.0,"reasoning_tokens":496,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T20:11:50.883958+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same Monte Carlo model with dust radial drift included (Stokes numbers of order 0.01–1) and ask whether crystalline fractions of 10–60% still reach 3–100 au in a compact, weakly turbulent disk; if they cannot, the gas-tracer premise and the derived protosolar disk parameters collapse. A robust observational determination that the protosolar disk's initial characteristic radius exceeded 5 au would also falsify the compact-disk inference.","supporting_citations":[],"review_version":1}