Optical Excitations of Chlorophyll $a$ and $b$ Monomers and Dimers

Ask Hjorth Larsen; Bruce Forbes Milne; Duncan John Mowbray; Keenan Lyon; Mar\'ia Rosa Preciado-Rivas

arxiv: 1907.09430 · v2 · pith:SGRUEDSEnew · submitted 2019-07-22 · ❄️ cond-mat.mes-hall · physics.comp-ph· quant-ph

Optical Excitations of Chlorophyll a and b Monomers and Dimers

Mar\'ia Rosa Preciado-Rivas , Duncan John Mowbray , Keenan Lyon , Ask Hjorth Larsen , Bruce Forbes Milne This is my paper

Pith reviewed 2026-05-24 17:59 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall physics.comp-phquant-ph

keywords chlorophyllTDDFToptical absorptionphotoinduced dissociationLCAOlight-harvestingdensity functional theory

0 comments

The pith

LCAO-TDDFT with GLLB-SC correction on KS eigenenergies semi-quantitatively reproduces experimental photoinduced dissociation spectra for chlorophyll a and b monomers and dimers.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper shows that linear combination of atomic orbitals time-dependent density functional theory, run in either frequency or real-time mode, produces optical absorption spectra for chlorophyll a and b that match more costly plane-wave calculations. Adding the GLLB-SC derivative discontinuity correction to the Kohn-Sham energies brings both methods into semi-quantitative agreement with measured photoinduced dissociation data. The reduced computational cost opens first-principles treatment of optical excitations in larger biomolecular assemblies such as light-harvesting complexes.

Core claim

LCAO-TDDFT-k-ω and LCAO-TDDFT-r-t calculations with the GLLB-SC Δx correction applied to KS eigenenergies semi-quantitatively reproduce experimentally measured photoinduced dissociation results for Chl a and b monomers and dimers, while reducing computational cost compared to plane-wave methods.

What carries the argument

The GLLB-SC derivative discontinuity correction Δx applied to the KS eigenenergies within LCAO-TDDFT frameworks, which shifts the energies to improve agreement with measured optical excitations.

Load-bearing premise

The GLLB-SC correction developed elsewhere applies without change to chlorophyll a and b, and the isolated monomer and dimer models capture the excitations that occur inside real protein matrices.

What would settle it

A direct comparison showing that the corrected LCAO-TDDFT spectra deviate strongly from new experimental PID measurements on larger chlorophyll aggregates or on chlorophylls inside actual LHC II proteins would falsify the central claim.

read the original abstract

A necessary first step in the development of technologies such as artificial photosynthesis is understanding the photoexcitation process within the basic building blocks of naturally-occurring light harvesting complexes (LHCs). The most important of these building blocks in biological LHCs such as LHC II from green plants are the chlorophyll $a$ (Chl $a$) and chlorophyll $b$ (Chl $b$) chromophores dispersed throughout the protein matrix. However, efforts to describe such systems are still hampered by the lack of computationally efficient and accurate methods that are able to describe optical absorption in large biomolecules. In this work we employ a highly efficient linear combination of atomic orbitals (LCAOs) to represent the Kohn--Sham (KS) wave functions at the density functional theory (DFT) level and perform time dependent density functional theory (TDDFT) in either the reciprocal space and frequency domain (LCAO-TDDFT-$k$-$\omega$) or real space and time (LCAO-TDDFT-$r$-$t$) calculations of the optical absorption spectra of Chl $a$ and $b$ monomers and dimers. We find our LCAO-TDDFT-$k$-$\omega$ and LCAO-TDDFT-$r$-$t$ calculations reproduce results obtained with a plane wave (PW) representation of the KS wave functions (PW-TDDFT-$k$-$\omega$), but with a significant reduction in computational effort. Moreover, by applying the GLLB-SC derivative discontinuity correction $\Delta_x$ to the KS eigenenergies, with both LCAO-TDDFT-$k$-$\omega$ and LCAO-TDDFT-$r$-$t$ methods we are able to semi-quantitatively reproduce the experimentally measured photoinduced dissociation (PID) results. This work opens the path to first principles calculations of optical excitations in macromolecular systems.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

LCAO-TDDFT with GLLB-SC correction reproduces PW spectra for Chl a/b monomers and dimers at lower cost and reaches semi-quantitative match to PID data.

read the letter

The paper shows that LCAO-TDDFT in both frequency-domain and real-time versions matches plane-wave TDDFT results for chlorophyll a and b monomers and dimers while cutting computational cost, and that adding the GLLB-SC derivative discontinuity shift to the KS eigenvalues brings the calculated spectra into semi-quantitative line with experimental photoinduced dissociation measurements. The dimer spectra are new numbers not in the cited prior work. The direct PW benchmark is the strongest part of the evidence; it confirms the basis set is adequate without hidden adjustments. The correction is applied uniformly before the response calculation, so there is no circular fitting to the target data inside the paper itself. Efficiency gains are stated plainly and follow from the LCAO representation. The isolated-monomer and dimer models are a clear limitation. The paper treats them as a necessary first step, but the protein matrix in real light-harvesting complexes will shift the excitations, and that environment is not tested here. The match to experiment remains semi-quantitative, which is useful for trends but leaves room for further checks on absolute positions. This is for groups that need scalable TDDFT routes to optical spectra of biomolecular chromophores. Readers working on computational photochemistry or light-harvesting simulations would get concrete benchmarks and the dimer results from it. The calculations rest on external comparisons rather than internal fitting, so the work deserves a serious referee.

