There's actually more than one way to calculate the dipole moment, and you can check whether or not the numbers you're obtaining from td=(singlets,nstates=5,root=1)
are making sense by comparing your result to what you'd get with a more "manual" approach in which you can see exactly what's going on.
The simplest formula is simply the expectation value of the dipole moment operator:
$$\tag{1}
\langle \psi | \mathbf{\hat{\mu} }| \psi\rangle.
$$
However, it can be much harder to get wavefunctions than energies, and a first-order error in the wavefunction can lead to a second-order error in the energy, so it would be convenient if we could calculate the dipole moment via only energy calculations (also, some electronic structure methods or implementations in software may be easily able to give energies but not wavefunctions). In fact we can calculate the dipole moment with only simple energy calculations.
It can be shown that for an electric field $\mathcal{E}$ we have (see section 4 of this PDF, for example):
$$\langle \psi | \mathbf{\hat{\mu} }| \psi\rangle = -\frac{\textrm{d}E}{\textrm{d}\mathcal{E}}\tag{2}.$$
This is related to perturbation theory and/or linear response theory, and it allows us to calculate Eq. 1 without ever needing the wavefunction, simply by calculating the energy multiple times and doing finite differences to approximate the derivative (if you want $n$ components of the dipole moment, you'll need to do $n+1$ extra energy calculations in addition to the zero-field energy calculation, but even for 3-components of the dipole moment, I often wouldn't mind doing 5 energy calculations if it could help me avoid requiring use of the wavefunction).
Whenever I hear "finite differences" I cringe a bit since I almost always prefer to do things analytically if possible, but this paper by Ernzerhof (the 'E' in PBE) et al. actually explains that using finite differences can be better than trying to calculate the dipole moment more directly. The same is said in section 3.1 of this paper by Lodi & Tennyson, along with more references that recommend the finite differences approach over the more "direct" approach. You can see the differences between numerical values obtained via finite differences versus DFPT here and the difference between finite differences and using the linear response here.
So if you're not sure if the direct approach is giving you sensible answers, you can just do a simple energy calculation $2 + n$ times if you want $n$ components of the dipole moment, and it might not be a bad idea to do this second method anyway, to see whether or not the more "direct" approach is giving sensible results!