When a quantum phase transition is crossed in finite time, critical slowing down leads to the breakdown of adiabatic dynamics and the formation of topological defects. The average density of defects scales with the quench rate following a universal power-law predicted by the Kibble-Zurek mechanism. We analyze the full counting statistics of kinks and report the exact kink number distribution in the transverse-field quantum Ising model. Kink statistics is described by the Poisson binomial distribution with all cumulants exhibiting a universal power-law scaling with the quench rate. In the absence of finite-size effects, the distribution approaches a normal one, a feature that is expected to apply broadly in systems described by the Kibble-Zurek mechanism.Across a quantum phase transition, the equilibrium relaxation time diverges. This phenomenon, known as critical slowing down, is responsible for the nonadiabatic character of critical dynamics. Preparing the ground state of the brokensymmetry phase, an ubiquitous task in quantum science and technology, is thus intrinsically challenging: traversing the phase transition in finite time leads to the formation of topological defects. The Kibble-Zurek mechanism (KZM) is the paradigmatic theory to describe this scenario [1][2][3]. Its origins are found in the pioneering insight by Kibble on the role of causality in structure formation in the early universe [4,5]. Soon after, it was pointed out by Zurek that condensed-matter systems offer a test-bed to study the dynamics of symmetry breaking [6][7][8]. The key prediction of the KZM is that the average density d of the resulting topological defects scales with the quench time τ Q in which the phase transition is crossed as a universal power-law, d ∝ τ −α Q . The power-law exponent α = Dν/(1 + νz) is set by a combination of the dimensionality of the system D, and the dynamic and correlation-length (equilibrium) critical exponents denoted by z and ν, respectively.The validity of the KZM is however not restricted to the classical domain. The paradigmatic Landau-Zener formula, describing excitation formation in two-level systems, was shown to capture the KZM for long quench times [9,10]. As a result, paradigmatic models exhibiting quantum phase transitions, such as the 1D Ising chain, could be shown to obey the KZM, establishing the validity of the mechanism in the quantum domain for thermally isolated systems [9][10][11][12]. Due to its broad applicability, the KZM stands out as a result in statistical mechanics describing nonequilibrium properties (density of defects) in terms of equilibrium quantities (critical exponents). On the applied side, it suggest the need to pursue adiabatic strategies in quantum simulations as well as in quantum annealing, where the mechanism provides useful heuristics.Under unitary dynamics the state of the system following the crossing of the phase transition is characterized by collective and coherent quantum excitations. One can thus expect that even for isolated quantum systems, the order para...