Distributed multipole analysis

Distributed multipole analysis

Distributed multipole analysis (DMA) is a compact and accurate way of describing the spatial distribution of electric charge within a molecule.

The DMA method was devised by Prof. Anthony Stone of Cambridge University to describe the charge distribution of a molecule in terms of a multipole expansion around a number of centers.[1] The idea of using a multi-center multipole expansion was earlier proposed by Robert Rein.[2] Typically, the centers correspond to the atoms constituting the molecule, though this is not a requirement. A multipole series, consisting of a charge, dipole, quadrupole and higher terms is located at each center. Importantly, the radius of convergence of this multipole series is sufficiently small that the relevant series will be convergent when describing two molecules in van der Waals contact.

The DMA series are derived from ab initio or density functional theory calculations using Gaussian basis sets. If the molecular orbitals are written as linear combinations of atomic basis functions the electron density takes the form of a sum of products of the basis functions, called density matrix elements. Boys (1950) showed that the product of two spherical Gaussian functions, centered at different points, can be expressed as a single Gaussian at an intermediate point known as the overlap center.

If a basis of Gaussian functions is used, the product of two s functions is spherically symmetric and can be represented completely just by a point charge at the ‘overlap center’ of the two Gaussian functions. The product of an s orbital and a p orbital has only charge and dipole components, and the product of two p functions has charge, dipole and quadrupole components.

If the overlap center is not at an atom, one can move the origin of the multipole expansion to the nearest distributed multipole site, re-expressing the series to account for the change of origin. The multipole expansion will no longer terminate, but the higher terms will be small. One may take the sites wherever one chooses, but they will usually be at the atoms. For small molecules one may wish to use additional sites at the centers of bonds; for larger molecules one may use a single site to describe a group of atoms such as a methyl group. The DMA procedure is exact and very fast, but for modern large basis sets with diffuse basis functions it has to be modified somewhat. When the basis functions have exponents that are small, the product function extends over several atoms, and it is better to calculate the distributed multipoles by numerical quadrature over a grid of points. The grid can be defined so that each point is associated with a particular site, and the multipoles for each site are obtained by quadrature over the points belonging to that site.

This description then includes at each site:

  • Charges, describing electronegativity effects in a chemically intuitive way;
  • Dipoles, arising from overlap of s and p orbitals and describing lone pairs and other atomic distortions;
  • Quadrupoles, arising from the overlap of p orbitals, and associated with pi bonds, for example;
  • Octopoles and hexadecapoles can be included if very high accuracy is required.

The DMA describes the potential at points outside the molecule with an accuracy which is essentially that of the wavefunction, so that its use entails no loss of precision. The DMA description gives the electrostatic energy of interaction between two molecules. It does not account for charge overlap effects and hence excludes the penetration energy.

DMA is inherently much more accurate than the commonly used partial charge methodologies for calculating intermolecular interaction energies, since it captures anisotropy of the atom-atom contributions to electrostatic interaction. It may therefore seem surprising that it has not been more widely used in molecular simulation. Possible reasons for this are:

  • Its non-inclusion in popular simulation codes;
  • The need to keep track of the orientation of a local axis system for each molecule;
  • The conformation-dependence of the DMA. As a consequence of its accuracy, the DMA captures features of the molecular charge distribution that depend strongly on molecular conformation. Thus, in a DMA-based simulation, the multipoles would have to be recalculated whenever a molecule underwent a conformational change.

DMA has found extensive use in crystal structure prediction for small organic molecules, where significant progress can often be made while using rigid molecular structures.

References

  1. ^ Stone, A.J.; Alderton, M. (1985). "Distributed multipole analysis methods and applications". Molecular Physics 56 (5): 1047–1064. doi:10.1080/00268978500102891. 
  2. ^ Rein, R. (1973). "On Physical Properties and Interactions of Polyatomic Molecules: With Application to Molecular Recognition in Biology". Advances in Quantum Chemistry 7: 335–396. doi:10.1016/S0065-3276(08)60568-7. 

Wikimedia Foundation. 2010.

Look at other dictionaries:

  • Distributed Multipole Analysis — (DMA) is a method devised by Anthony Stone to describe the charge distribution of a molecule in terms of a multipole expansion around a number of centers. References * …   Wikipedia

  • List of numerical analysis topics — This is a list of numerical analysis topics, by Wikipedia page. Contents 1 General 2 Error 3 Elementary and special functions 4 Numerical linear algebra …   Wikipedia

  • DMA — can refer to: DMA (magazine), a defunct dance music magazine Dallas Museum of Art, an art museum in Texas, USA Danish Music Awards, an award show held in Denmark since 1989 Database Marketing Agency, a financial services company in Beaverton,… …   Wikipedia

  • Crystal structure prediction — (CSP) is the calculation of the crystal structures of solids from first principles. Reliable methods of predicting the crystal structure of a compound, based only on its molecular structure, has been a goal of the physical sciences since the… …   Wikipedia

  • Cosmic microwave background radiation — CMB and Cosmic background radiation redirect here. For other uses see CMB (disambiguation) and Cosmic background (disambiguation). Physical cosmology …   Wikipedia

  • radioactivity — /ray dee oh ak tiv i tee/, n. Physics, Chem. the phenomenon, exhibited by and being a property of certain elements, of spontaneously emitting radiation resulting from changes in the nuclei of atoms of the element. Also called activity. [1895… …   Universalium

  • Gravitational wave — In physics, a gravitational wave is a fluctuation in the curvature of spacetime which propagates as a wave, traveling outward from a moving object or system of objects. Gravitational radiation is the energy transported by these waves. Important… …   Wikipedia

  • List of mathematics articles (S) — NOTOC S S duality S matrix S plane S transform S unit S.O.S. Mathematics SA subgroup Saccheri quadrilateral Sacks spiral Sacred geometry Saddle node bifurcation Saddle point Saddle surface Sadleirian Professor of Pure Mathematics Safe prime Safe… …   Wikipedia

  • Discrete element method — A discrete element method (DEM), also called a distinct element method is any of family of numerical methods for computing the motion of a large number of particles of micrometre scale size and above. Though DEM is very closely related to… …   Wikipedia

  • Boundary element method — The boundary element method is a numerical computational method of solving linear partial differential equations which have been formulated as integral equations (i.e. in boundary integral form). It can be applied in many areas of engineering and …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”