Acrylamide Conformers: A Revision of Published Density Functional Theory Studies

This paper revises previous literature on acrylamide conformers by using high-precision DFT calculations to clarify that the molecule possesses three stable structures (one planar "sys/trans" and two degenerate "skew" enantiomers) rather than the previously reported two or three, providing comprehensive vibrational and structural data for these isomers.

The Gist

Imagine a molecule called Acrylamide as a tiny, flexible dancer. For years, scientists have been trying to figure out exactly how many different poses this dancer can hold without falling over (becoming unstable).

Here is the story of a new study that cleared up a confusing mix-up about this dancer's poses.

The Confusing Backstory

For a long time, the scientific community couldn't agree on the number of stable poses Acrylamide could strike.

• Some researchers said there were two poses.
• Others said there were three.
• A giant online database (PubChem) listed four different poses, but it was messy. It included some poses that were actually just "mid-dance" movements (unstable) rather than final poses.

It was like a group of photographers trying to describe a gymnast, but some were counting the flips as final poses, while others missed a pose entirely.

The New Investigation

The authors of this paper, William Scott and Estela Blaisten-Barojas, decided to settle the score. They used a super-powerful computer simulation method called Density Functional Theory (DFT). Think of this as a "molecular microscope" that can see the energy of every single atom with extreme precision.

They took the four poses listed in the database and ran them through their high-precision simulation to see what actually happens.

The Big Discovery: Three Stable Poses

After running the numbers, they found that there are actually three distinct, stable poses (conformers) that the molecule can hold. They named them S1, S2, and S3.

Here is what they found, using a simple analogy:

1. The "Flat" Pose (S1) – The Most Stable

• What it is: This is the molecule lying completely flat, like a pancake on a table.
• Status: This is the "champion." It has the lowest energy, meaning it is the most comfortable and stable position. The molecule prefers to stay here.
• Literature name: Scientists have called this syn, trans, or conformer 1.

2. The "Twisted Mirror" Poses (S2 and S3) – The Runners-Up

• What they are: These are two poses where the molecule twists slightly out of the flat plane, like a person leaning to the side.
• The Twist: S2 and S3 are mirror images of each other. Imagine your left hand and right hand; they look the same but are flipped. They are "energy-degenerate," meaning they have the exact same energy level. They are slightly less stable than the flat pancake (S1), but they are still very stable.
• Literature name: These are the ones often called skew.

The "Fake" Pose and the Bridges

The study also clarified what happened to the other poses listed in the database:

• The "Fake" Stable Pose (Conformer 2): The database listed a second flat pose (called cis or anti) as a stable structure. The new study proved this is wrong. When you try to build this shape, it doesn't stay put; it immediately slides down into a transition state. It's not a resting pose; it's a bridge (a transition state) between the two twisted mirror poses (S2 and S3).
• The Bridges (Transition States): To get from one stable pose to another, the molecule has to climb a small energy hill. The study mapped out these three "hills" (T12, T13, and T23).
  • Going from the Flat Pose (S1) to a Twisted Pose (S2 or S3) requires climbing one hill.
  • Going from one Twisted Pose (S2) to its mirror twin (S3) requires climbing a different, flatter hill.

Why This Matters

Before this study, scientists were arguing about whether the "twisted" pose was one thing or two, and whether the "second flat" pose was real.

This paper acts like a referee with a slow-motion camera. It confirms:

1. There are three real, stable resting spots for the molecule.
2. The "second flat" spot is actually just a bridge between the two twisted spots.
3. The two twisted spots are mirror images, not the same thing.

The Takeaway

The authors didn't just guess; they provided a complete "map" of the molecule's energy landscape. They gave the exact coordinates (where every atom sits), the "vibrational fingerprint" (how the molecule shakes, which helps identify it in a lab), and the energy costs to move between these poses.

In short: Acrylamide has three stable homes, not two or four. One is a flat pancake, and the other two are twisted mirrors of each other. The "fourth" option people thought existed is actually just a bridge between the mirrors.

Technical Summary

1. Problem Statement

The study addresses a significant discrepancy in the scientific literature regarding the number and nature of stable conformers of acrylamide ($CH_2=CH-CONH_2$).

• Literature Conflict: Previous studies using various Density Functional Theory (DFT) methods (e.g., B3LYP/6-311+G**, B3LYP/cc-pVTZ) have reported conflicting results, identifying either two or three stable conformers. Some studies reported a planar "syn/trans" structure and a 3D "skew" structure, while others reported planar "trans" and "cis" structures.
• Database Ambiguity: The PubChem database (CID 6579) lists four distinct conformer geometries: one planar (syn/trans), a second planar structure (potentially cis/anti), and two 3D mirrored structures with the amide group rotated 90° out of plane.
• Objective: The authors aim to resolve these contradictions by performing high-precision DFT calculations to definitively determine the number of stable energy minima (conformers), the transition states (energy barriers) between them, and their intrinsic reaction coordinates (IRC).

