{"id":"f03d560f-50a4-4d00-89b3-6873898ec4bd","arxiv_id":"2608.10748","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"First relativistic MHD turbulence simulations with self-consistent synchrotron cooling find a cooling-driven two-phase plasma and frequency-dependent polarization and Faraday-rotation variability resembling high-energy source observations.","lead":"Simulations show that synchrotron cooling in relativistic magnetized turbulence splits the plasma into hot, dilute and cold, dense regions, which changes how the turbulence emits and rotates radio waves. The results produce synthetic spectra, polarization maps, and Faraday-rotation measures that resemble observations of blazars, pulsar wind nebulae, and fast radio bursts.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Thermal-instability attribution is not yet established: the driven runs lack a control that removes local cooling feedback, and the paper's own diagnostics point to turbulent compression instead.","rationale":"I focused on the causal attribution of the two-phase structure because it is the headline physical claim and because the emission, polarization, and variability predictions are framed as consequences of it. The numerical machinery is credible: the cooling term is benchmarked in Appendix B, the driving scheme conserves momentum with a clean temperature check (Appendix A), and a resolution/magnetization study is included (Appendix C). The weak point is interpretive. The linear thermal-instability analysis assumes a uniform equilibrium with prescribed heating, while the driven turbulence has no such equilibrium; moreover, the simulations' own density-anisotropy and β–ρ diagnostics indicate that turbulent compression dominates. The decisive missing experiment is a control in which local cooling feedback is switched off while the total radiated power is held fixed. Such a run would settle whether the bimodal phases are a genuine instability product or a consequence of compressible turbulence with cooling. Because this concern affects the explanation of the results rather than the descriptive numerical findings, I retain the reader's CONDITIONAL verdict and recommend no change; the revision should add the control run and soften the mechanism claim if the control preserves the phase structure.","tokens_in":27853,"tokens_out":7768,"duration_ms":91584,"concrete_test":"Run a control simulation identical to the η_syn = 6.4×10^-2 case but with the synchrotron cooling term (Eq. 6) replaced by a spatially uniform cooling rate per unit volume matched to the fiducial run's instantaneous volume-averaged power, so that the total energy loss is the same but the local dependence of cooling on B', ρ, and Θ_t is removed. If the (γρ, Θ_t) PDF and the hot-dilute/cold-dense filling fractions in Figures 3–4 remain bimodal, then turbulent compression plus bulk heat loss, not the local thermal-instability feedback, produces the phase structure; if the bimodality disappears, the instability interpretation is supported. This single 1024^3 run isolates the only term that distinguishes the claimed instability from turbulence with cooling.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim in the Abstract and Section 7—that synchrotron cooling triggers the thermal instability and bifurcates the plasma into hot-dilute and cold-dense phases—rests on applying the uniform-heating linear analysis of Appendix B (Eq. B3) to driven turbulence, in which the heating is time-dependent, spatially non-uniform, and supplied by turbulent dissipation rather than by a prescribed H. The instability criterion therefore does not directly govern the simulations. The manuscript itself weakens the attribution: Section 3.1 reports that density-fluctuation anisotropy 'resembles the one expected for uncooled turbulence rather than a thermal-instability pattern,' and Appendix E states that 'turbulent compression still dominates, because the β–γρ relation does not fully reverse, as would be expected for pure thermal instability.' No uncooled or uniformly cooled control using the same phase-identification thresholds is presented, so the non-zero filling fractions in Figure 4 do not demonstrate a feedback instability; they are also consistent with compressive turbulence whose dense, strongly magnetized cells cool faster. If phase coexistence is primarily a turbulent-compression-plus-cooling effect, the asserted mechanism and the claimed instability-enhanced high-frequency variability are unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript presents three-dimensional driven relativistic MHD turbulence simulations with synchrotron cooling, and computes synthetic synchrotron spectra, linear polarization maps, and Faraday rotation measures from the simulated turbulence. The central claim is that synchrotron cooling triggers the thermal instability, bifurcating the plasma into hot-dilute and cold-dense phases, which enhances high-frequency emission and polarization variability. The authors benchmark their stirring force (Appendix A), validate the synchrotron cooling implementation against 1D and 3D thermal-instability tests (Appendix B), and show convergence with resolution and initial magnetization (Appendix C). They also compare the mean-temperature scaling with the theoretical prediction of Uzdensky (2018). The synthetic diagnostics show plausible qualitative resemblance to blazar, pulsar wind nebula, and FRB observations.","tokens_in":28032,"tokens_out":5074,"duration_ms":48899,"significance":"If the central claim holds, this would be the first relativistic MHD simulation suite to couple driven turbulence with self-consistent synchrotron cooling and to derive observable radiative signatures. The paper's strengths are its careful benchmarking of the driving force (temperature preserved to ~2e-9), the explicit 1D and 3D linear/nonlinear thermal-instability tests, and the resolution/magnetization convergence checks. The qualitative results on frequency-dependent polarization and RM variability are potentially useful for interpreting multi-wavelength observations. The main weakness is that the attribution of the phase structure to the thermal instability is not supported by a control experiment and is internally qualified by the authors' own diagnostics.","major_comments":[{"comment":"The thermal-instability interpretation is not established for the driven runs because the linear instability criterion in Appendix B (Eq. B3) assumes a uniform heating H that balances cooling in a homogeneous equilibrium, whereas in the driven turbulence runs the heating is provided by turbulent dissipation and is spatially non-uniform and time-dependent. The manuscript's own diagnostics weaken the causal claim: Section 3.1 states that the density-fluctuation anisotropy 'resembles the one expected for uncooled turbulence rather than a thermal-instability pattern,' and Appendix E states that 'turbulent compression still dominates, because the β–γρ relation does not fully reverse, as would be expected for pure thermal instability.' Since the abstract and Section 7 attribute the phase coexistence to the thermal instability, this is a load-bearing point. Please add a control simulation without cooling (or with a cooling term that does not depend on local density and temperature), processed with the same phase-identification thresholds, and quantify the filling fractions relative to that control.","section":"Section 3.1, Appendix B (Eq. B3), Appendix E"},{"comment":"The statement that non-zero volume filling fractions for the hot-dilute and cold-dense phases 'confirm the presence of the synchrotron-cooling-induced thermal instability' is overstated because no null-hypothesis control is presented. The thresholds in Figure 4 are defined per run from the percentile p of the temperature and density PDFs, and the filling fraction is the intersection of the two selected tails. Compressive turbulence with density-dependent cooling can also populate the anti-correlated temperature-density quadrants without an instability. A control run with identical driving but without cooling, or with a cooling function that is artificially decoupled from local density/temperature, is needed to separate the instability contribution from the compressive-cooling contribution.","section":"Section 3.2, Figure 4"},{"comment":"The scaling relation S.D.(RM) ∝ sqrt(⟨σ⟩_V) K0(⟨Θ_t⟩_V^{-1})/K2(⟨Θ_t⟩_V^{-1}) in Eq. (12) is calibrated to the three simulations listed in Table 1; the proportionality constant is not quoted, and no independent run is used to validate the relation. The FRB temperature estimate in Section 6.3 (⟨Θ_t⟩_V ∼ 3×10^3 for FRB 121102) relies on this uncalibrated relation, which weakens the quantitative inference. Please provide the fitted proportionality constant with an estimate of scatter, or validate Eq. (12) against an additional run with different η_syn or σ0.","section":"Section 4.2, Table 1, Eq. (12)"}],"minor_comments":[{"comment":"The caption contains a duplicated phrase: 'In all runs, In all runs,'. Please remove the repetition.","section":"Figure 5 caption"},{"comment":"The expression for the Q Stokes emissivity contains ambiguous inline fractions '7Θ24/25_t + 35 / 10Θ24/25_t + 75'; please format these as explicit fractions with parentheses or display equations so the intended ratio is clear.","section":"Eq. (8)"},{"comment":"The derivation of the FRB temperature estimate should state explicitly how Eq. (12) is inverted (including the numerical value of the proportionality constant) so that the reader can reproduce the result.","section":"Section 6.3"},{"comment":"The caveat that no Lorentz-invariant kinetic PSD exists is placed only in an appendix; since the main text reports the kinetic PSD slope (Section 3.3), consider mentioning this frame-dependence caveat there to avoid over-interpretation.","section":"Appendix D"},{"comment":"The note that the forcing is non-causal because perturbations are applied simultaneously in Fourier space would be more informative if the range of forcing wavenumbers (kL/(2π) between 1 and 4) is stated here, as it quantifies the scale separation responsible for the non-causality.","section":"Section 5"}],"recommendation":"major_revision","confidential_remarks":"The paper presents a useful simulation framework and a rich set of synthetic diagnostics. The main concern is the thermal-instability attribution, which is central to the abstract and conclusions but is not backed by a control run and is explicitly qualified in the authors' own Appendix E. This is fixable within the manuscript's scope by adding an uncooled or locally-decoupled-cooling control and adjusting the language, so I do not recommend rejection. The RM scaling relation in Eq. (12) also needs a stated calibration constant or an independent validation before it is used for astrophysical inference. The paper is within the scope of the journal and the simulations appear carefully executed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my read of arXiv:2608.10748. The genuinely new thing is the setup: driven relativistic MHD turbulence with self-consistent synchrotron cooling, run in 3D at 1024^3, plus synthetic spectra, polarization maps, and Faraday rotation measures. That combination hasn't been done before, and it's a useful tool for the community. The implementation looks careful: the stirring force is benchmarked to not heat the plasma, the 1D and 3D thermal-instability tests reproduce earlier work, and a 512^3 run plus a lower-sigma run show the main statistics don't shift much. The mean temperature scaling <Theta_t> ~ eta_syn^{-1/2} matches Uzdensky (2018), so the thermodynamics has independent support.\n\nThe soft spot is the central claim that synchrotron cooling triggers the thermal instability and bifurcates the plasma into hot-dilute/cold-dense phases. The evidence isn't as clean as the abstract and Section 7 suggest. The paper applies the uniform-heating linear analysis from Appendix B to driven turbulence, where heating comes from turbulent dissipation and is neither spatially uniform nor prescribed. More importantly, the manuscript's own diagnostics weaken the attribution: Section 3.1 notes the density-fluctuation anisotropy resembles uncooled turbulence, and Appendix E states that the beta-gamma-rho relation does not fully reverse as would be expected for pure thermal instability. The nonzero filling fractions in Figure 4 are also consistent with compressive turbulence whose dense, strongly magnetized cells cool faster. There is no control run--uncooled or uniformly cooled--using the same phase-identification thresholds. So the mechanism attribution is shaky. If phase coexistence is mostly turbulent compression plus cooling, the claimed instability-enhanced high-frequency variability is unsupported.\n\nThe Faraday rotation part is mostly solid, but Eq. (12) is calibrated to three runs and then used to estimate a temperature for the FRB 121102 screen without error bars or a stated normalization. That is a minor issue if presented as a rough estimate, but the text gives it a quantitative flavor. Separately, no code or data release means the numbers can't be independently checked without contacting the authors.\n\nWho is this for? People doing GRMHD or turbulence-based emission modeling of blazars, PWNe, and FRBs will get value from the diagnostics and the frequency-dependent polarization behavior. It deserves a serious referee. I'd send it out, but the referee should press hard on the thermal-instability attribution and ask for a control run or softened language.","headline":"First driven RMHD turbulence with synchrotron cooling, well executed, but the thermal-instability attribution is not yet established and needs a control run.","tokens_in":765,"tokens_out":960,"would_cite":true,"duration_ms":27253,"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":"Synchrotron cooling, the paper argues, actively regulates relativistic turbulence by triggering a thermal instability that splits the plasma into hot dilute and cold dense phases, shaping emission, polarization, and Faraday rotation.","keywords":["relativistic magnetohydrodynamics","turbulence","synchrotron cooling","thermal instability","Faraday rotation","linear polarization","fast radio bursts","pulsar wind nebulae"],"falsifier":"Repeat the highest-cooling run at twice the linear resolution and check whether the two peaks in the lab-frame density PDF, near $\\sim\\rho_0/5$ and $\\sim 3\\rho_0$, persist with unchanged volume filling fractions; the paper's own appendix shows that the pure thermal instability is resolution-limited, so if the phase structure is numerically set rather than converged, the thermal-instability interpretation would not stand.","tokens_in":27547,"feed_emoji":"🌀","tokens_out":13496,"duration_ms":109166,"temperature":0.7,"pith_summary":"This paper tries to establish that synchrotron radiation acts as an active thermodynamic regulator of relativistic magnetized turbulence, not just an output. In three-dimensional driven-turbulence simulations with self-consistent synchrotron cooling, the balance between turbulent energy injection and radiative losses sets the mean plasma temperature, while the cooling-induced thermal instability pushes the plasma into coexisting hot-dilute and cold-dense phases. The paper then computes synthetic synchrotron spectra, linear polarization maps, and Faraday rotation measures from the simulated volume and shows that the two-phase structure makes high-frequency emission more intermittent and more strongly polarized than low-frequency emission. If correct, a single turbulent-cooling mechanism reproduces qualitative observational trends in blazar flares, pulsar wind nebulae, and repeating fast radio bursts.","feed_headline":"Synchrotron cooling splits turbulent plasma into hot and cold phases","feed_subtitle":"Cooling makes high-frequency emission and polarization vary more, matching blazars, nebulae, and FRBs","key_machinery":"The load-bearing control parameter is the dimensionless synchrotron cooling efficiency $\\eta_{\\rm syn}=4(n_e d_e^3)^{-1}(L/d_e)$, physically the Thomson optical depth of the box, which enters the fluid-frame cooling power $P_{\\rm syn}$ through the gas pressure, the magnetic energy density, and a Maxwell–Jüttner temperature factor $K_3(\\Theta_t^{-1})/K_2(\\Theta_t^{-1})$. The argument advances through a feedback loop: flux freezing makes overdense regions carry stronger magnetic fields, stronger fields increase $P_{\\rm syn}$, the local thermal pressure drops, and surrounding plasma is drawn inward, amplifying density and temperature contrasts until turbulent mixing balances the growth. This thermal-instability loop, together with the equality between turbulent energy injection and radiative losses, sets both the mean temperature and the two-phase structure that the synthetic spectra, polarization maps, and rotation measures inherit.","core_discovery":"The central discovery is the first relativistic MHD simulation of driven turbulence in which synchrotron cooling is included self-consistently as a momentum and energy sink, allowing the plasma thermodynamics to be regulated by the competition between turbulent heating and radiation. In quasi-steady state the total cooling power equals the energy injection rate, and the volume-averaged temperature falls as the inverse square root of the cooling efficiency, $\\langle\\Theta_t\\rangle_V \\propto \\eta_{\\rm syn}^{-1/2}$. At high cooling efficiency, the plasma bifurcates into a low-density, high-temperature phase and a high-density, low-temperature phase; the paper attributes this bifurcation to the synchrotron-cooling-induced thermal instability, in which overdense regions carry stronger magnetic fields, cool more rapidly, lose pressure, and draw in surrounding material until turbulent mixing arrests the runaway. The phase structure broadens the integrated synchrotron spectrum, depolarizes low-frequency emission along the line of sight because magnetic fields are tangled, and makes high-frequency emission come from rare, hot, strongly magnetized columns whose polarization approaches the theoretical maximum. Faraday rotation measures fluctuate in space and time, with the mean-field contribution comparable to the turbulent contribution and with fluctuations that weaken as the plasma becomes hotter. These synthetic diagnostics qualitatively match trends reported for blazars, pulsar wind nebulae, and repeating FRBs.","pith_inferences":["Editorial inference: replacing the Maxwell–Jüttner electron distribution with a nonthermal power-law tail would change the polarization ceiling and the temperature dependence of the rotation measure, so the qualitative trends are likely robust but the quantitative limits are not.","Editorial inference: the cold-dense phase is where the strongest magnetic fields concentrate, so if reconnection or current sheets live there, nonthermal particle acceleration would be spatially concentrated in exactly the regions that dominate high-frequency emission; kinetic simulations could test this.","Editorial inference: the rotation-measure scaling could be inverted as an observational tool; multi-frequency RM monitoring of a repeating FRB over days should map the cooling efficiency and temperature evolution of the screen rather than just its column density.","Editorial inference: adding synchrotron self-absorption or inverse-Compton losses would effectively renormalize the cooling efficiency, so the same two-phase structure should reappear in a broader class of radiatively cooled turbulent plasmas beyond the regimes simulated here."],"forward_implications":["Blazar flares that are faster and more strongly polarized at higher frequencies need no separate emission component; a single cooled turbulent region produces this frequency-dependent intermittency.","In pulsar wind nebulae, the observed rise of X-ray polarization degree with photon energy follows from line-of-sight integration through cooled turbulence, because high-energy emission is dominated by rare hot magnetized columns.","Faraday rotation measures around repeating fast radio bursts can fluctuate with Laplace-like distributions, and the paper's scaling ties the fluctuation amplitude to the mean electron temperature, so a measured RM constrains the temperature of the screen.","Because the magnetic power spectrum is Kolmogorov-like and nearly independent of cooling efficiency, the radiative predictions are stable to changes in magnetization and resolution at fixed cooling efficiency.","Stronger cooling drives the mean temperature down as $\\langle\\Theta_t\\rangle_V\\propto\\eta_{\\rm syn}^{-1/2}$, so systems with higher optical depth are naturally cooler and exhibit larger rotation-measure fluctuations."],"supporting_citations":[{"why":"supplies the theoretical Thomson-optical-depth scaling of mean temperature with cooling efficiency that the simulations reproduce.","marker":"Uzdensky 2018"},{"why":"introduces the synchrotron-cooling thermal instability that the paper invokes for the phase bifurcation.","marker":"Simon & Axford 1967"},{"why":"develops the hot-dilute and cold-dense equilibrium picture that the simulations are compared with.","marker":"Eilek & Caroff 1979"},{"why":"provides the nonlinear evolution of the instability used as the numerical benchmark in Appendix B.","marker":"Bodo et al. 1992"},{"why":"supplies the relativistic MHD source-term formulation for the cooling drag and energy loss.","marker":"Noble et al. 2009"},{"why":"provides the synchrotron emissivity formula for a Maxwell–Jüttner electron distribution used in the spectra.","marker":"Mahadevan et al. 1996"},{"why":"gives the emissivity and polarization implementation for thermal relativistic plasma used in the ray tracing.","marker":"Pandya et al. 2016"},{"why":"derives the Faraday rotation coefficient for relativistic thermal electrons entering the rotation measure.","marker":"Shcherbakov 2008"},{"why":"is the Vela nebula X-ray polarization observation used as a qualitative comparison target.","marker":"Xie et al. 2022"},{"why":"is the blazar multi-wavelength polarization variability observation the synthetic statistics are compared with.","marker":"Liodakis et al. 2022"}],"fun_headline_variants":["Synchrotron cooling splits relativistic plasma into two phases","Cooling drives phase separation in turbulent relativistic plasma","Relativistic turbulence bifurcates under synchrotron cooling","Synchrotron cooling triggers thermal instability in relativistic plasma","Turbulent plasma phase split matches blazar and FRB observations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The interpretation that the phase structure is caused by the classical thermal instability assumes that the instability criterion derived for a uniformly heated, quiescent plasma still applies when heating comes from patchy turbulent dissipation rather than a uniform source.","fun_headline_variants_meta":{"raw":{"variants":["Synchrotron cooling splits relativistic plasma into two phases","Cooling drives phase separation in turbulent relativistic plasma","Relativistic turbulence bifurcates under synchrotron cooling","Synchrotron cooling triggers thermal instability in relativistic plasma","Turbulent plasma phase split matches blazar and FRB observations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000659,"raw_usage":{"total_tokens":3039,"prompt_tokens":994,"completion_tokens":2045,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":610,"completion_tokens_details":{"reasoning_tokens":1962}},"tokens_in":610,"tokens_out":2045,"duration_ms":13833,"temperature":1.0,"reasoning_tokens":1962,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T18:08:10.838633+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the highest-cooling run at twice the linear resolution and check whether the two peaks in the lab-frame density PDF, near $\\sim\\rho_0/5$ and $\\sim 3\\rho_0$, persist with unchanged volume filling fractions; the paper's own appendix shows that the pure thermal instability is resolution-limited, so if the phase structure is numerically set rather than converged, the thermal-instability interpretation would not stand.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the theoretical Thomson-optical-depth scaling of mean temperature with cooling efficiency that the simulations reproduce."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"introduces the synchrotron-cooling thermal instability that the paper invokes for the phase bifurcation."},{"cited_title":"A., & Caroff, L","cited_arxiv_id":null,"evidence_quote":"develops the hot-dilute and cold-dense equilibrium picture that the simulations are compared with."},{"cited_title":"1996, ApJ, 465, 327, doi: 10.1086/177422","cited_arxiv_id":null,"evidence_quote":"provides the synchrotron emissivity formula for a Maxwell–Jüttner electron distribution used in the spectra."}],"review_version":1}