QUILD: QUantum-regions interconnected by local descriptions

A new program for multilevel (QM/QM and/or QM/MM) approaches is presented that is able to combine different computational descriptions for different regions in a transparent and flexible manner. This program, designated QUILD (for QUantum-regions Interconnected by Local Descriptions), uses adapted delocalized coordinates (Int J Quantum Chem 2006, 106, 2536) for efficient geometry optimizations of equilibrium and transition-state structures, where both weak and strong coordinates may be present. The Amsterdam Density Functional (ADF) program is used for providing density functional theory and MM energies and gradients, while an interface to the ORCA program is available for including RHF, MP2, or semiempirical descriptions. The QUILD optimization setup reduces the number of geometry steps needed for the Baker test-set of 30 organic molecules by approximately 30% and for a weakly-bound test-set of 18 molecules by approximately 75% compared with the old-style optimizer in ADF, i.e., a speedup of roughly a factor four. We report two examples of using geometry optimizations with numerical gradients, for spin-orbit relativistic ZORA and for excited-state geometries. Finally, we show examples of its multilevel capabilities for a number of systems, including the multilevel boundary region of amino acid residues, an S(N)2 reaction in the gas-phase and in solvent, and a DNA duplex.