DAMQT: A package for the analysis of electron density in molecules

Abstract DAMQT is a package for the analysis of the electron density in molecules and the fast computation of the density, density deformations, electrostatic potential and field, and Hellmann–Feynman forces. The method is based on the partition of the electron density into atomic fragments by means of a least deformation criterion. Each atomic fragment of the density is expanded in regular spherical harmonics times radial factors, which are piecewise represented in terms of analytical functions. This representation is used for the fast evaluation of the electrostatic potential and field generated by the electron density and nuclei, as well as for the computation of the Hellmann–Feynman forces on the nuclei. An analysis of the atomic and molecular deformations of the density can be also carried out, yielding a picture that connects with several concepts of the empirical structural chemistry. Program summary Program title: DAMQT1.0 Catalogue identifier: AEDL_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEDL_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GPLv3 No. of lines in distributed program, including test data, etc.: 278 356 No. of bytes in distributed program, including test data, etc.: 31 065 317 Distribution format: tar.gz Programming language: Fortran90 and C++ Computer: Any Operating system: Linux, Windows (Xp, Vista) RAM: 190 Mbytes Classification: 16.1 External routines: Trolltech's Qt (4.3 or higher) ( http://www.qtsoftware.com/products ), OpenGL (1.1 or higher) ( http://www.opengl.org/ ), GLUT 3.7 ( http://www.opengl.org/resources/libraries/glut/ ). Nature of problem: Analysis of the molecular electron density and density deformations, including fast evaluation of electrostatic potential, electric field and Hellmann–Feynman forces on nuclei. Solution method: The method of Deformed Atoms in Molecules, reported elsewhere [1], is used for partitioning the molecular electron density into atomic fragments, which are further expanded in spherical harmonics times radial factors. The partition is used for defining molecular density deformations and for the fast calculation of several properties associated to density. Restrictions: The current version is limited to 120 atoms, 2000 contracted functions, and l max = 5 in basis functions. Density must come from a LCAO calculation (any level) with spherical (not Cartesian) Gaussian functions. Unusual features: The program contains an OPEN statement to binary files (stream) in file GOPENMOL.F90. This statement has not a standard syntax in Fortran 90. Two possibilities are considered in conditional compilation: Intel's ifort and Fortran2003 standard. This latter is applied to compilers other than ifort (gfortran uses this one, for instance). Additional comments: The distribution file for this program is over 30 Mbytes and therefore is not delivered directly when download or e-mail is requested. Instead a html file giving details of how the program can be obtained is sent. Running time: Largely dependent on the system size and the module run (from fractions of a second to hours). References: [1] J. Fernandez Rico, R. Lopez, I. Ema, G. Ramirez, J. Mol. Struct. (Theochem) 727 (2005) 115.

[1]  R. Bader,et al.  Quantum Theory of Atoms in Molecules–Dalton Revisited , 1981 .

[2]  L Laaksonen,et al.  A graphics program for the analysis and display of molecular dynamics trajectories. , 1992, Journal of molecular graphics.

[3]  R. López,et al.  Analysis of the molecular density , 1999 .

[4]  A. Laaksonen,et al.  Visualization of solvation structures in liquid mixtures. , 1997, Journal of molecular graphics & modelling.

[5]  R. Bader,et al.  Quantum topology of molecular charge distributions. II. Molecular structure and its change , 1979 .

[6]  Peter Politzer,et al.  Properties of atoms in molecules. The position of the center of electronic charge of an atom in a molecule , 1971 .

[7]  F. L. Hirshfeld.  Bonded-atom fragments for describing molecular charge densities , 1977 .

[8]  R. López,et al.  Deformed atoms in molecules: analytical representation of atomic densities for Gaussian type orbitals , 2005 .

[9]  R. Stewart.  V. One‐Electron Density Functions and Many‐Centered Finite Multipole Expansions , 1977 .

[10]  Rafael López,et al.  Analytical method for the representation of atoms‐in‐molecules densities , 2004, J. Comput. Chem..

[11]  Peter Politzer,et al.  Properties of atoms in molecules: III. Atomic charges and centers of electronic charge in some heteronuclear diatomic molecules , 1971 .

[12]  Extracting atoms from molecular electron densities via integral equations. , 2004, The Journal of chemical physics.

[13]  M. Paniagua,et al.  1.1.1 electrostatic description of molecular systems , 1985 .

[14]  Á. M. Pendás,et al.  Ions in crystals: The topology of the electron density in ionic materials. I. Fundamentals , 1997 .

[15]  J. D. L. Vega,et al.  Atomic partitioning of long-range two-center potentials , 1985 .

[16]  R. Feynman.  Forces in Molecules , 1939 .

[17]  R. Parr,et al.  Regional stationary principles and regional virial theorems , 1973 .

[18]  Peter Politzer,et al.  Properties of atoms in molecules. I. Proposed definition of the charge on an atom in a molecule , 1970 .

[19]  M. Paniagua,et al.  Transferability of atomic descriptions in molecules , 1986 .

[20]  J. Maruani.  Molecules in Physics, Chemistry, and Biology , 1988 .

[21]  R. Bader.  Atoms in molecules , 1990 .

[22]  Paul L. A. Popelier,et al.  Convergence of the multipole expansion for electrostatic potentials of finite topological atoms , 2000 .

[23]  I. S. Gradshteyn,et al.  Table of Integrals, Series, and Products , 1976 .

[24]  R. Constanciel,et al.  Aspects of the Localizability of Electrons in Atoms and Molecules: Loge Theory and Related Methods , 1972 .

[25]  William T. Simpson,et al.  Coherent X‐Ray Scattering for the Hydrogen Atom in the Hydrogen Molecule , 1965 .

[26]  Rafael López,et al.  Electrostatic potentials and fields from density expansions of deformed atoms in molecules , 2004, J. Comput. Chem..

[27]  R. Daudel.  Quantum Theory of the Chemical Bond , 1974 .

[28]  P. Politzer,et al.  Properties of atoms in molecules. a proposed method for calculating the extent of distortion of an atom , 1973 .

[29]  R. López,et al.  Analysis of the molecular density: STO densities , 2002 .

[30]  Eduardo V. Ludeña,et al.  The loge theory as a starting point for variational calculations. I. General formalism , 1974 .

[31]  M. Paniagua,et al.  Atomic partitioning of two-center potentials for slater basis , 1986 .

[32]  L. Burnelle,et al.  Study of the Electron Density Distribution in the Simpler Hydrocarbons , 1962 .