The MRCC program system: Accurate quantum chemistry from water to proteins.

MRCC is a package of ab initio and density functional quantum chemistry programs for accurate electronic structure calculations. The suite has efficient implementations of both low- and high-level correlation methods, such as second-order Møller-Plesset (MP2), random-phase approximation (RPA), second-order algebraic-diagrammatic construction [ADC(2)], coupled-cluster (CC), configuration interaction (CI), and related techniques. It has a state-of-the-art CC singles and doubles with perturbative triples [CCSD(T)] code, and its specialties, the arbitrary-order iterative and perturbative CC methods developed by automated programming tools, enable achieving convergence with regard to the level of correlation. The package also offers a collection of multi-reference CC and CI approaches. Efficient implementations of density functional theory (DFT) and more advanced combined DFT-wave function approaches are also available. Its other special features, the highly competitive linear-scaling local correlation schemes, allow for MP2, RPA, ADC(2), CCSD(T), and higher-order CC calculations for extended systems. Local correlation calculations can be considerably accelerated by multi-level approximations and DFT-embedding techniques, and an interface to molecular dynamics software is provided for quantum mechanics/molecular mechanics calculations. All components of MRCC support shared-memory parallelism, and multi-node parallelization is also available for various methods. For academic purposes, the package is available free of charge.

An efficient linear-scaling CCSD(T) method based on local natural orbitals.

An improved version of our general-order local coupled-cluster (CC) approach [Z. Rolik and M. Kállay, J. Chem. Phys. 135, 104111 (2011)] and its efficient implementation at the CC singles and doubles with perturbative triples [CCSD(T)] level is presented. The method combines the cluster-in-molecule approach of Li and co-workers [J. Chem. Phys. 131, 114109 (2009)] with frozen natural orbital (NO) techniques. To break down the unfavorable fifth-power scaling of our original approach a two-level domain construction algorithm has been developed. First, an extended domain of localized molecular orbitals (LMOs) is assembled based on the spatial distance of the orbitals. The necessary integrals are evaluated and transformed in these domains invoking the density fitting approximation. In the second step, for each occupied LMO of the extended domain a local subspace of occupied and virtual orbitals is constructed including approximate second-order Mo̸ller-Plesset NOs. The CC equations are solved and the perturbative corrections are calculated in the local subspace for each occupied LMO using a highly-efficient CCSD(T) code, which was optimized for the typical sizes of the local subspaces. The total correlation energy is evaluated as the sum of the individual contributions. The computation time of our approach scales linearly with the system size, while its memory and disk space requirements are independent thereof. Test calculations demonstrate that currently our method is one of the most efficient local CCSD(T) approaches and can be routinely applied to molecules of up to 100 atoms with reasonable basis sets.