Archive pour mai 2008


Note in particular that GAMESS’ source can be compiled for Apple systems running OS X, or for PC systems running Linux. These are not just desktop screens, but run full-fledged Unix operating systems. Therefore, they are even capable of the parallel execution of GAMESS. However, if you do not feel comfortable about compiling source code, please look here for information about precompiled binary executables for desktop systems:  Macintosh

MacMolPlt freeware

MacMolPlt is:

  • A modern graphics program for plotting 3-D molecular structures and normal modes (vibrations). Modern means:
    • Mouse driven interface for real-time rotation and translation.
    • copy and paste functionality for interfacing to other programs such as word processors or other graphics programs (like ChemDraw).
    • simple printing to color or black and white printers (publication quality).
    • multiple files open at once.
  • It reads a variety of file formats including any GAMESS input, log or IRC file directly to create animations of IRC’s, DRC’s, and optimizations. You may also import a $VEC group from any file (such as a GAMESS .DAT file). In addition xMol XYZ files, MolDen format files and Chemical Markup Language (CML) files are supported. Also some PDB file support and MDL MolFile support is included.
  • Molecular point group symmetry is supported.
  • You may also paste GAMESS and Gaussian-92 style cartesian coordinates directly into the program.
  • Animation of Normal Modes.
  • Animation of IRC’s, and DRC’s including orbitals.
  • Simple Energy Plots (including geometrical parameters).
  • Simple frequency line graph of frequency versus infrared or Raman intensity.
  • Append multiple GAMESS files together to create a single animation.
  • Build or modify molecules using the graphical molecule builder.
    • Quickly build realistic 3D structures.
    • rotate selected atoms about bonds, change bond or dihedral angles.
    • translate and rotate selected subgroups.
  • Build molecules from scratch using cartesian or internal coordinates
  • 2D orbital, total electron density contour map display
  • 3D molecular orbital, total electron density display
  • density difference maps
  • Molecular Electrostatic Potential Maps
  • Simple GAMESS input (.inp) builder
  • 3D color display with lighting and shading using OpenGL under MacOS X.

Below are a few samples created directly with MacMolPlt:
Simple Molecule Display with labels and displaying the bond angle:

Water with the C2v operators.

Quinone with a total electron density surface colorized with the molecular electrostatic potential value.


MacOS X:

MacMolPlt 7.2.1 for MacOS X 10.3 and later on PowerPC or Intel processors (Universal binary) ( 5.8M )
MacMolPlt 7.2.1 for MacOS X 10.4 and later on PowerPC or Intel processors (Universal binary) supporting quad-buffered OpenGL ( 6.9M )
The only functional difference between these two binaries is that the latter supports hardware stereo display, but does not run on OS X 10.3.

Here is a source distribution primarily targeted at Linux: wxmacmolplt 7.2.1 ( 1.5M)
For further information on building and prebuilt binaries please refer to our Linux support page.

  • Most users should use this Windows installer (.msi) package MacMolPlt 7.2.1 for Windows ( 4.7M)
    You should be able to download it and double click it to install MacMolPlt.
  • If you can not run an installer on your system here is a plain zip archive. If the wxMacMolPlt.exe file fails to launch with a message saying the configuration is invalid you need to install the VS8 runtime libraries. Just run the file linked here.
  • The Windows build should support quad-buffered OpenGL.

These files contain the current version of MacMolPlt and the manual (in html format).

Check here for a list of known issues.

Click to subscribe to an email announcement list for new versions of MacMolPlt. Note this list is moderated and will only be used for announcing new releases.

Have a question or comment? Just include it in the comment section of the survey form. Make sure to include your email address if you want a reply.

If you need a version of MacMolPlt for MacOS 9 or MacOS X 10.2 and earlier you can grab the original MacMolPlt.


MacMolPlt is free, but if you use it I request that you send me a note (E-Mail or Snail mail) telling me what you like and what you don’t like (What’s not to like ;-).

If you use MacMolPlt for published work an appropriate reference is:

Bode, B. M. and Gordon, M. S. J. Mol. Graphics Mod., 16, 1998, 133-138.

Two-photon absorption strength



Two-photon absorption strength: A new tool for the quantification of two-photon absorption

Rémy Fortrie

