{"id":"cb1122de-5230-4d95-ae5d-5c0a9d3b1737","arxiv_id":"2504.20133","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In FIRE-2 cosmological simulations, the quadrupole (elongation) of Milky Way-mass dark matter halos stays aligned with cosmic filaments for gigayears, and a Large Magellanic Cloud-mass satellite produces a tenfold stronger quadrupole response than in spherical halo models.","lead":"This paper uses cosmological simulations of Milky Way-like galaxies to show that the dark matter halo's elongated shape persistently points along cosmic filaments, while satellite encounters add temporary disturbances. A new spectral decomposition separates these two effects and finds that the Large Magellanic Cloud's gravitational wake is about ten times stronger in a realistic, asymmetric halo than in simplified spherical models.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'order-of-magnitude' LMC quadrupole claim rests on an uncontrolled comparison of m12b with a spherical-halo simulation; a matched shape-control run is needed.","rationale":"The paper's qualitative core—dipole transient, quadrupole persistent, and filament alignment retained even at ~30 kpc—is supported by consistent behavior across seven FIRE-2 halos (Figs. 8-9) and by prior literature (e.g., Vera-Ciro et al. 2011; Bailin & Steinmetz 2005; Aubert et al. 2004). I do not see an internal inconsistency in the BFE/mSSA machinery; the heuristic PC grouping in Appendix C is a judgment call but is documented and the decomposition is internally coherent. The load-bearing weak point is exactly the quantitative factor-of-ten claim singled out by the reader. I agree with the reader's weakest assumption: comparing m12b with Garavito-Camargo et al. (2021) conflates initial halo shape with host mass, satellite mass and orbit, baryonic physics, and resolution. The paper frames this as evidence that pre-existing triaxiality boosts the satellite response, but that causal attribution requires a matched control. Without it, the abstract's 'order of magnitude larger' statement is not yet established. This does not undermine the filament-memory result, so the appropriate verdict remains CONDITIONAL, matching the reader's report.","tokens_in":32032,"tokens_out":6165,"duration_ms":71605,"concrete_test":"Run a matched shape-control experiment using the m12b snapshot just before the LMC-analog's first infall (e.g., ~6 Gyr): resimulate (A) the original triaxial halo plus the LMC-analog on its actual orbit, and (B) the same halo with its quadrupole artificially removed (e.g., re-sampling particles to a spherically symmetrized density profile with the same radial profile and velocity dispersion) plus the same LMC-analog orbit, at the same resolution and with the same baryonic content. Apply the same debris-removal, BFE, and mSSA pipeline to both runs and compare the LMC-group quadrupole power at the first pericenter. If the A/B ratio is not ~10, the enhancement is not specifically due to pre-existing halo shape; if it is ~10, the claim is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's most novel quantitative result is that an LMC-mass satellite induces a quadrupolar response an order of magnitude larger in cosmological FIRE-2 halos than in 'spherically symmetric MW models' (Sec. 5.1, Fig. 10; also Sec. 6). The evidence is one curve from m12b, selected as the best MW-LMC analog, plotted against the constrained MW-LMC simulation of Garavito-Camargo et al. (2021) after multiplying that curve by ten. The two simulations differ not only in initial halo shape but also in host mass and concentration, satellite mass and orbital parameters (m12b: mass ratio 1:8, pericenter ~40 kpc, eccentricity 0.83), baryonic treatment, numerical resolution, and in which multipoles and radial truncation enter the quoted power. Because no control run varies only the initial shape, the factor of ten cannot be attributed to pre-existing triaxiality as the text argues; it could be driven by, e.g., the orbit's deeper pericenter, a more concentrated host, baryonic contraction, or resolution-dependent wake growth. The qualitative claim that a triaxial host responds differently from a spherical one is plausible and independently motivated by prior work, but the quantitative factor of ten is not yet established from this comparison.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes seven FIRE-2 Latte zoom-in cosmological hydrodynamical Milky Way-mass halos using spherical-harmonic basis function expansions (BFE), focusing on the dipole and quadrupole gravitational power. It reports that the quadrupole orientation aligns with the dominant cosmic filament, that this alignment persists into the inner halo (~30 kpc), and that the quadrupole distortion can reach 1–2 times the spherical density near the virial radius. Using multichannel singular spectrum analysis (mSSA) on the m12b halo, the authors separate a slowly declining filament-associated quadrupole from an LMC-analog-associated response, and they argue that the latter is an order of magnitude stronger than in spherically symmetric MW–LMC constrained simulations. The paper then discusses implications for disk warps, the LMC-induced wake, and stellar-stream modeling.","tokens_in":32290,"tokens_out":7342,"duration_ms":74844,"significance":"If the quantitative claims hold, the paper makes a valuable contribution by connecting cosmic-web accretion to present-day Milky Way halo shape and by providing a method to separate filamentary and satellite-induced quadrupole distortions. The study's strengths include the use of public FIRE-2 simulations and open-source pyEXP/EXP software, explicit validation of the BFE density reconstruction (Fig. 2), tests of radial-basis stability (Appendix A), and an external spatial-alignment check of the mSSA decomposition via the satellite and filament position angles (Fig. 13). The qualitative picture—that triaxial cosmological halos respond more strongly to an LMC-mass satellite than spherical models—is plausible and well motivated by earlier work. The paper is clearly written and the methodology is mostly sound, but the headline order-of-magnitude comparison and the mSSA grouping procedure need additional support before the central claim can be considered established.","major_comments":[{"comment":"The central quantitative claim—that an LMC-mass satellite induces a quadrupolar response an order of magnitude larger in cosmological halos than in spherically symmetric MW models—rests on an uncontrolled comparison. The m12b halo differs from the Garavito-Camargo et al. (2021) constrained simulation not only in initial halo shape but also in host mass and concentration, satellite mass and orbital parameters (mass ratio 1:8, pericenter ~40 kpc, eccentricity 0.83), baryonic treatment, numerical resolution, and in which multipoles and radial truncation enter the quoted power. Moreover, the m12b analysis removes the satellite's debris (Sec. 2.3), whereas the comparison simulation may include different particle content. The factor of ten therefore cannot currently be attributed to pre-existing triaxiality as stated in the text. A matched control run that varies only the initial halo shape (e.g., sphericalizing the same host and rerunning the same LMC orbit), or a systematic comparison across several FIRE-2 halos with different triaxialities, is needed to isolate the shape effect. Without such a control, the order-of-magnitude statement should be softened to a qualitative enhancement.","section":"Sec. 5.1, Fig. 10"},{"comment":"The mSSA decomposition is the basis for separating filament and LMC quadrupole contributions, but the PC grouping (PC0–PC3 vs. PC4–PC11) is a judgment call based on qualitative inspection of time evolution, DFT spectra, and the w-correlation matrix. The paper does not report sensitivity to the mSSA window length L (fixed at half the series), the number of PCs retained, or the grouping rule. Because the reconstructed filament and LMC amplitudes, and hence the factor-of-ten comparison, depend on this grouping, the authors should provide robustness tests—for example, varying L, using an automated clustering of the w-correlation matrix, or showing that the LMC group continues to track the satellite position when the grouping is perturbed. Without such tests, the disentanglement claim is not uniquely determined by the data.","section":"Sec. 2.4, Appendix C"},{"comment":"The order-of-magnitude statement is demonstrated for a single halo, m12b, although the text says '(and in other FIRE-2 analogs)'. No other halo is shown in the mSSA comparison, and the selection of m12b as the MW–LMC analog is based on prior work rather than on an independent criterion established here. Given that the paper analyzes a sample of seven halos, the authors should either present the same decomposition for other halos hosting LMC-mass satellites or explicitly limit the quantitative claim to m12b and quantify the expected halo-to-halo scatter. This is particularly important because the paper draws general conclusions about the Milky Way's present-day response in Sec. 6.","section":"Sec. 6, Fig. 10"}],"minor_comments":[{"comment":"The bullet 'Haloresponseatpericenter. ,andthisbulkmotionmust...' is garbled and appears to be a duplicated or incomplete sentence; please fix the text.","section":"Sec. 6"},{"comment":"The text states that reconstructed profiles are 'accurate to within 5%' out to the virial radius, but the figure shows 8–10% deviations in the inner 2–3 kpc and beyond Rvir at early times; please restate the accuracy claim to match the figure.","section":"Sec. 2.3, Fig. 2"},{"comment":"The paper states in Sec. 5.1 that the first six PCs account for ~95% of power and the first eleven ~99%, while Appendix C says the first 12 PCs account for 98%; these numbers should be reconciled or more precisely defined.","section":"Sec. 5.1 and Appendix C"},{"comment":"The phrase 'quadrupole distortion can reach 1–2 times the spherical density' could be misread as a total density contrast; please define it explicitly as max|ρ_l=2|/ρ_l=0 at a given radius, as used in Sec. 4.2.","section":"Sec. 4.2, Abstract, Sec. 6"},{"comment":"There are several typos and grammatical slips, including 'filamentry memory' (Sec. 4.2 and Sec. 6), 'within within 5%' (Sec. 2.3), 'an shown in TNG' (Sec. 6), 'overalltheitremains' (Sec. 5.2), and 'color-coded curvesa' (Fig. 8 caption).","section":"Throughout"},{"comment":"The paper attributes the secular decline in quadrupole power to halo relaxation, but does not separate this from the growth of the monopole normalization over time; a brief discussion of this normalization effect would improve interpretation of Figs. 4 and 5.","section":"Sec. 3.3"}],"recommendation":"major_revision","confidential_remarks":"The paper is promising and the qualitative conclusions are likely correct, but the headline order-of-magnitude claim needs either a matched control simulation or a more qualified statement, and the mSSA grouping needs robustness checks. I would support publication after these issues are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The qualitative core of this paper is solid: across seven FIRE-2 halos, the quadrupole orientation tracks the dominant filament, persists at a fractional density of ~0.5 even at 30 kpc, and decays as the filamentary environment diffuses. That is consistent with earlier work but shown here more systematically and with a nice visualization framework. The BFE reconstruction is well-tested (within 5% out to the virial radius, with basis stability checked in Appendix A), and the mSSA decomposition into a slow filament mode (~4 Gyr) and a fast satellite mode (~1.5 Gyr) is a genuinely useful addition: the spatial reconstructions in Figs. 11-12 show the two components behaving differently, and the alignment diagnostics in Fig. 13 give the decomposition some independent support. The paper is clearly written and honestly flags its choices, including debris removal and the heuristic PC grouping.\n\nThe soft spots are concentrated in the quantitative headline. The claim that an LMC-mass satellite induces a quadrupole response an order of magnitude larger in cosmological halos than in spherical models rests entirely on comparing one curve from m12b to the Garavito-Camargo et al. (2021) constrained simulation, multiplied by ten. Those simulations differ in host mass, satellite mass, orbit, baryonic treatment, and resolution; no control varies only the initial halo shape. So the number is not established. The physical argument that a pre-existing quadrupole allows a stronger response is plausible, and it is supported by earlier idealized work, but this comparison alone cannot pin down the factor of ten. Also, the text says \"in m12b (and in other FIRE-2 analogs)\" but only m12b is shown; that claim is currently unquantified. The mSSA grouping is a judgment call; the Appendix C diagnostics are helpful, but there is no uncertainty estimate on the group assignment, so the separation should be treated as suggestive rather than exact.\n\nWho gets value from this? Anyone modeling the MW potential from stellar streams, disk warp, or LMC-wake observations, and anyone building time-dependent halo models. The qualitative result—that halos retain a filament-locked quadrupole deep into the inner regions—is important for that community. The factor-of-ten claim, if it survives a matched shape-control run, would be an important result; until then it should be reported as a single-halo case study. I would send this to review: the framework is reproducible, the BFE/mSSA tools are public, and the core qualitative science deserves referee time. The referee should ask for a matched control simulation or, failing that, for the paper to explicitly reframe the order-of-magnitude claim as a motivating hypothesis. That is a revision, not a rejection.","headline":"The filament-quadrupole memory result is solid and worth engaging, but the headline factor-of-ten LMC quadrupole claim rests on an uncontrolled comparison and should be treated as a hypothesis, not an established number.","tokens_in":32922,"tokens_out":2365,"would_cite":true,"duration_ms":26607,"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":"Cosmic filaments imprint a long-lived quadrupole shape on Milky Way-mass dark halos, and an LMC-like satellite stirs a far stronger response in those halos than spherical models predict.","keywords":["dark matter halo shape","halo triaxiality","quadrupole moment","cosmic filaments","Large Magellanic Cloud","spherical harmonic expansion","singular spectrum analysis","Milky Way dynamics"],"falsifier":"Run matched N-body simulations with identical host and satellite masses and orbital parameters, changing only the initial halo shape from spherical to a realistic triaxial configuration; the claim fails if the peak quadrupole power does not increase by roughly an order of magnitude. Observationally, the claim would also be tested by comparing the amplitude of the LMC-induced stellar wake in high-precision astrometric data with predictions from spherical versus triaxial time-evolving halo models.","tokens_in":31826,"feed_emoji":"🌌","tokens_out":8465,"duration_ms":75834,"temperature":0.7,"pith_summary":"The paper sets out to separate the two processes that sculpt the dark matter halos of Milky Way-mass galaxies: steady accretion of material along cosmic filaments and the sudden gravitational kicks from massive satellites such as the Large Magellanic Cloud. Using spherical-harmonic decompositions of zoomed cosmological hydrodynamical simulations, it shows that the dipole moment is a transient satellite-driven disturbance while the quadrupole moment is a persistent feature that mirrors the halo's triaxial shape. The claimed discovery is that the quadrupole's orientation stays aligned with the dominant feeding filament, imprinting a long-lived memory in the halo's morphology even in the inner ~30 kiloparsecs, where the quadrupole distortion reaches about half the spherical density. The paper further reports that an LMC-mass satellite produces a quadrupolar response roughly an order of magnitude larger in these realistic halos than in models that start from a spherical Milky Way halo, which matters because most dynamical models of the Milky Way assume spherical symmetry.","feed_headline":"Filaments leave a long-lived shape memory inside Milky Way halos","feed_subtitle":"An LMC-scale satellite drives a tenfold stronger quadrupole response in realistic triaxial halos than spherical models predict.","key_machinery":"The load-bearing tool is the basis function expansion: the halo density is written as $\\rho(\\mathbf{x},t)=\\sum_{\\ell m n} C_{\\ell m n}(t)\\,\\varrho_{\\ell n}(r)Y_\\ell^m(\\theta,\\phi)$, with an empirical radial basis fixed by the time-averaged spherical profile, and the gravitational power in harmonic $\\ell$ is $P_\\ell=\\sum_{n,m} C_{\\ell,n,m}^2$. The quadrupole position angle is read off the lowest-order $\\ell=2$ coefficients, and the orientation of the dominant filament is estimated from infalling particles outside the virial radius. Multichannel singular spectrum analysis on the coefficient time series separates slow (~4 Gyr) filament-driven modes from faster (~1.5 Gyr) satellite-driven modes, letting the paper reconstruct the density fields of each component separately.","core_discovery":"The paper decomposes the dark-matter density of Milky Way-mass halos into spherical harmonic basis functions and tracks the time-dependent gravitational power per harmonic. It finds that the dipole mode is a transient response to infalling satellites, whereas the quadrupole mode is a persistent, slowly declining feature whose orientation stays aligned with the dominant cosmic filament feeding the halo. The quadrupole's imprint reaches inward to roughly 30 kiloparsecs, and at the virial radius the quadrupole distortion can be one to two times the spherical density. Applying multichannel singular spectrum analysis to the quadrupole coefficients, the paper separates a slow, filament-driven component from a faster, satellite-driven component, and reports that an LMC-mass satellite induces a quadrupolar response about an order of magnitude larger in these cosmological halos than in constrained models that assume a spherical Milky Way halo.","pith_inferences":["If the tenfold enhancement holds for lower-mass satellites as well, even smaller infalling galaxies could leave detectable quadrupole signatures in a triaxial halo, a possibility the paper does not test.","If triaxiality boosts the LMC-induced quadrupole, then LMC mass estimates derived from wake measurements in spherical models may need rescaling, since the same wake amplitude would require a different satellite mass in a triaxial potential.","Extending the same decomposition to the octupole and higher harmonics could reveal whether residual angular structure at small radii is dominated by the live halo response or by stripped satellite debris, a distinction the paper leaves open.","A direct observational test would compare the inner halo's quadrupole orientation with the symmetry plane of the Magellanic Stream or the orbital plane of the Sagittarius stream, though the paper does not perform such a comparison."],"forward_implications":["Static, spherical Milky Way halo models will underestimate the gravitational response to the LMC, including the induced wake and the torque that drives disk warps.","The orientation of the inner ~30 kpc quadrupole can be read as a fossil record of the direction of the cosmic filament that fed the halo.","Because the satellite-induced quadrupole has stayed near its peak since the LMC's first pericenter, present-day Milky Way observables should be modeled with time-evolving asymmetric potentials.","The dipole must be measured and subtracted before the intrinsic triaxial shape can be inferred from stellar kinematics.","The mSSA decomposition gives a general method to compare how environment versus mergers set halo shape in different cosmological halos."],"supporting_citations":[{"why":"Provides the constrained Milky Way-LMC simulation with a spherical halo whose quadrupole power is the baseline for the factor-of-ten comparison.","marker":"N. Garavito-Camargo et al. 2021"},{"why":"Another spherical-halo Milky Way-LMC model used as a comparison for the satellite-induced quadrupole.","marker":"S. Lilleengen et al. 2023"},{"why":"Introduces the multichannel singular spectrum analysis technique used to separate filament and satellite components.","marker":"M. D. Weinberg & M. S. Petersen 2021"},{"why":"Supplies the basis function expansion formalism and code used to expand the halos.","marker":"M. S. Petersen et al. 2022b"},{"why":"Identifies the specific simulated halo used as the representative Milky Way-LMC analog, defining the case study.","marker":"A. Arora et al. 2024a"},{"why":"Describes the physics model used in the simulations that produce the analyzed halos.","marker":"P. F. Hopkins et al. 2018"},{"why":"Details the simulation suite analyzed here, including resolution, particle masses, and snapshot cadence.","marker":"A. Wetzel et al. 2023"},{"why":"Establishes the connection between filamentary accretion and halo triaxiality that this paper quantifies.","marker":"C. A. Vera-Ciro et al. 2011"},{"why":"Provides the empirical radial basis construction used to fit arbitrary spherical profiles.","marker":"M. D. Weinberg 1999"}],"fun_headline_variants":["Filaments leave a long-lived quadrupole memory in Milky Way halos","LMC-mass satellite triggers tenfold quadrupole response in realistic halos","Cosmic filaments mold Milky Way halo shape with a persistent quadrupole","Filament alignment wins over satellite kicks in shaping dark matter halos","Quadrupole imprint from filaments persists in Milky Way-scale halos"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The factor-of-ten comparison rests on two simulations that differ in host mass, satellite mass, orbit, resolution, and baryonic treatment; without a control run that varies only the halo's initial shape, part of the enhancement could come from those other differences.","fun_headline_variants_meta":{"raw":{"variants":["Filaments leave a long-lived quadrupole memory in Milky Way halos","LMC-mass satellite triggers tenfold quadrupole response in realistic halos","Cosmic filaments mold Milky Way halo shape with a persistent quadrupole","Filament alignment wins over satellite kicks in shaping dark matter halos","Quadrupole imprint from filaments persists in Milky Way-scale halos"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00104,"raw_usage":{"total_tokens":4424,"prompt_tokens":1041,"completion_tokens":3383,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":657,"completion_tokens_details":{"reasoning_tokens":3287}},"tokens_in":657,"tokens_out":3383,"duration_ms":21350,"temperature":1.0,"reasoning_tokens":3287,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:36:53.291980+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run matched N-body simulations with identical host and satellite masses and orbital parameters, changing only the initial halo shape from spherical to a realistic triaxial configuration; the claim fails if the peak quadrupole power does not increase by roughly an order of magnitude. Observationally, the claim would also be tested by comparing the amplitude of the LMC-induced stellar wake in high-precision astrometric data with predictions from spherical versus triaxial time-evolving halo models.","supporting_citations":[{"cited_title":"S., Erkal, D., et al","cited_arxiv_id":null,"evidence_quote":"Another spherical-halo Milky Way-LMC model used as a comparison for the satellite-induced quadrupole."},{"cited_title":"D., & Petersen, M","cited_arxiv_id":null,"evidence_quote":"Introduces the multichannel singular spectrum analysis technique used to separate filament and satellite components."},{"cited_title":"C., Sanderson, R","cited_arxiv_id":null,"evidence_quote":"Details the simulation suite analyzed here, including resolution, particle masses, and snapshot cadence."},{"cited_title":"A., Sales, L","cited_arxiv_id":null,"evidence_quote":"Establishes the connection between filamentary accretion and halo triaxiality that this paper quantifies."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the empirical radial basis construction used to fit arbitrary spherical profiles."}],"review_version":1}