Coupled Cluster Method Research Articles

Photoinduced hydrogen dissociation in thymine predicted by coupled cluster theory

The fate of thymine upon excitation by ultraviolet radiation has been the subject of intense debate. Today, it is widely believed that its ultrafast excited state gas phase decay stems from a radiationless transition from the bright ππ* state to a dark nπ* state. However, conflicting theoretical predictions have made the experimental data difficult to interpret. Here we simulate the early gas phase ultrafast dynamics in thymine at the highest level of theory to date. This is made possible by performing wavepacket dynamics with a recently developed coupled cluster method. Our simulation confirms an ultrafast ππ* to nπ* transition (τ = 41 ± 14 fs). Furthermore, the predicted oxygen-edge X-ray absorption spectra agree quantitatively with experiment. We also predict an as-yet uncharacterized πσ* channel that leads to hydrogen dissociation at one of the two N-H bonds. Similar behavior has been identified in other heteroaromatic compounds, including adenine, and several authors have speculated that a similar pathway may exist in thymine. However, this was never confirmed theoretically or experimentally. This prediction calls for renewed efforts to experimentally identify or exclude the presence of this channel.

Linear Response pCCD-Based Methods: LR-pCCD and LR-pCCD+S Approaches for the Efficient and Reliable Modeling of Excited State Properties.

In this work, we derive working equations for the linear response pair coupled cluster doubles (LR-pCCD) ansatz and its extension to singles (S), LR-pCCD+S. These methods allow us to compute electronic excitation energies and transition dipole moments based on a pCCD reference function. We benchmark the LR-pCCD+S model against the linear response coupled-cluster singles and doubles method for modeling electronic spectra (excitation energies and transition dipole moments) of the BH, H2O, H2CO, and furan molecules. We also analyze the effect of orbital optimization within pCCD on the resulting LR-pCCD+S transition dipole moments and oscillator strengths and perform a statistical error analysis. We show that the LR-pCCD+S method can correctly reproduce the transition dipole moments features, thus representing a reliable and cost-effective alternative to standard, more expensive electronic structure methods for modeling electronic spectra of simple molecules. Specifically, the proposed models require only mean-field-like computational cost, while excited-state properties may approach the CCSD level of accuracy. Moreover, we demonstrate the capability of our model to simulate electronic transitions with non-negligible contributions of double excitations and the electronic spectra of polyenes of various chain lengths, for which standard electronic structure methods perform purely.

High-Level Calculation for Assessing the Atmospheric Reactivity of Pentachlorophenol with Hydroxyl Radical: Mechanism and Kinetics.

This work aims to investigate the reactivity of the pentachlorophenol (PCP) compound, C6Cl5OH, with a hydroxyl (OH) radical in the gas phase. Geometry optimizations and vibrational frequency calculations were performed using DFT-M06-2X/6-311++G(d,p). Single-point energy calculations were carried out using the single-reference coupled cluster method, specifically CCSD(T), and the multiconfigurational CASPT2 level of theory to characterize the atmospheric degradation processes of PCP with OH radicals and to identify which chemical species resulting from their decomposition could remain in the gas phase. The performance of various families of basis sets was tested. The widely used augmented-correlation-consistent basis set family aug-cc-pVXZ-DK (X = D, T, and Q) was compared to the atomic natural orbital basis sets ANO-RCC-VXZP (X = D, T, and Q). The energy profile at 298 K showed that the Cl- and OH-abstractions are not energetically favorable. H-abstraction and OH-addition/Ck (k = 1-6) are characterized by the lowest Gibbs energies of activation and are strongly exothermic. The canonical transition state theory (TST) with a simple Wigner tunneling correction is used to predict the rate constants over the 220-400 K temperature range for each reaction channel. The overall rate constant at 298 K based on our calculations is about 3.25 × 10-13, 3.10 × 10-15, and 2.21 × 10-15 cm3 molecule-1 s-1 at M06-2X, CASPT2/ANO-RCC-VQZP, and CASPT2/aug-cc-pVQZ-DK levels of theory, respectively. The PCP atmospheric lifetime at 298 K is predicted to be in the range from 1 to 12-16 years at 0 km in the presence of variable OH concentration (9.00 × 105 to 1.50 × 107 molecules cm-3) based on CASPT2 data.

