We report on the experimental realization of a photonic system that simulates the one-dimensional two-particle Hubbard model. This analogy is realized by means of two-dimensional arrays of coupled optical waveguides, fabricated using femtosecond laser inscription. By tuning the analogous "interaction strength", we reach the strongly-interacting regime of the Hubbard Hamiltonian, and demonstrate the suppression of standard tunneling for individual "particles". In this regime, the formation of bound states is identified through the direct observation of pair tunneling. We then demonstrate the coherent destruction of tunneling (CDT) for the paired particles in the presence of an engineered oscillating force of high frequency. The precise control over the analogous "interaction strength" and driving force offered by our experimental system opens an exciting route towards quantum simulation of few-body physics in photonics.Introduction. Elucidating the physics of interacting electrons in crystalline solids constitutes one of the most challenging problems in modern physics, with direct implications for our understanding of quantum magnetism and superconductivity. Theoretical toy models, such as the Hubbard model [1], can be used to examine these problems, but challenging issues can arise, especially in the intermediate to strong coupling regimes where perturbation theory fails. However, the one-dimensional Hubbard model is exactly solvable by means of Bethe Ansatz techniques [2,3]. In particular, the two-particle solution in the singlet sector -the triplet sector is noninteracting -shows both scattering and bound states for any value of the interaction strength. Remarkably, bound state solutions, known as doublons, exist even in the presence of repulsive interactions. This phenomenon can be understood in terms of the band gap, or the boundedness of the spectrum in the Hubbard model, which implies that two repulsively interacting particles, initially occupying the same lattice well, will have no available scattering energies to dissociate into. This can be further explained by an exact symmetry between attractive and repulsive interactions reported in [4][5][6]. These repulsively bound states were experimentally observed using both bosonic [7] and fermionic [8] particles in optical lattices. These experiments triggered an intense activity exploring the physics of few particles in optical lattices, including the two-body [9][10][11][12][13][14][15][16][17][18][19] and three-body [20,21] problems, in which certain phenomena that have no analogy in free space occur. For example, the Mattis-Gallinar effect [22], due to the absence of Galilean invariance, states that the effective mass of a lattice bound pair is higher than its free-space analogue, and was observed for excitons [23]. Three-body composites in bosonic [21] or mass-imbalanced one-dimensional fermionic systems [20] can also exist due to an effective exchange mechanism