{"id":"74263c5d-1e8a-442d-8e1c-c1b9b01b6705","arxiv_id":"2501.12026","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Radiative torque disruption destroys weak aggregate grains up to 500 micrometers within less than 2 years in luminous protostellar outflows, while strong composite grains survive and migration continues.","lead":"This paper simulates whether fast-spinning dust grains, spun up by starlight, shatter as they travel through a young star's jet and outflow. It shows weak fluffy grains break apart quickly, which changes what ALMA polarization observations can tell us about grain sizes and migration.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The '<2 yr' aggregate-grain disruption claim rests on compact-grain RAT efficiencies, while the paper's own Footnote 2 cites aggregate Q_Gamma values 10–100x smaller, which can shift disruption timescales by one to two orders of magnitude.","rationale":"The paper is a competent and transparent extension of prior RATD work: it couples grain transport to an MHD jet/outflow simulation, incorporates RATD into POLARIS, explicitly compares the dynamic and rest-grain approaches, and admits several limitations, including the aggregate-Q_Gamma issue in Footnote 2, the artificial L_center = 100 L_sun choice, and the single MCRT iteration. The reader's weakest assumption identifies exactly the load-bearing issue I see: the headline aggregate-grain destruction timescale is computed with compact-grain RAT efficiencies. This is not an internal mathematical error in the compact-grain calculation; it is a mismatch between the grain type featured in the central claim and the radiative-torque model applied to it. A 10–100x reduction in Q_Gamma changes disruption timescales by one to two orders of magnitude and, because grains keep moving, can remove them from the high-radiation region before disruption occurs. This lands directly on the abstract, not merely on the quantitative polish. I agree with the reader's characterization, and I do not see a reason to move the verdict: CONDITIONAL remains appropriate because the qualitative weak-vs-strong tensile-strength separation is plausibly robust for the lowest S_max values, and the paper already shows that composite grains survive when transport is included. However, the S_max <= 1e5 erg/cm^3 boundary and the <2 yr statement should not be treated as quantitatively settled until the calculation is redone with aggregate-specific Q_Gamma values or the abstract is explicitly caveated as a compact-grain upper-limit scenario.","tokens_in":46897,"tokens_out":5619,"duration_ms":63525,"concrete_test":"Recompute the upper-row panels of Figure 2 with aggregate-grain RAT efficiencies from Jager et al. (2024), or from an independent DDSCAT/ADDA calculation for porous/ballistic aggregates, inserted into Eqs. A1–A4 using the same radiation field and gas density at the jet/outflow base. Evaluate t_disr for a_eff = 1, 10, 100, and 500 um with S_max = 1e5 erg/cm^3. If t_disr exceeds ~2 yr, or exceeds the advection time to leave the region where u_rad/u_ISRF is high, then the abstract's '<2 yr' claim and the disruption-size maps in Figures 3–4 are overestimates and should be restated with an explicit aggregate-Q_Gamma caveat. Repeat the same test at L_center = 20 L_sun to determine whether the quoted >= 20 L_sun activation threshold survives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract's central quantitative claim is that for bolometric luminosity >= 20 L_sun, RATD destroys aggregate grains of 1–500 um with S_max <= 1e5 erg/cm^3 in the jet/outflow base in <2 yr. This depends directly on the spin-up rate in Eqs. A1–A4, which adopts compact oblate-grain RAT efficiencies from Lazarian & Hoang (2007a), Hoang & Lazarian (2008), and Herranen et al. (2021). The paper's own Footnote 2 cites Jager et al. (2024), which finds aggregate-grain Q_Gamma is 10–100 times smaller at lambda/a_eff ~1–100, and then states that spin-up is 'slower by twice times'—an internal inconsistency with the cited factor. Since Eq. A4 gives Omega(t) ~ Gamma_RAT t / I at early times, the time to reach Omega_disr scales approximately as 1/Q_Gamma. A 10–100x reduction in Q_Gamma would turn a claimed <2 yr disruption time into roughly 20–200 yr, or no disruption at all if the saturated Omega falls below Omega_disr. During that longer time, grains moving at 6–300 km/s traverse hundreds to thousands of au and leave the strong-radiation base, so the RATD-active region and the migration-suppression picture in Figures 3–4 would shrink substantially. The qualitative weak-vs-strong tensile-strength hierarchy may survive for very weak grains with S_max ~1e3–1e4 erg/cm^3, but the specific <2 yr timescale and the S_max <= 1e5 erg/cm^3 boundary in the abstract are not supported by the RAT model actually adopted for aggregate grains.","agreement_with_reader":"agree"},"referee_report":null,"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, here's my read of Giang et al. (arXiv:2501.12026).\n\nThe genuinely new piece is the dynamic coupling: they solve grain spin-up along outflow streamlines using the MHD gas structure, instead of the static analytical RATD used in Hoang et al. (2021) and Le Gouellec et al. (2023b). They also implement RATD in POLARIS and compare the static POLARIS picture against the dynamic one. That comparison is worthwhile and honestly done: they show POLARIS works for weak aggregate grains but overestimates disruption for composite grains, and they quantify the discrepancy as roughly a factor of two. They also test several assumptions the previous papers ignored: initial rotation, one-round vs iterated MCRT, and luminosity cases down to 5 L_sun. That is real progress, and the qualitative conclusion — weak fluffy aggregates get destroyed near the outflow base, strong composite grains survive — is plausible and internally consistent.\n\nThe soft spot is the one the stress-test flagged, and it is load-bearing for the abstract. The calculations use compact oblate-grain RAT efficiencies (Eq. A1, Lazarian & Hoang 2007a, Hoang & Lazarian 2008, Herranen et al. 2021) for aggregate grains. Footnote 2 cites Jäger et al. (2024) saying aggregate Q_Gamma is 10–100 times smaller at lambda/a_eff ~1–100, but then says spin-up is only \"slower by twice times.\" That does not parse. Since early-time spin-up goes as Gamma_RAT t / I, a 10–100x lower Q_Gamma pushes the <2 yr disruption timescale to roughly 20–200 yr, during which grains at 6–300 km/s traverse hundreds to thousands of au and leave the strong-radiation region. So the specific claims — \"after <2 yrs\" and the S_max <= 1e5 erg/cm^3 boundary in the abstract — are not supported by the RAT model actually used for aggregates. The qualitative hierarchy likely survives for very weak grains, but the headline number should be rederived with aggregate Q_Gamma or restated with this caveat.\n\nMinor additional notes: Lcenter = 100 L_sun is an artificial choice to maximize RATD, though they do run 20 and 5 L_sun cases in the appendix; no code or data are released; and the MCRT iteration is acknowledged as one-round only. None of these sink the paper on its own. The paper's own discussion sections are honest about several limitations, which I read as a sign of care rather than a deflection.\n\nBottom line: this deserves a serious referee and likely publication after major revision. The referee should push the authors to fix the aggregate RAT efficiency inconsistency, rerun or clearly qualify the headline timescale, and ideally release the POLARIS RATD module. I would cite this work if the quantitative claim gets corrected; as it stands I'd cite it for the dynamic-vs-static comparison and the POLARIS implementation.","headline":"A solid, useful step forward in modeling RATD in protostellar outflows, but the headline <2 yr aggregate disruption time rests on compact-grain RAT efficiencies that the paper's own footnote undercuts, so the quantitative claim needs a fix or a caveat before publication.","tokens_in":47849,"tokens_out":1357,"would_cite":true,"duration_ms":17292,"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":"Radiative torque disruption can shatter weak aggregate grains up to 500 µm inside protostellar jets and outflows within two years, blocking the migration of very large grains from the disk to the inner envelope.","keywords":["radiative torque disruption","dust grains","protostellar jets","protostellar outflows","polarized dust emission","grain growth","young stellar objects","very large grains"],"falsifier":"Measure or compute the radiative-torque efficiency $Q_\\Gamma$ for realistic porous 1–500 µm aggregate grains at wavelengths from 0.1 µm to 3 mm; if the values are indeed 10–100 times below the compact-grain efficiencies used here, then the disruption timescales, disruption sizes, and the conclusion that RATD blocks large-grain migration must be revised. Observationally, detecting grains larger than about 100 µm inside the outflow or inner envelope of a Class 0/I protostar with bolometric luminosity near or above 100 $L_\\odot$ during an accretion burst would contradict the predicted dominance of submicron grains.","tokens_in":2314,"feed_emoji":"🌟","tokens_out":3220,"duration_ms":107533,"temperature":0.7,"pith_summary":"This paper asks whether radiative torque disruption (RATD) — the spinning-up and shattering of dust grains by starlight — can stop very large grains from travelling from the disk of a young star out through its jet and outflow and then settling into the inner envelope. Using the gas density and velocity structure of an MHD simulation of an intermediate-mass Class 0 protostar, the authors follow grains as they accelerate and spin up in the outflow. They find that if the protostar's bolometric luminosity is at least 20 $L_\\odot$, RATD destroys loose aggregate grains from 1 to 500 µm with tensile strength up to $10^5$ erg cm$^{-3}$ within less than two years in the jet/outflow base. That leaves submicron grains dominating the outflow and partially blocks the migration of large grains from the disk to the inner envelope, while compact composite grains with tensile strength above $10^7$ erg cm$^{-3}$ survive and continue migrating.","feed_headline":"Spinning starlight shatters large dust grains in young-star outflows","feed_subtitle":"Weak fluffy grains break apart, stopping large dust from migrating from disk to envelope.","key_machinery":"The central object is radiative torque disruption (RATD): an anisotropic radiation field exerts a net torque $\\Gamma_{\\rm RAT}$ on an irregular grain, spinning it up until the centrifugal stress exceeds the grain's maximum tensile strength $S_{\\max}$, at which point the grain breaks. The argument is carried by the comparison of two angular velocities: the spin $\\Omega(t,a)$ obtained by solving the rotational equation of motion with gas-drag damping, and the disruption threshold $\\Omega_{\\rm disr} = (2/a)(S_{\\max}/\\rho_{\\rm grain})^{1/2}$. The dynamic approach advects grains with the outflowing gas, updating the local radiation field and gas density along each trajectory and thereby producing a disruption size range $[a_{\\rm disr,dynamic}, a_{\\rm disr,max,dynamic}]$ as a function of position and time. A second, static-grain approach embeds RATD in the POLARIS radiative-transfer code, treating disruption as a local comparison of $\\Omega_{\\rm RAT}$ with $\\Omega_{\\rm disr}$ and then modifying the grain size distribution and polarization cross-sections.","core_discovery":"Inside the jet and outflow of an intermediate-mass Class 0 protostar, the authors claim, radiative torque disruption is the decisive grain-destruction channel. They solve for the angular velocity $\\Omega(t,a)$ gained by grains as they are accelerated by the outflowing gas, using the radiation field from a POLARIS post-processed MHD simulation. Comparing $\\Omega$ with the disruption threshold $\\Omega_{\\rm disr} = (2/a)(S_{\\max}/\\rho_{\\rm grain})^{1/2}$, they find that aggregate grains with $S_{\\max} \\leq 10^5$ erg cm$^{-3}$ and sizes $1\\!-\\!500\\,\\mu$m are shattered in less than two years in the jet/outflow base when the protostar's bolometric luminosity is at least 20 $L_\\odot$. Submicron fragments then dominate the outflow, partially preventing very large grains from migrating from the inner disk to the inner envelope. Composite grains with $S_{\\max} \\geq 10^7$ erg cm$^{-3}$ resist disruption and continue migrating. When RATD is inserted into POLARIS with grains held at rest, it reproduces the dynamic disruption pattern for slow aggregate grains but overestimates disruption for fast composite grains by about a factor of two; in the synthetic maps the polarization degree falls by a factor of two when weak aggregate grains are removed from the outflow cavity wall and inner envelope, although iron inclusions matter more than RATD for polarization.","pith_inferences":["If the cited result that aggregate radiative-torque efficiencies are 10–100 times smaller than compact-grain values at $\\lambda/a \\sim 1\\!-\\!100$ is correct, then the claimed sub-two-year destruction and the disruption-size maps for aggregates are likely optimistic; the migration-blocking conclusion for fluffy grains would weaken unless the luminosity or burst duration is larger than assumed.","A decisive observational test would compare outflow dust populations in the same high-luminosity Class 0/I protostar during an accretion burst and in quiescence: RATD predicts a temporary switch from very large grains to a submicron-dominated population within a few years of the burst turning on.","Because RATD converts large grains into submicron fragments, it shifts dust extinction toward UV–optical wavelengths and enriches the small-grain population; this should strengthen shock-produced SiO and other molecular tracers in jets, an effect the paper discusses qualitatively but has not yet folded into quantitative molecular-line predictions."],"forward_implications":["If RATD operates during an accretion burst with $L \\geq 20$ $L_\\odot$, aggregate grains up to 500 µm with $S_{\\max} \\leq 10^5$ erg cm$^{-3}$ are destroyed in the jet/outflow base in under two years, leaving submicron grains as the dominant outflow dust population for the burst lifetime of a few to a few hundred years.","The migration of weak, fluffy very large grains from the inner disk to the inner envelope is partially suppressed, while composite grains with $S_{\\max} \\geq 10^7$ erg cm$^{-3}$ survive and keep migrating.","Including RATD in POLARIS reduces the predicted polarization degree roughly twofold along the outflow cavity wall and inner envelope for aggregate grains with $S_{\\max} \\leq 10^4$ erg cm$^{-3}$, but leaves polarization unchanged for stronger grains.","The static-grain POLARIS implementation matches the dynamic disruption picture for slow-moving aggregate grains with gas velocity below 60 km/s, but overestimates disruption for fast-moving composite grains by roughly a factor of two.","Iron inclusions inside grains control the observed polarization more than RATD does, so ALMA polarization observations do not require composite or compact grain structures to explain the alignment efficiency in Class 0/I protostars."],"supporting_citations":[{"why":"Introduces radiative torque disruption and the disruption criterion $\\Omega_{\\rm disr}$ used throughout the paper.","marker":"Hoang et al. (2019)"},{"why":"Supplies the grain transport and terminal-velocity model, and the disk-to-envelope migration scenario this paper tests.","marker":"Wong et al. (2016)"},{"why":"Establishes the Ashfall dust cycle by which very large grains can migrate from the inner disk to the inner envelope.","marker":"Tsukamoto et al. (2021)"},{"why":"Prior RATD study in the same MHD outflow geometry that this paper extends with grain motion and POLARIS modeling.","marker":"Le Gouellec et al. (2023b)"},{"why":"Provides the radiative-torque formalism used to compute grain spin-up.","marker":"Lazarian & Hoang (2007a)"},{"why":"Provides the theory of grain alignment and radiative torques that underlies the RATD calculation.","marker":"Hoang & Lazarian (2008)"},{"why":"Supplies the compact-grain radiative-torque efficiencies adopted for aggregate grains in Appendix A.","marker":"Herranen et al. (2021)"},{"why":"Reports that aggregate radiative-torque efficiencies are 10–100 times smaller than compact-grain values at the relevant wavelengths, the main caveat on spin-up.","marker":"Jäger et al. (2024)"},{"why":"Provides the POLARIS radiative-transfer code into which RATD is incorporated.","marker":"Reissl et al. (2016)"},{"why":"Supplies the standard MRN grain size distribution used as the initial dust population.","marker":"Mathis et al. (1977)"}],"fun_headline_variants":["Radiative torques shatter large dust in young-star outflows","Weak fluffy grains break apart in protostellar jets","Starlight spins and shatters dust grains in outflows","Radiative torque disruption blocks large-grain migration","Outflow grain disruption halves polarized dust emission"],"cache_read_input_tokens":49792,"weakest_assumption_plain":"The calculation assumes that fluffy aggregate grains spin up under radiative torques as fast as compact grains, although a cited study reports that aggregate torque efficiencies can be 10–100 times lower at the relevant wavelengths; if aggregate grains spin up that much more slowly, the claimed sub-two-year destruction and the blocking of large-grain migration would be substantially weakened.","fun_headline_variants_meta":{"raw":{"variants":["Radiative torques shatter large dust in young-star outflows","Weak fluffy grains break apart in protostellar jets","Starlight spins and shatters dust grains in outflows","Radiative torque disruption blocks large-grain migration","Outflow grain disruption halves polarized dust emission"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000563,"raw_usage":{"total_tokens":2797,"prompt_tokens":1194,"completion_tokens":1603,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":810,"completion_tokens_details":{"reasoning_tokens":1523}},"tokens_in":810,"tokens_out":1603,"duration_ms":12698,"temperature":1.0,"reasoning_tokens":1523,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T17:34:47.377501+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure or compute the radiative-torque efficiency $Q_\\Gamma$ for realistic porous 1–500 µm aggregate grains at wavelengths from 0.1 µm to 3 mm; if the values are indeed 10–100 times below the compact-grain efficiencies used here, then the disruption timescales, disruption sizes, and the conclusion that RATD blocks large-grain migration must be revised. Observationally, detecting grains larger than about 100 µm inside the outflow or inner envelope of a Class 0/I protostar with bolometric luminosity near or above 100 $L_\\odot$ during an accretion burst would contradict the predicted dominance of submicron grains.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the grain transport and terminal-velocity model, and the disk-to-envelope migration scenario this paper tests."}],"review_version":1}