An explicitly correlated potential energy surface for N2-OCS complex: Out-of-plane motion and tunneling dynamics.

A four-dimensional potential energy surface (4D-PES) has been constructed for the N2-OCS complex. The PES is achieved by applying the explicitly correlated coupled cluster method, which incorporates single, double, and perturbative triple excitations [CCSD(T)-F12a], along with the augmented correlation consistent triple zeta (aug-cc-pVTZ) basis set. The rovibrational levels are precisely determined and assigned through bound state calculations and wavefunction analysis. The calculated transition frequencies reproduce the experimental observations accurately, achieving an RMSE of 0.0005cm-1 for the 23 rotational transitions (J ≤ 6, Ka ≤ 2). The R-φ contour plot of the wave function clearly demonstrates the unambiguous delocalization of the dihedral angle, and the averaged geometry of the ground vibrational state is determined to be non-planar with φ = 90°. To obtain a quantitative analysis of this phenomenon, we expanded the 3H-solution model [Guo et al., J. Quant. Spectrosc. Radiat. Transfer 309 (2023) 108711] from a three-dimensional system (Ar-AgF) to a nine-dimensional system (N2-OCS). Based on this model, the tunneling splitting was calculated to be 0.0822cm-1, which excellently matches the experimental result of 0.0817cm-1. The excellent agreement between the theoretical and experimental results suggests that the wavefunction delocalization and out-of-plane motion can be attributed to the tunneling effects in the ground vibrational state.

Quantum Chemistry Study on Cl-Initiated Reactions of 2-Chloropropane and 2-Methylpropanoyl Halogen (Cl, Br, F): Mechanism, Kinetics, and Atmospheric Implications.

Halogenated volatile organic compounds (HVOCs) pose significant bioaccumulative and toxicological risks, necessitating effective strategies for their removal. Here, we show, through a computational study employing density functional theory and coupled cluster methods, the detailed mechanism and kinetic properties of Cl-initiated degradation reactions of 2-chloropropane (2-CP, (CH3)2CHCl) and 2-methylpropanoyl halide ((CH3)2CHCOX, X = Cl, Br, F). The reaction rate constants of all the channels were calculated by the canonical variational transition state theory (CVT) with the correction of the small curvature tunneling effect (SCT) at 200-1000 K. The subsequent transformation pathways of the major radical products of (CH3)2CHCl and (CH3)2CHCOCl in the presence of O2, NO, and HO2 radical were investigated. The results elucidate the reaction pathways and rate constants, which are in excellent agreement with the experimental data at 296 K. We further explore the atmospheric implications of these reactions by assessing the atmospheric lifetime (τ) and ozone depletion potential (ODP). Additionally, we delve into the aquatic toxicity and bioaccumulation potential of the reactants and their transformation products. This study not only advances our knowledge of the atmospheric fate of halogenated hydrocarbons but also underscores the importance of considering the environmental and toxicological impacts in the development of HVOC mitigation strategies.

Accurate Thermochemistry with Multireference Methods: A Stress Test for Internally Contracted Multireference Coupled-Cluster Theory.

