{"id":"ec67fdb9-f8fd-4aa1-a534-333cdb39f463","arxiv_id":"1908.05350","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Random sequential adsorption of spheres on a cylinder yields a coverage that depends on the sphere-to-cylinder radius ratio, lies below planar RSA coverage, and shows weak short-range chiral ordering.","lead":"This paper simulates and measures how randomly packed spheres cover a cylindrical wire, as a function of the ratio of sphere radius to wire radius. It finds that wire curvature lowers the maximum random coverage below the flat-surface value and reveals a weak chiral ordering at short scales, which is relevant to colloidal coatings and biological structures like dental plaque.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper's high-curvature ρ(2D)∞ (Fig. 4) depends on extrapolating finite-time simulations with df=2, yet no quantitative evidence is shown that the kinetics exponent does not cross over to 1D for small r̃.","rationale":"I read the full text in good faith and verified the central mapping. The reduction of 3D sphere adsorption to a 2D RSA of the angular-envelope shapes in Eq. (1) is exact: for two spheres whose centers lie on the cylinder of radius R+r, their radial projections onto the φ–z strip overlap if and only if the 3D center distance is below 2R. This was initially a possible concern but is not the weak point. The genuine soft spot is the extrapolation of N∞ using the df=2 kinetic law. The reader's weakest assumption identified exactly this issue, and I agree. The paper's assertion that no crossover to 1D kinetics occurs is unsupported by any quantitative analysis, and the non-monotonic ρ(2D)∞ is presented without error bars or alternative estimation. Because the high-curvature regime is inaccessible experimentally, the simulation extrapolation is the sole basis for the non-monotonicity. This does not overturn the experimental validation at weak-to-moderate curvature, but it does mean the high-curvature prediction should be treated as provisional. The reader's CONDITIONAL verdict already captures this, so my independent pass does not move the verdict. The concrete test I propose would settle the concern by measuring the acceptance probability, which yields the jamming density without assuming a particular kinetic exponent, and by checking the local exponent directly.","tokens_in":9551,"tokens_out":30241,"duration_ms":291839,"concrete_test":"For r̃ = 0.05, 0.02, and 0.01, run the effective 2D RSA simulation to very long times (e.g., 10^9 attempts) and measure the acceptance probability φ(ρ) = dρ/dτ as a function of the instantaneous coverage ρ. Extrapolate φ(ρ) to zero to obtain an exponent-free estimate of ρ∞, and compare with the τ^{-1/2}-extrapolated ρ(2D)∞ from Fig. 4. In addition, compute the local logarithmic slope d ln(ρ∞_est−ρ(τ))/d ln τ in successive time decades; if the slope drifts from -1/2 toward -1 as τ increases for small r̃, the df=2 assumption fails and the reported non-monotonicity is an extrapolation artifact. A simpler internal check: re-extrapolate N∞ using only the last part of the time series and compare with the full-series extrapolation; systematic shifts with the fitting window indicate that the assumed power law is not valid.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claims at high curvature rest on the assumption that the asymptotic RSA kinetics obeys ρ∞−ρ(τ)∼τ^{-1/2} (df=2) for every r̃, as stated in the kinetics paragraph and Eq. (4). The non-monotonic ρ(2D)∞ reported for r̃≲0.2 is obtained by extrapolating N∞ from this law. As the authors note, the effective 2D strip becomes increasingly narrow in the φ direction as r̃→0, so a crossover to quasi-1D kinetics could change the exponent to df=1. The only defense offered is the sentence 'we do not observe a crossover to 1D asymptotic kinetics in our simulations as r̃→0', with no exponent measurement, no time-window analysis, and no supporting figure. If the true exponent drifts toward -1 at small r̃, fitting with τ^{-1/2} systematically underestimates the asymptotic density: a τ^{-1} decay fitted with a τ^{-1/2} model pulls the intercept below the true asymptote. Thus the values of ρ(2D)∞ at small r̃, and hence the claimed non-monotonicity, could be an artifact of the extrapolation. The same assumption underlies the reported monotonic increase of ρ(3D)∞ at small r̃. Because the high-curvature regime is not experimentally accessible in this work, the entire prediction in this regime hinges on this unverified kinetic assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies random sequential adsorption (RSA) of hard spheres on a long cylinder of radius r, with sphere radius R, parameterized by r̃ = r/R. It maps the process to RSA of two-dimensional shapes on a strip in the φ–z plane, derives coverage formulas ρ(2D) and ρ(3D), and obtains asymptotic coverage from simulations extrapolated using the τ^{-1/2} law. Experiments with electrostatic and DNA-mediated binding of colloidal spheres to silica wires validate the predicted linear density and 3D coverage for r̃ ≳ 0.2. For high curvature (r̃ ≲ 0.2), simulations predict non-monotonic ρ(2D)∞ and monotonic ρ(3D)∞. The paper also reports weak chiral ordering, with tan θmax ≈ 3 r̃.","tokens_in":9802,"tokens_out":5635,"duration_ms":53294,"significance":"If the results hold, the paper provides a parameter-free prediction of curvature-dependent RSA coverage, validated by two independent experimental chemistries, and identifies a surprising non-monotonic 2D coverage. The geometric mapping is elegant, the simulation methodology is standard, and the experimental protocols are detailed, including error bars that account for particle-counting uncertainty. However, the high-curvature predictions rest on an unverified kinetic extrapolation, and the chiral ordering law is reported without statistical uncertainty. The experimental validation covers only the weak-to-moderate curvature regime, so the high-curvature claims remain purely predictive.","major_comments":[{"comment":"The asymptotic values ρ(2D)∞ and ρ(3D)∞ at r̃ ≲ 0.2 are obtained by extrapolating finite-time simulation data with the law ρ∞−ρ(τ) ∼ τ^{-1/2}, which assumes df = 2 for all r̃. The paper's only defense is the sentence 'we do not observe a crossover to 1D asymptotic kinetics in our simulations as r̃ → 0,' but no exponent measurement, no local-slope analysis, and no supporting figure is provided. If the effective dimension crosses to df = 1 at small r̃, the τ^{-1/2} fit systematically biases the intercept below the true asymptote, and the claimed non-monotonicity of ρ(2D)∞ could be an artifact. Please provide a quantitative test, for example a plot of the logarithmic local slope d ln(ρ∞−ρ)/d ln τ versus 1/τ for the smallest r̃ values, or a comparison of extrapolations using df = 1 and df = 2.","section":"Results, kinetics paragraph and Eq. (4)"},{"comment":"The relation tan θmax ≈ 3 r̃ is presented as a quantitative result, but the paper does not report the number of simulation runs, the standard error of θmax, the range of r̃ over which the fit is performed, or the goodness of fit. Without these, the claimed linear law is not falsifiable. Please provide fit details and uncertainty, and state how θmax depends on the correlation cutoff (the paper notes peaks disappear when the range is increased).","section":"Discussion and Fig. 5C"}],"minor_comments":[{"comment":"The citation to packing spheres inside cylinders contains an unresolved '[?]' placeholder (the text reads '[?,35,41]'). Please complete the reference.","section":"Results, packing density paragraph"},{"comment":"The variables z and φ are introduced, but it would help to define the domain of φ and the periodicity explicitly (0 ≤ φ < 2π) before the strip of width 2π is used.","section":"Results, Eq. (1)"},{"comment":"There is a typo in the subsection heading: 'data anaysis' should be 'data analysis'.","section":"Methods, 'data anaysis' heading"},{"comment":"The caption says '⟨n(r)n(r',θ)⟩ as a function of angle θ' but the text defines it as counting pairs aligned at angle θ within tolerance Δy; the precise normalization of the correlation function (e.g., whether it is divided by the strip area) is not given.","section":"Results, Fig. 5B"},{"comment":"Several experimental points at small r̃ fall outside the two-sigma band, and the authors attribute this to weak binding on highly curved wires. This is plausible, but the statement 'these results validate our understanding' is too strong for the r̃ ≲ 0.2 region; the subsequent paragraph already restricts the validation to r̃ ≳ 0.2, so consider moving the validation phrase to that regime or adding a quantitative model of the undersaturation effect.","section":"Results, Fig. 3C and 3D"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope and the central geometric mapping is sound. The main concern is the unverified kinetic exponent in the high-curvature regime, which is load-bearing for the non-monotonic 2D coverage claim. The unresolved reference placeholders suggest the manuscript may not yet be in final form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What you should know: this paper gives the first systematic map of sphere random sequential adsorption on a cylinder over the whole radius ratio r̃, and the experimentally tested part is convincing. The new stuff is the full r̃ dependence of the asymptotic coverage, the non-monotonic ρ(2D)∞, the chiral ordering observation, and validation with two independent binding chemistries (electrostatic and DNA). The geometric mapping to a 2D strip is sound, the simulations are standard RSA, and the agreement at r̃ ≳ 0.2 is solid. The claim that cylinder curvature lowers long-time coverage relative to a plane holds up where the experiments actually reach.\n\nThe soft spots are real and concentrated in the high-curvature regime. The asymptotic density at small r̃ is extracted by assuming ρ∞−ρ(τ) ~ τ^{-1/2} (df=2), with no quantitative check of the kinetic exponent. As the unwrapped strip narrows, quasi-1D behavior is a plausible risk, and if the exponent drifts, the extrapolated N∞ and the non-monotonic ρ(2D)∞ could be biased. The authors say they see no crossover, but there is no supporting figure or exponent measurement. That is the load-bearing assumption for all r̃ ≲ 0.2 predictions, and it deserves a hard look before the paper is taken as settled. The chiral law tanθmax≈3r̃ is a fit with no uncertainty intervals, though it is a side result. Minor issues: the arXiv text contains ‘?’ placeholders for a citation and the fiber-pulling method, and there is no code/data deposit. None of these affect the experimentally validated middle range.\n\nOverall, this is a useful paper for the RSA and colloid science community, and the high-curvature part is an honest prediction, not an overclaim. A serious referee should probe the kinetic exponent assumption for r̃→0, but the paper deserves refereeing rather than desk rejection.","headline":"Solid RSA-on-cylinder paper with believable experiments and a clean mapping, but the high-curvature predictions rest on an unchecked kinetic exponent and should be treated as predictions, not settled results.","tokens_in":10361,"tokens_out":3700,"would_cite":true,"duration_ms":39924,"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":"The paper pinpoints how a cylinder's radius ratio governs the jamming coverage of spheres, which stays below the flat-plane value and makes the 2D density peak at intermediate curvature.","keywords":["random sequential adsorption","cylinder","wire curvature","jamming coverage","colloidal spheres","chiral ordering","non-monotonic coverage","irreversible adsorption"],"falsifier":"A direct simulation or experiment at very high curvature (r̃ ≲ 0.05) that measures the approach to jamming over many decades of time: if the exponent in ρ∞ − ρ(τ) ~ $τ^{{−1/d_f}}$ drifts away from d_f = 2 toward d_f = 1 as the wire thins, the extrapolated N∞ and hence the reported coverage curves are not the true asymptotic values.","tokens_in":9327,"feed_emoji":"🧪","tokens_out":6504,"duration_ms":59070,"temperature":0.7,"pith_summary":"The paper asks how tightly equal spheres can randomly pack onto a long cylindrical wire when they stick irreversibly where they first touch, the process known as random sequential adsorption. By unrolling the cylinder into a strip and treating each sphere as an oblong shape in the angle–axis plane, the authors reduce the problem to a 2D parking problem controlled by one number, the ratio r̃ = r/R. They find that the long-time (jamming) coverage in 3D increases monotonically with r̃ and always lies below the familiar flat-plane value, while the 2D area coverage is non-monotonic, peaking at an intermediate curvature. Colloidal experiments on silica wires confirm the predicted linear particle density and the below-plane coverage for weak to moderate curvature. The result matters because it gives a parameter-free prediction for a basic deposition geometry seen in dental plaque, sea-grapes, and coated wires.","feed_headline":"Wires pack adsorbed spheres looser than flat surfaces do","feed_subtitle":"Radius ratio alone sets how tightly spheres jam on a wire, and the 2D density peaks mid-curvature.","key_machinery":"The load-bearing object is the exact mapping of the 3D adsorption problem onto a 2D random sequential adsorption process in the φ–z plane. A sphere centered at (z0, φ0) becomes a lens-shaped exclusion zone with boundary φ(z) = φ0 ± arcsin[√(R² − (z − z0)²)/(R + r)], and the cylinder becomes a periodic strip of width 2π. Simulations deposit these shapes randomly on the strip, and the asymptotic particle number is extracted from the universal kinetics ρ∞ − ρ(τ) ~ $τ^{{−1/2}}$ for two-degree-of-freedom parking. The conversion between the 2D and 3D coverages is carried by the one-parameter formulas ρ(2D)(r̃) = λ(r̃)Ã(r̃)/(2π) and ρ(3D)(r̃) = λ(r̃)/(3(1 + r̃)).","core_discovery":"The central discovery is that random sequential adsorption of spheres on a cylinder is not a small perturbation of flat-surface parking: the asymptotic coverage depends on the ratio of sphere radius to cylinder radius, with ρ(3D)_∞ decreasing monotonically as the wire gets thinner and approaching ≈0.3647 on a plane, while the unwrapped 2D coverage ρ(2D)_∞ rises from the thin-wire limit, overshoots the planar value, and then falls back to ≈0.5471. The non-monotonicity follows from the product of two competing factors: the per-sphere angular area shrinks with r̃, while the saturated linear density grows because the maximum angular extent of a sphere decreases. Simulations show ρ(3D)_∞ increasing monotonically with r̃, and experiments with electrostatically and DNA-bound colloidal spheres match the predicted linear density over the accessible range r̃ ≳ 0.2. At short scales the adsorbed structures show weak chiral ordering with the preferred alignment angle obeying tan θ_max ≈ 3r̃.","pith_inferences":["The τ^{−1/2} extrapolation assumes two degrees of freedom down to arbitrarily thin wires; if quasi-1D kinetics emerge for r̃ → 0, the extrapolated N∞ and the reported non-monotonic peak of ρ(2D)_∞ could shift, so the high-curvature predictions remain the least secure part of the paper.","The experimental undersaturation at small r̃ hints that binding-energy decrease with curvature sets a practical irreversibility limit; using stronger bonds (covalent or longer DNA) could push experiments into the high-curvature regime and test the simulated coverage directly.","The method transfers to other convex particles (ellipsoids, rods) on a cylinder: one only needs the angular envelope shape, so the same equation (1) input could predict coverage families for different particle shapes.","The linear pitch-angle law tan θ_max ≈ 3r̃ offers a cheap experimental route to verifying the chiral ordering: count neighbor-center angles in 2D projection images of adsorbed spheres, no 3D reconstruction needed."],"forward_implications":["The 3D jamming coverage ρ(3D)_∞ increases monotonically with r̃ and never reaches the flat-plane limit of ≈0.3647, so any cylindrical curvature loosens the adsorbed monolayer in 3D.","The 2D coverage ρ(2D)_∞ is non-monotonic: it overshoots the planar disk value ≈0.5471 at intermediate r̃, a distinctive signature of the competition between per-sphere area and angular packing.","For very thin wires the random-parking density sits at about 61% of the densest possible packing, a ratio that stays near 0.62 across r̃ values.","Short-range order in the jammed state is chiral, with the preferred pitch angle obeying tan θ_max ≈ 3r̃, so wire curvature biases the local alignment of neighboring spheres.","The mapping to 2D RSA with shape (1) provides a parameter-free recipe to compute the coverage for any wire–sphere radius ratio without simulating the full 3D deposition."],"supporting_citations":[{"why":"Renyi's exact solution of 1D car parking, the benchmark for the one-dimensional limit of the problem.","marker":"[12]"},{"why":"Supplies the accepted asymptotic coverage of disks on a plane, ρ ≈ 0.5471, the flat-surface value the cylinder results are compared against.","marker":"[37]"},{"why":"Provides the universal kinetic law ρ∞ − ρ(τ) ~ τ^{−1/d_f} used to extrapolate the asymptotic particle number from simulations.","marker":"[40]"},{"why":"Establishes the τ^{−1/2} extrapolation method for two-degree-of-freedom RSA that the simulations rely on to obtain N∞.","marker":"[20]"},{"why":"Gives the earlier unrolling approximation ρ^cyl_∞ = ρ^plane_∞(1 + R/r) that this paper generalizes beyond the weak-curvature regime.","marker":"[5]"},{"why":"Applies the same unrolling approximation to sphere-on-cylinder adsorption, serving as the baseline prediction that fails as R/r grows.","marker":"[33]"},{"why":"Provides close-packed densities for spheres inside cylinders, used to compare random and close-packed coverage ratios.","marker":"[35]"}],"fun_headline_variants":["Cylinder curvature reshapes sphere parking limits","Adsorbed spheres pack densest mid-curvature wire","Wire radius alone sets how spheres jam and pack","Sphere coverage peaks on curved surfaces, not flat"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The asymptotic particle number is obtained by extrapolating finite-time simulations using the law ρ∞ − ρ(τ) ~ $τ^{{−1/2}}$, which assumes the effective number of degrees of freedom stays at two for every wire curvature; if the kinetics slow to one-dimensional behavior for very thin wires, the extrapolated coverages would be biased.","fun_headline_variants_meta":{"raw":{"variants":["Cylinder curvature reshapes sphere parking limits","Adsorbed spheres pack densest mid-curvature wire","Wire radius alone sets how spheres jam and pack","Sphere coverage peaks on curved surfaces, not flat"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000188,"raw_usage":{"total_tokens":1287,"prompt_tokens":853,"completion_tokens":434,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":469,"completion_tokens_details":{"reasoning_tokens":372}},"tokens_in":469,"tokens_out":434,"duration_ms":5431,"temperature":1.0,"reasoning_tokens":372,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:17:24.925080+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct simulation or experiment at very high curvature (r̃ ≲ 0.05) that measures the approach to jamming over many decades of time: if the exponent in ρ∞ − ρ(τ) ~ $τ^{{−1/d_f}}$ drifts away from d_f = 2 toward d_f = 1 as the wire thins, the extrapolated N∞ and hence the reported coverage curves are not the true asymptotic values.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Renyi's exact solution of 1D car parking, the benchmark for the one-dimensional limit of the problem."},{"cited_title":"and Ziff, R.M., 2018","cited_arxiv_id":null,"evidence_quote":"Supplies the accepted asymptotic coverage of disks on a plane, ρ ≈ 0.5471, the flat-surface value the cylinder results are compared against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the universal kinetic law ρ∞ − ρ(τ) ~ τ^{−1/d_f} used to extrapolate the asymptotic particle number from simulations."},{"cited_title":"J.; Ziff, R","cited_arxiv_id":null,"evidence_quote":"Establishes the τ^{−1/2} extrapolation method for two-degree-of-freedom RSA that the simulations rely on to obtain N∞."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the earlier unrolling approximation ρ^cyl_∞ = ρ^plane_∞(1 + R/r) that this paper generalizes beyond the weak-curvature regime."},{"cited_title":"Langmuir 2014, 30, 692-699","cited_arxiv_id":null,"evidence_quote":"Applies the same unrolling approximation to sphere-on-cylinder adsorption, serving as the baseline prediction that fails as R/r grows."},{"cited_title":"K., Weaire, D., and Hutzler, S","cited_arxiv_id":null,"evidence_quote":"Provides close-packed densities for spheres inside cylinders, used to compare random and close-packed coverage ratios."}],"review_version":1}