{"id":"8323b270-90ac-47f1-8005-2a8bb34e30aa","arxiv_id":"2608.12530","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Kelvin-Helmholtz-unstable shear flows generate large-scale, quasi-periodically reversing magnetic fields through jets that realize Yoshizawa's cross-helicity (Upsilon) dynamo effect.","lead":"This paper shows, with analytic theory and simulations using up to 4096 by 4096 by 8192 grid points, that turbulence from an unstable shear flow organizes into large-scale jets, and those jets stretch magnetic fluctuations to build a large-scale, oscillating magnetic field. It is a concrete realization of the long-postulated mean-vorticity, or Upsilon, dynamo, with applications to neutron star mergers and the Sun.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (S29), the self-consistent evolution equation for Upsilon that makes the dynamo exponential, is explicitly deferred to a 'forthcoming publication'; the central growth-rate claim is therefore uncheckable as written.","rationale":"The reader's conditional verdict identifies exactly the load-bearing weak point: the analytic derivation of Eq. (S29) is explicitly deferred, while the quantitative claim of exponential growth driven by the Upsilon effect depends on it. My reading of the full manuscript confirms this. The simulation campaign is substantial and the auxiliary diagnostics—especially Extended Data Fig. 3, which shows large-scale jets transferring energy to the mean field—provide genuine support for the qualitative mechanism. Methods VI is internally consistent as far as it goes: the quasilinear EMF coefficients in Eqs. (S11)-(S12) follow from the assumed closure with straining times tau_ZF and tau_ZM. However, those coefficients only parameterize the mean EMF; they do not provide the evolution equation for Upsilon. Without Eq. (S29), the step from 'the EMF is dominated by Upsilon in the regression' to 'the dynamo grows exponentially because of Upsilon' is not demonstrated. A missing derivation is not a falsification, and the numerical evidence may be correct, so a conditional verdict with a request for the derivation is appropriate. I do not see an internal inconsistency strong enough to reject the paper; nor is the central claim sufficiently established from the text alone to accept unconditionally. The proposed concrete test, an independent derivation of Eq. (S29) from the stated quasilinear system and a comparison of the resulting growth rate with the simulation value, would settle whether the concern lands.","tokens_in":18788,"tokens_out":7077,"duration_ms":72219,"concrete_test":"Independently derive Eq. (S29) from the quasilinear equations (S7)-(S10) and definitions (S12a)-(S12c) in Methods VI, without invoking any future publication, and compute the growth rate of the coupled (S28)+(S29) system; if the derivation requires an empirical closure or yields a rate inconsistent with the reported 0.04 U0/a, the Upsilon-attribution claim is not supported by the manuscript.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section XII states Eq. (S28) as the load-bearing induction equation and then asserts that Upsilon evolves as (∂t + c)Upsilon = f(Bx), with c constant and f a linear operator, so that (S28) and (S29) yield exponential growth. The very next sentence defers the derivation: 'we will report the details of the calculations in a forthcoming publication.' Neither c nor f(Bx) is given anywhere in Methods VI; the quasilinear calculation there derives the EMF coefficients (S12a)-(S12c) but stops short of closing the cross-helicity dynamics. The paper also claims Upsilon grows exponentially 'as we have observed in simulations' but provides no time series of Upsilon or growth-rate comparison. Consequently, the stated growth rate (0.04 U0/a) and the BNS e-folding estimate of ~8 µs are not checkable from the manuscript. The simulation data are real evidence for amplification, but the central claim that this is an ab initio Upsilon dynamo—rather than a large-scale shear-driven instability with a different cause—quantitatively depends on an absent derivation. This is a missing-support flag raised by the paper itself, not an outside objection.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper claims an ab initio realization of Yoshizawa's 1990 mean-vorticity (Upsilon) dynamo in driven, Kelvin-Helmholtz-unstable shear flows. The authors develop a quasilinear mean-field model in which the turbulent EMF is written as alpha B - beta curl B + Upsilon curl U, and they report three-dimensional incompressible MHD simulations (two independent codes, periodic and non-periodic domains, up to 4096x4096x8192 grid points, parameter scans in Re, Rm, and Pm) showing amplification of the mean field by three orders of magnitude, quasi-periodic polarity reversals, and generation even from infinitesimal initial magnetic fluctuations. They attribute the dynamo to large-scale, x-invariant 3D jets that are exact nonlinear MHD solutions, and they propose applications to binary neutron star mergers, the solar dynamo, galaxy clusters, and laboratory experiments. The central analytic statement is Eq. (S28), where the mean field evolves under a Upsilon-driven source, together with Eq. (S29), a stated self-consistent evolution equation for Upsilon that would produce exponential growth; however, the derivation of Eq. (S29) is explicitly deferred to a forthcoming publication.","tokens_in":19025,"tokens_out":4838,"duration_ms":47111,"significance":"If the central claim is substantiated, this would be a notable step: the first systematic numerical demonstration of a large-scale dynamo driven by the mean-vorticity effect, generated through self-organized jets rather than through the traditional alpha or Omega effects. The numerical campaign is extensive and well documented: two independent pseudospectral solvers, benchmarked against each other, a very large highest-resolution run, a wide parameter scan, and public data on Zenodo. The spectral transfer-function analysis in Extended Data Figs. 2 and 3 is a genuine independent diagnostic showing that large-scale jets transfer energy to the mean field and that the dynamo is not produced by an inverse cascade. If the mechanism survives scrutiny, it would have plausible implications for BNS mergers, the solar meridional-flow dynamo, and laboratory dynamo experiments, and it would help motivate cross-helicity-based mean-field closures. The significance is therefore high, but it is conditional on closing the gap between the numerical evidence and the analytic identification of the growth mechanism.","major_comments":[{"comment":"The load-bearing step of the analytic theory is deferred. Equations (S28) and (S29) are stated to be linearly coupled and to produce exponential growth, but neither the constant c nor the linear operator f(Bx) is specified, and the derivation is explicitly postponed to a forthcoming publication. The stated growth rate 0.04 U0/a and the BNS e-folding time of ~8 microseconds are therefore not checkable from the manuscript. This is a self-identified missing proof in the paper's own text, and it directly affects the central claim that the dynamo is an ab initio Upsilon dynamo rather than a shear-driven instability with a different cause.","section":"Methods Sec. XII, Eq. (S29)"},{"comment":"The Upsilon-dominance diagnostic is circular as presented. The coefficients alpha, beta, and Upsilon are recovered by fitting the assumed ansatz E = alpha B - beta curl B + Upsilon curl U to the simulated EMF using spatial and temporal regression, and the same fitted model is then used to conclude that Upsilon drives the dynamo. This confirms the ansatz only if the ansatz is already complete. The independent transfer-function diagnostic in Extended Data Fig. 3 supports the role of large-scale jets, but it measures a different quantity, the nonlinear transfer to the mean field, and does not by itself identify the EMF coefficient Upsilon. Please separate the regression-based model recovery from an out-of-sample or otherwise independent verification of the closure ansatz.","section":"Main text Fig. 2c and Methods Sec. VI"},{"comment":"The quasilinear closure depends on two free coherent-straining times, tau_ZF and tau_ZM, and the equality tau_ZF = tau_ZM is assumed without a derived basis or a sensitivity study. Since the paper motivates itself by criticizing traditional mean-field dynamos for containing parameters that are not justified from first principles, the present analytic model does not yet deliver a first-principles prediction of the growth rate; it provides a closure-based interpretation of the simulations. Please either derive or bound tau_ZF and tau_ZM, or state explicitly that the analytic component is interpretive and that the quantitative growth-rate claim rests on simulation measurement.","section":"Methods Sec. VI, Eqs. (S12a)-(S12c)"}],"minor_comments":[{"comment":"There is a typo: 'Following thar order' should be 'Following that order.'","section":"Main text, paragraph after Fig. 3a"},{"comment":"The caption says 'rms stands for room mean square'; this should be 'root mean square.'","section":"Fig. 4 caption"},{"comment":"The text describes the reversals as 'quasi-periodic,' but the intervals between reversals in Fig. 1a (around t = 3000, 5300, 6000, 6350) do not appear periodic; 'quasi-cyclic' or 'intermittent reversals' would be more accurate unless a periodicity analysis is provided.","section":"Main text, Fig. 1a"},{"comment":"The resolution range from 32^2 x 512 to 4096^2 x 8192 is stated and the two codes are benchmarked, but no explicit resolution-convergence study of the mean-field growth rate or saturation level is reported. A brief convergence statement would strengthen the quantitative claims.","section":"Methods Sec. IV"},{"comment":"The double shear-layer profile in Eq. (S3) is written with a trailing '1' after the two tanh terms; please clarify whether this is a normalization constant and define the domain extent consistently with the figure captions.","section":"Methods Sec. II, Eq. (S3)"}],"recommendation":"major_revision","confidential_remarks":"The numerical work appears substantial and the phenomenology of jet-driven mean-field generation is likely real. The main obstacle is not the simulations but the analytic identification: the central evolution equation for Upsilon is deferred, and the regression-based diagnostic is partly circular. The authors should be asked either to provide the full derivation of Eq. (S29), including c and f(Bx), and a simulated time series of Upsilon with a growth-rate comparison, or to explicitly weaken the claim from 'ab initio Upsilon dynamo' to 'large-scale jet-driven mean-field generation consistent with the Upsilon effect.' Given the paper's own acknowledgement of the omitted derivation, this is a fixable major revision rather than a rejection. The astrophysical extrapolation to BNS mergers is reasonable as a prediction but should be clearly separated from the verified numerical result."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a serious paper that likely delivers the first ab initio realization of Yoshizawa's mean-vorticity (Upsilon) dynamo. The simulation campaign is genuinely impressive: two independent pseudospectral codes, periodic and non-periodic domains, up to 4096x4096x8192 grid points, parameter sweeps in Re, Rm, and Pm, and a transfer-function analysis showing that x-invariant large-scale jets with kx=0, ky~0.2 are the dominant energy source for the mean field. Data and code are on Zenodo. That is real evidence, and identifying the 3D jets as the physical agent behind the previously postulated Upsilon effect is new.\n\nThe soft spots are in the theory, not the numerics. The load-bearing claim for exponential growth is Eq. (S29), the self-consistent evolution equation for Upsilon. The paper states that (dt+c)Upsilon = f(Bx) with c constant and f linear, says this closed system yields exponential solutions, and then defers the derivation to a forthcoming publication. Neither c nor f is given; no time series of Upsilon is shown. As a result, the stated growth rate 0.04 U0/a and the BNS e-folding estimate of ~8 microseconds are not checkable from the manuscript. The EMF coefficient recovery by regression is also partly circular, since the fitted Upsilon-dominated model is then used to confirm Upsilon dominance. That circularity is partially broken by the transfer-function analysis, which is independent. The \"topologically protected\" claim for the jets is asserted rather than proven; the vanishing pressure-gradient argument is a plausibility statement, not a topology proof. None of these flaws, taken together, falsifies the qualitative phenomenon: the simulations robustly show large-scale, reversing mean fields generated from small seed fields in forced KH turbulence. But the quantitative mechanism attribution rests on an absent derivation.\n\nThe applications sections are speculative but clearly framed as such; the MDE comparison is suggestive. The self-flagging of the deferred derivation is honest, and the paper is coherent on its own terms.\n\nMy recommendation: send it to peer review. A serious referee can check the numerics, ask for the Upsilon derivation or a softening of the growth-rate claims, and probe the regression validation. The core result deserves refereeing; the gaps are large but repairable.","headline":"Strong numerical evidence for a jet-driven Upsilon dynamo, but the exponential-growth theory is explicitly deferred, so the mechanism attribution runs ahead of the proof.","tokens_in":759,"tokens_out":1034,"would_cite":true,"duration_ms":36279,"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":"This paper claims that externally sustained shear turbulence spontaneously forms large-scale jets that drive a mean-vorticity dynamo, producing quasi-periodically reversing large-scale magnetic fields from near-zero initial fields.","keywords":["dynamo theory","mean-field electrodynamics","mean-vorticity effect","Upsilon dynamo","Kelvin-Helmholtz instability","shear-flow turbulence","large-scale magnetic fields","cross-helicity"],"falsifier":"Run the same externally forced Kelvin-Helmholtz setup with the $k_x=0, k_y\\neq 0$ jet modes artificially damped, or with the seed vertical fluctuation removed, and check whether the large-scale mean field still grows; if it does, the Upsilon mechanism is not the sole driver. Independently, a closure-free measurement of the time-resolved relation between $B_x$ and the cross-helicity coefficient Upsilon, fitting the coupled evolution equations directly to simulation data, would settle whether the exponential growth is real or a regression artifact.","tokens_in":2129,"feed_emoji":"🧲","tokens_out":5832,"duration_ms":95880,"temperature":0.7,"pith_summary":"The paper tries to establish that a large-scale magnetic dynamo can arise from scratch in an externally maintained, Kelvin-Helmholtz-unstable shear flow, with no initial large-scale field and no reliance on the classical alpha-effect. The driver is a different mean-field term, the mean-vorticity (Upsilon) effect postulated in 1990, in which turbulent cross-helicity converts the curl of the mean flow into mean magnetic field. The mechanism works through large-scale three-dimensional jets that form spontaneously in the turbulence, are exact nonlinear solutions of the ideal MHD equations, and survive over a wide range of Reynolds numbers. In simulations up to 4096 x 4096 x 8192 grid points, the mean field amplifies by three orders of magnitude and reverses polarity quasi-periodically. If correct, this is the first ab initio realization of the Upsilon dynamo and would change how shear-dominated systems from neutron-star mergers to the Sun are modeled.","feed_headline":"Shear-flow jets drive a new large-scale dynamo","feed_subtitle":"Simulations show reversing magnetic fields grow via the 1990 mean-vorticity effect, no initial field needed.","key_machinery":"The central object is the Upsilon coefficient, a turbulence closure parameter measuring turbulent cross-helicity $\\langle \\mathbf{u}'\\cdot\\mathbf{b}'\\rangle$. In the generalized mean-field ansatz $\\mathcal{E} = \\alpha B - \\beta \\nabla\\times B + \\Upsilon \\nabla\\times U$, Upsilon couples the mean vorticity of the shear flow to the mean magnetic field. The mechanical carriers are large-scale three-dimensional jets: velocity fluctuations directed along x, independent of x, varying in y and z, produced by stretching of Kelvin-Helmholtz-driven seed perturbations by the mean flow. They are exact nonlinear solutions of the ideal MHD equations with vanishing pressure and magnetic-pressure gradients, which is why they persist. The jets stretch $b_z$ fluctuations to produce a mean $B_x$ reversed across the shear layer; the load-bearing equation is $[\\partial_t - (\\eta+\\beta)\\nabla^2]B_x = \\mathbf{e}_x\\cdot[\\nabla\\times(\\Upsilon\\nabla\\times U)]$, with a companion evolution equation for Upsilon supplying exponential growth.","core_discovery":"The paper claims that the turbulent electromotive force in a driven shear flow contains a term $\\mathcal{E} = \\alpha B - \\beta \\nabla\\times B + \\Upsilon \\nabla\\times U$, and that the Upsilon term, not $\\alpha$ or $\\Omega$, generates the x-directed, z-reversed mean magnetic field. The source is the large-scale x-invariant jets $u_x(k_x=0,k_y)$, which stretch the magnetic fluctuation $b_z$; the resulting mean-field source is $B_x \\propto \\mathbf{e}_x\\cdot[\\nabla\\times(\\Upsilon\\nabla\\times U)]$. The jets form when the mean flow stretches a seed vertical fluctuation generated by the Kelvin-Helmholtz instability, and because they are exact nonlinear solutions with zero pressure-gradient forces, they resist turbulent and Lorentz destruction. The mean field grows exponentially because Upsilon itself grows with the fluctuations; the paper derives a coupled evolution equation for Upsilon, though the detailed derivation is deferred to a forthcoming publication. Simulations show quasi-periodic polarity reversals, dominance of cross-helicity over kinetic helicity by two orders of magnitude, and robustness to domain periodicity, initial-field presence, and magnetic Prandtl number above a threshold.","pith_inferences":["Editorial extension: if the Upsilon dynamo is confirmed, Alfvenization--the alignment of velocity and magnetic fluctuations--turns from a dynamo suppressor into the engine itself, because cross-helicity rather than kinetic helicity becomes the fuel; this could explain why strongly aligned, sheared astrophysical flows still produce large-scale fields.","Editorial extension: the same mechanism should appear in other vortex-dominated shear flows, such as accretion disks, galaxy-cluster merger shocks, and planetary atmospheres, whenever a large-scale vortical flow is externally sustained for longer than the dynamo growth time.","Editorial extension: a direct testable extension is to measure the Upsilon coefficient in liquid-sodium spherical Couette experiments with strong imposed shear, checking whether the turbulent electromotive force aligns with curl U rather than with B, as the paper's comparison with one laboratory experiment suggests.","Editorial extension: the deferred derivation of the Upsilon evolution equation is the logical place to stress-test the claim; until it appears, an independent, closure-free diagnostic of the time-resolved relation between $B_x$ and cross-helicity would settle whether the exponential growth is robust."],"forward_implications":["In binary neutron star mergers, the reported growth rate of about $0.04\\,U_0/a$ gives an e-folding time near 8 microseconds, so the mechanism could build $10^{16}$--$10^{17}$ G fields within the millisecond merger and shift gravitational-wave frequencies by 200--300 Hz, in reach of current and next-generation observatories.","In the Sun, the azimuthal vorticity of the meridional circulation would drive the poloidal field through the Upsilon effect, offering a parameter-light alternative to tuned alpha-profiles in flux-transport dynamo models.","The relaxed state of this dynamo has B nearly parallel or antiparallel to U, with B and curl B orthogonal, rather than the force-free state of alpha-dynamos; the sign of Upsilon sets the polarity, naturally producing cyclic reversals.","Because the jets are exact nonlinear MHD solutions and the mechanism needs no initial large-scale field, it should operate in any externally maintained shear layer above the Kelvin-Helmholtz threshold, regardless of boundary conditions and over a broad range of Reynolds numbers.","The dynamo is not driven by an inverse cascade: the energy flux is forward, from large scales to small scales, while the large-scale jets feed energy directly into the mean field."],"supporting_citations":[{"why":"Supplies the 1990 mean-vorticity-effect postulate that the paper confirms and generalizes.","marker":"4"},{"why":"Provides the laboratory turbulent-EMF measurement with field and EMF nearly orthogonal, which the paper interprets as Upsilon dominance.","marker":"37"},{"why":"Documents the large-scale radial vorticity in the same liquid-sodium experiment, the geometric ingredient the Upsilon effect needs.","marker":"38"},{"why":"Gives the binary-neutron-star merger shear-layer width and Kelvin-Helmholtz relevance used for the microsecond e-folding estimate.","marker":"5"},{"why":"Supplies the Newtonian Kelvin-Helmholtz treatment of neutron-star mergers and the estimate of ultrastrong amplified fields.","marker":"7"},{"why":"Provides the double-shear-layer benchmark geometry used in the highest-resolution simulation.","marker":"35"},{"why":"Supplies the external mean-flow forcing method that maintains the shear flow against turbulent depletion.","marker":"30"},{"why":"Provides the exact-solution framework for ideal MHD nonlinearities against which the jets are checked.","marker":"33"},{"why":"Defines the standard alpha-effect mean-field theory that the paper's generalized EMF ansatz extends.","marker":"24"}],"fun_headline_variants":["Shear-flow jets drive reversing magnetic fields from scratch","Mean-vorticity effect powers jet-driven dynamo without seed","Jets from shear flows switch on large-scale dynamo ab initio","Jet-driven dynamo creates reversing fields via mean vorticity"],"cache_read_input_tokens":21632,"weakest_assumption_plain":"The claim that the Upsilon term, rather than an artifact of the turbulence model, is what drives the growth depends on a closure with two unmeasured straining-time parameters, on assuming the three-term EMF formula is complete, and on a derivation of the Upsilon growth equation that is promised but not yet published.","fun_headline_variants_meta":{"raw":{"variants":["Shear-flow jets drive reversing magnetic fields from scratch","Mean-vorticity effect powers jet-driven dynamo without seed","Jets from shear flows switch on large-scale dynamo ab initio","Jet-driven dynamo creates reversing fields via mean vorticity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000705,"raw_usage":{"total_tokens":3243,"prompt_tokens":1071,"completion_tokens":2172,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":687,"completion_tokens_details":{"reasoning_tokens":2102}},"tokens_in":687,"tokens_out":2172,"duration_ms":14551,"temperature":1.0,"reasoning_tokens":2102,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:06:40.697792+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same externally forced Kelvin-Helmholtz setup with the $k_x=0, k_y\\neq 0$ jet modes artificially damped, or with the seed vertical fluctuation removed, and check whether the large-scale mean field still grows; if it does, the Upsilon mechanism is not the sole driver. Independently, a closure-free measurement of the time-resolved relation between $B_x$ and the cross-helicity coefficient Upsilon, fitting the coupled evolution equations directly to simulation data, would settle whether the exponential growth is real or a regression artifact.","supporting_citations":[],"review_version":1}