{"id":"1a6e7350-ebb4-464c-a92f-f776d1175706","arxiv_id":"2608.07644","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Hard-sphere nucleation times increase by roughly 6-fold per 0.1% volume fraction below about 53%, making spontaneous coexistence dynamically inaccessible over most of the coexistence window.","lead":"This reply argues that hard-sphere nucleation times grow exponentially as volume fraction moves a few tenths of a percent below about 53%, so spontaneous fluid-crystal coexistence is essentially unreachable over most of the coexistence region. It corrects a Comment's overgeneralization from a single near-melting simulation point.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Table 1's constant 6-fold/0.1% slope is fit at the melting end and cannot be extended to φ≈0.50; CNT curvature makes extrapolated times even longer, so the qualitative claim survives but the table is not a reliable benchmark.","rationale":"The reader's weakest-assumption assessment identifies the same load-bearing point: Table 1 extrapolates a linear fit over 11 orders of magnitude without error bars. My reading sharpens this by noting that the constant-slope extrapolation is not merely unquantified; it contradicts the functional form expected from the very same nucleation theory the reply invokes. Near φ_F, the nucleation barrier diverges as 1/(φ-φ_F)^2, so the log-rate versus volume-fraction curve steepens rapidly. Consequently the tabulated numbers below φ≈0.52 are unsupported, and if used as predictive benchmarks they could mislead in either direction depending on what one assumes—though the CNT correction makes times longer, not shorter. The qualitative central claim—that waiting times grow astronomically just a few tenths of a percent below 53%—is independently supported by the Auer–Frenkel data, the companion near-melting simulations, and Smallenburg's own sparse crystallization statistics. Therefore the right verdict remains CONDITIONAL: the qualitative argument is compelling, but Table 1 and the supersaturation mapping need correction or explicit caveats before the quantitative entries are used. My concern does not move the verdict from the reader's CONDITIONAL assessment.","tokens_in":4202,"tokens_out":12828,"duration_ms":126716,"concrete_test":"Recover the raw data points behind Fig. 11 of Ref. [6] and fit them two ways: log I* = a - b/(φ-φ_F)^2 (CNT form) and log I* = a - bφ (the linear form used in Table 1). Then compute the predicted waiting time at φ=0.520 from both fits for the 2,048,000-particle box. If the two predictions differ by more than an order of magnitude, Table 1's 16,000-year entry is not a reliable consequence of the cited data, and the table should be relabeled or restricted to the fitted range.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim rests on Table 1, which multiplies the Auer–Frenkel waiting time at φ≈0.5342 by a constant factor of 6 per 0.1% volume-fraction decrease, all the way down to φ=0.495. The caption says this is a linear fit to Fig. 11 of Ref. [6], but no fit details, error bars, or residuals are given. More fundamentally, classical nucleation theory for hard spheres predicts log I* ∝ -A/(φ-φ_F)^2 (with φ_F≈0.4918), so the logarithmic slope is not constant: it grows roughly as (φ-φ_F)^{-3} as the freezing point is approached. A constant slope fitted near 0.534 will therefore misrepresent rates at 0.52 and below. The direction of the error matters: the CNT curvature makes the true waiting times even longer than Table 1 reports, so the headline claim that spontaneous coexistence is dynamically inaccessible through most of the coexistence region is strengthened, not undermined. But the specific tabulated values—e.g., 16,000 years at φ=0.520 and 'never' at φ=0.495—are not a valid extrapolation of the cited linear fit. The 20%-supersaturation mapping to φ=0.541/0.536 is also asserted without derivation, although this matters less for the table than for locating the onset of accessibility.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript is a Reply to a Comment by Smallenburg on the authors' recent Perspective. The Reply argues that Smallenburg's conclusion that equilibrium hard-sphere coexistence is 'readily achievable' in simulation is an overgeneralization from a single high-density state point (φ=0.5325). The authors reiterate their central claim that nucleation times in monodisperse, purely repulsive hard spheres grow exponentially as the volume fraction decreases from about φ=0.53, based on Auer and Frenkel's measured nucleation rates. They present Table 1, which lists waiting times from 1.5 seconds at φ=0.535 to 'never' at φ=0.495, obtained by extrapolating a 6-fold increase per 0.1% volume-fraction decrease. They also propose a mapping between atomic supersaturation and colloidal volume fraction, identifying 20% supersaturation with φ≈0.541 or 0.536, and they critique Smallenburg's simulation statistics as not establishing the theoretical 80% crystal fraction.","tokens_in":4501,"tokens_out":7998,"duration_ms":70729,"significance":"If the exponential sensitivity claim is correct, it has practical consequences for the design of hard-sphere simulations: unbiased spontaneous nucleation would be feasible only in a narrow window very close to melting. The Reply credibly emphasizes that even a modest shift in volume fraction changes nucleation times by many orders of magnitude, a point grounded in established data (Auer and Frenkel, Ref. [6]). The manuscript also honestly acknowledges that the original Perspective was unclear about the range over which the long-time claim applies. However, the quantitative predictions in Table 1 and the supersaturation mapping are not adequately supported, and these are the load-bearing elements of the Reply's quantitative argument.","major_comments":[{"comment":"Table 1 is presented as the central quantitative evidence for 'astronomically long' nucleation times, but it is based on a linear fit to Fig. 11 of Ref. [6] with no documented fit range, no residuals, and no error bars. The table extrapolates a constant 6-fold-per-0.1% factor over 11 orders of magnitude, from φ=0.5342 to φ=0.495. Classical nucleation theory predicts that the logarithmic slope should increase as the freezing point is approached, so a constant-slope fit from the melting-side data is not a reliable predictor at lower volume fractions. The specific entries in Table 1 should therefore be labeled as an illustrative extrapolation, or the table should be restricted to the volume-fraction range over which the cited linear fit is actually valid.","section":"Table 1 and the paragraph beginning 'Now, Frenkel's study further shows...'"},{"comment":"The mapping %supercooling = (φ−φ_F)/(φ_HCP−φ_F) is asserted without derivation or comparison to any direct measurement. The physical analogy is not self-evident: Brownian motion does not vanish at φ_HCP, and the choice of φ_HCP as the analogue of absolute zero is not justified. This mapping is used to identify 20% supersaturation with φ≈0.541 (or 0.536) and underlies the statement that spontaneous phase separation is dynamically accessible only down to about 53% volume fraction. In its current form the mapping is an unsupported assertion; the Reply should either provide a derivation and justification or rely on the direct Auer–Frenkel data instead.","section":"Paragraph beginning 'It is straightforward to translate...'"},{"comment":"The Reply gives two widely different laboratory-time estimates for the authors' simulation box: about 2 years using the atomic nucleation rate of 1 nucleus cm⁻³ s⁻¹ at 20% supersaturation, and 9.1 seconds using the Auer–Frenkel colloidal rate at φ=0.5342. The discrepancy spans several orders of magnitude and is not reconciled. The phrase 'faster for colloids' is insufficient to explain the gap, and it is not clear which estimate is appropriate for which volume fraction. The Reply should clarify the relationship between the atomic and colloidal rates and avoid juxtaposing the two estimates without explicit reconciliation.","section":"Paragraphs 2 and 4 (laboratory-time estimates)"},{"comment":"The Reply dismisses the significance of Smallenburg's simulation results by noting that only 8 of 50 slab runs and 11 of 50 cubic runs produced crystallization, with final crystal fractions ranging from 56% to 100%. However, no statistical context is provided: there is no definition of what constitutes 'any crystallization,' no information on simulation duration, and no error bars on the fractions. Without such context, the statement that the results 'do not establish that the 80% theoretical value ... has been achieved consistently or spontaneously' is itself unsupported. This critique is a secondary point in the Reply but, as written, it does not meet the rigor expected for dismissing a commentator's simulation evidence.","section":"Penultimate paragraph (critique of Smallenburg's data)"}],"minor_comments":[{"comment":"The manuscript contains several typographical errors and formatting glitches, including 'reportsnear-melting pointsimulations' (missing spaces), '10 −16cm3' (missing superscripts), 'ϕ= 0.5342%of' (stray percent sign), and 'Table inFigure 1shows' (no Figure 1 is actually included; only Table 1 appears). These should be corrected before publication.","section":"Throughout"},{"comment":"The text says 'Independent calculations cited by ten Wolde et al.' but the reference given is to Kelton (Ref. [5]), not to ten Wolde et al. This should be clarified to indicate the original source of the φ=0.536 value.","section":"Paragraph 2, reference [5]"},{"comment":"The caption states that the table is 'calculated based on a linear fit to the data from Figure 11 in Ref. [6]' but does not give the fit parameters, the range of the fit, or the goodness of fit. Even if the fit details remain in the text, the caption should at least state the confidence interval or the range of validity.","section":"Table 1 caption"},{"comment":"The Reply uses the term 'supercooling' for colloids, but the relevant variable is volume-fraction supersaturation. Consider using a consistent term such as 'supersaturation' throughout, and number the defining equation for clarity.","section":"Terminology"},{"comment":"The phrase 'spotty results' in the penultimate paragraph is informal and should be replaced with a more neutral description in a formal journal reply.","section":"Language"}],"recommendation":"major_revision","confidential_remarks":"This is a reply to a Comment, so the scope is inherently narrow and the manuscript's fit with the journal depends on the journal's policy for exchange letters. The Reply's qualitative message is likely correct, but the quantitative table and the supersaturation mapping are not sufficiently supported. The manuscript's heavy reliance on the authors' own prior work (Refs. [1] and [3]) is expected in a reply, but the presentation would benefit from a more neutral framing. The editor may wish to consider whether the Reply should be allowed to include an extrapolated table without error analysis, or whether it should be reframed as an illustrative estimate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The reply makes a solid qualitative point: Smallenburg's single simulation at φ=0.5325, right next to melting, does not license the claim that coexistence is \"readily achievable\" throughout the coexistence region. The exponential growth of nucleation times as you move away from 53% is real—Auer and Frenkel measured it—and the paper is right to push back on that overgeneralization. It also credits its own earlier Perspective for being unclear on this point, which is honest.\n\nWhat's useful here: the paper frames φ≈0.53 as the practical departure point for accessible spontaneous nucleation, translates Auer and Frenkel's reduced rates into dimensional lab times for a colloidal system, and works through what that means for simulation design. That interpretive work is a service, even if the underlying law is not new.\n\nNow the soft spots, in proportion. Table 1 is the main problem. It takes a constant \"6-fold per 0.1%\" slope, fit near φ=0.534, and extends it linearly down to φ=0.495—eleven orders of magnitude. No fit details, no error bars, no residuals. Worse, the constant-slope assumption is inconsistent with classical nucleation theory, which gives log I* ∝ -A/(φ-φ_F)^2, so the logarithmic slope grows as you approach freezing. That means the true waiting times at φ=0.52 and below are even longer than the table says. So the qualitative conclusion—that spontaneous phase separation is dynamically inaccessible through most of the coexistence region—survives and is actually strengthened. But the specific numbers in Table 1 are not a reliable benchmark. They should be presented as illustrative, not predictive, or the fit should be replaced with a proper CNT-based extrapolation with uncertainties.\n\nTwo smaller issues. The supersaturation mapping \"% supercooling = (φ-φ_F)/(φ_HCP-φ_F)\" is asserted without derivation; it may be reasonable, but it needs a citation or a few lines of justification. And for a Reply, the paper never actually cites Smallenburg's Comment in the reference list. That's a noticeable omission, even if the journal version will include it.\n\nWho is this for? People working on hard-sphere or colloidal phase behavior and anyone designing coexistence simulations. It's a short, argumentative reply, not a major new result. But the caution it delivers is important, and the qualitative argument is grounded in prior published rates.\n\nRecommendation: it deserves peer review as a reply, because it corrects a potential misconception and the central claim holds up. A serious referee should, however, insist that Table 1 be reframed or recomputed before publication. The paper is fine to publish after that revision.","headline":"A mostly correct reply that overreaches in its table: the qualitative exponential-growth claim stands, but the constant-slope extrapolation is not a valid predictive benchmark.","tokens_in":5021,"tokens_out":2142,"would_cite":false,"duration_ms":21711,"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":"Hard-sphere crystallization waits grow 6-fold per 0.1% density drop, making spontaneous coexistence reachable only near melting.","keywords":["hard-sphere colloids","nucleation barrier","volume fraction","entropy-exchange mechanism","crystal nucleation","fluid-crystal coexistence","supersaturation","simulation timescales"],"falsifier":"A pristine, unbiased hard-sphere or nearly hard-sphere colloid simulation at φ=0.52 that produces a stable crystal nucleus within, say, a month of simulated time would contradict the predicted 16,000-year waiting time; conversely, rate measurements at φ=0.525, 0.52, and 0.515 that show the 6-fold-per-0.1% factor persisting would support the extrapolation.","tokens_in":3980,"feed_emoji":"⏳","tokens_out":6127,"duration_ms":53558,"temperature":0.7,"pith_summary":"This reply to a Comment defends a sharp claim about pristine hard-sphere simulations: spontaneous fluid-crystal coexistence is dynamically accessible only in a narrow window very close to the melting volume fraction, around 53%. The author argues that the Comment's successful nucleation run at a single near-melting packing fraction is the easy limiting case, not evidence that coexistence is readily achievable across the coexistence region. The central evidence is an exponential growth law: the nucleation time increases by roughly a factor of six for every 0.1% decrease in volume fraction, rising from about 1.5 seconds at φ=0.535 to 16,000 years at φ=0.520 and effectively infinite by φ≈0.505. If this is right, the long-observed absence of spontaneous coexistence in simulations is the expected outcome of the entropy-exchange mechanism, not a numerical artifact.","feed_headline":"Crystallization waits jump 6-fold per 0.1% density drop","feed_subtitle":"Hard-sphere simulations can nucleate spontaneously only in a narrow sliver near the melting point.","key_machinery":"The load-bearing object is Frenkel's entropy-exchange mechanism, defined as the competition between configurational and vibrational entropy that sets the nucleation barrier and makes it rise sharply as supersaturation increases. The reply quantifies this with a mapping from the atomic Lennard-Jones system to colloidal hard spheres, treating the hexagonal-close-packed volume fraction φ_HCP=0.74 as the analog of absolute zero, with melting at φ_M=0.543 and freezing at φ_F=0.4918; this maps 20% supersaturation to φ≈0.541, or φ≈0.536 by an independent estimate. The exponential growth law itself comes from a linear fit to the reduced nucleation rates measured in the 2004 hard-sphere study, giving roughly a 6-fold increase in nucleation time per 0.1% decrease in volume fraction, which the reply extrapolates across the coexistence region.","core_discovery":"The paper's central claim is that the nucleation time of monodisperse, purely repulsive hard spheres grows exponentially as the system moves away from the melting point: a 0.1% reduction in volume fraction multiplies the waiting time by about six. Using the absolute crystallization rate measured near φ=0.534 as a baseline, the author converts that growth law into a table of laboratory waiting times for a colloidal-scale simulation, showing that nucleation takes about 1.5 seconds at φ=0.535, 2.9 hours at φ=0.53, 16,000 years at φ=0.52, and becomes practically unattainable at lower densities. The same entropy-exchange mechanism that produces this barrier also explains why the Comment's simulation at φ=0.5325 succeeded: that state point sits just inside the narrow near-melting window where the barrier is low enough for spontaneous nucleation on feasible timescales. The reply concludes that the Comment's broader takeaway, that equilibrium coexistence is 'readily achievable,' is an overgeneralization from a single favorable state point.","pith_inferences":["The 6-fold-per-0.1% law is asserted from a linear fit over an 11-order-of-magnitude extrapolation; a direct measurement of nucleation rates at intermediate volume fractions, say φ=0.525 to 0.53, would test whether the growth rate itself changes with supersaturation.","The supersaturation mapping crucially depends on assigning HCP as the zero-motion reference for colloids; if that reference were revised, the nominal 20% supersaturation point would shift, moving the predicted fast-nucleation window without changing the exponential structure.","If the entropy-exchange argument is correct, reported spontaneous coexistence in other entropy-dominated systems should likewise cluster near their melting points, and any claim of rapid coexistence far from melting would warrant scrutiny.","The reply's table assumes a fixed system volume; because nucleation rates scale with volume, a larger simulation box could in principle bring lower-volume-fraction nucleation into reach, a testable prediction."],"forward_implications":["Spontaneous, unbiased coexistence simulations of hard spheres are feasible only above roughly φ=0.53; below that, waiting times exceed what any current simulation can reach.","A single near-melting state point cannot be used to infer that equilibrium coexistence is generally easy to reach; the Comment's own data (8 of 50 slab runs and 11 of 50 cubic runs crystallized) support the narrow-window picture.","Colloidal-scale simulations nucleate far faster than atomic-scale ones near melting, about 1.5 seconds versus years, but retain the same exponential sensitivity to volume fraction.","The lever-rule prediction of about 80% final crystal fraction requires simulation durations much longer than used in the Comment, so scattered final fractions likely reflect incomplete convergence rather than the equilibrium value."],"supporting_citations":[{"why":"Supplies the entropy-exchange mechanism that sets the nucleation barrier.","marker":"[2]"},{"why":"Supplies the atomic nucleation rate of 1 nucleus per cm3 per second and the 20% supersaturation reference point.","marker":"[4]"},{"why":"Supplies an independent estimate placing 20% supersaturation at φ=0.536.","marker":"[5]"},{"why":"Supplies the measured reduced nucleation rate at φ=0.5342 and the data from which the 6-fold-per-0.1% growth factor is fitted.","marker":"[6]"},{"why":"Reports the companion simulation showing spontaneous phase separation at φ=0.535, which serves as the baseline for the waiting-time table.","marker":"[3]"},{"why":"Provides the original Perspective's framing of the missing coexistence state that this reply defends.","marker":"[1]"}],"fun_headline_variants":["Nucleation waits grow 6x per 0.1% density fall","Hard-sphere nucleation feasible only near melting point","Crystallization times explode away from melting point"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The exponential waiting-time table assumes that the growth rate measured near φ=0.534 continues unchanged all the way down to φ=0.495, over an 11-order-of-magnitude range, with no error bars on the fit.","fun_headline_variants_meta":{"raw":{"variants":["Nucleation waits grow 6x per 0.1% density fall","Hard-sphere nucleation feasible only near melting point","Crystallization times explode away from melting point"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000196,"raw_usage":{"total_tokens":1346,"prompt_tokens":919,"completion_tokens":427,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":535,"completion_tokens_details":{"reasoning_tokens":373}},"tokens_in":535,"tokens_out":427,"duration_ms":4404,"temperature":1.0,"reasoning_tokens":373,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T00:26:42.204464+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A pristine, unbiased hard-sphere or nearly hard-sphere colloid simulation at φ=0.52 that produces a stable crystal nucleus within, say, a month of simulated time would contradict the predicted 16,000-year waiting time; conversely, rate measurements at φ=0.525, 0.52, and 0.515 that show the 6-fold-per-0.1% factor persisting would support the extrapolation.","supporting_citations":[{"cited_title":"Order through disorder: entropy strikes back.Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the entropy-exchange mechanism that sets the nucleation barrier."},{"cited_title":"Numerical calculation of the rate of crystal nu- cleation in a lennard-jones system at moderate undercool- ing.The Journal of Chemical Physics, 104(24):9932– 9947, 1996","cited_arxiv_id":null,"evidence_quote":"Supplies the atomic nucleation rate of 1 nucleus per cm3 per second and the 20% supersaturation reference point."},{"cited_title":"Crystal nucleation in liquids and glasses","cited_arxiv_id":null,"evidence_quote":"Supplies an independent estimate placing 20% supersaturation at φ=0.536."},{"cited_title":"Numerical prediction of absolute crystallization rates in hard-sphere colloids.The Journal of Chemical Physics, 120(6):3015–3029, 2004","cited_arxiv_id":null,"evidence_quote":"Supplies the measured reduced nucleation rate at φ=0.5342 and the data from which the 6-fold-per-0.1% growth factor is fitted."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the companion simulation showing spontaneous phase separation at φ=0.535, which serves as the baseline for the waiting-time table."},{"cited_title":"The elusive fluid-and-crystal coexis- tence state in simulations of monodisperse, hard-sphere colloids.AIChE Journal, 72(6):e70275, 2026","cited_arxiv_id":null,"evidence_quote":"Provides the original Perspective's framing of the missing coexistence state that this reply defends."}],"review_version":1}