Quantum Mechanical Programs

Ab initio calculations can give a variety of molecular properties besides energy. By changing the SCAN command to OPT in the above list of input statements, the value of the HCl bond length can be optimized to give the lowest energy and the bond length can then be compared with the value determined experimentally in Exp. 37. Addition of FREQ and POLAR commands to the second line yields properties such as the harmonic vibrational frequency (vj, the dipole moment ijl, and the polarizabihty tensor elements a , etc. The calculated dipole moment can be contrasted with that determined 

In all such computations, it can be expected that the accuracy of the calculated results will generally improve with an increase in the size of the basis set and in the level of the calculation. However, this improvement often comes at considerable increase in computation time, so explorations of basis set and level effects are probably best done as class projects, with each student responsible for one or two calculations. The grid below gives a 

Another method which takes into account electron correlation is density functional theory (DFT), which is implemented in Gaussian with the label B3LYP. The basic idea of DFT is to calculate the energy of an atom or molecule in terms of electron density rather than the wavefunction tf/ The method has become increasingly popular because it can do reasonably accurate calculations on large molecules in significantly less time than that for the HF method. 

Although the emphasis in this book is on the measurement of molecular properties, the above theoretical experiment on HCl is strongly recommended to enable you to judge, by comparison with your laboratoiy results, the quality of such calculations in at least this relatively simple case. For such small gas-phase molecules, one might well question today the validity of the famous 1929 declaration of P. A. M. Dirac, Nobel laureate in physics  

The nnderlying physical laws necessaiy for the mathematical theory of a large part of physics and the whole of ehemistiy ate thus completely known, and the difficulty is only that the exact application of these laws leads to equations much too complicated to be soluble. 

---

It is always possible to convert internal to Cartesian coordinates and vice versa. However, one coordinate system is usually preferred for a given application. Internal coordinates can usefully describe the relationship between the atoms in a single molecule, but Cartesian coordinates may be more appropriate when describing a collection of discrete molecules. Internal coordinates are commonly used as input to quantum mechanics programs, whereas calculations using molecular mechanics are usually done in Cartesian coordinates. The total number of coordinates that must be specified in the internal coordinate system is six fewer... 

The traditional way to provide the nuclear coordinates to a quantum mechanical program is via a Z-matrix, in which the positions of the nuclei are defined in terms of a set of intei ii.il coordinates (see Section 1.2). Some programs also accept coordinates in Cartesian formal, which can be more convenient for large systems. It can sometimes be important to choow an appropriate set of internal coordinates, especially when locating rninima or transitinn points or when following reaction pathways. This is discussed in more detail in Section 5.7. 

It is now possible to "see" the spatial nature of molecular orbitals (10). This information has always been available in the voluminous output from quantum mechanics programs, but it can be discerned much more rapidly when presented in visual form. Chemical reactivity is often governed by the nature of the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO). Spectroscopic phenomena usually depend on the HOMO and higher energy unoccupied states, all of which can be displayed and examined in detail. 

SemiChem products are available at http //www.semichem.com/prods.html. AMPAC , available as a stand-alone product with Windows-based and work-station-level interfaces, is a semiempirical quantum mechanical program featuring SAMI, AMI, MNDO, MNDO/d, PM3, MNDO/C, and MINDO/3 semiempirical... 

You can conduct your reaction experiment in a solution-like environment. One of the options in the set-up window of a program is to let you to turn on an implicit solvation model (48). Such models attempt to account for the mutual interaction of the dipoles of the water molecules and the electron distribution in the solute molecule, as well as to account for the energy required to create a cavity in the bulk solvent to accommodate the solute. A very dilute solution is assumed, so the solute molecules do not see each other. Several implicit solvation modeling schemes have been proffered in the literature and are incorporated in various quantum mechanical programs. 

The electrostatic contribution is modelled using a Coulombic potential. The electrostatic energy is a function of the charge on the non-bonded atoms, their interatomic distance, and a molecular dielectric expression that accounts for the attenuation of electrostatic interaction by the environment (e.g. solvent or the molecule itself). Partial atomic charges can be calculated for small molecules using an ab initio or semiempirical quantum mechanics program. 

The RHF (Restricted Flartree-Fock) and B3LYP-DFT (Density Functional Theory) calculations were performed with 6-3IG basis sets using the Spartan. 02 Quantum Mechanics Program package. 

In summary, Miller s 1970 papers [1, 2] on classical 5-matrix theory had a profound influence on the theory of molecular collisions and related topics such as photodissociation. Following earlier work on elastic scattering, they demonstrated how the results of classical mechanics can be built into a quantum mechanical framework of inelastic collisions. In my view the greatest asset of the classical 5-matrix theory is its interpretative power. The general shape of transition probabilities or collisional cross sections can be easily understood in terms of classical trajectories and their quantum mechanical interference. Exact quantum mechanical programs are like black boxes and the results are often difficult to understand without the help of classical mechanics or semiclassical analyses. The new developments such as the IVR are likely to become major tools for systems consisting of many atoms. 

The availability of TM and lanthanide elements in semiempirical SCF MO methods is shown in Table 2. See Chapter 2.58 for information about different quantum mechanical program packages. Note that atomic parameters and element availability may differ in different program packages. 

The bond order is calcnlated for any pair of atoms, but has high values only if the atoms are neighbors. The smaller the valne of B for two neighboring atoms, the more ionic is the character of the chemical bond. On the other hand, the larger the value of B, the more multiple is the character (aromatic, double, triple) of the chemical bond. The calculation methods used by the quantum mechanics programs are parameterized so that they calculate the bond order B 1 for the single chemical bond, B 1.5 for the aromatic chemical bond, B 2 for the double chemical bond, and B 3 for the triple chemical bond. 

---