{"id":"01cf6306-1004-40eb-934f-78ad02576def","arxiv_id":"2412.15345","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"On a cylinder, the 2D one-component plasma shows boundary density oscillations whose wavelength approaches the triangular lattice spacing as freezing is approached, and an oriented correlation function matches these oscillations better than the radial one.","lead":"Monte Carlo simulations show how charged particles in a 2D plasma on a cylinder organize into dense layers that ripple near the edges, with the ripple spacing set by the triangular crystal the plasma forms at low temperature. The results tie these edge ripples to bulk order and suggest a new orientational correlation function for studying freezing.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Commensurate cylinder radius (Eq. 10) may pin the boundary oscillation wavelength to d0 and suppress defects, so the claimed universality is not yet established.","rationale":"The reader's weakest_assumption identifies exactly the load-bearing concern: the cylinder radius is fixed to the commensurate value R = m a0/(2π) (Eq. 10), which removes geometric frustration by construction. This is not a minor detail. The paper's headline observations—(1) boundary density oscillation wavelength converges to the triangular lattice plane spacing d0 near freezing (Sec. 4, Fig. 4), and (2) topological defects are exponentially suppressed in the crystal phase (Sec. 5)—are both quantities that the commensurate geometry directly predisposes. At zero temperature the chosen radius makes the perfect triangular lattice the exact ground state, so the wavelength of the density profile along x is d0 by construction. At finite temperature near Γm, the same geometry favors the commensurate layer spacing. Without an incommensurate control, the claim that these are universal thermodynamic properties of the 2D OCP on a cylinder is under-supported. The disk data in Fig. 4 provide partial independent evidence, but the disk has curvature-induced frustration and does not isolate the commensuration effect. The proposed test is computationally straightforward and directly settles the issue: if incommensurate radii give the same λ, ξ, and defect suppression, the concern is resolved; if not, the conclusions must be restricted. Because this is an addressable numerical check rather than a fundamental flaw, the existing CONDITIONAL verdict is appropriate and should be kept. The reader's call to require this test before full acceptance is correct, and no change to the verdict is needed.","tokens_in":10642,"tokens_out":4855,"duration_ms":40498,"concrete_test":"Repeat the Monte Carlo at the same N (e.g., N=900) and temperatures Γ=100, 120, 130, 150, but set the cylinder circumference in Eq. (10) to (m±0.3)a0 and (m±0.5)a0, keeping ρ0=1. Extract λ and ξ from the boundary density profile using the same fit (12), and compute the defect density z from Delaunay triangulation as in Sec. 5. If λ converges to d0 within ~1% and z follows the same exponential-in-Γ trend for all incommensurate offsets, the commensurate choice is harmless; if λ shifts or z is enhanced by an order of magnitude, the universal-boundary-oscillation claim fails and must be restricted to commensurate cylinders.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central quantitative claim—λ→d0 near Γm and exponential defect suppression—rests entirely on simulations with R = m a0/(2π) (Eq. 10), a circumference that is exactly m triangular-lattice spacings. At Γ→∞ the ground state is then an unstrained triangular lattice with m rows separated by d0, so the boundary oscillation wavelength and the absence of free disclinations are both enforced by the simulation box rather than emergent. The authors motivate this choice as removing geometric frustration (Sec. 2.1) and use it to claim boundary-condition independence (Fig. 5), but they never test an incommensurate radius. The disk comparison (Fig. 4) is suggestive, yet the disk is frustrated and does not isolate the effect. If the commensuration pins the oscillation wavelength to d0 or suppresses defect nucleation, then the statements 'λ rapidly converges to the triangular lattice value' (Sec. 4) and 'exponentially suppressed' defects (Sec. 5) are properties of the chosen box, not generic thermodynamic features; the oriented-correlation improvement (Sec. 6) would inherit the same bias. This is the load-bearing weak point: the universality conclusion rests on a single, specially chosen circumference.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports Monte Carlo simulations of the two-dimensional one-component plasma (2D OCP) on a cylinder, focusing on the boundary density profile, its temperature dependence, and its relation to bulk correlations. The authors find damped oscillations in the density profile whose wavelength approaches the triangular-lattice row spacing d0 as the inverse temperature Γ approaches the melting value Γm ≈ 140, while the damping length increases sharply. They argue that the cylindrical geometry removes geometric frustration and that boundary conditions do not affect the universal shape of the profile. They also quantify topological defects via Delaunay triangulation, report exponential suppression of dislocations in the crystalline phase, and propose an oriented correlation function gx(x) that better matches the boundary density oscillations than the radial g(r). The results are interpreted through a phase-field crystal (PFC) model with a single-peak approximation for the liquid structure factor.","tokens_in":10938,"tokens_out":4361,"duration_ms":41512,"significance":"If the main quantitative claims hold, the paper provides a useful step toward connecting boundary density oscillations of the 2D OCP with the crystallization transition and with anisotropic generalizations of PFC models. The cylindrical geometry is a natural setup for separating boundary-parallel and boundary-normal fluctuations, and the parameter-free droplet-squeezing identity in Eq. (13) is a clean analytic result. The paper also contains an interesting proposal for an oriented correlation function that may be more relevant than the radial distribution for wall-bounded plasmas. However, the central quantitative statements rely on fits without reported uncertainties and on simulations performed at a single, specially commensurate cylinder radius; these limitations currently prevent the claims from being fully established.","major_comments":[{"comment":"The wavelength λ and damping length ξ extracted from density profiles are presented without any error bars, fit ranges, or statistical uncertainties. Given that the claimed differences between disk and cylinder data are of order a few percent in λ and that ξ varies by about a factor of four over the plotted Γ range, the absence of uncertainties makes the statements 'rapidly converges' and 'sharper increase' impossible to assess quantitatively. Please provide error bars, the fitting procedure, the number of independent runs, and the systematic uncertainty associated with the choice of fit window.","section":"Section 4, Fig. 4"},{"comment":"All cylinder simulations use R = m a0/(2π), a circumference that is exactly commensurate with m triangular-lattice spacings. At low temperature the crystalline ground state then has row spacing d0 by construction, so the wavelength of boundary density oscillations and the suppression of disclinations are geometrically favored by the simulation box. The paper claims these features are generic for a cylinder (Section 4.1 and Discussion), but it never tests an incommensurate radius or varies m at fixed Γ. A test with incommensurate R, or at least a quantitative argument for why Eq. (10) does not pin λ and defect densities, is needed to support the universality claim.","section":"Section 2.1, Eq. (10)"},{"comment":"The claim that dislocation density is exponentially suppressed for Γ > Γm is supported only by an inset without error bars, a functional fit, or finite-size analysis. The text also states that no strong divergence is observed near the transition, but no statistical power or system-size scaling is given. Please quantify z(Γ) with fits and uncertainties, and discuss how the commensurate boundary could affect the nucleation of defects.","section":"Section 5, Fig. 7"},{"comment":"The oriented correlation function gx(x) is described only verbally and in a figure caption; no defining equation is given. The agreement between its oscillation wavelength and that of the boundary density profile is assessed visually, without quantitative values or uncertainties. Since the proposed anisotropic PFC generalization is a stated motivation, a quantitative comparison of the wavelengths extracted from gx(x), g(r), and ρ(x) is necessary to support the claim of 'much more closely' matching.","section":"Section 6, Figs. 9 and 10"}],"minor_comments":[{"comment":"The left panel axis label 'Γh0' appears to be a typo and should likely be simply Γ.","section":"Section 4, Fig. 4"},{"comment":"The sentence comparing 'a0 ≈ 1.07 and d0 ≈ 0.93' in the q0 fit is confusing: a0 is the triangular lattice spacing and d0 the row spacing, and it is unclear which reference value 2π/q0 is being compared to. Please clarify.","section":"Section 6, Eq. (20)"},{"comment":"The sentence 'ρ(x) = ρ(x) represents a one-dimensional density profile' contains a typo; it should state that the density depends only on the coordinate x under symmetric boundary conditions.","section":"Section 2, Eq. (7)"},{"comment":"No data or code availability statement is included; given that the central results are numerical, making the simulation data and analysis scripts available would substantially improve reproducibility.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the scope of the journal and the research direction is timely, but the quantitative claims need stronger support. In particular, the commensuration issue raised in Eq. (10) should be addressed with additional simulations, and all fitted quantities need error bars. If the authors provide these, the paper could become acceptable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid numerical paper with a real new idea—the oriented correlation function gx—and a cleaner geometry for the 2D OCP edge. The central picture, that boundary density oscillations have wavelength approaching the triangular spacing d0 and damping length growing near freezing, is almost certainly correct. The weak spots are quantitative: fits without error bars, an asserted exponential law that never gets fitted, and an untested commensurate radius that could bias the universality claim.\n\nWhat's new: the cylinder with the frustration-free radius R = m a0/(2π), mixed boundary conditions that let you separate normal and parallel fluctuations, and the oriented correlation function that aligns snapshots to the nearest particle. The droplet-squeezing shift δ = d0/2 in Eq. (13) is a neat parameter-free geometric identity. The defect-density curves from Delaunay triangulation are a useful complement to the density data.\n\nWhat's done well: the paper compares soft, hard, and mixed boundary conditions, and disk versus cylinder, and shows that the density profiles largely collapse. The thermalization discussion is honest. The PFC treatment is used as a fitting framework, not over-sold as a derivation. The text is careful to mark numerical results as numerical.\n\nSoft spots, in order: (1) No error bars or fit uncertainties anywhere. For a paper whose headline claims are quantitative, this is a real gap. (2) The exponential suppression of topological defects is stated without showing a fit; a slope and range would make it a claim rather than an impression. (3) The commensurate radius is a legitimate confound. At large Γ the ground state is an unstrained triangular lattice by construction, so the absence of free disclinations is partially built in. An incommensurate-radius simulation, even at one Γ, would settle whether the defect suppression and the λ→d0 convergence are generic. I don't think this is fatal—local packing near a wall should still prefer d0 spacing—but the paper should test it. (4) The conclusion favoring weak first-order melting from the absence of a defect divergence is, as the authors admit, numerically elusive. That claim should be soft-pedaled.\n\nWho it's for: people working on 2D Coulomb systems, FQHE edge structure, and phase-field crystal models. It deserves a serious referee. My recommendation: send to review, ask for error bars, data release, and at least one incommensurate radius run.","headline":"Worth refereeing: a genuinely cleaner cylinder setup and a new oriented correlation function, but the fitted claims need error bars and a test of the commensurate radius.","tokens_in":11414,"tokens_out":3623,"would_cite":true,"duration_ms":34314,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82B21","82B26","82B80"],"pacs":[],"model":"deepseek-v4-flash","headline":"On a cylinder, the boundary density oscillations of the 2D one-component plasma have a wavelength that locks to the triangular lattice-plane spacing $d_0 \\approx 0.93$ as $\\Gamma$ approaches the freezing transition, while the damping…","keywords":["two-dimensional one-component plasma","boundary density oscillations","cylindrical geometry","triangular lattice","freezing transition","phase-field crystal","oriented correlation function","topological defects"],"falsifier":"Run the same Monte Carlo simulations on cylinders with incommensurate radii, for example a value midway between two consecutive commensurate values, at $\\Gamma = 130$; if the fitted wavelength $\\lambda$ shifts away from $d_0 \\approx 0.93$ or the dislocation density changes sharply, then the commensurate choice is pinning the observed boundary physics rather than revealing a thermodynamic property.","tokens_in":10431,"feed_emoji":"🌊","tokens_out":8995,"duration_ms":46630,"temperature":0.7,"pith_summary":"This paper uses Monte Carlo simulations of a two-dimensional one-component plasma wrapped on a cylinder to establish a sharp picture of boundary density oscillations in the strongly coupled fluid. It argues that, as the inverse temperature $\\Gamma$ approaches the freezing transition $\\Gamma_m \\approx 140$, the wavelength $\\lambda$ of the ripples near the boundary converges rapidly to $d_0 \\approx 0.93$, the spacing between triangular lattice planes, while the damping length $\\xi$ increases sharply. Because the cylinder circumference is chosen to fit an integer number of lattice spacings, the geometry removes frustration and disclinations, so these oscillations can be attributed to freezing at the edge rather than to geometric artifacts. The paper also proposes an oriented correlation function $g_x(x)$, defined by aligning each snapshot with the particle closest to the test particle, which matches the boundary-density wavelength much better than the radial $g(r)$. This connects boundary profiles to bulk correlations and to the melting transition through a phase-field crystal description.","feed_headline":"2D Coulomb plasma edge ripples lock to lattice spacing near freezing","feed_subtitle":"Their wavelength converges to d0≈0.93 and damping length surges as Γ approaches freezing near 140.","key_machinery":"The load-bearing object is the unfrustrated triangular lattice on the cylinder, obtained by fixing $R = m a_0/(2\\pi)$ so that $m$ spacings $a_0$ fit around the circumference; the lattice-plane spacing $d_0 = \\sqrt{3}a_0/2 \\approx 0.93$ is the wavelength that boundary density oscillations converge to. The oscillations are quantified by fitting $\\rho(x)$ to $e^{-x/\\xi}\\sin(2\\pi x/\\lambda)$, producing the damping length $\\xi$ and wavelength $\\lambda$. The phase-field crystal free energy $F[\\psi] = \\int d^2x\\,( \\frac{a}{2}\\psi^2 + \\frac{b}{2}\\psi(\\Delta + q_0^2)^2\\psi + \\frac{u}{4}\\psi^4 + V\\psi )$, with parameters linked to the bulk direct correlation function through the Ornstein-Zernike relation, connects the boundary profile to the static structure factor peak. The paper's new instrument is the oriented correlation function $g_x(x)$, defined by rotating each snapshot so the particle closest to a test particle lies along the $y$-axis; the density along $x$ then tracks the boundary profile more faithfully than $g(r)$. Together these objects carry the argument that edge oscillations are freezing signatures and that an anisotropic PFC action is needed.","core_discovery":"The central claim is that damped oscillations in the boundary density profile of the two-dimensional one-component plasma on a cylinder are a thermodynamic signature of a crystalline layer freezing at the edge. With the radius fixed to $R = m a_0/(2\\pi)$, exactly $m$ triangular-lattice spacings $a_0$ fit around the circumference, so a perfect triangular lattice can form without disclinations; the lattice-plane spacing $d_0 = \\sqrt{3} a_0/2 \\approx 0.93$ sets the oscillation wavelength. Fitting $\\rho(x)$ to $e^{-x/\\xi}\\sin(2\\pi x/\\lambda)$, the paper finds $\\lambda$ converging to $d_0$ as $\\Gamma \\to \\Gamma_m$ and $\\xi$ increasing sharply, consistent with a diverging correlation length at freezing. Comparing soft-wall, hard-wall, and mixed boundary conditions shows the profile shape is universal once shifted to the classical boundary, and the droplet squeezing shift is $\\delta = d_0/2$. In the crystal phase, Delaunay triangulation shows the dislocation density $z$ is exponentially suppressed for $\\Gamma > \\Gamma_m$, with no strong divergence near the transition. Finally, the paper defines an oriented correlation function $g_x(x)$ whose oscillation wavelength agrees with the boundary-density wavelength far better than the radial $g(r)$, pointing to an anisotropic generalization of the phase-field crystal model.","pith_inferences":["Editorial inference: if the commensurate-radius setup is not generic, then incommensurate cylinder radii $R \\neq m a_0/(2\\pi)$ should introduce phase slips or extra disclinations; repeating the simulation at such radii would show whether the clean $\\lambda \\to d_0$ convergence and the exponential defect suppression survive.","Editorial inference: the oriented correlation function suggests a quantitative test for anisotropic phase-field crystal models, where direction-dependent $q_0$ parameters could be fitted to the same data and checked against the residual five-to-twenty percent discrepancies in $q_0$ and $a/b$ reported for the isotropic model.","Editorial inference: because the cylinder suppresses disclinations by construction, the exponential defect suppression above $\\Gamma_m$ does not by itself discriminate between weakly first-order melting and BKTHNY; measuring the hexatic order parameter or defect correlation length in the same geometry would sharpen that distinction.","Editorial inference: the droplet squeezing shift $\\delta = d_0/2$ derived from lattice mismatch on the cylinder could be searched for in disk-geometry profiles as a boundary-dependent offset, providing a quantitative bridge between the two geometries."],"forward_implications":["The wavelength of boundary density oscillations in the strong-coupling fluid is set by the triangular lattice-plane spacing $d_0 \\approx 0.93$ and is essentially independent of $\\Gamma$, so measuring edge ripples gives the lattice constant of the incipient crystal.","The damping length $\\xi$ increases sharply as $\\Gamma$ approaches $\\Gamma_m$, consistent with a correlation length that diverges at freezing; the cylinder geometry makes this trend cleaner than in disk geometry.","The density-profile shape is universal across soft-wall, hard-wall, and mixed boundary conditions once shifted to the classical boundary, and the squeezing shift $\\delta = d_0/2$ follows from the lattice mismatch.","The dislocation density is exponentially suppressed for $\\Gamma > \\Gamma_m$, and the absence of a strong divergence near the transition is read as support for a weakly first-order melting picture, with the caveat that defect-divergence signatures can be numerically elusive.","The oriented correlation function $g_x(x)$ reproduces the boundary oscillation wavelength much better than $g(r)$, indicating that an anisotropic generalization of the phase-field crystal model would improve the calculation of density profiles."],"supporting_citations":[{"why":"Supplies the freezing-at-the-edge interpretation of boundary density oscillations that the cylinder simulations are designed to test and sharpen.","marker":"[40]"},{"why":"Provide exact and large-$N$ results showing the boundary density features are singular in the large-$N$ limit, motivating the cylindrical setup.","marker":"[34,41]"},{"why":"Defines the phase-field crystal free energy and saddle-point equation used to connect bulk correlations to boundary profiles.","marker":"[45]"},{"why":"Gives the direct correlation function and Ornstein-Zernike relation through which the PFC parameters $a$, $b$, $q_0$ are extracted from $g(r)$.","marker":"[46]"},{"why":"Provide the BKTHNY melting scenario whose defect-unbinding predictions the defect-density data are compared against.","marker":"[18–22]"},{"why":"Documents that BKTHNY defect-divergence signatures can be elusive in simulations, used to qualify the absence of a divergence in $z$ near $\\Gamma_m$.","marker":"[43]"},{"why":"Explains the non-perturbative droplet squeezing in the large-$N$ limit, supporting the interpretation of the shift $\\delta = d_0/2$.","marker":"[42]"}],"fun_headline_variants":["Cylinder plasma edge ripples lock to lattice spacing at freeze","2D plasma boundary ripples converge to lattice constant near melt","Plasma edge ripples mirror lattice spacing at freezing transition","Anisotropic correlation ties plasma boundary ripples to crystal order"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that choosing the cylinder radius as $R = m a_0/(2\\pi)$, so the circumference holds exactly $m$ lattice spacings, does not bias the physics; if this commensuration is what pins the wavelength to $d_0$ or suppresses defects, the universal-profile and exponential-defect-suppression claims would not survive for generic radii.","fun_headline_variants_meta":{"raw":{"variants":["Cylinder plasma edge ripples lock to lattice spacing at freeze","2D plasma boundary ripples converge to lattice constant near melt","Plasma edge ripples mirror lattice spacing at freezing transition","Anisotropic correlation ties plasma boundary ripples to crystal order"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000288,"raw_usage":{"total_tokens":1694,"prompt_tokens":957,"completion_tokens":737,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":573,"completion_tokens_details":{"reasoning_tokens":666}},"tokens_in":573,"tokens_out":737,"duration_ms":7103,"temperature":1.0,"reasoning_tokens":666,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T11:30:57.491094+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same Monte Carlo simulations on cylinders with incommensurate radii, for example a value midway between two consecutive commensurate values, at $\\Gamma = 130$; if the fitted wavelength $\\lambda$ shifts away from $d_0 \\approx 0.93$ or the dislocation density changes sharply, then the commensurate choice is pinning the observed boundary physics rather than revealing a thermodynamic property.","supporting_citations":[{"cited_title":"The boundary density profile of a coulomb droplet","cited_arxiv_id":null,"evidence_quote":"Supplies the freezing-at-the-edge interpretation of boundary density oscillations that the cylinder simulations are designed to test and sharpen."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the phase-field crystal free energy and saddle-point equation used to connect bulk correlations to boundary profiles."},{"cited_title":"Kosterlitz–thouless physics: a review of key issues.Reports on Progress in Physics, 79(2):026001, 2016","cited_arxiv_id":null,"evidence_quote":"Documents that BKTHNY defect-divergence signatures can be elusive in simulations, used to qualify the absence of a divergence in $z$ near $\\Gamma_m$."},{"cited_title":"Anomalous hydrodynamics of two-dimensional vortex fluids","cited_arxiv_id":null,"evidence_quote":"Explains the non-perturbative droplet squeezing in the large-$N$ limit, supporting the interpretation of the shift $\\delta = d_0/2$."}],"review_version":1}