{"id":"0c1971cc-66e3-4f05-8c95-0491da80d3e1","arxiv_id":"2608.01484","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"An asymmetric loss cut in a Bloch-oscillating waveguide array filters N=1 vs N=2 N00N states, with phase-controlled conditional probabilities of 80% and 90%.","lead":"A simulation shows that a photonic waveguide array can filter quantum light into mostly single-photon or two-photon states by cutting off half the waveguides at a specific point. By tuning an input phase, the device switches between outputting 80% two-photon states and 90% single-photon states, a step toward better photon state preparation in future photonic circuits.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Inconsistent N=1 phase factor in Eq. (9) vs Eq. (5) may invalidate the reported phase-tuning results; verify which initial state was simulated.","rationale":"I chose the phase inconsistency in Eq. (9) over the cut idealization because it is a concrete internal error that directly affects the central quantitative claim. The cut model is a known approximation that the authors acknowledge and bound (finite ϵ, sufficiently long Bloch period), and its effect is a matter of experimental realization rather than a flaw in the theoretical argument. In contrast, the phase factor error is an internal inconsistency: if the simulations used Eq. (9) as written, then the reported phase-switching results were obtained for an input state different from the one the paper defines in Eq. (5), and the claimed optimal phases (φ=π, φ≈π/2,3π/2) would not be correct for the stated state. Even if the error is only typographical and the simulations used the correct phase, the manuscript as written is not reproducible without resolving this discrepancy. This warrants the same CONDITIONAL verdict the Reader assigned: the idea is plausible, but the paper must correct the phase factor and verify that the results are unchanged. The Reader's weakest assumption (cut model) is a legitimate limitation, but the phase issue is more load-bearing because it can invalidate the headline numbers outright. Thus I partially agree with the Reader: we both flag the phase problem, but the Reader identifies the cut model as the weakest point, whereas I see the phase inconsistency as the primary concern. I recommend no change to the CONDITIONAL verdict, hence UNCHANGED.","tokens_in":8894,"tokens_out":5497,"duration_ms":47157,"concrete_test":"Re-implement the simulation with the N=1 sector phase factor e^{iφ} (as in Eq. 5), not e^{2iφ} (as written in Eq. 9), using the same parameters: C/B=0.3, α=β=1, 6 waveguides, cut on left half at z_c=λ_B/2, output at z=λ_B, with photon-number truncation as in the paper. Compute P1(φ) and P2(φ) and compare to Figs. 3–5. If the optimal phases shift (e.g., 2-photon filtering at φ=π/2 instead of φ=π) or the 80%/90% values change materially, then the paper's simulations used the wrong initial state and the central quantitative claim is invalid as stated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The manuscript's central quantitative claim—phase-controlled switching between 2-photon dominated (80%) and 1-photon dominated (90%) outputs—rests on the correct phase dependence of the N=1 N00N state Bloch oscillations. However, Eq. (9) writes the N=1 term as α(â†_μ(−z) + e^{2iφ} â†_ν(−z)), whereas the declared input state Eq. (5) has the N=1 term α(|1_μ0_ν⟩ + e^{iφ}|0_μ1_ν⟩). This means the evolution expression does not represent the stated initial state. If the MATLAB simulations were coded from Eq. (9), then all phase-dependent results, including the claimed optimal phases (φ=π for 2-photon filtering, φ≈π/2,3π/2 for 1-photon filtering) and the associated 80%/90% efficiencies, were computed for a different input state (N=1 phase 2φ instead of φ). A simple substitution shows that the N=1 fringe localization at the half Bloch period, which the cut exploits, would occur at φ=π/2 rather than φ=π when the phase is doubled, so the reported switching points would be off. The paper does not clarify whether the simulation used Eq. (5) or Eq. (9). This is a concrete internal inconsistency, not merely a philosophical concern about the cut idealization. The cut model is an acknowledged approximation; the phase error, if real, directly undermines the headline numbers.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies Bloch oscillations of N=1 and N=2 photonic N00N states in an array of evanescently coupled waveguides with a linearly graded propagation constant. It proposes placing total loss on one half of the waveguide array at the half Bloch period and simulates the resulting evolution with a reduced density matrix. The central claim is that tuning the input phase phi dynamically switches the output between a 2-photon-dominated subspace (about 80% conditional probability) and a 1-photon-dominated subspace (about 90% conditional probability). The paper includes an analytic Green's-function framework, density-matrix equations for the cut, and appendices with partial analytical checks.","tokens_in":9242,"tokens_out":9305,"duration_ms":90401,"significance":"If correct, the proposed scheme is a concrete, passive approach to photon-number filtering and routing that does not rely on destructive detection. The physical mechanism, based on the phase-dependent transverse localization of the N=1 N00N sector at the half Bloch period while the N=2 sector remains symmetric, is clear and falsifiable. The manuscript also provides useful analytical scaffolding and a plausible experimental parameter range. However, the quantitative claims currently rest on an internal phase-factor inconsistency that must be resolved before the results can be accepted.","major_comments":[{"comment":"The initial state Eq. (5) defines the N=1 sector as alpha(|1_mu 0_nu> + e^{i phi}|0_mu 1_nu>), but Eq. (9) writes the evolved N=1 term as alpha( a_dagger_mu(-z) + e^{2 i phi} a_dagger_nu(-z) )|0>. These differ: at phi=pi, Eq. (5) gives e^{i phi}=-1 while Eq. (9) gives e^{2 i phi}=+1. Since the cut simulations are stated to use Eq. (9), the reported phi=pi optimum for 2-photon filtration and the phi ~ pi/2, 3pi/2 optimum for 1-photon filtration may have been computed for a state whose N=1 phase is doubled. The authors must state which expression was implemented, correct Eq. (9) if it is a typo, and rerun or adjust Figs. 3-4 and the derived efficiencies accordingly. If the code used the erroneous factor, the headline 80%/90% numbers do not follow from the declared input state.","section":"III.1, Eq. (9); Figs. 3-4"},{"comment":"The cut model traces out the left modes at z_c and then reintroduces them as vacuum modes for subsequent evolution under the same Hamiltonian. This is not equivalent to a 'total loss' boundary: after z_c, light can tunnel back into the supposedly cut waveguides and return to the output modes. The paper acknowledges the idealized nature of the cut, but the repopulation effect can change the conditional probabilities in Figs. 3-5. Please provide an estimate or a separate simulation with the cut modes removed (or with explicit loss terms) to show that the filtering numbers are robust, or clearly state this limitation in the abstract and conclusions.","section":"III.1, Eqs. (12)-(13)"}],"minor_comments":[{"comment":"The sum over nu from -N/2 to N/2 contains N+1 terms for even N, inconsistent with the text 'sum over N waveguides' and with the 6-waveguide simulation. Please clarify the indexing used.","section":"II, Eq. (1)"},{"comment":"The analytical N=1 checks effectively treat the phase-symmetric state (phi=0) and do not include the e^{i phi} factor of Eq. (5). Since the phase dependence is the central mechanism, it would be valuable to show the analytical formula for general phi, or at least state the phi=0 restriction explicitly.","section":"Appendix B and C"},{"comment":"The notation sum_{n_i<1} is unclear; presumably the trace should sum over all Fock states of the cut modes. Please write the sum explicitly over n_{-N/2},...,n_0.","section":"III.1, Eq. (12)"},{"comment":"Reference [15] has a malformed DOI ('10.1103/grwm-1kj3'); please verify the bibliographic data.","section":"References"},{"comment":"The statement 'N=2 transmission is currently ~40%' can be confused with the 80% conditional filtering claim. Please define whether the 40% is absolute transmission efficiency and clarify the relation to the conditional probabilities P1 and P2.","section":"Conclusion"}],"recommendation":"major_revision","confidential_remarks":"The phase-factor issue in Eq. (9) is likely a typographical error, but it sits in the equation that is explicitly used for the cut simulations. I would ask the authors to confirm the exact form implemented in MATLAB and to make the corrected equation and figures consistent. If the numerical code already used the correct phase, the manuscript could become acceptable after a careful revision; if not, the quantitative conclusions need to be recomputed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nQuick take: the filtering idea is sound and the paper is worth a serious look, but there is a load-bearing phase-factor inconsistency between Eq. (5) and Eq. (9) that needs fixing before the quantitative claims can be trusted.\n\nWhat is new: previous work (Bromberg et al. and others) characterized the phase-dependent asymmetry of N=1 vs N=2 N00N Bloch oscillations. This paper converts that asymmetry into a concrete photon-number filter: an asymmetric loss cut at the half Bloch period, with phase-controlled switching. The density-matrix formalism, the conditional probability definitions, and the SiN waveguide parameters are all reasonable. The authors also acknowledge the input-state fidelity challenge and the idealized nature of the cut. That honesty is welcome.\n\nThe problem: Eq. (5) defines the N=1 input as α(|1μ0ν⟩ + e^{iφ}|0μ1ν⟩), but Eq. (9) evolves that sector as α(â†_μ(-z) + e^{2iφ} â†_ν(-z)). Doubling the phase in the evolution shifts the interference fringes: if the MATLAB code used Eq. (9), then the phases reported in Fig. 4 (φ=π for 2-photon filtering, φ≈π/2,3π/2 for 1-photon filtering) are not the physical input phases of Eq. (5); they correspond to φ=π/2 and π/4,3π/4. The paper never says which expression was implemented. This is not a nitpick: the 80%/90% switching claims depend on the correct phase mapping. The appendices add to the confusion by omitting the phase entirely in the N=1 analytical expressions (B1, C1).\n\nThe cut model is idealized—instantaneous total loss, no reflection, vacuum re-injection—but the authors acknowledge this and suggest future work with photonic crystals and partial loss. Minor.\n\nSummary: this is a promising proposal for photonic state routing and number filtering, with standard math and honest limitations. The phase inconsistency is fixable, but it is load-bearing. It deserves peer review, with a request to correct the phase factor and state explicitly which initial state was simulated. I would not cite it in its current form.\n\nBest","headline":"The filtering idea is sound, but the phase inconsistency between Eq. (5) and Eq. (9) undercuts the headline switching claims until clarified.","tokens_in":9722,"tokens_out":8336,"would_cite":false,"duration_ms":71475,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V80"],"pacs":["42.50.-p"],"model":"deepseek-v4-flash","headline":"A Bloch-oscillating waveguide array with an asymmetric loss cut can sort N=1 and N=2 N00N states, switching between ~80% two-photon and ~90% one-photon output by tuning input phase.","keywords":["N00N states","Bloch oscillations","photon number filtering","asymmetric loss","waveguide arrays","non-Hermitian photonics","phase-dependent switching","quantum state preparation"],"falsifier":"Measure the output photon-number statistics of an integrated waveguide array with a physically implemented cut (for example a Bragg scatterer or siphoned waveguides) at $z=\\lambda_B/2$; if the $P_1/P_2$ ratio does not track the input phase with peaks near $\\phi=\\pi$ and $\\phi\\approx\\pi/2,3\\pi/2$, or if reflected light from the cut repopulates the supposedly removed sector, the idealized trace-and-vacuum-reinsertion model fails.","tokens_in":1796,"feed_emoji":"💡","tokens_out":1727,"duration_ms":64232,"temperature":0.7,"pith_summary":"The paper proposes a way to separate the one- and two-photon parts of a photonic state without destroying it. In an array of waveguides whose effective index grows linearly, a one-photon N00N state (an entangled superposition of |1,0> and |0,1>) Bloch-oscillates with its bright fringe on one side fixed by input phase, while a two-photon N00N state spreads symmetrically on both sides. Removing half the waveguides at half a Bloch period deletes whichever sector sits on the lossy side; changing input phase switches between about 80% two-photon and 90% one-photon conditional output. Such a reconfigurable, non-destructive photon-number filter would be useful for preparing and routing states in photonic circuits.","feed_headline":"Bloch oscillations filter photon numbers: 80% two, 90% one","feed_subtitle":"Siphoning half a waveguide array at half a Bloch period makes input phase pick the 1- or 2-photon output.","key_machinery":"The argument is carried by the Green's function $U_{\\mu,\\nu}(z)=e^{i\\pi(\\mu-\\nu)/2}e^{iBz(\\mu+\\nu)/2}J_{\\mu-\\nu}\\big((4C/B)\\sin(Bz/2)\\big)$, which gives the single-photon propagation amplitudes and hence the photon-density profiles for N00N inputs. The cut is modeled as a partial trace over the modes on the lossy half of the array at $z_c=\\lambda_B/2$, followed by reintroducing those waveguides as vacuum modes and evolving under the unchanged Hamiltonian to $z_{\\text{out}}=\\lambda_B$. The ratio $C/B$, proportional to $\\lambda_B/L_t$, controls the transverse spread of the bright fringes and is the main optimization parameter.","core_discovery":"The paper claims that at the half Bloch period $z_c = \\lambda_B/2$, the N=1 sector of a N00N input localizes on one side of the waveguide array according to the input phase $\\phi$, while the N=2 sector spreads evenly across both sides. Imposing total loss on one half of the array at that point removes whichever sector sits on the lossy side. With equal N=1 and N=2 inputs, $\\phi=\\pi$ yields a conditional output that is more than 80% two-photon, while $\\phi\\approx\\pi/2,3\\pi/2$ yields roughly 90% one-photon output. The output is read from the two central waveguides at the full Bloch period $z_{\\text{out}}=\\lambda_B$, and the filtering is optimized by the ratio $C/B$ of coupling to index gradien","pith_inferences":["Because the filter exploits phase coherence between N=1 and N=2 sectors, a mixed initial state will degrade it; quantifying the tolerable phase noise or mixing would be a direct extension of the paper's density-matrix method.","The same trace-and-reinsert cut idea could be applied to staggered or partial loss profiles and to N>2 photon-number sectors, where numerical simulation is already the route the paper anticipates.","The appendix's beam-splitter coherent-state example implies this is a quantum-state filter, not a classical intensity filter; a natural application is cleaning or routing heralded non-classical photon sources.","A testable extension is to place output couplers on both central waveguides and measure photocount correlations; the predicted phase-dependent switch between one- and two-photon statistics should appear as a corresponding change in coincidence rates."],"forward_implications":["For equal N=1/N=2 inputs with $C/B=0.3$ and $\\phi=\\pi$, the conditional output is more than 80% two-photon; tuning $\\phi$ near $\\pi/2$ or $3\\pi/2$ gives roughly 90% one-photon output.","The filter is most effective when the input is N=1-dominated: with 90% N=1 input the relative two-photon efficiency stays high while the one-photon sector is strongly suppressed, at the cost of about 94% total intensity loss.","Only $C/B$ matters for the normalized dynamics, so the same design curves apply across different $B$ and $C$ values; a SiN/SiO2 platform with $L_t=200\\ \\mu$m, $C=0.008\\ \\mu$m$^{-1}$, $B=0.027\\ \\mu$m$^{-1}$ and $\\Delta n_{\\text{eff}}\\sim0.007$ is within a feasible parameter range.","The scheme is a phase-controlled switch between photon-number subspaces distributed over two output waveguides, offering demultiplexing and non-destructive routing rather than detection-triggered state preparation."],"supporting_citations":[{"why":"Supplies the Green's function solution and the key observation that N=1 N00N Bloch oscillations behave as coherent states with phase-localized fringes while N>1 sectors spread symmetrically.","marker":"[13]"},{"why":"Establishes photon-number dependent interference profiles in waveguide Bloch oscillations, which the paper extends to a filtering scheme.","marker":"[12]"},{"why":"Provides the concentric-arc waveguide geometry and the parameter mapping $\\Delta n_{\\text{eff}}=n_{\\text{eff}}d/R$ used to argue experimental feasibility.","marker":"[14]"},{"why":"Provides the context of nontrivial N00N-state interference in periodic structures that motivates the Bloch-oscillation platform.","marker":"[11]"},{"why":"Supplies the low-intensity coherent-state input model used in appendix A to show that the filter does not separate photon numbers from classical state components.","marker":"[16]"}],"fun_headline_variants":["Phase switches N00N output: 80% two-photon or 90% one","Siphoning half array at half Bloch period filters photons","N00N Bloch oscillations filter 1 vs 2 photons with phase","Selective photon number output via Bloch oscillations","Two-photon 80% or one-photon 90% by input phase"],"cache_read_input_tokens":11392,"weakest_assumption_plain":"The filtering effect is contingent on modeling the cut as instantaneous total loss in half the waveguides: the lost modes neither reflect nor re-couple, and the surviving array evolves under the same Hamiltonian; a real siphoning or Bragg-scattering cut of finite length could behave differently.","fun_headline_variants_meta":{"raw":{"variants":["Phase switches N00N output: 80% two-photon or 90% one","Siphoning half array at half Bloch period filters photons","N00N Bloch oscillations filter 1 vs 2 photons with phase","Selective photon number output via Bloch oscillations","Two-photon 80% or one-photon 90% by input phase"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000806,"raw_usage":{"total_tokens":3348,"prompt_tokens":690,"completion_tokens":2658,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":434,"completion_tokens_details":{"reasoning_tokens":2565}},"tokens_in":434,"tokens_out":2658,"duration_ms":17141,"temperature":1.0,"reasoning_tokens":2565,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T00:06:58.479397+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the output photon-number statistics of an integrated waveguide array with a physically implemented cut (for example a Bragg scatterer or siphoned waveguides) at $z=\\lambda_B/2$; if the $P_1/P_2$ ratio does not track the input phase with peaks near $\\phi=\\pi$ and $\\phi\\approx\\pi/2,3\\pi/2$, or if reflected light from the cut repopulates the supposedly removed sector, the idealized trace-and-vacuum-reinsertion model fails.","supporting_citations":[{"cited_title":"Bromberg, Y","cited_arxiv_id":null,"evidence_quote":"Supplies the Green's function solution and the key observation that N=1 N00N Bloch oscillations behave as coherent states with phase-localized fringes while N>1 sectors spread symmetrically."},{"cited_title":"Bromberg, Y","cited_arxiv_id":null,"evidence_quote":"Establishes photon-number dependent interference profiles in waveguide Bloch oscillations, which the paper extends to a filtering scheme."},{"cited_title":"Lebugle, M","cited_arxiv_id":null,"evidence_quote":"Provides the concentric-arc waveguide geometry and the parameter mapping $\\Delta n_{\\text{eff}}=n_{\\text{eff}}d/R$ used to argue experimental feasibility."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the context of nontrivial N00N-state interference in periodic structures that motivates the Bloch-oscillation platform."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the low-intensity coherent-state input model used in appendix A to show that the filter does not separate photon numbers from classical state components."}],"review_version":1}