{"id":"335e3b2b-97f4-4795-bafa-5492e94fdeca","arxiv_id":"2507.10858","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Kinetic coefficients for sublimation and deposition at a solid sphere in a binary vapor-noncondensable gas mixture are computed, revealing strong dependence on rarefaction, composition, and intermolecular potential.","lead":"This paper computes how the rate of sublimation or deposition at a solid argon sphere changes when helium is present in the surrounding vapor, using a kinetic model of molecular collisions. The result is a set of mass and heat transfer coefficients that could be used to design freeze-drying, vapor deposition, and aerosol processes with spherical particles.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The sign-change and >100% deviation claims for ΛPT rely on the unverified complete-condensation boundary condition (Eqs. 26–27); a partial sticking coefficient could shift or eliminate the zero-crossing.","rationale":"The reader's weakest assumption is correctly identified: the complete phase transition boundary condition, Eqs. (26)–(27), is the single most load-bearing premise for the central quantitative claims. The tables of ΛPP, ΛPT, ΛTT are moments of the solution to that boundary-value problem, so any departure from σ=1 or from Maxwellian emission changes the reported numbers. This is not merely an outside-consensus disagreement; it is an explicit modeling choice that the paper makes without sensitivity analysis. The reader's verdict of CONDITIONAL is therefore appropriate. I do not see an internal inconsistency that would overturn the calculation under the stated assumptions: the free-molecular coefficients (62)–(64) reproduce the table values, the Onsager symmetry is used as a check, and the governing equations are standard for this class of problems. The sign typo in Eq. (47) is an error in an auxiliary density moment and does not directly enter the kinetic coefficients; the numerical convergence data are absent but the paper claims 0.1% accuracy and reciprocity checks. Consequently, the central qualitative finding should stand as a conditional result, with the boundary-condition sensitivity as the key uncertainty to resolve before applying the coefficients to real argon–helium systems.","tokens_in":22944,"tokens_out":26121,"duration_ms":318412,"concrete_test":"Modify the species-2 boundary condition to a partial condensation model, e.g. h(P)_2 = σ + (1−σ)h_reflected and h(T)_2 = σ(c2^2−5/2) + (1−σ)h_reflected, with h_reflected a diffuse Maxwellian at Ts, and recompute the free-molecular limit analytically and the δ=0.1, 1, 10 cases numerically for σ=0.5 and 0.1 at C0=0.1, 0.5, 0.9. If ΛPT still changes sign with δ and still deviates by more than 100% from its free-molecular value for σ<1, the boundary-condition concern is not load-bearing for the central qualitative claim; if the zero-crossing shifts or disappears, the conclusions must be restated as conditional on σ=1.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central results are linear functionals of the boundary-value problem posed in Secs. 3–4. For the vapor species, Eqs. (26)–(27) impose complete phase transition: every incident argon/krypton atom is absorbed, and the surface emits a Maxwellian at Ts and p2s. This is an assumption, not a derived consequence of the kinetic equation, and it enters every entry in Tables 2–4. If the true condensation/sticking coefficient σ is below unity, or if desorption is not exactly Maxwellian, all three coefficients ΛPP, ΛPT, ΛTT change. The cross coefficient ΛPT is the most fragile: it is often an order of magnitude smaller than the diagonal coefficients and, in the tables, its sign change occurs at values very close to zero (e.g., Table 2, HS, C0=0.5, δ=5: ΛPT=0.0040 versus the free-molecular value −0.0523). A partial-condensation boundary condition does not merely rescale the coefficients uniformly; it changes the balance between emitted and reflected molecules, so the zero of ΛPT could move or disappear. The paper offers no sensitivity analysis and no experimental benchmark for σ, so the headline quantitative conclusions are not yet established for real argon–helium conditions. The numerical convergence issue noted by the reader is real but subsidiary: even with perfect numerics, the boundary-condition assumption remains load-bearing.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies sublimation/deposition at a solid argon sphere in a helium-argon mixture, using the linearized Boltzmann equation with the McCormack collision model. The authors compute the Onsager kinetic coefficients ΛPP, ΛPT, ΛTT that determine mass and energy flow rates at the interface, over a range of rarefaction parameter δ, molar fractions C0, and two temperatures (50 K and 70 K), using both hard-sphere and ab-initio potentials. They also give analytic free-molecular results, solve the transitional regime numerically by a discrete velocity method, report flow fields and temperature/pressure jump behavior, and include examples of mass and energy flow rates. The main claims are that the cross coefficient ΛPT is highly sensitive to rarefaction, can deviate by more than 100% from its free-molecular value, and can change sign with δ.","tokens_in":23268,"tokens_out":9047,"duration_ms":111483,"significance":"If the results hold, the paper provides a useful tabulation of spherical-geometry kinetic coefficients for sublimation/deposition in a binary mixture, based on ab-initio interaction data rather than fitted transport parameters. The free-molecular formulas are explicit and check against the tables, and the numerical reciprocity check ΛPT=ΛTP is a valuable internal consistency test. The comparison of hard-sphere and ab-initio potentials, and of argon and krypton, gives a concrete picture of where the interaction model matters. However, the central quantitative claims are conditional on the assumed complete phase transition boundary condition, and no sensitivity analysis for that assumption is provided; this is the main scientific weakness.","major_comments":[{"comment":"Equation (47) for ν(P)_2 has the wrong sign. With h(P)_2=1 on the cone 0≤θ≤θ0 and zero outside (Eq. (43)), the angular average is (1−cos θ0)/2 = (1−sqrt(1−(r0/r)^2))/2, which is positive and tends to 0 as r→∞. The printed expression −1/2 − sqrt(1−(r0/r)^2)/2 is negative everywhere and tends to −1 at infinity; it also contradicts the positive ν(P)_2 profiles in Figure 1. Please correct Eq. (47) and check the companion moments (48)–(51) for consistency.","section":"Sec. 6, Eq. (47)"},{"comment":"The complete phase transition boundary condition (26)–(27) is an assumption that enters every entry of Tables 2–4. The headline conclusions—that ΛPT can deviate by more than 100% from its free-molecular value and can change sign with δ—are statements about the solution of this particular boundary-value problem. A partial sticking coefficient or non-Maxwellian emission could shift or remove the zero of ΛPT. The paper offers no sensitivity analysis and no experimental benchmark for the sticking coefficient. Please add a parametric study with a partial-condensation boundary condition (e.g., a Maxwell-type condition) or, at minimum, a quantitative discussion of how the tabulated coefficients would change; otherwise the conclusions should be presented as conditional on the complete-transition assumption.","section":"Sec. 4, Eqs. (26)–(27); Tables 2–4"},{"comment":"The claimed numerical error of 0.1% for the kinetic coefficients is stated but not demonstrated. The sign change of ΛPT occurs at values as small as 0.0040 (Table 2, HS, C0=0.5, δ=5), and the reported deviations from free-molecular values exceed 100%. The accuracy statement needs to be supported by a convergence study in Nr, Nθ, Nc and rmax, or by a supplementary table reporting the grid-variation results. Please include this information so that the sign changes in the tables can be assessed as genuine physical effects rather than numerical artifacts.","section":"Sec. 7; Tables 2–4"}],"minor_comments":[{"comment":"There is a typo: “emmited” should be “emitted”.","section":"Sec. 4"},{"comment":"The sentence about rmax (“while the distance rmax was set as δ varied so that the increment in the radial distance ∆r ∼ 10−3”) is unclear and should be rewritten; please state the actual rmax values used for each δ.","section":"Sec. 7"},{"comment":"The text after Eq. (85) swaps the symbols for pressure and mole fraction deviations: ξ is defined as the pressure deviation in Eq. (84) and ζ as the mole fraction deviation in Eq. (85), but the paragraph says “The local mole fraction deviation, ξ, is always negative” and “the pressure deviation, ζ, is always positive.” This should be corrected in the text and checked against the panels in Figures 7–9.","section":"Sec. 8.7"},{"comment":"Please state in the captions that the rows with δ=0 are obtained from the analytic free-molecular formulas, not from the discrete velocity method.","section":"Tables 2–4"},{"comment":"The choice γα=p0α/µ0α is introduced by analogy with the single-gas Shakhov model; a brief justification or a citation to a previous mixture validation of this specific choice would help the reader assess its effect on the transitional-regime results.","section":"Appendix A, Eq. (A.9)"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the numerical methodology is solid, but the sign error in Eq. (47), the symbol swap in Sec. 8.7, and the unsupported 0.1% error claim should be fixed. The main scientific risk is the absence of sensitivity analysis for the complete phase transition boundary condition, which underpins the paper's headline claims about ΛPT. I do not see grounds for rejection, but the revision should be substantive rather than cosmetic."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague — this is a competent, incremental but useful paper. The new content is the spherical geometry for sublimation/deposition of a vapor species in a non-condensable background gas within the McCormack model, with ab initio potentials, plus an analytic free-molecular solution for the mixture. If you work in rarefied gas dynamics or freeze-drying modeling, the tabulated Onsager coefficients (ΛPP, ΛPT, ΛTT) across rarefaction, concentration, and potential are a solid addition; the reciprocity check ΛPT=ΛTP is a good numerical control.\n\nThe paper is honest about its method. The numerical approach is standard (discrete velocity method), the parameter space is well chosen, and the comparison between hard-sphere and ab initio potentials is meaningful. The free-molecular limits are a useful benchmark.\n\nThe soft spots are real but manageable. First, Eq. (47) for ν(P)_2 in the free-molecular regime is wrong as printed: it gives -1 at infinity instead of 0, and it is negative everywhere, contradicting the positive density deviations shown in Figures 1–3. It looks like a sign error in the printed formula; the tables of Λ do not depend on it directly, but the paper must be corrected. Second, the convergence study is described only in words (error 0.1%, grid variation) with no table of grid refinement; given the accuracy claim, that should be shown or at least put in a supplementary table. Third, and more substantive: the boundary condition of complete phase transition for the vapor (sticking coefficient 1, Maxwellian emission at Ts, p2s) is an idealization that is load-bearing for every coefficient. Real argon/helium surfaces may not stick perfectly, and the paper gives no sensitivity analysis. That does not make the results wrong as a kinetic-theory calculation, but it does mean the quantitative conclusions, especially the sign change of ΛPT near zero, are conditional on that assumption. The authors should state this more clearly as a limitation, not just as an assumption.\n\nFor a reader in this niche, the paper is worth engaging with. I would send it to a competent referee — it will be a useful contribution after minor revisions. I wouldn't cite it in my own work, but I can see colleagues in kinetic theory or aerosol physics doing so.","headline":"Solid incremental study of sublimation at a sphere in a binary mixture; useful tables, but Eq. (47) has a sign error and the complete-condensation boundary condition deserves a caveat.","tokens_in":23738,"tokens_out":3378,"would_cite":false,"duration_ms":40788,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["47.45.-n","51.10.+y","05.70.Ln"],"model":"deepseek-v4-flash","headline":"This paper computes the three Onsager kinetic coefficients for sublimation and deposition at a solid argon sphere in a helium background gas, and shows that the cross-coupling term can change sign and deviate by more than 100% from its…","keywords":["sublimation","deposition","mass transfer","Onsager coefficients","linearized Boltzmann equation","McCormack model","binary gas mixture","ab initio potential"],"falsifier":"Measure the steady mass flow from a small argon sphere suspended in a helium-argon mixture at $T_0=50$ K with $C_0=0.5$ while sweeping the pressure so that the rarefaction parameter runs from about 0.1 to 10, and applying a small temperature difference of about 0.1 K. The predicted sign change of $\\Lambda_{PT}$ would appear as a reversal of the net argon flow direction relative to the temperature gradient, and the particular value of $\\delta$ at which the reversal occurs would distinguish the hard-sphere from the ab initio potential prediction.","tokens_in":22752,"feed_emoji":"🧊","tokens_out":8019,"duration_ms":85129,"temperature":0.7,"pith_summary":"This paper computes the three Onsager kinetic coefficients $\\Lambda_{PP}$, $\\Lambda_{PT}$ and $\\Lambda_{TT}$ that determine the mass and energy flow rates when a solid argon sphere sublimes or deposits in a mixture with non-condensable helium. The coefficients are obtained from the linearized Boltzmann equation with the McCormack collision model, for both hard-sphere and ab initio interatomic potentials, across free-molecular, transitional, and near-continuum regimes. The central finding is that the cross coefficient $\\Lambda_{PT}$, which couples the temperature difference to the mass flow, is far more sensitive than the diagonal ones: it can deviate by more than 100% from its free-molecular value, depends strongly on the interaction potential and the helium fraction, and can even change sign as the rarefaction parameter, temperature, or species mass ratio varies. These tabulated coefficients let engineers predict the interface flow rates for arbitrary small driving forces without re-solving the kinetic equation.","feed_headline":"Cross term flips sign as argon sphere sublimates in helium","feed_subtitle":"Three tabulated kinetic coefficients let engineers compute mass and energy flow for any small driving force.","key_machinery":"The load-bearing object is the symmetric Onsager matrix built from the fluxes $J_P=(1-C_0)u_2$ and $J_T=C_0q_1+(1-C_0)q_2$, evaluated with the linearized Boltzmann equation for the binary mixture, where $u_2$ is the mean bulk velocity of the sublimating species and $q_1$, $q_2$ are the heat fluxes of the two species. The collision integral is replaced by the McCormack model, which respects conservation laws and the H-theorem and gives correct mixture transport coefficients, while the interatomic forces enter through the $\\Omega$-integrals for either hard spheres or ab initio potentials. To handle the discontinuity of the distribution function at the surface of a convex body, the solution is split into an analytically known free-streaming part and a numerical part solved with the discrete velocity method, and the reciprocity relation $\\Lambda_{PT}=\\Lambda_{TP}$ is used as a check of the numerical error. The free-molecular regime is solved analytically and serves as the reference for all deviations.","core_discovery":"For sublimation and deposition at a solid argon sphere surrounded by a helium-argon mixture, the interface mass and energy flow rates are linear functions of two thermodynamic forces, the vapor pressure difference and the temperature difference, with the proportionality encoded in a symmetric matrix of three kinetic coefficients: $\\Lambda_{PP}$, $\\Lambda_{PT}(=\\Lambda_{TP})$ and $\\Lambda_{TT}$. The paper establishes how these coefficients behave: $\\Lambda_{PP}$ and $\\Lambda_{TT}$ stay positive, while $\\Lambda_{PT}$ can be positive or negative and its sign depends on the rarefaction parameter $\\delta$, the helium molar fraction $C_0$, the temperature, the interatomic potential, and the mass ratio of the two gases. Unlike the planar case, where $\\Lambda_{PP}$ decreases monotonically with rarefaction, for the sphere $\\Lambda_{PP}$ can increase, peak, or decrease depending on the molar fraction. The cross coefficient is so sensitive that in some conditions the hard-sphere model predicts a sign opposite to the ab initio potential at the same state point, and the inverted temperature gradient appears near the sphere, with its presence controlled by the helium fraction.","pith_inferences":["A sign change in $\\Lambda_{PT}$ implies a regime where a temperature difference alone drives argon mass flow in the direction opposite to the saturated-pressure gradient; this could be searched for by levitating an argon particle and imposing a small thermal asymmetry while sweeping the background pressure.","If the real surface condensation or sticking coefficient is not exactly unity, the published tables would shift quantitatively, but the qualitative structure (sign change of $\\Lambda_{PT}$ and non-monotonic $\\Lambda_{PP}$) likely persists, so the free-molecular reference values could be rescaled by an effective sticking fraction as a first-order correction.","The same Onsager-matrix methodology extends naturally to dust grains, aerosol droplets, and sublimating particles in vacuum systems and planetary atmospheres, where the sphere is the natural geometry; the strong potential sensitivity warns that realistic collision cross-sections are needed in those applications too."],"forward_implications":["Any small driving force (pressure or temperature difference) at an argon sphere in helium now yields the mass and energy flow rates directly from the relations $\\dot M=4\\pi R_0^2 n_0 v_0 m_2 J_P$ and $\\dot E=4\\pi R_0^2 v_0 p_0\\left(J_T+\\tfrac52 J_P\\right)$ using the tabulated coefficients.","In the transitional and continuum regimes the hard-sphere model is not reliable for the cross effect: at $C_0=0.5$ and $\\delta=10$ the ab initio and hard-sphere potentials give $\\Lambda_{PT}$ values that differ by more than 100% and even have opposite signs.","Adding helium as a background gas suppresses the mass flow coefficient $\\Lambda_{PP}$ strongly, by about 60% when the helium fraction grows from 0.1 to 0.5 at $\\delta=1$, so even a modest amount of non-condensable gas slows sublimation.","The inverted temperature gradient near the sphere is confirmed and shown to depend on the helium molar fraction and the rarefaction parameter, not on the interatomic interaction potential.","The helium-krypton comparison shows that the mass ratio of the sublimating species strongly shifts the heat-flow coefficient $\\Lambda_{TT}$, with deviations up to about 43%, so results cannot be freely transferred from one sublimating species to another."],"supporting_citations":[{"why":"Supplies the McCormack collision model for binary gas mixtures that replaces the Boltzmann collision integral in all numerical solutions.","marker":"[27]"},{"why":"Gives the planar-geometry sublimation and deposition results that the sphere results are qualitatively compared against.","marker":"[10]"},{"why":"Provides the ab initio Omega-integrals for the helium-argon mixture that define the AI-potential collision data.","marker":"[48]"},{"why":"Supplies the transport coefficients of argon and its helium mixtures based on ab initio potentials, used for the viscosity in the mean free path and for the dimensional flow rates.","marker":"[43]"},{"why":"Provides the sublimation curves of argon and krypton that fix the two reference temperatures and the saturation pressures defining the thermodynamic forces.","marker":"[28]"},{"why":"Establishes the Onsager-Casimir reciprocal relations from the Boltzmann equation for gas mixtures, justifying the flux definitions and the symmetry of the coefficient matrix.","marker":"[37]"},{"why":"Gives the heat and mass transfer from a spherical particle in a rarefied gas that serves as the reference for the single-vapor behavior of $\\Lambda_{PP}$.","marker":"[42]"},{"why":"Provides the split-of-solution method that handles the distribution-function discontinuity at the sphere surface in the numerical scheme.","marker":"[40]"}],"fun_headline_variants":["Argon sphere sublimation cross term flips with helium fraction","Cross coefficient sign sensitive to rarefaction and potential","Sublimation at sphere: kinetic coefficients vary with gas mix","Helium background tunes sublimation cross term at argon sphere"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The surface is assumed to be a perfect phase-change interface: every incident argon or krypton atom is absorbed and the surface re-emits vapor atoms with a Maxwellian distribution at the surface temperature and saturation pressure, while helium is fully accommodated; if the real condensation or sticking coefficient is less than one, the tabulated coefficients change.","fun_headline_variants_meta":{"raw":{"variants":["Argon sphere sublimation cross term flips with helium fraction","Cross coefficient sign sensitive to rarefaction and potential","Sublimation at sphere: kinetic coefficients vary with gas mix","Helium background tunes sublimation cross term at argon sphere"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000674,"raw_usage":{"total_tokens":3065,"prompt_tokens":937,"completion_tokens":2128,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":553,"completion_tokens_details":{"reasoning_tokens":2060}},"tokens_in":553,"tokens_out":2128,"duration_ms":18742,"temperature":1.0,"reasoning_tokens":2060,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T17:24:53.912764+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the steady mass flow from a small argon sphere suspended in a helium-argon mixture at $T_0=50$ K with $C_0=0.5$ while sweeping the pressure so that the rarefaction parameter runs from about 0.1 to 10, and applying a small temperature difference of about 0.1 K. The predicted sign change of $\\Lambda_{PT}$ would appear as a reversal of the net argon flow direction relative to the temperature gradient, and the particular value of $\\delta$ at which the reversal occurs would distinguish the hard-sphere from the ab initio potential prediction.","supporting_citations":[{"cited_title":"Construction of linearized kinetic mode ls for gaseous mixture and molecular gases","cited_arxiv_id":null,"evidence_quote":"Supplies the McCormack collision model for binary gas mixtures that replaces the Boltzmann collision integral in all numerical solutions."},{"cited_title":"Sublimation and d eposition in gaseous mixtures","cited_arxiv_id":null,"evidence_quote":"Gives the planar-geometry sublimation and deposition results that the sphere results are qualitatively compared against."},{"cited_title":"Transport coefficients of heli um-argon mixture based on ab initio potential","cited_arxiv_id":null,"evidence_quote":"Provides the ab initio Omega-integrals for the helium-argon mixture that define the AI-potential collision data."},{"cited_title":"Transport coefficients of argo n and its mixtures with helium and neon at low density based ab initio potentials","cited_arxiv_id":null,"evidence_quote":"Supplies the transport coefficients of argon and its helium mixtures based on ab initio potentials, used for the viscosity in the mean free path and for the dimensional flow rates."},{"cited_title":"The sublimation of argon, kryp ton, and xenon","cited_arxiv_id":null,"evidence_quote":"Provides the sublimation curves of argon and krypton that fix the two reference temperatures and the saturation pressures defining the thermodynamic forces."},{"cited_title":"Onsager-Casimir reciprocal r elations based on the Boltz- mann equation and gas-surface interaction","cited_arxiv_id":null,"evidence_quote":"Establishes the Onsager-Casimir reciprocal relations from the Boltzmann equation for gas mixtures, justifying the flux definitions and the symmetry of the coefficient matrix."},{"cited_title":"The kinetic theory of heat and mass transfer from a spherical particle in a rarefied gas","cited_arxiv_id":null,"evidence_quote":"Gives the heat and mass transfer from a spherical particle in a rarefied gas that serves as the reference for the single-vapor behavior of $\\Lambda_{PP}$."},{"cited_title":"The driven cavity flow over th e whole range of the Knudsen number","cited_arxiv_id":null,"evidence_quote":"Provides the split-of-solution method that handles the distribution-function discontinuity at the sphere surface in the numerical scheme."}],"review_version":1}