Wilson's Secular Equation: A Guide to Molecular Vibrations

The atoms within a molecule are engaged in a constant, intricate dance of vibration. While essential to understanding chemical properties and reactivity, describing this complex motion presents a significant challenge. A simple list of atomic positions in Cartesian coordinates is cumbersome, mixing the interesting internal vibrations with the molecule's overall translation and rotation. To truly understand the "symphony" of the molecule, a more elegant and physically meaningful framework is required. This article introduces the Wilson secular equation, a cornerstone of vibrational spectroscopy that provides exactly such a framework. It deciphers the molecular dance by focusing on internal motions like bond stretches and angle bends. The following sections will first delve into the ​​Principles and Mechanisms​​, unpacking how the potential (F) and kinetic (G) energy matrices are combined to predict vibrational frequencies. Subsequently, the ​​Applications and Interdisciplinary Connections​​ section will demonstrate how this powerful tool is used to interpret experimental spectra, analyze isotopic effects, and reveal deep connections across different areas of chemical physics.

Imagine trying to describe a dance by listing the precise $x, y, z$ coordinates of the dancer's fingertips, elbows, and knees at every instant. It would be a nightmare of numbers, and you would completely miss the elegance of the motion—the graceful bend of an arm, the powerful leap. You would see the data, but you would miss the dance. The same challenge faces us when we try to describe the dance of atoms within a molecule. Every molecule, at any temperature above absolute zero, is in a state of constant, intricate vibration. This isn't just random, chaotic jiggling; it's a symphony of highly specific, quantized motions. Our mission is to become the conductors of this molecular orchestra, to understand its score and identify its fundamental notes.

The most straightforward way to describe a molecule with $N$ atoms is to list the $3N$ Cartesian coordinates ($x_1, y_1, z_1, x_2, y_2, z_2, \ldots$) for all its atoms. This is the "list of body parts" approach. While complete, it's terribly inconvenient. These coordinates mix up the truly interesting internal motions—the vibrations—with the boring overall translation and rotation of the entire molecule through space. It's like trying to listen to a symphony while the entire concert hall is rattling down a bumpy road.

A much more natural and intuitive way is to use ​​internal coordinates​​: changes in bond lengths, bond angles, and torsional angles. For a water molecule, we wouldn't track the nine Cartesian coordinates of its atoms; instead, we'd talk about the stretching of the two O-H bonds and the bending of the H-O-H angle. This simple switch in perspective is profound. It focuses our attention on the actions that actually store and release potential energy—the very essence of vibration. By choosing a set of $3N-6$ (for non-linear molecules) of these internal coordinates, we effectively tell the molecule, "I don't care where you are in the room or how you're spinning, just show me your dance". This choice cleans up the problem immensely, but as we shall see, it introduces a delightful subtlety into how we handle the physics.

At its heart, any oscillation, from a swinging pendulum to a vibrating guitar string, is a dance between two partners: a ​​restoring force​​ trying to pull the system back to equilibrium and ​​inertia​​ causing it to overshoot. In molecular vibrations, these two partners are represented by two remarkable mathematical objects: the $\mathbf{F}$ matrix and the $\mathbf{G}$ matrix.

The $\mathbf{F}$ Matrix: The 'Stiffness' of Bonds

Imagine the atoms are connected by a network of springs. Some springs are stiff (like a strong triple bond), others are floppy (like a wide bond angle). The ​​F matrix​​, or the force-constant matrix, is the master list of all these stiffness values. Its elements, $F_{ij}$, tell us how the potential energy of the molecule changes when we distort its internal coordinates. A diagonal element $F_{ii}$ represents the restoring force for stretching a single bond or bending a single angle, while an off-diagonal element $F_{ij}$ describes how stretching one bond affects the stiffness of its neighbor.

Here we encounter a beautiful and deep principle rooted in the ​​Born-Oppenheimer approximation​​. The potential energy landscape—the "hills and valleys" the nuclei move on—is created by the fast-moving electrons. The arrangement of these electrons depends on the nuclear charges (the number of protons) and positions, but not their masses. Now, consider what happens when we replace an atom with one of its isotopes, say, hydrogen with deuterium in a water molecule. We've doubled the mass of the nucleus, but we haven't touched the number of protons or electrons. The electronic cloud is blissfully unaware of the change! Therefore, the potential energy surface, and consequently the entire $\mathbf{F}$ matrix of force constants, remains unchanged. The "springs" of the molecule have exactly the same stiffness.

The $\mathbf{G}$ Matrix: The 'Inertia' of Geometry

While the stiffness is independent of mass, the inertia most certainly is not. A heavy weight on a spring will oscillate more slowly than a light one. For molecules, however, the situation is more complex than just a collection of masses. The kinetic energy of a vibration depends not only on the atomic masses but also on the molecule's geometry.

