{"id":"bdfc966d-c618-4ca2-93a4-0dc5f61b0b35","arxiv_id":"2507.08925","paper_version":1,"verdict":"UNVERDICTED","confidence":"HIGH","novelty_score":0.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"A comprehensive introductory review of cosmological galaxy simulation methods, covering initial conditions, numerical solvers, star formation and feedback, analysis, and validation.","lead":"This paper is a textbook-style review of how scientists build computer simulations of galaxy formation, from the Big Bang's leftover seeds to the galaxies we see today. It is a useful entry point for anyone who wants to know what modern galaxy simulations contain and how their results are tested against observations.","discovery_kind":"review","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: the review is accurate, well-hedged, and its declared exclusions are explicit rather than misleading.","rationale":"The paper is a review article, not a research contribution, so the Pith verdict semantics correctly mark it UNVERDICTED. The reader's weakest assumption—representativeness of the field—is indeed the only plausible point of failure for the central claim. I examined whether the taxonomy of numerical solvers, the catalog of sub-grid baryonic physics, the discussion of post-processing, and the validation strategies omit or misrepresent major approaches. The review includes the dominant modern codes and methods, names the main simulation suites (EAGLE, Illustris-TNG, Simba, FIRE/FIREbox, NewHorizon, Flamingo, etc.), and explicitly lists excluded physics. It also consistently hedges the status of sub-grid models as effective, observationally calibrated prescriptions, noting degeneracies and open questions. The historical and technical descriptions are consistent with standard practice as presented in other recent reviews. No load-bearing objection emerged; the only worthwhile check is a systematic coverage audit against the existing review literature, which would confirm or refute representativeness in a falsifiable way. The verdict should remain UNCHANGED.","tokens_in":30966,"tokens_out":7533,"duration_ms":87508,"concrete_test":"Cross-check the code and method lists in §2.3–§2.4 against independent recent reviews (e.g., Vogelsberger et al. 2020, Naab & Ostriker 2017, Crain & van de Voort 2023): if any widely-used cosmological galaxy simulation code or baryonic physics module is absent from both the main text and the cited literature, the representativeness assumption would need revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is pedagogical: the review must accurately represent the essential components and methods of cosmological galaxy simulations. That could fail only if the solver taxonomy or sub-grid catalog omitted or misrepresented a major approach. I checked the relevant sections: §2.3 covers tree, particle-mesh, Tree-PM, SPH, AMR, and moving-mesh/mesh-free methods, citing Gadget-4, Gasoline, ChaNGa, Swift, Arepo, Gizmo, Enzo, Ramses, and ART; §2.4 covers cooling, ISM pressure floors, star formation criteria, SN feedback variants (thermal, kinetic, decoupled winds, superbubbles), BH seeding, Bondi/torque accretion, and AGN feedback modes; §2.6 covers convergence, parameter variation, cross-simulation comparisons, and observational validation. Scope exclusions (magnetohydrodynamics, cosmic rays, radiation hydrodynamics, thermal conduction, viscosity) are stated explicitly in §2.4 and revisited in the outlook. The paper repeatedly acknowledges calibration degeneracies and unresolved coupling of feedback, so a newcomer would not be misled into overtrusting sub-grid models. No internal inconsistency or factual error that would undermine the review's stated purpose was found.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript is a review of cosmological galaxy simulations. It begins with a brief historical account (Section 2.1) from Holmberg's analog experiments to modern cosmological hydrodynamics, then describes the generation of CMB-constrained initial conditions (Section 2.2), the main numerical solvers for dark matter and gas (Section 2.3), sub-grid baryonic physics including cooling, star formation, stellar feedback, black hole seeding and AGN feedback (Section 2.4), standard post-processing analysis (Section 2.5), and strategies for validating simulations through convergence tests, parameter variations, cross-simulation comparisons, and observational benchmarks (Section 2.6). The final section discusses next-generation directions: higher resolution, additional physics, and machine-learning tools. The abstract's central assertion is that the review 'provides an introductory overview'; the manuscript delivers on this claim as a synthesis.","tokens_in":31040,"tokens_out":6494,"duration_ms":68415,"significance":"The review's value is pedagogical and organizational rather than containing new results. It is technically careful: it repeatedly flags uncertainties and unresolved issues such as the overcooling problem, the resolution dependence and calibration degeneracies of sub-grid models, and the limited understanding of feedback coupling. It explicitly lists scope exclusions (magnetohydrodynamics, cosmic rays, radiation hydrodynamics, thermal conduction, viscosity) and revisits them in the outlook. The bibliography is extensive and up-to-date (e.g., Flamingo 2025, EDGE-INFERNO 2025, FIREbox HR 2025, Rose et al. 2025), and the figures (Figs. 2-4) effectively illustrate key concepts. If the review is correct, it provides a trustworthy entry point for newcomers.","major_comments":[],"minor_comments":[{"comment":"The sentence 'Dark matter is a fundamental component of the universe, comprising approximately 85% of its total mass' is imprecise: dark matter constitutes roughly 85% of the matter content, but only about 26% of the total energy density of the universe. Please rephrase to avoid confusing 'matter' with 'mass-energy'.","section":"Sec. 2.3.1"},{"comment":"In the 'Star Formation and Evolution' paragraph, 'star formation generally only occurs when certain gas conditions as met' contains a typo ('as' should be 'are').","section":"Sec. 2.4"},{"comment":"The statement that AGN simulations 'have successfully reproduced the AGN luminosity function, which is dominated by black holes with masses around 10^8 M☉' is an overgeneralization; the luminosity function is a population-level statistic and is not simply dominated by a single black-hole mass. I recommend softening the claim or citing a specific study that makes this point.","section":"Sec. 2.4"},{"comment":"In the list of dust radiative transfer codes, 'and and POWDERDAY' contains a doubled conjunction; please delete one 'and'.","section":"Sec. 2.5"},{"comment":"Equation (1) writes the collisionless Boltzmann equation with df/dt = 0, but f is a function of (r, v, t); for clarity, I suggest writing the partial differential equation explicitly as ∂f/∂t + v·∂f/∂r − ∇φ·∂f/∂v = 0.","section":"Sec. 2.3.2"},{"comment":"The sentence beginning 'Addressing these discrepancies is one of the central goals of cross-simulation comparison projects such as Aquila (Scannapieco et al., 2012), AGORA (Kim et al., 2014; Roca-Fabrega et al., 2024).' would read more smoothly with a colon or semicolon after 'projects'.","section":"Sec. 2.6"}],"recommendation":"minor_revision","confidential_remarks":"I see no controversial claims or internal inconsistencies in this review. The editorial issues listed are straightforward to fix. The manuscript is well within the scope of a review journal and will be useful to the intended audience."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a review, and it does not pretend to be anything else. The reader's UNVERDICTED call is right. What it does well: the taxonomy of solvers (tree/PM/Tree-PM, SPH/AMR/moving-mesh), sub-grid models (cooling, star formation, stellar and AGN feedback), and validation approaches (convergence, parameter variation, cross-simulation comparison, mock observations) is accurate and carefully hedged. Scope exclusions - MHD, cosmic rays, radiation hydrodynamics, thermal conduction, viscosity - are stated up front in Section 2.4, so a newcomer would not be misled. The discussion of the overcooling problem and calibration degeneracies is honest, and the citations are broad and balanced; the self-citations lean toward FIRE-style work, but EAGLE, Illustris-TNG, Simba, and other major suites get fair coverage.\n\nThe soft spots are minor and mostly intrinsic to the genre. There is no critical comparison of evidence for different sub-grid choices; a reader gets the landscape but not the empirical basis for preferring one approach over another. The exclusions are sensible, but they are the physics many in the field would name as the next frontier, and the outlook gives them only a gesture. The ML/emulator section is superficial, though appropriate for an introductory review. One slightly bigger gap: the text says sub-grid models are calibrated to observations but never dwells on how degenerate that calibration is, beyond a passing mention. None of this undercuts the paper's stated purpose.\n\nIf this were submitted as a review article, I would send it to a serious referee and expect it to survive with minor comments. I would cite it as a methods reference in the next year, and it would be a fine reading-group orientation for students or non-specialists. Not a deep dive, but a trustworthy map.","headline":"A careful, well-hedged introductory review of galaxy simulation methods; no new science, but a solid orientation piece that deserves a real referee if submitted.","tokens_in":579,"tokens_out":845,"would_cite":true,"duration_ms":30418,"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":"This review claims that modern cosmological galaxy simulations can be understood as a single pipeline running from CMB-calibrated initial conditions, through gravity and hydrodynamics solvers, to sub-grid baryonic feedback and mock…","keywords":["cosmological simulations","galaxy formation","initial conditions","N-body methods","hydrodynamics","sub-grid physics","stellar feedback","AGN feedback"],"falsifier":"Run two hydrodynamics solvers with identical initial conditions, cooling tables, and feedback prescriptions, and check whether their predicted galaxy stellar mass functions by $z=0$ agree within observational uncertainties; a disagreement larger than those uncertainties would call into doubt the review's claim that the pipeline is stable and adequately validated.","tokens_in":30666,"feed_emoji":"🌌","tokens_out":7351,"duration_ms":78759,"temperature":0.7,"pith_summary":"The paper is an introductory review that sets out to show that modern cosmological galaxy simulations are built from a small set of standard ingredients: initial conditions seeded by measured cosmic microwave background fluctuations, numerical solvers for gravity and gas dynamics, sub-grid models for star formation and feedback, and a post-processing and validation stage. It argues that these pieces fit together as a pipeline that lets simulations reproduce many observed galaxy properties, and that unresolved baryonic physics is intentionally treated with 'effective' recipes calibrated to observations. If the review is right, a newcomer can acquire a working map of the field from one chapter, and the key open problems, such as resolution limits, sub-grid degeneracies, and the circumgalactic medium, are clearly identifiable.","feed_headline":"A full pipeline for simulating galaxies, from CMB to feedback","feed_subtitle":"An introductory review lays out the ingredients—gravity solvers, gas physics, star formation, AGN feedback—and how simulations are…","key_machinery":"The organizing object is the simulation pipeline itself, from initial conditions to mock observations. Initial conditions are drawn from a Gaussian random field with power spectrum $P(k)=A_s k^{n_s} T^2(k)$, with the amplitude and shape fixed by CMB measurements and the transfer function $T(k)$ computed by Boltzmann solvers. The load-bearing mechanism within the pipeline is the baryonic cycle: gas cools and accretes onto galaxies, forms stars, and is ejected back into the circumgalactic medium by stellar and AGN feedback, which regulates further star formation. Because the relevant scales are unresolved, this cycle is implemented through effective sub-grid prescriptions, including cooling tables, density-threshold star formation, and thermal or kinetic feedback injection, whose parameters are calibrated to reproduce observed scaling relations.","core_discovery":"The central claim, stated in the authors' own terms, is that cosmological galaxy simulations are numerical experiments that follow dark matter, gas, stars, and black holes in an expanding Universe, and that their essential structure is now standardized. The review lays out the chain: a Gaussian random density field with a CMB-calibrated power spectrum is evolved from high redshift using gravity solvers (tree, particle-mesh, tree-PM) and hydrodynamics solvers (smoothed particle hydrodynamics, adaptive mesh refinement, moving mesh); star formation, stellar feedback, and AGN feedback are inserted as sub-grid models; halos and galaxies are identified in post-processing; and the results are checked by convergence tests, parameter variations, cross-code comparisons, and comparison with observed scaling relations. The review maintains that this pipeline reproduces key observed galaxy properties and that remaining discrepancies concentrate in the circumgalactic medium, faint low-mass galaxies, and low-surface-brightness features. It closes by arguing that next-generation simulations will push resolution, add physics such as magnetic fields, cosmic rays, and non-equilibrium cooling, and incorporate machine-learning emulators.","pith_inferences":["If this pipeline description is correct, the field's reproducibility would be improved by standardizing initial conditions and validation metrics across codes, since the review shows the ingredient list is already shared.","Because the review deliberately excludes magnetohydrodynamics, cosmic rays, radiation hydrodynamics, and conduction, current flagship predictions for circumgalactic gas and high-redshift galaxies may shift once these processes become standard; a reader should treat those predictions as provisional.","A natural testable extension is to run identical feedback models in different hydrodynamics solvers at matched resolution and measure how much of the spread in predicted galaxy properties is numerical rather than physical; the review's own cross-code discussion suggests this spread is largest in the circumgalactic medium.","The review's emphasis on observation-calibrated effective models implies that predictions for observables outside the calibration set, such as detailed interstellar-medium phase structure or faint low-surface-brightness features, are the most likely place for model failures to show up."],"forward_implications":["Any modern galaxy simulation can be understood by identifying its initial-condition generator, its gravity and hydrodynamics solvers, and its sub-grid feedback choices.","Simulation success is judged by reproduction of observed scaling relations such as the stellar mass function, the Kennicutt-Schmidt relation, and the mass-metallicity relation, rather than by resolving every microphysical process.","The missing-satellites problem is presented as largely resolved within the standard cosmological model once stellar feedback and environmental effects are included.","Cross-code comparison projects that fix initial conditions and physical models expose where predictions are stable and where they are code-dependent, with the circumgalactic medium highlighted as the main site of disagreement.","Next-generation simulations are expected to extend resolution, add magnetic fields, cosmic rays, non-equilibrium cooling, and thermal conduction, and accelerate analysis with machine-learning emulators."],"supporting_citations":[{"why":"Sets the review's framing of simulations as tools for galaxy formation theory and surveys the field.","marker":"Somerville & Davé, 2015"},{"why":"Provides the companion review of numerical methods and galaxy formation results that this chapter distills.","marker":"Naab & Ostriker, 2017"},{"why":"Is the stated reference for cross-simulation comparisons and for where persistent discrepancies, especially in the circumgalactic medium, remain.","marker":"Crain & van de Voort, 2023"},{"why":"Supplies the cosmological parameters and fluctuation amplitude used to initialize Gaussian random density fields.","marker":"Planck Collaboration et al., 2020"},{"why":"Is the canonical adaptive mesh refinement hydrodynamics solver cited for Eulerian methods.","marker":"Teyssier, 2002"},{"why":"Introduces the moving-mesh Voronoi hydrodynamics technique used by Arepo, representing arbitrary Lagrangian-Eulerian solvers.","marker":"Springel, 2010b"},{"why":"Describes EAGLE, the large-volume simulation whose sub-grid feedback model and convergence discussion anchor the review's treatment of baryonic physics.","marker":"Schaye et al., 2015"},{"why":"Defines the FIRE zoom-in simulation framework that resolves the interstellar medium and stellar feedback, anchoring the high-resolution end of the surveyed landscape.","marker":"Hopkins et al., 2018"},{"why":"Presents Illustris, a flagship cosmological-volume simulation demonstrating baryonic feedback in a full cosmological context.","marker":"Vogelsberger et al., 2014"}],"fun_headline_variants":["The galaxy sim pipeline: CMB to feedback","Simulating galaxies: the full recipe explained","From CMB to galaxy formation: a simulation guide","The standard pipeline for galaxy simulations","Inside galaxy simulations: from initial conditions to stars"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The review's pedagogical claim collapses if its condensed taxonomy of solvers and sub-grid models, together with the deliberately excluded physics, misrepresents the essential ingredients enough to mislead a newcomer.","fun_headline_variants_meta":{"raw":{"variants":["The galaxy sim pipeline: CMB to feedback","Simulating galaxies: the full recipe explained","From CMB to galaxy formation: a simulation guide","The standard pipeline for galaxy simulations","Inside galaxy simulations: from initial conditions to stars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000565,"raw_usage":{"total_tokens":2671,"prompt_tokens":930,"completion_tokens":1741,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":546,"completion_tokens_details":{"reasoning_tokens":1673}},"tokens_in":546,"tokens_out":1741,"duration_ms":16602,"temperature":1.0,"reasoning_tokens":1673,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T18:08:56.629658+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run two hydrodynamics solvers with identical initial conditions, cooling tables, and feedback prescriptions, and check whether their predicted galaxy stellar mass functions by $z=0$ agree within observational uncertainties; a disagreement larger than those uncertainties would call into doubt the review's claim that the pipeline is stable and adequately validated.","supporting_citations":[],"review_version":1}