The phenomenon of pairing in atomic and nuclear many-body systems gives rise to a great number of different physical properties of matter, from areas as seemingly diverse as the shape of stable nuclei to superconductivity in metals and superfluidity in neutron stars. With the experimental realisation of the long sought BCS-BEC crossover observed in trapped atomic gases - where it is possible to fine tune the s-wave scattering length a of a many-fermion system between a dilute, correlated BCS-like superfluid of Cooper pairs and a densely packed BEC of composite bosons - pairing problems in atomic physics have found renewed interest in recent years. Given the high precision techniques involved in producing these trapped gas condensates, we would like to employ a suitably accurate many-body method to study such systems, preferably one which goes beyond the simple mean-field picture.The Coupled Cluster Method (CCM) is a widely applied and highly successful ab initio method in the realm of quantum many-body physics and quantum chemistry, known to be capable of producing extremely accurate results for a wide variety of different many-body systems. It has not found many applications in pairing problems however, at least not in a general sense. Our aim, therefore, is to study various models of pairing using a variety of CCM techniques - we are interested in studying the generic features of pairing problems and in particular, we are especially interested in probing the collective modes of a system which exhibits the BCS-BEC crossover, in either the BCS or BEC limit. The CCM seems a rather good candidate for the job, given the high precision results it can produce.

In this thesis, the implementation of the incremental scheme for open-shell systems with unrestricted Hartree-Fock reference wave functions is described. The implemented scheme is tested within robustness and performance with respect to the accuracy in the energy and the computation times. New approaches are discussed to implement a fully automated incremental scheme in combination with the domain-specific basis set approximation. The alpha Domain Partitioning and Template Equalization are presented to handle unrestricted wave functions for the local correlation treatment. Both orbital schemes are analyzed with a test set of structures and reactions. As a further goal, the DSBSenv orbital basis sets and auxiliary basis sets are optimized to be used as environmental basis in the domain-specific basis set approach. The performance with respect to the accuracy and computation times is analyzed with a test set of structures and reactions. In another project, a scheme for the optimization of auxiliary basis sets for uranium is presented. This scheme was used to optimize the MP2Fit auxiliary basis sets for uranium. These auxiliary basis enable density fitting in quantum chemical methods and the application of the incremental scheme for systems containing uranium. Another project was the systematical analysis of the binding energies of four water dodecamers. The incremental scheme in combination with the CCSD(T) and CCSD(T)(F12*) method were used to calculate benchmark energies for these large clusters.