Laboratoire de Chimie, UMR 5182 du CNRS, École Normale Supérieure de Lyon, 46 Allée d’Italie, F-69364 Lyon Cedex 07, France

Henry Chermettea

Université Claude Bernard Lyon 1 and UMR 5182 du CNRS Chimie Physique Théorique, Bât. Paul Dirac (210), 43 Boulevard du 11 Novembre 1918, F-69622 Villeurbanne Cedex, France

Received 21 December 2005; accepted 29 March 2006; published online 22 May 2006􏰀

In this paper, we define the two-photon absorption strength, a new characterization tool, similar to the oscillator strength, but for two-photon absorption. It allows the quantification of the two-photon absorption properties of molecular systems which are one-photon transparent. Its definition is such that the corresponding numerical values are around 100 for small molecules. We also show that this new theoretical tool allows the direct comparison of experimental and theoretical data without requiring the introduction of any arbitrary band width. As an example, the experimental and theoretical 􏰁AM1 + CNDO / S and HF + CIS / 3-21G􏰀 two-photon absorption properties of the 2 , 2􏰒-bi􏰁9,9-dihexylfluorene􏰀 molecule are compared. © 2006 American Institute of Physics. 􏰃DOI: 10.1063/1.2198530􏰄 THE JOURNAL OF CHEMICAL PHYSICS 124, 204104 􏰁2006􏰀


MOLCAS-7 is the latest version of the MOLCAS quantum chemistry software program. It has been released for distribution in the summer of 2007. MOLCAS-7 is a quantum chemistry code that allows studies of molecular systems using a variety of quantum chemical models ranging from SCF/DFT to coupled cluster and multiconfigurational SCF (RASSCF) with dynamical electron correlation treated with multi-reference CI or second order perturbation theory (MS-CASPT2). MOLCAS emphasizes the multiconfigurational approach for quantum chemical calculations, which allows studies of systems where a single configuration does not give a good representation of the electronic structure. Examples are excited states, transition states for chemical reactions, heavy element systems (transition metals, lanthanides, actinides), and much more.

Molcas is distributed as a source code with the set of configuration files for different hardware and compilers. In fact, Molcas is running under almost all Unix and Linux platforms, and also under MS Windows and Mac OS 10 operating systems.

Excited States and Electronic Spectra

MOLCAS is in particular designed to study potential surfaces for excited states.

  • Energies may be obtained using all the wave function methods. Geometry optimization is possible also for state average RASSCF energies.
  • Transition properties are computed at the RASSCF level using the RASSCF State Interaction Method, which is unique to the MOLCAS program.
  • The same code can also be used to compute spin-orbit coupling using a effective one-electron SO Hamiltonian and so called Atomic Mean Field Integrals (AMFI).
  • Automatic search at the RASSCF level for energy barriers, conical intersections, etc on excited state surfaces.
  • Vibrationally resolved electronic spectra may be obtained using the MULA code for computing transition dipole moments between harmonic vibrational levels of two electronic states.

lectures on Quantum mechanical methods

III. Semiempirical Molecular Orbital Methods

Lecture Notes

Quantum mechanical methods for the study of molecules can be divided into two categories: ab initio and semiempirical models. Ab initio methods refer to quantum chemical methods in which all the integrals are exactly evaluated in the course of a calculation. Ab initio methods include Hartree-Fock (HF) or molecular orbital (MO) theory, configuration interaction (CI) theory, perturbation theory (PT), and density functional theory (DFT). Ab initio methods that include correlation can have an accuracy comparable with experiment in structure and energy predictions. However, a drawback is that ab initio calculations are extremely demanding in computer resources, especially for large molecular systems. Semiempirical quantum chemical methods lie between ab initio and molecular mechanics (MM). Like MM, they use experimentally derived parameters to strive for accuracy; like ab initio methods, they are quantum-mechanical in nature. Semiempirical methods are computationally fast because many of the difficult integrals are neglected. The error introduced is compensated through the use of parameters. Thus, semiempirical procedures can often produce greater accuracy than ab initio calculations at a similar level.

1. Hartree-Fock Theory

We start with the stationary state Schrodinger equation:


The nonrelativistic, time independent, fixed nuclei Hamiltonian is:


where A and B designate nuclei, and I and j electrons, Z are the atomic numbers, and atomic units are used in eq (1). The many-electron wave-function Ψ is approximated as a product of one-electron functions φ(I) – the orbital approximation.


