  We present a general approach within the relativistic coupled-cluster theory
framework to calculate exactly the first order wave functions due to any rank
perturbation operators. Using this method, we calculate the static dipole and
quadrupole polarizabilities in some alkali atoms and alkaline earth-metal ions.
This may be a good test of the present theory for different rank and parity
interaction operators. This shows a wide range of applications including
precise calculations of both parity and CP violating amplitudes due to rank
zero and rank one weak interaction Hamiltonians. We also give contributions
from correlation effects and discuss them in terms of lower order many-body
perturbation theory.
