Jump to content

Z-matrix (chemistry)

From Wikipedia, the free encyclopedia
(Redirected from Internal coordinates)

In chemistry, the Z-matrix is a way to represent a system built of atoms. A Z-matrix is also known as an internal coordinate representation. It provides a description of each atom in a molecule in terms of its atomic number, bond length, bond angle, and dihedral angle, the so-called internal coordinates,[1][2] although it is not always the case that a Z-matrix will give information regarding bonding since the matrix itself is based on a series of vectors describing atomic orientations in space. However, it is convenient to write a Z-matrix in terms of bond lengths, angles, and dihedrals since this will preserve the actual bonding characteristics. The name arises because the Z-matrix assigns the second atom along the Z axis from the first atom, which is at the origin.

Z-matrices can be converted to Cartesian coordinates and back, as the structural information content is identical, the position and orientation in space, however is not meaning the Cartesian coordinates recovered will be accurate in terms of relative positions of atoms, but will not necessarily be the same as an original set of Cartesian coordinates if you convert Cartesian coordinates to a Z matrix and back again. While the transform is conceptually straightforward, algorithms of doing the conversion vary significantly in speed, numerical precision and parallelism.[1] These matter because macromolecular chains, such as polymers, proteins, and DNA, can have thousands of connected atoms and atoms consecutively distant along the chain that may be close in Cartesian space (and thus small round-off errors can accumulate to large force-field errors.) The optimally fastest and most numerically accurate algorithm for conversion from torsion-space to cartesian-space is the Natural Extension Reference Frame method.[1] Back-conversion from Cartesian to torsion angles is simple trigonometry and has no risk of cumulative errors.

They are used for creating input geometries for molecular systems in many molecular modelling and computational chemistry programs. A skillful choice of internal coordinates can make the interpretation of results straightforward. Also, since Z-matrices can contain molecular connectivity information (but do not always contain this information), quantum chemical calculations such as geometry optimization may be performed faster, because an educated guess is available for an initial Hessian matrix, and more natural internal coordinates are used rather than Cartesian coordinates. The Z-matrix representation is often preferred, because this allows symmetry to be enforced upon the molecule (or parts thereof) by setting certain angles as constant. The Z-matrix simply is a representation for placing atomic positions in a relative way with the obvious convenience that the vectors it uses easily correspond to bonds. A conceptual pitfall is to assume all bonds appear as a line in the Z-matrix which is not true. For example: in ringed molecules like benzene, a z-matrix will not include all six bonds in the ring, because all of the atoms are uniquely positioned after just 5 bonds making the 6th redundant.

Example

[edit]

The methane molecule can be described by the following Cartesian coordinates (in Ångströms):

C     0.000000     0.000000     0.000000
H     0.000000     0.000000     1.089000
H     1.026719     0.000000    -0.363000
H    -0.513360    -0.889165    -0.363000
H    -0.513360     0.889165    -0.363000

Reorienting the molecule leads to Cartesian coordinates that make the symmetry more obvious. This removes the bond length of 1.089 from the explicit parameters.

C     0.000000     0.000000     0.000000
H     0.628736     0.628736     0.628736
H    -0.628736    -0.628736     0.628736
H    -0.628736     0.628736    -0.628736
H     0.628736    -0.628736    -0.628736

The corresponding Z-matrix, which starts from the carbon atom, could look like this:

C
H   1 1.089000
H   1 1.089000  2  109.4710
H   1 1.089000  2  109.4710  3  120.0000
H   1 1.089000  2  109.4710  3 -120.0000

Only the 1.089000 value is not fixed by tetrahedral symmetry.

References

[edit]
  1. ^ a b c Parsons, Jerod; Holmes, J. Bradley; Rojas, J. Maurice; Tsai, Jerry; Strauss, Charlie E. M. (2005). "Practical conversion from torsion space to Cartesian space for in silico protein synthesis". Journal of Computational Chemistry. 26 (10): 1063–1068. CiteSeerX 10.1.1.83.8235. doi:10.1002/jcc.20237. PMID 15898109. S2CID 2279574.
  2. ^ Gordon, M. S.; Pople, J. A. (1968). "Approximate Self-Consistent Molecular-Orbital Theory. VI. INDO Calculated Equilibrium Geometries". The Journal of Chemical Physics. 49 (10): 4643–4650. Bibcode:1968JChPh..49.4643G. doi:10.1063/1.1669925.
[edit]