Cited by 13Quantum Mechanics and Non-Hermitian PhysicsQuantum Mechanics and ApplicationsQuantum Computing Algorithms and Architecture

Photonic quantum simulations of coupled PT -symmetric Hamiltonians

Open access
Maraviglia, Nicola, Yard, Patrick, Wakefield, Ross, Carolan, Jacques, Sparrow, Chris, Chakhmakhchyan, Levon, Harrold, Chris, Hashimoto, Toshikazu +4 more·2022·Physical Review Research

Abstract

Parity-time-symmetric ($\mathcal{PT}$-symmetric) Hamiltonians are generally non-Hermitian and give rise to exotic behavior in quantum systems at exceptional points, where eigenvectors coalesce. The recent realization of $\mathcal{PT}$-symmetric Hamiltonians in quantum systems has ignited efforts to simulate and investigate many-particle quantum systems across exceptional points. Here, we use a programmable integrated photonic chip to simulate a model composed of twin pairs of $\mathcal{PT}$-symmetric Hamiltonians, with each the time reverse of its twin. We simulate quantum dynamics across exceptional points including two- and three-particle interference, and a particle-trembling behavior that arises due to interference between subsystems undergoing time-reversed evolutions. These results show how programmable quantum simulators can be used to investigate foundational questions in quantum mechanics.

AI snapshot

Written from the paper’s full text.

TL;DR

The paper uses a programmable six-mode integrated photonic chip to simulate coupled PT-symmetric Hamiltonians and their time-reversed twins via unitary dilation. It observes dynamics across exceptional points, including two- and three-particle interference and a Zitterbewegung-like trembling from interference between forward- and reverse-time subsystems.

Key findings

  • The authors simulate two-mode PT-symmetric systems with eigenvalues \(\pm\sqrt{1-\gamma^2}\) and an exceptional point at \(\gamma_c=1\); they observe periodic evolution for \(\gamma=0.25\), nonperiodic steady-state behavior for \(\gamma=1.1\), and divergent oscillation period at the EP \(\gamma=1\).
  • For an initially mixed PT qubit \(\rho_0=0.5|1_F\rangle\langle1_F|+0.5|2_F\rangle\langle2_F|\), the von Neumann entropy is nonmonotonic and periodic for \(\gamma<1\), but decays to 0 for \(\gamma>1\); a no-signaling violation test is maximal at the exceptional point \(\gamma=1\).
  • Three-mode PT-symmetric systems have eigenvalues \(\{0,\sqrt{2-\gamma^2},-\sqrt{2-\gamma^2}\}\) and an EP at \(\gamma_c=\sqrt{2}\); two-photon and three-photon interference measurements show opposite-time evolution and parity-symmetry-related bunching patterns, with the exchange symmetry between modes 1 and 3 broken when \(\gamma>0\).
  • Coherent superpositions across forward and reverse subsystems produce oscillatory expectation values of \(S_{F(R)}\), interpreted as a discrete-variable analogue of Zitterbewegung; the oscillation amplitude depends on the initial coherence and vanishes for a statistical mixture.

Methods

  • The simulation embeds each nonunitary evolution operator \(G_N(\gamma,t)=e^{-iH_N(\gamma)t}\) into a \(2N\times2N\) unitary \(U_{2N}(\gamma,t)\) using Halmos unitary dilation, with a defect operator \(D_N\); the forward and reverse subsystems are the first and last \(N\) modes of the interferometer.
  • Experiments use a programmable integrated photonic chip comprising a triangular mesh of 15 Mach-Zehnder interferometers, 30 thermo-optic phase shifters, and six waveguides, implementing arbitrary six-mode unitaries via Reck decomposition.
  • Single photons are generated by spontaneous parametric downconversion in a bismuth borate crystal pumped at 404 nm, filtered at 808 nm, and detected with silicon avalanche photodiodes; one-, two-, and three-photon input states are used, with tomography, detector calibration, and probabilistic number-resolving detection for bunching events.
  • Experimental data are compared with theory curves; error bars for three-particle data are obtained via bootstrapping.

Limitations

  • The paper does not include a dedicated limitations section.
  • Evident from the text: the setup is capable of simulating coupled PT-symmetric Hamiltonians with \(N=2\) and \(N=3\) modes in each subsystem and multiparticle evolution up to three photons.
  • Evident from the text: three-photon simulations require probabilistic number-resolving detection of up to two photons in the same mode, and final probabilities are normalized over output states with three photons in the forward modes but less than three photons bunched in the same mode.
  • Evident from the text: higher-order pair generation is neglected, occurring at a two-order-of-magnitude-lower rate.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.1103/physrevresearch.4.013051