where A is the antisymmetrizer, ensuring the wavefunction obeys the Pauli exclusion principle, and O(s) is a spin projection operator that ensures that the wavefunction remains an eigenfunction of the spin-squared operator, S2.

Each molecular orbital is then expanded as a linear combination of atom orbitals, or basis functions, χu:


Utilizing the variational principle <Ψ| H |Ψ>/<Ψ|Ψ> ³ E(exp), Ψ is varied with respect to C, and an eigenvalue equation that yields the molecular orbitals and energies is obtained.


where f is an effective one-electron Fock operator with matrix elements given by


The one-electron matrix elements are give by


and two-electron integrals by


P is the first-order Fock-Dirac (or one-particle) density matrix


where ni is occupation number, 2, 1, or 0.

The Fock equations can be solved by matrix diagonalization of the matrix equations


where S is the overlap matrix and C is a square matrix with the ith column being the MO coefficients of the ith molecular orbital.

The following steps are common to all MO procedures:

1. Calculate the integrals needed to form the Fock matrix F.

2. Calculate the overlap matrix S.

3. Diagonalize S. W+SW = D

4. Form S-1/2 = WD-1/2W+.

5. Form the F matrix from equation (6).

6. Form F’ = S-1/2FS-1/2.

7. Diagonalize F’ for the MO eigenvalues E, V+F’V = E.

8. Back transform V to obtain the MO coefficients C, C = S-1/2V.

9. Form the density matrix P.

    1. Check P for convergence. If not converged, repeat steps 5-10 until self-consistent electronic field (SCF) is obtained.

2. Approximate MO theories.

Pople and coworkers introduced a series of zero-differential overlap approximations in 1965. The hierarchy of integral approximations and the effect on the Fock matrix formation are given below.

Complete Neglect of Differential Overlap (CNDO) method.

In this approximation, all integrals involving different atomic orbitals, χu, are ignored; <uv|st> = δ(u,v) δ(s,t)<uv|st> = <uu|ss>. Thus, the overlap matrix becomes the unit matrix, S = 1, and the Fock matrix elements:



Parameterization and implementation scheme of the CNDO method was also proposed by Pople:

* All two-center two-electron integrals between a pair of atoms are set equal: <uv|st> (= <uu|ss>) = γAB, where γAB is a function of atoms A and B, depending on the interatomic distance RAB.

* The off-diagonal one-electron, or resonance integrals, are set proportional to S:


where βA is a parameter that depends only on the nature of atom A.

* Electron-core attraction interactions for a given pair of atoms are set equal: <u| VB|v> = δ(u,v) VAB, with VAB = <uA| ZB/R | vA >.

* The one-center core integral Uuu is approximated by assuming that the same atomic orbital for the atom are appropriate for the positive ion: Uuu = -Iu – (ZA -1)γAA. This represents the energy required to remove an electron from atomic orbital χu in the fully ionized atom. An alternative would have been to derive Uuu from the electron affinity – the energy gained when a fully ionized atom captured an electron. CNDO/2 utilized the average of the two.

* A minimum basis set of valence orbitals was used chosen, using Slater type orbitals (STO).


where Y(θ,φ) is the real spherical harmonics.

These approximations reduce the Fock equations from the full Roothaan form to




However, it was soon realized that the electron-core attraction and the electron-electron repulsion are imbalanced in the CNDO ZDO approximation, and the appropriate two-center two-electron integral was used to approximate the electron-core attraction.


The total electronic energy was calculated from


Intermediate Neglect of Differential Overlap (INDO).

In INDO, the constraint in CNDO that monocentric two-electron integrals be equal was removed. Consequently, for atoms with s and p orbitals, there five unique two-electron one-center integrals: <ss|ss> = gss, <ss|pp> = gsp, <pp|pp> = gpp, <pp|p’p’> = gp2, p ¹ p’, and <sp|sp> = hsp. The diagonal Fock matrix element in INDO is


in which eq (15) was used in place of the electron-core attraction integral.


