As one of the essential features in non-Hermitian systems coupled with environment, the exceptional point has attracted much attention in many physical fields. The phenomena that eigenvalues and eigenvectors of the system simultaneously coalesce at the exceptional point are also one of the important properties to distinguish from Hermitian systems. In non-Hermitian systems with parity-time reversal symmetry, the eigenvalues can be continuously adjusted in parameter space from all real spectra to pairs of complex-conjugate values by crossing the phase transition from the parity-time reversal symmetry preserving phase to the broken phase. The phase transition point is called an exceptional point of the system, which occurs in company with the spontaneous symmetry broken and many novel physical phenomena, such as sensitivity-enhanced measurement and loss induced transparency or lasing. Here, we focus on a two-qubit quantum system with parity-time reversal symmetry and construct an experimental scheme, prove and verify the features at its third-order exceptional point, including high-order energy response induced by perturbation and the coalescence of eigenvectors. We first theoretically study a two-qubit non-Hermitian system with parity-time reversal symmetry, calculate the properties of eigenvalues and eigenvectors, and prove the existence of a third-order exceptional point. Then, in order to study the energy response of the system induced by perturbation, we introduce an Ising-type interaction as perturbation and quantitatively demonstrate the response of eigenvalues. In logarithmic coordinates, three of the eigenvalues are indeed in the cubic root relationship with perturbation strength, while the fourth one is a linear function. Moreover, we study the eigenvectors around exceptional point and show the coalescence phenomenon as the perturbation strength becomes smaller. The characterization of the response of eigenvalues at high-order exceptional points is a quite difficult task as it is in general difficult to directly measure eigenenergies in a quantum system composed of a few qubits. In practice, the time evolution of occupation on a particular state is used to indirectly fit the eigenvalues. In order to make the fitting of experimental data more reliable, we want to determine an accurate enough expressions for the eigenvalues and eigenstates. To this aim, we employ a perturbation treatment and show good agreement with the numerical results of states occupation obtained by direct evolution. Moreover, we find that after the system evolves for a long enough time, it will end up to one of the eigenstates, which gives us a way to demonstrate eigenvector coalescence by measuring the density matrix via tomography and parity-time reversal transformation. To show our scheme is experimentally applicable, we propose an implementation using trapped 171Yb+ ions. We can map the parity-time reversal symmetric Hamiltonian to a purely dissipative two-ion system: use microwave to achieve spin state inversion, shine a 370 nm laser to realize dissipation of spin-up state, and apply Raman operation for Mølmer-Sørensen gates to implement Ising interaction. By adjusting the corresponding microwave and laser intensity, the spin coupling strength, the dissipation rate and the perturbation strength can be well controlled. We can record the probability distribution of the four product states of the two ions and measure the density matrix by detecting the fluorescence of each ion on different Pauli basis. © 2022 Chinese Physical Society