Cites

50 of 50 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Claims the paper makes

0 of 37 cited sources checked · 0 supported · 0 flagged

  • Findingown data

    The probabilities of the two alternative outputs for qubit B are conditioned by a local unitary transformation applied to the PT qubit A prior to its nonunitary evolution.

    “Figure 2(d) shows that the probabilities of the two alternative outputs for qubit B are conditioned by the a local unitary transformation applied to the PT qubit A prior to its nonunitary evolution.”
    Quote in the paperSubstantial claim
  • Findingown data

    In the unbroken symmetry phase with γ=0.25, the time evolution of the system is periodic.

    “In the unbroken symmetry phase, shown in Fig. 2(a) with γ = 0.25, the time evolution of the system is periodic, and the finite probability of detecting |vac F indicates that the excitation can tunnel from the forward system to its time-reversed twin.”
    Quote in the paperSubstantial claim
  • Findingown data

    In the unbroken symmetry phase, the von Neumann entropy of the initially mixed PT qubit is a nonmonotonic function of time, with a minimum at t=0.5τ when time is in units of τ=π/|1-γ^2|.

    “As opposed to a closed Hermitian system whose entropy is a constant of motion, in the unbroken symmetry phase, the entropy is a nonmonotonic function of time.By reporting the time parameter in units of τ = π/ |1 -γ 2 |, the entropy has a minimum when t = 0.5τ .”
    Quote in the paperSubstantial claim
  • Findingcitation

    The authors experimentally simulate multiparticle dynamics in two- and three-mode PT-symmetric Hamiltonians using a programmable photonic chip and ensembles of one-, two-, and three-photon input states.

    “We experimentally simulate multiparticle dynamics [26] in two-and three-mode PT -symmetric Hamiltonians using a programmable photonic chip [27] and ensembles of one-, two-, and three-photon input states.”
    Quote in the paperSubstantial claim0 of 2 sources checked
    • cites [26]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites [27]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • Findingreasoning

    In the Hermitian limit, due to the underlying parity symmetry, the probability of two photons in mode 1 and one photon in mode 2 (pattern |2p,1p,0p_F) is equal to that of its parity-symmetric pattern |0p,1p,2p_F.

    “Additionally, in the Hermitian limit, due to the underlying parity symmetry, the probability of two photons in mode 1 and one photon in mode 2 (i.e., pattern |2 p , 1 p , 0 p F ) is equal to that of its paritysymmetric pattern |0 p , 1 p , 2 p F .”
    Quote in the paperNot a substantial claim
  • May combine two claims. 30
  • Classified as probably not a substantial claim. 26
  • May combine two claims.

    Findingreasoning

    At the exceptional point, the probabilities evolve towards a never-reached turning point, and the system effectively tends to a steady state with coefficients that follow a polynomial power law.

    “The probabilities evolve towards a never-reached turning point, and the systems effectively tends to a steady state with coefficients that follows a polynomial power law.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    In the broken symmetry phase, the system approaches its steady state with an exponential power law whose characteristic time τ=π/√(γ^2-1) decreases for increasing γ.

    “In contrast, in the broken symmetry phase, the system approaches its steady state with an exponential power law whose characteristic time, τ = π/ γ 2 -1, decreases for increasing γ .”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    The no-signaling condition is violated for all γ=0 and is maximal for γ=γ_c.

    “The no-signaling condition is violated for all γ = 0 and is maximal for γ = γ c .”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    The authors reproduce dynamics in the PT-symmetric regime and across the exceptional point into the PT-symmetry-broken regime, including effects of mutual coherence and interference between forward-in-time and backward-in-time subsystems.

    “We reproduce the dynamics in the PT -symmetric regime (defined by real eigenvalues) and across the EP into the PT -symmetry-broken regime (defined by complex-conjugate eigenvalues), including the effects of mutual coherence and interference between the forward-in-time and backward-in-time subsystems.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    In the broken symmetry phase with γ=1.1, the system rapidly tends to a steady state with no oscillations and with a nonzero probability of finding the excitation in either the forward or reverse systems.

    “Figure 2(b) shows, instead, the evolution in the broken symmetry phase with γ = 1.1.The system rapidly tends to a steady state with no oscillations and with a nonzero probability of finding the excitation in either the forward or reverse systems.”
    Quote in the paperSubstantial claim
  • Classified as probably not a substantial claim.

    Findingown data

    The data in Figs. 4(d)-4(f) show the evolution of the probability of different photon bunching events in the forward system modes.

    “The data in Figs.4(d)-4(f) show the evolution of the probability of different photon bunching events in the forward system modes.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    The amplitude of these oscillations is reduced if the coherence between |1_F and |1_R decreases, and it vanishes if the input is a statistical mixture of |1_F and |1_R.

    “The amplitude of these oscillations is reduced if the coherence between the states |1 F and |1 R decreases, and it vanishes if the input is a statistical mixture of |1 F and |1 R [Fig.3(a)].”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingnone given

    By increasing γ, the system undergoes the transition from the symmetric to the broken symmetry phase characterized by an increase of the period of evolution for γ<√2 and by a nonperiodic evolution for γ>√2.

    “By increasing the value of γ , the system undergoes the transition from the symmetric to the broken symmetry phase characterized by an increase of the period of the evolution, for γ < √ 2, and by a nonperiodic evolution, for γ > √ 2.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    In the broken symmetry phase, S_F(R) no longer exhibits oscillatory behavior because the system loses its periodicity, but the final configuration changes greatly depending on the initial coherence between the two subsystems.

    “We finally note that in the broken symmetry phase [Fig.3(b)], S F(R) no longer exhibits oscillatory behavior since the system loses its periodicity, but the final configuration of the system changes greatly depending on the initial coherence between the two subsystems.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingreasoning

    When γ>0, the exchange symmetry between modes 1 and 3 is broken; the weight at smaller times t<T(γ)/2 shifts to the gain mode, and the weight at times t>T(γ)/2 shifts to the loss mode.

    “When γ > 0, the exchange symmetry between modes 1 and 3 is broken; the weight at smaller times t < T (γ )/2 shifts to the gain mode, and the weight at times t > T (γ )/2 shifts to the loss mode.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    The entropy of the system is periodic for γ<1 and goes asymptotically to 0 for γ>1.

    “(c) The entropy of the system is periodic for γ < 1 and goes asymptotically to 0 for γ > 1.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    The data for |ψ2 follow the opposite temporal evolution as |ψ1, while the probabilities of observing |1 and |3 are exchanged.

    “Conversely, the data for |ψ 2 follow the opposite temporal evolution as |ψ 1 , while the probabilities of observing |1 and |3 are exchanged.”
    Quote in the paperNot a substantial claim
    • The authors state this as an assumption: The production of the desired input state is assumed every time a coincidence detection happens at two detectors, neglecting higher-order pair generation occurring with a two-order-of-magnitude-lower rate.
    • Rests on an assumption: For N=3, the third-order exceptional point occurs at γ_c=√2.
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Findingown data

    The Hong-Ou-Mandel visibility of the interference of two indistinguishable single photons is 0.97.

    “The Hong-Ou-Mandel visibility of their interference is measured to be 0.97.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingown data

    In the Hermitian regime (γ=0), the vacuum probability is constant at 0 and there are sinusoidal oscillations between the two populated levels.

    “the inset shows the Hermitian regime (γ = 0) with a constant vacuum probability of 0 and sinusoidal oscillations between the two populated levels.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingreasoning

    The partial trace over the reverse (forward) system sums the diagonal elements of the bottom (top) N modes of ρ(t) and places this value in the vacuum state of the forward (reverse) system.

    “Similar results are obtained when the partial trace is performed over the forward system, showing that the partial trace operation over the reverse (forward) system sums the diagonal elements of the bottom (top) N modes of ρ(t ) and places this value in the vacuum state of the forward (reverse) system.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingown data

    In the broken symmetry phase (γ=1.1), the evolution is not periodic, and the system tends rapidly to a steady state.

    “(b) In the broken symmetry phase (γ = 1.1), the evolution is not periodic, and the system tends rapidly to a steady state; the inset shows the evolution of the system at the exceptional point.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Findingnone given

    The spikes appearing in H_1,3 and H_2,4 represent asymptotes.

    “The spikes appearing in H 1,3 and H 2,4 represent asymptotes.”
    Quote in the paperNot a substantial claim
  • Classified as probably not a substantial claim.

    Findingreasoning

    The Fock spaces of the forward and reverse systems are decomposed into particle number subspaces as F(F/R)=H(F/R)_0 ⊕ H(F/R)_1.

    “That is, the Fock spaces of the systems are decomposed into particle number subspaces as F (F/R) = H (F/R) 0 ⊕ H (F/R) 1 .”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingreasoning

    When H_N is Hermitian, G_N(t) is unitary and the singular values are 1, so D_N=0, giving two independently evolving systems.

    “In the case when H N is Hermitian, G N (t ) is unitary and the singular values are 1, so D N = 0, giving two independently evolving systems.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Findingown data

    The unplotted elements in Fig. 5 are all numerically 0, including all diagonal elements H_k,k for 1≤k≤4.

    “The unplotted elements in Fig. 5 are all (numerically) 0 including all diagonal elements H k,k for 1 k 4.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

    Premisecitation

    In classical wave systems, PT-symmetric systems have shown unidirectional invisibility, single-mode and topological lasing, topological energy transfer, mode switching, and enhanced sensitivity.

    “In the classical domain, PT -symmetric systems have shown remarkable properties such as unidirectional invisibility [5], single-mode [6,7] and topological [8,9] lasing, topological energy transfer [10], mode switching [11], and enhanced sensitivity [12,13].”
    Quote in the paperNot a substantial claim0 of 5 sources checked
    • cites [5]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites [6]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites [7]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites [8]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites [9]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • Classified as probably not a substantial claim.

    Premisecitation

    The role of thermal and quantum noise on exceptional points in semiclassical PT-symmetric systems is not yet understood.

    “For all these advances, the role of thermal and quantum noise on EPs in these semiclassical systems is not yet understood [14,15].”
    Quote in the paperSubstantial claim0 of 2 sources checked
    • cites [14]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites [15]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • May combine two claims.

    Premisecitation

    PT-symmetric Hamiltonians are generally non-Hermitian and exhibit exceptional points at which two or more eigenvalues become degenerate and the corresponding eigenstates coalesce.

    “While generally non-Hermitian, PT -symmetric Hamiltonians exhibit exceptional points (EPs) at which two (or more) eigenvalues become degenerate and the corresponding eigenstates coalesce.”
    Quote in the paperSubstantial claim0 of 2 sources checked
    • cites [1]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites [2]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • May combine two claims.

    Interpretationreasoning

    Due to PT symmetry between the input states |ψ1 and |ψ2, their evolution proceeds in opposite directions in time and presents an exchange of probabilities between the patterns |1p1p0p and |0p1p1p.

    “We note that, due to the PT symmetry between these input states, we expect their evolution to proceed in opposite directions in time and to present an exchange of probabilities between the observation of the patterns |1 p 1 p 0 p and |0 p 1 p 1 p .”
    Quote in the paperNot a substantial claim
  • Classified as probably not a substantial claim.

    Premisecitation

    Passive PT symmetry is based on the equivalence, after normalization, of a PT-symmetric Hamiltonian and a dissipative Hamiltonian with mode-selective losses.

    “The first approach, known as passive PT symmetry, is based on the equivalence (after normalization) of a PT -symmetric Hamiltonian and a dissipative Hamiltonian with mode-selective losses [17,18].”
    Quote in the paperSubstantial claim0 of 2 sources checked
    • cites [17]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites [18]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • May combine two claims.

    Interpretationnone given

    Whenever γ=0, the two systems are coupled by the defect matrix D_N, and excitations can tunnel between these two systems.

    “Whenever γ = 0, the two systems are coupled by the defect matrix D N , and excitations can tunnel between these two systems.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Study design and conductnone given

    Calibration of the on-chip thermo-optic phase shifters allows mapping a desired phase shift to a voltage bias to be set.

    “Calibration of the on-chip thermo-optic phase shifters allows us to map a desired phase shift to a voltage bias to be set.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisenone given

    With N=2, the eigenvalues of the Hamiltonian are ±√(1-γ^2).

    “With N = 2, the eigenvalues of the Hamiltonian are = ± 1 -γ 2 .”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

    Interpretationnone given

    As γ increases, the coupling between the systems iD(γ,t) increases, correspondingly decreasing the probability of a photon remaining in the desired subsystem.

    “As we increase γ , we also increase the coupling between the systems, iD(γ , t ), which correspondingly decreases the probability of a photon remaining in the desired subsystem.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Interpretationreasoning

    G†_N(γ,t)=G*_N(γ,t)=T G_N(γ,t) T^{-1}, meaning a system initialized in the last N modes transforms with a backwards-in-time evolution.

    “Indeed, G † N (γ , t ) = G * N (γ , t ) = T GN (γ , t )T -1 , meaning that a system initialized in the last N modes of the interferometer transforms with a backwards-in-time (reverse) evolution.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Premisecitation

    The spectrum of H_N(γ) is purely real for γ<γ_c, and the corresponding eigenvectors are simultaneous eigenvectors of the antilinear PT operator.

    “The spectrum of the Hamiltonian H N (γ ) is purely real for γ < γ c , and the corresponding eigenvectors are simultaneous eigenvectors of the antilinear PT operator (unbroken symmetry phase).”
    Quote in the paperSubstantial claim0 of 2 sources checked
    • cites [28]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites [29]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • May combine two claims.

    Premisecitation

    Hamiltonian dilation places the PT-symmetric Hamiltonian in a larger Hilbert space by introducing an ancillary two-level system, and the coherent nonunitary dynamics are recovered from the unitary dynamics in the larger Hilbert space.

    “The second method uses Hamiltonian dilation that places the PT -symmetric Hamiltonian in a larger Hilbert space by introducing an ancillary two-level system [20][21][22], and the coherent, nonunitary dynamics of interest are recovered from the unitary dynamics in the larger Hilbert space.”
    Quote in the paperSubstantial claim0 of 3 sources checked
    • cites [20]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites [21]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites [22]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • May combine two claims.

    Study design and conductreasoning

    All values of γ can be implemented in the same way, allowing investigation of system dynamics in four regimes: uncoupled Hermitian evolution, coupled PT-symmetric evolution, evolution at the EP, and evolution in the broken symmetry phase.

    “In this way, all values of γ can be implemented in the same way, allowing us to investigate the system dynamics in four regimes: uncoupled Hermitian evolution, coupled PT -symmetric evolution, evolution at the EP, and evolution in the broken symmetry phase.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Study design and conductcitation

    The photonic chip comprises a triangular mesh of 15 interconnected Mach-Zehnder interferometers, controlled by 30 metallic thermo-optic phase shifters, and is capable of implementing any unitary linear optical transformation over its six waveguides.

    “Our photonic chip comprises a triangular mesh of 15 interconnected Mach-Zehnder interferometers, controlled by 30 metallic thermo-optic phase shifters, and is capable of implementing any unitary linear optical transformation over its six waveguides (see Ref. [27] for additional device details).”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites [27]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • May combine two claims.

    Premisereasoning

    In passive PT symmetry, evolution time is limited by exponentially decaying signal or postselection success probability, while Hamiltonian dilation requires a time-dependent Hermitian Hamiltonian for the larger Hilbert space.

    “In the first approach, the evolution time of the system is limited by the exponentially decaying signal or postselection success probability, while the second method requires a time-dependent Hermitian Hamiltonian for the larger Hilbert space.”
    Quote in the paperNot a substantial claim
  • Classified as probably not a substantial claim. May combine two claims.

    Study design and conductreasoning

    The first section of the chip implements the state evolution described by {U_4(γ,t_i)}_i, and the remaining part performs projective measurements allowing reconstruction of density matrices for evolved single-photon states.

    “The first section of the chip implements the state evolution described by {U 4 (γ , t i )} i .The remaining part of the chip is used to perform projective measurements allowing the reconstruction of the density matrices for evolved single-photon states.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Premisecitation

    γ_c/J equals 1 whenever N is even and equals √((N+1)/(N-1)) when N is odd.

    “γ c /J = 1 whenever N is even and γ c /J = √ (N + 1)/(N -1) when N is odd.”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites [28]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • May combine two claims.

    Premisecitation

    A dissipative Hamiltonian arises naturally from a Lindblad approach for open quantum systems when certain quantum jumps are ignored and allows coherent nonunitary dynamics in a minimal quantum system.

    “Such a dissipative Hamiltonian arises naturally from a Lindblad approach for open quantum systems when certain quantum jumps are ignored and allows the implementation of coherent, nonunitary dynamics in a minimal quantum system [19].”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites [19]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    The dilation process has already been used to simulate Lindbladian dynamics.

    “This process has already been used to simulate Lindbladian dynamics [32] and could be extended to simulate anti-PT systems [33].”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites [32]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • Probably not a substantial claim: may not pass S2 (checkable).

    Interpretationown data

    These demonstrations point towards the inclusion of non-Hermitian simulations in the realm of quantum technologies.

    “Combined with the current advances in photonic technologies, these demonstrations point towards the inclusion of non-Hermitian simulations in the realm of quantum technologies.”
    Quote in the paperNot a substantial claim
  • Classified as probably not a substantial claim.

    Study design and conductown data

    The relative detector efficiency at different output modes for single-photon events is obtained by testing 500 random phase configurations and comparing experimental statistics with expected photon statistics.

    “The relative detector efficiency at the different output modes, for single-photon events, is obtained by testing 500 random phase configurations of our interferometer and comparing the experimental statistic with the expected photon statistics due to the phase settings.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Study design and conductnone given

    The final probabilities are normalized over all output states characterized by three photons in the forward system modes but with less than three photons bunched in the same mode.

    “The final probabilities are normalized over all the output states characterized by three photons in the forward system modes but with less than three photons bunched in the same mode.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    In boson sampling, classical intractability arises because the detection statistics are given by permanents of complex matrices.

    “In the case of boson sampling [46], in which photons propagate in a circuit that is ideally described by a unitary matrix, classical intractability arises because the detection statistics are given by permanents of complex matrices.”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites [46]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • Classified as probably not a substantial claim.

    Study design and conductnone given

    The efficiency correction of each mode that maximizes statistical fidelity between experimental statistics and theoretical distribution is chosen as numerical correction for single-photon data.

    “The efficiency correction of each mode that maximizes the statistical fidelity between the experimental statistics and the theoretical distribution is chosen as a numerical correction for the data collected during single-photon experiments.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    Photon loss reduces the computational complexity in experimental implementations of boson sampling.

    “Yet photon loss reduces the computational complexity in experimental implementations of boson sampling [47].”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites [47]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • Classified as probably not a substantial claim.

    Study design and conductnone given

    For two-photon antibunched events, the occurrence of each pattern is multiplied by the product of the correction factor of the individual output modes populated with photons.

    “For two-photon antibunched events, the occurrence of each pattern is multiplied by the product of the correction factor of the individual output modes populated with photons.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Study design and conductcitation

    Error bars are obtained via bootstrapping methods.

    “Error bars are obtained via bootstrapping methods [45].”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites [45]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    The same device can be used as a cascade of 2×2 interferometers to sequentially implement state preparation, evolution, or projection in physically different parts of the chip.

    “Alternatively, the same device can be used as a cascade of 2 × 2 interferometers in order to sequentially implement state preparation, evolution, or projection in physically different parts of the chip.”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites [27]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • Probably not a substantial claim: may not pass S3 (load-bearing), S4 (based).

    Study design and conductnone given

    The number of events associated with an output state |N1p N2p N3p with N_kp = n_k,1 + n_k,2 is given as the average of the occurrences of equivalent detection patterns.

    “The number of events associated with an output state |N 1p N 2p N 3p with N k p = n k,1 + n k,2 is given as the average of the occurrences of equivalent detection patterns.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    It is computationally complex to simulate the evolution of quantum Hamiltonians that produce particle statistics governed by intractable matrix functions.

    “It is computationally complex to simulate the evolution of quantum Hamiltonians [32] that produce particle statistics which are governed by intractable matrix functions.”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites [32]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Interpretationreasoning

    A photonic simulation of a traditional model for an open system, such as the Lindblad approach, could become classically tractable via the loss that is necessarily built into the model.

    “Similarly, a photonic simulation of a traditional model for an open system, such as the Lindblad approach, could become classically tractable via the loss that is necessarily built into the model.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Study design and conductnone given

    In the pseudo-number-resolving configuration, the detection calibration factor is obtained by using the single-photon detector efficiency ν_i recorded at each output of the FBS.

    “In the pseudo-number-resolving configuration, the detection calibration factor is obtained by using the single-photon detector efficiency ν i recorded at each output of the FBS.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Interpretationnone given

    In the authors' simulation model for twin PT-symmetric Hamiltonians, while each subsystem is open and lossy, the overall model is unitary and an ideal implementation is lossless; generally, statistics provide information about both the forward and backward subsystems and the relation between them.

    “However, in our simulation model for twin PT -symmetric Hamiltonians, while each subsystem is open and lossy, the overall model is unitary and an ideal implementation is lossless; generally, statistics provide information about both the forward and backward subsystems, and the relation between them.”
    Quote in the paperSubstantial claim