Referee Report

2 major / 3 minor

Summary. The paper claims that LCAO-based TDDFT calculations (both LCAO-TDDFT-k-ω and LCAO-TDDFT-r-t) reproduce plane-wave TDDFT results for the optical absorption spectra of isolated chlorophyll a and b monomers and dimers at significantly reduced computational cost. Applying the GLLB-SC derivative discontinuity correction Δx to the KS eigenenergies prior to the response calculation yields semi-quantitative agreement with experimental photoinduced dissociation spectra.

Significance. If validated, the work demonstrates a computationally efficient LCAO route to TDDFT spectra for chlorophyll chromophores, with explicit cross-validation against independent PW-TDDFT calculations and two distinct TDDFT implementations (frequency- and time-domain). The uniform, literature-derived GLLB-SC shift introduces no additional free parameters within the manuscript, which is a methodological strength. The isolated-monomer/dimer focus is appropriate to the stated claim and does not require the protein-matrix embedding to be addressed for the present results to stand.

major comments (2)

[Results] The central experimental claim rests on 'semi-quantitative' reproduction of PID spectra after the Δx shift. A table or figure panel that tabulates the positions and relative intensities of the main Q and Soret bands (monomer and dimer) against both PW-TDDFT and experiment, together with a quantitative measure of agreement (e.g., mean absolute deviation in eV), is needed to substantiate the assertion.
[Methods] § on computational methods: the manuscript must state the numerical value(s) of Δx applied to Chl a and Chl b and confirm that these values are taken unchanged from the cited GLLB-SC literature without any system-specific re-fitting or adjustment.

minor comments (3)

[Figures] Figure captions should explicitly label which curves correspond to LCAO-k-ω, LCAO-r-t, PW-TDDFT, and experiment so that the claimed reproduction is immediately verifiable.
[Abstract] The abstract states a 'significant reduction in computational effort' but provides no timing or scaling data; a brief statement or supplementary table comparing wall-clock times or basis sizes would be useful.
[Abstract] Notation: the symbol Δx is introduced without a preceding definition in the abstract; a parenthetical reminder of its origin (GLLB-SC) would improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive evaluation and the recommendation for minor revision. We address the two major comments below and will incorporate the requested changes into the revised manuscript.

read point-by-point responses

Referee: [Results] The central experimental claim rests on 'semi-quantitative' reproduction of PID spectra after the Δx shift. A table or figure panel that tabulates the positions and relative intensities of the main Q and Soret bands (monomer and dimer) against both PW-TDDFT and experiment, together with a quantitative measure of agreement (e.g., mean absolute deviation in eV), is needed to substantiate the assertion.

Authors: We agree that a tabulated quantitative comparison would strengthen the presentation of our results. In the revised manuscript we will add a new table (or expanded figure panel) listing the Q- and Soret-band positions and relative intensities for both Chl a and Chl b monomers and dimers. The table will compare the LCAO-TDDFT-k-ω and LCAO-TDDFT-r-t results (with and without the GLLB-SC shift) against the corresponding PW-TDDFT values and the experimental PID peak positions. We will also report the mean absolute deviation (in eV) between the calculated and experimental band positions to provide a clear numerical measure of agreement. revision: yes
Referee: [Methods] § on computational methods: the manuscript must state the numerical value(s) of Δx applied to Chl a and Chl b and confirm that these values are taken unchanged from the cited GLLB-SC literature without any system-specific re-fitting or adjustment.

Authors: We will revise the computational-methods section to state explicitly the numerical values of Δx applied to Chl a and Chl b. These values are taken directly from the cited GLLB-SC literature references and are used without any system-specific re-fitting or adjustment, consistent with our description of a uniform, literature-derived correction that introduces no additional free parameters. revision: yes

Circularity Check

0 steps flagged

No circularity detected; results are independent numerical benchmarks

full rationale

The paper computes optical spectra via LCAO-TDDFT (k-ω and r-t) for isolated Chl a/b monomers/dimers, benchmarks the basis against separate PW-TDDFT calculations, and applies the GLLB-SC Δx correction (taken from external prior literature) uniformly to KS eigenvalues before the response calculation. Target observables are direct comparisons to independent experimental PID spectra. No equations, parameters, or self-citations reduce the claimed spectra to quantities fitted or defined inside the paper itself. The isolated-model limitation is acknowledged as outside the explicit claim scope.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard DFT and TDDFT approximations drawn from prior literature together with the GLLB-SC correction also taken from earlier work. No new free parameters, ad-hoc axioms, or invented entities are introduced in the abstract.

axioms (2)

domain assumption Kohn-Sham DFT with standard approximate functionals yields sufficiently accurate ground-state wavefunctions for chlorophyll molecules.
Invoked by the use of DFT to generate the KS states that enter the TDDFT step.
domain assumption Linear-response or real-time TDDFT with the chosen functionals and correction produces excitation energies that can be compared to experimental PID spectra.
Core premise of the optical absorption calculations.

pith-pipeline@v0.9.0 · 5901 in / 1357 out tokens · 31883 ms · 2026-05-24T17:59:01.711438+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

by applying the GLLB-SC derivative discontinuity correction Δx to the KS eigenenergies... semi-quantitatively reproduce the experimentally measured photoinduced dissociation (PID) results
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

LCAO-TDDFT-k-ω and LCAO-TDDFT-r-t calculations

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.