Photo-induced electronic excitations drive polymerization of carbon monoxide: A first-principles study

Basic Information

• Paper ID: 2505.16926
• Title: Photo-induced electronic excitations drive polymerization of carbon monoxide: A first-principles study
• Authors: Rasool Ahmad, Jonathan C. Crowhurst, Stanimir A. Bonev (Lawrence Livermore National Laboratory)
• Classification: cond-mat.mtrl-sci (Condensed Matter Physics - Materials Science)
• Publication Date: May 22, 2025 (arXiv preprint)
• Paper Link: https://arxiv.org/abs/2505.16926

Abstract

Carbon monoxide (CO) can be converted to a polymeric phase under high pressure, and this polymer remains stable at ambient pressure. Although this transformation can be achieved through compression alone, experimental observations demonstrate that laser irradiation can induce similar transformations at significantly lower pressures. Using time-dependent density functional theory (TDDFT) and Born-Oppenheimer molecular dynamics (BOMD) simulations, this study reveals the mechanism by which electronic excitations promote CO polymerization. Computational results demonstrate that electronic excitations enhance carbon-carbon bonding, enabling polymerization reactions to occur at pressures substantially lower than required by conventional compression methods. These findings provide a photo-assisted synthesis pathway for synthesizing novel high-energy materials under milder conditions.

Research Background and Motivation

Problem Statement

1. Core Issue: Carbon monoxide can be converted to a polymer under high pressure (~5-7 GPa), but this requires extreme pressure conditions. Experiments show that laser irradiation can reduce the transition pressure to 3.2 GPa, yet the microscopic mechanism remains unclear.
2. Significance:
   • CO polymers are metastable high-energy-density materials capable of releasing energy through decomposition reactions (2CO → CO₂ + C)
   • Understanding the photo-induced polymerization mechanism facilitates development of novel energetic material synthesis methods
   • Reducing synthesis pressure makes material preparation more practical
3. Limitations of Existing Research:
   • Extensive experimental characterization of polymerized CO structure and energetic properties exists
   • Theoretical simulations have primarily focused on phase transitions under pure compression
   • Critical Gap: Systematic first-principles investigation of the microscopic mechanism of laser-induced photochemical transformation is lacking
4. Research Motivation:
   • Evans et al. experimentally demonstrated that 488 nm and 514 nm laser radiation can reduce the transition pressure from 5.2 GPa to 3.2 GPa
   • Theoretical clarification is needed on how electronic excitations alter chemical bonding and phase transition dynamics
   • Theoretical guidance is required for photo-assisted material synthesis

Core Contributions

1. First Systematic Study: Systematic investigation using first-principles methods of the mechanism by which electronic excitations affect the CO molecular-polymer transition
2. Multi-scale Theoretical Approach: Integration of three complementary computational methods:
   • FPMD simulations reveal macroscopic phase transition behavior
   • TDDFT/TDA calculations elucidate molecular-level bonding mechanisms
   • LR-TDDFT determines optical absorption spectra and pressure effects
3. Key Findings:
   • Electronic excitations significantly reduce the polymerization transition pressure (from ~10 GPa to ~5 GPa)
   • Excited states promote carbon-carbon bond formation and suppress CO molecular stability
   • Compression shifts the optical absorption band gap from 4.85 eV to 3.51 eV, facilitating visible light excitation
4. Theory-Experiment Comparison: Qualitatively explains the experimental observations of Evans et al.:
   • Reduced transition pressure
   • Progressive nature of photo-assisted transition with broader coexistence range of molecular and polymeric phases

Methodology Details

Task Definition

Research Objective: Clarify how electronic excitations affect the CO molecular → polymeric phase transition, specifically including:

• Input: CO systems at various electronic excitation levels, pressure/density conditions
• Output: Phase transition pressure, structural characteristics, bonding mechanisms, optical absorption spectra
• Constraints: Constant temperature at 300 K, quasi-static compression process

Methodological Framework

This study employs a three-level computational strategy:

1. FPMD Simulation (Macroscopic Dynamics)

Implementation Details:

• Software: VASP
• System Scale: 64 CO molecules, periodic boundary conditions
• Exchange-Correlation Functional: PBE-GGA
• Pseudopotentials: PAW method, core radii of 1.50 and 1.52 Bohr for C and O respectively
• Plane Wave Cutoff Energy: 550 eV
• k-point Sampling: Single Γ point
• Ensemble: NVT (Langevin thermostat, T = 300 K)
• Time Step: 1 fs
• Compression Protocol: Starting from 1.12 g/cm³, volume reduced by 1% per step, 0.1 ps equilibration per step

Electronic Excitation Simulation:

• Innovative Approach: Artificially transfer electrons from highest occupied states to lowest unoccupied states
• Three Excitation Levels:
  • Level 1: Excitation energy 2.9 eV
  • Level 2: Excitation energy 2.8 eV
  • Level 3: Excitation energy 8.0 eV
• Physical Validity: Excitation energies correspond to visible light range, achievable via laser

2. TDDFT/TDA Calculation (Molecular Bonding Mechanism)

Implementation Details:

• Software: ORCA
• System: Two CO molecules (in vacuum)
• Configuration: C-C relative orientation, O-O back-to-back
• Functional: PBE-GGA
• Basis Set: def2-tzvpp all-electron basis set
• Computational Content: Ground state and first 9 excited state energies as function of C-C distance (1.45-2.0 Å)
• Purpose: Direct verification whether excited states promote C-C bond formation

3. LR-TDDFT Calculation (Optical Response)

Implementation Details:

• Software: SALMON (real-space TDDFT code)
• System: 64 CO molecules
• Density: 1.12 g/cm³ and 1.70 g/cm³ (pre-polymerization)
• Grid: 80×80×80 real-space grid
• Functional: PBE-GGA (adiabatic approximation)
• Time Step: 10⁻³ fs, total duration 20 fs
• Method: Apply pulsed electric field, time-evolve Kohn-Sham equations, Fourier transform to obtain frequency-dependent dielectric function

Technical Innovations

1. Novel Excited-State FPMD Scheme:
   • Traditional BOMD considers only ground state
   • This work simulates sustained excitation effects through occupation number control
   • While not strictly excited-state dynamics, effectively captures qualitative trends
2. Multi-method Cross-Validation:
   • FPMD: Statistical behavior and phase transition thermodynamics
   • TDDFT/TDA: Precise excited-state energy surfaces
   • LR-TDDFT: Actual optical absorption characteristics
   • Three methods complement each other, enhancing conclusion reliability
3. Quantitative-Qualitative Integration with Experiment:
   • Acknowledges PBE-GGA overestimation of transition pressure (~10 GPa vs experimental 5-7 GPa)
   • References Bonev et al. work showing HSE hybrid functional can improve results
   • Emphasizes correctness of qualitative trends

Experimental Setup

Computational Parameters

Evaluation Metrics

1. Structural Indicators:
   • Radial distribution function (RDF): C-C, C-O, O-O pairs
   • Evolution of chain/molecule count
   • Ring structure statistics (4-membered, 5-membered rings)
2. Thermodynamic Indicators:
   • Pressure-density relationship curves
   • Transition pressure (P_trans) definition: pressure drop point
3. Electronic Structure Indicators:
   • Excited-state energy curves
   • Optical absorption coefficient
   • Band gap variation with pressure

Comparison Schemes

• Baseline: Ground-state molecular dynamics (GSMD), no electronic excitation
• Experimental Comparison: Three excitation levels
• Density Variable: Continuous compression from 1.12 → 2.5 g/cm³

Experimental Results

Main Results

1. Significant Reduction in Transition Pressure (Figure 2a)

Key Observations:

• Electronic excitations reduce transition pressure by approximately 50%
• Consistent with experimental trends (5.2 GPa → 3.2 GPa)

2. Transition Mechanism Differences (Figure 2b)

• GSMD: Sharp drop in chain count at 1.8 g/cm³, sharp transition
• Excited States:
  • C-C bond formation begins at initial density (1.12 g/cm³)
  • Smooth transition with broad coexistence of molecular and polymeric phases
  • Chain count begins decreasing from <64 (partial pre-polymerization)