Since INDO and CNDO execute on a computer at about the same speeds and INDO contains some important integrals neglected in the CNDO method, INDO performs much better than CNDO especially in prediction of molecular spectral properties. One most successful program is that developed by Ridley and Zerner in 1973. The model, CNDO/S, was parameterized by carrying out CI-singles (CIS) calculations, and based one atomic spectroscopic data. In general, the INDO/S model reproduces the excitation energies of transitions below 40,000 cm-1 within 2000 cm-1. Another well-known program at this level is Dewar’s MINDO/3 model (1975). Here, all quantities that entered the Fock matrix and the energy expression were treated as free parameters. As a result, the MINDO/3 model achieved an impressive predictive power.

Neglect of Diatomic Differential Overlap (NDDO).

The NDDO model can be derived from the replacement:


This leads to the following Fock matrix elements:




Here, all two-electron two-center integrals involving charge clouds arising from pairs of orbitals on an atom were retained. This model allowed lone-pair lone-pair repulsions to be represented. Since there are four valence orbitals (for an sp atom), there are 10 unique pairs, giving rise to 100 integrals. For (spd) basis, 2025 two-electron integrals are needed compared to INDO’s 4.

The NDDO model is the only model so far that really relates to an actual basis set. The first practical NDDO model was introduced by Dewar and Thiel in 1977. Thiel has continued recently to expand the original MNDO model to include d orbitals. The model was parameterized on molecular geometries, heats of formation, dipole moments, and ionization potentials. A feature in the MNDO model is that the core-core repulsion term was made a function of the electron-electron repulsion integrals:


A major fatal problem in the MNDO model is its inability to reproduce hydrogen bonding interactions due to a spurious repulsion at just outside chemical bonding distances. The solution to this problem was to assign a number of spherical Gaussians to mimic the correlation effects. This is the Austin Model 1 (AM1):


This increased the number of parameters from the original 7 to 13-16 per atom. The AM1 model is undoubtedly a much improvement over the MNDO model for a wide range of properties. Stewart’s Parametric Method Number 3 (PM3) model treat all quantities that enter the Fock matrix as free parameters, while MNDO and AM1 derive the one-center two-electron integrals from atomic spectroscopy.

3. Integral evaluation

While in CNDO and INDO, most electron repulsion integrals were neglected, and those that are left had a common value (γAB) for a given pair of atoms. On the other hand, in NDDO, there are 22 non-vanishing integrals between first-row atoms. The NDDO electron repulsion integrals are determined in semiempirical programs in terms of multipole-multipole interactions, a method introduced by Dewar and Thiel in 1977. The atomic orbitals used in this treatment are the Slater-Zener orbitals, which are products of a radial function Rnl(r) and a normalized real spherical harmonic Ylm(θ,φ), with quantum numbers n, l, and m.


where ξ is the orbital exponent of the Slater AO, and Pl|m|(cos θ) is an associated Legendre function. Within MNDO, AM1, and PM3, only s and p basis AO were used (l = 0, 1). The multipole moments Mlm of a charge distribution ρ(r, θ, φ) is defined by


The notation of multipole moments may be compared with the common notations: M00 = q;

M10 = μz; M11 = μx; M1-1 = μy; etc.

The semiempirical formalism for the NDDO repulsion integrals were derived by first assuming that the two interacting charge distributions do not overlap. The solution to this problem is to introduce normalized real spherical harmonics and a bipolar expansion for 1/r12.

After substituting the 1/r12 expansion into the two-electron integral expression and comparison with the definition of multipoles, the two electron repulsion integrals can be expressed by


where the coefficients depends on l, m and dlm. Eq (24) has the correct asymptotic behavior for RAB ® ¥. For short distances, the assumption of zero diatomic overlap is no longer valid. Therefore, eq (23) must be modified semiempirically to ensure the behavior at short distances. The modification was provided by applying the Klopman formula for monopole-monopole interactions between charge distributions uu and vv:


The nonvanishing multipoles Mlm in the NDDO integrals are represented by configurations of 2l point charges of magnitude e/2l. The interactions between the multipoles are calculated by applying the Klopman formula to each point charge pair.


There are two values needed to be determined in this treatment, the distance separating the point charges for each multipole configuration and the parameters ρ0. The first was determined analytically by comparing the exact expression for the multipoles defined by eq (23), and the results from the point charges. These distances depend on the orbital exponents. ρ0 were chosen to yield the correct one-center limit for the monopole-monopole, dipole-dipole, and quadrupole-quadrupole interactions.