The internally contracted multireference coupled-cluster method with single, double and perturbative triple excitations, icMRCCSD(T), was tested for its performance in the context of computational high-accuracy thermochemistry. The results were gauged against the standard single-reference coupled-cluster hierarchy with up to 5-fold excitations. The test set comprised of a selection of first-row dinuclear compounds and the three 3d-transition metal compounds MnH, FeH, and CoH. The results revealed two problems with the current formulation of icMRCCSD(T). First, the choice of the Dyall Hamiltonian as the zeroth-order Hamiltonian, which leads to a biased description of the different orbital subspaces and particularly poor results for the atomic correlation energies, and second, the tendency to overestimate the perturbative correction for triply excited clusters, in particular in the presence of open shells and correspondingly low orbital-energy gaps. The two problems could be solved by resorting to the effective Fock operator as zeroth-order Hamiltonian and by adopting a modified amplitude equation that includes terms quadratic in the pair clusters. A similar modification was recently proposed by Masios et al. (Phys. Rev. Lett. 2023, 131, 186401) in the context of applying single-reference coupled-cluster theory to systems with small or vanishing band gaps and we chose the acronym '(cT*) correction' in analogy to that work. In contrast to the work of Masios et al., additional terms including single excitation clusters were omitted, as these again lead to an overestimation of correlation effects in more difficult cases. We also tested another alternative for the zeroth-order Hamiltonian and additional higher-order corrections for the correlation energy. These extensions did not significantly improve the results and were also computationally more demanding. The improved icMRCCSD(cT*)F method yields very accurate results with errors, relative to accurate benchmarks, better than 2 kJ/mol for total energies and atomization energies for the entire set of examples considered in this work.

An "ultimate" coupled cluster method based entirely on T2.

Electronic structure methods built around double-electron excitations have a rich history in quantum chemistry. However, it seems to be the case that such methods are only suitable in particular situations and are not naturally equipped to simultaneously handle the variety of electron correlations that might be present in chemical systems. To this end, the current work seeks a computationally efficient, low-rank, "ultimate" coupled cluster method based exclusively on T2 and its products that can effectively emulate more "complete" methods that explicitly consider higher-rank, T2m, operators. We introduce a hierarchy of methods designed to systematically account for higher, even order cluster operators, such as T4, T6, …, T2m, by invoking tenets of the factorization theorem of many-body perturbation theory (MBPT) and expectation-value coupled cluster theory. It is shown that each member within this methodological hierarchy is defined such that both the wavefunction and energy are correct through some order in MBPT and can be extended up to arbitrarily high orders in T2. The efficacy of such approximations are determined by studying the potential energy surface of several closed and open-shell molecules. We find that the proposed hierarchy of augmented T2 methods essentially reduces to standard CCD for problems where dynamic electron correlations dominate but offer improvements in situations where non-dynamic and static correlations become relevant. A notable highlight of this work is that the cheapest methods in this hierarchy-which are correct through fifth-order in MBPT-consistently emulate the behavior of the O(N10) CCDQ method, yet only require a O(N6) algorithm by virtue of factorized intermediates.

Automatic Orbital Pair Selection for Multilevel Local Coupled-Cluster Based on Orbital Maps.

We present an automatic, orbital-map based orbital-pair selection scheme for multilevel local coupled-cluster approaches that exploits the locality of chemical reactions by focusing on the part of the molecule directly involved in the reaction. The previously introduced pair-selected multilevel extension to domain-based local pair natural orbital coupled-cluster with singles, doubles, and semicanonical perturbative triples [DLPNO-CCSD(T0)] partitions the orbital pairs according to relative changes in pair correlation energies [Bensberg, M.; Neugebauer, J. Orbital pair selection for relative energies in the domain-based local pair natural orbital coupled-cluster method. J. Chem. Phys. 2022, 157, 064102. 10.1063/5.0100010]. To this end, maps between localized orbitals are required which in turn require maps between the atoms of structures along reaction paths. So far, these atom maps have been manually determined, which can be a (human) time-consuming procedure. Here, we present an automatic atom mapping algorithm based on the principle of minimum chemical distance that incorporates orientation dependence through dihedral angles. A similar strategy is then introduced to obtain orbital maps, which proves advantageous over the previously used direct orbital selection. Along with a modified orbital pair prescreening, this results in an improved variant of the pair-selected multilevel DLPNO-CCSD(T0) method. The performance of this approach is demonstrated for various reaction types showing a significant efficiency gain and accurate results due to beneficial, systematic error cancellation. The presented method operates in a black-box manner due to its fully automatized algorithms with only the need to specify a single target-accuracy parameter. Additionally, we demonstrate that basis set extrapolation techniques can be applied. In this context, the approach shows deficiencies for the use of large basis sets, especially with diffuse basis functions, which can be traced back to the semicanonical triples correction.

A static quantum embedding scheme based on coupled cluster theory.

We develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller-Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller-Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C-C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.

Computational Studies of Nucleophilic Substitution at Nitrogen Center: Reactions of NH2Cl with HO-, CH3O- and C2H5O.

The atomic-level mechanisms of the nucleophilic substitution reactions at the nitrogen center (SN2@N) were investigated for the reactions of chloramine (NH2Cl) with the alkoxide ions (RO-, where R=H, CH3, and C2H5) using DFT and MP2 methods. The computed potential energy profiles for the SN2@N pathways involving the back-side attack of the nucleophiles show the typical double-well potential with submerged barriers similar to the SN2 reactions at the carbon center (SN2@C). However, the pre-reaction and post-reaction complexes are, respectively, the N-H⋅⋅⋅O and N-H⋅⋅⋅Cl hydrogen-bonded intermediates, which are different from those generally seen in SN2@C reactions. The SN2@N pathways involving front-side attack of the nucleophiles have high-energy barriers. The potential energy surfaces (PESs) along the proton-transfer pathways were flat. In addition to the proton-transfer and SN2 pathways, we also observed a new path for the methoxide and ethoxide nucleophiles where a hydride-transfer from the nucleophile to chloramine resulted in the products Cl-+R'CHO+NH3, (R'=H, CH3), and was the most exoergic. A comparison of the energetics obtained used different DFT and MP2 methods with that of the benchmark coupled-cluster methods reveals that CAM-B3LYP best describes the PESs.

Rotational excitation of methyl mercaptan (CH3SH) in collisions with molecular hydrogen

