{"id":"598a0599-1542-47cd-b0cf-f89fd52cb3c4","arxiv_id":"1908.03467","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A two-parameter formula with anisotropic tidal terms reproduces simulated halo spin-spin correlations out to 20 Mpc/h on dwarf to Milky Way mass scales.","lead":"This paper measures how the spin directions of neighboring dark matter halos align at different separations, using two N-body simulations, and proposes a new formula for the spin-spin correlation function. The formula, with two fitted parameters, matches the simulated correlations out to 20 Mpc/h, including a large-scale tail that earlier models missed, and the author suggests it could help probe dark matter.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (16) is not a derived two-point tidal statistic: in a statistically isotropic universe the J3,J5 terms cancel in ⟨T_ij T'_ij⟩, so g3,g5 are empirical fit parameters and the anisotropic-tidal origin of the large-scale tail is not established.","rationale":"The reader's weakest assumption is essentially correct, and the sharper version is that the assumed non-vanishing of J3,J5 in the two-point contraction is inconsistent with statistical isotropy of the cosmological ensemble. This does not invalidate the paper's empirical fit, but it removes the physical interpretation that anisotropic tides, rather than non-Gaussianity or halo selection, produce the large-scale tail. The Section 2.2 validation of Eq. (5) with no fitting parameters is genuine supporting evidence, but it constrains the conditional spin alignment model, not the new η(r) formula. The proposed check would settle the issue by measuring the tidal two-point function directly. Because the empirical formula may still be useful and the paper is explicit about the fitting parameters, the verdict remains CONDITIONAL rather than REJECT.","tokens_in":33,"tokens_out":17383,"duration_ms":361744,"concrete_test":"From the same ν2GC-H2 and SMDPL outputs, compute C(r)=⟨\\tilde T_ij \\tilde T'_ij⟩/⟨|\\tilde T'|^2⟩ directly from the reconstructed tidal tensors. Fit C(r) to \\tilde ξ(r)+g3 \\tilde J3(r)−g5 \\tilde J5(r) and compare the fitted coefficients with Table 1. If the directly fitted coefficients are consistent with zero (so C(r)= \\tilde ξ(r)) or differ significantly from Table 1, then the J3,J5 terms are not present in the tidal two-point correlation and Eq. (16) is not a tidal prediction. As a second check, fit g3,g5 to η(r) only for r<10 h^-1 Mpc and test whether the tail at r≥10 h^-1 Mpc is predicted without re-fitting.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is Section 3.2's assertion that for 'anisotropic T' the J3 and J5 terms survive the contraction used to write Eq. (16). Equation (14) is the general two-point function of a traceless tensor field in a statistically homogeneous and isotropic universe; Gaussianity is not required for that form. Contracting i=k, j=l makes the J3 and J5 terms cancel algebraically, leaving exactly ξ(r), as the paper itself notes. That individual tidal tensors are anisotropic (the cosmic web) is a property of a realization; the ensemble average over a statistically isotropic cosmological model is still isotropic, so the cancellation cannot be evaded by citing local filamentary structure. Therefore the g3 J3 and g5 J5 terms in Eq. (16) are an empirical basis extension, not a consequence of the anisotropic tidal effect at the two-point level. Since g3 and g5 are fitted per mass/redshift sample and dt is measured from the same halos, the excellent agreement with η(r) does not establish that anisotropic tides generate the r≥10 h^-1 Mpc tail; it is equally consistent with curve fitting or with non-Gaussianity. A further red flag is that printed Eq. (11) has η→−1/3 as r→∞ (because ⟨T_ij T'_ij⟩→0), so the derivation from Eq. (3) to Eq. (11) is not transparent; Eq. (16) is an ansatz rather than a derived consequence.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper investigates the halo spin-spin correlation function η(r) for low-mass dark matter halos using the ν2GC-H2 and SMDPL N-body simulations. It finds that η(r) decreases with separation much more slowly than the quadratic model η∝ξ² and retains a statistically significant tail at r≥10 h−1 Mpc on dwarf-galaxy scales. Building on the extended tidal-torque model of Lee (2019), which adds a linear tidal term controlled by parameter dt to the conditional spin covariance, the paper proposes Eq. (16), expressing η(r) in terms of the rescaled density correlation ξ̃(r) and the integral quantities J̃3, J̃5, with two adjustable parameters g3 and g5. The formula is fitted to the simulated η(r) for several mass ranges and redshifts and is reported to agree excellently with the measurements, especially at large separations. The paper interprets this agreement as evidence that anisotropic tidal fields produce the large-scale tail and discusses the possibility of using dwarf-halo spin-spin correlations as a complementary probe of dark matter.","tokens_in":13327,"tokens_out":16488,"duration_ms":162654,"significance":"If Eq. (16) were a derived consequence of the extended tidal-torque model, the result would be significant: it would provide a physical explanation for the large-scale spin-spin tail and a new connection between spin alignments and the linear density correlation function. The paper has genuine strengths: the validation of Eq. (5) on dwarf-galaxy scales in Section 2.2 is parameter-free and tests the extended model against a high-resolution simulation; the numerical measurements cover a broad mass range and several redshifts; and the empirical formula (16) does appear to describe the simulated η(r) well in Figures 3–8. However, the physical interpretation is not currently established. Equation (16) is introduced as a fitting formula with two free parameters and rests on an explicit assumption about anisotropic tidal tensors that is not derived. The fit is performed on the same data that supply ξ̃, J̃3, J̃5, and dt, and no out-of-sample prediction or goodness-of-fit statistic is provided.","major_comments":[{"comment":"The central claim that anisotropic tidal effects generate the r≥10 h−1 Mpc tail rests on the assertion that for anisotropic T the J3 and J5 terms survive the contraction ⟨T̃_ij T̃'_ij⟩. This is an assumption, not a derivation. In a statistically homogeneous and isotropic cosmological model, the ensemble-averaged two-point function of the tidal tensor retains the isotropic form given by Eq. (14), and the anisotropy of individual tidal tensors in a particular realization (the cosmic web) does not alter the symmetry of the ensemble average. The paper's own wording in Section 3.2 (\"Assuming here that for the anisotropic T, the terms containing J3 and J5 ... would not vanish, we propose the following fitting formula\") acknowledges this gap. Since g3 and g5 are fitted to each η(r) curve (Table 1) and the mean dt is measured from the same halos, the excellent agreement shown in Figures 3–8 is equally consistent with curve fitting and does not by itself establish that anisotropic tides cause the large-scale tail. A derivation of the contraction in a controlled anisotropic model, or an explicit demonstration that the J3/J5 terms appear with nonzero coefficients in the ensemble average, is required.","section":"Section 3.2, Eq. (16)"},{"comment":"The step from Eq. (3) to Eq. (11) is compressed and appears to contain an incorrect constant. For r→∞, the two-point function ⟨T̃_ij T̃'_ij⟩ of unit traceless tensors vanishes, so Eq. (11) gives η→−1/3, whereas the definition Eq. (6) requires η→0. The expansion of Eq. (3) into Eq. (11) also omits terms involving the product of the T̃² and T̃ terms; the derivation should be written out in full and the constant offset checked. This issue is load-bearing because Figure 2 uses Eq. (11) to validate the model before Eq. (16) is introduced.","section":"Section 3.2, Eqs. (11)–(12)"},{"comment":"The best-fit parameters g3 and g5 are tabulated without uncertainties or a goodness-of-fit statistic such as reduced χ², and the claim of \"excellent agreement\" is assessed visually. Moreover, the fit is performed on the same data that supply ξ̃, J̃3, J̃5, and the mean dt; no out-of-sample prediction (for example, a mass or redshift split not used in the fit) is presented. Consequently, the paper does not currently demonstrate that Eq. (16) has predictive content beyond being an interpolation formula for the simulated η(r).","section":"Section 3.3, Table 1"}],"minor_comments":[{"comment":"The phrase \"between the neighbor halos\" should be \"between neighboring halos\" or \"between neighbor halos\".","section":"Abstract"},{"comment":"\"snapshopts\" is a typo for \"snapshots\".","section":"Section 2.2"},{"comment":"The differentials \"dr'^2\" and \"dr'^4\" are likely typos; they should presumably be r'^2 dr' and r'^4 dr'.","section":"Eq. (15)"},{"comment":"The text states that the range 0≤r/(h−1 Mpc)≤20 is divided into bins of length Δr, but the value of Δr is never specified.","section":"Section 3.3"},{"comment":"The notation for the mass ranges is confusing; explicit ranges such as \"0.01≤M<0.5\" would be clearer than brackets followed by a comma.","section":"Table 1"},{"comment":"The figures would be easier to interpret if the numerical η(r) points included error bars, since the text says the one-standard-deviation errors are calculated.","section":"Figures 3–8"}],"recommendation":"major_revision","confidential_remarks":"The main technical gap is between the assumed anisotropic-tidal mechanism and the fitted formula: Eq. (16) is presented as a consequence of the extended model but is in fact an ansatz with two free parameters, and the physical interpretation outruns the derivation. The parameter-free validation in Section 2.2 shows that the author can perform careful numerical work, so the core issue is fixable in revision by either deriving the J3/J5 terms from the two-point tidal statistics or by substantially softening the causal claim and presenting Eq. (16) as an empirical fitting formula. I would also encourage the author to add out-of-sample tests and error bars on the fitted parameters before the paper can support the dark-matter-probe discussion in Section 4."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, the dwarf-scale measurements of η(r) out to r ≥ 10 h⁻¹ Mpc are new and worth having; the paper shows convincingly that neither the quadratic nor the linear ξ scaling describes that tail. Second, the new fitting formula Eq. (16) fits the simulations well, but I don't think it establishes the paper's causal claim. The J3 and J5 terms in the contracted two-point function ⟨T_ij T'_ij⟩ cancel in a statistically isotropic ensemble; local filamentary structure doesn't rescue that because the ensemble average is still isotropic. The authors acknowledge this and simply assume the terms survive for anisotropic T, so including them is an empirical basis extension, not a derivation.\n\nWhat the paper does well: Section 2.2 extends the Lee (2019) alignment model to dwarf scales with no fitting, and that part looks solid. The comparison of Eqs. (9), (10), and (16) across mass and redshift samples is thorough, and Table 1 gives a clear picture. The discussion of observational issues—baryonic spins, redshift-space distortions, MaNGA kinematics—is honest and useful.\n\nWhere it's soft: Eq. (16) has three free parameters per sample—g3, g5, and the mean dt—all determined from the same data being fit, so the excellent agreement is not a strong test. The step from Eq. (3) to Eq. (11) is compressed, and Eq. (11) as written has η → −1/3 as r → ∞, which is inconsistent with the definition of η; some of that may be absorbed into the later normalization, but it's not shown. The paper drops the -1/3 offset in Eq. (16) without comment. For the causal claim, an out-of-sample test or an independent derivation of g3 and g5 is needed; without it, the large-scale tail is just as consistent with non-Gaussianity or with curve fitting.\n\nThat said, the empirical part stands. This paper is for people working on intrinsic alignments and spin correlations who want the cleanest current measurement of dwarf-scale behavior. It deserves a serious referee, with a request for clarification of Eq. (11) and an explicit toning down of the anisotropic-tidal interpretation—or better, a prediction on an independent sample.","headline":"New dwarf-scale spin-correlation tail is worth having, but the anisotropic-tidal origin claim rests on an unproven assumption about ensemble averages.","tokens_in":13833,"tokens_out":3966,"would_cite":false,"duration_ms":42001,"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 claims that anisotropic tidal fields, not just non-Gaussianity, produce the large-scale halo spin-spin correlation tail.","keywords":["halo spin-spin correlation","tidal torque theory","anisotropic tidal field","intrinsic alignments","dwarf galaxies","dark matter halos","N-body simulation","large-scale structure"],"falsifier":"Measure $\\langle \\tilde T_{ij}(x)\\tilde T_{ij}(x+r)\\rangle$ directly from the same simulated tidal-field reconstructions and check whether, after subtracting $\\tilde{\\xi}(r)$, the residual is proportional to $\\tilde J_3(r)$ and $\\tilde J_5(r)$ with constant coefficients across all $r$ bins. A second decisive test is to fit Eq. (16) using only $r<10\\,h^{-1}$Mpc and see whether the extrapolation reproduces the $r\\ge10\\,h^{-1}$Mpc tail without refitting.","tokens_in":12747,"feed_emoji":"🌌","tokens_out":10429,"duration_ms":87483,"temperature":0.7,"pith_summary":"The paper argues that the long-range halo spin-spin correlation signal, which earlier tidal-torque models could not explain, is produced by the anisotropy of the tidal field rather than by non-Gaussianity alone. It proposes a new formula for $\\eta(r)$, the spin-spin correlation function, built from the rescaled linear density correlation $\\tilde{\\xi}(r)$ and two of its integrals, $\\tilde J_3(r)$ and $\\tilde J_5(r)$, multiplied by two fitted coefficients. With those coefficients, the formula reproduces N-body measurements of spin correlations of low-mass and dwarf halos across the mass range $0.05 \\le M/(10^{11}\\,h^{-1}M_\\odot) \\le 50$ at redshifts $z=0$, $0.2$, and $0.4$, including the tail at $r \\ge 10\\,h^{-1}$Mpc where earlier models failed. If the formula is correct, the large-scale spin alignment of dwarf halos becomes a practical complement to other probes of dark matter and of the linear density field.","feed_headline":"Tidal anisotropy explains long-range halo spin correlations","feed_subtitle":"Two-parameter formula ties dwarf-halo spin alignments to integrals of the density correlation, matching simulations.","key_machinery":"The load-bearing object is the extended conditional covariance for tidally induced spin alignments, $\\langle \\hat s_i\\hat s_j|\\hat T\\rangle = (\\frac13+\\frac35 c_t+\\frac35 d_t)\\delta_{ij}-\\frac35 c_t \\hat T_{ik}\\hat T_{kj}-\\frac35 d_t \\hat T_{ij}$, where the linear-in-$\\hat T$ term with the second spin-correlation parameter $d_t$ encodes the alignment of low-mass halo spins with the third eigenvector of the tidal field. Substituting this covariance into the definition $\\eta(r)=\\langle|\\hat s(x)\\cdot\\hat s'(x')|^2\\rangle-\\frac13$ and applying the same Wick-contraction approximations as earlier tidal-torque work yields Eq. (16). The integral functions $\\tilde J_3(r)$ and $\\tilde J_5(r)$ are rescaled versions of $J_3(r)=3r^{-3}\\int_0^r\\xi(r')r'^2\\,dr'$ and $J_5(r)=5r^{-5}\\int_0^r\\xi(r')r'^4\\,dr'$; the entire claim is that these terms survive when the tidal field is anisotropic and that they carry the large-scale tail.","core_discovery":"The central claim is that $\\eta(r)$ decays far more slowly with separation than the old $\\eta \\propto \\xi^2$ prediction, and that its statistically significant tail at $r \\ge 10\\,h^{-1}$Mpc cannot be explained by the linear-scaling $\\eta \\propto \\xi$ model either. The reason, the paper argues, is the isotropy assumption hidden in earlier derivations: when the tidal field is allowed to be anisotropic, the two-point correlation of traceless tidal tensors retains extra terms involving $J_3(r)$ and $J_5(r)$. This leads to Eq. (16), $\\eta(r) \\approx (18/25)d_t^2\\big[\\tilde{\\xi}(r)+g_3\\tilde J_3(r)-g_5\\tilde J_5(r)\\big]$, where $d_t$ is the linear tidal coupling of the spin alignment and $g_3,g_5$ are fitted coefficients absorbing the degree of anisotropy. Fitting those two coefficients to the $\\nu^2$GC-H2 and SMDPL simulations gives close agreement over the whole range of $r$ probed, with the $g_5$ term dropping sharply as mass and redshift increase, marking it as the most sensitive indicator of the anisotropic tidal effect.","pith_inferences":["The same $\\tilde J_3,\\tilde J_5$ decomposition could be tried on galaxy shape-shape (intrinsic alignment) correlations, where anisotropic tides may leave a comparable large-scale signature.","A sharper diagnostic would be to split the fitted coefficients by cosmic-web environment (filament, wall, void) and check whether the anisotropic term grows where filaments dominate.","The redshift-space caveat in the paper points to a concrete extension: convolve Eq. (16) with velocity-dispersion kernels and test whether the large-scale tail survives in redshift-space galaxy samples.","The severe drop of $g_5$ with mass and redshift suggests $g_5$ itself could be used as a measurable index of tidal anisotropy, if it can be calibrated against cosmic-web tracers."],"forward_implications":["If the formula holds, the large-scale tail of $\\eta(r)$ is a direct signature of anisotropic tides, so measuring dwarf-halo spin alignments at $r\\ge10\\,h^{-1}$Mpc can constrain tidal anisotropy rather than only non-Gaussianity.","Because $\\eta(r)$ is expressed through integrals of the linear density correlation $\\xi(r)$, a measured spin-spin correlation could in principle be inverted to reconstruct features of $\\tilde{\\xi}(r)$ without higher-order statistics.","The rescaled form of Eq. (16) is independent of the amplitude of $\\xi(r)$, giving it the potential to break degeneracies between the power-spectrum amplitude and the dark-matter abundance.","The model predicts that the anisotropic term weakens with increasing halo mass and redshift, so the linear-scaling approximation should become accurate for massive halos and at $z\\gtrsim0.4$."],"supporting_citations":[{"why":"supplies the original tidal-torque formula $\\eta\\propto\\xi^2$ that the new model must improve on.","marker":"Pen et al. 2000"},{"why":"establishes the Gaussian-isotropic Wick contraction and the parameter relation $a=3c_t/5$ used in the derivation.","marker":"Lee & Pen 2001"},{"why":"introduces the rival non-Gaussian explanation $\\eta\\propto\\xi$ that the paper argues is insufficient for the large-scale tail.","marker":"Hui & Zhang 2002"},{"why":"provides the linear-plus-quadratic fitting formula and reports the unexplained $r\\ge10\\,h^{-1}$Mpc signal.","marker":"Lee & Pen 2008"},{"why":"supplies the extended spin-alignment model with the second parameter $d_t$ and the linear-in-$\\hat T$ term.","marker":"Lee 2019"},{"why":"provides the high-resolution $\\nu^2$GC-H2 simulation dataset used for the dwarf-halo spin correlations.","marker":"Ishiyama et al. 2015"},{"why":"provides the SMDPL simulation used to test the formula for higher-mass halos and larger smoothing scales.","marker":"Klypin et al. 2016"}],"fun_headline_variants":["Anisotropic tides stretch halo spin correlations to 10 Mpc","Long-range halo spin alignments hint at anisotropic tides","New formula ties dwarf-halo spin alignments to density integrals","Two-parameter formula matches halo spin correlations across scales","Anisotropic tidal fields rework halo spin correlation law"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument stands on the assumption that the anisotropy of the tidal field changes the tensor two-point correlation exactly in the form of two extra terms proportional to $J_3$ and $J_5$ with coefficients that do not depend on separation; if that shape is wrong, the agreement with simulations is just curve fitting.","fun_headline_variants_meta":{"raw":{"variants":["Anisotropic tides stretch halo spin correlations to 10 Mpc","Long-range halo spin alignments hint at anisotropic tides","New formula ties dwarf-halo spin alignments to density integrals","Two-parameter formula matches halo spin correlations across scales","Anisotropic tidal fields rework halo spin correlation law"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000726,"raw_usage":{"total_tokens":3324,"prompt_tokens":1086,"completion_tokens":2238,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":702,"completion_tokens_details":{"reasoning_tokens":2156}},"tokens_in":702,"tokens_out":2238,"duration_ms":16589,"temperature":1.0,"reasoning_tokens":2156,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:12:50.230072+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure $\\langle \\tilde T_{ij}(x)\\tilde T_{ij}(x+r)\\rangle$ directly from the same simulated tidal-field reconstructions and check whether, after subtracting $\\tilde{\\xi}(r)$, the residual is proportional to $\\tilde J_3(r)$ and $\\tilde J_5(r)$ with constant coefficients across all $r$ bins. A second decisive test is to fit Eq. (16) using only $r<10\\,h^{-1}$Mpc and see whether the extrapolation reproduces the $r\\ge10\\,h^{-1}$Mpc tail without refitting.","supporting_citations":[{"cited_title":"2000, ApJ, 543, L107 Planck Collaboration, Ade, P","cited_arxiv_id":null,"evidence_quote":"supplies the original tidal-torque formula $\\eta\\propto\\xi^2$ that the new model must improve on."},{"cited_title":"2001, ApJ, 555, 106","cited_arxiv_id":null,"evidence_quote":"establishes the Gaussian-isotropic Wick contraction and the parameter relation $a=3c_t/5$ used in the derivation."},{"cited_title":"Intrinsic/Extrinsic Density-Ellipticity Correlations and Galaxy-Galaxy Lensing","cited_arxiv_id":"astro-ph/0205512","evidence_quote":"introduces the rival non-Gaussian explanation $\\eta\\propto\\xi$ that the paper argues is insufficient for the large-scale tail."},{"cited_title":"2008, ApJ, 681, 798","cited_arxiv_id":null,"evidence_quote":"provides the linear-plus-quadratic fitting formula and reports the unexplained $r\\ge10\\,h^{-1}$Mpc signal."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the high-resolution $\\nu^2$GC-H2 simulation dataset used for the dwarf-halo spin correlations."},{"cited_title":"2016, MNRAS, 457, 4340","cited_arxiv_id":null,"evidence_quote":"provides the SMDPL simulation used to test the formula for higher-mass halos and larger smoothing scales."}],"review_version":1}