pith. sign in

arxiv: 2604.03011 · v1 · submitted 2026-04-03 · ⚛️ physics.chem-ph · cond-mat.mtrl-sci

Regio-Connectivity and Torsional Angle Effects on Singlet Fission and SOCT-ISC in Aza-BODIPY Dimers

Sophiya Goyal , S. Rajagopala Reddy This is my paper

Pith reviewed 2026-05-13 18:15 UTC · model grok-4.3

classification ⚛️ physics.chem-ph cond-mat.mtrl-sci
keywords Aza-BODIPY dimersIntramolecular singlet fissionSOCT-ISCTorsional angleRegioisomersTriplet generationQuantum chemical calculations
0
0 comments X

The pith

Torsional angle between monomer units in aza-BODIPY dimers controls whether triplet generation occurs mainly through singlet fission or spin-orbit charge transfer intersystem crossing.

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

This study uses quantum chemical calculations to show how the shape of aza-BODIPY dimers affects their ability to generate triplet states. The key factor is the twist angle between the two monomer halves, which influences the energy levels and couplings needed for singlet fission or for charge-transfer assisted intersystem crossing. Different connection patterns between the units change these angles and lead to varying efficiencies. Dimer versions connected at positions 1 and 3 favor fission, while the 2-2 version boosts the alternative crossing route. These results highlight geometry as a design handle for molecules that produce triplets efficiently for applications like photovoltaics or photocatalysis.

Core claim

Triplet formation is strongly governed by the torsional angle (Φ) between monomer units, with regio-connectivity exerting a secondary influence. Dimers D[1,1] and D[1,3] exhibit favorable iSF energetics and coupling magnitudes, whereas D[2,2] displays low iSF rate constant but enhanced SOCT-ISC activity. The D[3,3] dimer shows exothermic multiexciton formation but reduced iSF efficiency due to destructive coupling interactions. The dominant ISC channel proceeds through the S1-T3 transition with large spin-orbit coupling and a small energy gap.

What carries the argument

Torsional angle Φ between monomer units, which sets the diabatic couplings for intramolecular singlet fission rates and the spin-orbit matrix elements for spin-orbit charge transfer intersystem crossing rates.

If this is right

  • D[1,1] and D[1,3] dimers achieve higher intramolecular singlet fission rates due to favorable energetics and couplings.
  • D[2,2] regioisomer shifts the dominant pathway to enhanced SOCT-ISC because of its geometry-induced properties.
  • D[3,3] shows reduced singlet fission efficiency from destructive interference in the couplings despite exothermic multiexciton states.
  • The S1-T3 transition dominates the intersystem crossing channel across the dimers owing to its large spin-orbit coupling and small energy gap.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Synthesizing dimers with locked or adjustable torsional angles could switch the preferred triplet-generation mechanism on demand.
  • The same geometry dependence may appear in other conjugated organic dimers and could be exploited to tune excited-state lifetimes.
  • Device-level testing of these isomers in thin films would reveal whether the calculated solution-phase preferences hold under operating conditions.

Load-bearing premise

The SA-XMCQDPT multireference calculations with the chosen active space and basis sets accurately capture the diabatic couplings, spin-orbit matrix elements, and relative energetics without significant errors from method limitations or missing dynamic correlation.

What would settle it

Experimental measurement of triplet formation rates and yields in synthesized samples of the D[1,1], D[1,3], D[2,2], and D[3,3] aza-BODIPY dimers, compared directly to the computed rate constants.

Figures

Figures reproduced from arXiv: 2604.03011 by Sophiya Goyal, S. Rajagopala Reddy.

