Third-Body Stabilization of Supercritical CO2 in CO Oxidation: Development and Application of a ReaxFF Force Field for the CO/O/CO2 System

This paper presents the development and validation of a novel ReaxFF force field for the CO/O/CO2 system, which reveals that a supercritical CO2 matrix acts as an efficient third body to stabilize the exothermic CO oxidation reaction by dissipating excess energy through molecular collisions, thereby preventing the immediate dissociation of the newly formed CO2 product.

The Gist

The Big Picture: A Molecular Dance Floor

Imagine you are trying to build a stable house (a molecule of Carbon Dioxide, or CO₂) by snapping together two Lego bricks: a Carbon monoxide (CO) and a floating Oxygen atom (O).

In a normal, empty room (a dilute environment), when you snap these bricks together, the connection is so violent and energetic that the house immediately explodes back apart. The energy released from the "snap" is so strong that it blows the bricks apart before they can settle.

Now, imagine that same reaction happening in a crowded, packed dance floor filled with thousands of other people (a dense supercritical CO₂ environment). When you try to snap those bricks together, the crowd immediately surrounds you. They bump into you, shake your hand, and absorb the excess energy from your excitement. Because the crowd helps calm you down, your new house stays built.

This paper is about proving that the "crowd" (supercritical CO₂) acts as a safety net, stabilizing chemical reactions that would otherwise fail.

---

1. The Problem: The "Explosive" Reaction

The scientists were studying how Carbon Monoxide (a toxic gas) turns into Carbon Dioxide (a greenhouse gas, but less toxic). This happens when CO meets a highly energetic Oxygen atom.

