{"id":"f36d0e36-259c-4ef3-b445-2ec395e79929","arxiv_id":"2608.10948","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A one-dimensional hybrid contact process shows a continuous absorbing-to-active phase transition with a critical exponent near 0.79, suggesting a new universality class.","lead":"This paper studies a one-dimensional lattice where particles spread and die through both quantum and classical processes, and finds two types of transitions between an empty absorbing state and a populated active state. It matters because the continuous transition appears to have a new, nonclassical critical exponent that quantum simulators could test.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed β*≈0.79 outside the directed-percolation class rests on a CAM fit to only five to six cluster sizes (L=5–10) with no error bars; the paper's own caveat that larger clusters are needed undermines the central claim.","rationale":"The paper's qualitative results—the MF phase diagram with absorbing, active, and bistable regions, the saddle-node bifurcation indicating a discontinuous transition, and the preservation of a continuous absorbing-to-active transition in CMF—are coherent and clearly presented. The Liouvillian spectrum analysis is also consistent with critical slowing down, though not conclusive. The load-bearing weakness is exclusively the extraction of β* in Sec. VI. The CAM method is a legitimate extrapolation tool, but its validity here rests on two unverified assumptions: (i) that L=5–10 lies in the asymptotic scaling window described by Eq. (23), and (ii) that the exponential form for κ_c(L) is the correct finite-size convergence. With only five or six data points and no error estimates, the reported difference between β*=0.7962 and 0.7865 is not meaningful, and the claim of a new universality class (i.e., a value far from the directed-percolation exponent 0.276) is not robust. The paper's own caveat in Sec. VI—that larger clusters are required—directly undercuts the headline result as it currently stands. A concrete computational test with clusters up to L=14, or a bootstrap uncertainty analysis, would settle whether the exponent stabilizes. For these reasons, the conditional verdict is appropriate: the paper should be published only after the exponent claim is either confirmed with larger clusters or independent methods, or softened to a preliminary estimate. The reader's weakest-assumption analysis correctly identified the same concern, so no adjustment to the verdict is needed.","tokens_in":13686,"tokens_out":5139,"duration_ms":40620,"concrete_test":"Compute CMF data for L=8,10,12,14 (and L=9,11,13 if parity matters), re-fit Eq. (23) to the subset L=8–14, and compare the resulting β* and κ* with the values obtained from L=5–10. If β* changes by more than 0.05 or the residuals worsen, the claimed exponent is an artifact of the fitting range. Additionally, bootstrap the fit (e.g., leave-one-out or synthetic noise) to obtain a confidence interval on β*; if the interval includes the directed-percolation value 0.276 or the mean-field value 1, the universality-class claim is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the HCP continuous transition defines a new universality class hinges entirely on the coherent anomaly method (CAM) analysis in Sec. VI. Equation (23) postulates the asymptotic scaling n0(L) ∼ C0[(κ*−κ_c(L))/κ*]^{-(β*−β_mf)} with β_mf fixed to the classical value 1, but the paper provides no evidence that L=5–10 lies in the asymptotic regime. Only five data points (L=5–10) are fitted for Ω/Γ=0.1 and four for Ω/Γ=0.4, and no residuals or confidence intervals are reported; the two extracted values β*=0.7962 and 0.7865 are therefore statistically indistinguishable from each other and from a range of other exponents. The companion exponential ansatz for κ_c(L) (Fig. 5(c)) is equally ad hoc: no physical derivation is given, and the authors themselves note that the QCP counterpart exhibits an even-odd size effect, which would invalidate a smooth exponential fit. Because the CAM exponent is highly sensitive to the extrapolated κ_c, a small bias in the κ_c(L) fitting form—e.g., if the convergence is actually power-law—could shift β* substantially. The paper explicitly concedes 'larger sizes are required to obtain a more precise β*', yet the abstract and conclusion present β*≈0.79 as the central result. Without a demonstrated plateau in β* as L grows, or an independent calculation (e.g., tensor-network or Monte Carlo simulation of the actual 1D HCP), the existence of a new universality class is not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the one-dimensional hybrid contact process (HCP), where coherent (Ω) and incoherent (κ) coagulation and branching coexist with local decay (Γ). It derives single-site mean-field (MF) Bloch equations, identifies stable absorbing, active, and bistable steady states, and shows that the absorbing-to-active transition along κ is continuous while the absorbing-to-bistable transition can be discontinuous via a saddle-node bifurcation. Finite-size Liouvillian spectra show a decreasing gap and a nonmonotonic κ dependence. The authors then use cluster mean-field (CMF) theory with self-consistent effective fields to compute critical points κ_c(L) for cluster sizes L=5–10, and the coherent anomaly method (CAM) to extract thermodynamic-limit critical points κ*/Γ≈0.988 and exponents β*≈0.7962 and 0.7865 for Ω/Γ=0.1 and 0.4. The paper concludes that the continuous transition is non-DP and likely defines a new universality class.","tokens_in":14082,"tokens_out":6740,"duration_ms":57411,"significance":"The strength of the paper is the transparent MF and CMF construction: the Bloch equations, stability analysis, phase diagram, Liouvillian spectrum, and self-consistent effective-field CMF are all clearly laid out and internally consistent. If the CAM result β*≈0.79 were established, it would be a significant finding because it would place the HCP continuous absorbing transition outside the directed-percolation universality class, with implications for how quantum coherence modifies absorbing-state criticality. However, the central exponent claim is not yet supported at the required level: it relies on five or six CMF cluster sizes, an ad hoc exponential-size scaling, no uncertainty quantification, and no independent check. The paper's own caveat that larger clusters are needed is appropriate, but it is in tension with the abstract and conclusion presenting β*≈0.79 as the main quantitative result.","major_comments":[{"comment":"The central claim that β*≈0.79 defines a new universality class rests entirely on a three-parameter fit of n0(L) to Eq. (23) using only L=5,...,10 (and L=6,...,10 for Ω/Γ=0.4). No residuals, confidence intervals, or χ² values are reported, and the fit is not tested for stability under removal of the smallest or largest L. With this number of points, the extracted β* is statistically indistinguishable from a wide range of values, and the statement 'the true critical exponents are estimated' overstates what the data can support. The authors themselves note that larger clusters are required; the abstract and conclusion should reflect this limitation, or the fits should be augmented with error bars and a convergence test.","section":"Section VI, Eq. (23), Fig. 6"},{"comment":"The exponential finite-size scaling κ_c(L)/Γ = κ_c/Γ − A exp(−bL) is assumed with no derivation or empirical justification. If the convergence is actually power-law, or if there is an even-odd size effect as the authors report for the QCP, the extrapolated κ_c will be biased. Since β* in Eq. (23) is very sensitive to the value of κ*, a small bias in κ_c can produce a large shift in β*. The agreement between the direct extrapolation and the CAM critical point is not an independent confirmation because both analyses use the same CMF data and related assumptions.","section":"Section V, Fig. 5(c)"},{"comment":"The analysis assumes that the CMF order parameter has the classical exponent β_mf=1 and that the coherent anomaly scaling form Eq. (23), taken from previous work, applies to this model. The paper gives no evidence that L=5–10 is in the asymptotic regime for either the CMF critical points or the amplitude n0(L). This is a load-bearing assumption: if Eq. (23) is not the correct scaling form for the HCP, the extracted β* is not the true critical exponent. The authors should either provide a numerical check of the scaling collapse or present the CAM result as a preliminary estimate only.","section":"Section VI, Eq. (22)"}],"minor_comments":[{"comment":"The word 'self-desctruction' should be 'self-destruction'.","section":"Section II"},{"comment":"The text contains 'Lioullian' in several places; this should be 'Liouvillian'.","section":"Section IV"},{"comment":"The phrase 'lower brach' should be 'lower branch'.","section":"Section III D"},{"comment":"The notation 2Ω⟨σx⟩⟨σx⟩ in the second Bloch equation is confusing; please write 2Ω(⟨σx⟩)^2.","section":"Eq. (7)"},{"comment":"The x-axis is labeled 1/L, but the exponential fit is performed in L; please clarify in the caption that the extrapolation to 1/L=0 is made after fitting κ_c(L) as a function of L.","section":"Fig. 5(c)"},{"comment":"The fit parameters A, b, C0 and the values of n0(L) and κ_c(L) should be provided in a table or in supplementary material so that the CAM extrapolation can be reproduced and assessed.","section":"Section VI"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely to be of interest to the community working on absorbing-state transitions and open quantum systems, but the current form is better suited to publication if the CAM claims are softened or, preferably, if the exponent extraction is strengthened with uncertainty quantification and a test of the assumed scaling forms. Independent numerical evidence, for example from tensor-network or quantum-jump Monte Carlo simulations of the actual 1D HCP, would substantially increase confidence in the central claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does a genuinely useful thing: it defines a hybrid contact process that interpolates between the classical contact process and the quantum contact process, and it maps out the mean-field phase diagram with absorbing, active, and bistable regions. The single-site and cluster mean-field calculations are internally consistent, and the Liouvillian spectrum analysis adds a reasonable check on relaxation dynamics. The model itself is new, and the observation that the absorbing-to-active transition stays continuous while the absorbing-to-bistable transition is discontinuous is worth having on record.\n\nThe soft spot is the central quantitative claim. The extracted exponent β*≈0.79, which would put the transition outside the directed percolation class, rests on a coherent anomaly method fit to five or six cluster sizes (L=5–10) with no error bars. The scaling form in Eq. (23) is assumed, not derived for this model, and the exponential finite-size scaling for κ_c(L) is likewise ad hoc. The paper itself concedes that larger clusters are needed for a precise β*, yet the abstract and conclusion still present β*≈0.79 as the main result. I also note a mild circularity: the CAM critical point and the direct extrapolation of κ_c(L) use the same underlying CMF data, so their agreement is not an independent confirmation. The stress-test note is fair on all these points.\n\nThat said, the paper is not careless. The authors flag their own limitations, the mean-field analysis is thorough, and the phase diagram is a solid contribution. The problem is that the headline claim about a new universality class is over-reached relative to the numerical evidence. This is a fixable issue: larger clusters, uncertainty quantification, and an independent method (tensor network or Monte Carlo on the actual quantum master equation) would either firm up the claim or kill it.\n\nFor a reader, the paper is useful as a model introduction and a demonstration of the CMF+CAM machinery, but I would not yet cite it for the exponent. It deserves a serious referee because the model is new and the analysis is honest, but the referee should demand the missing evidence before accepting the universality-class conclusion. I would send it to peer review with a request for major revision rather than desk-reject it.","headline":"A well-executed but numerically fragile study of a new hybrid contact process; the claimed new universality class is plausible but not established with the current cluster sizes.","tokens_in":14572,"tokens_out":1157,"would_cite":false,"duration_ms":11957,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The one-dimensional hybrid contact process shows a continuous absorbing-to-active transition whose critical exponent, estimated at about 0.79, lies outside the directed percolation universality class.","keywords":["hybrid contact process","absorbing-state phase transition","coherent anomaly method","cluster mean-field approximation","directed percolation universality","Liouvillian spectrum","bistability","open quantum many-body systems"],"falsifier":"Run a direct steady-state simulation of the one-dimensional hybrid contact process at $\\Omega/\\Gamma = 0.1$ and fit the order parameter against $\\kappa$ near $\\kappa^* \\approx 0.9884\\,\\Gamma$; if the measured exponent is the directed percolation value $\\beta \\approx 0.276$ rather than about $0.79$, the coherent anomaly extrapolation is wrong. Equivalently, compute cluster mean-field results for $L = 12, 14, 16$ and check whether the CAM estimates of $\\beta^*$ remain stable near $0.79$.","tokens_in":53,"feed_emoji":"⚛️","tokens_out":12506,"duration_ms":314708,"temperature":0.7,"pith_summary":"This paper studies the one-dimensional hybrid contact process, a spin-chain model in which coherent (quantum) and incoherent (classical) branching and coagulation compete with local decay. Using single-site and cluster mean-field approximations together with the coherent anomaly method, it maps the steady-state phase diagram and identifies a discontinuous absorbing-to-bistable transition and a continuous absorbing-to-active transition. The paper's central quantitative claim is that the continuous transition has a nonclassical true critical exponent: $\\beta^* \\approx 0.7962$ for $\\Omega/\\Gamma = 0.1$ and $\\beta^* \\approx 0.7865$ for $\\Omega/\\Gamma = 0.4$, values that are close to each other and far from the directed percolation value. If the extraction is correct, adding quantum coherence changes the universality class of the one-dimensional absorbing-state transition.","feed_headline":"Transition in hybrid contact process leaves directed-percolation class","feed_subtitle":"If confirmed, quantum coherence puts this absorbing transition in a universality class of its own.","key_machinery":"The central object is the Lindblad master equation for a spin-1/2 chain, with a coherent Hamiltonian $\\hat H = \\Omega \\sum_j (\\hat\\sigma^x_j \\hat n_{j+1} + \\hat n_j \\hat\\sigma^x_{j+1})$, correlated incoherent jump operators at rate $\\kappa$, and local decay at rate $\\Gamma$. The load-bearing analysis is the coherent anomaly method (CAM): cluster mean-field (CMF) solutions of increasing size $L$ give pseudo-critical points $\\kappa_c(L)$ and order-parameter amplitudes $n_0(L)$, and CAM assumes $n_0(L) \\sim C_0 [(\\kappa^* - \\kappa_c(L))/\\kappa^*]^{-(\\beta^* - \\beta_{\\rm mf})}$ with $\\beta_{\\rm mf}=1$, then fits this relation to extract the true $\\kappa^*$ and $\\beta^*$. The CMF self-consistent effective-field condition captures both stable and unstable steady states, and the Liouvillian spectrum is used to characterize the relaxation dynamics.","core_discovery":"On its own terms, the paper establishes that the steady state of the one-dimensional hybrid contact process has three regimes: an absorbing phase, an active phase, and a bistable region. The boundary between the absorbing phase and the bistable region is a saddle-node bifurcation, so that transition is discontinuous, while the absorbing-to-active boundary is continuous. Combining cluster mean-field results for cluster sizes $L=5$ through $10$ with the coherent anomaly ansatz, the paper estimates the true critical points $\\kappa^*/\\Gamma = 0.9884$ and $0.9862$ for $\\Omega/\\Gamma = 0.1$ and $0.4$, respectively, and the true order-parameter exponents $\\beta^* = 0.7962$ and $0.7865$. Since those exponents are very close to each other and clearly differ from both the classical mean-field value $1$ and the one-dimensional directed percolation value, the paper concludes that the continuous transition likely defines its own universal critical behavior, while noting that larger clusters are needed for a more precise exponent.","pith_inferences":["A decisive test would be a direct large-system simulation of the steady-state density (quantum-jump Monte Carlo or tensor network) at $\\Omega/\\Gamma = 0.1$; if the measured exponent returns to the directed percolation value rather than staying near $0.79$, the coherent anomaly extrapolation is the source of the apparent new class.","If the exponent stays near $0.79$ as $\\Omega/\\Gamma \\to 0$, the universality change would persist in the weak-coherence limit; if it instead approaches the DP value, coherence must exceed a threshold to alter the class.","The absence of metastability and the nonmonotonic Liouvillian gap suggest the continuous transition is genuinely second-order, making the HCP a clean testbed for coherent anomaly extrapolations in open quantum systems."],"forward_implications":["For small coherent coupling, sweeping the incoherent rate $\\kappa$ at fixed $\\Omega$ crosses a continuous absorbing-to-active transition, while sweeping $\\Omega$ at fixed $\\kappa$ crosses a discontinuous absorbing-to-bistable boundary.","The two estimated exponents, $\\beta^* = 0.7962$ and $0.7865$, are close enough that the continuous transition appears to belong to a single universality class across coherent coupling strengths.","That universality class is not directed percolation: the extracted $\\beta^*$ differs from the one-dimensional DP value and from the classical mean-field value $1$.","The CAM critical points agree with the direct exponential finite-size extrapolation of the cluster mean-field critical points, giving consistent locations for the transition."],"supporting_citations":[{"why":"Compiles the universality classes for absorbing-state transitions, providing the directed percolation exponent that the reported $\\beta^*$ is compared against.","marker":"[27]"},{"why":"The quantum contact process study this model interpolates from, whose phase diagram and method this paper extends.","marker":"[39]"},{"why":"Introduces the cluster mean-field approximation and effective-field parameterization used to obtain the size-dependent steady states.","marker":"[41]"},{"why":"Derives the coherent anomaly scaling relation for mean-field singularities that underlies Eq. (23).","marker":"[46]"},{"why":"Adapts the coherent anomaly method to dissipative phase transitions, giving the specific extrapolation used to extract $\\kappa^*$ and $\\beta^*$.","marker":"[48]"}],"fun_headline_variants":["Hybrid contact process hints at new universality class","1D absorbing transition leaves mean-field and DP universality","Continuous absorbing transition in 1D defines novel universality class","Hybrid contact process: two transitions, one in new universality class","Nonclassical exponent found in 1D absorbing transition"],"cache_read_input_tokens":16640,"weakest_assumption_plain":"The extraction of $\\beta^*$ assumes the coherent anomaly scaling relation (Eq. 23) and the exponential finite-size form for $\\kappa_c(L)$ correctly describe how the cluster mean-field results approach the true transition, and the fits use only five or six cluster sizes from $L=5$ to $10$, so the reported exponents inherit that assumption's validity.","fun_headline_variants_meta":{"raw":{"variants":["Hybrid contact process hints at new universality class","1D absorbing transition leaves mean-field and DP universality","Continuous absorbing transition in 1D defines novel universality class","Hybrid contact process: two transitions, one in new universality class","Nonclassical exponent found in 1D absorbing transition"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000692,"raw_usage":{"total_tokens":3087,"prompt_tokens":857,"completion_tokens":2230,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":473,"completion_tokens_details":{"reasoning_tokens":2146}},"tokens_in":473,"tokens_out":2230,"duration_ms":16960,"temperature":1.0,"reasoning_tokens":2146,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:39:03.166847+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a direct steady-state simulation of the one-dimensional hybrid contact process at $\\Omega/\\Gamma = 0.1$ and fit the order parameter against $\\kappa$ near $\\kappa^* \\approx 0.9884\\,\\Gamma$; if the measured exponent is the directed percolation value $\\beta \\approx 0.276$ rather than about $0.79$, the coherent anomaly extrapolation is wrong. Equivalently, compute cluster mean-field results for $L = 12, 14, 16$ and check whether the CAM estimates of $\\beta^*$ remain stable near $0.79$.","supporting_citations":[{"cited_title":"Shang, S","cited_arxiv_id":null,"evidence_quote":"The quantum contact process study this model interpolates from, whose phase diagram and method this paper extends."},{"cited_title":"Suzuki, M","cited_arxiv_id":null,"evidence_quote":"Derives the coherent anomaly scaling relation for mean-field singularities that underlies Eq. (23)."},{"cited_title":"Jin, W.-B","cited_arxiv_id":null,"evidence_quote":"Adapts the coherent anomaly method to dissipative phase transitions, giving the specific extrapolation used to extract $\\kappa^*$ and $\\beta^*$."}],"review_version":1}