Molecular vibration

From Wikipedia, the free encyclopedia

A molecular vibration occurs when atoms in a molecule are in periodic motion while the molecule as a whole has constant translational and rotational motion. The frequency of the periodic motion is known as a vibration frequency, and the typical frequencies of molecular vibrations range from less than 10¹² to approximately 10¹⁴ Hz.

In general, a molecule with N atoms has 3N – 6 normal modes of vibration, but a linear molecule has 3N – 5 such modes, as rotation about its molecular axis cannot be observed.^[1] A diatomic molecule has one normal mode of vibration. The normal modes of vibration of polyatomic molecules are independent of each other but each normal mode will involve simultaneous vibrations of different parts of the molecule such as different chemical bonds.

A molecular vibration is excited when the molecule absorbs a quantum of energy, E, corresponding to the vibration's frequency, ν, according to the relation E = hν (where h is Planck's constant). A fundamental vibration is excited when one such quantum of energy is absorbed by the molecule in its ground state. When two quanta are absorbed the first overtone is excited, and so on to higher overtones.

To a first approximation, the motion in a normal vibration can be described as a kind of simple harmonic motion. In this approximation, the vibrational energy is a quadratic function (parabola) with respect to the atomic displacements and the first overtone has twice the frequency of the fundamental. In reality, vibrations are anharmonic and the first overtone has a frequency that is slightly lower than twice that of the fundamental. Excitation of the higher overtones involves progressively less and less additional energy and eventually leads to dissociation of the molecule, as the potential energy of the molecule is more like a Morse potential.

The vibrational states of a molecule can be probed in a variety of ways. The most direct way is through infrared spectroscopy, as vibrational transitions typically require an amount of energy that corresponds to the infrared region of the spectrum. Raman spectroscopy, which typically uses visible light, can also be used to measure vibration frequencies directly. The two techniques are complementary and comparison between the two can provide useful structural information such as in the case of the rule of mutual exclusion for centrosymmetric molecules.

Vibrational excitation can occur in conjunction with electronic excitation (vibronic transition), giving vibrational fine structure to electronic transitions, particularly with molecules in the gas state.

Simultaneous excitation of a vibration and rotations gives rise to vibration-rotation spectra.

Vibrational coordinates

The coordinate of a normal vibration is a combination of changes in the positions of atoms in the molecule. When the vibration is excited the coordinate changes sinusoidally with a frequency ν, the frequency of the vibration.

Internal coordinates

Internal coordinates are of the following types, illustrated with reference to the planar molecule ethylene,

• Stretching: a change in the length of a bond, such as C-H or C-C
• Bending: a change in the angle between two bonds, such as the HCH angle in a methylene group
• Rocking: a change in angle between a group of atoms, such as a methylene group and the rest of the molecule.
• Wagging: a change in angle between the plane of a group of atoms, such as a methylene group and a plane through the rest of the molecule,
• Twisting: a change in the angle between the planes of two groups of atoms, such as a change in the angle between the two methylene groups.
• Out-of-plane: a change in the angle between any one of the C-H bonds and the plane defined by the remaining atoms of the ethylene molecule. Another example is in BF₃ when the boron atom moves in and out of the plane of the three fluorine atoms.

In a rocking, wagging or twisting coordinate the bond lengths within the groups involved do not change. The angles do. Rocking is distinguished from wagging by the fact that the atoms in the group stay in the same plane.

In ethene there are 12 internal coordinates: 4 C-H stretching, 1 C-C stretching, 2 H-C-H bending, 2 CH₂ rocking, 2 CH₂ wagging, 1 twisting. Note that the H-C-C angles cannot be used as internal coordinates as the angles at each carbon atom cannot all increase at the same time.

Vibrations of a methylene group (-CH₂-) in a molecule for illustration

The atoms in a CH₂ group, commonly found in organic compounds, can vibrate in six different ways: symmetric and asymmetric stretching, scissoring, rocking, wagging and twisting as shown here:

Symmetrical stretching	Asymmetrical stretching	Scissoring (Bending)
Rocking	Wagging	Twisting

(These figures do not represent the "recoil" of the C atoms, which, though necessarily present to balance the overall movements of the molecule, are much smaller than the movements of the lighter H atoms).

Symmetry-adapted coordinates

Symmetry-adapted coordinates may be created by applying a projection operator to a set of internal coordinates.^[2] The projection operator is constructed with the aid of the character table of the molecular point group. For example, the four(un-normalised) C-H stretching coordinates of the molecule ethene are given by

where are the internal coordinates for stretching of each of the four C-H bonds.

Illustrations of symmetry-adapted coordinates for most small molecules can be found in Nakamoto.^[3]

Normal coordinates

The normal coordinates, denoted as Q, refer to the positions of atoms away from their equilibrium positions, with respect to a normal mode of vibration. Each normal mode is assigned a single normal coordinate, and so the normal coordinate refers to the "progress" along that normal mode at any given time. Formally, normal modes are determined by solving a secular determinant, and then the normal coordinates (over the normal modes) can be expressed as a summation over the cartesian coordinates (over the atom positions). The advantage of working in normal modes is that they diagonalize the matrix governing the molecular vibrations, so each normal mode is an independent molecular vibration, associated with its own spectrum of quantum mechanical states. If the molecule possesses symmetries, it will belong to a point group, and the normal modes will "transform as" an irreducible representation under that group. The normal modes can then be qualitatively determined by applying group theory and projecting the irreducible representation onto the cartesian coordinates. For example, when this treatment is applied to CO₂, it is found that the C=O stretches are not independent, but rather there is an O=C=O symmetric stretch and an O=C=O asymmetric stretch.

• symmetric stretching: the sum of the two C-O stretching coordinates; the two C-O bond lengths change by the same amount and the carbon atom is stationary. Q = q₁ + q₂
• asymmetric stretching: the difference of the two C-O stretching coordinates; one C-O bond length increases while the other decreases. Q = q₁ - q₂

When two or more normal coordinates belong to the same irreducible representation of the molecular point group (colloquially, have the same symmetry) there is "mixing" and the coefficients of the combination cannot be determined a priori. For example, in the linear molecule hydrogen cyanide, HCN, The two stretching vibrations are

1. principally C-H stretching with a little C-N stretching; Q₁ = q₁ + a q₂ (a << 1)
2. principally C-N stretching with a little C-H stretching; Q₂ = b q₁ + q₂ (b << 1)

The coefficients a and b are found by performing a full normal coordinate analysis by means of the Wilson GF method.^[4]

Newtonian mechanics

The HCl molecule as an anharmonic oscillator vibrating at energy level E₃. D₀ is dissociation energy here, r₀ bond length, U potential energy. Energy is expressed in wavenumbers. The hydrogen chloride molecule is attached to the coordinate system to show bond length changes on the curve.

Perhaps surprisingly, molecular vibrations can be treated using Newtonian mechanics to calculate the correct vibration frequencies. The basic assumption is that each vibration can be treated as though it corresponds to a spring. In the harmonic approximation the spring obeys Hooke's law: the force required to extend the spring is proportional to the extension. The proportionality constant is known as a force constant, k. The anharmonic oscillator is considered elsewhere.^[5]

By Newton’s second law of motion this force is also equal to a reduced mass, μ, times acceleration.

Since this is one and the same force the ordinary differential equation follows.

The solution to this equation of simple harmonic motion is

A is the maximum amplitude of the vibration coordinate Q. It remains to define the reduced mass, μ. In general, the reduced mass of a diatomic molecule, AB, is expressed in terms of the atomic masses, m_A and m_B, as

The use of the reduced mass ensures that the centre of mass of the molecule is not affected by the vibration. In the harmonic approximation the potential energy of the molecule is a quadratic function of the normal coordinate. It follows that the force-constant is equal to the second derivative of the potential energy.

