An efficient finite difference/Hermite-Galerkin spectral method for time-fractional coupled sine-Gordon equations on multidimensional unbounded domains and its application in numerical simulations of vector solitons.