This is where Edgar Bright Wilson Jr.'s genius shines through. He devised the ​​G matrix​​ to perfectly capture this intricate relationship between mass and geometry. Unlike the simple diagonal mass matrix in Cartesian coordinates, the $\mathbf{G}$ matrix accounts for how the motion of one internal coordinate is kinematically coupled to others. For example, in a bent molecule, bending the central angle will force the terminal atoms to move in a way that is far from simple. The $\mathbf{G}$ matrix elements are calculated from the atomic masses and the vectors describing the equilibrium geometry.

Because the $\mathbf{G}$ matrix elements are explicit functions of the atomic masses (specifically, their reciprocals, $1/m$), it must change upon isotopic substitution. Swapping H for D in water doesn't alter the $\mathbf{F}$ matrix, but it dramatically alters the $\mathbf{G}$ matrix. We have the same springs, but we've attached different weights to them. We should, therefore, expect the vibrational frequencies—the notes of the molecular song—to change. This is the theoretical origin of the isotopic shift observed in vibrational spectroscopy, a powerful tool for assigning spectral features.

So we have the potential energy information in $\mathbf{F}$ and the kinetic energy information in $\mathbf{G}$. How do we combine them to predict the molecule's dance? Wilson provided the master equation, a cornerstone of vibrational spectroscopy, known as the ​​Wilson secular equation​​:

This compact and elegant equation is a thing of beauty. It states that we must find special values, $\lambda$, for which the matrix product $\mathbf{G} \mathbf{F}$ has a non-trivial solution. This is an eigenvalue problem. The product $\mathbf{G} \mathbf{F}$ represents the full vibrational dynamics, a marriage of the kinetic and potential energy aspects of the system. The solutions, or ​​eigenvalues​​, $\lambda_k$, are directly related to the vibrational frequencies $\omega_k$ by the simple relation $\lambda_k = \omega_k^2$. Solving this single equation gives us all the fundamental frequencies of our molecular instrument! The problem of dimension $3N \times 3N$ in Cartesian coordinates, fraught with translations and rotations, is reduced to a clean $(3N-6) \times (3N-6)$ problem in internal coordinates that yields only the true vibrations.

Solving the secular equation gives us more than just the frequencies (the eigenvalues). For each eigenvalue $\lambda_k$, we also find a corresponding ​​eigenvector​​, $\mathbf{L}_k$. This vector is the "recipe" for the vibration. It tells us precisely how much each internal coordinate—each stretch and bend—participates in that specific vibration.

This eigenvector describes a ​​normal mode​​ of vibration. A normal mode is a special, synchronized motion where every atom in the molecule oscillates at the same frequency and passes through its equilibrium position at the same time. Any seemingly chaotic vibration a molecule undergoes is, in reality, just a linear combination—a chord, if you will—of these fundamental, pure tones.

These normal modes are beautifully independent. The mathematics ensures that when the molecular motion is expressed in the basis of normal coordinates, both the kinetic and potential energy matrices become diagonal. This means there are no cross-terms; energy does not flow from one normal mode to another. They are the true, orthogonal, fundamental vibrations of the system.

Let's see how this works. For a simple bent triatomic molecule like water, which has $C_{2v}$ symmetry, we can set up the $3 \times 3$ $\mathbf{G}$ and $\mathbf{F}$ matrices for the two stretches and one bend. But applying the principles of symmetry allows us to simplify the problem even before we solve it. We can combine our simple internal coordinates into ​​Symmetry-Adapted Linear Combinations (SALCs)​​. For water, this means creating a "symmetric stretch" coordinate (where both bonds stretch in unison), an "antisymmetric stretch" coordinate (where one stretches as the other compresses), and the bending coordinate.

When we transform our $\mathbf{F}$ and $\mathbf{G}$ matrices into this new, symmetry-adapted basis, a miracle happens: they become ​​block-diagonal​​. The antisymmetric stretch, belonging to the $B_2$ symmetry species, mathematically decouples entirely from the symmetric stretch and the bend, which both belong to the $A_1$ species. The $3 \times 3$ problem splits into a trivial $1 \times 1$ problem for the antisymmetric stretch and a manageable $2 \times 2$ problem for the two symmetric modes.

For a more complex molecule like a trigonal planar $XY_3$ molecule with $D_{3h}$ symmetry, this simplification is even more dramatic. A daunting problem involving many coupled stretches and bends fractures into small, independent blocks corresponding to the different symmetry species ($A_1'$, $E'$, etc.), which can be solved one by one. Symmetry is not just an aesthetic quality of molecules; it is a powerful computational tool that makes complex problems tractable.

