{"id":"0e8c9c2f-7fae-4d71-8021-c2584a610c7d","arxiv_id":"1908.01494","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In a fully connected open Ising model driven by classical colored noise, pink noise speeds up relaxation, blue noise slows it down, and the strong-coupling regime shows metastable nonzero magnetization.","lead":"An open Ising model made from superconducting qubits is simulated under white, pink, and blue classical noise. Pink noise accelerates the system's relaxation while blue noise slows it down, and strong qubit couplings produce a metastable plateau in the magnetization.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Pink-noise master equation is used outside its stated weak-coupling regime: normalization at πf0≈1570Γ makes the effective pink dissipative rate ~150Γ, exceeding ϵ=10Γ.","rationale":"I identified a different concern from the reader's finite-size worry. The finite-size N=8 mainly threatens the metastable part of the abstract; the colored-noise acceleration/slowing is a single-qubit-level spectral effect and would likely survive at larger N. The more load-bearing issue is whether Eq. (6) is valid at all for the pink-noise parameters: the noise power at the system's transition frequencies is ~150 times larger than the nominal Γ, so the perturbative small parameter is ~15 instead of 0.1. If the trajectory test confirms Eq. (6), the concern does not land and the paper's qualitative claims are supported; if it does not, the central quantitative comparison is an artifact. Since this is a testable revision rather than a demonstrated contradiction, the appropriate outcome remains conditional acceptance pending the check, matching the reader's verdict category; hence UNCHANGED.","tokens_in":15130,"tokens_out":28620,"duration_ms":283840,"concrete_test":"Compute m(t) by exact ensemble averaging of the stochastic Schrödinger equation for Hamiltonian (2): generate an ensemble of pink/blue noise realizations with the same PSDs and parameters (N=8, ϵ/Γ=10, λ/ϵ=1 and 10), integrate the time-dependent Schrödinger equation numerically (e.g., Euler–Maruyama with a fine time step), average the resulting trajectory magnetizations, and compare with Fig. 3 (middle and right). If the trajectory-averaged curves disagree with Eq. (6) or if pink no longer accelerates when the PSD is renormalized at a transition frequency so that the effective coupling is genuinely weak (S_pink(2ϵ) ≪ ϵ), the central claim is not established. This check directly tests whether the second-order TCL master equation is valid in the simulated parameter regime.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is obtained by solving the second-order time-convolutionless master equation (6), whose derivation is explicitly restricted to ϵ ≫ Γ (Section II, Fig. 2). The numerical simulations set ϵ/Γ = 10, 2f0/Γ = 10^3, and normalize all noise PSDs to the same value at the Nyquist frequency πf0 ≈ 1570Γ. But the system's relevant Bohr frequencies are ω ~ 2ϵ = 20Γ (and λ-scale splittings up to ~100Γ). For pink noise, Sη(2πf) ∝ 1/f, so the noise power at ω = 20Γ is larger than the white-noise reference by a factor ≈ (πf0)/(ω/2) ≈ 157. Consequently the actual dissipative rate entering Eq. (6) in the pink case is of order 100Γ, i.e., an order of magnitude larger than the unperturbed energy scale ϵ = 10Γ. The weak-coupling condition under which the master equation was derived is therefore violated precisely in the pink-noise case that is claimed to accelerate relaxation; the blue case is suppressed by the same normalization. The reported pink/blue asymmetry may be an artifact of extrapolating a perturbative TCL equation far outside its radius of validity, not a genuine memory effect. The paper does not provide a non-perturbative check (e.g., exact trajectory averages of the original stochastic Hamiltonian (2)).","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies an N-qubit fully connected Ising model realized in a superconducting charge-qubit circuit, with each qubit subject to classical white, pink, or blue noise. Starting from a circuit Hamiltonian derived in Appendix A, the authors use a time-convolutionless master equation from Ref. [14] to compute the relaxation of the magnetization m(t). For white noise they find a metastable plateau m_ms != 0 in the strong-interaction regime lambda/epsilon >> 1, supported by a Liouvillian spectral analysis in Appendix C. For colored noise, they report that pink noise accelerates relaxation relative to the white-noise Markovian limit, while blue noise slows it down, and they attribute this to positive versus negative memory in the noise correlation function. The paper also examines ground-state weight and von Neumann entropy for the three noise colors.","tokens_in":15447,"tokens_out":14342,"duration_ms":158277,"significance":"If the colored-noise results were robust, the paper would provide a useful, experimentally motivated step toward engineering the environment of superconducting many-body simulators. The manuscript has clear strengths: the circuit derivation in Appendix A is explicit, the white-noise metastability is diagnosed through the Liouvillian spectrum rather than asserted, and the numerical methods are standard and described sufficiently for reproduction. However, the central pink/blue comparison is not currently established: the pink-noise operating point lies outside the stated validity regime of the master equation, and the normalization convention used for the three spectra may predetermine the reported ordering of relaxation rates. The white-noise metastability result appears sound and could stand independently; the colored-noise claims need additional work, including a non-perturbative check or a properly controlled normalization.","major_comments":[{"comment":"The pink-noise simulations violate the assumption epsilon >> Gamma used to derive Eq. (6). The paper sets epsilon/Gamma = 10 and 2f0/Gamma = 10^3, so the Nyquist frequency is pi f0 ≈ 1570 Gamma. Because all noise spectra are normalized to the same value at pi f0 and S_pink ∝ 1/f, the pink spectral density at a Bohr frequency Delta omega = 2 epsilon = 20 Gamma is S_pink ≈ (f0 / (Delta omega / 2 pi)) S_white ≈ (500 / 3.18) S_white ≈ 157 S_white. The dissipative rate entering Eq. (6) is therefore of order 100 Gamma for pink noise, an order of magnitude larger than epsilon = 10 Gamma. This is precisely the regime in which the paper states that Eq. (6) is no longer valid, because the noise term becomes comparable to Hs. The reported acceleration by pink noise may thus be an artifact of using a second-order time-convolutionless equation outside its radius of validity, while the blue case is suppressed by the same normalization. I ask for a non-perturbative check, e.g., exact averaging over trajectories of the stochastic Hamiltonian (2), or a reparameterization for which the effective dissipative rates are at most of order Gamma for all three colors.","section":"Section II and IV, Eq. (6) and Fig. 3"},{"comment":"The representativeness of the N = 8 simulations is asserted, not demonstrated. The paper states that N = 8 is used 'to capture the general features of larger systems', but the Liouvillian spectral analysis in Fig. 6 is actually performed at N = 5 because of matrix-size limits, and no finite-size scaling or large-N argument is provided. For the all-to-all Hamiltonian (3), the spectrum and the Liouvillian gap structure change with N, so the metastable plateau and the relaxation rates extracted from N = 8 may not transfer to the large-scale processors mentioned in the abstract. Please add N-dependence scans (e.g., N = 4, 6, 8, 10) or a mean-field/large-N analysis to support the 'general features' claim.","section":"Section III, after Eq. (14), and Appendix C"},{"comment":"Equation (19) as written is not equivalent to Eq. (6). For Hs eigenstates, the relevant matrix element of U_s(t,t') sigma_x U_s^dagger(t,t') is <alpha|sigma_x|alpha'> exp[-i Delta omega_{alpha alpha'}(t-t')], so after the time integral the coefficient should contain exp[-i Delta omega_{alpha alpha'} t] times an integral over exp[+i Delta omega_{alpha alpha'} t'] K(t,t'), or an equivalent rotating-frame expression. The present definition K^{alpha alpha'}_k(t) = integral_0^t K(t,t') exp[-i Delta omega t'] dt' has the wrong phase and no compensating factor. Since the colored-noise results are obtained by directly integrating Eq. (19), the phase convention affects which frequency components of the noise are sampled. Please correct the equation and repeat the numerical calculation, or explicitly specify the rotating frame in which Eq. (19) is intended to hold.","section":"Section IV, Eq. (19)"},{"comment":"The central interpretation of the colored-noise results does not isolate memory effects from the imposed spectral shape. The authors normalize all PSDs to the same value at pi f0 and then observe that for all f < f0 the pink PSD lies above and the blue PSD lies below the white value. With the stated slopes S ∝ f^-1 and S ∝ f, this ordering is fixed by construction, so the resulting hierarchy of relaxation rates is largely predetermined by the normalization. To make the claim that positive memory enhances and negative memory suppresses relaxation non-tautological, the authors should fix a common physical noise strength, for instance equal spectral density at the system's characteristic transition frequency or equal integrated variance over the relevant frequency band, and show that the asymmetry persists under that comparison.","section":"Section IV, paragraph after Eq. (20)"}],"minor_comments":[{"comment":"The Fourier transform in Eq. (18) should be S_sigma(omega) = integral C_a(tau) exp(-i omega tau) d tau; as written it integrates over t while C_a is a function of tau.","section":"Section III, Eq. (18)"},{"comment":"The text says the eigenstates are labeled alpha = 0, 1, ..., 2N; since the Hilbert space of N qubits has dimension 2^N, this should read alpha = 0, 1, ..., 2^N - 1.","section":"Section IV, basis enumeration"},{"comment":"The manuscript alternates between 'pink' and 'red' noise: the abstract and most of Section IV use 'pink', but the text before Eq. (8) and the Fig. 3 caption use 'red-colored'; please harmonize the terminology.","section":"Section IV and Fig. 3 caption"},{"comment":"There is a typo in the Introduction: 'have alreadly provided' should be 'have already provided'.","section":"Introduction"},{"comment":"References [32] and [47] are the same Nakamura, Pashkin, and Tsai paper; one duplicate should be removed or cross-referenced.","section":"References"},{"comment":"The statement that all transition frequencies |Delta omega| are smaller than pi f0 'due to pi f0 > epsilon and lambda' is not correct as written, because the spectral width of Hs scales as O(N(epsilon + lambda)) for N qubits; the condition should be stated for the actual spectral support of the transitions coupled by sigma_x.","section":"Section IV, spectral support statement"}],"recommendation":"major_revision","confidential_remarks":"The colored-noise central claim depends on a noise normalization that may predetermine the result, and the pink-noise operating point appears to violate the stated weak-coupling assumption; the white-noise metastability part is more solid. I have no additional remarks beyond the referee report."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a clean, well-written numerical study of a dissipative fully-connected Ising model under classical colored noise, using the time-convolutionless master equation of Chenu et al. The white-noise section is genuinely useful: the metastable plateau at strong coupling is a nice demonstration, and the spectral-gap argument in the appendix is plausible. The pink/blue comparison, however, is built on an arbitrary normalization and—in the pink case, at least—the equations are used outside the weak-coupling regime they were derived for. That makes the central 'pink accelerates, blue slows' claim less solid than the abstract suggests.\n\nWhat's new: applying this master equation to a fully-connected SC-qubit Ising model is a legitimate new application. The metastability in the strong-coupling regime, supported by the relaxation-mode analysis, is a nice addition. The white-noise MCWF results, including jump statistics and correlation functions, are a solid piece of numerics.\n\nWhere it gets shaky. First, the normalization. All PSDs are set equal at the Nyquist frequency πf0≈1570Γ. Since the system has transition energies around 2ϵ=20Γ, the pink PSD at those frequencies is about 80–150 times larger than the white reference. The effective dissipative rate in Eq. (6) for pink is then of order 100Γ, not Γ. That is not a weak perturbation on ϵ=10Γ. The TCL equation (6) is explicitly derived under ϵ≫Γ; solving it for pink is like using lowest-order perturbation theory past the point where the perturbation is strong. To your earlier question: I think the stress-test is right, and it is not answered by the paper. A quick fix would be to normalize all noises to the same PSD at the actual qubit frequency, or to show that a non-perturbative trajectory simulation of Eq. (2) gives the same asymmetry. The authors do note the PSD interpretation in Sec. IV, so they are aware of the effect of normalization; but they don't discuss the validity boundary.\n\nSecond, the metastable result is only for N=8 with no finite-size scaling or large-N analysis, and no code/data released. The gap structure could change with N.\n\nMinor: the rough approximation in Sec. IV (setting the phase factors to 1) is very crude; it gives a nice intuition but also makes the 'positive memory' language sound more causal than it is.\n\nWho's this for? People working on noisy superconducting circuits and open many-body systems will find the white-noise part worth a look. The colored-noise comparison needs more care. I'd send it to review, but the referee should ask for a check of the pink-noise regime and a finite-size study, or a reframing of the claims. I'd accept for peer review with a request for major revision, not a desk reject. The metastability result is worth publishing on its own if the colored-noise part is fixed. That's my take.","headline":"Pink-noise comparison is built on a normalization that pushes the effective dissipation far outside the weak-coupling limit; the white-noise metastability result is the more solid part.","tokens_in":15945,"tokens_out":10354,"would_cite":false,"duration_ms":101995,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["85.25.-j","42.50.-p","06.20.-f"],"model":"deepseek-v4-flash","headline":"Noise color controls relaxation of an open Ising model: pink noise accelerates decay, blue noise slows it, and strong qubit coupling yields a metastable nonzero magnetization.","keywords":["open quantum systems","non-Markovian dynamics","Ising model","superconducting qubits","colored noise","metastability","quantum simulation","magnetization relaxation"],"falsifier":"Measure the magnetization relaxation in a superconducting qubit array with gate voltage filtered to pink and blue noise and compare it with the white-noise case; alternatively, compute the spectral gap and metastable magnetization for N = 10, 12, 16 or via a mean-field treatment and check whether the plateau and the pink/blue relaxation ordering survive.","tokens_in":14917,"feed_emoji":"⚛️","tokens_out":4362,"duration_ms":43262,"temperature":0.7,"pith_summary":"The paper claims that the color of environmental noise changes how fast a fully connected Ising spin system relaxes. It models eight superconducting charge qubits weakly coupled to classical white, pink, or blue noise, using a non-Markovian master equation. In the strong-interaction regime, the magnetization first relaxes to a metastable nonzero value before eventually decaying to zero. Pink noise, whose correlations are positive, accelerates relaxation; blue noise, with negative correlations, slows it down. If true, this gives a practical handle for designing qubit environments in large quantum processors.","feed_headline":"Pink noise relaxes a qubit Ising array faster, blue slower","feed_subtitle":"In a superconducting-circuit model, noise memory sets the relaxation speed, and strong coupling creates a metastable plateau.","key_machinery":"The central object is the non-Markovian master equation with a time-convolution integral over the noise correlation function $K_{k,k'}(t,t')$, in the local-in-time form derived from a perturbation formalism. The sign of the noise memory enters through the effective time-dependent decay rate $\\tilde{\\Gamma}_k(t)=\\Gamma\\int_0^t K_{k,k}(t,t')\\,dt'$, which is larger than $\\Gamma$ for pink noise and smaller for blue noise. Metastability is explained by diagonalizing the Markovian Liouvillian: in the strong-coupling regime, several relaxation modes have decay rates $\\gamma_\\mu$ close to $\\gamma_0=0$, producing a long-lived transient plateau in the magnetization.","core_discovery":"The dissipative dynamics of a fully connected quantum Ising model, simulated with superconducting charge qubits and perturbed by classical colored noise, depends on the sign of the noise memory. For pink noise, the effective decay rate is enhanced, so the system reaches its stationary state faster than with white noise; for blue noise, the decay rate is suppressed. In the strong inter-qubit coupling limit ($\\lambda/\\epsilon \\gg 1$), the Markovian white-noise system exhibits metastability: the magnetization $m(t)$ settles to a nonzero metastable value $m_{\\rm ms}\\neq 0$ before decaying to the true stationary state $m_{\\rm ss}=0$.","pith_inferences":["The sign-of-memory rule suggests a predictable monotonic ordering for any $1/f^\\alpha$ noise: as the spectral exponent $\\alpha$ is swept from negative (blue-like) to positive (pink-like), relaxation should interpolate between the blue and pink extremes, which could be tested on the same circuit by changing the filter.","If the metastable plateau survives finite-size scaling, the plateau value and lifetime as a function of $\\lambda/\\epsilon$ could serve as a sensitive probe of the environment's noise spectrum.","Because the derivation assumes weak coupling ($\\epsilon \\gg \\Gamma$), the claimed acceleration or slowing might change in the ultrastrong-coupling regime; extending the treatment beyond perturbation would test whether the memory-sign rule is universal.","For the fully connected model, a mean-field or large-$N$ analysis could confirm whether the metastable gap structure seen at $N=8$ persists or is a finite-size artifact."],"forward_implications":["In the strong-coupling regime ($\\lambda/\\epsilon \\gg 1$), the magnetization first relaxes to a nonzero metastable value $m_{\\rm ms}<0$ before decaying to the stationary state.","Pink-colored noise accelerates relaxation toward the true stationary state relative to white noise, shortening both the metastable period and the ground-state lifetime.","Blue-colored noise slows relaxation, so within the same observation window the system may not even reach its metastable plateau.","The single-spin power spectrum shows a narrow central peak plus side peaks, with side-peak positions set by $\\pm(\\lambda/4N - 2\\epsilon)$ and integer multiples.","The acceleration or suppression of relaxation is governed by the sign of the noise correlation: positive memory enhances the effective decay rate, negative memory reduces it."],"supporting_citations":[{"why":"Supplies the perturbation formalism that yields the non-Markovian master equation used throughout the paper.","marker":"[14]"},{"why":"Provides the digital-filtering algorithm that generates pink and blue colored noise fields with the specified power spectra.","marker":"[46]"},{"why":"Gives the Monte Carlo wave-function method used to simulate the Markovian white-noise limit and record quantum-jump trajectories.","marker":"[51]"},{"why":"Establishes the spectral criterion for metastability that the paper uses to interpret the slow relaxation modes.","marker":"[55]"},{"why":"Defines the charge-qubit (Cooper-pair box) building block that forms the simulated Ising system.","marker":"[47]"}],"fun_headline_variants":["Pink noise relaxes qubit arrays faster, blue slower","Noise memory controls relaxation in superconducting Ising model","Colored noise alters decay rate in open quantum Ising chain","Strong coupling yields metastable plateau in qubit relaxation","Pink vs blue noise: which relaxes a qubit chain faster?"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central results are computed for an eight-qubit system, and the paper assumes, without finite-size scaling, that these eight qubits capture the behavior of much larger qubit arrays.","fun_headline_variants_meta":{"raw":{"variants":["Pink noise relaxes qubit arrays faster, blue slower","Noise memory controls relaxation in superconducting Ising model","Colored noise alters decay rate in open quantum Ising chain","Strong coupling yields metastable plateau in qubit relaxation","Pink vs blue noise: which relaxes a qubit chain faster?"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000908,"raw_usage":{"total_tokens":3827,"prompt_tokens":789,"completion_tokens":3038,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":405,"completion_tokens_details":{"reasoning_tokens":2954}},"tokens_in":405,"tokens_out":3038,"duration_ms":22610,"temperature":1.0,"reasoning_tokens":2954,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:11:22.464890+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the magnetization relaxation in a superconducting qubit array with gate voltage filtered to pink and blue noise and compare it with the white-noise case; alternatively, compute the spectral gap and metastable magnetization for N = 10, 12, 16 or via a mean-field treatment and check whether the plateau and the pink/blue relaxation ordering survive.","supporting_citations":[{"cited_title":"Zhivomirov, Romanian Journal of Acoustics and Vi- bration 15, 14 (2018)","cited_arxiv_id":null,"evidence_quote":"Provides the digital-filtering algorithm that generates pink and blue colored noise fields with the specified power spectra."},{"cited_title":"Macieszczak, M","cited_arxiv_id":null,"evidence_quote":"Establishes the spectral criterion for metastability that the paper uses to interpret the slow relaxation modes."},{"cited_title":"Nakamura, Yu","cited_arxiv_id":null,"evidence_quote":"Defines the charge-qubit (Cooper-pair box) building block that forms the simulated Ising system."}],"review_version":1}