3. Structural Characteristics (Figure 3 RDF Analysis)

GSMD at 1.12 g/cm³:

• C-O peak: 1.12 Å (pure molecular phase)
• C-C: No peak within 2.5 Å

Excited States at 1.12 g/cm³:

• C-C peak: Already appears at ~1.5 Å (C-C bond formation)
• C-O peak: Slightly weakened

All Cases at 2.5 g/cm³:

• C-C peak: Significantly enhanced at ~1.5 Å
• C-O peak: Markedly weakened
• Polymerization degree: Excited states > GSMD

4. Ring Structure Formation (Figures 2c-d)

• GSMD: 4- and 5-membered rings appear after 1.8 g/cm³
• Excited States: Ring formation begins at 1.3-1.5 g/cm³
• Physical Significance: Polymer is not purely linear, contains complex topology

TDDFT/TDA Results (Figure 4)

Ground State Energy Curve:

• Decreasing C-C distance → Energy monotonically increases
• Reflects intrinsic CO molecular stability

First Excited State Energy Curve:

• Decreasing C-C distance → Energy decreases (1.45-2.0 Å range)
• Key Finding: Excited state favors C-C bond formation!

Higher-Order Excited States:

• 4th excited state: Energy minimum at ~1.5 Å
• 6th excited state: Energy minimum at ~1.8 Å
• Confirms multiple excited states promote C-C bonding

Theoretical Significance:

• Directly proves excited states alter potential energy surface
• Explains enhanced bonding observed in FPMD

Optical Absorption Spectrum Results (Figure 5)

Key Findings:

1. Compression Reduces Band Gap: 4.85 eV → 3.51 eV (27% reduction)
2. Absorption Threshold Red-Shift: Lower energy photons can be absorbed at high pressure
3. Visible Light Feasibility: At 1.70 g/cm³, approaching visible light range (1.6-3.1 eV)

Practical Implications:

• Initially requires UV light excitation
• After pre-compression, visible light becomes viable (e.g., experimental 488/514 nm)
• Synergistic effect of pressure and light

Ablation Studies

While not traditional ablation, systematic investigation of excitation level effects through three levels:

Trend: Higher excitation energy → Lower transition pressure → More pronounced pre-polymerization

Experimental Findings Summary

1. Qualitative Consistency:
   • Reduced transition pressure ✓
   • Smoothed transition ✓
   • Molecular-polymer coexistence ✓
2. Quantitative Discrepancy:
   • Calculated transition pressure higher (PBE-GGA limitation)
   • Relative pressure reduction magnitude approaches experiment
3. New Predictions:
   • Compression causes optical absorption red-shift
   • Multiple excited states promote polymerization
   • Pre-compression + visible light potentially more efficient

Related Work

Experimental Research Timeline

1. Early Discovery (1980s):
   • Katz et al. (1984): First report of high-pressure CO polymerization
   • Mills et al. (1985): X-ray characterization of phase diagram
2. Photochemical Transformation (1990s-2000s):
   • Lipp et al. (1998, 2005): Spectroscopic characterization of polymeric phase
   • Evans et al. (2006): Key experiment discovering laser-reduced transition pressure
   • Ceppatelli et al. (2009): FTIR study of pressure-induced reactivity
3. High-Pressure Extension (2010s):
   • Ryu et al. (2016): Progressive polymerization at 160 GPa

Theoretical Research Progress

1. Structure Prediction:
   • Bernard et al. (1998): Decomposition and polymerization mechanisms
   • Sun et al. (2011): Pressure-controlled bonding and band gap
   • Xia et al. (2017): Ground-state structure prediction
2. Energetics:
   • Bonev et al. (2021): FPMD simulation with good experimental agreement
   • Authors' prior work validating PBE-GGA applicability
3. This Work's Position:
   • First systematic investigation of photo-induced mechanism
   • Fills gap between experiment (Evans 2006) and theory
   • Multi-method cross-validation

Methodological Comparison

Conclusions and Discussion

Main Conclusions

