REVIEW 3 major objections 6 minor 34 references
Quantum Phase Transition of Non-Hermitian Systems using Variational Quantum Techniques
T0 review · 3 major / 6 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read The non-Hermitian variational quantum eigensolver, applied to the transverse Ising model with real and imaginary transverse fields, reproduces exact-diagonalization magnetization on five qubits and gives a finite-size susceptibility peak…
desk verdict A modest but honest proceedings paper: known algorithm applied to a new non-Hermitian Ising model, with a real but fixable gap in specifying which eigenstate the variational search targets. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the operator pair $M_\pm(E; H) = (H^\dagger - E^*)(H - E)$ and $(H - E)(H^\dagger - E^*)$, which are Hermitian and positive semi-definite even when $H$ is not. Their expectation values $C_\pm$ vanish exactly when the parametrized state $|\phi(\theta)\rangle$ and the complex shift $E$ form a right- or left-hand eigenpair, respectively. Minimizing $C_+$ with a fully expressive circuit built from single-qubit rotations and CNOTs turns the non-Hermitian eigenproblem into a VQE-style optimization, with a modified multi-stage gradient strategy adapted from Ref. [25]. The paper also uses the decomposition of any matrix into Pauli strings to feed $H$ and $H^\dagger$ into the cost function, and uses the transverse magnetization $\langle M_x \rangle$ and the zero-temperature susceptibility $\chi_x \approx \langle M_x^2 \rangle - \langle M_x \rangle^2$ as the order parameters for phase-transition searches.
What would settle it
Run the same non-Hermitian VQE on actual five-qubit hardware (or a shot-level simulator) for the real-field Ising chain near $\Gamma = 1$ and compare the resulting $\langle M_x \rangle$ and $\chi_x$ with the paper's exact-diagonalization curves; agreement outside the $\pm 0.04$ uniform band would show the noise model is not representative. Alternatively, compute the shot-noise variance of $\langle M_x \rangle$ under the cost function; if it is inconsistent with a uniform bound, the noiseless-plus-uniform model is falsified.
Extended reading notes
Core claim
On its own terms, the paper's central claim is that the non-Hermitian variational quantum eigensolver of Ref. [25] can serve as a practical probe of quantum phase transitions in non-Hermitian spin systems. The authors construct cost functions $C_\pm$ from the positive semi-definite operators $M_\pm = (H^\dagger - E^*)(H - E)$ and $(H - E)(H^\dagger - E^*)$, minimize $C_+$ with an over-expressive ansatz, and add a uniform noise of $\pm 0.04$ to every expectation value to mimic a roughly 1000-shot measurement. With five spins, the resulting $\langle M_x \rangle$ values agree with explicit diagonalization for both the real-field Ising chain and the non-Hermitian chain with an imaginary transverse field. The susceptibility $\chi_x$ computed from exact diagonalization develops a peak that moves with system size and is consistent with the known transition at $\Gamma = 1$, while for the imaginary-field model the size trend of $\langle M_x \rangle$ points toward zero in the thermodynamic limit, suggesting no ground-state quantum phase transition. The paper explicitly leaves the imaginary-field conclusion tentative, noting that larger systems are needed.
Load-bearing premise
The load-bearing premise is that adding a uniform random error in $(-0.04, 0.04)$ to every expectation value adequately represents the noise of a real roughly-1000-shot quantum measurement, even though the paper concedes this does not capture the complexity of true noise.
Editorial extensions
If this is right
- If the five-qubit agreement carries over to larger circuits, the non-Hermitian VQE gives a concrete quantum route to the ground-state properties of sign-problem Hamiltonians, not just the two Ising variants tested here.
- The susceptibility peak in the real-field model reproduces the known $\Gamma = 1$ transition from finite-size data, so $\chi_x$ can be used to extrapolate transition locations in larger non-Hermitian systems.
- The even-odd spin alternation in $\langle M_x \rangle$ for the imaginary-field model is a finite-size artifact that appears to vanish in the thermodynamic limit, supporting a zero-order-parameter ground state.
- Adapting the optimization strategy (adding a parameter stage, using a stochastic optimizer) preserves accuracy under the paper's noise model, which is a direct requirement for any NISQ implementation.
Reading between the lines
- This suggests a sharper test: rerun the same five-qubit pipeline with actual shot-noise sampling instead of uniform perturbations; if the cost-function noise is not approximately uniform, the reported agreement will not transfer to hardware.
- The absence of a transition is established only for the $\langle M_x \rangle$ order parameter; a transition might still appear in the complex spectrum itself, such as a gap closing or exceptional-point structure, which the paper does not scan.
- The same $M_+$ construction applies to $\mathcal{PT}$-symmetric chains with staggered complex fields, so the method could be extended to the class of models studied in earlier numerical diagonalization work.
- A definitive thermodynamic-limit answer for the imaginary-field model would come from pushing exact diagonalization or the statevector simulator to roughly $N = 20$ and extrapolating the odd-spin $\langle M_x \rangle$ curve; the paper's current five-spin quantum data alone cannot separate a true zero from a slow decay.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates the use of the variational quantum algorithm of Xie, Xue, and Zhang (Ref. [25]) for computing eigenpairs of non-Hermitian Hamiltonians, with the goal of studying quantum phase transitions. Two one-dimensional spin models are considered: the transverse Ising model with a real transverse field (Eq. 2) and a non-Hermitian variant with an imaginary transverse field (Eq. 3). The authors use a statevector simulator of a five-qubit circuit with a uniform noise perturbation on expectation values, compare the resulting <M_x> values with exact diagonalization, and then use exact diagonalization for larger system sizes to examine finite-size behavior. They report agreement for N=5 and use the susceptibility peak to identify the real-field transition near Gamma=1, while for the imaginary-field model the order parameter appears to vanish with increasing odd system size, suggesting no quantum phase transition.
Significance. If the central feasibility claim were fully established, this would be a valuable step toward using quantum computers for non-Hermitian systems and sign-problem-affected lattice field theories. The paper's strengths include the direct comparison against exact diagonalization as an external benchmark, the use of the established algorithm from Ref. [25], and a conservative interpretation of the imaginary-field results. The physics conclusion for the real-field Ising model is consistent with known results. However, the lack of a defined eigenstate-selection rule and the ad hoc noise model leave the main claim only partially supported.
major comments (3)
- [Sec. 4, Eq. (7)] The cost function C+(theta,E) vanishes for every right eigenpair (E,|phi>) of the non-Hermitian Hamiltonian, not just for the ground state. Since non-Hermitian eigenvalues are generally complex, there is no natural variational ordering analogous to the minimum energy in Hermitian VQE. The manuscript states in Sec. 4 that a 'slightly modified version' of the two-stage optimization from Ref. [25] is used, but it never specifies how the algorithm selects a particular eigenpair or how it tracks the physical ground state as Gamma or Gamma_I is varied. The order parameter <M_x> in Eq. (4) is defined with respect to 'the ground state |phi0>', but for the non-Hermitian model it is not defined which eigenstate this is (e.g., lowest real part, largest overlap with the Hermitian Gamma->0 ground state, or some other rule). Without this definition and a demonstrated selection rule, the N=5 agreement with exact diagonalization could reflect convergence to different eigenpairs at different parameter values, and the central claim that the algorithm can 'quantify and study' quantum phase transitions is not yet supported. Please specify the selection mechanism and provide evidence that a single physical state is tracked across the transition.
- [Sec. 5] The noisy-quantum-simulation results rest on an ad hoc noise model: each expectation value is perturbed by a uniform random variable in (-0.04,0.04), claimed to be comparable to about 1000 shots. This model does not reproduce the binomial statistics of actual measurement shots, does not include gate errors or decoherence, and the manuscript itself concedes that it 'does not capture the complexity of true noise.' Moreover, no error bars are shown on the noisy data points in Fig. 1, so the reader cannot judge whether the observed agreement with exact diagonalization is statistically meaningful. As a result, the statement in Sec. 7 that the method works 'in a noisy environment' is not substantiated. Please replace this with shot-based sampling or at least show statistical uncertainties, and temper the noisy-environment claim accordingly.
- [Sec. 6.2] The conclusion that the imaginary-field Ising model has no quantum phase transition is based on the observation that <M_x> for odd numbers of spins appears to tend to zero as N increases (Fig. 3). This is a qualitative reading of the plot; there is no quantitative finite-size scaling analysis, no extrapolation to the thermodynamic limit, and no discussion of the even-odd discrepancy beyond citing boundary effects. Since a negative claim (absence of a transition) requires more care than a positive one, the conclusion 'there is no quantum phase transition' is premature. Please add a quantitative scaling analysis or explicitly reframe the conclusion as a provisional statement.
minor comments (6)
- [Abstract and Sec. 1] The abstract says PT-symmetry is discussed, but the two models studied (Eqs. 2-3) do not respect PT-symmetry; consider making the connection between the introductory discussion and the actual models more explicit.
- [Sec. 4] The 'slightly modified version' of the optimization strategy is not described in enough detail for reproduction. Please give the number of stages, the parameter update rules, and the convergence criteria used.
- [Sec. 5] Typo: 'enviroment' should be 'environment'.
- [Sec. 6.1] The phrase 'on the hardware we used' is misleading because the results were obtained from a statevector simulation, not quantum hardware.
- [Fig. 2] The caption says quantum-simulator results are also reported for five spins, but the figure does not distinguish them from the exact-diagonalization points; add distinct markers and a legend.
- [Eq. (5)] The notation 'approx' for the susceptibility is vague; it should be stated that this is a zero-temperature fluctuation measure rather than a standard thermodynamic response function.
Circularity Check
No circularity: the variational results are checked against independent exact diagonalization, and the QPT conclusions rest on exact diagonalization, not on fitted inputs or self-citations.
full rationale
The paper's derivation chain is self-contained in the relevant sense. The non-Hermitian variational algorithm is taken from independent prior work (Ref. [25]); the cost functions C_+ and C_- (Eq. 7) are constructed from the Hamiltonian and a complex parameter E, and the paper checks the resulting eigenpairs against explicit diagonalization with the C++ Eigen library. The N=5 <M_x> curves in Fig. 1 are benchmarked against exact diagonalization, an external standard that does not depend on any parameter fitted in this paper. The central physics claims—the susceptibility peak near Gamma=1 for the real-field model and the vanishing <M_x> for the imaginary-field model—are drawn from exact diagonalization at larger system sizes (Secs. 6.1 and 6.2), not from the variational algorithm's output. There are no self-citations by the present authors, no imported uniqueness theorems, and no renamed known result. The paper explicitly acknowledges its ad hoc noise model ('does not capture the complexity of true noise') and the five-qubit resource limit, and it leaves unspecified how the 'ground state' is selected for the non-Hermitian Hamiltonian; but those are correctness and scope limitations, not circular reductions. No equation in the paper is defined in terms of the quantity it is used to predict, and no fitted parameter is relabeled as a prediction. Therefore no circular step is exhibited.
Assumptions & free parameters
free parameters (3)
- Complex eigenvalue search parameter E =
Optimized value (not reported per data point)
- Ansatz circuit parameters theta =
Optimized per field value (values not reported)
- Noise amplitude epsilon =
0.04
assumptions (4)
- domain assumption The variational ansatz G(theta) is expressive enough to represent the ground-state eigenvector, and the optimizer reaches the global minimum of C+.
- ad hoc to paper The uniform noise perturbation (epsilon=0.04) on each expectation value approximates a ~1000-shot quantum measurement.
- standard math Exact diagonalization with the C++ Eigen library provides correct reference eigenpairs.
- domain assumption The finite-size trend of <M_x> for odd spin chains indicates the thermodynamic limit value.
Cite this review
Pith. "Pith review of Quantum Phase Transition of Non-Hermitian Systems using Variational Quantum Techniques." pith.science (2026). https://pith.science/paper/BTM45K7A
@misc{pith2026250117003,
author = {Pith},
title = {Pith review of: Quantum Phase Transition of Non-Hermitian Systems using Variational Quantum Techniques},
year = {2026},
howpublished = {\url{https://pith.science/paper/BTM45K7A}},
note = {Machine review of arXiv:2501.17003}
}
abstract
The motivation for studying non-Hermitian systems and the role of $\mathcal{PT}$-symmetry is discussed. We investigate the use of a quantum algorithm to find the eigenvalues and eigenvectors of non-Hermitian Hamiltonians, with applications to quantum phase transitions. We use a recently proposed variational algorithm. The systems studied are the transverse Ising model with both a purely real and a purely complex transverse field.
Figures
Reference graph
Works this paper leans on
-
[25]
Variational quantum algorithms for scanning the complex spectrum of non-hermitian systems,
X.-D. Xie, Z.-Y. Xue, and D.-B. Zhang, “Variational quantum algorithms for scanning the complex spectrum of non-hermitian systems,”Frontiers of Physics, vol. 19, no. 4, Feb. 2024. [Online]. Available: http://dx.doi.org/10.1007/s11467-023-1382-3
-
[1]
Review on Quantum Computing for Lattice Field Theory,
L. Funcke, T. Hartung, K. Jansen, and S. Kühn, “Review on Quantum Computing for Lattice Field Theory,”PoS, vol. LATTICE2022, p. 228, 2023
work page 2023
-
[2]
Quantum Computing for High-Energy Physics: State of the Art and Challenges,
A. Di Meglioet al., “Quantum Computing for High-Energy Physics: State of the Art and Challenges,” PRX Quantum, vol. 5, no. 3, p. 037001, 2024
work page 2024
-
[3]
Overview of external electromagnetism and rotation in lattice qcd,
A. Yamamoto, “Overview of external electromagnetism and rotation in lattice qcd,” The European Physical Journal A , vol. 57, no. 6, Jun. 2021. [Online]. Available: http://dx.doi.org/10.1140/epja/s10050-021-00530-8
-
[4]
S. Sachdev, “Quantum phase transitions,”Quantum Phase Transitions, 05 2011
work page 2011
-
[5]
Quantum-classical correspondence in quantum channels
B. Vijaywargia and A. Lakshminarayan, “Quantum-classical correspondence in quantum channels,” 2024. [Online]. Available: https://arxiv.org/abs/2407.14067
work page Pith review arXiv 2024
-
[6]
Non-hermitianopticsandphotonics: from classical to quantum,
C.Wang,Z.Fu,W.Mao,J.Qie,A.Stone,andL.Yang,“Non-hermitianopticsandphotonics: from classical to quantum,”Advances in Optics and Photonics, vol. 15, 06 2023
work page 2023
-
[7]
Non-Hermitian Topological Phenomena: A Review,
N. Okuma and M. Sato, “Non-Hermitian Topological Phenomena: A Review,”Ann. Rev. Condensed Matter Phys., vol. 14, pp. 83–107, 2023
work page 2023
Show all 34 references
-
[8]
Coherent control of nonreciprocal optical properties of the defect modes in 1d defective photonic crystals with atomic doping,
N. Ghangas and S. Dasgupta, “Coherent control of nonreciprocal optical properties of the defect modes in 1d defective photonic crystals with atomic doping,” 2024. [Online]. Available: https://arxiv.org/abs/2408.07446
2024 arXiv
-
[9]
Skin effect in non-hermitian systems with spin,
W. Zhang, Y. Hu, H. Zhang, X. Liu, G. Veronis, Y. Shen, Y. Huang, W. Luo, and A. Alu‘, “Skin effect in non-hermitian systems with spin,” 2024. [Online]. Available: https://arxiv.org/abs/2408.07406
2024 arXiv
-
[10]
Moiseyev, Non-Hermitian Quantum Mechanics
N. Moiseyev, Non-Hermitian Quantum Mechanics . Cambridge: Cambridge University Press, 2011. [Online]. Available: https://www.cambridge.org/gb/universitypress/ subjects/physics/quantum-physics-quantum-information-and-quantum-computation/ non-hermitian-quantum-mechanics
2011
-
[11]
A review: Rise ofPT-symmetry for laser applications,
P. Sampath and K. Senthilnathan, “A review: Rise ofPT-symmetry for laser applications,” Optik, vol. 289, p. 171260, 08 2023
2023
-
[12]
Enhanced energy harvesting near exceptionalpointsinsystemswith(pseudo-) PT-symmetry,
L. J. Fernández-Alcázar, R. Kononchuk, and T. Kottos, “Enhanced energy harvesting near exceptionalpointsinsystemswith(pseudo-) PT-symmetry,” Communications Physics,vol.4, no. 1, p. 79, Apr 2021. [Online]. Available: https://doi.org/10.1038/s42005-021-00577-5
2021 doi
-
[13]
Enhancing the sensitivity of frequency and energy splitting detection by using exceptional points: Application to microcavity sensors for single-particle detection,
J. Wiersig, “Enhancing the sensitivity of frequency and energy splitting detection by using exceptional points: Application to microcavity sensors for single-particle detection,” Phys. Rev. Lett. , vol. 112, p. 203901, May 2014. [Online]. Available: https://link.aps.org/doi/10...
2014 doi
-
[14]
Oscillation probabilities for a PT-symmetric non-Hermitian two-state system,
J. Alexandre, M. Dale, J. Ellis, R. Mason, and P. Millington, “Oscillation probabilities for a PT-symmetric non-Hermitian two-state system,” 2 2023
2023
-
[15]
Numerical simulations ofPT-symmetric quantum field theories,
C. W. Bernard and V. M. Savage, “Numerical simulations ofPT-symmetric quantum field theories,” Phys. Rev. D, vol. 64, p. 085010, 2001
2001
-
[16]
Introduction to PT-symmetric quantum theory,
C. M. Bender, “Introduction to PT-symmetric quantum theory,” Contemporary Physics, vol. 46, no. 4, p. 277–292, Jul. 2005. [Online]. Available: http: //dx.doi.org/10.1080/00107500072632
2005 doi
-
[17]
Making sense of non-hermitian hamiltonians,
——, “Making sense of non-hermitian hamiltonians,” Reports on Progress in Physics, vol. 70, no. 6, p. 947–1018, May 2007. [Online]. Available: http: //dx.doi.org/10.1088/0034-4885/70/6/R03
2007 doi
-
[18]
Computational complexity and fundamental limitations to fermionic quantum monte carlo simulations,
M. Troyer and U.-J. Wiese, “Computational complexity and fundamental limitations to fermionic quantum monte carlo simulations,”Physical Review Letters, vol. 94, no. 17, May
-
[19]
Onpropertiesoftheisingmodelforcomplexenergy/temperature and magnetic field,
V.MatveevandR.Shrock,“Onpropertiesoftheisingmodelforcomplexenergy/temperature and magnetic field,”Journal of Physics A: Mathematical and Theoretical , vol. 41, no. 13, p. 135002, 2008
2008
-
[20]
AntiferromagneticIsingmodel in an imaginary magnetic field,
V.Azcoiti,G.DiCarlo,E.Follana,andE.Royo-Amondarain,“AntiferromagneticIsingmodel in an imaginary magnetic field,”Phys. Rev. E, vol. 96, no. 3, p. 032114, 2017
2017
-
[21]
R. A. Horn and C. R. Johnson,Matrix Analysis, 2nd ed. New York, NY, USA: Cambridge University Press, 2013
2013
-
[22]
Sachdev, Quantum Phase Transitions , 2nd ed
S. Sachdev, Quantum Phase Transitions , 2nd ed. Cambridge: Cambridge University Press, 2011. [Online]. Available: https://www.cambridge.org/core/books/ quantum-phase-transitions/33C1C81500346005E54C1DE4223E5562
2011
-
[23]
Quantum phase transitions in non-hermitian PT-symmetric transverse-field ising spin chains,
G. A. Starkov, M. V. Fistul, and I. M. Eremin, “Quantum phase transitions in non-hermitian PT-symmetric transverse-field ising spin chains,”Annals of Physics , vol. 456, p. 169268, Sep. 2023. [Online]. Available: http://dx.doi.org/10.1016/j.aop.2023.169268
2023
-
[24]
TheVariationalQuantumEigensolver: Areviewofmethodsandbestpractices,
J.Tilly et al.,“TheVariationalQuantumEigensolver: Areviewofmethodsandbestpractices,” Phys. Rept., vol. 986, pp. 1–128, 2022
2022
-
[26]
Adam: Amethodforstochasticoptimization,
D.P.Kingma,“Adam: Amethodforstochasticoptimization,” arXiv preprint arXiv:1412.6980, 2014
2014 arXiv
-
[27]
Exact ising model simulation on a quantum computer,
A. Cervera-Lierta, “Exact ising model simulation on a quantum computer,”Quantum, vol. 2, p. 114, Dec. 2018. [Online]. Available: http://dx.doi.org/10.22331/q-2018-12-21-114 9 Quantum Phase Transition of Non-Hermitian Systems using Variational Quantum T echniques James Hancock
2018 doi
-
[28]
A universal variational quantum eigensolver for non- Hermitian systems,
H. Zhao, P. Zhang, and T.-C. Wei, “A universal variational quantum eigensolver for non- Hermitian systems,”Sci. Rep., vol. 13, no. 1, p. 22313, 2023
2023
-
[29]
Variational principles for complex eigenvalues,
C. Detar and F. E. Low, “Variational principles for complex eigenvalues,”Astrophysical Journal, Vol. 227, pp. 349-358 (1979)., vol. 227, pp. 349–358, 1979
1979
-
[30]
Modifiedstochasticvariationalapproachtonon-Hermitianquantum systems,
D.KraftandW.Plessas,“Modifiedstochasticvariationalapproachtonon-Hermitianquantum systems,” J. Phys. Conf. Ser., vol. 738, no. 1, p. 012029, 2016
2016
-
[31]
Guennebaud, B
G. Guennebaud, B. Jacobet al., Eigen: A C++ template library for linear algebra , 2010, version 3.4.0. [Online]. Available: http://eigen.tuxfamily.org
2010
-
[32]
The boundary effects of transverse field ising model,
Y. He and H. Guo, “The boundary effects of transverse field ising model,”Journal of Statistical Mechanics: Theory and Experiment , vol. 2017, no. 9, p. 093101, Sep. 2017. [Online]. Available: http://dx.doi.org/10.1088/1742-5468/aa85b0
2017 doi
-
[33]
Cuaoa: A novel cuda-accelerated simulation framework for the qaoa,
J. Stein, J. Blenninger, D. Bucher, P. J. Eder, E. Çetiner, M. Zorn, and C. Linnhoff-Popien, “Cuaoa: A novel cuda-accelerated simulation framework for the qaoa,” 2024. [Online]. Available: https://arxiv.org/abs/2407.13012 10
2024 arXiv
-
[2005]
Available: http://dx.doi.org/10.1103/PhysRevLett.94.170201
[Online]. Available: http://dx.doi.org/10.1103/PhysRevLett.94.170201
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.