Many-body characterization of particle-conserving topological superfluids

Determination of electronic excitation energies within the doubly occupied configuration interaction space by means of the Hermitian operator method

In this work, we formulate the equations of motion corresponding to the Hermitian operator method in the framework of the doubly occupied configuration interaction space. The resulting algorithms turn out to be considerably simpler than the equations provided by that method in more conventional spaces, enabling the determination of excitation energies in N-electron systems under an affordable polynomial computational cost. The implementation of this technique only requires to know the elements of low-order reduced density matrices of an N-electron reference state, which can be obtained from any approximate method. We contrast our procedure against the reduced Bardeen-Cooper-Schrieffer and Richardson-Gaudin-Kitaev integrable models, pointing out the reliability of our proposal.

Solving the nuclear pairing model with neural network quantum states

We present a variational Monte Carlo method that solves the nuclear many-body problem in the occupation number formalism exploiting an artificial neural network representation of the ground-state wave function. A memory-efficient version of the stochastic reconfiguration algorithm is developed to train the network by minimizing the expectation value of the Hamiltonian. We benchmark this approach against widely used nuclear many-body methods by solving a model used to describe pairing in nuclei for different types of interaction and different values of the interaction strength. Despite its polynomial computational cost, our method outperforms coupled-cluster and provides energies that are in excellent agreement with the numerically exact full configuration-interaction values.

Variational determination of ground and excited-state two-electron reduced density matrices in the doubly occupied configuration space: A dispersion operator approach

This work implements a variational determination of the elements of two-electron reduced density matrices corresponding to the ground and excited states of N-electron interacting systems based on the dispersion operator technique. The procedure extends the previously reported proposal [Nakata et al., J. Chem. Phys. 125, 244109 (2006)] to two-particle interaction Hamiltonians and N-representability conditions for the two-, three-, and four-particle reduced density matrices in the doubly occupied configuration interaction space. The treatment has been applied to describe electronic spectra using two benchmark exactly solvable pairing models: reduced Bardeen-Cooper-Schrieffer and Richardson-Gaudin-Kitaev Hamiltonians. The dispersion operator combined with N-representability conditions up to the four-particle reduced density matrices provides excellent results.

Floquet Gauge Pumps as Sensors for Spectral Degeneracies Protected by Symmetry or Topology

We introduce the concept of a Floquet gauge pump whereby a dynamically engineered Floquet Hamiltonian is employed to reveal the inherent degeneracy of the ground state in interacting systems. We demonstrate this concept in a one-dimensional XY model with periodically driven couplings and transverse field. In the high-frequency limit, we obtain the Floquet Hamiltonian consisting of the static XY and dynamically generated Dzyaloshinsky-Moriya interaction (DMI) terms. The dynamically generated magnetization current depends on the phases of complex coupling terms, with the XY interaction as the real and DMI as the imaginary part. As these phases are cycled, the current reveals the ground-state degeneracies that distinguish the ordered and disordered phases. We discuss experimental requirements needed to realize the Floquet gauge pump in a synthetic quantum spin system of interacting trapped ions.

Two-Fluid Coexistence in a Spinless Fermions Chain with Pair Hopping

We show that a simple one-dimensional model of spinless fermions with pair hopping displays a phase in which a Luttinger liquid of paired fermions coexists with a Luttinger liquid of unpaired fermions. Our results are based on extensive numerical density-matrix renormalization-group calculations and are supported by a two-fluid model that captures the essence of the coexistence region.

Rigorous Results on Topological Superconductivity with Particle Number Conservation

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.

Certified Quantum Measurement of Majorana Fermions

We present a quantum self-testing protocol to certify measurements of fermion parity involving Majorana fermion modes. We show that observing a set of ideal measurement statistics implies anti-commutativity of the implemented Majorana fermion parity operators, a necessary prerequisite for Majorana detection. Our protocol is robust to experimental errors. We obtain lower bounds on the fidelities of the state and measurement operators that are linear in the errors. We propose to analyze experimental outcomes in terms of a contextuality witness W, which satisfies 〈W〉 ≤ 3 for any classical probabilistic model of the data. A violation of the inequality witnesses quantum contextuality, and the closeness to the maximum ideal value 〈W〉 = 5 indicates the degree of confidence in the detection of Majorana fermions.

Benchmarking the Variational Reduced Density Matrix Theory in the Doubly Occupied Configuration Interaction Space with Integrable Pairing Models

The variational reduced density matrix theory has been recently applied with great success to models within the truncated doubly occupied configuration interaction space, which corresponds to the seniority zero subspace. Conservation of the seniority quantum number restricts the Hamiltonians to be based on the SU(2) algebra. Among them there is a whole family of exactly solvable Richardson-Gaudin pairing Hamiltonians. We benchmark the variational theory against two different exactly solvable models, the Richardson-Gaudin-Kitaev and the reduced BCS Hamiltonians. We obtain exact numerical results for the so-called [Formula: see text] N-representability conditions in both cases for systems that go from 10 to 100 particles. However, when random single-particle energies as appropriate for small superconducting grains are considered, the exactness is lost but still a high accuracy is obtained.