53 of 53
  • The paper relies on this without support.

    Most approximate methods for many-body systems (variational principle, tensor networks, perturbation theory) are designed for and work only with Hermitian Hamiltonians.

    “most approximate methods-variational principle, tensor networks, and perturbation theory-developed to reduce the exponential complexity of a many-body system, are designed for and work only with Hermitian Hamiltonians.”
  • The paper relies on this without support.

    Non-Hermitian quantum many-particle systems are an emerging and challenging frontier.

    “non-Hermitian quantum many-particle systems are an emerging and challenging frontier.”
  • The authors state this as an assumption.

    The production of the desired input state is assumed every time a coincidence detection happens at two detectors, neglecting higher-order pair generation occurring with a two-order-of-magnitude-lower rate.

    “The production of the desired input state is assumed every time a coincidence detection happens at two of the detectors connected to the output modes of our device neglecting higher-order pair generation occurring with a twoorder-of-magnitude-lower rate.”
  • The paper relies on this without support.

    For N=3, the third-order exceptional point occurs at γ_c=√2.

    “Since, for N = 3, the third-order exceptional point occurs at γ c = √ 2, to analyze unbroken and broken symmetry regimes, we perform simulations for γ = √ 2/2 and γ = 1.1 √ 2, respectively.”
4 of 4

