Relativistic coupled-cluster calculations of nuclear spin-dependent parity non-conservation in Cs, Ba$^+$ and Ra$^+$

We have developed a relativistic coupled-cluster theory to incorporate nuclear spin-dependent interaction Hamiltonians perturbatively. This theory is ideal to calculate parity violating nuclear spin-dependent electric dipole transition amplitudes, $E1_{\rm PNC}^{\rm NSD}$, of heavy atoms. Experimental observation of which is a clear signature of nuclear anapole moment, the dominant source of nuclear spin-dependent parity violation in atoms and ions. We apply the theory to calculate $E1_{\rm PNC}^{\rm NSD}$ of Cs, which to date has provided the best atomic parity violation measurements. We also calculate $E1_{\rm PNC}^{\rm NSD}$ of Ba$^+$ and Ra$^+$, candidates of ongoing and proposed experiments.