{"id":"e7b029c3-5580-42b6-854a-5090cb20fc24","arxiv_id":"2501.10349","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Single-photon Langevin simulations show chiral state transfer near a Liouvillian exceptional point is transient, with chirality collapsing onto a fitted scaling curve C=f(γ2 T^(1/ν)) as a function of encircling time and dephasing.","lead":"This paper reports a photonic experiment that simulates an open two-level quantum system and shows that chiral state transfer near a Liouvillian exceptional point appears only for intermediate encircling times and vanishes as the system relaxes to its steady state. It matters because it tests exceptional-point dynamics in honest open systems and offers a single-photon way to simulate open-system master equations.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (12) renormalizes each trajectory by the largest singular value of U_j rather than the actual survival probability for the initial state; if this is wrong, the measured chirality and the universal scaling collapse are reconstruction artifacts.","rationale":"The reader's weakest assumption identifies the density-matrix reconstruction as the key risk; my read agrees and sharpens it into a concrete mathematical inconsistency. Eq. (12) weights each tomographically normalized trajectory by the largest singular value of U_j, whereas the exact unnormalized contribution is U_jρ(0)U_j†, whose trace is the initial-state-specific survival probability. These two weights differ for the nonnormal, non-scalar U_j produced by the stroboscopic evolution (5), especially for the initial state |1⟩⟨1| used in Fig. 4. Because the chirality C is a trace-distance between ensemble-averaged density matrices, any systematic weighting error changes both the peak position and the long-time decay, which are the two phenomena that support the transient-chirality and universal-scaling claims. The n=10 ensemble is a secondary but related issue: even with correct weights, ten realizations of a stochastic unraveling may be insufficient to resolve a trace-distance observable to the precision claimed. I do not reject the paper outright, because the qualitative claim that dissipative gap shrinks chirality at long times is physically plausible and the experimental platform is sophisticated; the correct response is to require the numerical test described above. If the test shows that Eq. (12) reproduces the exact unravelling and master-equation results, the conditional acceptance remains appropriate. If it does not, the universal scaling result and the extracted exponent ν=1.7221 should be withdrawn or re-derived with the correct normalization. Thus the reader's CONDITIONAL verdict is unchanged pending this check.","tokens_in":10790,"tokens_out":6254,"duration_ms":72878,"concrete_test":"Recompute the final-time density matrix for the Fig. 4 parameters using the exact unnormalized ensemble average ρ_exact(T) = (1/n) Σ_j U_j |1⟩⟨1| U_j†, with the same stochastic noise schedule and U_j as in Eq. (5), and compare the resulting chirality C(T,γ2) with Eq. (12). Also compare both against the direct Lindblad master-equation solution of Eq. (1) at the same parameters. If the Eq. (12) chirality surface or the fitted scaling exponent ν differs from the exact-unravelling or master-equation result beyond the quoted error bars, the reconstruction formula is invalid and the universal scaling claim is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing assumption is the density-matrix reconstruction formula in Eq. (12). The paper claims the ensemble average of Eq. (12) reproduces the Lindblad dynamics of Eq. (1), but this is not the standard unravelling average. In the quantum Langevin equation (4), the j-th unnormalized state is U_j|ψ0⟩, whose contribution to the ensemble density matrix has trace ⟨ψ0|U_j†U_j|ψ0⟩. Equation (12) instead multiplies the tomographically normalized ρ_j(T) by max|ξ_j|, the largest singular value squared of U_j. These weights coincide only if |ψ0⟩ is the principal right singular vector of every U_j, which is not true for the general two-level operators generated by Eqs. (4)-(10) and the initial state |1⟩⟨1| used in Fig. 4. No argument is given that the largest singular value equals the survival probability of the actual initial state. Since the chirality C in Eq. (11) and the universal scaling C = f(γ2 T^(1/ν)) in Fig. 4 are computed from this reconstructed ρ(T), an incorrect normalization would make the central quantitative claims artifacts of the reconstruction procedure. The small ensemble size n=10 compounds the problem: even with correct trace weights, no convergence analysis is shown for the trace-distance observable C.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a single-photon interferometric experiment that simulates the Lindblad dynamics of a two-level open system via a quantum Langevin equation. By parametrically encircling a Liouvillian exceptional point with different total encircling times, the authors observe chiral state transfer at an intermediate time (T=90) and its disappearance at long times (T=600), where the system relaxes to the steady state. They quantify the chirality with the trace-distance definition in Eq. (11) and claim a universal scaling C = f(γ2 T^{1/ν}) with an exponent ν=1.7221 obtained by fitting the experimental data.","tokens_in":11101,"tokens_out":5001,"duration_ms":53899,"significance":"If the central claims hold, the paper provides a valuable experimental demonstration that chiral state transfer near a Liouvillian exceptional point is a transient effect, in contrast to the long-time chiral transfer familiar from non-Hermitian Hamiltonian dynamics. The experimental scheme, based on reconstructing density-matrix evolution from stochastic wave-function realizations, is original and the data in Figs. 2 and 3 are shown to be in reasonable agreement with master-equation curves. However, the reconstruction formula in Eq. (12), the absence of a convergence analysis for the n=10 ensemble, and the fitting of the scaling exponent to the same data used to demonstrate collapse leave the quantitative claims—especially the universal scaling—insufficiently supported.","major_comments":[{"comment":"The reconstruction weight in Eq. (12) is not the standard unravelling weight and is not derived. For the linear quantum Langevin equation (4), the unnormalized contribution of trajectory j is U_j|ψ0⟩, whose correct weight in the ensemble density matrix is ⟨ψ0|U_j† U_j|ψ0⟩. Equation (12) instead multiplies the tomographically normalized ρ_j(T) by max|ξ_j|, the largest eigenvalue of U_j U_j† (i.e., the squared largest singular value). These weights coincide only if |ψ0⟩ is the principal right singular vector of every U_j. The initial state in Fig. 4 is |1⟩⟨1|, and no argument is given that this state is the principal singular vector of the operators generated by Eqs. (4)–(10). Because the chirality C in Eq. (11) and the scaling collapse in Fig. 4 are computed from this reconstructed ρ(T), the central quantitative claims rest on an unjustified normalization. Please derive Eq. (12) from Eq. (4), or explicitly benchmark both reconstruction rules against a direct numerical solution of the Lindblad master equation (1) for the same parameters.","section":"Eq. (12), together with Eqs. (4) and (11)"},{"comment":"The paper states that n=10 'provides a good enough estimation of the density-matrix dynamics,' but no convergence analysis is shown. For nonunitary trajectories with state-dependent norm fluctuations, ten realizations can bias the ensemble average of ρ(T) and therefore bias the measured chirality C. Please present the ensemble average or the trace-distance observable C as a function of n for representative parameters (for example, T=90, γ2=2), including statistical error bars, and justify why n=10 is sufficient for the scaling collapse in Fig. 4(b).","section":"Section 'Simulation of open systems' and Fig. 4"},{"comment":"The universal scaling C = f(γ2 T^{1/ν}) is demonstrated with an exponent ν=1.7221 obtained by fitting the same experimental data that are then collapsed. As written, this is an empirical fit with a free parameter rather than an independent prediction of a scaling law. To substantiate the universality claim, the exponent should be derived from the Liouvillian spectral structure or from an independent numerical dataset, and the fit uncertainty should be reported.","section":"Fig. 4(b) and the paragraph 'Chirality'"}],"minor_comments":[{"comment":"Please define ξ_j unambiguously: is it an eigenvalue of U_j U_j†, and is the maximum taken over the eigenvalues of a single trajectory or over trajectories? The current notation 'max|ξ_j|' is unclear.","section":"Eq. (12)"},{"comment":"The caption contains a grammatical error: 'States |1⟩ is coupled' should read 'State |1⟩ is coupled'.","section":"Fig. 1(a) caption"},{"comment":"The abstract states that chirality scales with 1/T^ν, while the main text and Fig. 4 use C = f(γ2 T^{1/ν}); please reconcile this notation.","section":"Abstract and 'Chirality' section"},{"comment":"The figure legend distinguishes hollow squares and solid dots, but the caption does not explain which symbol corresponds to clockwise and which to counterclockwise encircling, nor what the solid curves represent.","section":"Fig. 2"},{"comment":"The data availability statement says the data are available 'within the Letter [49]', but reference [49] is a supplemental materials file; please provide a persistent repository link or a stable DOI for the raw data.","section":"Data Availability Statement"}],"recommendation":"major_revision","confidential_remarks":"The manuscript builds heavily on the authors' own prior framework (Refs. [11,19,44]), and the central scaling claim is weakened by the fact that the exponent is fitted to the same data. The editor may wish to consider whether the incremental advance over these prior works is sufficiently articulated. The reconstruction issue in Eq. (12) is the most serious technical point; it is potentially fixable by a re-analysis or by explicit benchmarking against master-equation simulations."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The core of this paper is the experimental demonstration that chiral state transfer near a Liouvillian exceptional point is a transient effect: at long encircling times the system relaxes to the steady state and chirality disappears, while at intermediate times it survives. That message is plausible, the data in Figs. 2–3 are consistent with numerical master-equation curves, and the single-photon simulation of the quantum Langevin equation with full density-matrix reconstruction is a genuinely useful technical step. The qualitative conclusion is probably right and worth broadcasting.\n\nThe soft spots are all in the quantitative claims. The most load-bearing is Eq. (12). The paper renormalizes each tomographically normalized trajectory by max|ξ_j|, the largest singular value squared of U(τ), but the correct unravelling weight for the initial state |1⟩⟨1| is its survival probability ⟨ψ0|U_j†U_j|ψ0⟩. Those coincide only when |ψ0| is the principal right singular vector of every U_j, which is not generally true for the operators generated by their encircling path. Without a derivation or a numerical check, the reconstructed ρ(T), the chirality C, and the scaling collapse in Fig. 4 could be artifacts of the weighting. This is not a fatal flaw in the qualitative story, but it is a real gap and the paper should address it.\n\nSecond, n=10 Langevin realizations is small and no convergence analysis is shown for the trace-distance observable C. The error bars from Monte Carlo photon-counting statistics do not capture the ensemble-size bias. Third, the \"universal scaling\" C = f(γ2 T^(1/ν)) rests on ν = 1.7221 fitted to the same experimental data, so the collapse is a fit, not a prediction. The paper says the exponent is path-independent but parameter-dependent; that is honest, but it lowers the strength of the claim. Fourth, the theoretical framework is largely from the authors' own prior work (Refs. [11,19]), so the new content is the experimental implementation, not the physics.\n\nNone of this undermines the main qualitative result, and the citation pattern looks fair. The paper deserves a serious referee: the experiment is real, the technique is interesting, and the reconstruction issue is fixable with a clear derivation and convergence checks. I would engage with it, but I would not cite the scaling exponent as a law until the weighting is sorted out.","headline":"A careful photonic experiment showing that chiral state transfer near a Liouvillian EP is transient, with a reconstruction-weighting issue in Eq. (12) that needs scrutiny before the quantitative scaling claim can be trusted.","tokens_in":11591,"tokens_out":1349,"would_cite":true,"duration_ms":16052,"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":"This paper claims that chiral state transfer near a Liouvillian exceptional point is a transient effect: it appears only for intermediate encircling times, disappears at long times as the system relaxes to its steady state, and obeys a…","keywords":["chiral state transfer","Liouvillian exceptional point","open quantum system","quantum Langevin equation","single-photon interferometry","universal scaling","density-matrix reconstruction","non-Hermitian dynamics"],"falsifier":"Run the same encircling protocol with the density matrix computed directly from the Lindblad master equation, or with a much larger stochastic ensemble, at the parameters of Fig. 4; if the chirality values and the collapse $C=f(\\gamma_2 T^{1/\\nu})$ with $\\nu\\approx1.7221$ do not survive, the observed chirality is an artifact of the ten-realization reconstruction. As a second check, encircle a loop of the same shape that does not enclose the Liouvillian exceptional point: if direction-dependent state flipping persists, the chirality is not caused by exceptional-point encircling.","tokens_in":10581,"feed_emoji":"⚛️","tokens_out":7762,"duration_ms":74735,"temperature":0.7,"pith_summary":"The authors set out to show that in a genuinely open quantum system, the direction-dependent state flip associated with encircling a Liouvillian exceptional point is not a long-time feature of the dynamics. Using single-photon interferometry to simulate the density-matrix evolution through the quantum Langevin equation, they find that with dephasing present the chirality survives only for intermediate encircling times, such as their $T=90$ case, and vanishes for long encircling times such as $T=600$. They also report that the chirality obeys the scaling relation $C=f(\\gamma_2 T^{1/\\nu})$ with a fitted exponent $\\nu\\approx1.7221$. The result matters because it clarifies when Liouvillian exceptional points, which sit away from the steady state, can produce observable dynamical consequences in open systems.","feed_headline":"Chiral state transfer near Liouvillian exceptional point is transient","feed_subtitle":"A photonic experiment shows chirality survives only at intermediate times and follows a universal scaling law.","key_machinery":"The central objects are the Liouvillian superoperator $\\mathcal{L}$ of a two-level open system and its stochastic unraveling through the quantum Langevin equation $i\\frac{d}{dt}|\\psi(t)\\rangle=\\left(H-\\frac{i}{2}L_{\\phi}^{\\dagger}L_{\\phi}+i l(t)L_{\\phi}\\right)|\\psi(t)\\rangle$, where $l(t)$ is white noise and $L_{\\phi}$ is the dephasing jump operator. The experiment reconstructs the density matrix as an ensemble average of $n=10$ stochastic wave-function evolutions, implemented as nonunitary single-photon interferometry, and quantifies directionality with the chirality $C=\\frac{1}{2}\\operatorname{Tr}\\sqrt{(\\tilde{\\rho}_{\\mathrm{cw}}-\\tilde{\\rho}_{\\mathrm{ccw}})^{\\dagger}(\\tilde{\\rho}_{\\mathrm{cw}}-\\tilde{\\rho}_{\\mathrm{ccw}})}$. This machinery connects the observable state flip to the spectral landscape of the Liouvillian and produces the scaling collapse in $\\gamma_2 T^{1/\\nu}$.","core_discovery":"The authors claim that when a Liouvillian exceptional point is parametrically encircled in the presence of dephasing, chiral state transfer is a transient phenomenon. At sufficiently long encircling times the Liouvillian gap makes the system relax toward its instantaneous steady state in both encircling directions, so the final states merge and the chirality disappears; at an intermediate time, one direction follows the Liouvillian spectral landscape nearly adiabatically while the other undergoes a non-adiabatic jump, leaving a mixed final state and a finite chirality. They further claim that this intermediate-time chirality is purely due to the encircling of the Liouvillian exceptional point and that the final-time chirality follows the universal scaling $C=f(\\gamma_2 T^{1/\\nu})$, with $\\nu=1.7221$ determined by fitting the experimental data.","pith_inferences":["A microscopic derivation of the fitted exponent $\\nu\\approx1.7221$ from the Liouvillian spectrum would test whether the scaling is truly universal or specific to the chosen encircling path; the paper reports the scaling form as universal but the coefficient as parameter-dependent.","Because each reported density matrix uses only ten stochastic realizations, benchmarking the same parameters against an exact Lindblad master-equation solution would show how much of the chirality magnitude depends on the ensemble size.","The same simulation scheme could be applied to bistable or nonlinear open systems, where the paper notes a recent exception to the transient rule, to search for steady-state chiral switching."],"forward_implications":["At long encircling times, any dephased open system will relax toward its instantaneous steady state, so protocols relying on Liouvillian-exceptional-point chirality must operate in an intermediate time window.","The measured scaling $C=f(\\gamma_2 T^{1/\\nu})$ gives a quantitative rule for choosing encircling time and dephasing rate to maximize the chiral signal.","Comparing the two encircling directions at intermediate times reveals near-adiabatic following in one direction and a non-adiabatic jump in the other, providing a signature for identifying Liouvillian exceptional-point dynamics.","The single-photon stochastic-unraveling scheme reconstructs density-matrix evolution from wave-function trajectories and can therefore simulate other Lindblad open-system dynamics beyond this two-level example."],"supporting_citations":[{"why":"Establishes the equivalence between the Lindblad master equation and the quantum Langevin equation that the experiment uses to reconstruct the density matrix.","marker":"[42]"},{"why":"Supplies the stroboscopic single-photon realization of nonunitary evolution on which the interferometric gates are based.","marker":"[44]"},{"why":"Gives the theory of chiral state transfer under dephasing and defines how the state trajectory is projected onto the Liouvillian spectral landscape.","marker":"[11]"},{"why":"Demonstrates quantum jumps in the non-Hermitian dynamics of a superconducting qubit, the dissipative-qubit context this experiment extends.","marker":"[18]"},{"why":"Provides the chirality measure used in Eq. (11) and the decoherence-induced exceptional-point comparison.","marker":"[19]"},{"why":"Supplies the eigenvalue renormalization used to convert nonunitary evolution into an ensemble-averaged density matrix in Eq. (12).","marker":"[45–47]"},{"why":"Provides the general quasi-adiabatic description of exceptional-point encircling dynamics that motivates the chirality analysis.","marker":"[5]"}],"fun_headline_variants":["Chiral transfer near Liouvillian EP is transient","Liouvillian EP chirality seen only mid-encirclement","Photonic test shows chiral transfer fades with time","Universal scaling law for transient chiral state transfer","Chirality near Liouvillian EP vanishes at long times"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that averaging ten noisy photon trajectories, rescaled by a normalization factor obtained from the evolution operator, faithfully reproduces the true Lindblad density matrix well enough that the measured chirality and scaling reflect the physical Liouvillian dynamics rather than the simulation procedure.","fun_headline_variants_meta":{"raw":{"variants":["Chiral transfer near Liouvillian EP is transient","Liouvillian EP chirality seen only mid-encirclement","Photonic test shows chiral transfer fades with time","Universal scaling law for transient chiral state transfer","Chirality near Liouvillian EP vanishes at long times"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000234,"raw_usage":{"total_tokens":1480,"prompt_tokens":909,"completion_tokens":571,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":525,"completion_tokens_details":{"reasoning_tokens":492}},"tokens_in":525,"tokens_out":571,"duration_ms":5551,"temperature":1.0,"reasoning_tokens":492,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T19:10:58.814398+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same encircling protocol with the density matrix computed directly from the Lindblad master equation, or with a much larger stochastic ensemble, at the parameters of Fig. 4; if the chirality values and the collapse $C=f(\\gamma_2 T^{1/\\nu})$ with $\\nu\\approx1.7221$ do not survive, the observed chirality is an artifact of the ten-realization reconstruction. As a second check, encircle a loop of the same shape that does not enclose the Liouvillian exceptional point: if direction-dependent state flipping persists, the chirality is not caused by exceptional-point encircling.","supporting_citations":[{"cited_title":"Kampen, and N","cited_arxiv_id":null,"evidence_quote":"Establishes the equivalence between the Lindblad master equation and the quantum Langevin equation that the experiment uses to reconstruct the density matrix."},{"cited_title":"Sun, and W","cited_arxiv_id":null,"evidence_quote":"Gives the theory of chiral state transfer under dephasing and defines how the state trajectory is projected onto the Liouvillian spectral landscape."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Demonstrates quantum jumps in the non-Hermitian dynamics of a superconducting qubit, the dissipative-qubit context this experiment extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the chirality measure used in Eq. (11) and the decoherence-induced exceptional-point comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the general quasi-adiabatic description of exceptional-point encircling dynamics that motivates the chirality analysis."}],"review_version":1}