When two or more normal vibrations have the same symmetry a full normal coordinate analysis must be performed (see GF method). The vibration frequencies,ν_i are obtained from the eigenvalues,λ_i, of the matrix product GF. G is a matrix of numbers derived from the masses of the atoms and the geometry of the molecule.^[4] F is a matrix derived from force-constant values. Details concerning the determination of the eigenvalues can be found in.^[6]

Quantum mechanics

In the harmonic approximation the potential energy is a quadratic function of the normal coordinates. Solving the Schrödinger wave equation, the energy states for each normal coordinate are given by

,

where n is a quantum number that can take values of 0, 1, 2 ... In molecular spectroscopy where several types of molecular energy are studied and several quantum numbers are used, this vibrational quantum number is often designated as v.^[7][8]

The difference in energy when n (or v) changes by 1 is therefore equal to , the product of the Planck constant and the vibration frequency derived using classical mechanics. For a transition from level n to level n+1 due to absorption of a photon, the frequency of the photon is equal to the classical vibration frequency (in the harmonic oscillator approximation).

See quantum harmonic oscillator for graphs of the first 5 wave functions, which allow certain selection rules to be formulated. For example, for a harmonic oscillator transitions are allowed only when the quantum number n changes by one,

but this does not apply to an anharmonic oscillator; the observation of overtones is only possible because vibrations are anharmonic. Another consequence of anharmonicity is that transitions such as between states n=2 and n=1 have slightly less energy than transitions between the ground state and first excited state. Such a transition gives rise to a hot band.

Intensities

In an infrared spectrum the intensity of an absorption band is proportional to the derivative of the molecular dipole moment with respect to the normal coordinate.^[9] The intensity of Raman bands depends on polarizability.

Quantum harmonic oscillator

From Wikipedia, the free encyclopedia

Some trajectories of a harmonic oscillator according to Newton's laws of classical mechanics (A–B), and according to the Schrödinger equation of quantum mechanics (C–H). In A–B, the particle (represented as a ball attached to a spring) oscillates back and forth. In C–H, some solutions to the Schrödinger Equation are shown, where the horizontal axis is position, and the vertical axis is the real part (blue) or imaginary part (red) of the wavefunction. C, D, E, F, but not G, H, are energy eigenstates. H is a coherent state—a quantum state that approximates the classical trajectory.

The quantum harmonic oscillator is the quantum-mechanical analog of the classical harmonic oscillator. Because an arbitrary potential can usually be approximated as a harmonic potential at the vicinity of a stable equilibrium point, it is one of the most important model systems in quantum mechanics. Furthermore, it is one of the few quantum-mechanical systems for which an exact, analytical solution is known.^[1][2][3]

One-dimensional harmonic oscillator

Hamiltonian and energy eigenstates

Wavefunction representations for the first eight bound eigenstates, n = 0 to 7. The horizontal axis shows the position x. Note: The graphs are not normalized, and the signs of some of the functions differ from those given in the text.

Corresponding probability densities.

The Hamiltonian of the particle is:

where m is the particle's mass, ω is the angular frequency of the oscillator, ∧x is the position operator (= x), and ∧p is the momentum operator, given by

The first term in the Hamiltonian represents the possible kinetic energy states of the particle, and the second term represents its respectively corresponding possible potential energy states.

One may write the time-independent Schrödinger equation,

where E denotes a yet-to-be-determined real number that will specify a time-independent energy level, or eigenvalue, and the solution |ψ⟩ denotes that level's energy eigenstate.

One may solve the differential equation representing this eigenvalue problem in the coordinate basis, for the wave function ⟨x|ψ⟩ = ψ(x), using a spectral method. It turns out that there is a family of solutions. In this basis, they amount to

The functions H_n are the physicists' Hermite polynomials,

The corresponding energy levels are

