{"id":"06ba9d12-2270-4c14-adc0-80d80e2ab5ef","arxiv_id":"2608.01632","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In isolated MOND simulations, the UGC 7321 model keeps h_z/R_D mostly below 0.1 for 5 Gyr, and higher-MOND-depth (more diffuse) models stay thinner than low-depth models.","lead":"Simulations show that a model of the superthin galaxy UGC 7321 can stay extremely flat for five billion years under Milgromian gravity, even after forming a bar and buckling early. It is the first systematic MOND test of why these rare razor-thin galaxies survive.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Fiducial UGC 7321 mass is ~2x too low relative to the paper's own MOND rotation-curve fit; the higher-mass regime (M2) thickens beyond h_z/R_D=0.1, so the 'observationally constrained' demonstration may not apply to the real galaxy.","rationale":"The paper's central claim is a dynamical-viability statement about real superthin galaxies, and the UGC 7321 model is the observational anchor. The strongest test of that anchor is whether the adopted baryonic mass is compatible with the observed rotation curve in MOND. It is not: the Fig. 2 fit requires Upsilon_B=2.99, giving M_star about twice the Table 1 value, and the authors explicitly note the DICE model underpredicts the outer rotation. Since MOND rotation curves are fixed by baryonic mass, a model that misses the rotation is not a valid representation of UGC 7321 in MOND. The same paper shows that doubling the mass (M2) produces substantial thickening beyond the superthin threshold, so the mass choice directly biases the headline result. This is an internal inconsistency between the adopted model and the paper's own best-fit mass, not a matter of outside consensus. A targeted re-run at the fitted mass would settle whether the 'observationally constrained' conclusion survives. The D_M comparison is also confounded with surface density, but that is a secondary issue; the mass normalization is the single load-bearing concern because it determines whether the fiducial model is the right test.","tokens_in":18401,"tokens_out":9346,"duration_ms":101803,"concrete_test":"Re-run the fiducial UGC 7321 setup with the stellar mass set to the best-fit MOND value, M_star = 2.99 × 1.1e9 Msun = 3.29e9 Msun, leaving M_gas, scale lengths, vertical structure, and resolution unchanged. Measure the radial h_z/R_D profile at t=5 Gyr. If the disc exceeds h_z/R_D=0.1 over a substantial radial range (similar to M2 in Fig. 9), the 'observationally constrained' demonstration fails for the real galaxy; if it remains below 0.1, the mass-normalization concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that superthin discs can remain vertically thin over 5 Gyr of isolated MOND evolution, demonstrated with a model labeled 'observationally constrained' for UGC 7321. However, the adopted stellar mass (M_star=1.6e9 Msun, Table 1) is inconsistent with the paper's own MOND rotation-curve fit (Fig. 2, blue), which requires Upsilon_B=2.99 and hence M_star~3.3e9 Msun for L_B=1.1e9 Lsun. The DICE model (red) underpredicts the observed outer rotation curve. This mass choice sets D_M=0.81 and h_z/R_D(0)=0.057, both favorable to thinness. Model M2, with double the baryonic mass (D_M=0.69), thickens beyond h_z/R_D=0.1 over a substantial radial range (Fig. 9). A model with the fitted stellar mass (M_bar~4.8e9 Msun) would lie between M1 and M2, so the trend implies appreciable thickening. Thus the simulation may only demonstrate that a hypothetical low-mass, high-D_M disc stays thin, not that a galaxy matching UGC 7321's observed kinematics does. The authors acknowledge the underpredicting rotation in Sec. 2 but retain the low-mass model; this unresolved tension weakens the observational grounding of the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses the Phantom of RAMSES (POR) code to run three-dimensional QUMOND hydrodynamical N-body simulations of an isolated, superthin, low-surface-brightness galaxy modeled on UGC 7321. The fiducial model consists of exponential stellar and gas discs with parameters chosen from photometric and HI observations, evolved for 5 Gyr. The paper reports that the disc develops a strong bar and a buckling episode near 1 Gyr, but the vertical heating is limited and most of the stellar disc remains below h_z/R_D < 0.1 at the end of the simulation. The paper also introduces the MOND depth index D_M (Eq. 14) and compares five models (M1–M5) that vary total baryonic mass or disc scale lengths, finding that lower-D_M models tend to heat and thicken more, while higher-D_M models remain thinner. The main conclusions are that superthin discs can remain vertically thin in isolated MOND evolution and that the long-term maintenance of thinness is at least partly related to the degree to which the galaxy lies in the low-acceleration regime.","tokens_in":18757,"tokens_out":4375,"duration_ms":54255,"significance":"If the central claim is correct, the paper would extend the empirical successes of MOND from rotation curves to the long-term vertical structure of the thinnest known disc galaxies, a regime in which standard ΛCDM simulations have difficulty producing such flat discs. The study has notable strengths: it uses a publicly available, well-tested code; the diagnostics (Fourier amplitudes, h_z/R_D profiles, vertical velocity dispersion, vertical heating, and restoring force) are standard; the model parameters are reported in enough detail to reproduce the runs; and the authors are transparent about numerical resolution and the caveats on the isothermal gas treatment. The five-model D_M suite is a reasonable first step. However, the observational anchoring of the fiducial model is undermined by the paper's own rotation-curve fit, which requires roughly twice the adopted stellar mass. Since the higher-mass model M2 thickens beyond the superthin threshold over a substantial radial range, the paper does not currently establish that a galaxy matching UGC 7321's observed kinematics would remain superthin in MOND.","major_comments":[{"comment":"The 'observationally constrained' fiducial model is inconsistent with the paper's own MOND rotation-curve fit. The fit in Fig. 2 requires Υ_B=2.99, which with L_B=1.1e9 L_sun implies M_star≈3.3e9 M_sun, whereas Table 1 uses M_star=1.6e9 M_sun and M_gas=1.5e9 M_sun. The DICE model (red curve) visibly underpredicts the observed outer rotation curve. This low mass sets D_M=0.81 and the favorable initial h_z/R_D=0.057. Model M2, with M_bar=6.2e9 M_sun and D_M=0.69, exceeds h_z/R_D=0.1 over a substantial radial range (Fig. 9). The mass implied by the rotation-curve fit lies between M1 and M2, so the trend suggests realistic UGC 7321 could thicken well beyond the fiducial result. The central claim that an observationally constrained UGC 7321 remains superthin is therefore not established. The authors should either rerun the fiducial model at the fitted mass normalization (or an intermediate va","section":"§2, Fig. 2, Table 1"},{"comment":"The D_M comparison does not isolate MOND depth as a causal variable. Models M2 and M3 have identical D_M=0.69 but differ in total baryonic mass and scale lengths (Table 2), and they show different vertical heating and morphology; M4 and M5 likewise share D_M=0.89 with different outcomes. Thus the apparent trend 'lower D_M implies stronger vertical heating and thickening' is confounded with total baryonic mass, surface density, and compactness. The text acknowledges that D_M alone does not uniquely determine evolution, but the Abstract and Conclusions still state that maintenance is influenced 'at least partly' by the degree of low-acceleration regime. To support that causal statement, the authors need a model set that varies D_M while holding mass and scale lengths fixed (e.g., via a0 or an external field) or a quantitative decomposition separating these effects. As presented, the eviden","section":"§3.2.2, Figs. 9–10"},{"comment":"The resolution comparison is reported only qualitatively for the vertical-thickness claim. The high-resolution fiducial and the lower-resolution M1 have the same physical parameters but different peak bar strengths, buckling times, and late-time m=2 amplitudes; the text says their edge-on morphologies and h_z/R_D profiles both indicate limited thickening, but no quantitative h_z/R_D comparison is shown. Because vertical heating in this study is driven by the bar and buckling, which are explicitly resolution-sensitive, the central claim 'the disc remains globally superthin after 5 Gyr' needs a quantitative resolution check, e.g., a plot of Δ(h_z/R_D) at late times for the two runs, or an error bar on the final thickness. Without this, the reader cannot judge whether the survival result is converged.","section":"§3.2.2"}],"minor_comments":[{"comment":"Minor grammatical issue: 'The radial gas scale-height profile constructed following the method of Banik et al. (2020)' should read '... profile was constructed following...'.","section":"§2 text"},{"comment":"The caption uses 'D = 10.00 Mpc, i = 88.0' while the text states D=10.0 Mpc and i=88°. Please use consistent notation and units throughout.","section":"Fig. 2 caption"},{"comment":"The angular-momentum flux tensor equation does not show the local average indicated by the overbar in the text, and the index contraction is implicit. Adding an explicit overbar on ρ_star v_k v_α or stating 'local average' directly in the equation would improve clarity.","section":"Eq. (11)"},{"comment":"The colorbar label 'log10 ( /10^3 M pc^-2)' is ambiguous; please specify the quantity explicitly, e.g., log10(Σ / (10^3 M_sun pc^-2)), and state the projected mass surface density of which component is shown.","section":"Fig. 4 colorbar"},{"comment":"The statement about |K_z|/(2πG) having dimensions of surface density but not being a direct baryonic surface density is useful; consider adding a one-sentence reminder when interpreting Fig. 7 in the text, since the shape of this profile is close to a surface-density profile.","section":"§3.1.2"}],"recommendation":"major_revision","confidential_remarks":"I do not see circularity in the use of D_M: it is computed from the input mass distribution via Eq. (14) before evolution, and the outcome h_z/R_D is not used to readjust D_M. The main issue is calibration: the paper's own algebraic MOND fit requires roughly twice the adopted stellar mass, and the model with that approximate mass (M2) does not stay superthin. The authors should address this directly, either with a new simulation or by clearly downgrading the claim to a hypothetical low-mass model. The paper also needs to justify the causal language about D_M despite the M2/M3 and M4/M5 degeneracies. These are fixable within the manuscript's scope, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper asks the right question—can MOND keep a superthin disc thin for 5 Gyr—and the answer it gets is plausible but softer than the abstract implies. The UGC 7321 model they call observationally constrained uses a stellar mass half of what their own MOND rotation-curve fit requires. That matters because the doubled-mass model M2 is the one that thickens past h_z/R_D=0.1.\n\nGenuinely new: first systematic 3D hydrodynamical MOND simulations of long-term vertical evolution of a superthin galaxy. The model sequence using D_M (MOND depth) is a useful construct, and the diagnostics—Fourier amplitudes, h_z/R_D, sigma_z, vertical restoring force—are standard and well presented. They are also transparent: they admit the DICE model underpredicts the outer rotation curve, note the resolution sensitivity of the bar, and flag the gas-temperature choice as a free parameter. That honesty is real.\n\nSoft spots, in order. First, the mass normalization. The adopted M_star=1.6e9 Msun comes from Bell & de Jong, but their own one-parameter MOND fit to the same photometry needs Upsilon_B=2.99, i.e., M_star~3.3e9. They keep the low-mass model because their aim is an idealized superthin disc, not a rotation-curve fit. Fair enough, but then \"observationally constrained\" overstates the case: a model that misses the observed kinematics by the amount shown in Fig. 2 is not constrained by those kinematics. Since M2 (double mass) thickens substantially, the actual UGC 7321 probably sits in an uncomfortable middle ground. The broad conclusion—diffuse, high-D_M discs can stay superthin—survives, but the specific demonstration for UGC 7321 does not.\n\nSecond, reproducibility. They use public POR and DICE, but ship no configuration files or initial conditions. With only single realizations per model, no seed variation is tested. The resolution check (levelmax 12 vs 13) shows the non-axisymmetric amplitudes are resolution-sensitive; vertical thickness is less affected, so that is minor, not fatal.\n\nThird, only one galactic archetype. They acknowledge this in the conclusions.\n\nOverall: the central mechanism—MOND's vertical restoring field can offset bar/buckling heating in low-surface-density discs—is a plausible and genuinely interesting result that deserves attention. The mass normalization inflates the strength of the UGC 7321-specific conclusion. This paper warrants a serious referee, and a good referee will send it back for a higher-mass fiducial model or a rewritten claim.\n\nRecommendation: yes, send to peer review.","headline":"First systematic MOND hydro sims of superthin survival; the broad trend with MOND depth is credible, but the UGC 7321-specific claim rests on a low-mass model that its own rotation-curve fit contradicts.","tokens_in":19285,"tokens_out":3267,"would_cite":true,"duration_ms":37068,"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":"Simulations show that a MOND model of UGC 7321 keeps most of its stellar disc superthin after 5 Gyr, despite a strong bar and a buckling episode.","keywords":["superthin galaxies","MOND","QUMOND","vertical disc heating","bar buckling","UGC 7321","MOND depth index","N-body simulations"],"falsifier":"Re-run the fiducial simulation with the stellar mass increased to the MOND fit value (M/L_B ≈ 2.99, roughly doubling the stellar mass) and measure h_z/R_D after 5 Gyr; if a substantial radial range exceeds 0.1, the claim that a mass-consistent MOND model of UGC 7321 stays superthin is falsified.","tokens_in":18263,"feed_emoji":"🌌","tokens_out":8306,"duration_ms":88595,"temperature":0.7,"pith_summary":"Superthin galaxies are edge-on discs with vertical scale height less than about a tenth of their radial scale length; keeping them that flat in ordinary Newtonian gravity usually requires unusually compact dark-matter haloes. This paper asks whether Milgromian dynamics (MOND) can do it instead, and simulates an observationally constrained model of the prototype UGC 7321 for 5 Gyr. The disc develops a strong bar and buckles at about 1 Gyr, yet most of the stellar disc still ends below h_z/R_D < 0.1. A model sequence spanning the MOND depth index D_M shows that more diffuse, deeper-MOND discs stay thinner, while more massive or compact systems heat and thicken more. The conclusion is that long-lived superthin discs are dynamically viable under MOND, so the special halo conditions invoked for them may not be necessary.","feed_headline":"Superthin discs survive 5 Gyr of isolated MOND evolution","feed_subtitle":"An observationally tuned model of UGC 7321 stays below h_z/R_D = 0.1 despite a strong bar and buckling.","key_machinery":"The argument is carried by the quasi-linear formulation of MOND called QUMOND. The Newtonian potential is first solved from the baryonic density, then a second Poisson equation with the interpolating function ν(y) = (1 + sqrt(1 + 4/y))/2 builds the Milgromian potential; this is the gravity law the simulations evolve. The second ingredient is the MOND depth index D_M = 1 − M_bar(< r_M)/M_bar, the fraction of baryonic mass outside the MOND radius r_M = sqrt(G M_bar/a0), which quantifies how much of the galaxy lies in the low-acceleration regime. The third is the simulated vertical restoring field K_z, measured near the disc plane; its persistence through the run shows that baryonic matter alon","core_discovery":"The central claim is that a superthin baryonic disc can survive long-term isolated evolution in MOND without being puffed up by its own instabilities. In the UGC 7321 model, the stellar disc forms a strong bar and undergoes a buckling episode at t ≈ 1 Gyr, yet after 5 Gyr the average h_z/R_D has increased by only about 1.5–1.8 times and most of the disc remains below the 0.1 superthin threshold. The authors interpret this as a balance between heating from non-axisymmetric structures and vertical confinement by the Milgromian gravitational field. Across the M1–M5 sequence, lower-D_M models (higher baryonic mass or smaller scale length) heat more and thicken more, while higher-D_M models (lowe","pith_inferences":["If the fiducial mass were raised to match the paper's own best-fit MOND rotation curve (M/L_B ≈ 2.99), the galaxy would sit near model M2, which exceeds h_z/R_D = 0.1 over a substantial radial range; so the superthin conclusion may hinge on the photometric mass scale.","The MOND-depth trend suggests a testable mapping between observable surface brightness and diffuseness and vertical thickness among edge-on LSB galaxies: high-D_M candidates should show the smallest h_z/R_D at fixed age.","Because these runs are isolated, the warped outer H I disc of UGC 7321 likely requires an external field or accretion; the paper's setup could be extended with static or time-varying external fields to test both warp and vertical heating.","Including a more realistic multiphase gas or star-formation feedback could change the early bar strength and hence the heating; the paper's isothermal gas treatment is a simplification that future simulations should vary."],"forward_implications":["If MOND is the correct description, superthin galaxies no longer require compact dark-matter haloes; the baryonic disc itself generates enough vertical restoring force to survive 5 Gyr of isolated evolution.","Bars and bar buckling do not by themselves destroy superthin discs in MOND: the simulated bar weakens after buckling and the disc remains thin, matching the boxy/peanut isophotes seen in UGC 7321.","The D_M trend predicts a population correlation: galaxies whose baryonic mass lies mostly inside the MOND radius (low D_M) should be thicker and warmer, while high-D_M galaxies should be the thinnest and most persistent.","D_M alone is not enough: models with equal D_M, like M2 and M3, evolve differently, so mass and scale length still matter for predicting vertical structure."],"supporting_citations":[{"why":"Supplies the observed H I rotation curve, H I mass, distance, and inclination used to build and test the UGC 7321 model.","marker":"Uson & Matthews (2003)"},{"why":"Provides the photometric structure of UGC 7321, including the stellar scale length and color used to set the initial stellar disc.","marker":"Matthews et al. (1999)"},{"why":"Provides the Newtonian mass model and central surface-density calibration that the initial conditions inherit.","marker":"Banerjee et al. (2010)"},{"why":"Supplies the B-band colour-to-mass-to-light calibration used to set the stellar mass of the fiducial model.","marker":"Bell & de Jong (2001)"},{"why":"Defines QUMOND, the modified-gravity formulation that the simulations solve.","marker":"Milgrom (2010)"},{"why":"Introduces Phantom of RAMSES, the MOND N-body and hydrodynamical code used for the evolution.","marker":"Lüghausen et al. (2015)"},{"why":"Supplies the method for initialising the gas scale height in MOND vertical equilibrium and the effective-temperature choice.","marker":"Banik et al. (2020)"},{"why":"Defines the MOND depth index D_M used to build the M1–M5 model sequence.","marker":"Eappen & Kroupa (2026)"}],"fun_headline_variants":["MOND keeps superthin discs flat despite buckling","Superthin galaxy stays thin after 5 Gyr in MOND","Bar and buckling don't ruin superthin galaxies in MOND","Superthin discs resist heating in MOND simulations"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"Everything rests on assuming the photometric stellar mass of 1.6e9 solar masses is right; the paper's own MOND rotation-curve fit requires about twice that, and the adopted model under-predicts the observed outer rotation, so if the galaxy really is heavier the disc may not stay superthin.","fun_headline_variants_meta":{"raw":{"variants":["MOND keeps superthin discs flat despite buckling","Superthin galaxy stays thin after 5 Gyr in MOND","Bar and buckling don't ruin superthin galaxies in MOND","Superthin discs resist heating in MOND simulations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000718,"raw_usage":{"total_tokens":3142,"prompt_tokens":908,"completion_tokens":2234,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":652,"completion_tokens_details":{"reasoning_tokens":2168}},"tokens_in":652,"tokens_out":2234,"duration_ms":16246,"temperature":1.0,"reasoning_tokens":2168,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T23:51:35.010300+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the fiducial simulation with the stellar mass increased to the MOND fit value (M/L_B ≈ 2.99, roughly doubling the stellar mass) and measure h_z/R_D after 5 Gyr; if a substantial radial range exceeds 0.1, the claim that a mass-consistent MOND model of UGC 7321 stays superthin is falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the observed H I rotation curve, H I mass, distance, and inclination used to build and test the UGC 7321 model."},{"cited_title":"D., Gallagher, J","cited_arxiv_id":null,"evidence_quote":"Provides the photometric structure of UGC 7321, including the stellar scale length and color used to set the initial stellar disc."},{"cited_title":"D., & Jog, C","cited_arxiv_id":null,"evidence_quote":"Provides the Newtonian mass model and central surface-density calibration that the initial conditions inherit."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the B-band colour-to-mass-to-light calibration used to set the stellar mass of the fiducial model."},{"cited_title":"2020, ApJ, 905, 135","cited_arxiv_id":null,"evidence_quote":"Supplies the method for initialising the gas scale height in MOND vertical equilibrium and the effective-temperature choice."},{"cited_title":"& Kroupa, P","cited_arxiv_id":null,"evidence_quote":"Defines the MOND depth index D_M used to build the M1–M5 model sequence."}],"review_version":1}