{"id":"1d2dfb56-b61e-46f1-9ba4-4fd055885fce","arxiv_id":"2412.18289","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Chiral active crystals are predicted to acquire a net angular momentum with a non-monotonic dependence on chirality, plus a non-dispersive spectral peak and a rotational entropy-production contribution.","lead":"This paper studies a two-dimensional crystal made of chiral active particles that swim in circles. It predicts the crystal develops a net angular momentum, strongest at intermediate chirality, and that this rotation shows up in particle vibrations and in entropy production.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (27) is derived in the overdamped limit, but for the paper's own parameters most phonon modes are underdamped; no simulation or calculation shows the nonmonotonic M survives finite inertia.","rationale":"The paper's formalism is clean and the derivation of Eq. (27) from Eq. (26) via the Debye approximation is internally consistent; I found no algebraic error in the small-Ω and large-Ω limits. The load-bearing weak point is the regime in which the central observable is evaluated. The model is introduced as underdamped (Eqs. (1)-(2)), and the full spectral density (11) is underdamped, but the angular-momentum integral leading to Eq. (27) is performed with the overdamped propagator iωγ+ω_q² (Eq. (77)). For the parameters used in the figures (K/(mγ²)=10^3, τγ=1) the modes that dominate the integral have ω_q ~ γ, so the overdamped condition is not met. This is not a matter of anharmonicity or defects; it is a mismatch within the stated model. A direct numerical integration of the exact underdamped spectral density or a Langevin simulation would settle the issue. If finite-inertia corrections only shift the peak, the qualitative conclusion stands; if they change the sign or monotonicity, the central claim fails. The reader's harmonic/defect concern is real but secondary; the overdamped reduction is a more immediate and testable threat.","tokens_in":30899,"tokens_out":19608,"duration_ms":182892,"concrete_test":"Simulate N≈4×10^4 particles obeying the full underdamped Langevin equations (1)-(2) on a triangular lattice with harmonic springs, using the Fig. 4 parameters (K/(mγ²)=10^3, τγ=1, mv0²/T=10^2). Measure the steady-state angular momentum M=(m/N)Σ_n r_n×v_n as a function of Ω/γ over 0.1-100 and compare with Eq. (27). Also integrate Eq. (11) over ω and q numerically to obtain the exact underdamped M_q without the overdamped reduction. If the simulated or exact M(Ω) deviates substantially from Eq. (27) or loses nonmonotonicity, the central claim needs revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative prediction, Eq. (27), follows from Eq. (26), which is obtained from the overdamped equation of motion (77). That reduction requires γ²/4 ≫ ω_q² for every mode contributing to the q-integral. For the representative parameters used in Fig. 4 (K/(mγ²)=10^3, τγ=1), the integrand of Eq. (26) peaks near qσ ≈ 1/√(τγ · (3/2)(K/(mγ²))) ≈ 0.026, where ω_q ≈ γ; the overdamped inequality is violated by roughly an order of magnitude (ω_q²/γ² ≈ 1, versus the required ≪ 1/4). Thus the dominant modes contributing to the angular momentum are not in the overdamped regime in which Eq. (26) was derived. The paper asserts in Sec. 4 that chirality effects are similar in both regimes, but this is not demonstrated for M; the exact underdamped spectral density Eq. (11) has phonon poles at ±ω_q that could change the q-dependence of the integrated M. Since the model itself is underdamped, Eqs. (1)-(2), the headline nonmonotonic curve Fig. 4(c) is not yet established by the stated derivation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"Marini Bettolo Marconi and Caprini study a two-dimensional triangular harmonic crystal composed of underdamped chiral active particles, with the active force following chiral active Ornstein--Uhlenbeck dynamics. Using a Fourier-space solution of the linear Langevin equations, they derive displacement and velocity correlation spectra, the spectral density of angular momentum, equal-time and time-dependent correlations, and a path-integral decomposition of the entropy production rate. The central prediction is Eq. (27): the total angular momentum per particle M is nonzero for nonzero chirality Ω, vanishes as Ω→0 and Ω→∞, and is non-monotonic in Ω. The paper also reports a non-dispersive peak at the chiral frequency in the displacement spectrum, a chirality-reduced velocity correlation length, and an additional rotational contribution to entropy production.","tokens_in":31190,"tokens_out":25306,"duration_ms":222384,"significance":"If the central prediction is correct, this is a valuable exactly solvable model of collective rotation emerging from individual chirality in a two-dimensional active solid. The analytic treatment is largely self-contained and the torque-balance relation in Appendix B (Eq. (74)) is exact, so the angular momentum is not inserted by hand; the non-monotonic M(Ω) of Eq. (27) is a concrete, falsifiable prediction. The paper also provides explicit closed-form results for the displacement spectrum, velocity correlations, mean-square displacement, and entropy-production decomposition, which can guide simulations and experiments with chiral colloids or granular particles. The harmonic, defect-free idealization is appropriate for an analytic theory, but it limits quantitative comparison with real active crystals, where anharmonicities and plastic events may modify the predicted angular momentum.","major_comments":[{"comment":"The time-domain cross-correlation in Eq. (3b) is written as −m²γ²v0² δ_{n,n′} e^{−|t−t′|/τ} sin(Ω|t−t′|). As written, this function is even in the time difference, so its Fourier transform is real and even in ω, in contradiction with the frequency-domain expression Eq. (63b) (and Eq. (69)), which is imaginary and odd in ω. The correct steady-state correlation from Eq. (2) is −m²γ²v0² δ_{n,n′} e^{−|t−t′|/τ} sin(Ω(t−t′)); the absolute value should appear only in the exponential decay factor. This is not a cosmetic issue: the oddness of the cross-correlation in the time difference is exactly what produces the imaginary odd spectrum and hence the nonzero angular momentum. Please correct Eq. (3b) and check any passage that quotes or relies on it.","section":"§4.2, Eq. (26)-(27), and Appendix D, Eq. (77)"}],"minor_comments":[{"comment":"The small- and large-Ω asymptotics stated after Eq. (27) appear to contain an extra factor 1/(2π) compared with the limiting behavior of Eq. (27); for example, the small-Ω limit of Eq. (27) is (vc/2π) m v0² τ [Ωτ q_D²/(1+τ/γ c²q_D²)], while Eq. (28) has (vc/2π) m v0² τ (1/(2π)) [Ωτ q_D²/(1+τ/γ c²q_D²)]. Please verify these prefactors.","section":"Eqs. (28)-(29)"},{"comment":"The sentence 'vertical dashed lines are used to denote the peak frequencies, i.e. Ω/γ and' is incomplete; it should finish with '±ω_q/γ' (or equivalent).","section":"Figure 2 caption"},{"comment":"The thermal noise is denoted ζ_n in Eq. (1) but ξ_n in Eq. (4) and in most of the rest of the paper; please unify the notation.","section":"Eq. (1) and throughout"},{"comment":"The phrase 'i.e., for time i.e., for time t → ∞' contains a duplication and should be reworded.","section":"Section 4, first paragraph"},{"comment":"Several typos remain, e.g., 'diferent' in Section 5.1 and 'interpretaion' in Section 4.3; a careful proofreading pass is needed.","section":"General"},{"comment":"The Debye cutoff q_D is introduced without a precise definition of the corresponding Debye frequency, and Eq. (27) depends explicitly on q_D; the statement in Section 1 that the method is 'without parameter fitting' should be softened or q_D should be fixed by the lattice, e.g., by the Brillouin-zone edge.","section":"Introduction and Eq. (27)"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the harmonic idealization is acceptable for an analytic theory, but the underdamped/overdamped gap and the sign error in Eq. (3b) need to be addressed before publication. The second issue is a simple fix, while the first requires additional calculation or simulation support."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper gives a clean exact solution for a chiral AOUP harmonic crystal and identifies a genuinely new observable, the steady-state angular momentum, with a plausible nonmonotonic dependence on chirality. But the headline curve is computed in an overdamped reduction that the paper's own parameters violate, so the central quantitative claim is not actually established.\n\nWhat's good: the model extends the authors' earlier AOUP crystal work to chiral active forces; the algebra is internally consistent (frequency-domain active-force correlations match time-domain expressions); it reproduces the Shee et al. velocity correlation length; and it exposes a non-dispersive spectral peak at the chirality frequency that appears in both the underdamped and overdamped forms of the spectral density. The entropy-production decomposition is careful, and the convention-dependence of the active-force time reversal is acknowledged. This is a useful reference calculation for chiral active solids.\n\nThe soft spot is load-bearing. Equation (26), from which the angular momentum M and the nonmonotonic curve Eq. (27) follow, is derived from the overdamped Langevin equation (77) in Appendix C. The overdamped condition is γ²/4 ≫ ω_q² for every contributing mode. For the representative parameters of Fig. 4 (K/(mγ²)=10³, τγ=1), the integrand of (26) peaks at qσ ≈ 0.026, where ω_q ≈ 0.58γ. That is not ≪ γ/2; the dominant modes are underdamped. The paper's Sec. 4 says chirality effects are similar in both regimes, but that statement concerns the presence of the non-dispersive peak, not the integrated angular momentum. The exact underdamped spectral density (11) has phonon poles at ±ω_q, which can change the q-integration and the shape of M(Ω). Without an underdamped calculation, an appropriate parameter scan, or simulations, the central prediction is unsupported for the model as stated.\n\nMinor issues: \"spontaneous\" is the wrong word—handedness is an input, so \"chirality-induced\" would be accurate. The harmonic defect-free lattice idealization limits quantitative relevance, but that is a standard trade-off for an exact solution.\n\nWho it's for: active matter theorists working on chiral solids, odd elasticity, and entropy production. It deserves a serious referee: the machinery is substantial and the new observable is worth vetting. I would send it out, but with a request to fix the overdamped inconsistency, either by performing the full underdamped integration or by restricting all claims and plots to genuinely overdamped parameters and checking a few cases numerically.","headline":"Clean exact solution, new angular-momentum prediction, but the headline curve rests on an overdamped reduction that the paper's own parameters violate.","tokens_in":31685,"tokens_out":4522,"would_cite":true,"duration_ms":39321,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A two-dimensional crystal of chiral active particles spontaneously acquires net angular momentum.","keywords":["chiral active matter","active crystals","angular momentum","harmonic lattice","active Ornstein-Uhlenbeck process","non-equilibrium steady state","entropy production","spatial velocity correlations"],"falsifier":"Simulate or construct a two-dimensional chiral active crystal and measure the per-particle angular momentum $M$ and the displacement power spectrum as the chirality rate $\\Omega$ is varied: the paper predicts $M\\neq 0$ for every finite $\\Omega$, $M\\to 0$ as $\\Omega\\to 0$ and $\\Omega\\to\\infty$, and a peak in the spectrum at frequency $\\Omega$ whose position does not change with wavevector, so finding $M=0$ at finite $\\Omega$ or no non-dispersive peak would settle the claim against it.","tokens_in":30725,"feed_emoji":"🌀","tokens_out":8476,"duration_ms":75762,"temperature":0.7,"pith_summary":"Marconi and Caprini study a two-dimensional triangular crystal in which every particle is chiral active: it self-propels, and its propulsion direction precesses at a fixed angular rate Ω. Working with linear harmonic springs between neighbors, they solve the coupled Langevin equations exactly by Fourier modes. They show that chirality generates a nonzero total angular momentum per particle whenever Ω≠0, with a magnitude that rises at small chirality, peaks, and falls back to zero at large chirality. The same mechanism puts a non-dispersive peak at frequency Ω in the displacement spectrum, producing oscillations in mean-square displacement and autocorrelations, and adds a chirality-dependent term to the entropy production rate. If correct, a chiral active crystal is a bulk rotating steady state whose handedness can be read from its dynamics without any external bias.","feed_headline":"Chiral active crystals spontaneously rotate","feed_subtitle":"Net angular momentum emerges from the chiral forces alone, peaking at intermediate chirality and vanishing in both extremes.","key_machinery":"The engine is the chiral active Ornstein-Uhlenbeck force, a self-propulsion whose direction both randomizes with persistence time $\\tau$ and precesses at angular rate $\\Omega$, together with the linearized harmonic lattice that makes the dynamics exactly solvable: Fourier transforming the coupled Langevin equations turns the crystal into independent damped oscillators driven by correlated chiral noise. The load-bearing object is the spectral angular momentum density Eq. (11), an even function of frequency whose peak at the non-dispersive frequency $\\Omega$ is produced by the chiral correlations of the active force; integrating it over the Brillouin zone and the spectrum yields the closed-form arctangent formula for the total angular momentum per particle, Eq. (27).","core_discovery":"The central claim is that the particles' chirality alone transforms a passive harmonic crystal into a rotating steady state. Specifically, the paper derives Eq. (27): the average angular momentum per particle $M \\approx \\frac{v_c}{2\\pi}\\frac{m\\gamma v_0^2}{c^2}\\arctan\\!\\left[\\frac{\\tau c^2 q_D^2/\\gamma}{|\\Omega|\\tau}\\frac{1}{1+(1+\\tau c^2 q_D^2/\\gamma)/(\\Omega^2\\tau^2)}\\right]$, which is nonzero for every finite chirality $\\Omega$, vanishes for $\\Omega=0$ and in the limit $\\Omega\\to\\infty$, and therefore has a maximum at an intermediate chirality. The angular momentum is entirely active in origin: it is independent of temperature, comes from the torque the chiral active force exerts on the lattice, and satisfies the torque balance $|\\mathcal{T}|=\\gamma M$ with the frictional torque. The same calculation shows that equal-time cross-correlations between $x$ and $y$ displacement and velocity components are nonzero and odd, that the displacement spectrum acquires a non-dispersive peak at the chirality frequency, and that the velocity correlation length shrinks as $1/\\sqrt{1+\\Omega^2\\tau^2}$, matching earlier continuum results.","pith_inferences":["A testable extension the paper does not pursue: in a dense monolayer of chiral colloids linked by soft springs, Eq. (27) implies that the measured mean circulation $\\langle r\\times v\\rangle$ should follow the same non-monotonic curve as $\\Omega$ is tuned by magnetic field or particle shape.","Because the $\\Omega$-peak is non-dispersive, the displacement spectrum offers a chirality fingerprint that could survive even where the net angular momentum is masked by boundaries or defects.","The torque balance $|\\mathcal{T}|=\\gamma M$ suggests that the microscopic chiral torque is directly readable from kinematic measurements, without force probes.","Coarse-graining the Gaussian correlations should yield an antisymmetric, odd elastic modulus, connecting this microscopic model to macroscopic chiral elastodynamics; the paper only notes the possibility."],"forward_implications":["A harmonic chiral active crystal in steady state carries a net angular momentum per particle whenever $\\Omega\\neq 0$, even though no external torque is applied.","The angular momentum is non-monotonic in chirality: it grows linearly for small $\\Omega\\tau$, peaks at an intermediate chirality, and decays as $1/(\\Omega\\tau)$ for large $\\Omega\\tau$.","The displacement spectrum acquires a non-dispersive peak at the chirality frequency $\\Omega$, present in both underdamped and overdamped regimes, which generates damped oscillations in mean-square displacement and in two-time displacement and velocity correlations.","Chirality shortens the spatial velocity correlation length according to $\\xi^2=\\xi^2(\\Omega=0)/(1+\\Omega^2\\tau^2)$, reproducing the prediction of the continuum theory.","The steady-state entropy production gains a rotational contribution tied to the angular momentum, so measuring dissipation can reveal the chiral torque."],"supporting_citations":[{"why":"Supplies the chiral active Ornstein-Uhlenbeck force dynamics with the angular-drift term $\\Omega\\times f^a$ that the model builds on.","marker":"48"},{"why":"Provides the harmonic-crystal Fourier-mode method for non-chiral active particles that this paper extends to chirality.","marker":"82"},{"why":"Continuum theory of chiral active crystals whose velocity-correlation length decreasing with chirality the paper's particle model reproduces.","marker":"76"},{"why":"Single chiral active particle in a harmonic trap; source of the fluctuation-reduction behavior used to interpret the non-monotonic angular momentum.","marker":"47"},{"why":"Shows the odd antisymmetric correlations of chiral active particles (odd diffusivity) that underlie the cross-correlation producing angular momentum.","marker":"25"},{"why":"Establishes spontaneous velocity alignment and the Ornstein-Zernike spatial velocity correlations in active crystals, the baseline that Eq. (20) extends.","marker":"15"},{"why":"Gives the non-chiral Ornstein-Zernike velocity-correlation profile and correlation length used as the $\\Omega=0$ limit of Eq. (21).","marker":"80"},{"why":"Provides the entropy-production calculation for non-chiral active crystals to which Eq. (33) reduces when $\\Omega\\to 0$.","marker":"87"},{"why":"Derives the analytical correlation length for active crystals from model parameters, referenced as the source of that result.","marker":"17"}],"fun_headline_variants":["Chirality alone spins active crystals","No torque needed: chiral crystals rotate","Chirality generates angular momentum in crystals","Active crystals rotate solely from chirality","Angular momentum emerges from chiral activity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analytical solution assumes an infinite, defect-free crystal with purely linear nearest-neighbor springs, so all predictions are exact only for Gaussian vibrations about an ideal triangular lattice; in a real active crystal with anharmonicities, vacancies, or plastic rearrangements, the predicted net angular momentum and non-dispersive peak could be weakened or destroyed.","fun_headline_variants_meta":{"raw":{"variants":["Chirality alone spins active crystals","No torque needed: chiral crystals rotate","Chirality generates angular momentum in crystals","Active crystals rotate solely from chirality","Angular momentum emerges from chiral activity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000894,"raw_usage":{"total_tokens":3867,"prompt_tokens":975,"completion_tokens":2892,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":591,"completion_tokens_details":{"reasoning_tokens":2830}},"tokens_in":591,"tokens_out":2892,"duration_ms":20338,"temperature":1.0,"reasoning_tokens":2830,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T04:50:40.264338+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate or construct a two-dimensional chiral active crystal and measure the per-particle angular momentum $M$ and the displacement power spectrum as the chirality rate $\\Omega$ is varied: the paper predicts $M\\neq 0$ for every finite $\\Omega$, $M\\to 0$ as $\\Omega\\to 0$ and $\\Omega\\to\\infty$, and a peak in the spectrum at frequency $\\Omega$ whose position does not change with wavevector, so finding $M=0$ at finite $\\Omega$ or no non-dispersive peak would settle the claim against it.","supporting_citations":[],"review_version":1}