1. Mechanism Clarification: Electronic excitations promote C-C bond formation and suppress CO molecular stability by altering the potential energy surface
2. Quantitative Prediction: Excited states can reduce transition pressure from ~10 GPa to ~5 GPa (50% reduction)
3. Optical Characteristics: Compression shifts band gap from 4.85 eV to 3.51 eV, facilitating visible light excitation
4. Experimental Validation: Qualitatively explains Evans et al. experimental observations, including reduced transition pressure and phase coexistence phenomena
5. Universal Significance: Photo-assisted synthesis strategy may be applicable to other energetic materials (e.g., high-pressure nitrogen)

Limitations

Methodological Level

1. PBE-GGA Functional:
   • Overestimates transition pressure (~10 GPa vs experimental 5-7 GPa)
   • Limited accuracy for optical absorption spectra
   • Authors recommend HSE hybrid functional for improvement
2. Excited-State Treatment:
   • Artificial occupation number control is not rigorous excited-state dynamics
   • Does not account for excited-state lifetime and relaxation
   • Cannot directly correlate with laser parameters (intensity, wavelength)
3. System Size:
   • TDDFT/TDA only for 2 CO molecules
   • May miss many-body effects

Physical Assumptions

1. Quasi-Static Process: 0.1 ps equilibration may be insufficient
2. Single Temperature: Only 300 K considered
3. Pure CO System: Actual samples may contain impurities

Future Directions

Authors explicitly propose:

1. Ehrenfest Dynamics:
   • Explicit coupling of ionic and excited electron dynamics
   • Direct simulation of laser-matter interaction
   • Requires further validation for large systems
2. Hybrid Functionals:
   • HSE and others to improve quantitative accuracy
   • More accurate band gaps and optical properties
3. Other Materials:
   • Extension to nitrogen, CO-N₂ mixtures, etc.
   • Explore universality of photo-assisted synthesis
4. Experimental Verification:
   • Systematic experiments with different laser wavelengths
   • Time-resolved spectroscopy tracking dynamics

In-Depth Evaluation

Strengths

Scientific Contribution

1. Fills Important Gap: First theoretical clarification of photo-induced CO polymerization mechanism, addressing theoretical void since Evans 2006 experiment
2. Methodological Innovation:
   • Practical excited-state FPMD scheme
   • Three complementary methods with mutual verification, enhancing credibility
   • Balances computational cost with physical insight
3. Experimental Relevance:
   • Qualitative predictions highly consistent with experiment
   • Quantitative discrepancies have reasonable explanations
   • Provides testable new predictions (compression-optical absorption relationship)

Technical Rigor

1. Reasonable Parameter Selection:
   • 64-molecule system sufficient for statistical sampling
   • 550 eV cutoff energy adequately converged
   • Three excitation levels systematically explore parameter space
2. Comprehensive Result Analysis:
   • Multi-angle characterization via RDF, chain count, ring structures
   • System evolution at different densities
   • Detailed excited-state energy surface analysis
3. Honest Discussion of Limitations:
   • Clearly identifies PBE-GGA shortcomings
   • Acknowledges simplifications in excited-state treatment
   • References prior work validating method reliability

Presentation Quality

1. Clear Figures: Figures 1-5 logically coherent, information-rich
2. Fluent Writing: Reasonable structure, sufficient technical detail
3. Comprehensive Literature Review: Covers experimental and theoretical progress

Weaknesses

Methodological Limitations

1. Simplified Excited-State Dynamics:
   • Artificial occupation numbers cannot reflect real excited-state lifetime
   • Neglects non-adiabatic coupling and excited-state relaxation
   • Impact: Limited reliability of quantitative predictions
2. Single Functional:
   • Only PBE-GGA employed
   • No comparison with HSE or others
   • Impact: Systematic overestimation of transition pressure
3. Small System TDDFT/TDA:
   • 2 molecules may be insufficient to represent condensed phase
   • Lacks solvation/periodicity effects
   • Impact: Universality of bonding mechanism needs verification

Experimental Verification

1. Lack of Direct Comparison:
   • No Raman spectrum calculation for experimental comparison
   • No simulation of specific experimental conditions (488 nm, 3.2 GPa)
   • Recommendation: Supplement with spectral calculations to enhance persuasiveness