ABSTRACT This paper presents the calculation of rate coefficients for transitions between rotational levels of the A-type and E-type levels of methyl mercaptan (CH$_3$SH), resulting from collisions with molecular hydrogen. Radiative transfer modelling requires both radiative and collisional rates to describe the rotational populations under the usual conditions in interstellar clouds where local thermodynamic equilibrium conditions do not apply. To compute the intermolecular interaction between CH$_3$SH and H$_2$, the explicitly correlated CCSD(T)-F12a coupled-cluster method that utilized a correlation-consistent aug-cc-pVTZ basis was employed. The computed energies were fit to a functional form suitable for use in scattering calculations. Rate coefficients were calculated over the temperature range from 5 to 100 K for transitions between the 110 lowest CH$_3$SH rotational levels (having energies less than 107 cm$^{-1}$ (ca. 150 K) within both the A-type and E-type manifolds caused by collisions with para- and ortho-H$_2$. The rate coefficients were obtained through time-independent quantum close coupling quantum scattering calculations utilizing the calculated potential energy surface.

DMRG-Tailored Coupled Cluster Method in the 4c-Relativistic Domain: General Implementation and Application to the NUHFI and NUF3 Molecules.

Heavy atom compounds represent a challenge for computational chemistry due to the need for simultaneous treatment of relativistic and correlation effects. Often such systems also exhibit strong correlation, which hampers the application of perturbation theory or single-reference coupled cluster (CC) methods. As a viable alternative, we have proposed externally correcting the CC method using the density matrix renormalization group (DMRG) wave functions, yielding the DMRG-tailored CC method. In a previous paper [J. Chem. Phys. 2020, 152, 174107], we reported a first implementation of this method in the relativistic context, which was restricted to molecules with real double group symmetry. In this work, we present a fully general implementation of the method, covering complex and quaternion double groups as well. The 4c-TCC method thus becomes applicable to polyatomic molecules, including heavy atoms. For the assessment of the method, we performed calculations of the chiral uranium compound NUHFI, which was previously studied in the context of the enhancement of parity violation effects. In particular, we performed calculations of a cut of the potential energy surface of this molecule along the stretching of the N-U bond, where the system exhibits strong multireference character. Since there are no experimental data for NUHFI, we have performed also an analogous study of the (more symmetric) NUF3 molecule, where the vibrational frequency of the N-U bond can be compared with spectroscopic data.

Two-Photon Absorption Strengths of Small Molecules: Reference CC3 Values and Benchmarks.

We present a large dataset of highly accurate two-photon transition strengths (δTPA) determined for standard small molecules. Our reference values have been calculated using the quadratic response implementation of the third-order coupled cluster method including iterative triples (Q-CC3). The aug-cc-pVTZ atomic basis set is used for molecules with up to five non-hydrogen atoms, while larger molecules are assessed with aug-cc-pVDZ; the differences due to the basis sets are discussed. This dataset, encompassing 82 singlet transitions of various characters (Rydberg, valence, and double excitations), enables a comprehensive benchmark of smaller basis sets and alternative wavefunction methods when Q-CC3 calculations become beyond reach as well as time-dependent density functional theory (TD-DFT) approaches. The evaluated wavefunction methods include quadratic response and equation-of-motion CCSD approximations, Q-CC2, and second-order algebraic diagrammatic construction in its intermediate state representation (I-ADC2). In the TD-DFT framework, a set of five commonly used exchange-correlation functionals are evaluted. This extensive analysis provides a quantitative assessment of these methods, revealing how different system sizes, response intensities, and types of transitions affect their performances.

Rank-Reduced Equation-of-Motion Coupled Cluster Triples: an Accurate and Affordable Way of Calculating Electronic Excitation Energies.

In the present work, we report an implementation of the rank-reduced equation-of-motion coupled cluster method with approximate triple excitations (RR-EOM-CC3). The proposed variant relies on tensor decomposition techniques in order to alleviate the high cost of computing and manipulating the triply excited amplitudes. In the RR-EOM-CC3 method, both ground-state and excited-state triple-excitation amplitudes are compressed according to the Tucker-3 format. This enables factorization of the working equations such that the formal scaling of the method is reduced to N6, where N is the system size. An additional advantage of our method is the fact that the accuracy can be strictly controlled by proper choice of two parameters defining sizes of triple-excitation subspaces in the Tucker decomposition for the ground and excited states. Optimal strategies of selecting these parameters are discussed. The developed method has been tested in a series of calculations of electronic excitation energies and compared to its canonical EOM-CC3 counterpart. Errors several times smaller than the inherent error of the canonical EOM-CC3 method (in comparison to FCI) are straightforward to achieve. This conclusion holds both for valence states dominated by single excitations and for states with pronounced doubly excited character. Taking advantage of the decreased scaling, we demonstrate substantial computational costs reductions (in comparison with the canonical EOM-CC3) in the case of two large molecules - l-proline and heptazine. This illustrates the usefulness of the RR-EOM-CC3 method for accurate determination of excitation energies of large molecules.

Complete Active Space Iterative Coupled Cluster Theory.

In this work, we investigate the possibility of improving multireference-driven coupled cluster (CC) approaches with an algorithm that iteratively combines complete active space (CAS) calculations with tailored CC and externally corrected CC. This is accomplished by establishing a feedback loop between the CC and CAS parts of a calculation through a similarity transformation of the Hamiltonian with those CC amplitudes that are not encompassed by the active space. We denote this approach as the complete active space iterative coupled cluster (CASiCC) ansatz. We investigate its efficiency and accuracy in the singles and doubles approximation by studying the prototypical molecules H4, H8, H2O, and N2. Our results demonstrate that CASiCC systematically improves on the single-reference CCSD and the externally corrected CCSD methods across entire potential energy curves while retaining modest computational costs. However, the tailored coupled cluster method shows superior performance in the strong correlation regime, suggesting that its accuracy is based on error compensation. We find that the iterative versions of externally corrected and tailored coupled cluster methods converge to the same results.

Low-lying [formula omitted] electronic states of BiF: Perturbative versus variational approaches of spin–orbit coupling

The low-lying Ω states of bismuth monofluoride (BiF) below 50000 cm−1 have been theoretically studied using the multi-reference configuration interaction and the equation-of-motion coupled-cluster methods, where the spin–orbit coupling (SOC) effects are respectively considered perturbatively or variationally. Since the perturbative treatment of SOC is not a good approximation for the heavy atom Bi, the latter method exhibits higher accuracy and better agreements with the experimental spectroscopic constants. Our results show that the second 3Π state of BiF is due to the occupation on the Rydberg 7s-shell of Bi, which gives rise to the so-called “triplet C” state and its sub-states reported in the early literatures. With the help of our theoretical results, some experimentally assigned spectral bands with considerable confusion have also been clarified and reassigned.

Triple Electron Attachments with a New Intermediate-Hamiltonian Fock-Space Coupled-Cluster Method.

The implementation of a new intermediate-Hamiltonian Fock-space coupled-cluster (IHFSCC) scheme for the (3,0) sector of the Fock space is reported. In this IHFSCC approach, the three-body contributions in the cluster operator S(3,0) corresponding to the (3,0) sector of the Fock space are considered, while S(1,0) and S(2,0) at the (1,0) and (2,0) level only include one- and two-body operators. By introducing a suitable partition of the wave operator, the intermediate Hamiltonian, which only depends on the two-body operator of S(1,0), is obtained. S(2,0) and S(3,0) are not required within this new IHFSCC scheme, and a large reference space can be possibly employed. The performance of this (3,0) IHFSCC method in calculating triple ionization potentials and excitation energies for atoms and cations with a ground p3 configuration as well as adiabatic excitation energies for some molecules is investigated. The effect of the number of active virtual orbitals and three different types of orbitals, i.e., reference orbitals, restricted open-shell Hartree-Fock orbitals (ROHF) of the target systems, and canonicalized ROHF orbitals, on IHFSCC results, is also studied. Our calculations indicate that reasonable results can be achieved with this (3,0) IHFSCC method when a minimal reference space is employed. Further increasing the number of active orbitals does not necessarily improve the results. In addition, the IHFSCC results using canonicalized ROHF orbitals generally agree well with reference values, and they are not very sensitive to the number of active orbitals compared with results using the reference orbitals. The new (3,0) IHFSCC method can be applied to open-shell systems with three unpaired electrons with reasonable accuracy at a relatively low computational cost.

Time-dependent orbital-optimized coupled-cluster methods families for fermion-mixtures dynamics.

Five time-dependent orbital optimized coupled-cluster methods, of which four can converge to the time-dependent complete active space self-consistent-field method, are presented for fermion-mixtures with arbitrary fermion kinds and numbers. Truncation schemes maintaining the intragroup orbital rotation invariance, as well as equations of motion of coupled-cluster (CC) amplitudes and orbitals, are derived. Present methods are compact CC-parameterization alternatives to the time-dependent multiconfiguration self-consistent-field method for systems consisting of arbitrarily different kinds and numbers of interacting fermions. Theoretical analysis of applications of present methods to various chemical systems is reported.

Benchmark Study of Core-Ionization Energies with the Generalized Active Space-Driven Similarity Renormalization Group.

X-ray photoelectron spectroscopy (XPS) is a powerful experimental technique for probing the electronic structure of molecules and materials; however, interpreting XPS data requires accurate computational methods to model core-ionized states. This work proposes and benchmarks a new approach based on the generalized active space-driven similarity renormalization group (GAS-DSRG) for calculating core-ionization energies and treating correlation effects at the perturbative and nonperturbative levels. We tested the GAS-DSRG across three data sets. First, the vertical core-ionization energies of small molecules containing first-row elements are evaluated. GAS-DSRG achieves mean absolute errors below 0.3 eV, which is comparable to high-level coupled cluster methods. Next, the accuracy of GAS-DSRG is evaluated for larger organic molecules using the CORE65 data set, with the DSRG-MRPT3 level yielding a mean absolute error of only 0.34 eV for 65 core-ionization transitions. Insights are provided into the treatment of static and dynamic correlation, the importance of high-order perturbation theory, and notable differences from density functional theory in the predicted energy ordering of core-ionized states for specific molecules. Finally, vibrationally resolved XPS spectra of diatomic molecules (CO, N2, and O2) are simulated, showing excellent agreement with experimental data.

Calculation of Some Low-Lying Electronic Excitations of Barium Monofluoride Using the Equation of Motion (EOM)-CC3 Method with an Effective Core Potential Approach.

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.