{"id":"5cd36615-c4f8-446f-8b05-efbdf3fd7da9","arxiv_id":"2607.09574","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":2,"one_line_summary":"Exponential models in Mach number, with coefficients as low-order polynomials in the specific-heat ratio, accurately approximate maximal-entropy-production evaporation curves for calorically perfect polyatomic gases.","lead":"Simple exponential formulas accurately describe pressure and temperature ratios versus Mach number for steady evaporation of polyatomic gases under maximal total entropy production. The same forms also match independent kinetic simulations and generalize cleanly via the specific-heat ratio.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified","rationale":"The Reader correctly isolates the maximal-entropy-production identification as the sole non-trivial assumption and correctly judges that the reported fitting accuracy stands independently of that assumption. Because the same elementary functional form also matches the independent kinetic data to within a few thousandths, the practical utility of the proposed models is already partially corroborated. No sharper technical flaw (e.g., in the derivation of the modified bound, the choice of polynomial degrees, or the least-squares procedure) is visible that would justify altering the ACCEPT verdict or lowering confidence.","tokens_in":16412,"tokens_out":485,"duration_ms":25298,"concrete_test":"Extract the numerical (pδ#(M),Tδ#(M)) pairs underlying Tables 1–2 (or re-optimize eΛ on a denser M-grid for δ=0 and δ=2 at M=0.5 and M=0.9) and compare them pointwise with the corresponding DSMC or moment-method states already tabulated by Morozov et al. and Frezzotti; if any absolute discrepancy exceeds 0.01 the claimed relevance of the maximizers would be weakened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper’s central accuracy claims (MAPE <0.01 % for T and ~0.10 % for p on the maximizers of eΛ, plus absolute differences <0.004 when the same functional form is refitted to independent DSMC/BGK/MM/Holway data) are directly supported by the tabulated coefficients, R^{2}/MAPE values, and Figures 2–3 and 6. The only substantive premise that could undermine practical utility as fluid-dynamic boundary conditions is the inherited maximal-entropy-production hypothesis that the maximizers of the explicit upper-bound functional eΛ(p,T,M) (Section 4) lie close to the true far-field states of the Boltzmann half-space problem. That premise is already flagged by the Reader; the present work mitigates it empirically by showing that the identical exponential/γ-polynomial ansatz also reproduces the independent kinetic-simulation curves to high accuracy. No internal inconsistency, algebraic error, or unsupported numerical claim appears in the fitting procedure itself.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper revisits the half-space evaporation problem for a calorically perfect polyatomic gas (rigid-rotor density of states) under complete absorption, using a modified total-entropy-production functional eD(f) = D(f)/(n∞u). It derives an explicit upper bound eΛ(p,T,M) and numerically maximizes it over a dense Mach-number grid to obtain the curves (pδ#(M), Tδ#(M)). These curves are then fitted by simple exponential models, first for fixed δ and then by a γ-dependent generalization whose coefficients are low-order polynomials in the ratio of specific heats. The same functional form, after re-fitting coefficients, also reproduces independent DSMC, BGK, S-model, Holway and moment-method evaporation data to high accuracy. Mixed diffuse-reflection/absorption boundary conditions are treated briefly, and the necessary conditions of Proposition 1 are extended accordingly.","tokens_in":16647,"tokens_out":716,"duration_ms":10901,"significance":"If the maximizers of the entropy bound remain close to the true far-field states of the Boltzmann half-space problem, the resulting closed-form expressions supply compact, ready-to-use fluid-dynamic boundary conditions for polyatomic evaporation that cover a continuous range of γ. The empirical observation that the identical ansatz also fits independent kinetic-simulation data sets (absolute differences <0.004) strengthens the practical utility of the formulas even if the maximal-entropy-production hypothesis is only approximate. The work therefore offers both a quantitative characterization of the entropy-production surface and immediately usable analytic models for continuum simulations.","major_comments":[],"minor_comments":[{"comment":"The abstract and introduction repeatedly use the phrase “slightly modified entropy functional” without an immediate pointer to the precise definition (11) and the relation eD = D/(n∞u). A single clarifying sentence early in Section 4 would help the reader.","section":null},{"comment":"Figure 1 caption and the surrounding text refer to “required corrections for a perfect fit,” yet the vertical scales of the two panels differ by an order of magnitude; a common scale or an explicit statement of the different ranges would improve readability.","section":null},{"comment":"In Table 3 the column header “n” is never defined in the caption; a parenthetical remark that n is the number of data points used for each fit would remove any ambiguity.","section":null},{"comment":"Equation (21) presents the optimized coefficient functions without stating the precise least-squares residual or the training set size; a short sentence quantifying the residual would make the claim of “excellent accuracy” fully self-contained.","section":null},{"comment":"A few typographical inconsistencies appear (e.g., “Alsoanextension” in the abstract, missing spaces after commas in several places). A careful copy-edit pass would eliminate them.","section":null}],"recommendation":"accept","confidential_remarks":"The manuscript is a natural and well-executed sequel to the authors’ earlier Nonlinearity paper. Its main novelty is the systematic construction of compact analytic models rather than a new existence theory; that is appropriate for the journal provided the fitting results are presented with the same care as the analytic bounds. I see no citation or novelty-disclosure issues."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The punchline is that Bernhoff and Wadbro give compact, ready-to-use exponential formulas for the pressure and temperature ratios versus Mach number that come from maximizing a slightly rescaled total-entropy-production functional, then show the same functional form (with refitted coefficients) also reproduces independent DSMC, BGK, Holway and moment-method evaporation curves to absolute differences under 0.004. They further package the coefficients as low-order polynomials in γ so one expression covers a range of internal degrees of freedom.\n\nWhat is new is the modified functional eD = D/(n∞u), the systematic least-squares fits (basic exponential, then improved with M^{2} and M^{3} terms, then the four-polynomial γ-model with explicit coefficients in (21)), the MAPE tables (sub-0.01 % for T, ~0.1 % for p on their maximizers), and the mixed-boundary extension with σe. The math is standard and transparent: the upper bound eΛ is written out, maximization is done on a dense Mach grid, and the nonlinear fits report near-unit R^{2}. Cross-checks against Frezzotti and Morozov et al. data are independent and look solid. Citation pattern is appropriate; self-citation to their earlier bound is necessary rather than decorative.\n\nThe soft spot is the inherited premise that the maximizers of eΛ sit close enough to the true far-field states of the Boltzmann half-space problem to serve as fluid-dynamic boundary conditions. That is already flagged in the literature they build on; the present paper mitigates it empirically by showing the same ansatz fits the kinetic-simulation curves well, but it does not prove the maximizers are the actual solutions. No public code or raw tables is a minor practical inconvenience, not a soundness issue. Nothing circular or algebraically broken appears.\n\nThis is for people who need simple analytic interface conditions for continuum or hybrid simulations of polyatomic evaporation/condensation. It is not a foundational theoretical advance, but it is careful, usable, and well-executed. I would send it to peer review without hesitation; a serious referee can check the fitting details and the physical relevance of the maximizers. Worth engaging if you work in the area.","headline":"Clean, usable exponential and γ-polynomial fits to polyatomic evaporation curves that also match independent kinetic data; solid incremental work.","tokens_in":17266,"tokens_out":549,"would_cite":true,"duration_ms":7813,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76P05","82C40","35Q20"],"pacs":[],"model":"grok-4.5","headline":"Simple exponential functions of Mach number and specific-heat ratio accurately reproduce maximal-entropy-production evaporation curves for polyatomic gases.","keywords":["kinetic theory","Boltzmann equation","evaporation","polyatomic gas","entropy production","half-space problem","Mach number","ratio of specific heats"],"falsifier":"Compute the exact asymptotic pressure and temperature ratios of the Boltzmann half-space problem for a few polyatomic gases (by DSMC or a high-resolution kinetic solver) and check whether they deviate systematically from the exponential maximizers by more than a few thousandths over the Mach-number interval (0,1].","tokens_in":17284,"feed_emoji":"☁️","tokens_out":993,"duration_ms":10939,"temperature":0.7,"pith_summary":"The paper studies the half-space problem of steady evaporation of a calorically perfect polyatomic gas whose molecules behave as rigid rotors. It maximizes a modified total-entropy-production functional over the relative pressure and temperature that appear as far-field parameters, then shows that the resulting curves of pressure and temperature versus Mach number are extremely well described by compact exponential expressions. These expressions, once their four coefficients are allowed to depend on the ratio of specific heats through low-order polynomials, remain accurate across a continuous range of internal degrees of freedom. The same functional form also fits independent numerical evaporation data obtained by other kinetic methods to within absolute differences of a few thousandths. The resulting closed-form models supply ready-to-use boundary conditions for the compressible Euler equations at an evaporating interface.","feed_headline":"Four numbers fit polyatomic evaporation curves to 0.1 percent","feed_subtitle":"Exponential models of Mach number and heat-capacity ratio replace costly kinetic simulations of the Knudsen layer.","key_machinery":"The modified total-entropy-production functional eD(f) and its explicit upper bound eΛ(p,T,M). Maximizing eΛ with respect to the relative pressure and temperature for each Mach number produces the target curves that the exponential models are fitted to.","core_discovery":"The pairs (pressure ratio, temperature ratio) that maximize the modified total-entropy-production upper bound for each fixed Mach number and number of internal degrees of freedom lie on curves that are captured, to MAPE below 0.01 percent for temperature and about 0.10 percent for pressure, by the two-parameter exponential ansätze p = exp(-β₁ M + β₂ M²) and T = exp(-α₁ M - α₂ M³). The four coefficients themselves admit smooth polynomial representations in the specific-heat ratio γ, yielding a single four-function model valid for the whole family of calorically perfect gases considered.","pith_inferences":["Because the temperature maximizer is independent of the accommodation coefficient while the pressure maximizer is not, mixed-boundary simulations can be reduced to a one-parameter family of pressure curves once temperature is fixed by the pure-absorption formula.","The low-dimensional structure (only four smooth coefficient functions of γ) suggests that a similar exponential reduction may exist for the condensation side of the phase-transition surface.","If the maximizers prove close to true Boltzmann states, the analytic models could be inserted directly into continuum CFD codes as interface conditions without any kinetic pre-computation."],"forward_implications":["Closed-form exponential boundary conditions for the Euler equations become available for any calorically perfect polyatomic gas once its specific-heat ratio is known.","The same functional shape can be reused, after a trivial rescaling of coefficients, for mixed diffuse-reflection/absorption interfaces.","Numerical evaporation data obtained by DSMC, BGK, S-model or moment methods can be replaced by the analytic fits with absolute error less than 0.004.","The continuous dependence on γ permits direct interpolation to any intermediate number of internal degrees of freedom without new kinetic simulations."],"fun_headline_variants":["Four coefficients fit polyatomic evaporation curves to 0.1%","Two exponentials capture pressure and temperature ratios vs Mach","γ-polynomials turn four numbers into full evaporation-curve model","Maximal entropy production yields exponential fits for any γ","Simple M-exponentials replace Knudsen-layer simulations"],"cache_read_input_tokens":128,"weakest_assumption_plain":"That the pressure and temperature values which maximize the entropy-production upper bound are the same values that actually appear in a true solution of the Boltzmann half-space problem.","fun_headline_variants_meta":{"raw":{"variants":["Four coefficients fit polyatomic evaporation curves to 0.1%","Two exponentials capture pressure and temperature ratios vs Mach","γ-polynomials turn four numbers into full evaporation-curve model","Maximal entropy production yields exponential fits for any γ","Simple M-exponentials replace Knudsen-layer simulations"]},"model":"grok-4.5","effort":"low","cost_usd":0.004966,"raw_usage":{"total_tokens":1427,"prompt_tokens":803,"num_sources_used":0,"completion_tokens":86,"cost_in_usd_ticks":49660000,"prompt_tokens_details":{"text_tokens":803,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":538,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":803,"tokens_out":86,"duration_ms":8944,"temperature":1.0,"reasoning_tokens":538,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-13T02:02:06.317988+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Compute the exact asymptotic pressure and temperature ratios of the Boltzmann half-space problem for a few polyatomic gases (by DSMC or a high-resolution kinetic solver) and check whether they deviate systematically from the exponential maximizers by more than a few thousandths over the Mach-number interval (0,1].","supporting_citations":[],"review_version":1}