We show that the one-particle density matrix ρ can be used to characterize the interaction-driven many-body localization transition in closed fermionic systems. The natural orbitals (the eigenstates of ρ) are localized in the many-body localized phase and spread out when one enters the delocalized phase, while the occupation spectrum (the set of eigenvalues of ρ) reveals the distinctive Fockspace structure of the many-body eigenstates, exhibiting a step-like discontinuity in the localized phase. The associated one-particle occupation entropy is small in the localized phase and large in the delocalized phase, with diverging fluctuations at the transition. We analyze the inverse participation ratio of the natural orbitals and find that it is independent of system size in the localized phase. Introduction. While the theory of noninteracting disordered systems is well developed [1,2], the possibility of a localization transition in closed interacting systems has only recently been firmly established . This manybody localization (MBL) transition occurs at finite energy densities and is not a conventional thermodynamic transition [24,25]. Instead, it can be understood as a dynamical phase transition, associated with the emergence of a complete set of local conserved quantities in the localized phase, which thus behaves as an integrable system [26][27][28][29][30]. This restricts the entanglement entropy of the eigenstates to an area law [31], in contrast to the volume law predicted by the eigenstate thermalization hypothesis for the ergodic delocalized phase [32][33][34]. At the localization transition, the fluctuations of the entanglement entropy diverge [16,35]. The effects of MBL are also observed in the dynamics following, for example, a global quench from a product state, wherein dephasing between the effective degrees of freedom leads to a characteristic logarithmic growth of the entanglement entropy [6,10,12]. These features comprise a much richer set of signatures than in the context of noninteracting systems, for which, in the spirit of one-parameter scaling, the notion of a localization length based on single-particle wave functions generally suffices [1,2].