We develop a technique to directly study spinons (emergent spin S = 1/2 particles) in quantum spin models in any number of dimensions. The size of a spinon wave packet and of a bound pair (a triplon) are defined in terms of wave-function overlaps that can be evaluated by quantum Monte Carlo simulations. We show that the same information is contained in the spin-spin correlation function as well. We illustrate the method in one dimension. We confirm that spinons are well defined particles (have exponentially localized wave packet) in a valence-bond-solid state, are marginally defined (with power-law shaped wave packet) in the standard Heisenberg critical state, and are not well defined in an ordered Néel state (achieved in one dimension using long-range interactions).PACS numbers: 75.10. Kt, 75.10.Jm, 75.40.Mg, 75.10.Pq Spinons are emergent spin S = 1/2 particles (fractional excitations) of quantum magnets [1][2][3] and potentially exist also in strongly-correlated electron systems such as the high-T c cuprate superconductors [4]. Their existence is well established in one-dimensional (1D) systems [1,2], where they correspond to kinks and antikinks (solitons). In higher dimensions, gapped magnons ("triplons") can be viewed as bound states of spinons. Under some conditions, in spin liquid states [3] and at certain quantumcritical points [5], these spinons may become deconfined (unbound). Even in cases where the spinons are not completely deconfined, such as in a valence-bond-solid (VBS) state of a two-dimensional (2D) system close to a phase transition into the antiferromagnetic (Néel) state, the bound state can become very large [5]. The spinons can then be viewed as deconfined below the length-scale of the pair size, and above a corresponding (relatively low) energy scale. This is analogous to quarks, which are the elementary constituent particles of the baryons although they are strictly speaking always confined.Observing deconfined or almost deconfined spinons in experiments is in general difficult [6]. In 1D systems, e.g., the Heisenberg chain, it is well understood (based on the exact Bethe ansatz solution and numerical calculations [2,7]) that spinons lead to a broad continuum in the dynamic spin structure factor S(q, ω). This continuum has been observed in neutron scattering experiments on quasi-1D quantum antiferromagnets [8]. In 2D systems, there is no known reference model with deconfined spinons in which S(q, ω) can be computed exactly. One nevertheless expects a broad continuum also in this case, and such experimental signatures have been claimed in some quasi-2D systems [9]. The issue is complicated, however, by the fact that a continuum is also expected due to multi-magnon processes [10].In this Letter, we discuss spinon detection in numerical model calculations. This has also been a challenging problem, the solution of which will greatly help to understand the conditions under which spinons can exist as independent elementary particles. Recently, signatures in thermodynamic properties were obser...