{"id":"2d3c3e6f-a835-428d-a560-85efe377f4b3","arxiv_id":"2502.02411","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"In the CEJSN scenario, the accretion disk can entrain up to ~0.01-0.03 M_sun of neutron star crust material via Kelvin-Helmholtz mixing, enhancing the r-process yield.","lead":"This paper argues that in a proposed site for making heavy elements, the common envelope jets supernova scenario, the dense disk around a neutron star can scrape off and mix in neutron-rich crust material, boosting the amount of r-process elements ejected. The mechanism could help explain where the heaviest elements in the universe come from, but it rests on several untested assumptions about how a solid neutron star crust behaves.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The KHI mixing analysis treats the NS crust as a fluid (Eqs. 11-14) while conceding it is solid; without a neutrino-cooled melting timescale, the 0.01 M_sun entrainment claim is not secured.","rationale":"The paper's argument chain is: (i) the disk density at the NS surface (Eq. 1) reaches ~10^12 g/cm^3; (ii) the disk penetrates the crust to the depth where densities match; (iii) KHI mixes deeper crust material, giving q_rho up to ~86 and Mcr,mix ~ 0.01 M_sun (Eq. 17); (iv) jets carry this material out, raising the r-process yield. The weakest link is (iii): the KHI derivation (Eqs. 9-14) is for a fluid-fluid shear layer, whereas the NS crust is a solid lattice at these densities. The paper itself concedes the crust is solid and does not provide a quantitative melting calculation under neutrino cooling. If the crust remains solid, the interface is elastic, the KHI criterion does not apply, and the entrained mass would be set by ablation rather than KHI, so the 0.01 M_sun enhancement is unsupported. The reader's verdict (CONDITIONAL) already captures this appropriately: the claim is plausible and well-posed, but requires verification of the fluidization assumption. My proposed test would settle the matter by evaluating whether the boundary layer actually reaches T_SL ~ 10^10 K; a null result would reduce the yield to negligible, while a positive result would validate the fluid assumption. Therefore I see no reason to change the reader's verdict, and I agree that the solid-crust fluidization is the weakest assumption.","tokens_in":12470,"tokens_out":12134,"duration_ms":109448,"concrete_test":"Compute the steady-state temperature of the disk-crust boundary layer by balancing the kinetic-energy dissipation rate per unit area (approximately rho_d v_Kep^3 (Delta_qv)^2, using values from Eqs. (1), (2), and (8)) against the neutrino emissivity from modified Urca and pair processes at rho ~ 10^12-10^13 g/cm^3. If the equilibrium temperature remains below the melting temperature T_SL ~ 10^10 K (e.g., from Carreau et al. 2020), the crust does not melt and the KHI entrainment of Section 3 fails; if T exceeds T_SL over a layer thicker than the unstable wavelength, the fluid assumption is justified. Alternatively, a 2D boundary-layer simulation with an elastic-plastic crust equation of state can directly test whether the crust melts or is ablated before KHI grows.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim, entrainment of up to ~0.01 M_sun of neutron-rich crust material, rests on the Kelvin-Helmholtz instability (KHI) analysis in Section 3. The KHI condition (Eqs. 7-9) and the mixing-depth derivation (Eqs. 11-14) assume the crust behaves as a fluid with a continuous density gradient, but at densities of 10^12-10^13 g/cm^3 the neutron star crust is a Coulomb solid. The paper explicitly acknowledges this in Section 3 ('which is a solid') but supports fluidization with only a qualitative statement: a simple calculation without neutrino cooling gives T > 10^11 K, and 'nonetheless, the dissipated energy can turn the solid mixed crust material to liquid.' No quantitative melting timescale, dissipation rate, neutrino cooling rate, or steady-state temperature is provided. If the crust remains solid, the interface is elastic, the KHI criterion (7) derived for two fluids does not apply, and the shear is accommodated in the fluid boundary layer rather than by entraining crust material. The entrained mass would then be limited to whatever can be ablated, collapsing the 0.01 M_sun yield. This is load-bearing because both Eq. (15) and the mass estimate Eq. (17) depend on the fluid assumption. The paper also assumes the mixed material is expelled in the jets (Section 4) without modeling, but the fluidization of the crust is the more fundamental unverified step.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes a new mechanism in the common-envelope jets supernova (CEJSN) r-process scenario: the high-density accretion disk around a neutron star (NS) that enters a massive stellar core can penetrate the NS crust, and Kelvin-Helmholtz instability (KHI) can mix neutron-rich crust material into the inner disk. Using scaling relations for the disk density (Eq. 1), a fitted crust pressure-density relation (Eq. 4), and a KHI criterion (Eqs. 7-15), the author estimates that up to ~0.01 M_sun of original NS crust material can be entrained, and possibly ~0.03 M_sun with replenishment, with this material then carried out by jets and contributing to r-process nucleosynthesis. The paper explicitly frames the result as strengthening, not replacing, the CEJSN r-process scenario.","tokens_in":12722,"tokens_out":10325,"duration_ms":93545,"significance":"If the mechanism operates, it would provide a new channel for supplying neutron-rich material to the r-process in CEJSNe and would help explain large per-event r-process yields. The paper is commendably transparent: it uses published scaling relations, fits no free parameters to the claimed yield, explicitly acknowledges key assumptions, and calls for numerical simulations. The algebraic derivation of the KHI condition (Eqs. 11-14) is internally consistent, and the manuscript is clearly written. However, the central quantitative claim rests on two unverified physical steps: treating the solid NS crust as a fluid for KHI, and assuming that the mixed boundary-layer material is ejected in jets rather than accreted. Both are acknowledged but not quantitatively supported, so the claimed 0.01-0.03 M_sun yield remains conditional.","major_comments":[{"comment":"The KHI analysis treats the NS crust as a fluid with a continuous density gradient, but at densities of rho_cr ~ 10^12-10^13 g cm^-3 the outer and inner crust is a Coulomb solid. The paper itself concedes this ('which is a solid') and supports fluidization only with the qualitative statement that a calculation without neutrino cooling gives T > 10^11 K and that 'nonetheless, the dissipated energy can turn the solid mixed crust material to liquid.' No dissipation rate, neutrino-cooling timescale, melting timescale, or steady-state temperature is provided. This is load-bearing because the KHI criterion in Eq. (7) is derived for two fluids and does not apply to an elastic solid; if the crust remains solid, the shear is accommodated in the disk's boundary layer and the entrained mass in Eq. (17) is not secured. The manuscript should provide a quantitative estimate of the local heating and cooling balance in the boundary layer, or demonstrate that KHI can operate on a partially molten or elastic interface, before the 0.01 M_sun claim is accepted.","section":"Section 3, Eqs. (7)-(15) and Eq. (17)"},{"comment":"The final r-process yield also depends on the assumption that the mixed crust-disk material is ejected in the jets. The text states 'The proposed scenario assumes that the crust-disk mixed material is ejected from the disk and is part of the jets' material,' but no model or estimate of the ejection fraction from the boundary layer is given. A boundary layer can accrete rather than eject, and the mass mixed into the disk (Eq. 17) need not be launched. Since the abstract's 0.01-0.03 M_sun yield depends on this step, the manuscript should either provide a quantitative basis for the ejection fraction or clearly label the 0.01-0.03 M_sun value as an upper limit conditional on efficient jet launching from the mixed layer.","section":"Section 4"}],"minor_comments":[{"comment":"The phrase 'at a density of 10^12 g cm^-2' should read '10^12 g cm^-3'.","section":"Section 3, after Eq. (17)"},{"comment":"The text '0.05 M_sun yr^-1' appears to be a typo; it should be '0.05 M_sun s^-1' to be consistent with Eq. (1) and the rest of the paper.","section":"Section 4, first paragraph"},{"comment":"The expression '2GMMN/(c^2 RNS)' contains a typographical error ('MMN'); it should be 'GMNS'.","section":"Section 2, relativistic-effects paragraph"},{"comment":"The relation between d(rho_cr)/rho_cr and (rho_cr - rho_d)/((rho_cr + rho_d)/2) is easy to misread because of missing parentheses; adding an explicit definition or parentheses would improve clarity.","section":"Section 3, Eqs. (12)-(13)"},{"comment":"The assumption that the mixing depth is of the order of the KHI wavelength is ad hoc, and the subsequent statement that a 'more accurate treatment should yield a deeper mixing length' is not substantiated; either justify this claim or remove it.","section":"Section 3, mixing-depth assumption"},{"comment":"The abstract's '0.01-0.03 M_sun' conflates the newly entrained mass with the total r-process yield from earlier CEJSN models; it would be clearer to state explicitly that the entrained mass is added to an already 0.01-0.03 M_sun scenario.","section":"Abstract and Section 4"}],"recommendation":"major_revision","confidential_remarks":"This is a speculative but clearly written scenario paper. The main barrier to acceptance is not the algebra but the physical applicability of fluid KHI to a solid crust; the author should be encouraged to add a quantitative melting argument or reframe the central claim as conditional. The manuscript's heavy reliance on the author's own previous scaling relations is acceptable in this context, but the novelty relative to prior CEJSN r-process papers could be stated more sharply. The paper fits the scope of a fast-turnaround astrophysics journal, provided the load-bearing assumptions are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hi X,\n\nQuick take on Soker's arXiv:2502.02411. The paper argues that in the common-envelope jets supernova (CEJSN) r-process scenario, the high-density accretion disk can penetrate the neutron star crust, and the Kelvin-Helmholtz instability can mix crust material into the disk, adding up to ~0.01 Msun of neutron-rich mass to the jets and raising the per-event r-process yield to 0.01-0.03 Msun. That mechanism is new for this scenario, and the algebra is internally consistent. The paper is appropriately hedged: it claims a strengthening of the CEJSN r-process site, not the main site, and it explicitly lists most of its own assumptions.\n\nWhat it does well: The derivation is transparent and checkable. It addresses magnetic field suppression of KHI, it gives a simple mass estimate from a crust pressure-density fit, and it highlights the dynamical-range problem for future simulations. It is a clean order-of-magnitude plausibility argument.\n\nThe soft spots are real, and the stress-test note puts its finger on the load-bearing one. The KHI mixing calculation treats the NS crust as a fluid, but at 10^12-10^13 g/cm^3 the crust is a Coulomb solid. The paper acknowledges this in Section 3 and asserts that dissipated kinetic energy can melt it, but it gives no quantitative melting timescale, no neutrino cooling estimate, and no steady-state temperature. If the crust stays solid, the two-fluid KHI criterion does not apply and the 0.01 Msun number collapses. The disk density at the NS surface is also extrapolated from Grichener & Soker (2019a) with a scaling law; if the true density is lower, the entrained mass drops. And the assumption that the mixed material ends up in the jets is asserted, not modeled. The paper flags the first two as limitations, but they are still load-bearing.\n\nNone of this kills the paper as a plausibility study. The argument is well-posed, the caveats are mostly explicit, and the author asks for the right follow-up (high-resolution simulations of the boundary layer). But the central yield should be read as an upper limit, not a benchmark prediction.\n\nWho gets value: people comparing r-process site candidates, and anyone modeling accretion onto NSs in common envelope events. It is a reasonable submission for a refereed journal, and I would send it to a referee. The referee should push for either a melting timescale estimate or a reduction in the claimed mass. I would not cite it as evidence for a specific yield, but I would cite it as a proposed mechanism.","headline":"Plausible order-of-magnitude argument that CEJSN disks can entrain NS crust, but the crust-fluidization assumption is unquantified, so the yield is an upper limit.","tokens_in":13339,"tokens_out":2650,"would_cite":false,"duration_ms":26398,"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":"In the common-envelope jets supernova r-process scenario, a dense accretion disk around a neutron star can penetrate the star's crust, and Kelvin-Helmholtz mixing feeds neutron-rich crust material into the jets, raising r-process yields…","keywords":["r-process nucleosynthesis","common envelope jets supernova","neutron star crust","Kelvin-Helmholtz instability","accretion disk","jets","heavy element nucleosynthesis","neutron star accretion"],"falsifier":"A three-dimensional simulation that resolves the mixing layer (scale height $H_\\rho \\approx 0.01 R_{\\rm NS}$), includes neutrino cooling and a realistic solid-crust equation of state, and shows that the crust stays solid and the Kelvin-Helmholtz instability stops at the surface would falsify the $0.01\\,M_\\odot$ entrainment claim; a softer observable falsifier is CEJSN candidates whose r-process ejecta mass comes out below $0.01\\,M_\\odot$.","tokens_in":12138,"feed_emoji":"💥","tokens_out":11528,"duration_ms":95053,"temperature":0.7,"pith_summary":"In the common-envelope jets supernova (CEJSN) r-process scenario, a neutron star plunges into the core of a massive evolved star, accretes at about 0.05 solar masses per second, and launches jets. This paper argues that the resulting high-density accretion disk does not stop at the neutron star's surface: it can penetrate the outer crust down to the density where disk and crust match, and the Kelvin-Helmholtz instability then mixes neutron-rich crust material from layers tens of times denser into the disk. The jets carry this material out, so each event could eject roughly 0.01-0.03 solar masses of r-process elements. The author presents this as strengthening the CEJSN scenario as one of several r-process sites, not as the main site.","feed_headline":"Accretion disk can shred a neutron star's crust into r-process jets","feed_subtitle":"Mixing neutron-rich crust into the jets could raise r-process yields to 0.01-0.03 solar masses per event.","key_machinery":"The load-bearing mechanism is the Kelvin-Helmholtz instability acting at the interface between the Keplerian accretion disk and the much slower neutron-star crust. The key derived scale is the density ratio $q_\\rho \\equiv \\rho_{\\rm cr}/\\rho_d$; the instability condition reduces to $q_\\rho \\lesssim 86 (\\Delta q_v/0.5)^2$, which sets how deep the mixing reaches. Combining that with the crust mass-density relation $M_{\\rm cr}(>\\rho) \\approx 10^{-4} (\\rho/10^{12}\\,\\mathrm{g\\,cm^{-3}})\\,M_\\odot$ gives the entrained mass, up to about $0.01\\,M_\\odot$. A supporting ingredient is the claim that magnetic fields cannot suppress the instability unless $B \\gtrsim 3\\times 10^{16}\\,\\mathrm{G}$, far above typical neutron-star fields.","core_discovery":"The paper's central claim is that the accretion disk in the CEJSN r-process scenario is dense enough to eat into the neutron star's crust, and that the shear between the disk and the crust drives Kelvin-Helmholtz mixing that pulls neutron-rich material from the deep inner crust into the disk. Using the disk density $\\rho_d \\approx 10^{12}\\,\\mathrm{g\\,cm^{-3}}$ from earlier scalings, the disk penetrates to the equal-density level, and the Kelvin-Helmholtz condition (with shear $\\Delta q_v \\sim 0.5$) keeps the instability alive up to density ratios $q_\\rho \\approx 86$, i.e. densities of a few $\\times 10^{13}\\,\\mathrm{g\\,cm^{-3}}$. The mass above that depth is about $0.01\\,M_\\odot$ of original cold crust; because accretion replenishes the crust on a timescale of about 0.002 s while the event lasts 10-100 s, the total entrained mass could reach about $0.03\\,M_\\odot$. Entraining this neutron-rich material lowers the electron fraction and enlarges the neutron reservoir, so the estimated r-process ejecta per event rises to $0.01$-$0.03\\,M_\\odot$. The paper does not claim CEJSNe are the main r-process site; it claims they are a contributing site, consistent with evidence that two or more sites are needed.","pith_inferences":["If crust entrainment operates, the abundance pattern of CEJSN r-process ejecta should carry a signature of neutron-star crust composition that differs from neutron-star merger ejecta; comparing predicted iridium-to-europium and lanthanide fractions with observed metal-poor stars could test this beyond what the paper computes.","The same shear-mixing argument could be applied to other configurations where a dense disk forms on a compact object with a solid crust, such as accreting white dwarfs; the paper cites numerical work there without developing the analogy.","Magnetic fields far below about 10^16 gauss do not suppress the instability in this estimate, but a CEJSN system with a magnetar-strength field in a stabilizing geometry could cut the entrained mass sharply; current r-process-rich star observations cannot yet distinguish this case."],"forward_implications":["Each CEJSN r-process event can eject 0.01-0.03 solar masses of r-process elements, making rare early-Universe events viable heavy-element sources.","Neutron-rich crust material guarantees a low electron fraction in the jet-launching region even for accretion rates below the previously quoted threshold.","The heaviest r-process nuclei, including third-peak and actinide tracers, are likely products of the inner, crust-enriched disk because that material has the lowest electron fraction and highest density.","The scenario supports the claim that at least two r-process sites contribute to galactic nucleosynthesis, with CEJSNe acting on short delay times in the young Universe.","Future simulations must resolve density scale heights of about 0.01 neutron-star radii and include both the disk and the neutron star, a dynamical range of about three orders of magnitude."],"supporting_citations":[{"why":"Supplies the CEJSN scenario, the disk density and scale-height scalings, and the 0.05 solar masses per second accretion baseline used throughout.","marker":"Grichener & Soker (2019a)"},{"why":"Provides the neutron-star crust density profile and the pressure-density fit used to compute entrained crust mass.","marker":"Chamel & Haensel (2008)"},{"why":"Gives the mass-above-radius approximation used to convert mixing depth into entrained mass and the magnetar-flare yield used for comparison.","marker":"Cehula et al. (2024)"},{"why":"Gives the magnetospheric truncation radius expression used to rule out disk truncation by typical neutron-star magnetic fields.","marker":"Long et al. (2005)"},{"why":"Provides a lower-accretion-rate formula for the mixing depth used as a consistency check on the few times 10^13 grams per cubic centimeter penetration estimate.","marker":"Fujimoto (1993)"},{"why":"Sets the accretion-rate threshold for electron fraction below 0.15 that the disk already exceeds.","marker":"Kohri et al. (2005)"},{"why":"Lowers the strong-r-process accretion-rate limit, making the CEJSN disk comfortably r-process-favorable.","marker":"Siegel et al. (2019)"},{"why":"Supplies the magnetic-field criterion used to argue that the Kelvin-Helmholtz instability is not suppressed for realistic neutron-star fields.","marker":"Prajapati et al. (2009)"},{"why":"Demonstrates experimentally that the Kelvin-Helmholtz instability mixes stratified shear layers.","marker":"Strang & Fernando (2001)"}],"fun_headline_variants":["Disk crunches neutron star crust to power r-process jets","Accretion disk shears off neutron star crust for r-process","Neutron-rich crust mixed into jets boosts r-process yields","Disk eats neutron star crust, enriches r-process jets","Crust material from neutron star enriches r-process jets"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation treats the neutron-star crust as a continuous fluid with a smooth density gradient, but the outer and inner crust are solid lattices at the relevant densities; the paper only asserts that dissipated kinetic energy can liquefy the mixed material and does not show that melting outruns neutrino cooling.","fun_headline_variants_meta":{"raw":{"variants":["Disk crunches neutron star crust to power r-process jets","Accretion disk shears off neutron star crust for r-process","Neutron-rich crust mixed into jets boosts r-process yields","Disk eats neutron star crust, enriches r-process jets","Crust material from neutron star enriches r-process jets"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000549,"raw_usage":{"total_tokens":2698,"prompt_tokens":1101,"completion_tokens":1597,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":717,"completion_tokens_details":{"reasoning_tokens":1513}},"tokens_in":717,"tokens_out":1597,"duration_ms":12153,"temperature":1.0,"reasoning_tokens":1513,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T12:13:13.342756+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A three-dimensional simulation that resolves the mixing layer (scale height $H_\\rho \\approx 0.01 R_{\\rm NS}$), includes neutrino cooling and a realistic solid-crust equation of state, and shows that the crust stays solid and the Kelvin-Helmholtz instability stops at the surface would falsify the $0.01\\,M_\\odot$ entrainment claim; a softer observable falsifier is CEJSN candidates whose r-process ejecta mass comes out below $0.01\\,M_\\odot$.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides a lower-accretion-rate formula for the mixing depth used as a consistency check on the few times 10^13 grams per cubic centimeter penetration estimate."},{"cited_title":"K., & Chhajlani, R","cited_arxiv_id":null,"evidence_quote":"Supplies the magnetic-field criterion used to argue that the Kelvin-Helmholtz instability is not suppressed for realistic neutron-star fields."},{"cited_title":"J., & Fernando, H","cited_arxiv_id":null,"evidence_quote":"Demonstrates experimentally that the Kelvin-Helmholtz instability mixes stratified shear layers."}],"review_version":1}