Cites

Flagged first, then by citations

1 of the 51 references couldn’t be matched to a known paper.

Research works printed but not linked: 4

  1. P. R. Halmos, 1950, Summa Bras. Math., 2, 125
  2. E. Schrödinger, 1930, Sitzungber. Preuss. Akad. Wiss. Phys.-Math. Kl., 24, 418
  3. G. Simpson, 1986, Astron. Astrophys., 162, 340
  4. 10.4086/toc.2013.v009a004 · DOI

Cited by

Sorted by citations

Authors
Maraviglia, N., Yard, P., Wakefield, R., Carolan, J., Sparrow, C., Chakhmakhchyan, L., Harrold, C., Hashimoto, T., Matsuda, N., Harter, A. K., Joglekar, Y. N., & Laing, A.
Affiliations
  • University of Bristol (GB)
  • University of Copenhagen (DK)
  • NTT (Japan) (JP)
  • Indiana University – Purdue University Indianapolis (US)
Venue
Physical Review Research
Published
2022-01-26
Volume
4
Issue
1
ISSN
2643-1564 (ISSN-L)
OpenAlex
W4210411560
Cited by
13
References
51
Citations by year
2026: 1 · 2025: 1 · 2024: 5 · 2023: 4 · 2022: 2
FWCI
1.21
Type
journal article
Language
en
Open access
Locations
  • Physical Review Research (publishedVersion, cc-by)
  • Research at the University of Copenhagen (University of Copenhagen) (publishedVersion, cc-by)
  • arXiv (Cornell University) (submittedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion)
  • Figshare (submittedVersion, cc-by)
  • Bristol Research (University of Bristol) (submittedVersion, cc-by)
  • IUScholarWorks (Indiana University) (submittedVersion, cc-by)
  • arXiv (Cornell University) (cc-by)
Funders
  • Engineering and Physical Sciences Research Council (EP/T001062/1, EP/N003470/1)
  • National Science Foundation (DMR-1054020)
  • Institute of Advanced Studies, University of Bristol
Topics
  • Physical Sciences › Physics and Astronomy › Atomic and Molecular Physics, and Optics › Quantum Mechanics and Non-Hermitian Physics
  • Physical Sciences › Physics and Astronomy › Atomic and Molecular Physics, and Optics › Quantum Mechanics and Applications
  • Physical Sciences › Computer Science › Artificial Intelligence › Quantum Computing Algorithms and Architecture
Keywords
PT-symmetric Hamiltonians, non-Hermitian quantum mechanics, exceptional points, time-reversal symmetry, quantum simulation, quantum foundations

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Claims read by deepseek-v4.1-flash under definition version 3 (extract-v3) on 2026-10-23.

Classified by jev-1.13-20260917 (gold-guideline-v0.2) on 2026-10-23.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review