In recent years, tools from quantum information theory have become indispensable in characterizing many-body systems. In this work, we employ measures of entanglement to study the interplay between disorder and the topological phase in 1D systems of the Kitaev type, which can host Majorana end modes at their edges. We find that the entanglement entropy may actually increase as a result of disorder, and identify the origin of this behavior in the appearance of an infinite-disorder critical point. We also employ the entanglement spectrum to accurately determine the phase diagram of the system, and find that disorder may enhance the topological phase, and lead to the appearance of Majorana zero modes in systems whose clean version is trivial.Entanglement is the quintessential characteristic of the quantum world. Recently, much attention has been devoted to manifestations of entanglement in quantum many-body systems [1,2]. Motivations include the prospect of using many-body ground and excited states as a resource for quantum communication and computation; the amount of entanglement controls the applicability of a powerful tensor-network-based numerical method; entanglement is an indicator of quantum correlations which is independent of the details of the system, allowing, in particular, one to identify quantum phase transitions even without knowing the relevant order parameter. In this work, we will concentrate on the latter property, in the context of disordered 1D quantum many-body systems.Majorana zero modes [3][4][5][6] and their more complicated relatives [7] have been studied extensively recently, both due to their predicted exotic non-abelian braiding statistics, which gives rise to the prospect of topologically-protected quantum computation [8], and to many concrete proposals for their realization in the laboratory. Leading candidates are semiconductor quantum wires with strong spin orbit interaction, which are rendered effectively spinless by the application of an appropriate magnetic field and gate voltage, and driven into a topological phase by proximity-coupling to a superconductor, leading to the formation of Majorana end modes [9,10]. Indications for these modes have recently been measured experimentally [11][12][13][14][15]. These Majorana zero modes show up not only in the presence of a real edge, but also in the entanglement spectrum, that is, in the spectrum of the reduced density matrix of a subsystem [16].Disorder naturally occurs in all these systems, and may hamper their topological characteristics, especially the Majorana edge modes of 1D systems . In this work, we will study the interplay of disorder and the topological phase from the entanglement point of view. After introducing the model and our method for calculating the entanglement spectrum and entropy in Section 2, we will