• The Issue: This reaction releases a massive amount of energy (it's exothermic).
• The Result: In a sparse environment, the new CO₂ molecule is born with so much "hype" (kinetic energy) and "wiggling" (vibrational energy) that it literally shakes itself apart. It's like trying to build a sandcastle in the middle of a hurricane; the moment the sand touches, the wind blows it away.
• The Mystery: Scientists knew this happened, but they couldn't see how it happened because the intermediate steps happen too fast for cameras (experiments) and are too complex for standard computer models.

2. The Solution: A New "Rulebook" for Atoms

To solve this, the researchers needed a new set of rules to simulate how atoms behave.

• Old Rules (Standard Models): Think of these like a rigid instruction manual. They say, "If atoms are close, they are bonded. If they are far, they aren't." They can't handle the messy middle ground where bonds are breaking and forming.
• The New Rulebook (ReaxFF): The team created a new, flexible rulebook called ReaxFF. Imagine this as a smart, adaptive game engine. It doesn't just say "bond" or "no bond"; it calculates the strength of the connection in real-time. It can handle the chaos of atoms breaking apart and snapping back together, just like real life.

They trained this new rulebook using super-accurate quantum physics calculations (the "gold standard" of science) to make sure it was right.

3. The Experiment: Empty Room vs. Packed Dance Floor

Once they had their new rulebook, they ran two simulations:

Scenario A: The Empty Room (Dilute Gas)

• They let CO and O collide in a sparse space.
• Outcome: The CO₂ formed, but it was so energetic it immediately fell apart. The reaction was inefficient. It was like trying to light a firecracker in a vacuum; it sparks, but nothing catches.

Scenario B: The Packed Dance Floor (Supercritical CO₂)

• They added thousands of other CO₂ molecules to act as the "crowd."
• Outcome: When CO and O collided, the new CO₂ molecule was born. But immediately, the surrounding "crowd" molecules bumped into it.
• The Magic: These collisions acted like a third body (a referee) that absorbed the excess energy. The new CO₂ molecule stopped shaking, cooled down, and stayed stable.

4. The Key Findings: Where did the energy go?

The researchers looked closely at the "newly born" CO₂ molecules to see how the energy was managed. They found two fascinating things:

1. The "Third-Body" Effect: The dense fluid didn't just sit there; it actively participated. It acted like a heat sink, soaking up the heat and vibration from the new molecule through billions of tiny collisions.
2. Internal vs. External Energy: They discovered that the excess energy wasn't making the molecule fly around the room (translational energy). Instead, 92% of the energy was stored inside the molecule, making it vibrate and spin wildly (rotational and vibrational energy).
   • Analogy: Imagine a spinning top. The "crowd" doesn't stop the top from spinning in place; they just keep hitting it until it stops wobbling and spins smoothly. The energy was dissipated by making the molecule "wiggle" less, not by moving it across the room.

5. Why Does This Matter?

This isn't just about chemistry homework; it has real-world applications:

• Cleaner Energy: Supercritical CO₂ is being used in advanced power plants to generate electricity more efficiently. Understanding how these reactions work helps engineers design better, safer systems.
• Carbon Capture: It helps us understand how to trap and store carbon dioxide effectively.
• Safety: It explains why certain chemical reactions behave differently under high pressure, which is crucial for preventing explosions or designing better industrial processes.

Summary

In short, this paper developed a new, super-smart computer model to watch atoms dance. They discovered that in a crowded, high-pressure environment (supercritical CO₂), the surrounding molecules act as a safety net. They catch the "hot" new molecules, absorb their excess energy, and help them settle down, turning a reaction that would normally fail into a successful, stable process.

Technical Summary

1. Problem Statement

Supercritical carbon dioxide (scCO₂) is a critical solvent in advanced power cycles, separation processes, and materials processing. While its macroscopic thermophysical properties are well understood, the atomistic mechanisms governing how the dense scCO₂ matrix influences fundamental chemical reactions—specifically the oxidation of carbon monoxide (CO)—remain poorly explored.

• The Specific Challenge: The reaction $CO + O \rightarrow CO_2$ is highly exothermic. In a dilute environment, the newly formed $CO_2$ molecule acquires excessive kinetic and potential energy, often exceeding its bond dissociation energy (~125 kcal/mol), causing it to immediately dissociate back into $CO + O$.
• The Knowledge Gap: Experimental methods and standard Molecular Dynamics (MD) simulations struggle to detect the transient, highly reactive intermediates (like atomic oxygen) and the rapid energy transfer events required to stabilize the product. Classical non-reactive force fields cannot model bond breaking/formation, while ab initio methods are too computationally expensive for the large system sizes and time scales (nanoseconds) required to observe solvent-mediated stabilization.

2. Methodology

To bridge this gap, the authors developed a novel ReaxFF (Reactive Force Field) specifically parameterized for the $CO/O/CO_2$ system.

A. Force Field Development & Training

The ReaxFF parameters were optimized against a comprehensive training set comprising nearly 500 data points derived from:

• Quantum Mechanical (QM) Calculations: Density Functional Theory (DFT) using B3LYP/TZ2P and high-level second-order Møller-Plesset (MP2) calculations.
• Training Targets Included:
  • Crystal Properties: Bulk modulus, stiffness tensor components, and equation of state (EOS) for cubic $CO_2$ crystals.
  • Intermolecular Interactions: $CO_2$ dimer energies (slid-parallel and T-shaped) to capture van der Waals forces and long-range interactions.
  • Reaction Pathways: Energy barriers for $CO_2$ dimer formation (to prevent unphysical $C_2O_4$ formation) and the formation of transient $CO_3$.
  • Bond Dissociation: Dissociation curves for $O=O$, $O-O$, and $C \equiv O$ bonds.
  • Cohesive Energy: Solid $CO_2$ cohesive energy.
  • Bulk Properties: Energy-density EOS generated using the validated Cygan potential for larger clusters (25–50 molecules) to ensure accurate pressure/density predictions in the supercritical regime.

B. Simulation Protocols

• Validation: The force field was validated against NIST data for pressure-temperature behavior, radial distribution functions (RDF), and C-O bond length compression under high pressure.
• Reaction Dynamics: Two scenarios were simulated using NVE (microcanonical) ensembles:
  1. Dilute System: 20 $CO$+ 20$O$ atoms in a small box (low density).
  2. Dense scCO₂ System: 20 $CO$+ 20$O$ atoms embedded in 400 $CO_2$ solvent molecules at supercritical conditions (330 K, ~25 MPa, density ~0.84 g/cm³).

3. Key Contributions

1. Novel ReaxFF Parameterization: The first ReaxFF force field specifically trained to accurately model the $CO/O/CO_2$ system, capable of reproducing both bulk thermodynamic properties (EOS, RDF) and reactive chemical events (bond breaking/formation).
2. Mechanistic Insight into Third-Body Stabilization: The study provides the first atomistic demonstration of how a dense scCO₂ matrix acts as an efficient "third body" to stabilize highly exothermic reaction products.
3. Energy Dissipation Analysis: Detailed decomposition of how excess energy is transferred from the nascent product to the solvent, distinguishing between translational and internal (rotational/vibrational) modes.

4. Key Results

A. Force Field Validation

• Thermodynamics: The ReaxFF model reproduced NIST pressure data within ~5–10% across various temperatures and densities.
• Structure: Radial Distribution Functions (RDF) for liquid and scCO₂ matched well with existing non-reactive force fields (MSM, TraPPE, Cygan) and BOMD data.
• High-Pressure Behavior: The model correctly captured the qualitative trend of C-O bond shortening under compression (0–20 GPa), though with slightly lower sensitivity than MP2 reference data.

B. Reaction Dynamics: Dilute vs. Dense Environment

• Dilute Environment: The reaction $CO + O \rightarrow CO_2$ was highly inefficient. The nascent $CO_2$ accumulated excess energy (~180 kcal/mol) from the collision, exceeding the bond dissociation energy. Consequently, the molecule rapidly dissociated back into $CO + O$. No stable $CO_2$ was retained.
• Dense scCO₂ Environment: The reaction became efficient and stable. The surrounding $CO_2$ matrix acted as a heat bath.
  • Stabilization: 11 stable $CO_2$ molecules were formed and retained over a 2 ns simulation.
  • Energy Dissipation: The nascent $CO_2$ molecules dissipated an average excess energy of 133.9 ± 3.6 kcal/mol through collisions with the solvent.
  • Timescale: The average stabilization time was 112.4 ± 17.9 ps.

C. Energy Transfer Mechanism

• Internal vs. Translational: Kinetic energy decomposition revealed that ~92% of the excess kinetic energy was stored in internal degrees of freedom (rotational and vibrational modes), causing intense structural distortion (bending) of the $CO_2$ molecule. Only ~8% was in translational motion.
• Collisional Cooling: The solvent molecules absorbed this energy via repeated collisions, quenching the vibrational excitation and allowing the $CO_2$ to relax to its stable, linear equilibrium geometry.
• Transient Intermediates: The transient species $CO_3$ was observed but was extremely short-lived (~20 fs) and unstable, dissociating back to $CO_2 + O$ because the $CO_2 + O$ state is energetically more favorable.

5. Significance

• Combustion and Power Cycles: The findings are critical for optimizing direct-fired supercritical CO₂ power cycles and oxy-fuel combustion systems. They confirm that scCO₂ is not merely an inert diluent but an active participant that enhances reaction efficiency by stabilizing reactive intermediates and products.
• Methodological Advancement: This work demonstrates that ReaxFF is a viable tool for studying solvent-mediated reaction kinetics in dense fluids, overcoming the limitations of both classical non-reactive force fields (no chemistry) and ab initio MD (too small/slow).
• Future Applications: The developed force field provides a foundation for investigating impurity effects (e.g., $H_2O$, $N_2$, $O_2$) in scCO₂ systems and designing more efficient high-pressure chemical reactors.

In conclusion, the paper establishes that the dense scCO₂ environment fundamentally alters the fate of exothermic reaction products by acting as a highly efficient third body, dissipating excess energy primarily through vibrational and rotational coupling, thereby enabling the survival of species that would otherwise dissociate in dilute conditions.