This energy spectrum is noteworthy for three reasons. First, the energies are quantized, meaning that only discrete energy values (integer-plus-half multiples of ħω) are possible; this is a general feature of quantum-mechanical systems when a particle is confined. Second, these discrete energy levels are equally spaced, unlike in the Bohr model of the atom, or the particle in a box. Third, the lowest achievable energy (the energy of the n = 0 state, called the ground state) is not equal to the minimum of the potential well, but ħω/2 above it; this is called zero-point energy. Because of the zero-point energy, the position and momentum of the oscillator in the ground state are not fixed (as they would be in a classical oscillator), but have a small range of variance, in accordance with the Heisenberg uncertainty principle. This zero-point energy further has important implications in quantum field theory and quantum gravity.

Note that the ground state probability density is concentrated at the origin. This means the particle spends most of its time at the bottom of the potential well, as we would expect for a state with little energy. As the energy increases, the probability density becomes concentrated at the classical "turning points", where the state's energy coincides with the potential energy. This is consistent with the classical harmonic oscillator, in which the particle spends most of its time (and is therefore most likely to be found) at the turning points, where it is the slowest. The correspondence principle is thus satisfied. Moreover, special nondispersive wave packets, with minimum uncertainty, called coherent states in fact oscillate very much like classical objects, as illustrated in the figure; they are not eigenstates of the Hamiltonian.

Ladder operator method

Probability densities |ψ_n(x)|² for the bound eigenstates, beginning with the ground state (n = 0) at the bottom and increasing in energy toward the top. The horizontal axis shows the position x, and brighter colors represent higher probability densities.

The spectral method solution, though straightforward, is rather tedious. The "ladder operator" method, developed by Paul Dirac, allows us to extract the energy eigenvalues without directly solving the differential equation. Furthermore, it is readily generalizable to more complicated problems, notably in quantum field theory. Following this approach, we define the operators a and its adjoint a^†,

This leads to the useful representation of ∧x and ∧p,

The operator a is not Hermitian, since itself and its adjoint a^† are not equal. Yet the energy eigenstates |n>, when operated on by these ladder operators, give

It is then evident that a^†, in essence, appends a single quantum of energy to the oscillator, while a removes a quantum. For this reason, they are sometimes referred to as "creation" and "annihilation" operators.
From the relations above, we can also define a number operator N, which has the following property:

The following commutators can be easily obtained by substituting the canonical commutation relation,

And the Hamilton operator can be expressed as

so the eigenstate of N is also the eigenstate of energy.

The commutation property yields

and similarly,

This means that a acts on |n⟩ to produce, up to a multiplicative constant, |n–1⟩, and a^† acts on |n⟩ to produce |n+1⟩.
For this reason, a is called a "lowering operator", and a^† a "raising operator". The two operators together are called ladder operators. In quantum field theory, a and a^† are alternatively called "annihilation" and "creation" operators because they destroy and create particles, which correspond to our quanta of energy.

Given any energy eigenstate, we can act on it with the lowering operator, a, to produce another eigenstate with ħω less energy. By repeated application of the lowering operator, it seems that we can produce energy eigenstates down to E = −∞. However, since

the smallest eigen-number is 0, and

In this case, subsequent applications of the lowering operator will just produce zero kets, instead of additional energy eigenstates. Furthermore, we have shown above that

Finally, by acting on |0⟩ with the raising operator and multiplying by suitable normalization factors, we can produce an infinite set of energy eigenstates

such that

which matches the energy spectrum given in the preceding section.

Arbitrary eigenstates can be expressed in terms of |0⟩,

Proof:

The ground state |0⟩ in the position representation is determined by a |0⟩ = 0,

and hence

and so on, as in the previous section.

Natural length and energy scales

The quantum harmonic oscillator possesses natural scales for length and energy, which can be used to simplify the problem. These can be found by nondimensionalization.

The result is that, if we measure energy in units of ħω and distance in units of √ħ/(mω), then the Hamiltonian simplifies to

