η-pairing in correlated fermion models with spin-orbit coupling 1,2, Kai Li ∗ 1School of Physics, Zhengzhou University, Zhengzhou 450001, China 2National Laboratory of Solid State Microstructures and School of Physics, Nanjing University, Nanjing 210093, China (Dated: October 30, 2020) We generalize the η-pairing theory in Hubbard models to the ones with spin-orbit coupling (SOC) and obtain the conditions under which the η-pairing operator is an eigenoperator of the Hamiltonian. The η pairing thus reveals an exact SU(2) pseudospin symmetry in our spin-orbit coupled Hubbard model, even though the SU(2) spin symmetry is explicitly broken by the SOC. In particular, these exact results can be applied to a variety of Hubbard models with SOC on either bipartite or non- bipartite lattices, whose noninteracting limit can be a Dirac semimetal, a Weyl semimetal, a nodal- line semimetal, and a Chern insulator. The η pairing conditions also impose constraints on the band topology of these systems. We then construct and focus on an interacting Dirac-semimetal model, which exhibits an exact pseudospin symmetry with fine-tuned parameters. The stability regions for the exact η-pairing ground states (with momentum π or 0) and the exact charge-density-wave ground states are established. Between these distinct symmetry-breaking phases, there exists an exactly solvable multicritical line. In the end, we discuss possible experimental realizations of our results. I. INTRODUCTION in these systems, while the SU(2) spin symmetry is ex- plicitly broken by the SOC. We then exemplify these Exact results or solutions of interacting many-electron exact results with four concrete spin-orbit coupled sys- models beyond one spatial dimension are rare and highly tems, including a Dirac semimetal, a Weyl semimetal, appreciated. One important advance in studying the a nodal-line semimetal, and a Chern insulator, in the Hubbard model is the so-called η-pairing theory, which noninteracting limit. Some consequences of the pseu- is an exact result due to Yang1. The η-pairing operator dospin symmetry and the relation between the present allows for the construction of exact eigenstates possess- study and recent works are discussed. To convert the ing off-diagonal long-range order (ODLRO) and reveals η-pairing eigenstates into ground states, additional in- an SU(2) pseudospin symmetry2 of the Hubbard model. teractions besides the Hubbard U are needed. In par- More importantly, the η-pairing states can serve as exact ticular, we construct and focus on an interacting Dirac- superconducting ground states in a class of strongly cor- semimetal model on the square lattice. By making the related electronic models3–6, which is beyond the frame- Hamiltonian frustration-free in some parameter regions, work of BCS theory. Studies inspired by the η-pairing we establish the exact η-pairing ground states and the theory have since been pursued, including the unifica- exact charge-density-wave (CDW) ground states, where 0 tion of various symmetry-breaking orders (and related the η-pairing ground states with momentum π and are multicritical behavior study)2,7–11, the construction of an identified with a pair-density-wave (PDW) phase and a SU(2) gauge theory of the Hubbard model12–14, the cal- uniform superconductor, respectively. Moreover, there culation of entanglement entropy of η-pairing states15,16, exists an exactly solvable multicritical line between the and the study of quantum thermalization physics and three symmetry-breaking phases. In the end, we discuss quantum many-body scars in Hubbard-like models17–22. possible experimental realizations of our results in optical lattices. It is now recognized that spin-orbit coupling (SOC) can produce a wealth of topological phases and phys- ical systems with SOC are currently under intense in- vestigation. In particular, a number of theoretical un- arXiv:1901.06914v4 [cond-mat.str-el] 28 Oct 2020 derstandings have been obtained for strongly correlated electronic systems with SOC23–25. Nevertheless, most of these theoretical investigations used various approximate methods such as mean field theory and more rigorous and This paper is organized as follows: In Sec. II we in- clearer physical understanding of these systems deserves troduce our model and generalize the η-pairing theory. further studies. Therefore, exact results or solutions for In Sec. III we exemplify the generalized η-pairing the- spin-orbit coupled interacting electronic systems would ory with four concrete models. An interacting Dirac- be valuable. semimetal model with exact η-pairing ground states is In this work, we first show that Yang’s η-pairing the- presented in Sec. IV and we identify an exactly solvable ory can be generalized to a class of Hubbard models with multicritical line in Sec. V. We discuss possible exper- SOC and derive the equations which make the η-pairing imental realizations in Sec. VI and present a summary operator an eigenoperator of the Hamiltonian. The ex- and discussion in Sec. VII. Certain details can be found act SU(2) pseudospin symmetry can be then preserved in the appendices. 2 II. MODEL AND GENERALIZATION OF η (note the different meaning of an eigenoperator intro- PAIRING duced in the context of the rigged Hilbert space26,27) of the total Hamiltonian, namely The generic Hamiltonian for our spin-orbit coupled Hubbard models in arbitrary spatial dimensions is given [H, ηQ† ] = (λ + U)ηQ† . (6) by Therefore, one can construct many exact eigenstates H = H0 + HU m m H(ηQ† ) |0i = m(λ + U)(ηQ† ) |0i possessing ODLRO. H0 = Ck† (αkσx + βkσy + γkσz + εkσ0)Ck The SU(2) spin symmetry is explicitly broken by the Xk SOC terms in Eq.(1). However, from Eq.(6), we see that our spin-orbit coupled Hubbard model possesses an exact HU = U ni ni . (1) ↑ ↓ SU(2) pseudospin symmetry when λ + U = 0. The term Xi pseudospin is attributed to the SU(2) algebra formed by The noninteracting part H0 is written in momentum η-pairing operators, i.e., [Jα,Jβ] = iǫαβγJγ (α,β,γ = space, where C† = (c† ,c† ). In H , the three 1 1 k k k 0 x,y,z), where J = (η† + ηQ),J = (η† − ηQ), J = ↑ ↓ x 2 Q y 2i Q z momentum-dependent real coefficients αk,βk,γk with 1 (Nˆe −N), Nˆe = niσ, and N is the number of lattice Pauli matrices σx,y,z represent the SOC, and the term εk 2 iσ sites. P with the identity matrix σ0 is the usual energy band with- out SOC. The interacting part H written in real space Physically, the pseudospin symmetry can be under- U stood via a particle-hole transformation for the spin- is the on-site Hubbard interaction, where niσ = ciσ† ciσ. iQ j down electrons: fj† = cj† ,fj† = e− · cj . The pseu- In the absence of SOC (αk = βk = γk = 0), Yang finds ↑ ↑ ↓ ↓ iπ j dospin operators can now be rewritten in terms of the f that the so-called η-pairing operator ηπ† = j e · cj† cj† 1 ↑ ↓ fermions as Jα = 2 j fj†ταfj, where τα are the Pauli obeys the commutation relations [H0, ηπ† ]P = −2µηπ† and matrices and fj† = (fPj† ,fj† ). Thus, the pseudospin can [HU , ηπ† ] = Uηπ† , where µ is the chemical potential and ↑ ↓ π = (π,π,π) or (π, π) for the cubic or square lat- be interpreted as the spin of f fermions. This is rem- iniscent of the SU(2) gauge structure in the fermionic- tice. The total Hamiltonian and η† thus obey [H, η† ]= π π spinon representation of spin operators28–30. It is impor- (U − 2µ)ηπ† , and one can construct many exact eigen- m m tant to emphasize that the η-pairing states do not have states H(ηπ† ) |0i = m(U − 2µ)(ηπ† ) |0i (m = 0, 1, ...) possessing ODLRO. the pseudospin symmetry. In fact, they correspond to In Yang’s observation, the key point that makes the the ferromagnetic states of f fermions with maximal total 2 m N N m ”spin” N/2, i.e., J (η† ) |0i = ( + 1)(η† ) |0i and possible η pairing is the fact that the energy band εk Q 2 2 Q m N m 2 2 (for the cubic or square lattice without SOC) has the Jz(ηQ† ) |0i = (m − 2 )(ηQ† ) |0i (where J = α Jα), property εk + επ k = −2µ (independent of k), which which break the f fermion’s spin symmetry and hence − P leads to the commutation relation [H0, ηπ† ]= −2µηπ† . break the pseudospin symmetry. Therefore, to generalize the η-pairing theory to Hub- bard models with SOC, e.g., Eq.(1), we need some new constraints on the coefficients αk,βk, and γk. To achieve III. APPLICATION OF THE GENERALIZED these new constraints, we first define the η-pairing oper- η-PAIRING THEORY ator with a generic momentum Q as Based on the above exact results for model (1), we find iQ j ηQ† = ck† cQ† k = e · cj† cj† . (2) that the the η pairing can exist in a variety of spin-orbit ↑ − ↓ ↑ ↓ Xk Xj coupled systems. In particular, we apply these exact re- sults to four concrete examples, which include three gap- It can then be shown that (see Appendix A) less and one gapped noninteracting topological phases due to SOC. [H0, η† ]= λη† , (3) Q Q The first example is a Dirac-semimetal Hubbard model if αk,βk,γk,εk, and Q satisfy the following equations on a square lattice. Its Hamiltonian is given by Eq.(1) with αk = 2t′ sin kx,βk = 2t′ sin ky,γk = 0, and αk = αQ k, βk = βQ k, − − εk = −2t(cos kx + cos ky) − µ. The lattice spacing εk + εQ k + γk − γQ k = λ, (4) has been set to unity throughout our paper.