We present the first formulation and application of relativistic unitary coupled cluster theory to atomic properties. The remarkable features of this theory are highlighted, and it is used to calculate the lifetimes of 5 2 D 3/2 and 6 2 P 3/2 states of Ba + and P b + respectively. The results clearly suggest that it is very well suited for accurate ab initio calculations of properties of heavy atomic systems.PACS numbers: 31.15. Ar, 31.15.Dv, 31.25.Jf, 32.10.Fn There have been a number of attempts to modify coupled-cluster (CC) theory [1], despite its spectacular success in elucidating the properties of a wide range of many-body systems [1,2,3,4,5]. One interesting case in point is unitary coupled-cluster (UCC) theory which was first proposed by Kutzelnigg [6]. In this theory, the effective Hamiltonian is Hermitian by construction and the energy which is the expectation value of this operator in the reference state is unlike in CC theory, an upper bound to the ground state energy [5]. Another attractive feature of this theory which we shall discuss later is that at a given level of approximation it incorporates certain higher order excitations that are not present in CC theory. Furthermore, it is well suited for the calculation of properties where core relaxation effects are important [7,8]. In spite of the aforementioned advantages, there have been relatively few studies based on this method [9,10]. This work is the first relativistic formulation of unitary coupled cluster (UCC) theory and also the first application of this theory to atomic properties. In this letter, we first present the formal aspects of relativistic UCC theory and then apply it to calculate the lifetimes of the 5 2 D 3/2 and 6 2 P 3/2 states of Ba + and P b + respectively; which depend strongly on both relativistic and correlation effects. The comparison of the results of these calculations with accurate experimental data would constitute an important test of this theory.The exact wave function for a closed shell state in CC theory is obtained by the action of the operator exp(T ) on the reference state |Φ . However, in UCC theory [10], it is written aswhere σ = T − T † ; and T and T † are the excitation and deexcitation operators respectively. σ is clearly antiHermitian, since σ † = −σ. Using this unitary ansatz for the correlated wave function, the relativistic UCC equation in the Dirac-Coulomb approximation can be written aswhere H is the Dirac-Coulomb HamiltonianUsing the normal ordered Hamiltonian, Eq.(2) can be rewritten aswhere the normal ordered Hamiltonian is defined as H N = H − Φ| H |Φ and ∆E = E − Φ| H |Φ . . The choice of the operator σ makes the effective Hamiltonian H N = exp(−σ)H N exp(σ) Hermitian.The effective Hamiltonian is expressed by the Hausdorff expansion in CC theory as