Relativistic DMRG calculations on the curve crossing of cesium hydride

Over the past few years, it has been shown in various studies on small molecules with only a few electrons that the density-matrix renormalization group (DMRG) method converges to results close to the full configuration-interaction limit for the total electronic energy. In order to test the capabilities of the method for molecules with complex electronic structures, we performed a study on the potential-energy curves of the ground state and the first excited state of 1sigma+ symmetry of the cesium hydride molecule. For cesium relativistic effects cannot be neglected, therefore we have used the generalized arbitrary-order Douglas-Kroll-Hess protocol up to tenth order, which allows for a complete decoupling of the Dirac Hamiltonian. Scalar-relativistic effects are thus fully incorporated in the calculations. The potential curves of the cesium hydride molecule feature an avoided crossing between the ground state and the first excited state, which is shown to be very well described by the DMRG method. Compared to multireference configuration-interaction results, the potential curves hardly differ in shape, for both the ground state and the excited state, but the total energies from the DMRG calculations are in general consistently lower. However, the DMRG energies are as accurate as corresponding coupled cluster energies at the equilibrium distance, but convergence to the full configuration-interaction limit is not achieved.