We've solved the equations. We have our list of frequencies and the eigenvectors for each normal mode. What does a chemist do with this? The eigenvector is the key to interpreting the music. For instance, in our $2 \times 2$ problem for the symmetric modes of water, the two eigenvectors will be mixtures of the symmetric stretch and the bend. One mode will be mostly stretching with a little bit of bending, and the other will be mostly bending with a little bit of stretching.

By calculating the ​​Potential Energy Distribution (PED)​​, we can quantify this. The PED tells us what percentage of the potential energy of a given normal mode is stored in each internal coordinate. This allows us to look at an experimental spectrum, see a peak at a certain frequency, and confidently declare, "That is the C=O symmetric stretch," or "That is the C-H wagging mode." The Wilson GF method, therefore, provides not just a set of numbers, but a bridge between an abstract quantum mechanical description and the tangible, chemical reality of bonds vibrating, bending, and twisting in a beautiful, predictable, and ultimately understandable molecular dance.

Having laid the groundwork of the Wilson secular equation, we might feel we've conquered a rather abstract mathematical mountain. We have our $\mathbf{F}$ matrix of forces and our $\mathbf{G}$ matrix of motions, and we know how to combine them to solve for the vibrational frequencies of a molecule. But the real magic, the true joy of physics, is not in the equation itself, but in what it tells us about the world. It’s a key that unlocks a hidden reality, the ceaseless, beautiful dance of atoms that underlies everything. Now, let’s leave the abstract peaks and descend into the vibrant world of its applications, where this equation becomes our lens for seeing the unseen.

At its heart, the Wilson FG method is a tool for interpretation and prediction. A chemist running a spectrometer sees a spectrum, a series of peaks that look like a cryptic message from the molecular world. Our equation is the cipher. It allows us to assign each peak to a specific, well-defined motion—a normal mode—and in doing so, to understand the molecule's structure and the nature of the bonds that hold it together.

Let’s start with a molecule we all know: carbon dioxide, $\mathrm{CO}_2$. This linear, symmetric molecule seems simple enough, but its atoms are constantly stretching and vibrating. How can we make sense of it? The brute-force approach, tracking every possible jiggle, is a mess. But here, symmetry comes to our rescue. By recognizing that the two oxygen atoms are identical, we can simplify the problem immensely. The motions must be either symmetric, with both oxygen atoms moving in or out together, or asymmetric, with one moving in while the other moves out. The Wilson equation, when applied in these "symmetry coordinates," splits into two independent, much simpler problems. It tells us not just that these two modes exist, but it provides a precise relationship between their frequencies. The ratio of their squared frequencies, $(\nu_{\text{as}}/\nu_{\text{s}})^2$, is elegantly expressed in terms of the atomic masses and the bond force constants. We've turned a complex dance into a solvable two-step.

This power isn't limited to linear molecules. Consider the bent water molecule, $\mathrm{H_2O}$. An experimentalist might measure the frequency of its symmetric stretching mode, where both H atoms move in concert away from the oxygen. We can then work backwards. Knowing the frequency, the atomic masses, and the H-O-H angle, the Wilson equation allows us to calculate the "effective force constant" for that motion. This number isn't just an abstract parameter; it is a direct measure of the stiffness of the O-H bonds, a fundamental property of the chemical bond itself. We are, in essence, using the vibration as a tiny probe to measure the forces inside the molecule.

The true beauty of the method is its scalability. As molecules get more complex, our intuition might fail, but the mathematics holds strong. For a long chain like diacetylene (H−C≡C−C≡C−H), a planar molecule like boron trifluoride ($\mathrm{BF_3}$), or a highly symmetric three-dimensional cage like sulfur hexafluoride ($\mathrm{SF_6}$), the principles remain the same. Symmetry continues to be our guiding light, breaking down a dauntingly large problem into smaller, manageable blocks. For $\mathrm{SF_6}$, a molecule with 7 atoms and 15 vibrations, we can isolate the single infrared-active stretching mode ($T_{1u}$) and calculate its frequency with a simple, one-dimensional equation, as if it were a simple diatomic molecule! This is the power of a good physical theory: it finds the simplicity hidden within the complexity.

We've talked a lot about frequencies—the "tempo" of the molecular dance. But what does the dance itself look like? Each frequency, or eigenvalue, that comes out of the secular equation has a corresponding eigenvector. This eigenvector is the choreography. It describes the precise, synchronized pattern of motion of every atom in the molecule for that specific normal mode.

