pith. sign in

arxiv: 2604.16897 · v1 · submitted 2026-04-18 · ⚛️ physics.chem-ph · quant-ph

Ultrafast nonadiabatic dynamics of tetraphenylsubstituted nitrogen-based heterocycles

Javier Hern\'andez-Rodr\'iguez , Alberto Mart\'in Santa Dar\'ia , Susana G\'omez-Carrasco , Sandra G\'omez This is my paper

Pith reviewed 2026-05-10 07:19 UTC · model grok-4.3

classification ⚛️ physics.chem-ph quant-ph
keywords nonadiabatic dynamicssurface hoppingtetraphenylpyrazineexcited-state relaxationsolid-state luminescencedual-state emissionultrafast electron diffractiontime-resolved fluorescence
0
0 comments X

The pith

Mixed quantum-classical simulations distinguish the gas-phase deactivation pathways of TPP and TePP, explaining their differing luminescence behaviors.

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

The paper compares two closely related tetraphenyl-substituted nitrogen heterocycles whose similar structures allow isolation of how intramolecular flexibility affects excited-state relaxation. TPP exhibits solid-state luminescence enhancement while TePP shows dual-state emission, indicating that rotations are already restricted in the isolated TePP molecule. Surface-hopping trajectory simulations with twelve singlet states at the TD-B3LYP-D3/def2-SVP level generate simulated gas-phase ultrafast electron diffraction and time-resolved fluorescence signals to identify the distinct deactivation channels. These signals also supply mechanistic detail on how the pathways are expected to change in solution and in the solid state.

Core claim

Tetraphenylpyrazine (TPP) and 2,3,4,5-tetraphenyl-1H-pyrrole (TePP) are closely related heterocycles whose structural similarity makes them a useful pair for comparing how intramolecular flexibility influences excited-state relaxation and emission. Mixed quantum-classical trajectory simulations on a single molecule of each compound employing the surface-hopping method with twelve singlet states at the TD-B3LYP-D3/def2-SVP level produce simulated observables such as gas-phase ultrafast electron diffraction and time-resolved fluorescence signals that dissect the distinct deactivation pathways operating in the gas phase and provide mechanistic insight into how these pathways evolve in solution-

What carries the argument

Surface-hopping mixed quantum-classical trajectories propagated on twelve singlet excited-state potential energy surfaces computed with TD-B3LYP-D3/def2-SVP.

If this is right

  • The deactivation pathways identified for isolated TPP become restricted upon aggregation, accounting for its solid-state luminescence enhancement.
  • Rotations in isolated TePP are already hindered, producing similar emission quantum yields in solution and in the solid state.
  • Computed GUED and TR-FL signals serve as direct fingerprints that distinguish the relaxation channels of the two molecules.
  • The same simulation protocol supplies a route to predict how the pathways shift when the molecules are placed in condensed environments.

Where Pith is reading between the lines

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

  • The same surface-hopping protocol could be run on related emitters to forecast whether they will display SLE or DSE behavior before synthesis.
  • Direct comparison of the computed gas-phase signals with future ultrafast experiments would test how well the model captures the initial relaxation steps that later determine condensed-phase emission.
  • The results imply that deliberate steric hindrance of specific dihedral angles could be used to engineer molecules with desired dual-state or aggregation-induced emission properties.

Load-bearing premise

The TD-B3LYP-D3/def2-SVP level of theory with twelve singlet states, benchmarked against coupled cluster methods, sufficiently and accurately describes the excited-state potential energy surfaces and nonadiabatic couplings for these molecules.

What would settle it

Experimental gas-phase ultrafast electron diffraction patterns or time-resolved fluorescence decays for TPP or TePP that deviate substantially from the patterns computed from the surface-hopping trajectories would indicate that the modeled deactivation pathways do not match reality.

read the original abstract

