{"id":"feb04e25-3174-427e-8ee5-7cba847f8e99","arxiv_id":"2504.14228","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The first empirical phosphorus budget for a protoplanetary disk shows that nearly all P in HD 100546 is in refractory dust, with gas-phase P depleted by factors of 129 to 513 relative to solar.","lead":"Using ultraviolet and sub-millimeter observations, astronomers find that phosphorus in the planet-forming disk around the young star HD 100546 is mostly locked in dust grains, not floating as gas. Because phosphorus is essential for life, this is one of the first direct clues to how a key biogenic element becomes available to forming planets.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The gas-phase P upper limits are derived from PO/PN non-detections only, while the observed HCP line was not modeled; HCP could hide unconstrained gas-phase P, weakening the refractory-P conclusion.","rationale":"The reader's weakest_assumption targets the stellar mixing fraction f_acc = 98%. That concern is relevant but not the most load-bearing: a simple two-component dilution argument shows that if the initial photospheric P/H were solar, f_acc would need to be above roughly 88% to produce the observed depletion, so the quantitative depletion factor is not highly sensitive to f_acc over its plausible range. The more direct vulnerability is the gas-phase upper limit, which is the key discriminator between 'P is in dust' and 'P is absent from the gas but perhaps hidden in an unobserved carrier.' The paper itself flags that HCP was observed but not modeled, and HCP is named in the introduction as one of the dominant gas-phase P species. Without modeling HCP or justifying its neglect, the reported gas-phase P/H limits are conditional on the PO/PN chemistry of one network. This concern does not overturn the paper; the PO/PN-based limits may well be correct, and the stellar abundance provides independent evidence of depletion. But the omission should be addressed for the central claim to be fully supported. Since the reader already returned CONDITIONAL, my concern reinforces that verdict rather than changing it, so verdict_should_be is UNCHANGED.","tokens_in":18334,"tokens_out":13118,"duration_ms":130582,"concrete_test":"Re-run the DALI grid of §4.2 with the HCP J=6–5 transition included and, separately, with an alternative sticking/desorption and branching treatment that maximizes HCP and atomic-P fractions. For each model, compute the disk-integrated HCP 239.694 GHz line flux at total P/H = 2×10^-9 and 1×10^-8. If the predicted flux exceeds the observed 26.4 mK km/s upper limit for P/H = 2×10^-9, the PO/PN-based limits are confirmed as total-gas-P limits. If not, the gas-phase P budget is unconstrained and the refractory-P conclusion requires revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In §4.2, the disk-integrated gas-phase P/H upper limits are obtained from APEX non-detections of PO and PN rotational lines using DALI with a modified UMIST P-network. Table 2 also lists an HCP J=6–5 3σ upper limit of 26.4 mK km/s, but §4.2 states 'The HCP line was not included in the modelling.' This matters because §1 identifies HCP as one of the three dominant gas-phase P carriers in the interstellar medium, alongside PO and PN. The quoted limits P/H_in ≲ 2×10^-9 and P/H_out ≲ 5×10^-10 are therefore upper limits on the gas-phase P inventory only if the DALI network partitions volatile P into PO/PN efficiently enough that their lines are the dominant tracers. If a significant fraction of gas-phase P is in HCP, or in atomic P or PH3 with different excitation or chemistry, the same non-detections would permit a substantially higher total gas-phase P/H. The central claim that 'almost all' P is in a refractory reservoir depends on the contrast between the stellar total P/H and the gas-phase limit; if the gas-phase limit is not a true total-gas-P limit, that contrast could shrink by an order of magnitude. The omission of HCP is explicitly acknowledged, but no argument or test is provided that HCP is negligible. This is a load-bearing, model-dependent gap in the gas-phase constraint.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constrains the phosphorus budget of the HD 100546 protoplanetary disk by combining a photospheric P abundance measurement from HST/STIS ultraviolet spectroscopy (log(P/H)_star = -7.50(+0.23/-0.28)) with APEX sub-millimetre non-detections of PO and PN rotational lines, modeled with the DALI thermochemical code to obtain gas-phase upper limits of (P/H)_in < 2e-9 inside the 22 au dust cavity and (P/H)_out < 5e-10 outside. Interpreting the stellar abundance as the total gas-plus-dust P/H of the accretion stream via the stellar mixing fraction f_acc = 98%, they conclude that phosphorus in the disk is predominantly locked in refractory solids, consistent with apatite/schreibersite, and discuss implications for the two protoplanet candidates and for chemical habitability.","tokens_in":18587,"tokens_out":10480,"duration_ms":88185,"significance":"This would be the first empirical determination of the volatile-versus-refractory partitioning of phosphorus in a planet-forming disk, a quantity relevant to planet formation, atmospheric chemistry, and prebiotic potential. The study uses two independent datasets and an established accretion-contamination method, and the stellar depletion factor (~8) is consistent with the depletion of other refractory elements. The paper also makes falsifiable predictions for future JWST/Ariel observations. The principal caveats are the dependence of the conclusions on a single stellar-mixing model and on the completeness of the disk phosphorus chemical network, in particular the unmodeled HCP line.","major_comments":[{"comment":"The observed HCP J=6-5 transition is listed in Table 2 with a 3-sigma upper limit of 26.4 mK km/s, but the text states that 'The HCP line was not included in the modelling.' Since Section 1 identifies HCP as one of the three dominant gas-phase P carriers in the interstellar medium, the gas-phase P/H upper limits derived solely from PO and PN (P/H_in < 2e-9, P/H_out < 5e-10) are not total gas-phase phosphorus limits unless HCP is explicitly shown to be negligible. The authors should either incorporate HCP into the DALI network and model its line, or provide a quantitative argument, for example based on photodissociation rates or chemical network predictions for protoplanetary disks, that HCP carries a negligible fraction of gas-phase P. Without this, the central claim in Conclusions (iii) that phosphorus is 'strongly depleted from the gas into a refractory reservoir' is not fully supported.","section":"§4.2 / Table 2"},{"comment":"The identification of the measured photospheric abundance log(P/H)_star = -7.50 with the total P/H of the inner-disk accretion stream relies on the mixing fraction f_acc = 98% from Jermyn & Kama (2022). This is a model-dependent quantity, and if f_acc were substantially lower, the stellar abundance would not trace the inner-disk P inventory, weakening the inferred factor ~8 total-P depletion and the contrast with the gas-phase limits. The authors should include a sensitivity test or an explicit error budget for f_acc, or cite independent validation of this mixing fraction, to make the central interpretation robust.","section":"§4 (stellar mixing fraction)"}],"minor_comments":[{"comment":"The text says the gas-phase chemistry was 'probed by our ALMA observations,' but the observations were taken with APEX; please correct the instrument name.","section":"§5.1"},{"comment":"The entries Jenkins 2009a and Jenkins 2009b refer to the same publication (ApJ 700, 1299) and should be merged or distinguished consistently.","section":"References"},{"comment":"The text states that the 234 GHz setting covers PO and HCP lines and the 242 GHz setting covers PN, but the observed PN line is at 234.936 GHz while the PO and HCP lines are near 240 GHz; the assignment of lines to local oscillator settings appears reversed.","section":"§3.2 / Table 2"},{"comment":"Section 4.2 states that the disk extends to R_out = 500 au, while Table A1 lists R_out = 1000 au; please resolve this inconsistency.","section":"Table A1"},{"comment":"The sentence 'The conclusion from above that most, or all, elemental P is locked in refractory solids.' is a sentence fragment and should be rephrased.","section":"§5.2"},{"comment":"There is a typo, 'phoshporous,' in the caption; it should be 'phosphorus.'","section":"Figure 5 caption"},{"comment":"The word 'phosphorous' is used in several places (e.g., abstract, Figure 5 caption); the correct spelling for the element is 'phosphorus.'","section":"Abstract / throughout"},{"comment":"The phrase 'constant volatile P/H ratio' could be misunderstood; consider using 'constant gas-phase elemental P/H ratio' to avoid confusion between volatility and gas phase.","section":"§4.2"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses an important question and the two-strand approach is elegant. The main substantive issue is the incomplete treatment of HCP in the gas-phase analysis, which is acknowledged in the table note but not addressed in the text; this is fixable with additional modeling or a quantitative argument. The stellar mixing fraction should also be given a sensitivity analysis. I believe the paper will be suitable for publication after these points are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, here's my read. The paper gives the first empirical handle on the gas/solid partitioning of phosphorus in a planet-forming disk. That's genuinely new and worth knowing. The method is inherited – stellar accretion contamination, previously used for S and refractories – but applying it to P with APEX non-detections and a P network yields a concrete result: the star is depleted in P by ~8 relative to solar, consistent with other refractory elements, and the inner disk gas P/H is at least ~16 times lower than the total P/H in the accreting stream, implying most P is locked in solids. The two strands (stellar UV and sub-mm upper limits) are independent and internally consistent. The mineral discussion is appropriately hedged, and the second stellar fit region is a good validation.\n\nSoft spots, in order. First, the HCP issue. The observed HCP line is not modeled; the gas-phase limits come only from PO and PN. Since HCP is a known gas-phase P carrier in some environments, the limits are only as strong as the network's claim that PO/PN dominate the volatile P. The paper acknowledges the omission but gives no argument that HCP is negligible. This is a real gap; it should be fixed before publication, either by modeling HCP or by showing its abundance is small. It doesn't break the central refractory conclusion, because the stellar depletion already implies refractory behavior, but it weakens the quantitative gas-phase limits.\n\nSecond, the stellar P/H sits on narrow, blended UV lines with gf values calibrated from Vega. The authors acknowledge that systematics from other stellar parameters are not quantified. That's a moderate weakness; the internal consistency of two regions helps, but a different Teff/log g could shift the abundance.\n\nThird, the 98% accretion mixing fraction is load-bearing for equating stellar P/H to inner disk total P/H. It comes from a published model, so it's not a new step, but I'd want a sensitivity check.\n\nFourth, the disk model and updated network are only available on request. That will slow independent reproduction.\n\nBottom line: a serious, useful paper. It deserves peer review, not desk reject. I'd ask the authors to address HCP and the stellar systematics before acceptance. For a reader in astrochemistry or planet formation, it's worth a close look; I'd bring it to reading group.","headline":"First empirical P volatile/refractory budget for a planet-forming disk; solid but conditional on HCP modeling and stellar abundance systematics.","tokens_in":19183,"tokens_out":3138,"would_cite":true,"duration_ms":29765,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Phosphorus in the HD 100546 planet-forming disk is almost entirely carried by refractory dust, with gas-phase P depleted by a factor of at least 129 inside the main dust trap and at least 513 outside it.","keywords":["phosphorus","protoplanetary disks","HD 100546","refractory reservoir","dust trap","Herbig Ae/Be stars","gas-phase depletion","chemical habitability"],"falsifier":"A deeper sub-millimetre or mid-infrared search for P-bearing molecules in the HD 100546 disk that detects a line flux above the 3-sigma upper limits used here, or an independent stellar abundance re-analysis that returns a phosphorus-to-hydrogen ratio within 0.1 dex of solar, would contradict the inferred depletion and settle the claim.","tokens_in":18111,"feed_emoji":"🪐","tokens_out":12315,"duration_ms":101945,"temperature":0.7,"pith_summary":"The paper sets out to determine where elemental phosphorus sits in a planet-forming disk: in gas, in volatile ices, or in refractory dust. It uses the Herbig Ae/Be star HD 100546, whose slowly mixing radiative envelope means that accreting disk material contaminates the photosphere, so the stellar phosphorus abundance measures the total gas-plus-dust phosphorus in the inner disk. The stellar abundance is $\\log(\\mathrm{P/H})_\\star=-7.50^{+0.23}_{-0.28}$, about eight times below the solar value, while non-detections of the P-bearing molecules PO and PN in the sub-millimetre put the inner-disk gas-phase P/H below $2\\times10^{-9}$ and the outer disk below $5\\times10^{-10}$. The paper concludes that almost all elemental phosphorus in the disk is locked in a refractory (dusty) reservoir that survives the roughly 180 K temperatures at the inner edge of the 22 au dust trap, with apatite or schreibersite minerals, or possibly ammonium phosphate salts, as plausible carriers. If true, a planet's phosphorus budget is set by its solid-accretion history, with consequences for the formation and habitability of planets around such stars.","feed_headline":"Most phosphorus in this planet-forming disk rides in dust, not gas","feed_subtitle":"Gas-phase P is depleted 129-513x; photospheric P shows the missing phosphorus arrives as refractory grains.","key_machinery":"The central diagnostic is the stellar accretion-contamination method: the star's slowly mixing radiative envelope preserves freshly accreted disk material in its photosphere, so the photospheric phosphorus-to-hydrogen ratio measures the total gas-plus-dust inventory of the material that crosses the last major dust trap. The dust trap is the 22-40 au ring where large grains are stopped, which removes the refractory fraction from the accretion stream. Comparing the stellar total with the gas-phase disk abundance isolates the fraction locked in solids. The gas-phase side is carried by a two-dimensional thermo-chemical disk model with a phosphorus chemical network, whose calculated PO and PN line fluxes are compared with the sub-millimetre upper limits, in a two-zone version that lets the inner and outer disk gas-phase P/H vary independently.","core_discovery":"On its own terms, the discovery is that phosphorus in the HD 100546 disk is strongly depleted from gas into dust. The gas-phase elemental abundance inside the main dust trap is at least a factor of 129 below the solar reference, outside it at least a factor of 513 below, and the total gas-plus-dust inventory crossing into the accretion stream is depleted by a factor of about 8. Because roughly 90 percent of the refractory element mass is held back by the dust trap while volatile carriers like water ice pass through, the phosphorus carrier must remain solid at temperatures up to roughly 180 K rather than being a volatile ice such as PO, PN, or $\\mathrm{PH_3}$. The preferred solid reservoir is a mineral such as apatite or schreibersite; ammonium phosphate remains a candidate but would require a special formation mechanism to avoid depleting nitrogen on the star.","pith_inferences":["Editorial extension: if the stellar mixing fraction is close to the assumed 98 percent, a survey of Herbig Ae/Be stars with dust traps should find broadly similar phosphorus behaviour, which would turn this single-object result into a general property of transitional disks.","Editorial extension: the current non-detections leave room for a phosphorus carrier that the chemical network does not include; a targeted search for HCP or $\\mathrm{PH_3}$ lines at the same sensitivity would distinguish apatite or schreibersite carriers from an unmodelled volatile reservoir.","Editorial extension: simultaneous abundance measurements of calcium, iron, and phosphorus in the accreting photosphere could test whether phosphorus rides on the same refractory grains as these metals, and whether the inferred 90 percent grain-trapping efficiency is consistent across elements."],"forward_implications":["Giant planets forming inside and outside the 22 au dust trap will accrete very different phosphorus inventories: the inner planet candidate receives dust-poor, low-phosphorus gas, while the outer planet can gather more of the refractory reservoir.","The P/H ratio in a giant planet's envelope becomes a tracer of its refractory-element accretion history, comparable to or better than sulfur, because almost all phosphorus is delivered in solids.","Chemical habitability arguments should treat phosphorus as a primarily solid-delivered element in transitional disks, so exoplanet-atmosphere phosphorus detections must be interpreted against a refractory delivery pathway.","The same comparison between a stellar photospheric abundance and disk gas-phase upper limits can map the volatile-versus-refractory budget of other elements in other transitional disks."],"supporting_citations":[{"why":"supplies the stellar mixing fraction f_acc=98% that lets the photospheric P/H stand in for the total inner-disk P inventory.","marker":"Jermyn & Kama (2022)"},{"why":"provides the high disk accretion rate whose value drives the high mixing fraction.","marker":"Fairlamb et al. (2015)"},{"why":"supplies the adopted stellar parameters, photospheric abundances, and disk structure for HD 100546.","marker":"Kama et al. (2016)"},{"why":"establishes the thermo-chemical disk model and the stellar UV spectrum used to compute P-bearing line fluxes.","marker":"Bruderer et al. (2012)"},{"why":"sets the solar reference abundance against which all depletion factors are quoted.","marker":"Asplund et al. (2009)"},{"why":"fixes the disk temperature and density structure, including the roughly 180 K temperature at the dust-ring inner edge.","marker":"Keyte et al. (2023)"},{"why":"demonstrates the same accretion-contamination diagnostic for sulfur, providing the methodological template.","marker":"Kama et al. (2019)"},{"why":"supplies the binding energies for PO and PN used in the chemical network.","marker":"Piacentino & Öberg (2022)"},{"why":"gives the equilibrium condensation temperatures used to argue that phosphorus condenses into schreibersite and phosphates.","marker":"Lodders (2003)"}],"fun_headline_variants":["Phosphorus in HD 100546 disk: gas depleted, dust holds the load","Gas-phase phosphorus depleted 129-513x in this planet-forming disk","Refractory phosphorus rules: dust traps the element, gas runs low","In HD 100546, missing phosphorus is in dust, not volatile ice","Phosphorus budget: HD 100546 disk gas depleted, dust carries the rest"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the phosphorus-to-hydrogen ratio measured on the star's surface equals the total phosphorus-to-hydrogen ratio of the gas and dust that cross the 22 au dust trap and fall onto the star, which requires a roughly 98 percent mixing fraction for freshly accreted material and no chemical sorting between gas and dust during accretion.","fun_headline_variants_meta":{"raw":{"variants":["Phosphorus in HD 100546 disk: gas depleted, dust holds the load","Gas-phase phosphorus depleted 129-513x in this planet-forming disk","Refractory phosphorus rules: dust traps the element, gas runs low","In HD 100546, missing phosphorus is in dust, not volatile ice","Phosphorus budget: HD 100546 disk gas depleted, dust carries the rest"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000747,"raw_usage":{"total_tokens":3357,"prompt_tokens":999,"completion_tokens":2358,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":615,"completion_tokens_details":{"reasoning_tokens":2255}},"tokens_in":615,"tokens_out":2358,"duration_ms":14830,"temperature":1.0,"reasoning_tokens":2255,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:54:12.658934+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A deeper sub-millimetre or mid-infrared search for P-bearing molecules in the HD 100546 disk that detects a line flux above the 3-sigma upper limits used here, or an independent stellar abundance re-analysis that returns a phosphorus-to-hydrogen ratio within 0.1 dex of solar, would contradict the inferred depletion and settle the claim.","supporting_citations":[{"cited_title":"S., Kama M., 2022, @doi [Research Notes of the American Astronomical Society] 10.3847/2515-5172/ac79b9 , https://ui.adsabs.harvard.edu/abs/2022RNAAS...6..131J 6, 131","cited_arxiv_id":null,"evidence_quote":"supplies the stellar mixing fraction f_acc=98% that lets the photospheric P/H stand in for the total inner-disk P inventory."},{"cited_title":"L., \\\"O berg K","cited_arxiv_id":null,"evidence_quote":"supplies the binding energies for PO and PN used in the chemical network."}],"review_version":1}