In this work, we present a result on the non-equilibrium dynamics causing equilibration and Gaussification of quadratic non-interacting fermionic Hamiltonians. Specifically, based on two basic assumptions -clustering of correlations in the initial state and the Hamiltonian exhibiting delocalizing transport -we prove that nonGaussian initial states become locally indistinguishable from fermionic Gaussian states after a short and well controlled time. This relaxation dynamics is governed by a power-law independent of the system size. Our argument is general enough to allow for pure and mixed initial states, including thermal and ground states of interacting Hamiltonians on and large classes of lattices as well as certain spin systems. The argument gives rise to rigorously proven instances of a convergence to a generalized Gibbs ensemble. Our results allow to develop an intuition of equilibration that is expected to be more generally valid and relates to current experiments of cold atoms in optical lattices.Despite the great complexity of quantum many-body systems out-of-equilibrium, local expectation values in such systems show the remarkable tendency to equilibrate to stationary values that do not depend on the microscopic details of the initial state, but rather can be described with few parameters using thermal states or generalized Gibbs ensembles [1][2][3]. Such behavior has been successfully studied in many settings theoretically and experimentally, most notably in instances of quantum simulations in optical lattices [2,4,5].By now, it is clear that, despite the unitary nature of quantum mechanical evolution, local expectation values equilibrate due to a dephasing between the eigenstates [3,[6][7][8][9][10][11]. So far it is, however, unclear why this dephasing tends to happen so rapidly. In fact, experiments often observe equilibration after very short times which are independent of the system size [5,12], while even the best theoretical bounds for general initial states of concrete systems diverge exponentially [2,11]. This discrepancy poses the challenge of precisely identifying the equilibration time, which constitutes one of the main open questions in the field [1][2][3].What is more, only little is known about how exactly the equilibrium expectation values emerge. Due to the exponentially many constants of motion present in quantum manybody systems, corresponding to the overlaps with the eigenvectors of the system, there seems to be no obvious reason why equilibrium values often only depend on few macroscopic properties such as temperature or particle number. In short: It is unclear how precisely the memory of the initial conditions is lost during time evolution.To make progress towards a solution of these two problems, it is instructive to study the behavior of non-interacting particles captured by so-called quadratic or free models. In these models the time evolution of so called Gaussian states, which are fully described by their correlation matrix, is particularly simple to describe. While studying...