2. Methodology

The authors employed a rigorous computational approach to re-evaluate the potential energy surface (PES) of acrylamide.

• Primary Method: High-precision DFT calculations using the $\omega$B97XD functional with the Def2TZVPP basis set. This method includes dispersion corrections, which are critical for accurately modeling non-covalent interactions and conformational stability.
• Comparative Method: Calculations were also performed using the B3LYP/6-311+G** method to compare results with historical literature.
• Software & Settings: The Gaussian 16 package was used with strict optimization criteria (opt=verytight, Int=Ultrafine, scf=QC) to ensure convergence to true minima.
• Analysis Tools:
  • VMD (Visual Molecular Dynamics) was used to calculate Root Mean Square Deviation (RMSD) to quantify structural differences.
  • QST3 algorithm was used to locate Transition State (TS) structures.
  • Intrinsic Reaction Coordinate (IRC) calculations were performed to trace the path from energy minima, through transition states, to adjacent minima.
  • Spectral Analysis: Vibrational frequencies were calculated and scaled (0.9552 for $\omega$B97XD, 0.967 for B3LYP) to generate IR spectra. Partial atomic charges were derived using Mulliken and Electrostatic Potential (ESP) analyses.

3. Key Contributions

The paper provides a definitive revision of acrylamide's conformational landscape:

1. Clarification of Stable Conformers: The study confirms the existence of three distinct stable conformers (energy minima), not two or four as previously debated.
2. Identification of Transition States: The authors identified and characterized three first-order saddle points (transition states) connecting these minima.
3. Resolution of PubChem Discrepancy: The study clarifies that the "second planar structure" (Conformer 2 in PubChem) is not a stable conformer but a transition state. Conversely, the two 3D mirrored structures (Conformers 3 and 4 in PubChem) are confirmed as stable, energy-degenerate isomers.

4. Key Results

A. Stable Conformers (Minima)

Three stable isomers were identified, all in singlet ground states:

• S1 (Lowest Energy): A planar structure with a dihedral angle ($\theta$) of $180^\circ$ between the C=C and C=O bonds. This corresponds to the literature terms "syn" or "trans."
  • Energy: Reference point (0.0 kcal/mol).
  • Dipole Moment: 3.6512 D (in-plane).
• S2 and S3 (Higher Energy): Two 3-dimensional, energy-degenerate structures that are mirror images of each other. They possess a dihedral angle $\theta \approx 26.15^\circ - 28.15^\circ$.
  • Energy: $\sim$1.22 kcal/mol ($\omega$B97XD) or 1.41 kcal/mol (B3LYP) above S1.
  • Dipole Moment: 3.9083 D (out-of-plane, opposite directions for S2 and S3).
  • Significance: These correspond to the "skew" structure mentioned in earlier literature.

B. Transition States (Saddle Points)

• T12 and T13: Transition states connecting the planar S1 to the 3D S2 and S3, respectively. These are mirrored structures with an energy barrier of $\sim$4.32 kcal/mol.
• T23: A planar transition state connecting the two 3D isomers (S2 $\leftrightarrow$ S3). This structure corresponds to the PubChem "Conformer 2" and the literature "cis/anti" structure, which was previously misidentified as a stable conformer.
  • Energy Barrier: $\sim$1.61 kcal/mol (relative to S2/S3).

C. Spectral and Structural Validation

• Vibrational Spectra: The IR spectra of S1 are distinct from the identical spectra of S2 and S3. The authors provide scaled frequencies, confirming that experimental IR spectroscopy could easily distinguish between the planar and 3D conformers.
• Structural Differences: RMSD analysis confirms that S2 and S3 are distinct 3D geometries despite having identical energies.
• Triplet States: The triplet states for all conformers are significantly higher in energy (>90 kcal/mol), confirming the ground state is a singlet.

5. Significance

This revision is critical for the accuracy of computational chemistry databases and future molecular modeling involving acrylamide.

• Correcting the Literature: It resolves the confusion between "cis/anti" and "skew" by proving that the planar "cis/anti" geometry is actually a transition state, not a stable minimum.
• Database Integrity: It corrects the interpretation of PubChem data, clarifying that the database contains two stable 3D isomers and one stable planar isomer, with the "fourth" geometry being a transition state.
• Experimental Guidance: By providing precise Cartesian coordinates, vibrational frequencies, and dipole moments, the study offers a reliable benchmark for experimentalists attempting to characterize acrylamide conformers via IR spectroscopy or microwave spectroscopy.
• Methodological Validation: The study demonstrates the necessity of using dispersion-corrected functionals (like $\omega$B97XD) and high-precision settings to correctly identify subtle energy differences between planar and non-planar conformers.