point-group-symmetry-quantum - SKILL.md Agent Skill

name: point-group-symmetry-quantum description: "Point-group symmetry analysis of many-electron wavefunctions on quantum computers. Ancilla-free hybrid method for abelian and non-abelian groups using orbital rotations from representation matrix eigenvectors, tensor-network encoding, and error mitigation for molecular simulation. Application: quantum chemistry, drug discovery, materials science." arxiv: "2605.24824" categories: ["quantum-computing", "quantum-chemistry", "molecular-simulation"] keywords: ["point-group symmetry", "many-electron wavefunction", "quantum computer", "molecular simulation", "tensor-network", "error mitigation", "drug discovery", "quantum chemistry"]

Point-Group Symmetry Analysis on Quantum Computers

Based on arXiv:2605.24824 — "Point-group symmetry analysis of many-electron wavefunctions on a quantum computer" (Sakuma et al., May 2026).

Core Problem

Point groups (spatial symmetry operations in molecular systems) are essential for analyzing molecular orbitals and spectroscopy in chemistry. While quantum algorithms to exploit symmetry exist, practical implementations of point-group symmetry operations and detailed symmetry analysis of realistic many-electron wavefunctions were missing — especially for non-abelian groups and arbitrary basis functions.

Key Innovation

Ancilla-free hybrid method for analyzing point-group symmetries of many-electron states:

Works for both abelian and non-abelian groups
Calculates projection weights of irreducible representations
Uses orbital rotations derived from eigenvectors of representation matrices
Applicable to arbitrary basis functions (not restricted to symmetry-adapted basis)
Combines tensor-network based encoding with error mitigation for hardware execution

Methodology

Step 1: Representation Matrix Construction

For a given point group $G$ and molecular wavefunction $|\psi\rangle$:

Construct representation matrices $D(g)$ for each symmetry operation $g \in G$
Compute eigenvectors of $D(g)$ to define orbital rotations

Step 2: Orbital Rotation Application

Apply unitary orbital rotations to transform the wavefunction into symmetry-adapted basis
No ancilla qubits needed — the rotation is absorbed into the state preparation circuit

Step 3: Projection Weight Calculation

For each irreducible representation (irrep) $\Gamma$: $$w_\Gamma = \frac{d_\Gamma}{|G|} \sum_{g \in G} \chi_\Gamma(g)^* \langle \psi | U(g) | \psi \rangle$$ where $d_\Gamma$ is the irrep dimension, $\chi_\Gamma$ is the character, and $U(g)$ is the unitary representation.

Step 4: Tensor-Network Encoding + Error Mitigation

Use tensor-network based state encoding to compress the wavefunction representation
Apply Zero-Noise Extrapolation (ZNE) and measurement error mitigation on hardware
Achieves faithful reproduction of irrep weights within a few percent error

Implementation Results

Benzene (C₆H₆) — D₆ₕ symmetry

Up to 32 qubits on IBM ibm_kawasaki device
Ground state and first excited state symmetry analysis
Irrep weights reproduced within few percent error after error mitigation

Ferrocene (Fe(C₅H₅)₂) — D₅d symmetry

Numerical simulation for larger molecular system
Demonstrates scalability to transition metal complexes

Reusable Patterns

Pattern 1: Symmetry-Adapted Wavefunction Analysis

# Pseudocode for symmetry analysis pipeline
def analyze_symmetry(wavefunction, point_group):
    # 1. Get representation matrices for the group
    rep_matrices = get_representation_matrices(point_group, basis)
    
    # 2. Compute eigenvectors for orbital rotations
    rotations = compute_eigenvector_rotations(rep_matrices)
    
    # 3. Apply rotations (ancilla-free)
    rotated_wf = apply_orbital_rotations(wavefunction, rotations)
    
    # 4. Calculate projection weights for each irrep
    weights = {}
    for irrep in point_group.irreps:
        weights[irrep] = calculate_projection_weight(rotated_wf, irrep, point_group)
    
    return weights

Pattern 2: Tensor-Network State Encoding

Encode molecular wavefunction using Matrix Product States (MPS) or similar TN formats
Reduces qubit count requirements for large molecular systems
Compatible with VQE and other quantum chemistry algorithms

Pattern 3: Error Mitigation Stack for Chemistry

Readout error mitigation — calibration-based correction
Zero-Noise Extrapolation (ZNE) — Richardson extrapolation across noise levels
Symmetry verification — post-select on correct particle number / spin symmetry
Tensor-network denoising — exploit low-entanglement structure

Application Domains

Drug Discovery: Symmetry analysis of drug-target molecular complexes
Materials Science: Point-group classification of crystal structures and defects
Spectroscopy: Predicting and interpreting molecular spectra from quantum simulations
Catalysis: Symmetry properties of transition metal catalyst active sites
NISQ Chemistry: Practical quantum advantage pathway for near-term devices

Activation Keywords

point-group symmetry, molecular symmetry, quantum chemistry, wavefunction analysis, irreducible representation, tensor-network encoding, error mitigation, drug discovery, molecular orbital, spectroscopy, abelian group, non-abelian group

Pitfalls

Basis set dependency: The method works for arbitrary basis functions, but convergence requires sufficiently large basis sets for accurate symmetry analysis
Non-abelian groups: Require careful handling of multi-dimensional irreps and character orthogonality relations
Hardware noise: Without error mitigation, irrep weights can deviate significantly from true values — always apply the full error mitigation stack
State preparation: The quality of symmetry analysis depends on accurate state preparation; poor VQE convergence leads to incorrect symmetry assignments