{"id":"11d8e39e-27ea-4dae-9cc0-8daae3e17c46","arxiv_id":"1908.05971","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Near-critical collapse with a scalar field and a Yang-Mills field shows the scalar field dominating at small scales, explained by a conjectured quasi-periodically self-similar critical solution with one unstable mode.","lead":"In simulations of gravitational collapse with two competing matter fields, the scalar field takes over and dominates at the smallest scales, regardless of the initial mix. The authors propose that a single quasi-periodically self-similar critical solution connects two known collapse behaviors, offering a template for how matter and gravitational waves might compete in more general collapse.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central scalar-dominance claim rests on unreported Floquet exponents; if either sign of Re lambda is wrong the universal conclusion fails.","rationale":"The paper has a real numerical discovery: in mixed scalar/YM collapse the near-critical evolution displays a crossover from YM-dominated to scalar-dominated scaling, with a shared accumulation point and a visible break in the mass/curvature exponent around |p-p*|~1e-10 for q=0.92. The single-field checks reproduce the known periods and exponents, and machine-precision fine-tuning to |p-p*|~1e-15 is a strong consistency test. My concern is not about the existence of the crossover but about the stronger universal statement that scalar dominance holds for sufficiently good fine-tuning for every nonzero scalar admixture. The q=0.95 run shows no crossover in the accessible range, and the only quantitative handle on the asymptotic regime is the pair of Floquet exponents after Eq. (44). Neither the boundary-value problem nor error bars are given, and the exactly periodic YM background used for the scalar mode is an idealized truncated-model solution rather than the actual asymptotically DSS YM critical solution (footnote 4). A positive test-field growth rate is necessary but not sufficient for the claimed nonlinear attractor; the QSS solution is explicitly conjectural. An independent spectral computation of the two exponents, with resolution control, directly tests the sign that the central claim depends on. This is the same weakness the reader identified; it warrants a conditional verdict but not rejection, because the observed crossover and the consistency of the single-field limits are solid evidence.","tokens_in":14312,"tokens_out":11821,"duration_ms":119252,"concrete_test":"Independently solve the Floquet eigenvalue problem for the scalar test field on the truncated YM DSS background and for the YM test field on the scalar DSS background, using a spectral method with Richardson extrapolation to estimate the error in Re lambda; check that Re(lambda_scalar-on-YM) is robustly positive (near 0.09) and Re(lambda_YM-on-scalar) is robustly negative (near -0.40). If either sign changes under resolution or under including the non-ideal YM self-coupling corrections at finite T, the universal dominance claim is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The universal statement that 'the scalar field always dominates on sufficiently small scales' is not directly observed for all mixtures: for q=0.95 the authors report a constant exponent 0.22 with YM dominance throughout the available fine-tuning range. The claim therefore rests on the test-field Floquet exponents reported after Eq. (44), in particular Re lambda = 0.09 for the scalar on the YM background and Re lambda = -0.40 for the YM on the scalar background. These values are stated without the boundary-value problem, without a specification of the exactly periodic YM background (which, as footnote 4 notes, is only a solution of the truncated model (37,38)), and without error estimates. A positive decoupled growth rate does not by itself establish that the fully nonlinear coupled evolution is attracted to the pure scalar critical solution; the QSS connection is explicitly conjectural and no construction is given. If either sign is wrong, the scalar-dominance conclusion fails.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper investigates type II critical collapse in spherical symmetry for a system containing two massless matter fields, a scalar field and a magnetic SU(2) Yang-Mills field, coupled only gravitationally. For each field separately the critical solution is discretely self-similar, with known periods and critical exponents. The authors evolve a two-parameter family of initial data (p,q) with a double-null code and adaptive mesh refinement down to |p-p*|≈10^-15. In the mixed cases q≈0.9–0.93 they find a single 'shared' critical solution in which the Yang-Mills field dominates at large scales and the scalar field takes over at small scales, with a corresponding break in the scaling exponents; for q=0.95 the available fine-tuning shows only Yang-Mills dominance. The scalar takeover is explained by test-field mode results reported after Eq. (44): the scalar grows on the Yang-Mills background while the Yang-Mills decays on the scalar background. On this basis the authors conjecture a quasi-periodically self-similar (QSS) solution connecting the two single-field critical solutions, with one unstable mode, and propose that it governs critical collapse for generic mixtures of the two fields.","tokens_in":14455,"tokens_out":8404,"duration_ms":81791,"significance":"This is a well-executed numerical study with a surprising central observation. The authors fine-tune down to |p-p*|≈10^-15, verify the single-field benchmarks (echo counts and exponents for the scalar and Yang-Mills critical solutions), and document a clean stress-energy switchover with a shared accumulation point in the mixed case q=0.92. The test-field mode analysis is a genuinely independent ingredient that, if correct, explains the scalar dominance in a simple way. The two proposed phase-space pictures (Figs. 12 and 13) are useful for framing the result and for guiding future work on competing massless degrees of freedom. If the missing eigenvalue documentation is supplied and the extrapolation to all q is either substantiated or made explicitly conjectural, the paper would be a valuable contribution to critical collapse.","major_comments":[{"comment":"The central claim that the scalar field always dominates on sufficiently small scales is load-bearing and is supported, beyond the observed q≈0.9–0.93 cases, by the two test-field eigenvalues reported in the paragraph after Eq. (44): λ=0.09+4.2i for the scalar on the Yang-Mills background and λ=-0.40+3.7i for the Yang-Mills on the scalar background. The manuscript does not describe the boundary-value problem used to compute these modes, the discretization, the boundary conditions at the regular centre and the past light cone, or any error estimates, and it does not specify whether the background is the numerical asymptotic Yang-Mills solution or the exactly DSS solution of the truncated model (37)-(38); footnote 4 states that an exact DSS Yang-Mills solution exists only in the truncated model. Because a change of sign of either real part would reverse the conclusion, these missing details are load-bearing and must be supplied, or the universal statement must be explicitly downgraded to a conjecture.","section":"Section VI, after Eq. (44)"},{"comment":"The universal statement is not directly observed for all mixtures, since for q=0.95 the authors report a constant exponent 0.22 with Yang-Mills dominance throughout the fine-tuning range down to |p-p*|≈10^-15. The claim that scalar dominance sets in 'for sufficiently good fine-tuning' is therefore an extrapolation based on the test-field eigenvalues and the conjectured QSS solution, rather than an observed result. Please either add an explicit test of the nonlinear evolution for a case with no observed break, or state clearly that the universal small-scale statement is a conjecture motivated by the test-field calculation.","section":"Section VI, q=0.95 paragraph"},{"comment":"The statement that the QSS solution 'has only one unstable mode, and so acts as the critical solution' is not established by the numerical evolutions presented, because no linear perturbation analysis of the QSS solution itself is given and the QSS solution is not constructed. The numerical results show a shared accumulation point and stress-energy switchover for an intermediate range of q, but the one-unstable-mode property and the resulting global attractor claim are part of the conjecture. Please mark this distinction explicitly in the abstract and conclusions, or provide a perturbation analysis of the QSS solution.","section":"Section VII, phase-space paragraph"}],"minor_comments":[{"comment":"The phrase 'at small q the YM field dominates' appears inconsistent with the definition q=0 (pure scalar) and q=1 (pure Yang-Mills) and with the reported q=0.95 behavior; please check whether 'small' and 'large' should be interchanged.","section":"Section VII, first paragraph of 'Our main evidence...'"},{"comment":"The abstract uses 'quasi-discretely self-similar (QSS)' while the conclusions define 'quasi-periodically self-similar (QSS)'; please use one consistent term throughout.","section":"Abstract and Section VII"},{"comment":"The axis label 'min/max(Ttilde1/2)' is ambiguous; write explicitly that the plotted curves are the minimum and maximum over x of \\tilde T_{(1)} and \\tilde T_{(2)}.","section":"Section VI, Fig. 10"},{"comment":"The quoted mixed exponents (0.25, 0.33, 0.37, 0.23, 0.27) are described as 'fitted by eye'; if these values are intended as quantitative results, describe the fitting procedure and give uncertainties.","section":"Section VI, paragraph after Fig. 8"}],"recommendation":"major_revision","confidential_remarks":"To the editor: The paper is interesting and likely publishable after revision. The main risk is the undocumented Floquet eigenvalue calculation, which carries the universal small-scale claim; I would like to see either a full description with convergence checks or a softening of the abstract and conclusions. There is also a small-q/large-q typo in Section VII that should be corrected before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things you should know. First, this paper reports a genuinely new numerical result: in spherical collapse with a scalar field plus Yang-Mills, the scalar takes over at small scales for intermediate mixtures, with a visible break in the mass/curvature scaling exponent. Second, the explanation—a conjectured quasi-periodically self-similar (QSS) solution connecting the YM critical solution at large scales to the Choptuik scalar solution at small scales—is honest about its status, but the linear-stability eigenvalues that support it are under-documented. That is where I would push back.\n\nWhat is new: this is the cleanest two-field, gravity-only-coupled critical collapse study I know. Kain's scalar triplet and Maliborski-Rinne's sphaleronic YM have direct couplings; here the two fields interact only through the Einstein equations. The crossover in exponents and the asymmetric perturbative stability are absent from those papers. The numerics are careful: fine-tuning to |p-p*| ~ 10^-15, consistent scaling breaks, stress-energy switchover, and validation against the single-field critical solutions.\n\nNow the soft spots. The Floquet exponents (lambda = 0.09 + 4.2i for the scalar on YM, -0.40 + 3.7i for YM on scalar) are stated without the boundary-value problem, without error estimates, and the YM background is only a solution of the truncated model (37,38) (see footnote 4). This matters because the universal statement \"scalar always dominates\" is not directly observed for q=0.95, where YM dominates throughout the accessible fine-tuning range. The asymptotic claim leans on the sign of Re(lambda), and that calculation is not reproducible from the text. The QSS solution itself is not constructed; it is inferred from intermediate-T behaviour. Mixed exponents are fitted by eye. Code and data are not public.\n\nNone of this breaks the central observation for q around 0.9-0.93, where the scalar takeover is directly seen. The paper is well written and the authors are clear about what is observed versus conjectured. The stress-test worry about the eigenvalues is real, but it mainly affects the universal conclusion, not the existence of the crossover.\n\nWho should read it: anyone working on critical collapse or on how matter and gravitational-wave degrees of freedom compete. It deserves peer review. I would send it out, but the referee should ask for the eigenvalue calculation to be described or the relevant code/tables made available. I would not block on the absence of a constructed QSS solution; that is a follow-up.","headline":"Solid numerical discovery of scalar-field takeover in two-field critical collapse, with an under-documented Floquet-eigenvalue argument behind the universal claim.","tokens_in":15010,"tokens_out":4903,"would_cite":true,"duration_ms":46851,"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":"For sufficiently fine-tuned collapse with two massless matter fields, the scalar field always wins at small scales, and the two fields share a single quasi-discretely self-similar critical solution.","keywords":["critical collapse","type II critical phenomena","discrete self-similarity","scalar field","Yang-Mills field","black hole threshold","quasi-discretely self-similar solution","gravitational wave toy model"],"falsifier":"Compute the dominant perturbation spectrum with an independent method and error estimates: if a scalar perturbation on the Yang-Mills critical background has $\\operatorname{Re}(\\lambda)<0$, or a Yang-Mills perturbation on the scalar background has $\\operatorname{Re}(\\lambda)>0$, the scalar-dominance conclusion collapses. A cheaper numerical test is to fine-tune mixed data beyond the current $|p-p_*|\\sim10^{-15}$ and check whether the mass-scaling exponent continues to break toward the scalar value; if the exponent stays at the Yang-Mills value at all accessible scales, the QSS conjecture is wrong.","tokens_in":14087,"feed_emoji":"🕳️","tokens_out":9255,"duration_ms":82318,"temperature":0.7,"pith_summary":"This paper asks what happens at the threshold of black hole formation when two independent massless matter fields—a scalar field and a Yang-Mills field—are coupled only through gravity. It claims that, for sufficiently good fine-tuning, the scalar field always dominates on sufficiently small scales, regardless of the initial mixture. To explain this, the authors conjecture a 'quasi-discretely self-similar' (QSS) solution that looks like the Yang-Mills critical solution at large scales and like the Choptuik scalar critical solution at small scales, with the scalar progressively taking over in between. Because this QSS solution would have only one unstable mode, it would act as the critical solution for any mixture of the two fields. The point of the toy model is to ask whether, in real axisymmetric collapse, gravitational waves might similarly lose a small-scale competition to matter.","feed_headline":"Scalar field always dominates in mixed critical collapse","feed_subtitle":"Two competing matter fields share one critical solution: Yang-Mills at large scales, scalar at small.","key_machinery":"The load-bearing object is the conjectured 'quasi-discretely self-similar' (QSS) solution: a nonlinear spacetime that is exactly discretely self-similar only in the two scale limits, with period $\\Delta_{\\mathrm{YM}}\\simeq0.6$ at large scales and $\\Delta_{\\mathrm{scal}}\\simeq3.44$ at small scales, and only quasiperiodic in between. The mechanism that selects the scalar field is the linear mode asymmetry computed on each single-field critical background: a scalar perturbation of the Yang-Mills solution grows like $e^{0.09T}$, while a Yang-Mills perturbation of the scalar solution decays like $e^{-0.40T}$. This one-growing-mode structure is what gives the QSS solution its single unstable direction and makes it the attractor for mixed initial data.","core_discovery":"On the paper's own terms, the central discovery is that the type II critical solution for the combined Einstein-scalar-Yang-Mills system is not a new two-field solution but a scale-dependent takeover. At large scales the spacetime approaches the known Yang-Mills DSS critical solution, with period $\\Delta_{\\mathrm{YM}}\\simeq0.6$; as $T\\to\\infty$ and scales shrink, the scalar field's stress-energy grows and the solution approaches the Choptuik scalar critical solution, with period $\\Delta_{\\mathrm{scal}}\\simeq3.44$. The evidence is an asymmetry in the linear mode problem: a scalar test field on the Yang-Mills critical background grows as $\\lambda=0.09+4.2i$, while a Yang-Mills test field on the scalar background decays as $\\lambda=-0.40+3.7i$. In mixed near-critical evolutions the two fields share one accumulation point $u_*$, and the mass scaling exponent breaks from a Yang-Mills-like value ($\\gamma\\simeq0.25$) at low fine-tuning to a larger scalar-like value (roughly $0.33$ to $0.37$) at high fine-tuning. The authors conjecture a one-parameter family of QSS solutions connecting the two single-field critical solutions, each with exactly one unstable mode.","pith_inferences":["If the QSS picture transfers to axisymmetric Einstein-matter collapse, matter may always dominate gravitational waves at the smallest scales, with a purely gravitational critical solution appearing only when matter is exactly absent; the paper leaves this as an open question.","The asymmetry between the two growth rates may be linked to the ratio of the two DSS periods: the faster-oscillating Yang-Mills solution, with smaller $\\Delta$, is the one that is destabilised. This could be tested by pairing other matter models with different periods.","A nonlinear boundary-value construction of the QSS solution itself would be a direct check: its eigenvalue spectrum should reproduce the one growing mode and predict the $q$-dependent effective critical exponents that the paper currently fits by eye.","The same numerical setup could measure whether the takeover scale shifts with the initial scalar amplitude in the way the one-growing-mode picture predicts, which would be a sharper test of the conjectured family than the present data."],"forward_implications":["Any one-parameter fine-tuning of mixed scalar/Yang-Mills data that contains a nonzero scalar component should produce the same QSS critical solution, so there is no separate two-field critical solution to find.","The pure Yang-Mills critical solution acquires a second unstable mode once the scalar field is allowed, while the pure scalar critical solution keeps one, so the scalar end is the generic small-scale endpoint.","The effective black-hole-mass exponent in mixed evolutions is not a constant but drifts from the Yang-Mills value toward the scalar value as fine-tuning improves, with the break location depending on the initial scalar amplitude.","In the QSS picture, curvature and mass scaling oscillations should become quasiperiodic in the intermediate regime, showing broadened frequency peaks rather than the sharp peaks of either single-field DSS solution."],"supporting_citations":[{"why":"Supplies the pure scalar-field critical solution, its DSS period and critical exponent against which the mixed results are compared.","marker":"[1]"},{"why":"Provides the linear perturbation and scaling analysis of the scalar critical solution, including the periodic wiggle used to validate the numerics.","marker":"[18]"},{"why":"Supplies the pure Yang-Mills critical solution and critical exponent used as the other endpoint of the QSS family.","marker":"[19]"},{"why":"Gives the Yang-Mills critical solution and exponent from a DSS ansatz, the large-scale limit of the conjectured QSS solution.","marker":"[20]"},{"why":"Supplies the double-null numerical method in spherical symmetry that this paper adapts to the two-field system.","marker":"[26]"},{"why":"Describes the closest previous two-field system, Yang-Mills plus triplet scalar, and provides the contrast for the mixed critical behaviour found here.","marker":"[21]"},{"why":"Studies the related two-degree-of-freedom Yang-Mills system and is interpreted as showing one field, the magnetic amplitude, dominating the critical dynamics.","marker":"[22]"},{"why":"Provides the analytic description of the periodic scaling wiggle that the paper uses to confirm discrete self-similarity in its evolutions.","marker":"[28]"}],"fun_headline_variants":["Scalar wins small scales in two-field collapse","Mixed collapse: one critical solution, two regimes","Choptuik wins small scales in mixed critical collapse","Quasi-self-similar solution unites two critical collapses"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the reported linear-mode rates are correct—specifically that the scalar mode on the Yang-Mills background grows with $\\operatorname{Re}\\lambda=0.09$ while the Yang-Mills mode on the scalar background decays with $\\operatorname{Re}\\lambda=-0.40$—but the paper gives these numbers without describing the boundary-value calculation or providing error estimates.","fun_headline_variants_meta":{"raw":{"variants":["Scalar wins small scales in two-field collapse","Mixed collapse: one critical solution, two regimes","Choptuik wins small scales in mixed critical collapse","Quasi-self-similar solution unites two critical collapses"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000805,"raw_usage":{"total_tokens":3585,"prompt_tokens":1044,"completion_tokens":2541,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":660,"completion_tokens_details":{"reasoning_tokens":2478}},"tokens_in":660,"tokens_out":2541,"duration_ms":17699,"temperature":1.0,"reasoning_tokens":2478,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:58:44.181510+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the dominant perturbation spectrum with an independent method and error estimates: if a scalar perturbation on the Yang-Mills critical background has $\\operatorname{Re}(\\lambda)<0$, or a Yang-Mills perturbation on the scalar background has $\\operatorname{Re}(\\lambda)>0$, the scalar-dominance conclusion collapses. A cheaper numerical test is to fine-tune mixed data beyond the current $|p-p_*|\\sim10^{-15}$ and check whether the mass-scaling exponent continues to break toward the scalar value; if the exponent stays at the Yang-Mills value at all accessible scales, the QSS conjecture is wrong.","supporting_citations":[{"cited_title":"Alcubierre, G","cited_arxiv_id":null,"evidence_quote":"Provides the linear perturbation and scaling analysis of the scalar critical solution, including the periodic wiggle used to validate the numerics."},{"cited_title":"Garfinkle and G.C","cited_arxiv_id":null,"evidence_quote":"Supplies the pure Yang-Mills critical solution and critical exponent used as the other endpoint of the QSS family."},{"cited_title":"We represent our fields on a grid at fixed values of v, and numerically advance in the retarded time u","cited_arxiv_id":null,"evidence_quote":"Supplies the double-null numerical method in spherical symmetry that this paper adapts to the two-field system."},{"cited_title":"Sorkin, Class","cited_arxiv_id":null,"evidence_quote":"Describes the closest previous two-field system, Yang-Mills plus triplet scalar, and provides the contrast for the mixed critical behaviour found here."},{"cited_title":"Maliborski and O","cited_arxiv_id":null,"evidence_quote":"Studies the related two-degree-of-freedom Yang-Mills system and is interpreted as showing one field, the magnetic amplitude, dominating the critical dynamics."},{"cited_title":"Bartnik and J","cited_arxiv_id":null,"evidence_quote":"Provides the analytic description of the periodic scaling wiggle that the paper uses to confirm discrete self-similarity in its evolutions."}],"review_version":1}