Figure 1
Figure 1. Figure 1: Structure of the aza-BODIPY monomer (M) with atom numbering following IU [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
read the original abstract

Aza-BODIPY dimers represent promising molecular systems for efficient triplet-state generation through either intramolecular-singlet fission (iSF) or spin-orbit charge transfer intersystem crossing (SOCT-ISC). In this work, we investigate the role of molecular geometry in governing these mechanisms across four regioisomeric aza-BODIPY dimers (D[1,1], D[1,3], D[3,3], and D[2,2]) using multireference quantum-chemical calculations. Ground- and excited-state properties were analyzed at the MP2 and SA-XMCQDPT levels of theory, while diabatic couplings and spin-orbit matrix elements were evaluated to estimate iSF and SOCT-ISC rate constants, respectively. Our results reveal that triplet formation is strongly governed by the torsional angle ({\Phi}) between monomer units, with regio-connectivity exerting a secondary influence. Dimers D[1,1] and D[1,3] exhibit favorable iSF energetics and coupling magnitudes, whereas D[2,2] displays low iSF rate constant (kSF ) but enhanced SOCT-ISC activity. The D[3,3] dimer shows exothermic multiexciton formation but reduced iSF efficiency due to destructive coupling interactions. The dominant ISC channel proceeds through the S1-T3 transition with large spin-orbit coupling and a small energy gap. These findings provide critical mechanistic insights into geometry-dependent triplet generation in aza-BODIPY dimers.

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.

Referee Report

2 major / 3 minor

Summary. The manuscript computationally examines four regioisomeric aza-BODIPY dimers (D[1,1], D[1,3], D[3,3], D[2,2]) to assess how regio-connectivity and the inter-monomer torsional angle Φ control intramolecular singlet fission (iSF) versus spin-orbit charge-transfer intersystem crossing (SOCT-ISC). Ground- and excited-state properties are obtained at MP2 and SA-XMCQDPT levels; diabatic couplings and spin-orbit matrix elements are evaluated to estimate rate constants k_SF and k_ISC. The central claim is that Φ dominates over regio-connectivity, with D[1,1] and D[1,3] exhibiting favorable iSF energetics and couplings, D[2,2] showing suppressed k_SF but enhanced SOCT-ISC, and D[3,3] displaying exothermic multiexciton formation yet reduced iSF efficiency due to destructive interference; the dominant ISC channel is identified as S1–T3.

Significance. If the computed ordering of rates and the dominance of Φ prove robust, the work supplies concrete mechanistic guidance for engineering aza-BODIPY dimers toward selective triplet generation, relevant to singlet-fission photovoltaics and triplet photosensitizers. The explicit evaluation of diabatic couplings and SOCMEs from multireference wavefunctions constitutes a methodological strength that moves beyond purely phenomenological arguments.

major comments (2)
  1. [Computational Methods] Computational Methods section: the SA-XMCQDPT active-space selection is presented without convergence tests or explicit justification for the inclusion/exclusion of π/π* and nitrogen lone-pair orbitals. Because diabatic couplings and spin-orbit matrix elements are known to vary strongly with active-space choice in BODIPY-like chromophores, this omission directly undermines in the reported ordering of k_SF values and the claim that torsional angle Φ dominates regio-connectivity (abstract and Results on D[1,1]/D[1,3] versus D[2,2]).
  2. [Results] Results section (rate-constant estimation): the formulas used to convert computed energies, couplings, and SOCMEs into numerical k_SF and k_ISC values are not stated explicitly (no reference to Fermi’s golden rule, Marcus-Levich-Dogonadze parameters, or the precise prefactors). Without these equations or a sensitivity analysis, it is impossible to verify whether the quantitative ranking of the four dimers is an artifact of the chosen rate model.
minor comments (3)
  1. [Abstract] Abstract: the symbol 'kSF' contains an extraneous space; consistent subscript notation (k_SF) should be used throughout.
  2. [Figure 1] Figure 1 (or equivalent geometry figure): the definition and measurement convention for the torsional angle Φ are not stated in the caption or main text, making it difficult to reproduce the reported geometries.
  3. [Introduction] Introduction: several recent experimental and computational studies on aza-BODIPY singlet fission or SOCT-ISC are not cited, weakening the placement of the present results in the literature.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major point below and have revised the manuscript to incorporate the requested clarifications and supporting information.

read point-by-point responses

  1. Referee: [Computational Methods] Computational Methods section: the SA-XMCQDPT active-space selection is presented without convergence tests or explicit justification for the inclusion/exclusion of π/π* and nitrogen lone-pair orbitals. Because diabatic couplings and spin-orbit matrix elements are known to vary strongly with active-space choice in BODIPY-like chromophores, this omission directly undermines in the reported ordering of k_SF values and the claim that torsional angle Φ dominates regio-connectivity (abstract and Results on D[1,1]/D[1,3] versus D[2,2]).

    Authors: We acknowledge that the original manuscript did not include explicit convergence tests or a detailed justification for the SA-XMCQDPT active space. The active space was selected to encompass all relevant π/π* orbitals of the aza-BODIPY units together with the nitrogen lone-pair orbitals, following established protocols for BODIPY and aza-BODIPY systems in the literature. To strengthen the presentation, we have added a dedicated paragraph in the revised Computational Methods section that explains this choice with supporting references and reports a brief active-space convergence test (expansion by two additional virtual orbitals). The test confirms that the diabatic couplings, SOCMEs, and resulting rate ordering remain stable, thereby supporting the conclusion that torsional angle Φ is the dominant factor over regio-connectivity. revision: yes

  2. Referee: [Results] Results section (rate-constant estimation): the formulas used to convert computed energies, couplings, and SOCMEs into numerical k_SF and k_ISC values are not stated explicitly (no reference to Fermi’s golden rule, Marcus-Levich-Dogonadze parameters, or the precise prefactors). Without these equations or a sensitivity analysis, it is impossible to verify whether the quantitative ranking of the four dimers is an artifact of the chosen rate model.

    Authors: We agree that the explicit rate expressions were omitted from the original Results section. The k_SF and k_ISC values were obtained via Fermi’s golden rule in the weak-coupling limit, employing the Marcus-Levich-Dogonadze form for the Franck-Condon weighted density of states with the computed energy gaps, electronic couplings, and SOCMEs; reorganization energies were estimated from the MP2 geometry relaxations. In the revised manuscript we now state these equations explicitly (including the precise prefactors and the value of the electronic coupling threshold), together with a short sensitivity analysis demonstrating that the relative ordering of the four dimers is robust against ±20 % variations in reorganization energy and prefactor. This addition removes any ambiguity regarding the quantitative ranking. revision: yes

Circularity Check

0 steps flagged

No circularity: rate constants derived directly from computed matrix elements

full rationale

The derivation computes ground- and excited-state properties at MP2/SA-XMCQDPT, extracts diabatic couplings and spin-orbit matrix elements, then estimates kSF and ISC rates from those quantities across fixed geometries. No parameters are fitted to the target observables, no equations define a result in terms of itself, and no self-citations or imported uniqueness theorems carry the central claim about torsional-angle dominance. The ordering of D[1,1]/D[1,3] versus D[2,2] follows from explicit comparison of the computed energetics and couplings, which remain independent of the final mechanistic conclusion.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on standard quantum-chemistry approximations and method-specific choices rather than new postulates or fitted parameters.

axioms (2)
  • standard math Born-Oppenheimer approximation separating nuclear and electronic degrees of freedom
    Invoked implicitly in all ground- and excited-state geometry optimizations and energy evaluations.
  • domain assumption Selected active space in SA-XMCQDPT is adequate to describe the relevant singlet, triplet, and charge-transfer states
    Multireference calculations require an appropriate active space; its sufficiency is assumed without explicit justification in the abstract.

pith-pipeline@v0.9.0 · 5589 in / 1418 out tokens · 23759 ms · 2026-05-13T18:15:36.000462+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Design Principles for Singlet Fission in Aza-BODIPY Dimers: Spacer-Controlled Electronic Structure and Energy Ordering

    physics.chem-ph 2026-05 unverdicted novelty 5.0

    Spacers in aza-BODIPY dimers reduce excessive state mixing, promote distinct multiexcitonic states, and tune singlet fission rates over several orders of magnitude via substitution and conjugation patterns.

Reference graph

Works this paper leans on

8 extracted references · 8 canonical work pages · cited by 1 Pith paper

  1. [1]

    R.; Le, C.; Arias-Rotondo, D

    (1) Welin, E. R.; Le, C.; Arias-Rotondo, D. M.; McCusker, J. K.; MacMillan, D. W. Pho- tosensitized, energy transfer-mediated organometallic catalysis through electronically excited nickel (II).Science2017,355, 380–385. (2) Kim, T.; McCarver, S. J.; Lee, C.; MacMillan, D. W. Sulfonamidation of aryl and heteroaryl halides through photosensitized nickel cat...

    work page 2018

  2. [2]

    Organoboron molecules and polymers for organic solar cell applications.Chem

    (6) Miao, J.; Wang, Y.; Liu, J.; Wang, L. Organoboron molecules and polymers for organic solar cell applications.Chem. Soc. Rev.2022,51, 153–187. (7) Fernández-Moreira, V.; Thorp-Greenwood, F. L.; Coogan, M. P. Application of d6 transition metal complexes in fluorescence cell imaging.Chem Commun2010,46, 186–202. (8) Zhao, J.; Xu, K.; Yang, W.; Wang, Z.; Z...

    work page 2022

  3. [3]

    26 Table 13: The rate constant(kISC ins -1 units) of triplet state formation through the SOCT-ISC by varying the parameters reorganization energy (λin eV units), Haung-Rhys factor (Sin dimensionless units) and effective coupling (¯hωin cm-1 units) for D[2,2]. Non-orthogonal Orthogonal λ S¯hω k S1→T1 kS1→T2 kS1→T3 kS1→T1 kS1→T2 kS1→T3 0.1 0.3 700 2.69x10 −...

    work page 1968

  4. [4]

    (37) Dunning Jr, T. H. Gaussian Basis Sets for Use in Correlated Molecular Calculations. I. the Atoms Boron through Neon and Hydrogen.J. Chem. Phys.1989,90, 1007–1023. 29 (38) Granovsky,A.A.ExtendedMulti-configurationQuasi-degeneratePerturbationTheory: The New Approach to Multi-state Multi-reference Perturbation Theory.J. Chem. Phys. 2011,134, 214113. (39...

    work page 1989

  5. [5]

    Nonadiabatic Dynamics Simulations of Sin- glet Fission in 2,5-bis(fluorene-9-ylidene)-2,5-dihydrothiophene Crystals.Phys

    30 (46) Wibowo, M.; Persico, M.; Granucci, G. Nonadiabatic Dynamics Simulations of Sin- glet Fission in 2,5-bis(fluorene-9-ylidene)-2,5-dihydrothiophene Crystals.Phys. Chem. Chem. Phys.2019,21, 692–701. (47) Hartzler,D.A.; Slipchenko,L.V.; Savikhin,S.Triplet–tripletCouplinginChromophore Dimers: Theory and Experiment.J. Phys. Chem. A2018,122, 6713–6723. (4...

    work page 2019

  6. [6]

    K.; Kim, D.; Coropceanu, V.; Brédas, J.-L

    (59) Samanta, P. K.; Kim, D.; Coropceanu, V.; Brédas, J.-L. Up-Conversion Intersystem Crossing Rates in Organic Emitters for Thermally Activated Delayed Fluorescence: Impact of the Nature of Singlet vs Triplet Excited States.J. Am. Chem. Soc.2017, 139, 4042–4051. (60) Schmidt, K.; Brovelli, S.; Coropceanu, V.; Beljonne, D.; Cornil, J.; Bazzini, C.; Caronn...

    work page 2017

  7. [7]

    The optimized orbitals are delocalized. To S-1 Table S1: The VEEs for aza-BODIPY Dimers with Different Active Spaces 4,4 8,8 12,10 Dimer S 1 T1 ∆SF S 1 T1 ∆SF S 1 T1 ∆SF D[1,1] 1.88 1.06 0.24 1.85 1.09 0.33 1.68 0.95 0.22 D[1,3] 1.66 0.93 0.21 1.63 0.94 0.25 1.52 0.95 0.37 D[3,3] 1.64 0.75 -0.13 1.54 0.69 -0.15 1.49 0.53 -0.43 D[2,2] 1.78 1.13 0.49 1.76 1...

    work page 2020

  8. [8]

    Additionally, two doubly excited (DE) diabatic states|1(DE)1⟩and| 1(DE)2⟩were included due to energy considerations

    is localized in one of the aza-BODIPY moieties, and the charge transfer (CT) states|1(CA)⟩and| 1(AC)⟩, where C and A denote the radical cation and radical anion forms of the aza-BODIPY moiety, respectively. Additionally, two doubly excited (DE) diabatic states|1(DE)1⟩and| 1(DE)2⟩were included due to energy considerations. The electronic Hamiltonian (Hel) ...

    work page 1949