[Paper Review] Environment-Induced Exciton Renormalization in the Photosystem II Reaction Center

Tucker Allen, Barry Y. Li+2 more|arXiv (Cornell University)|Feb 23, 2026

Photosynthetic Processes and Mechanisms0 citations

TL;DR

The paper presents an ab initio, many-body Bethe–Salpeter equation study of the PSII reaction center with explicit protein environment using TDHF@vW, showing polarization-dependent exciton renormalization and changes in exciton delocalization.

ABSTRACT

Protein electrostatics tune excitation energies in the Photosystem II reaction center (PSII-RC), yet a fully quantum-mechanical many-body description of how the surrounding protein environment renormalizes excitons has remained computationally inaccessible. The Bethe-Salpeter equation (BSE) within many-body perturbation theory accurately describes excitonic physics through an explicit electron-hole interaction, but is prohibitively expensive for systems containing thousands of valence electrons. Here, we show that for sufficiently large systems the BSE becomes simpler to solve when treated with modern stochastic sampling techniques, as atomistic interactions self-average. In this regime, the effective electron-hole interaction mediated by the environment is governed by collective $k$-dependent polarization. These insights enable an ab initio study of the PSII-RC in which all six chlorins forming the hexameric dye core are treated explicitly together with a roughly seven Angstrom local protein environment. We directly compare the low-lying optical excitations of the isolated chromophore hexamer (1276 valence electrons) and the protein-dye cluster (3238 valence electrons). For $Q_y$ excitations near 680 nm, inclusion of the protein environment induces polarization-dependent energy shifts, redistributes spectral weight, and alters exciton delocalization and pigment character. Lateral and transverse asymmetries in the low-lying excited states are captured at the BSE level of theory. These results establish that we now have the tools for many-body calculations of biological nanostructures.

Motivation & Objective

• Motivated by how protein electrostatics tune excitation energies in PSII-RC.
• Aim to provide a fully quantum-mechanical many-body description including the protein environment.
• Demonstrate feasibility of large-scale BSE-like calculations with stochastic sampling.
• Compare isolated chromophore hexamer to protein-embedded system to quantify environmental effects.

Proposed method

• Use TDHF@v_W to replace explicit W with a translationally invariant screened exchange kernel v_W learned from stochastic fits.
• Apply real-time stochastic time-dependent Hartree propagation to capture k-dependent polarization without building the dielectric matrix.
• Treat long-wavelength exchange deterministically and high-k space via sparse-stochastic basis to compress exchange integrals.
• Diagonalize the two-particle Hamiltonian iteratively to go beyond Tamm-Dancoff approximation and include resonant-antiresonant coupling.
• Compute polarization-resolved optical spectra, transition densities, and participation ratios for isolated and embedded PSII-RC.

Experimental results

Research questions

• RQ1 How does explicit protein environment renormalize low-lying excitons in PSII-RC?
• RQ2 What polarization- and geometry-dependent shifts arise when embedding the pigment core in protein surroundings?
• RQ3 How does environment affect exciton delocalization and pigment character in PSII-RC?
• RQ4 Is there evidence of resonant-antiresonant coupling contributing to environment-induced renormalization?
• RQ5 Do results align with observed lateral and transverse asymmetries in PSII-RC excitations?

Key findings

• Protein embedding induces polarization-dependent energy shifts in low-lying excitons near 680 nm.
• The environment redistributes spectral weight and alters exciton delocalization and pigment character.
• Lateral and transverse asymmetries in excited states are captured at the BSE level and reflect environment effects.
• The bright Qy excitation shifts from 1.79 eV (isolated) to 1.82 eV (embedded), with a reduced participation ratio.
• Exciton characters show reduced delocalization and more localized density on the D1 branch in the embedded system.
• Tamm-Dancoff approximation misses the environment-induced shifts, highlighting the importance of resonant-antiresonant coupling.