while the energy eigenfunctions and eigenvalues simplify to

where H_n(x) are the Hermite polynomials.

To avoid confusion, we will not adopt these "natural units" in this article. However, they frequently come in handy when performing calculations, by bypassing clutter.

For example, the fundamental solution (propagator) of H−i∂_t, the time-dependent Schroedinger operator for this oscillator, simply boils down to the Mehler kernel,^[4][5]

where K(x,y;0) =δ(x − y). The most general solution for a given initial configuration ψ(x,0) then is simply

Phase space solutions

In the phase space formulation of quantum mechanics, solutions to the quantum harmonic oscillator in several different representations of the quasiprobability distribution can be written in closed form. The most widely used of these is for the Wigner quasiprobability distribution, which has the solution

where

and L_n are the Laguerre polynomials.

This example illustrates how the Hermite and Laguerre polynomials are linked through the Wigner map.

N-dimensional harmonic oscillator

The one-dimensional harmonic oscillator is readily generalizable to N dimensions, where N = 1, 2, 3, ... . In one dimension, the position of the particle was specified by a single coordinate, x. In N dimensions, this is replaced by N position coordinates, which we label x₁, ..., x_N. Corresponding to each position coordinate is a momentum; we label these p₁, ..., p_N. The canonical commutation relations between these operators are

The Hamiltonian for this system is

As the form of this Hamiltonian makes clear, the N-dimensional harmonic oscillator is exactly analogous to N independent one-dimensional harmonic oscillators with the same mass and spring constant. In this case, the quantities x₁, ..., x_N would refer to the positions of each of the N particles. This is a convenient property of the potential, which allows the potential energy to be separated into terms depending on one coordinate each.
This observation makes the solution straightforward. For a particular set of quantum numbers {n} the energy eigenfunctions for the N-dimensional oscillator are expressed in terms of the 1-dimensional eigenfunctions as:

In the ladder operator method, we define N sets of ladder operators,

By a procedure analogous to the one-dimensional case, we can then show that each of the a_i and a^†_i operators lower and raise the energy by ℏω respectively. The Hamiltonian is

This Hamiltonian is invariant under the dynamic symmetry group U(N) (the unitary group in N dimensions), defined by

where is an element in the defining matrix representation of U(N).

The energy levels of the system are

As in the one-dimensional case, the energy is quantized. The ground state energy is N times the one-dimensional energy, as we would expect using the analogy to N independent one-dimensional oscillators. There is one further difference: in the one-dimensional case, each energy level corresponds to a unique quantum state. In N-dimensions, except for the ground state, the energy levels are degenerate, meaning there are several states with the same energy.

The degeneracy can be calculated relatively easily. As an example, consider the 3-dimensional case: Define n = n₁ + n₂ + n₃. All states with the same n will have the same energy. For a given n, we choose a particular n₁. Then n₂ + n₃ = n − n₁. There are n − n₁ + 1 possible pairs {n₂, n₃}. n₂ can take on the values 0 to n − n₁, and for each n₂ the value of n₃ is fixed. The degree of degeneracy therefore is:

Formula for general N and n [g_n being the dimension of the symmetric irreducible n^th power representation of the unitary group U(N)]:

The special case N = 3, given above, follows directly from this general equation. This is however, only true for distinguishable particles, or one particle in N dimensions (as dimensions are distinguishable). For the case of N bosons in a one-dimension harmonic trap, the degeneracy scales as the number of ways to partition an integer n using integers less than or equal to N.

This arises due to the constraint of putting N quanta into a state ket where and , which are the same constraints as in integer partition.

Example: 3D isotropic harmonic oscillator

The Schrödinger equation of a spherically-symmetric three-dimensional harmonic oscillator can be solved explicitly by separation of variables; see this article for the present case. This procedure is analogous to the separation performed in the hydrogen-like atom problem, but with the spherically symmetric potential

