Accurately reproducing ab initio electrostatic potentials with multipoles and fragmentation.

In this work, we show that our energy based fragmentation method (Bettens, R. P. A.; Lee, A. M. J. Phys. Chem. A 2006, 110, 8777) accurately reproduces the electrostatic potential for a selection of peptides, both charged and uncharged, and other molecules of biological interest at the solvent accessible surface and beyond when compared with the full ab initio or density functional theory electrostatic potential. We also consider the ability of various point charge models to reproduce the full electrostatic potential and compare the results to our fragmentation electrostatic potentials with the latter being significantly superior. We demonstrate that our fragmentation approach can be readily applied to very large systems and provide the fragmentation electrostatic potential for the neuraminidase tetramer (ca. 24,000 atom system) at the MP2/6-311(+)G(2d,p) level. We also show that by using at least distributed monopoles, dipoles, and quadrupoles at atomic sites in the fragment molecules an essentially identical electrostatic potential to that given by the fragmentation electrostatic potential at and beyond the solvent accessible surface can be obtained.

First principles NMR calculations by fragmentation.

The nuclear magnetic shielding tensor is a molecular property that can be computed from first principles. In this work we show that by utilizing the fragmentation approach, one is able to accurately compute this property for a large class of molecules. This is of great significance because the computational expense required in the evaluation of the shielding tensor for all nuclei in a large molecule is now subject to near linear scaling. On the basis of previous studies and this work, it is also very likely that all molecular properties that can be expressed as derivatives of the total energy of the system are also amenable to accurate evaluation via fragmentation. If only the chemical shifts for nuclei in a small part of a large molecule are of interest, then only those molecular fragments containing those nuclei need to have their shielding tensors evaluated. Further, the fragmentation approach allows one to construct a database of molecular fragments that could, in principle, be used in the NMR characterization of molecules and at the same time provide possible three-dimensional representations of these molecules.

A new algorithm for molecular fragmentation in quantum chemical calculations.

In this study, we present a "black-box" method for fragmenting a molecule with a well-defined Kekulé or valence-bond structure into a significant number of smaller fragment molecules that are more amenable to high level quantum chemical calculations. By taking an appropriate linear combination of the fragment energies, we show that it is possible in many cases to obtain highly accurate total energies when compared to the total energy of the full molecule. Our method is derived from the approach reported by Deev and Collins, but it contains significant unique elements, including an isodesmic approach to the fragmentation process. Using a method such as that described in this work it is in principle possible to obtain very accurate total energies of systems containing hundreds, if not thousands, of atoms as the approach is subject to massive parallelization.