Crystalline Symmetry-Protected Majorana Mode in Number-Conserving Dirac Semimetal Nanowires

One of the cornerstones for topological quantum computations is the Majorana zero mode, which has been intensively searched in fractional quantum Hall systems and topological superconductors. Several recent works suggest that such an exotic mode can also exist in a one-dimensional (1D) interacting double-wire setup even without long-range superconductivity. A notable instability in these proposals comes from interchannel single-particle tunneling that spoils the topological ground state degeneracy. Here we show that a 1D Dirac semimetal (DSM) nanowire is an ideal number-conserving platform to realize such Majorana physics. By inserting magnetic flux, a DSM nanowire is driven into a 1D crystalline-symmetry-protected semimetallic phase. Interaction enables the emergence of boundary Majorana zero modes, which is robust as a result of crystalline symmetry protection. We also explore several experimental consequences of Majorana signals.

Polynomial description of inhomogeneous topological superconducting wires

We present the universal features of the topological invariant for p-wave superconducting wires after the inclusion of spatial inhomogeneities. Three classes of distributed potentials are studied, a single-defect, a commensurate and an incommensurate model, using periodic site modulations. An analytic polynomial description is achieved by splitting the topological invariant into two parts; one part depends on the chemical potential and the other does not. For the homogeneous case, an elliptical region is found where the topological invariant oscillates. The zeros of these oscillations occur at points where the fermion parity switches for finite wires. The increase of these oscillations with the inhomogeneity strength leads to new isolated non-topological phases. We characterize these new phases according to each class of spatial distributions. Such phases could also be observed in the XY model, to which our model is dual.

Many-Body Topological Invariants for Fermionic Symmetry-Protected Topological Phases

We define and compute many-body topological invariants of interacting fermionic symmetry-protected topological phases, protected by an orientation-reversing symmetry, such as time-reversal or reflection symmetry. The topological invariants are given by partition functions obtained by a path integral on unoriented spacetime which, as we show, can be computed for a given ground state wave function by considering a nonlocal operation, "partial" reflection or transpose. As an application of our scheme, we study the Z_{8} and Z_{16} classification of topological superconductors in one and three dimensions.

Majorana Quasiparticles Protected by Z_{2} Angular Momentum Conservation

We show how angular momentum conservation can stabilize a symmetry-protected quasitopological phase of matter supporting Majorana quasiparticles as edge modes in one-dimensional cold atom gases. We investigate a number-conserving four-species Hubbard model in the presence of spin-orbit coupling. The latter reduces the global spin symmetry to an angular momentum parity symmetry, which provides an extremely robust protection mechanism that does not rely on any coupling to additional reservoirs. The emergence of Majorana edge modes is elucidated using field theory techniques, and corroborated by density-matrix-renormalization-group simulations. Our results pave the way toward the observation of Majorana edge modes with alkaline-earth-like fermions in optical lattices, where all basic ingredients for our recipe-spin-orbit coupling and strong interorbital interactions-have been experimentally realized over the last two years.

Exact Solution of Quadratic Fermionic Hamiltonians for Arbitrary Boundary Conditions

We present a procedure for exactly diagonalizing finite-range quadratic fermionic Hamiltonians with arbitrary boundary conditions in one of D dimensions, and periodic in the remaining D-1. The key is a Hamiltonian-dependent separation of the bulk from the boundary. By combining information from the two, we identify a matrix function that fully characterizes the solutions, and may be used to construct an efficiently computable indicator of bulk-boundary correspondence. As an illustration, we show how our approach correctly describes the zero-energy Majorana modes of a time-reversal-invariant s-wave two-band superconductor in a Josephson ring configuration, and predicts that a fractional 4π-periodic Josephson effect can only be observed in phases hosting an odd number of Majorana pairs per boundary.

Localized Majorana-Like Modes in a Number-Conserving Setting: An Exactly Solvable Model

In this Letter we present, in a number conserving framework, a model of interacting fermions in a two-wire geometry supporting nonlocal zero-energy Majorana-like edge excitations. The model has an exactly solvable line, on varying the density of fermions, described by a topologically nontrivial ground state wave function. Away from the exactly solvable line we study the system by means of the numerical density matrix renormalization group. We characterize its topological properties through the explicit calculation of a degenerate entanglement spectrum and of the braiding operators which are exponentially localized at the edges. Furthermore, we establish the presence of a gap in its single particle spectrum while the Hamiltonian is gapless, and compute the correlations between the edge modes as well as the superfluid correlations. The topological phase covers a sizable portion of the phase diagram, the solvable line being one of its boundaries.