{"id":"ac41377a-4c1d-41ed-8d3d-5b7e440bfe80","arxiv_id":"2608.05603","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":11,"one_line_summary":"Argon preplating of MCM-41 nanopores screens the silica surface corrugation by up to about 34% and yields a helium confinement potential whose radial average reproduces the smooth cylindrical model used in earlier quantum Monte Carlo simulations.","lead":"This paper simulates argon-coated MCM-41 nanopores atom by atom and computes the energy landscape a single helium atom would feel inside them. It finds that the argon layer smooths the rough silica surface and shifts helium's preferred position into an annular ring, matching the simplified cylindrical model that earlier quantum simulations assumed.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The load-bearing weakness is that the smooth Eq. (11) potential is fit to an arithmetic average over θ and z (Eq. 10); residual corrugation (Fig.","rationale":"Good-faith reading: the paper is careful, the TPI computation is reasonable, the atomistic model is benchmarked against isotherm and S(q), and the fitted σ, R, and nε are consistent with He-Ar Lennard-Jones parameters and the BET-derived Ar density, which is genuine supporting evidence. However, the step from TPI data to Eq. (11) is an arithmetic average, and the validity of that average for a quantum fluid is not demonstrated. Because the paper explicitly frames Eq. (11) as microscopic support for the QMC potential, the averaging operation is load-bearing. The reader's CONDITIONAL verdict already flags validation weaknesses (four-parameter fit, model fidelity); my concern is a distinct but related threat to the support claim. It does not warrant rejection because it is testable and addressable with the provided data, so I leave the verdict UNCHANGED as CONDITIONAL. Agreement is partial: the reader identified related validation concerns, but not the averaging issue as the primary weakness.","tokens_in":16418,"tokens_out":12407,"duration_ms":145194,"concrete_test":"Perform a single-4He path-integral ground-state (PIGS) or diffusion Monte Carlo calculation in the full frozen TPI potential U(r,θ,z) and in the fitted Eq. (11) potential, using the available code and data [51], and compare the ground-state energy and radial density. Alternatively, compute the exponential (free-energy) radial potential W(r) = -k_B T ln[(1/(Nθ Nz)) Σ_{θ,z} exp(-U(r,θ,z)/k_B T)] at T ≈ 0.5 K from the TPI grid; if |W(r_min) - U(r_min)| exceeds about 1 K, the arithmetic average in Eq. (10) is not a faithful effective potential for low-temperature helium and the microscopic-support claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the mean radial potential U(r) from Eq. (10) validates the smooth cylindrical model used in Ref. [26]. The concern is not the fit itself but the physical relevance of the averaging operation. U(r) is an unweighted arithmetic mean over θ and z of the TPI potential. A helium atom at low temperature does not uniformly sample that mean: its density is exponentially weighted toward the deepest angular and axial channels. The paper reports that corrugation is only partially screened (S(Δr) max ≈34%, Fig. 9), that δU spans hundreds of kelvin in parts of the pore (Fig. 11), and that Fourier modes A_n(r) reach about 25 K near the wall (Fig. 11 insets), implying local deviations of order 10 K at the annular minimum. Replacing the full potential by its arithmetic mean therefore produces a smoother, systematically shallower effective potential than the one a helium atom actually experiences. A fit of Eq. (11) to this mean is thus not by itself microscopic evidence that the QMC potential is faithful; the averaging step assumes the very smoothness the paper claims to establish. The authors acknowledge residual corrugation and provide a GP surrogate, but they do not use it to compute the correct coarse-grained radial potential or to test the QMC approximation against the full atomistic potential. This gap is load-bearing because the support claim is the paper's central message.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reports an atomistic simulation study of Ar-preplated MCM-41. GCMC simulations of Ar adsorption are compared with a measured 90 K isotherm and neutron S(q), MD cooling produces frozen monolayer configurations, and He test-particle insertion yields the spatially resolved one-body confinement potential. The authors find an annular minimum at r≈11.66 Å with depth −44.16 K, fit the radially averaged potential to the continuum cylindrical model of Eq. (11), and quantify residual angular corrugation via Fourier analysis and a Gaussian process surrogate.","tokens_in":16811,"tokens_out":5844,"duration_ms":63232,"significance":"If the central claim holds, the paper provides a valuable bottom-up connection between a specific experimental pore preparation and the one-body potential used in quantum many-body simulations of confined helium, and the open code and data would benefit the field. The work is careful in reporting uncertainties, acknowledges discrepancies in S(q), and does not claim that the potential is featureless. The main limitation is that the validation of the smooth cylindrical potential rests on an arithmetic average of a corrugated energy landscape, so the microscopic support is not yet established at the quantitative level claimed.","major_comments":[{"comment":"The radial potential U(r) is an unweighted arithmetic mean over θ and z. For a helium atom at low temperature, the relevant one-body radial potential is the Boltzmann-weighted average, U_eff(r) = -k_B T ln⟨exp[-U(r,θ,z)/k_B T]⟩_{θ,z}, because the density is exponentially enhanced in the deepest angular and axial channels. The paper's own data show residual corrugation with Fourier amplitudes of order 10–25 K near the wall and δU spanning hundreds of kelvin (Fig. 11), and only partial screening S≈34% near the minimum region (Fig. 9). The arithmetic mean is therefore systematically shallower than the potential actually sampled, and the excellent fit of Eq. (11) to U(r) does not by itself demonstrate that the smooth potential of Ref. [26] faithfully represents helium confinement. I request a quantitative test: compute the Boltzmann-weighted radial potential from the TPI or GP data at a relevant temperature (e.g., 1–4 K) and compare its minimum depth and shape with Eq. (11).","section":"Section IV B, Eq. (10) and Fig. 8"},{"comment":"The comparison with the experimental isotherm uses two reference parameters, µ_ref and P_ref/P0, chosen to align the simulation and experimental inflection points. This calibration removes the main horizontal and vertical offsets, so the 'excellent agreement' in Fig. 2 is partly by construction. Please state explicitly which features of the isotherm (shape, slope, plateau level) are parameter-free predictions, or re-fit with only one global parameter if possible. Without this, the isotherm comparison provides weaker experimental validation of the atomistic model than implied.","section":"Section IV A, Eq. (8) and Fig. 2"},{"comment":"The calculated S(q) reproduces the main peak near 1.8 Å⁻¹ and the secondary peak, but misses fine features just above 2.5 Å⁻¹ and persistent modulations at higher q. Since the TPI potential is computed from the same atomistic model, this discrepancy indicates a possible mismatch in surface structure or surface chemistry that could affect the reported He minimum position and depth. The acknowledgment in the text is helpful, but the central claim needs an assessment of sensitivity, for example by varying silanol density or pore roughness and re-computing U(r).","section":"Section IV A, Fig. 5"},{"comment":"The Gaussian process surrogate is trained only on the z-averaged transverse potential U(r,θ), not on the full U(r,θ,z). The abstract and discussion state that residual corrugation is 'accurately captured by a Gaussian process surrogate,' but the axial corrugation documented in Fig. 10 is not captured, and the text states that a full z-dependent GP is hindered by matrix size. Please either restrict the GP claim to the angular corrugation of the z-averaged potential or provide a full 3D surrogate.","section":"Section III C/D and IV C"},{"comment":"The He test-particle insertion uses a single Lennard-Jones parameter set (σ=2.640 Å, ε/k_B=10.9 K) without specifying how cross interactions with each host species (Ar, O(1), O(2), Si(2), H) are formed. Since the position and depth of the He minimum are central quantitative results, the manuscript must state the mixing rule and the per-species parameters, or point to the exact table in the released code.","section":"Section III C"}],"minor_comments":[{"comment":"There is a typo: 'dependnece' should be 'dependence', and 'Gaussian Process' should be capitalized consistently.","section":"Section IV C"},{"comment":"The normalization of the cross term in S(q) should be stated more carefully: the denominator differs from the first term's normalization, and the convention for the loaded-minus-empty subtraction should be justified.","section":"Eq. (5)"},{"comment":"The color scale for U(r,z) is difficult to read; please label the color axis and clarify whether the values are multiplied by k_B⁻¹ in the same way as the curves.","section":"Fig. 10"},{"comment":"Please state how the parameter uncertainties were obtained (for example, from the covariance of the least-squares fit) and note that U0 is consistent with zero within the reported uncertainty.","section":"Table II"},{"comment":"The phrase 'areal number density of 0.059 Å⁻²' should use explicit units such as atoms/Ų; the current notation is ambiguous.","section":"Section II"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope for cond-mat.mtrl-sci and provides useful atomistic data and open code. The main concern is that the headline validation of the continuum potential is based on an arithmetic average of a strongly corrugated energy landscape; I believe the requested Boltzmann-average test is feasible with existing TPI/GP data and would decide whether the support is quantitative or only qualitative. If the test shows a significant difference, the central claim should be softened to a statement about the mean radial profile rather than the effective potential experienced by helium."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper gives a genuinely new atomistic picture of the He confinement potential in Ar-preplated MCM-41 — annular minimum at 11.66 Å, 44 K depth, screening of the silica corrugation by up to 34%, plus a GP surrogate that reproduces the residual landscape to ~0.8 K. Code and data are available. That is real, useful output for people doing QMC in confined helium.\n\nWhat I’d push back on is the selling line. The claim that the atomistic calculation provides 'microscopic support' for the smooth cylindrical potential used in Ref. [26] rests on a four-parameter fit of Eq. (11) to U(r), where U(r) is the arithmetic average of the TPI potential over θ and z. A helium atom at low temperature does not experience that arithmetic mean; its density is exponentially weighted toward the deepest angular and axial channels. The paper’s own numbers show the corrugation is not tiny: δU spans hundreds of kelvin near the wall, Fourier modes run to ~25 K, and the residual RMS is only partially screened. So fitting the smooth cylinder to the mean demonstrates that the mean is smooth — which the atomistic picture already encoded — but it does not by itself show that a quantum helium drop sees the smooth potential. The authors are honest in the Discussion that the fit 'should not be interpreted as evidence that the microscopic potential is featureless,' but the abstract and introduction still sell it as support. That gap is load-bearing for the paper’s central message.\n\nOther soft spots are minor. The isotherm comparison is calibrated with two reference parameters in Eq. (8), so it is consistency checking, not a sharp test. The neutron S(q) misses fine features above 2.5 Å⁻¹; the authors flag this. The atomistic model comes from Ref. [3], and its fidelity to the real Sigma-Aldrich sample is assumed.\n\nAll that said, this is a careful, well-documented paper and the TPI potential itself is a useful object. The right fix in revision is to compute a coarse-grained radial potential that respects the corrugation — e.g., a Boltzmann-weighted average at the relevant T or the ground-state transverse projection — and then see whether Eq. (11) still describes that. If it does, the microscopic support claim gets teeth. If it doesn’t, the paper still stands as a characterization of the atomistic potential, just with a weaker punchline.\n\nVerdict: send it to peer review. A serious referee will want the coarse-graining addressed, but the raw results deserve publication.","headline":"Useful atomistic TPI potential for helium in Ar-plated MCM-41, but the 'microscopic support' for the smooth continuum model is a four-parameter fit to an arithmetic mean, not a coarse-graining that honors the corrugation.","tokens_in":17362,"tokens_out":3118,"would_cite":true,"duration_ms":32355,"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":"An argon monolayer screens the rough silica wall of MCM-41 and reshapes helium confinement into a smooth cylindrical trap.","keywords":["argon preplating","MCM-41","helium confinement","test-particle insertion","confinement potential","Gaussian process surrogate","neutron scattering","grand-canonical Monte Carlo"],"falsifier":"Run the same helium test-particle insertion on an alternative atomistic MCM-41 model whose silanol density or surface roughness is adjusted to reproduce the fine features near $q \\simeq 2.5$ Å$^{-1}$ and at higher $q$ in the measured static structure factor; if the resulting radial minimum shifts by more than the reported fit uncertainties ($\\sigma_R \\simeq 0.19$ Å and $\\sigma_{U_0} \\simeq 0.65$ K), the claim that the current model fixes the confinement potential is falsified.","tokens_in":16221,"feed_emoji":"🧊","tokens_out":13289,"duration_ms":116854,"temperature":0.7,"pith_summary":"The paper aims to show that preplating MCM-41 with one argon monolayer is a controlled way to engineer the confinement seen by helium: the argon fills the most strongly attractive crevices of the silica surface, screens much of its atomic-scale roughness, and shifts the helium adsorption minimum to an annular region at $r_{\\min}=11.66$ Å with depth $-44.16$ K. It then demonstrates that the radially and axially averaged helium potential is accurately captured by a smooth, continuum cylindrical potential, exactly the effective form used in earlier quantum Monte Carlo simulations of argon-plated MCM-41. If correct, this validates the coarse-grained description of preplated pores and gives an experimentally constrained microscopic potential ready for predictive many-body studies of confined quantum fluids. The result matters because it connects a concrete laboratory preparation, an argon monolayer on a real mesoporous silica, to the one-body term in the Hamiltonian that controls whether confined helium can behave as a quasi-one-dimensional quantum liquid.","feed_headline":"Argon lining turns rough silica pores into smooth helium traps","feed_subtitle":"Helium's minimum moves inward to 11.66 Å with depth −44 K, backing the smooth-potential model.","key_machinery":"The load-bearing object is the helium test-particle insertion potential $U_\\alpha(\\mathbf{r})$ computed on a cylindrical grid ($200\\times144\\times50$ points per frozen configuration) against ten low-temperature Ar-preplated host configurations. Averaging this one-body energy over angle and axial position produces the radial confinement $U(r)$; fitting $U(r)$ to the continuum cylindrical potential of Eq. (11), a helium atom interacting with a perfect cylindrical cavity carved in a continuous medium, is what turns the atomistic model into a smooth effective description. The Gaussian-process surrogate with summed Matérn kernels $K_{5/2}+K_{3/2}$ is the machinery that captures the residual multi-scale corrugation after the smooth radial part is removed.","core_discovery":"On its own terms, the paper's central discovery is that the Ar monolayer does not merely narrow the pore: it transforms the heterogeneous silica adsorption landscape into a nearly cylindrical one. Helium test-particle insertion into frozen, low-temperature Ar-preplated configurations gives a transverse potential whose most attractive part is an annulus at $r_{\\min}\\simeq 11.66$ Å with depth $U(r_{\\min})/k_B = -44.16$ K, compared with the bare-pore minimum at $r_{\\min}\\simeq 14.17$ Å. Averaging over angle and along the pore axis yields a radial potential that fits the continuum cylindrical form of Eq. (11) with parameters $\\sigma = 3.01\\pm0.27$ Å, $R = 14.28\\pm0.19$ Å, $n\\varepsilon/k_B = 0.64\\pm0.21$ K Å$^{-3}$, and $U_0/k_B = 0.14\\pm0.65$ K. The same analysis quantifies the screening: the RMS corrugation of the potential is reduced by up to about 34% near $\\Delta r \\simeq -0.8$ Å. Residual angular and axial corrugation remains, is spread over many Fourier modes, and is reproduced by a Gaussian-process surrogate with a two-Matérn-kernel covariance and a mean absolute error of $0.77\\pm0.03$ K.","pith_inferences":["The same test-particle-insertion-plus-Gaussian-process pipeline could be applied to other preplating agents (neon, molecular hydrogen, cesium) and other probe fluids, producing a one-body confinement potential for each preparation without new experimental input.","The measured maximum screening factor of about 34% suggests an optimization target: preplating coverage or species that raises screening near the adsorption minimum should push the confinement closer to ideal cylindrical geometry.","Because the atomistic model misses fine high-$q$ features in the measured structure factor (for example the peak just above 2.5 Å$^{-1}$ in Fig. 5), the true silanol chemistry or surface roughness could shift the reported potential beyond the fit uncertainties; quantum-simulation predictions using this potential should carry a model-error band."],"forward_implications":["The smooth cylindrical potential used in earlier quantum Monte Carlo work is microscopically justified for argon-preplated MCM-41, not merely assumed.","Helium at low temperature in the preplated pore will first populate an annular shell centered near $r=11.66$ Å rather than adsorbing directly on the silica wall.","The atomistic model reproduces the measured adsorption isotherm and the dominant neutron-scattering peak, anchoring its energetics to experiment.","A quantitative many-body description should retain residual corrugation, which persists over many length scales and can pin defects or renormalize the sound velocity of an emergent one-dimensional quantum liquid."],"supporting_citations":[{"why":"Supplies the experiment-guided atomistic MCM-41 framework and Lennard-Jones parameters that define the pore surface.","marker":"[3]"},{"why":"Provides the experimental sample characterization (hexagonal lattice constant, surface area, pore size) that fixes the pore geometry.","marker":"[16]"},{"why":"Defines the smooth cylindrical effective potential and the earlier quantum Monte Carlo treatment that this paper microscopically validates.","marker":"[26]"},{"why":"Supplies the pore-size-distribution analysis used to set the nominal pore radius of about 3.0 nm.","marker":"[31]"},{"why":"Documents the neutron-scattering instrument used to measure the argon-preplated sample for comparison with simulation.","marker":"[32]"},{"why":"Supplies the grand-canonical Monte Carlo and molecular-dynamics engine used to build and cool the argon monolayer.","marker":"[34]"},{"why":"Supplies the Debye formula used to compute the powder-averaged static structure factor from simulation.","marker":"[39]"},{"why":"Supplies the Gaussian-process regression framework used for the corrugated potential-energy surrogate.","marker":"[40]"},{"why":"Supplies the continuum cylindrical wall-fluid potential of Eq. (11) to which the radial helium potential is fit.","marker":"[49]"}],"fun_headline_variants":["Argon lining gives silica pores a smooth radial helium trap","Ar preplating suppresses silica corrugation, pulls helium inward","Smooth cylindrical potential from Ar-preplated MCM-41 nanopores","Ar monolayer transforms rough pore into smooth quantum-fluid trap"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The atomistic MCM-41 model with its table of Lennard-Jones parameters faithfully represents the real silicon-oxygen surface, including its silanol chemistry and roughness; if the real surface differs, the computed helium minimum, depth, and screening factor shift.","fun_headline_variants_meta":{"raw":{"variants":["Argon lining gives silica pores a smooth radial helium trap","Ar preplating suppresses silica corrugation, pulls helium inward","Smooth cylindrical potential from Ar-preplated MCM-41 nanopores","Ar monolayer transforms rough pore into smooth quantum-fluid trap"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000228,"raw_usage":{"total_tokens":1529,"prompt_tokens":1051,"completion_tokens":478,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":667,"completion_tokens_details":{"reasoning_tokens":406}},"tokens_in":667,"tokens_out":478,"duration_ms":6015,"temperature":1.0,"reasoning_tokens":406,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T05:40:22.764050+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same helium test-particle insertion on an alternative atomistic MCM-41 model whose silanol density or surface roughness is adjusted to reproduce the fine features near $q \\simeq 2.5$ Å$^{-1}$ and at higher $q$ in the measured static structure factor; if the resulting radial minimum shifts by more than the reported fit uncertainties ($\\sigma_R \\simeq 0.19$ Å and $\\sigma_{U_0} \\simeq 0.65$ K), the claim that the current model fixes the confinement potential is falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the experimental sample characterization (hexagonal lattice constant, surface area, pore size) that fixes the pore geometry."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the smooth cylindrical effective potential and the earlier quantum Monte Carlo treatment that this paper microscopically validates."},{"cited_title":"Schlumberger and M","cited_arxiv_id":null,"evidence_quote":"Supplies the pore-size-distribution analysis used to set the nominal pore radius of about 3.0 nm."},{"cited_title":"Thommes, K","cited_arxiv_id":null,"evidence_quote":"Documents the neutron-scattering instrument used to measure the argon-preplated sample for comparison with simulation."},{"cited_title":"Copley and J","cited_arxiv_id":null,"evidence_quote":"Supplies the grand-canonical Monte Carlo and molecular-dynamics engine used to build and cool the argon monolayer."},{"cited_title":"Furukawa, T","cited_arxiv_id":null,"evidence_quote":"Supplies the Debye formula used to compute the powder-averaged static structure factor from simulation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Gaussian-process regression framework used for the corrugated potential-energy surrogate."},{"cited_title":"Akram, S","cited_arxiv_id":null,"evidence_quote":"Supplies the continuum cylindrical wall-fluid potential of Eq. (11) to which the radial helium potential is fit."}],"review_version":1}