The amplitude ("Higgs") mode is a ubiquitous collective excitation related to spontaneous breaking of a continuous symmetry. We combine quantum Monte Carlo (QMC) simulations with stochastic analytic continuation to investigate the dynamics of the amplitude mode in a three-dimensional dimerized quantum spin system. We characterize this mode by calculating the spin and dimer spectral functions on both sides of the quantum critical point, finding that both the energies and the intrinsic widths of the excitations satisfy field-theoretical scaling predictions. While the line width of the spin response is close to that observed in neutron scattering experiments on TlCuCl3, the dimer response is significantly broader. Our results demonstrate that highly non-trivial dynamical properties are accessible by modern QMC and analytic continuation methods.The spontaneous breaking of a continuous symmetry allows collective excitations of the direction and amplitude of the order parameter; for O(N ) symmetry, there are N−1 massless directional (Goldstone) modes and one massive amplitude mode [1][2][3][4]. In loose analogy with the Standard Model, the latter is often called a Higgs mode. A strongly damped amplitude mode has been reported in two dimensions (2D) at the Mott transition of ultracold bosons [5] and at the disorder-driven superconductorinsulator transition [6,7]. In 3D, the amplitude mode is expected on theoretical grounds to be more robust, and indeed the cleanest observation to date of a "Higgs boson" in condensed matter is at the pressure-induced magnetic quantum phase transition (QPT) in the dimerized quantum antiferromagnet TlCuCl 3 [8][9][10].Below the upper critical number of space-time dimensions, which for an O(N ) model is D c = 4, the amplitude mode is unstable, decaying primarily into pairs of Goldstone bosons [11][12][13]. In both 2D and 3D, the longitudinal dynamic susceptibility exhibits an infrared singularity due to the Goldstone modes [14], whose consequences for the visibility of the amplitude mode have been investigated extensively in 2D [15][16][17]. It was noted [14] that the scalar O(N )-symmetric susceptibility remains uncontaminated by infrared contributions, which should permit the amplitude mode to be observed as a welldefined peak. The (3+1)D O(3) case of TlCuCl 3 is at D c and the amplitude mode is critically damped, meaning that its width is proportional to its energy at the meanfield level [9,[18][19][20]. This mode can be probed through the spin response (longitudinal susceptibility) by neutron spectroscopy, and measurements over a wide range of pressures reveal a rather narrow peak width of just 15% of the excitation energy [10]. The value of this nearconstant width-to-energy ratio is the key to the mode visibility, thus calling for unbiased numerical calculations in suitable model Hamiltonians.FIG. 1. Schematic representation of ground states, excitation processes, and corresponding gaps in a dimerized antiferromagnet. The ratio g = J /J of the intra-and inter-dimer coupling constants ...