Barium monofluoride (BaF) is a polar molecule of interest in measurements of the electron electric dipole moment. For this purpose, efforts are underway to investigate this molecule embedded within cryogenic matrices, e.g., in solid Ne. For a theoretical understanding of the electronic structure of such an embedded molecule, the need arises for efficient methods which are accurate but also able to handle a number of atoms which surround the molecule. The calculation for gas-phase BaF can be reduced to involve only outer electrons by representing the inner core of Ba with a pseudopotential, while carrying out a non-relativistic calculation with an appropriate basis set. Thus, the method is effectively at a scalar-relativistic level. In this work, we demonstrate to which extent this can be achieved using coupled-cluster methods to deal with electron correlation. As a test case, the SrF(X2Σ+→B2Σ+) transition is investigated, and excellent accuracy is obtained with the EOM-CC3 method. For the BaF(X2Σ+→A'2Δ, X2Σ+→A2Π, X2Σ+→B2Σ+) transitions, various coupled-cluster approaches are compared with very good agreement for EOM-CC3 with experimentally derived spectroscopic parameters, at the level of tens of cm-1. An exception is the excitation to the A'2Δ state, for which the energy is overestimated by 230cm-1. The poor convergence behavior for this particular state is demonstrated by providing results from calculations with basis sets of n = 3, 4, 5)-zeta quality. The calculated excitation energy for the B2Σ+ state agrees better with a deperturbation analysis than with the effective spectroscopic value, with a difference of 120cm-1.
Read full abstract