What distinguishes trivial superfluids from topological superfluids in interacting many-body systems where the number of particles is conserved? Building on a class of integrable pairing Hamiltonians, we present a number-conserving, interacting variation of the Kitaev model, the Richardson-Gaudin-Kitaev chain, that remains exactly solvable for periodic and antiperiodic boundary conditions. Our model allows us to identify fermion parity switches that distinctively characterize topological superconductivity (fermion superfluidity) in generic interacting many-body systems. Although the Majorana zero modes in this model have only a power-law confinement, we may still define many-body Majorana operators by tuning the flux to a fermion parity switch. We derive a closed-form expression for an interacting topological invariant and show that the transition away from the topological phase is of third order.
What distinguishes trivial superfluids from topological superfluids in interacting many-body systems where the number of particles is conserved? Building on a class of integrable pairing Hamiltonians, we present a number-conserving, interacting variation of the Kitaev model, the Richardson-Gaudin-Kitaev chain, that remains exactly solvable for periodic and antiperiodic boundary conditions. Our model allows us to identify fermion parity switches that distinctively characterize topological superconductivity (fermion superfluidity) in generic interacting many-body systems. Although the Majorana zero modes in this model have only a power-law confinement, we may still define many-body Majorana operators by tuning the flux to a fermion parity switch. We derive a closed-form expression for an interacting topological invariant and show that the transition away from the topological phase is of third order.
We investigate the spectral and topological properties of a zigzag Kitaev chain constructed from two diagonally coupled one dimensional Kitaev chains with a zero and finite superconducting pairing phase difference. Using a Bogoliubov-de Gennes formulation, we analyze the energy spectrum, distribution of Majorana zero modes (MZMs), the quasi-particle dispersion, and the winding number, respectively. For a zero phase difference, the resulting energy spectrum shows topological phases with two, four MZMs, and trivial regions. The phases of gap closure determine the topological phase boundaries. In particular, for the phase difference between φ=π, the degeneracy of MZMs is partially lifted, leading to modified topological phases compared to the case φ=0. The topological and trivial phase boundaries are further confirmed by evaluating the quasi-particle dispersion and the topological invariant, namely the winding number. We show that the zigzag Kitaev chain contributes independently to the total winding number ν=1 and 2, giving rise to distinct topological phases that support two and four MZMs. The ν=0 gives rise to a trivial region. The energy spectrum of systems corroborates the analytical phase boundaries and reveals characteristics associated with hybridization, enabling us to obtain the complete phase diagram of the zigzag model. Our results establish the zigzag Kitaev chain as a minimal platform for engineering MZM quantum computations, with potential applications in the study of topological phases and Majorana based qubit physics.
Topological superconductors are believed to host exotic quasiparticle excitations known as Majorana zero modes (MZMs), with much of the evidence based on BCS mean-field theory. The direct application of mean-field arguments is tenuous in finite, isolated systems relevant in some experiments. Here, we develop a new correlation-based method for identifying MZMs in interacting, number-conserving systems. Using the density matrix renormalization group, we study fermion number-conserving models with long-range interactions, which under periodic boundary conditions exhibit robust topological and nontopological superconductivity, tuned by the strength of interaction [Ortiz et al., Phys. Rev. Lett. 113, 267002 (2014)PRLTAO0031-900710.1103/PhysRevLett.113.267002]. We find evidence that, on the topological side, Majorana edge modes appear in open chains, manifesting as the vanishing of the energy splitting between odd- and even-parity ground states with increasing system size. Additionally, off-diagonal two-point correlation functions show nonlocal, parity-dependent edge effects. These correlations reveal the spatial structure of Majorana modes in the many-body wave function. We show that the correlation diagnostic applies broadly, including to short-range interacting models, where topological superconductivity is more fragile due to the absence of a bulk excitation gap.
Most theoretical studies of topological superconductors and Majorana-based quantum computation rely on a mean-field approach to describe superconductivity. A potential problem with this approach is that real superconductors are described by number-conserving Hamiltonians with long-range interactions, so their topological properties may not be correctly captured by mean-field models that violate number conservation and have short-range interactions. To resolve this issue, reliable results on number-conserving models of superconductivity are essential. As a first step in this direction, we use rigorous methods to study a number-conserving toy model of a topological superconducting wire. We prove that this model exhibits many of the desired properties of the mean-field models, including a finite energy gap in a sector of fixed total particle number, the existence of long-range Majorana-like correlations between the ends of an open wire, and a change in the ground state fermion parity for periodic vs antiperiodic boundary conditions. These results show that many of the remarkable properties of mean-field models of topological superconductivity persist in more realistic models with number-conserving dynamics.
This thesis is devoted to the topic of one- and two-dimensional models of topological superconductivity. We study the Kitaev chain subject to closed and open boundary conditions. In the closed chain, we derive the energy spectrum, the ground state, and the topological invariant viewed as a certain Berry phase. We study two-point correlation functions and find enlarged values close to the topological phase transitions. The open chain is studied with focus on describing the degenerate ground states for a simple parameter choice in the topological phase. Then, the open system Hamiltonian is diagonalized numerically, and we model the order parameter with a spatial dependency. Additional Majorana zero modes appear if the order parameter changes sign. We also consider the p+ip model. The localization of a Majorana zero mode, bound to a vortex that is described by Ginzburg-Landau theory, is found numerically. We propose an argument that results in a non-Abelian exchange transformation for a system of several vortices. The p+ip model is also studied on an annulus; we approximate the ground state as a combination of boundary states and calculate its energy. The results are compared to a numerical implementation, and the agreement is convincing as the boundary separation becomes large compared to other length scales.
A time-reversal invariant Kitaev-type model is introduced in which spins (Dirac matrices) on the square lattice interact via anisotropic nearest-neighbor and next-nearest-neighbor exchange interactions. The model is exactly solved by mapping it onto a tight-binding model of free Majorana fermions coupled with static Z2 gauge fields. The Majorana fermion model can be viewed as a model of time-reversal-invariant superconductor and is classified as a member of symmetry class DIII in the Altland-Zirnbauer classification. The ground-state phase diagram has two topologically distinct gapped phases which are distinguished by a Z2 topological invariant. The topologically nontrivial phase supports both a Kramers' pair of gapless Majorana edge modes at the boundary and a Kramers' pair of zero-energy Majorana states bound to a 0-flux vortex in the $\ensuremath{\pi}$-flux background. Power-law decaying correlation functions of spins along the edge are obtained by taking the gapless Majorana edge modes into account. The model is also defined on the one-dimension ladder, in which case again the ground-state phase diagram has Z2 trivial and nontrivial phases.
We investigate the robustness of Majorana edge modes under disorder and interactions. We exploit a recently found mapping of the interacting Kitaev chain in the symmetric region ($\ensuremath{\mu}=0, t=\mathrm{\ensuremath{\Delta}}$) to free fermions. Extending the exact solution to the disordered case allows us to calculate analytically the topological phase boundary for all interaction and disorder strengths, which has been thought to be only accessible numerically. We discover a regime in which moderate disorder in the interaction matrix elements enhances topological order well into the strongly interacting regime U>t. We also derive the explicit form of the many-body Majorana edge wave function, revealing how it is dressed by many-particle fluctuations from interactions. The qualitative features of our analytical results are valid beyond the fine-tuned integrable point, as expected from the robustness of topological order and as corroborated here by an exact diagonalization study of small systems.