In this talk, I will present our recent work on scattering of two boson stars by taking into account three effects: point-mass gravitational, tidal, and short-range scalar-field interactions. We first derive the analytical expressions of the scattering angle using PM-EFT techniques, providing the first analytical treatment of boson stars as a two-body problem. We then compare analytic results to the scattering angle extracted from four sequences of numerical-relativity simulations at fixed energy and varying impact parameter. The very good agreement exhibits the attractive (repulsive) effect of in-phase (out-of-phase) binaries. For small impact parameters, where the stars approach more closely before separating to infinity, the short-range scalar-field interaction is found to dominate.