Nodal Andreev spectra in multi-Majorana three-terminal Josephson junctions

We investigate the Andreev-bound-state (ABS) spectra of three-terminal Josephson junctions which consist of 1D topological superconductors (TSCs) harboring multiple zero-energy edge Majorana bound states (MBSs) protected by chiral symmetry. Our theoretical analysis relies on the exact numerical diagonalization of the Bogoliubov-de Gennes (BdG) Hamiltonian describing the three interfaced TSCs, complemented by an effective low-energy description solely based on the coupling of the interfacial MBSs arising before the leads get contacted. Considering the 2D synthetic space spanned by the two independent superconducting phase differences, we demonstrate that the ABS spectra may contain either point or line nodes and identify Z_2 topological invariants to classify them. We show that the resulting type of nodes depends on the number of preexisting interfacial MBSs, with nodal lines necessarily appearing when two TSCs harbor an unequal number of MBSs. Specifically, the precise number of interfacial MBSs determines the periodicity of the spectrum under 2π slidings of the phase differences and, as a result, also controls the shape of the nodal lines in synthetic space. When chiral symmetry is preserved, the lines are open and coincide with high-symmetry lines of synthetic space, while when it is violated, the lines can also transform into loops and chains. The nodal spectra are robust by virtue of the inherent particle-hole symmetry of the BdG Hamiltonian, and give rise to distinctive experimental signatures that we identify.