The quantum critical properties of interacting fermions in the presence of disorder are still not fully understood. While it is well known that for Dirac fermions, interactions are irrelevant to the noninteracting infinite randomness fixed point (IRFP), the problem remains largely open in the case of Majorana fermions which further display a much richer disorder-free phase diagram. Here, pushing the limits of density matrix renormalization group simulations, we carefully examine the ground state of a Majorana chain with both disorder and interactions. Building on appropriate boundary conditions and key observables such as entanglement, energy gap, and correlations, we strikingly find that the noninteracting Majorana IRFP is very stable against finite interactions, in contrast with previous claims.

We study Majorana chain with the shortest possible interaction term and in the presence of hopping alternation. When formulated in terms of spins the model corresponds to the transverse field Ising model with nearest-neighbor transverse and next-nearest-neighbor longitudinal repulsion. The phase diagram obtained with extensive DMRG simulations is very rich and contains six phases. Four gapped phases include paramagnetic, period-2 with broken translation symmetry, Z₂ with broken parity symmetry and the period-2-Z₂ phase with both symmetries broken. In addition there are two floating phases: gapless and critical Luttinger liquid with incommensurate correlations, and with an additional spontaneously broken Z₂ symmetry in one of them. By analyzing an extended phase diagram we demonstrate that, in contrast with a common belief, the Luttinger liquid phase along the self-dual critical line terminates at a weaker interaction strength than the end point of the Ising critical line that we find to be in the tri-critical Ising universality class. We also show that none of these two points is a Lifshitz point terminating the incommensurability. In addition, we analyzed topological properties through Majorana zero modes emergent in the two topological phases, with and without incommensurability. In the weak interaction regime, a self-consistent mean-field treatment provides a remarkable accuracy for the description of the spectral pairing and the parity switches induced by the interaction.