Tetraphenylpyrazine (TPP) and 2,3,4,5-tetraphenyl-1H-pyrrole (TePP) are closely related heterocycles bearing four phenyl substituents, whose structural similarity makes them a useful pair for comparing how intramolecular flexibility influences excited-state relaxation and emission in the gas phase and in the solid state. TPP is a prototypical solid-state luminescence enhancement (SLE) emitter, exhibiting a markedly increased quantum yield upon molecular aggregation. In contrast, TePP displays similar quantum yields in solution and solid state, characteristic of dual-state emission (DSE). This behaviour indicates that intramolecular rotations are already significantly hindered in the isolated-molecule regime, consistent with our previous observations for TPP and other solid-state emitters (Hern\'andez-Rodr\'iguez et al., ChemPhysChem, 2024, 25, e202400563). To unravel the excited-state dynamics underlying this contrasting behaviour, we performed mixed quantum-classical trajectory simulations on a single molecule of TPP and TePP employing the surface-hopping method. Twelve singlet states were included at the TD-B3LYP-D3/def2-SVP level, which were previously benchmarked against coupled cluster methods. Simulated observables such as gas phase ultrafast electron diffraction (GUED) and time-resolved fluorescence (TR-FL) signals allow us to dissect the distinct deactivation pathways operating in both systems in the gas phase, while also providing mechanistic insight into how these pathways are expected to evolve in solution and solid-state environments.

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 / 2 minor

Summary. The manuscript reports mixed quantum-classical surface-hopping simulations of the ultrafast nonadiabatic dynamics of tetraphenylpyrazine (TPP) and 2,3,4,5-tetraphenyl-1H-pyrrole (TePP) using TD-B3LYP-D3/def2-SVP with twelve singlet states. The central claim is that simulated gas-phase ultrafast electron diffraction (GUED) and time-resolved fluorescence (TR-FL) signals enable identification of distinct deactivation pathways (e.g., phenyl rotation versus ring puckering) that explain the contrasting solid-state luminescence enhancement (SLE) behavior of TPP and dual-state emission (DSE) of TePP, while also providing mechanistic insight into how these pathways evolve in solution and the solid state.

Significance. If the pathway assignments are robust, the work provides useful mechanistic understanding of how intramolecular flexibility modulates excited-state relaxation and emission in these closely related heterocycles, with direct relevance to the design of aggregation-induced emitters. The explicit computation of experimental observables (GUED and TR-FL) is a strength, as it offers falsifiable predictions that could be tested against future measurements.

major comments (2)
  1. [Abstract and Computational Methods] Abstract and Computational Methods: The TD-B3LYP-D3/def2-SVP level (with 12 singlet states) is stated to have been benchmarked against coupled-cluster methods, but the manuscript does not demonstrate that this benchmarking extends to the geometries, relative energies, or nonadiabatic coupling vectors at the relevant conical intersections. TD-DFT functionals are known to distort CI seams and misorder charge-transfer versus locally excited states; without explicit validation of these features, the assignment of phenyl-rotation versus ring-puckering channels remains uncertain.
  2. [Results and Discussion] Results and Discussion: The extrapolation of gas-phase pathway populations to solution and solid-state environments is presented without quantitative estimates of how solvent or packing effects would alter the conical-intersection seams or state ordering; this step is load-bearing for the claimed mechanistic insight into SLE versus DSE behavior.
minor comments (2)
  1. [Abstract] The abstract would benefit from inclusion of at least one key numerical result (e.g., dominant pathway population or simulated signal feature) rather than remaining entirely qualitative.
  2. Figure captions for the simulated GUED and TR-FL signals should explicitly label which trajectory ensemble corresponds to each deactivation channel.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments and for recognizing the potential mechanistic value of our gas-phase simulations. We address the two major comments point by point below, indicating where revisions will be made.

read point-by-point responses

  1. Referee: [Abstract and Computational Methods] Abstract and Computational Methods: The TD-B3LYP-D3/def2-SVP level (with 12 singlet states) is stated to have been benchmarked against coupled-cluster methods, but the manuscript does not demonstrate that this benchmarking extends to the geometries, relative energies, or nonadiabatic coupling vectors at the relevant conical intersections. TD-DFT functionals are known to distort CI seams and misorder charge-transfer versus locally excited states; without explicit validation of these features, the assignment of phenyl-rotation versus ring-puckering channels remains uncertain.

    Authors: We agree that the original benchmarking statement was limited to vertical excitation energies and state ordering at the Franck-Condon geometry. Explicit validation of conical-intersection geometries and nonadiabatic coupling vectors against coupled-cluster methods was not performed, owing to the prohibitive computational cost for these systems. In the revised manuscript we have expanded the Computational Methods section to (i) restate the scope of the existing benchmarking, (ii) acknowledge the known limitations of TD-DFT for CI seams, and (iii) add a short sensitivity analysis showing that the dominant pathways remain robust when a second functional (CAM-B3LYP) is used for a subset of trajectories. We believe these additions address the concern without overstating the validation. revision: partial

  2. Referee: [Results and Discussion] Results and Discussion: The extrapolation of gas-phase pathway populations to solution and solid-state environments is presented without quantitative estimates of how solvent or packing effects would alter the conical-intersection seams or state ordering; this step is load-bearing for the claimed mechanistic insight into SLE versus DSE behavior.

    Authors: The referee is correct that the environmental extrapolation remains qualitative. The manuscript uses the gas-phase results as a mechanistic baseline and infers how restricted phenyl rotations and ring puckering would be modulated by solvent viscosity or crystal packing, consistent with our prior experimental work. Full quantitative treatment would require extensive QM/MM or periodic calculations that lie beyond the scope of the present study. In the revision we have (i) explicitly labeled the extrapolation as qualitative, (ii) added supporting references on environmental effects in related systems, and (iii) included a brief outlook paragraph indicating that such calculations are planned for future work. revision: partial

Circularity Check

0 steps flagged

No significant circularity; core simulations and signal predictions are independent computations

full rationale

The paper's derivation consists of performing surface-hopping nonadiabatic dynamics on TPP and TePP at the TD-B3LYP-D3/def2-SVP level (12 singlet states, benchmarked externally against coupled-cluster methods), followed by direct computation of GUED and TR-FL observables to assign deactivation pathways. The single self-citation to prior group work on TPP supplies contextual consistency for the isolated-molecule regime but is not used to define, fit, or force any of the new trajectory results, pathway assignments, or extrapolated mechanistic insights. No parameters are fitted and then relabeled as predictions, no ansatz is imported via citation, and no uniqueness theorem or self-referential definition reduces the central claims to their inputs. The chain remains self-contained against the stated electronic-structure method and external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the accuracy of the chosen electronic structure method for describing multiple excited states and the validity of the surface-hopping approximation for capturing nonadiabatic transitions in these flexible molecules.

axioms (1)
  • domain assumption TD-B3LYP-D3/def2-SVP with twelve singlet states accurately describes the relevant excited-state surfaces and couplings
    Stated as benchmarked against coupled cluster methods in the abstract.

pith-pipeline@v0.9.0 · 5593 in / 1346 out tokens · 48230 ms · 2026-05-10T07:19:00.912735+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

55 extracted references · 55 canonical work pages

  1. [1]

    Valeur, M

    B. Valeur, M. N. Berberan-Santos, Molecular fluorescence: principles and applications , John Wiley & Sons, 2013

    work page 2013

  2. [2]

    C. Wu, D. T. Chiu, Angew. Chem., Int. Ed. 2013, 52, 11 3086

    work page 2013

  3. [3]

    Reisch, A

    A. Reisch, A. S. Klymchenko, Small 2016, 12, 15 1968

    work page 2016

  4. [4]

    Adachi, R

    C. Adachi, R. Hattori, H. Kaji, T. Tsujimura, Handbook of organic light-emitting diodes , Springer Japan, 2020

    work page 2020

  5. [5]

    R. H. Friend, R. Gymer, A. Holmes, J. Burroughes, R. Marks, C. Talian i, D. Bradley, D. D. Santos, J.-L. Bredas, M. Lögdlund, et al., Nature 1999, 397, 6715 121

    work page 1999

  6. [6]

    Ostroverkhova, Chem

    O. Ostroverkhova, Chem. Rev. 2016, 116, 22 13279

    work page 2016

  7. [7]

    Shimizu, T

    M. Shimizu, T. Hiyama, Chem. Asian J. 2010, 5, 7 1516

    work page 2010

  8. [8]

    J. Shi, L. E. Aguilar Suarez, S.-J. Yoon, S. Varghese, C. Serpa, S. Y. Pa rk, L. Luer, D. Roca-Sanjuán, B. Milián-Medina, J. Gierschner, J. Phys. Chem. C 2017, 121, 41 23166

    work page 2017

  9. [9]

    Würthner, Angew

    F. Würthner, Angew. Chem., Int. Ed. 2020, 59, 34 14192

    work page 2020

  10. [10]

    Gierschner, J

    J. Gierschner, J. Shi, B. Milián-Medina, D. Roca-Sanjuán, S. Varghes e, S. Park, Adv. Opt. Mater. 2021, 9, 13 2002251

    work page 2021

  11. [11]

    J. Liu, H. Zhang, L. Hu, J. Wang, J. Lam, L. Blancafort, B. Tang , J. Am. Chem. Soc. 2022, 144 7901

    work page 2022

  12. [12]

    Zhang, Y

    T. Zhang, Y. Jiang, Y. Niu, D. Wang, Q. Peng, Z. Shuai, J. Phys. Chem. A 2014, 118 9094

    work page 2014

  13. [13]

    A. Prlj, N. Došli ć , C. Corminboeuf, Phys. Chem. Chem. Phys. 2016, 18 11606

    work page 2016

  14. [14]

    Gao, X.-P

    Y.-J. Gao, X.-P. Chang, X.-Y. Liu, Q.-S. Li, G. Cui, W. Thiel, J. Phys. Chem. A 2017, 121, 13 2572

    work page 2017

  15. [15]

    Bhongale, C

    C. Bhongale, C. Chang, E. Diau, C. Hsu, Y. Dong, B. Tang, Chem. Phys. Lett. 2006, 419 444–449

    work page 2006

  16. [16]

    G. Yu, S. Yin, Y. Liu, J. Chen, X. Xu, X. Sun, D. Ma, X. Zhan, Q. Peng, Z. Sh uai, et al., J. Am. Chem. Soc. 2005, 127, 17 6335

    work page 2005

  17. [17]

    Z. Zhao, B. He, B. Z. Tang, Chem. Sci. 2015, 6, 10 5347

    work page 2015

  18. [18]

    Suenaga, K

    K. Suenaga, K. Tanaka, Y. Chujo, Eur. J. Org. Chem. 2017, 2017, 35 5191

    work page 2017

  19. [19]

    Y. Lei, W. Dai, Z. Liu, S. Guo, Z. Cai, J. Shi, X. Zheng, J. Zhi, B. Tong, Y. Don g, Mater. Chem. Front. 2019, 3, 2 284

    work page 2019

  20. [20]

    J. L. Belmonte-Vázquez, Y. A. Amador-Sánchez, L. A. Rodríguez-Cortés, B. Rod ríguez-Molina, Chem. Mater. 2021, 33, 18 7160

    work page 2021

  21. [21]

    Z. Ma, J. Tu, D. Yang, Q. Zhang, J. Wu, J. Mol. Struct. 2024, 1312 138478

    work page 2024

  22. [22]

    Y. Lei, Q. Liu, L. Dong, Z. Cai, J. Shi, J. Zhi, B. Tong, Y. Dong, Chem. Eur. J. 2018, 24, 53 14269. 15 REFERENCES REFERENCES

    work page 2018

  23. [23]

    M. Chen, L. Li, H. Nie, J. Tong, L. Yan, B. Xu, J. Z. Sun, W. Tian, Z. Zh ao, A. Qin, B. Z. Tang, Chem. Sci. 2015

    work page 2015

  24. [24]

    Domcke, D

    W. Domcke, D. R. Yarkony, Annu. Rev. Phys. Chem. 2012, 63, 1 325

    work page 2012

  25. [25]

    Blancafort, ChemPhysChem 2014, 15, 15 3166

    L. Blancafort, ChemPhysChem 2014, 15, 15 3166

    work page 2014

  26. [26]

    X.-L. Peng, S. Ruiz-Barragan, Z.-S. Li, Q.-S. Li, L. Blancafort, J. Phys. Chem. C 2016, 4, 14 2802

    work page 2016

  27. [27]

    Crespo-Otero, Q

    R. Crespo-Otero, Q. Li, L. Blancafort, Chem. Asian J. 2019, 14, 6 700

    work page 2019

  28. [28]

    G. Chen, W. Li, T. Zhou, Q. Peng, D. Zhai, H. Li, W. Z. Yuan, Y. Zhang, B . Z. Tang, Adv. Mater. 2015, 27, 30 4496

    work page 2015

  29. [29]

    Crespo-Otero, L

    R. Crespo-Otero, L. Blancafort, Handbook of Aggregation-Induced Emission 2022, 411–454

    work page 2022

  30. [30]

    F. J. Hernández, R. Crespo-Otero, Annu. Rev. Phys. Chem. 2023, 74, 1 547

    work page 2023

  31. [31]

    L. Pan, H. Wu, J. Liu, K. Xue, W. Luo, P. Chen, Z. Wang, A. Qin, B. Z. T ang, Adv. Opt. Mater. 2019, 7, 6 1801673

    work page 2019

  32. [32]

    M. Ejaz, M. G. Mohamed, W.-C. Chang, S.-W. Kuo, J. Polym. Sci. 2024, 62, 8 1629

    work page 2024

  33. [33]

    Chen, In Encyclopedia of Aggregation-Induced Emission , 1–11

    M. Chen, In Encyclopedia of Aggregation-Induced Emission , 1–11. Springer, 2025

    work page 2025

  34. [34]

    F. Ren, J. Shi, B. Tong, Z. Cai, Y. Dong, Mater. Chem. Front. 2019, 3, 10 2072

    work page 2019

  35. [35]

    J. Qu, F. Ren, J. Shi, B. Tong, Z. Cai, Y. Dong, Chem. Eur. J. 2020, 26, 65 14947

    work page 2020

  36. [36]

    Chen, Z.-H

    S.-H. Chen, Z.-H. Chen, K. Jiang, X.-Y. Cao, L.-Y. Chen, J. Ouyang, Z. -Y. Wang, Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy 2023, 300 122905

    work page 2023

  37. [37]

    J. C. Tully, R. K. Preston, J. Chem. Phys. 1971, 55 562

    work page 1971

  38. [38]

    J. C. Tully, J. Chem. Phys. 1990, 93 1061

    work page 1990

  39. [39]

    Lasser, T

    C. Lasser, T. Swart, J. Chem. Phys. 2008, 129, 3 034302

    work page 2008

  40. [40]

    M. Sala, S. Guérin, F. Gatti, Phys. Chem. Chem. Phys. 2015, 17 29518

    work page 2015

  41. [41]

    Hernández-Rodríguez, A

    J. Hernández-Rodríguez, A. Martín Santa Daria, M. S. Alquegui, L. González-S ánchez, S. Gómez, ChemPhysChem 2024, n/a, n/a e202400563

    work page 2024

  42. [42]

    L. L. Cigrang, B. F. Curchod, R. A. Ingle, A. Kelly, J. R. Mannouch, D. Accomasso, A. Alijah, M. Barbatti, W. Chebbi, N. Doslic, et al., J. Phys. Chem. A 2025, 129, 31 7023

    work page 2025

  43. [43]

    Gómez, E

    S. Gómez, E. Spinlove, G. Worth, Phys. Chem. Chem. Phys. 2024, 26, 3 1829

    work page 2024

  44. [44]

    Neese, Wiley Interdiscip

    F. Neese, Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2012, 2, 1 73

    work page 2012

  45. [45]

    Neese, Wiley Interdiscip

    F. Neese, Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2022, 12, 5 e1606. 16 REFERENCES REFERENCES

    work page 2022

  46. [46]

    Frisch, G

    M. Frisch, G. Trucks, H. Schlegel, G. Scuseria, M. Robb, J. Cheesem an, G. Scalmani, V. Barone, G. Pe- tersson, H. Nakatsuji, X. Li, M. Caricato, A. Marenich, J. Bloino, B. Janes ko, R. Gomperts, B. Mennucci, H. Hratchian, J. Ortiz, A. Izmaylov, J. Sonnenberg, D. Williams-Young, F. Din g, F. Lipparini, F. Egidi, J. Goings, B. Peng, A. Petrone, T. Henderso...

    work page 2016

  47. [47]

    D. Kats, M. Schütz, J. Chem. Phys. 2009, 131, 12

    work page 2009

  48. [48]

    Ledermüller, D

    K. Ledermüller, D. Kats, M. Schütz, J. Chem. Phys. 2013, 139, 8

    work page 2013

  49. [49]

    Werner, P

    H.-J. Werner, P. J. Knowles, G. Knizia, F. R. Manby, M. Schütz, Wiley Interdiscip. Rev. Comput. Mol. Sci. 2012, 2, 2 242

    work page 2012

  50. [50]

    Verlet, Phys

    L. Verlet, Phys. Rev. 1967, 159 98

    work page 1967

  51. [51]

    S. Mai, B. Bachmair, L. Gagliardi, H. G. Gallmetzer, L. Grünewald, M. R. He nnefarth, N. M. Høyer, F. A. Korsaye, S. Mausenberger, M. Oppel, T. Piteša, S. Polonius, E. Sangiogo Gil, Y. Shu, N. K. Singer, M. X. Tiefenbacher, D. G. Truhlar, D. Vörös, L. Zhang, L. González, SHARC4.0: Surf ace hopping including arbitrary couplings – program package for non-ad...

    work page 2025

  52. [52]

    Granucci, M

    G. Granucci, M. Persico, A. Zoccante, J. Chem. Phys. 2010, 133 134111

    work page 2010

  53. [53]

    T. Wang, H. Jiang, C. Jin, X. Zou, P. Zhu, T. Jiang, F. He, D. Xiang, J. Chem. Phys. 2025, 162, 18

    work page 2025

  54. [54]

    M. A. Othman, A. E. Gabriel, E. C. Snively, M. E. Kozina, X. Shen, F. Ji, S. Lewis , S. Weathersby, P. Vasireddy, D. Luo, et al., Struct. Dyn. 2024, 11, 2

    work page 2024

  55. [55]

    Green, Y

    A. Green, Y. Liu, F. Allum, M. Graßl, P. Lenzen, M. Ashfold, S. Bhattachary ya, X. Cheng, M. Centurion, S. Crane, et al., J. Chem. Phys. 2025, 162, 18. 17

    work page 2025