Recently developed explicitly correlated local coupled-cluster methods [PNO-LCCSD(T)-F12] are reviewed. Extensive benchmarks for reaction energies and intermolecular interaction energies are presented, in which the convergence of the results with respect to all local approximations is studied. The explicit correlation treatment (F12) is shown to be essential to minimize basis set incompleteness errors, as well as errors caused by domain approximations. Generally, the errors of relative energies due to local approximations can be reduced to below 1 kcal/mol. The methods are well parallelized, and using small computer clusters with 100-200 computing cores, calculations for systems with 100-200 atoms using augmented triple-ζ basis sets can be carried out within a few hours of elapsed time. Recommendations are made on how such calculations should be carried out, how the accuracy can be tested, and which computational resources are required. This article is categorized under: Electronic Structure Theory> Ab Initio Electronic Structure Methods Software> Quantum Chemistry K E Y W O R D S explicit correlation, coupled cluster, local correlation, pair natural orbitals 1 | INTRODUCTIONThe coupled-cluster method with single and double excitations and a perturbative treatment of triple excitations [CCSD(T)] is considered the gold standard of quantum chemistry for computing energies and other properties of molecules. Often, the accuracy is comparable or even better than that of experimental data, and chemical accuracy (1 kcal/mol or better for relative energies) can be achieved, provided that the system under consideration is of single-reference character, that is, the Hartree-Fock (HF) Slater determinant dominates the total wave function expansion. However, in its traditional canonical form, the application of the CCSD(T) method is limited to rather small systems (20-30 atoms) due to the steep increase of the computational effort with system size: The central processing unit (CPU) time scales as O N 7 el