{"id":"6a4fb525-ed08-4f5b-a7bb-904d50d8c18e","arxiv_id":"1908.03989","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"In subsonic, super-Alfvénic, high-beta MHD turbulence, the density-magnetic field anticorrelation weakens and the density-fluctuation versus Mach number slope steepens as the adiabatic index increases from isothermal to 5/3.","lead":"This paper uses 30 MHD simulations to test how the gas equation of state and cooling change the statistics of magnetized turbulence. It finds that the known anticorrelation between density and magnetic field weakens as the gas becomes less isothermal, which complicates inferring Mach numbers from observations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The steep gamma=5/3 free-free slopes in Table 2 may be inflated by a runaway-cooling run that violates the paper's own stationarity criterion.","rationale":"The reader's weakest assumption concerns whether the idealized cooling functions approximate real optically thin cooling well enough for the gamma-dependences to be generic. My concern is narrower and more internal: one simulation included in the all-simulation fits is explicitly non-stationary, undergoing runaway cooling at t = 4T, despite the paper's stated rule that only stationary snapshots from 5T to 10T enter the statistics. Because the gamma=5/3 free-free branch is also called thermally marginally stable at Ms ~ 0.5, the distinctive results for that branch, the steepest density-slope (0.664) and the weakest rho-B anticorrelation (-0.55), could be consequences of proximity to thermal instability rather than of gamma and cooling per se. This is the most load-bearing issue because the paper's headline quantitative claims are slope and correlation differences between branches, and a single non-stationary high-Ms point can materially change a linear fit over an Ms range of only 0.2-0.6. The paper deserves credit for flagging the marginal and runaway behavior in Section 3.4, but it does not say whether those runs were excluded from the fits. A focused reanalysis can settle this quickly, and the central trend may well survive; but until that is shown, the acceptance should be conditional on this robustness check.","tokens_in":18906,"tokens_out":5811,"duration_ms":66904,"concrete_test":"Recompute the Table 2 regressions and the Fig. 5 moments for the gamma=5/3 free-free branch with M0.70P0.37gamma=5/3 Cool:ff A1.00H removed, and as a sensitivity check also remove the marginally stable runs M0.54P0.58gamma=5/3 Cool:ff A1.00H and M0.50P0.74gamma=5/3 Cool:ff A1.00H, restricting all fits to stationary snapshots with 5T <= t <= 10T. Then compare the resulting sigma_ln rho versus Ms slope and the rho-B correlation coefficients with the published values. If the slope remains near 0.66 and corr[rho,B] remains near -0.55, the concern is resolved; if the slope drops toward the linear-cooling value 0.575 or the anticorrelation strengthens appreciably, the central gamma/cooling trend is partly a thermal-instability artifact.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Section 2.3 restricts all statistical results to stationary snapshots in 5T <= t <= 10T, with statistics calculated over 51 snapshots. However, Section 3.4 describes run M0.70P0.37gamma=5/3 Cool:ff A1.00H as undergoing runaway cooling at t = 4T, and Table 1 labels its values as final values at t = 4T. Section 3.3 nevertheless states that moments are computed for all 30 simulations, and Section 3.3 and Table 2 present linear fits over all outputs of all simulations for each EOS/cooling combination. If that non-stationary run was not excluded from the regressions, the gamma=5/3 free-free density-slope m = 0.664(5) and the weakened rho-B anticorrelation could be controlled by transient high-density tails from thermal instability rather than by the thermodynamics of balanced turbulence. The paper does not state that this point was excluded. This matters because the runaway run anchors the high-Ms end of the fit, and the gamma=5/3 free-free branch at Ms ~ 0.5 is itself described as thermally marginally stable, so the reported cooling-function dependence may be entangled with thermal stability instead of being a generic thermodynamic effect.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a suite of 30 ideal MHD simulations of driven, thermally balanced turbulence in the subsonic, super-Alfvenic, high-beta regime, with the equation of state varied between approximately isothermal, gamma = 7/5, and gamma = 5/3, and with no cooling, linear cooling, or free-free-like cooling used to balance turbulent dissipation. The central results are that the energy spectra are insensitive to thermodynamics, that the thermal-magnetic pressure anticorrelation is essentially independent of gamma and cooling (about -0.8), while the density-magnetic field anticorrelation weakens systematically with increasing gamma (from about -0.8 isothermal to about -0.55 for gamma = 5/3 with free-free cooling), and that the slopes of the linear relations between density/pressure fluctuation amplitudes and sonic Mach number steepen with gamma. The paper argues this creates a degeneracy between thermodynamics and Mach number in interpreting observations, which may be broken by higher-order moments.","tokens_in":19077,"tokens_out":11706,"duration_ms":111195,"significance":"If the reported gamma-dependences are robust, this is a useful and timely quantification of a usually neglected thermodynamic axis in MHD turbulence simulations, with direct implications for interpreting Faraday-rotation-based and density-fluctuation-based diagnostics in ICM-like plasmas. The study is carefully set up: it uses a fixed numerical scheme with resolution pairs at 512^3 and 1024^3, an explicit stationarity interval (5T-10T), a published code fork, and a candid limitations section that acknowledges the idealized cooling functions and the narrow parameter range. The comparison to previous isothermal and hydrodynamic work (Nolan et al. 2015; Mohapatra & Sharma 2019) is useful. The main caveat is that the quantitative slopes and the cooling-function dependence rest on a small number of simulations and, as detailed below, on the treatment of one non-stationary run.","major_comments":[{"comment":"The manuscript does not state whether the non-stationary, runaway-cooling run M0.70P0.37gamma=5/3 Cool:ff A1.00H is included in the regressions of Table 2 and in the correlation and moment plots. Section 2.3 restricts all statistical results to the stationary interval 5T <= t <= 10T, but Table 1 gives only final values at t = 4T for this run, while Sec. 3.3 says moments are computed for all snapshots of all 30 simulations and that regressions use all outputs of all simulations. If the runaway run contributes, the high-Ms endpoint that anchors the gamma=5/3 free-free slope m = 0.664(5) is a transient thermal-instability tail rather than balanced turbulence, so the steepening with gamma and the weakened rho-B anticorrelation would not be clean thermodynamic statements. Please state explicitly whether this run is excluded, and if so, report the fitted slopes and correlations with the run excluded; the gamma=5/3 free-free branch is also described as only marginally stable at Ms ~ 0.5, so the sensitivity of the fits to this branch's stability should be quantified.","section":"Secs. 2.3, 3.3, and Tables 1-2"},{"comment":"The text states that higher resolution leads to slightly weaker rho-B anticorrelations of about the same order as the cooling functions, yet no numerical estimate of this resolution effect is given, and Appendix A demonstrates convergence only for PDFs, not for the correlation coefficients. Because the cooling-function difference at fixed gamma (e.g., -0.64(3) for linear versus -0.55(3) for free-free at gamma = 5/3) is one of the paper's quantitative results, please report the resolution-pair differences from the gray-shaded simulations in Fig. 3 and show explicitly that the cooling trend survives after accounting for resolution.","section":"Sec. 3.2 and Appendix A"},{"comment":"The slopes in Table 2 are fit to 51 temporally correlated snapshots per simulation, with only 3-4 simulations per EOS/cooling combination (and 15-16 for gamma = 5/3 free-free) over the narrow range 0.2 < Ms < 0.6. The quoted standard errors, e.g., 0.664(5), therefore reflect within-run scatter rather than simulation-to-simulation variance, and the comparisons of slopes across EOS and cooling functions may not be statistically robust. Please provide a fit or bootstrap that treats each simulation as an independent sample, for example using time-averaged moments per run, or otherwise quantify the inter-run uncertainty.","section":"Sec. 3.3 and Table 2"}],"minor_comments":[{"comment":"The abstract and Sec. 4.2 contain awkward or ungrammatical phrasings such as the physical properties turbulence and Similarly, with the linear relation between variations in density and thermal pressure with sonic Mach number becomes steeper with increasing gamma; these should be corrected.","section":"Abstract and Sec. 4.2"},{"comment":"Section 3.3 states that linear fits are performed only when the absolute correlation exceeds 0.9, but Table 2 includes a fit for the isothermal pth kurtosis with Corr = 0.86; please reconcile the criterion with the reported fits.","section":"Sec. 3.3 and Table 2"},{"comment":"The caption of Fig. 7 is internally inconsistent about line styles, referring to solid, dashed, and dotted lines in ways that do not agree with the described convergence comparison; the visual distinction between 512^3 and 1024^3 runs should be made unambiguous.","section":"Fig. 7 caption"},{"comment":"The runaway run M0.70P0.37gamma=5/3 Cool:ff A1.00H is listed without temporal standard deviations; a brief footnote explaining that these are single-epoch final values would help avoid confusion with the stationary rows.","section":"Table 1"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this is a careful numerical study of how the equation of state and cooling change the statistics of subsonic, high-beta MHD turbulence. The main trends—weaker density–magnetic-field anticorrelation with increasing gamma, and steeper density-fluctuation versus Mach-number slopes—are visible in the stationary runs and are probably robust. The paper deserves peer review, not desk rejection.\n\nWhat is new: a systematic sweep across isothermal, gamma=7/5, and gamma=5/3 with two idealized cooling functions in the MHD case. The gamma-dependence of the rho–B correlation and the steepening of sigma_rho with Ms are new relative to the isothermal MHD and hydrodynamic-adiabatic literature. The finding that the thermal–magnetic pressure anticorrelation is pinned near -0.8 regardless of thermodynamics is clean and useful. The authors are also honest about limitations: idealized cooling, no multiphase behavior, subsonic-only, and they include resolution checks.\n\nThe soft spots are real but not fatal. The stress-test about run M0.70P0.37gamma=5/3 Cool:ff A1.00H lands. The paper says all statistics come from the stationary window 5T–10T, but Section 3.3 says moments are computed for “all snapshots of all 30 simulations,” and Table 2 regressions use all outputs per EOS/cooling combination. The runaway run has no stationary snapshots, so it should have been explicitly excluded from the fits. If it entered, the gamma=5/3 free-free density slope of 0.664(5) is inflated by thermal-instability tails, and the paper’s secondary claim that free-free cooling steepens the slope more than linear cooling is shaky. Even without that run, the free-free slopes appear steeper than linear based on the other gamma=5/3 runs, so the central story survives, but the authors need to state the exclusion and rerun the fits. The narrow Mach range (0.2–0.6) is another soft spot; the authors themselves tell readers not to over-interpret the fitted slopes. And the resolution effect on the correlations is comparable to the cooling-function effect, which tempers the finer distinctions.\n\nWho this is for: anyone inferring Mach numbers from density or Faraday-rotation statistics in clusters, and turbulence simulators who need a non-isothermal MHD reference point. The paper is readable and the code is public, though the analysis scripts are not—releasing them would settle the runaway-run question quickly.\n\nRecommendation: send it to peer review. A serious referee should ask for the explicit exclusion of the non-stationary run and a re-fit, but the core is solid and the trend is likely to stand.","headline":"Solid, systematic MHD turbulence study on thermodynamics; the main trends hold, but the runaway-cooling run must be explicitly excluded from the fits.","tokens_in":19662,"tokens_out":6486,"would_cite":true,"duration_ms":58211,"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":"Thermodynamics weakens density–magnetic field link in turbulence","keywords":["MHD turbulence","equation of state","adiabatic index","optically thin cooling","sonic Mach number","density–magnetic field correlation","Faraday rotation","intracluster medium"],"falsifier":"Run the same turbulence suite with a realistic, tabulated optically thin cooling curve for a specific plasma such as the intracluster medium, and compare the density–magnetic field correlation coefficients and the density-fluctuation versus Mach number slopes: if the ordering by $\\gamma$ weakens, reverses, or vanishes, the idealized cooling functions are responsible for the claimed degeneracy.","tokens_in":18656,"feed_emoji":"🌪️","tokens_out":10876,"duration_ms":98258,"temperature":0.7,"pith_summary":"This paper argues that in subsonic, super-Alfvénic, high-beta magnetohydrodynamic turbulence, the equation of state of the gas changes the statistical fingerprints used to interpret astrophysical observations. Through 30 simulations that balance turbulent dissipation with idealized cooling, the authors show that as the adiabatic index rises from isothermal to monatomic, the anticorrelation between density and magnetic field weakens (from about $-0.8$ to about $-0.55$), while the linear slope connecting density fluctuations to sonic Mach number steepens. The result is a degeneracy: the same measured relation between density fluctuations and Mach number can arise from different combinations of equation of state and Mach number. The paper proposes that higher-order moments, such as the skewness of the density distribution, depend mainly on the equation of state and may break that degeneracy.","feed_headline":"Thermodynamics weakens density–magnetic field link in turbulence","feed_subtitle":"Same density–Mach slope can come from different gas states, so observers need higher moments to tell them apart.","key_machinery":"The load-bearing tool is a suite of 30 ideal MHD simulations of stationary, thermally balanced turbulence, run with two idealized optically thin cooling functions—linear cooling $L\\propto\\rho e$ and free-free-like cooling $L\\propto\\rho^2 e^{1/2}$—with the cooling coefficient chosen to balance turbulent dissipation. The analysis rests on correlation coefficients and probability density functions (mean, standard deviation, skewness, kurtosis) of density, thermal pressure, total pressure, and a Faraday-rotation-derived line-of-sight magnetic field strength, computed over 51 snapshots in the stationary regime ($5T\\le t\\le10T$). The key relations are the $\\rho$–$B$ anticorrelation, whose dependence on $\\gamma$ and cooling encodes the shift in the compressible mode balance, and the constant $p_{\\rm th}$–$p_B$ anticorrelation that signals total-pressure equilibrium.","core_discovery":"In the regime studied—subsonic ($M_s\\approx0.2$–$0.6$), super-Alfvénic ($M_a\\approx1.8$), and high plasma $\\beta$ ($10\\lesssim\\beta_p\\lesssim100$)—the thermal–magnetic pressure anticorrelation is essentially fixed near $-0.8$ regardless of thermodynamics, indicating total-pressure equilibrium. The density–magnetic field anticorrelation, however, is set primarily by the adiabatic index: about $-0.81$ for isothermal gas, $-0.71$ for $\\gamma=7/5$, and $-0.64$ (linear cooling) or $-0.55$ (free-free cooling) for $\\gamma=5/3$. The linear relation between the standard deviation of logarithmic density and sonic Mach number steepens with $\\gamma$, from a slope near $0.48$ in the isothermal case to about $0.66$ for $\\gamma=5/3$ with free-free cooling. Because Mach number and thermodynamics shift these relations in the same direction, the authors conclude that inferring Mach numbers from such statistics alone is degenerate, but that higher-order moments may separate the two effects.","pith_inferences":["We infer that published $M_s$ estimates from intracluster-medium-like observations should carry an additional systematic uncertainty for unknown $\\gamma$; the slopes reported here provide a first table for quantifying it.","A testable extension would be to repeat the analysis with a realistic, tabulated cooling curve: if the ordering by $\\gamma$ persists, the idealized cooling functions are not the controlling factor.","The weakening $\\rho$–$B$ anticorrelation implies a shift in the mix of slow and fast magnetosonic modes with $\\gamma$; synthetic polarization or Faraday rotation maps could search for that shift in observed clusters.","We infer that the degeneracy should strengthen in the supersonic regime, where density fluctuations are larger, so extending this suite to higher $M_s$ would sharpen or refute the higher-moment diagnostic."],"forward_implications":["Observational estimates of sonic Mach number from density fluctuation amplitudes carry a hidden thermodynamic dependence: the same slope can correspond to different combinations of $\\gamma$ and $M_s$.","Faraday-rotation-based line-of-sight field strengths stay within about 10 percent of the true field in this regime, so the thermodynamic effect is unlikely to dominate measurement error there.","The $\\rho$–$B$ anticorrelation is not a clean slow-mode diagnostic outside isothermal turbulence, because its strength changes with $\\gamma$ and cooling.","Higher-order moments, especially the skewness of the density distribution, are nearly independent of $M_s$ and mainly track $\\gamma$, offering a path to break the degeneracy.","Kinetic, magnetic, and internal energy spectra are practically identical across equations of state, so spectral shape alone cannot reveal the thermodynamics."],"supporting_citations":[{"why":"introduced the linear relation between density fluctuations and sonic Mach number that this paper extends.","marker":"Padoan et al. (1997)"},{"why":"derived and tested the lognormal density distribution and the density–Mach number relation used as the baseline.","marker":"Passot & Vázquez-Semadeni (1998)"},{"why":"showed that the proportionality constant depends on solenoidal versus compressive forcing, the baseline for comparing slopes.","marker":"Federrath et al. (2008)"},{"why":"predicted steeper slopes for larger adiabatic index in subsonic hydrodynamic turbulence, which the MHD results support.","marker":"Nolan et al. (2015)"},{"why":"supplied the plasma-beta correction to the density–Mach number relation, negligible in the high-beta regime here.","marker":"Molina et al. (2012)"},{"why":"reported the isothermal density–magnetic field anticorrelation used as the reference case.","marker":"Yoon et al. (2016)"},{"why":"provided the isothermal density–field correlation and forcing-correlation-time baseline used in this study.","marker":"Grete et al. (2018)"},{"why":"offered isothermal MHD higher-order moment statistics, though with a different forcing model.","marker":"Kowal et al. (2007)"},{"why":"provided the finite-volume MHD code that this study modified and ran.","marker":"Stone et al. (2008)"},{"why":"described the numerical scheme used in the simulations.","marker":"Stone & Gardiner (2009)"}],"fun_headline_variants":["Gas state mimics Mach number in turbulence relations","Turbulence relations hide gas-state–Mach degeneracy","Higher moments needed to untangle gas state from Mach","Isothermal assumption distorts turbulence-derived Mach numbers","Thermodynamics weakens density–field link in turbulence"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the two idealized cooling functions and the dissipation-balance condition lose energy the way real optically thin plasmas do, closely enough that the reported dependence on $\\gamma$ is generic rather than an artifact of the chosen cooling forms.","fun_headline_variants_meta":{"raw":{"variants":["Gas state mimics Mach number in turbulence relations","Turbulence relations hide gas-state–Mach degeneracy","Higher moments needed to untangle gas state from Mach","Isothermal assumption distorts turbulence-derived Mach numbers","Thermodynamics weakens density–field link in turbulence"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000253,"raw_usage":{"total_tokens":1641,"prompt_tokens":1100,"completion_tokens":541,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":716,"completion_tokens_details":{"reasoning_tokens":466}},"tokens_in":716,"tokens_out":541,"duration_ms":6738,"temperature":1.0,"reasoning_tokens":466,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:55:45.791671+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same turbulence suite with a realistic, tabulated optically thin cooling curve for a specific plasma such as the intracluster medium, and compare the density–magnetic field correlation coefficients and the density-fluctuation versus Mach number slopes: if the ordering by $\\gamma$ weakens, reverses, or vanishes, the idealized cooling functions are responsible for the claimed degeneracy.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"introduced the linear relation between density fluctuations and sonic Mach number that this paper extends."},{"cited_title":"1998, Phys","cited_arxiv_id":null,"evidence_quote":"derived and tested the lognormal density distribution and the density–Mach number relation used as the baseline."},{"cited_title":"S., & Schmidt, W","cited_arxiv_id":null,"evidence_quote":"showed that the proportionality constant depends on solenoidal versus compressive forcing, the baseline for comparing slopes."},{"cited_title":"A., Federrath, C., & Sutherland, R","cited_arxiv_id":null,"evidence_quote":"predicted steeper slopes for larger adiabatic index in subsonic hydrodynamic turbulence, which the MHD results support."},{"cited_title":"Z., Glover, S","cited_arxiv_id":null,"evidence_quote":"supplied the plasma-beta correction to the density–Mach number relation, negligible in the high-beta regime here."},{"cited_title":"2016, The Astrophysical Journal, 831, 85","cited_arxiv_id":null,"evidence_quote":"reported the isothermal density–magnetic field anticorrelation used as the reference case."},{"cited_title":"W., & Beckwith, K","cited_arxiv_id":null,"evidence_quote":"provided the isothermal density–field correlation and forcing-correlation-time baseline used in this study."},{"cited_title":"2007, The Astrophysical Journal, 658, 423","cited_arxiv_id":null,"evidence_quote":"offered isothermal MHD higher-order moment statistics, though with a different forcing model."},{"cited_title":"M., & Gardiner, T","cited_arxiv_id":null,"evidence_quote":"described the numerical scheme used in the simulations."}],"review_version":1}