To get a feel for this, let's imagine a simple, hypothetical linear chain of atoms X-Y-Z, with masses $m$, $2m$, and $3m$ respectively, all connected by identical springs. Solving the equations of motion reveals two distinct stretching frequencies. For the higher frequency mode, what is the pattern of motion? The solution to the secular equation tells us more than just the frequency; it tells us the relative amplitudes of the atoms' movements. We find that the lightest atom, X, moves with a much larger amplitude than the heaviest atom, Z. For this specific mode, the ratio of their amplitudes, $A_X/A_Z$, is a precise number, $2+\sqrt{13}$. This isn't random; it's a direct consequence of the laws of mechanics. The normal mode is a collective, holistic motion of the entire system, a perfectly choreographed dance where the role of each atom is rigidly defined.

One of the most profound ideas in chemistry is the Born-Oppenheimer approximation, which states that the light, zippy electrons rearrange themselves so quickly that for the slow, heavy nuclei, the electronic potential energy landscape is essentially fixed. A direct consequence is that the force constants—the stiffness of the chemical bonds—depend on the electron configuration, not the nuclear masses.

What happens if we swap an atom for one of its heavier isotopes, for example, replacing the hydrogen atoms in water ($\mathrm{H_2O}$) with deuterium to make heavy water ($\mathrm{D_2O}$)? Chemically, nothing changes. The molecule's shape and the stiffness of its bonds remain almost perfectly identical. But the masses have changed. It's like having the same set of springs but attaching heavier weights to them. And just as you'd expect, the vibrational frequencies drop. The stretches and bends in $\mathrm{D_2O}$ are slower, a lower-pitched hum compared to $\mathrm{H_2O}$. The Wilson $\mathbf{G}$ matrix, which depends on mass, changes, while the $\mathbf{F}$ matrix does not.

This simple idea leads to a remarkably powerful and general principle: the ​​Teller-Redlich product rule​​. If we take the product of all vibrational frequencies within a given symmetry species for one isotopologue, and divide it by the same product for a second isotopologue, the force constants—all the messy details of the chemical bonds—completely cancel out! The resulting ratio depends only on the atomic masses and the molecular geometry. For a simple case like the out-of-plane vibrational mode of a planar $\mathrm{XY_3}$ molecule, the ratio of the frequency of an isotopologue $\mathrm{X'Y'_3}$ to the original is just a function of the masses $m_X, m_Y, m_{X'}, m_{Y'}$. This is an incredibly useful tool for spectroscopists. If they believe they have identified all the vibrational modes of a certain symmetry, they can check their work by seeing if the frequencies of an isotopic cousin obey this elegant rule. It's a consistency check provided by nature itself.

The influence of the Wilson equation extends beyond just identifying vibrations. It builds bridges to other areas of physics, revealing the deep unity of the principles governing molecular behavior.

A molecule in the gas phase doesn't just vibrate; it also rotates. A simple model treats it as a rigid rotor. But real molecules are not perfectly rigid. As a molecule spins faster, centrifugal force pulls the atoms apart, slightly increasing the bond lengths and thus the moment of inertia. This effect is measurable in a high-resolution rotational spectrum as a small correction term, characterized by the centrifugal distortion constant, $D_J$. Where does this constant come from? It's a battle between the centrifugal force trying to rip the molecule apart and the chemical bonds trying to hold it together. The restoring force is described by the very same force constants in our vibrational $\mathbf{F}$ matrix. It is therefore no surprise that we can derive a direct relationship between the rotational constant $D_J$ and the vibrational frequencies. This beautifully connects two seemingly different types of spectroscopy—rotational and vibrational—showing they are just different manifestations of the same underlying mechanics of masses and springs.

Finally, let's step back and admire the mathematical elegance one more time. Sometimes, we don't need to know every single frequency to test a model. The structure of the matrix equation $| \mathbf{GF} - \lambda \mathbf{I} | = 0$ offers us some wonderful shortcuts. From linear algebra, we know that the sum of the eigenvalues of a matrix is equal to its trace (the sum of its diagonal elements), and the product of the eigenvalues is equal to its determinant. Applying this, the sum of the squared frequencies ($\sum_i \omega_i^2$) for a set of vibrations is nothing more than the trace of the $\mathbf{GF}$ matrix. The product of the squared frequencies ($\prod_i \omega_i^2$) is its determinant. These relationships allow us to compute global properties of a spectrum directly from the fundamental masses and force constants, without the tedious process of finding each root. It’s another example of how the right mathematical framework doesn't just give answers; it reveals deeper structure and elegance.

From the greenhouse effect driven by the vibrations of $\mathrm{CO_2}$ and $\mathrm{H_2O}$, to the identification of complex molecules in interstellar clouds from their spectral fingerprints, the principles captured by Wilson's equation are universal. It is far more than a formula; it is a language for conversing with the molecular world, a testament to the fact that even the most complex phenomena can be understood through the harmonious interplay of fundamental laws.