2. Missing Kinetic Information:
   • No polymerization reaction rate provided
   • No time-scale estimates
   • Impact: Difficult to guide experimental optimization

Analysis Depth

1. Electronic Structure Details:
   • No analysis of which specific orbitals participate in excitation
   • Lacks charge density/orbital diagrams
   • Recommendation: Supplement with orbital analysis for enhanced physical picture
2. Temperature Effects:
   • Only single temperature (300 K)
   • No exploration of temperature-pressure-light synergy
   • Impact: Limited application condition range

Impact Assessment

Academic Value

1. Theoretical Breakthrough: ★★★★☆
   • First systematic theoretical investigation
   • Clear mechanism elucidation
   • Method still has improvement potential
2. Methodological Contribution: ★★★☆☆
   • Practical but non-rigorous excited-state FPMD scheme
   • Multi-method combination strategy worth emulating
   • Ehrenfest dynamics outlook important
3. Experimental Guidance: ★★★★☆
   • Testable predictions of compression-optical absorption relationship
   • Provides theoretical basis for experimental parameter optimization
   • Quantitative accuracy needs improvement

Practical Value

1. Material Synthesis:
   • Reduced synthesis pressure has industrial significance
   • Photo-assisted strategy extensible to other materials
   • Potential Applications: High-energy materials, extreme-condition chemistry
2. Technical Feasibility:
   • Visible light + moderate pressure (<5 GPa) achievable
   • More practical than pure compression method (>10 GPa)
   • Requires further engineering optimization

Reproducibility

1. Computational Reproducibility: ★★★★☆
   • Detailed parameter description
   • Open-source (SALMON) or commercial (VASP, ORCA) software
   • Excited-state occupation scheme requires VASP scripting
2. Experimental Verifiability: ★★★★☆
   • Clear predictions (transition pressure, optical absorption)
   • Requires diamond anvil cell + laser equipment
   • High cost

Applicable Scenarios

Direct Applications

1. CO Polymer Synthesis:
   • Laboratory-scale synthesis
   • Parameter optimization (pressure, laser wavelength)
2. High-Energy Material Research:
   • Model system for understanding photo-assisted polymerization
   • Guides similar material (N₂, CO-N₂ mixture) research

Method Transfer

1. Photo-Chemical High-Pressure Research:
   • Excited-state FPMD scheme applicable to other systems
   • Multi-method combination strategy universal
2. Energetic Material Design:
   • Universal mechanism of electronic excitation enhancing bonding
   • Provides insights for new material screening

Theoretical Extension

1. Excited-State Dynamics:
   • Advances Ehrenfest method development
   • Large-system non-adiabatic dynamics methodology
2. Extreme-Condition Physics:
   • Multi-field coupling of pressure-light-temperature
   • Metastable material stability

Comprehensive Scoring

Key References

1. Bonev et al., J. Chem. Phys. 155, 54501 (2021)
   • Authors' prior work validating FPMD method
2. Evans et al., Chem. Mater. 18, 2520 (2006)
   • Key experiment discovering laser-reduced transition pressure
3. Lipp et al., Nat. Mater. 4, 211 (2005)
   • High-energy-density CO polymer characterization
4. Heyd et al., J. Chem. Phys. 118, 8207 (2003)
   • HSE hybrid functional, authors' recommended improvement direction
5. Noda et al., Comput. Phys. Commun. 235, 356 (2019)
   • SALMON software, LR-TDDFT implementation

Summary Commentary: This is a high-quality computational materials science paper that systematically clarifies the microscopic mechanism of photo-induced CO polymerization for the first time. The authors skillfully combine three complementary methods, obtaining physical insights at reasonable computational cost. While excited-state treatment involves simplifications, qualitative conclusions are reliable and highly consistent with experiment. The paper honestly discusses limitations and clearly identifies improvement directions (Ehrenfest dynamics, hybrid functionals). The work not only addresses a specific scientific question but also provides universally applicable insights for photo-assisted material synthesis, possessing significant academic and practical value. Subsequent work should supplement spectral calculations and implement more rigorous excited-state dynamics simulations.