Quantum quench in a p plus ip superfluid: Winding numbers and topological states far from equilibrium

We study the nonadiabatic dynamics of a two-dimensional p + ip superfluid following an instantaneous quantum quench of the BCS coupling constant. The model describes a topological superconductor with a nontrivial BCS (trivial BEC) phase appearing at weak-(strong-) coupling strengths. We extract the exact long-time asymptotics of the order parameter Delta(t) by exploiting the integrability of the classical p-wave Hamiltonian, which we establish via a Lax construction. Three different types of asymptotic behavior can occur depending upon the strength and direction of the interaction quench. We refer to these as the nonequilibrium phases {I, II, III}, characterized as follows. In phase I, the order parameter asymptotes to zero due to dephasing. In phase II, Delta -> Delta(infinity), a nonzero constant. Phase III is characterized by persistent oscillations of Delta(t). For quenches within phases I and II, we determine the topological character of the asymptotic states. We show that two different formulations of the bulk topological winding number, although equivalent in the BCS or BEC ground states, must be regarded as independent out of equilibrium. The first winding number Q characterizes the Anderson pseudospin texture of the initial state; we show that Q is generically conserved. For Q not equal 0, this leads to the prediction of a gapless topological state when Delta asymptotes to zero. The presence or absence of Majorana edge modes in a sample with a boundary is encoded in the second winding number W, which is formulated in terms of the retarded Green's function. We establish that W can change following a quench across the quantum critical point. When the order parameter asymptotes to a nonzero constant, the final value of W is well defined and quantized. We discuss the implications for the (dis)appearance of Majorana edge modes. Finally, we show that the parity of zeros in the bulk out-of-equilibrium Cooper-pair distribution function constitutes a Z(2)-valued quantum number, which is nonzero whenever W not equal Q. The pair distribution can in principle be measured using rf spectroscopy in an ultracold-atom realization, allowing direct experimental detection of the Z(2) number. This has the following interesting implication: topological information that is experimentally inaccessible in the bulk ground state can be transferred to an observable distribution function when the system is driven far from equilibrium.