Consider, for example, the integral <spz|spz>. The only nonvanishing multipole for a charge distribution of sp is the dipole term μz. Therefore, <spz|spz> = [ μz, μz ].

In hybrid QM/MM application, the one-electron integral between a charge distribution uv in the QM region and a point charge qmm in the MM region is approximated, consistent with the MNDO/AM1/PM3 procedure, by the two-electron integrals: <u|qmm/Rmm|v> = <uv|smmsmm>. The relevant integrals, thus, include <ss|ss>, <pπ pπ|ss>, <pzpz|ss>, and <spz|ss> types. In ab initio/DFT hybrid QM/MM models, the one-electron integrals are, of course, computed analytically as do all other integrals (Freindorf, J. Comput. Chem. 17, 386 (1996)).


4. Computer programs.

The MOPAC package developed by J. J. P. Stewart was a public domain program up to MOPAC 7, and MOPAC 93. It contains the MINDO/3, MNDO, AM1, and PM3 models, which allows a variety of properties to be computed. There is also a commercial version of MOPAC as well as AMPAC. The latter contains the newest SAM1 model. The most familiar ab initio package is the Gaussian series of programs. More than a dozen of other similar programs are available. A public domain ab initio program, the GAMESS program, is also available. For spectroscopy calculations, Zerner’s INDO/S would be an excellent choice. The most sophisticated ab initio package on spectroscopy would be the MOLCAS from Roos.


Sir John Anthony Pople, FRS, (October 31, 1925March 15, 2004) was a theoretical chemist. Born in Burnham on Sea, Somerset, England, he attended Bristol Grammar School. He won a scholarship to Trinity College, Cambridge in 1943. He received his B. A. in 1946. Between 1945 and 1947 he work at the Bristol Aeroplane Company. He then returned to Cambridge University and was awarded his doctorate degree in mathematics in 1951. He moved to the United States of America in 1964, where he lived the rest of his life, though he retained British citizenship. Pople considered himself more of a mathematician than a chemist, but theoretical chemists consider him one of the most important of their number.

He received the Nobel Prize in Chemistry in 1998.[11]

Statistical mechanics of water

His early paper[3] on the statistical mechanics of water, according to Michael J. Frisch, « remained the standard for many years.[2] This was his thesis topic for his Ph D at Cambridge supervised by John Lennard-Jones.[1]

Ab Initio Electronic Structure Theory

He pioneered the development of more sophisticated computational methods, called ab initio quantum chemistry methods, that use basis sets of either Slater type orbitals or Gaussian orbitals to model the wave function. While in the early days these calculations were extremely expensive to perform, the advent of high speed microprocessors has made them much more feasible today. He was instrumental in the development of one of the most widely used computational chemistry packages, the « GAUSSIAN« (tm) suite of programs, including coauthorship of the first version, Gaussian 70.[7] One of his most important original contributions is the concept of a model chemistry whereby a method is rigorous evaluated across a range of molecules.[2] [8] He instigated the quantum chemistry composite methods such as Gaussian-1 (G1) and Gaussian-2 (G2). He was a founder of the Q-Chem computational chemistry program.[9]

ZINDO is a semi-empirical quantum chemistry method used in computational chemistry. It is a development of the INDO method. It stands for Zerner’s Intermediate Neglect of Differential Overlap, as it was developed by Michael Zerner. [1] Unlike INDO which was really restricted to organic molecules and those containing the atoms B to F, ZINDO covers a wide range of the periodic table, even including the rare earth elements. There are two distinct versions of the method:-

  • ZINDO/1 – for calculating ground state properties such as bond lengths and bond angles.
  • ZINDO/S (sometimes just called INDO/S) – for calculating excited states and hence electronic spectra.

The original program from the Zerner group is not widely available but the method is implemented in HyperChem and, in part, in Gaussian.

A major complaint against ZINDO in particular is that while it could reproduce the low lying spectra of larger polyenes and some organo-metallic compounds, it lacked a consistent parameterization. To obtain good results, it had been frequently necessary to fit the parameters to a given molecule, thereby reducing its generality and predictive capacity. In contrast, ab initio packages, such as GAUSSIAN, gained popularity because they would always produce the same results for a given input, even if the results were sometimes inaccurate.



Blog Stats

  • 223 070 hits


Flickr Photos

mai 2008