where μ is the mass of the problem. (Because m will be used below for the magnetic quantum number, mass is indicated by μ, instead of m, as earlier in this article.)

The solution reads

where

is a normalization constant; ;

are generalized Laguerre polynomials; The order k of the polynomial is a non-negative integer;

is a spherical harmonic function;

ħ is the reduced Planck constant:  

The energy eigenvalue is

The energy is usually described by the single quantum number

Because k is a non-negative integer, for every even n we have ℓ = 0, 2, ...,n − 2, n and for every odd n we have ℓ =1,3,...,n − 2,n . The magnetic quantum number m is an integer satisfying −ℓ ≤ m ≤ ℓ, so for every n and ℓ there are 2ℓ + 1 different quantum states, labeled by m . Thus, the degeneracy at level n is

where the sum starts from 0 or 1, according to whether n is even or odd. This result is in accordance with the dimension formula above, and amounts to the dimensionality of a symmetric representation of SU(3), the relevant degeneracy group.

Harmonic oscillators lattice: phonons

We can extend the notion of a harmonic oscillator to a one lattice of many particles. Consider a one-dimensional quantum mechanical harmonic chain of N identical atoms. This is the simplest quantum mechanical model of a lattice, and we will see how phonons arise from it. The formalism that we will develop for this model is readily generalizable to two and three dimensions.

As in the previous section, we denote the positions of the masses by x₁,x₂,..., as measured from their equilibrium positions (i.e. x_i = 0 if the particle i is at its equilibrium position.) In two or more dimensions, the x_i are vector quantities. The Hamiltonian for this system is

where m is the (assumed uniform) mass of each atom, and x_i and p_i are the position and momentum operators for the i th atom and the sum is made over the nearest neighbors (nn). However, it is customary to rewrite the Hamiltonian in terms of the normal modes of the wavevector rather than in terms of the particle coordinates so that one can work in the more convenient Fourier space.

We introduce, then, a set of N "normal coordinates" Q_k, defined as the discrete Fourier transforms of the xs, and N "conjugate momenta" Π defined as the Fourier transforms of the ps,

The quantity k_n will turn out to be the wave number of the phonon, i.e. 2π divided by the wavelength. It takes on quantized values, because the number of atoms is finite.

This preserves the desired commutation relations in either real space or wave vector space

From the general result

it is easy to show, through elementary trigonometry, that the potential energy term is

where

The Hamiltonian may be written in wave vector space as

Note that the couplings between the position variables have been transformed away; if the Qs and Πs were hermitian(which they are not), the transformed Hamiltonian would describe N uncoupled harmonic oscillators.
The form of the quantization depends on the choice of boundary conditions; for simplicity, we impose periodic boundary conditions, defining the (N + 1)th atom as equivalent to the first atom. Physically, this corresponds to joining the chain at its ends. The resulting quantization is

The upper bound to n comes from the minimum wavelength, which is twice the lattice spacing a, as discussed above.

The harmonic oscillator eigenvalues or energy levels for the mode ω_k are

If we ignore the zero-point energy then the levels are evenly spaced at

So an exact amount of energy ħω, must be supplied to the harmonic oscillator lattice to push it to the next energy level. In comparison to the photon case when the electromagnetic field is quantised, the quantum of vibrational energy is called a phonon.

All quantum systems show wave-like and particle-like properties. The particle-like properties of the phonon are best understood using the methods of second quantization and operator techniques described later.^[6]

Applications

• The vibrations of a diatomic molecule are an example of a two-body version of the quantum harmonic oscillator. In this case, the angular frequency is given by

where μ = m₁m₂/(m₁ + m₂) is the reduced mass and is determined by the masses m₁, m₂ of the two atoms.^[7]

• The Hooke's atom is a simple model of the helium atom using the quantum harmonic oscillator
• Modelling phonons, as discussed above
• A charge, q, with mass, m, in a uniform magnetic field, B, is an example of a one-dimensional quantum harmonic oscillator: the Landau quantization.