Pith. sign in

REVIEW 3 major objections 7 minor 49 references

Hidden nonlinear optical susceptibilities in linear polaritonic spectra

T0 review · 3 major / 7 minor · reviewed 2026-08-12 · deepseek-v4-flash

Pith's one-line read Molecular polariton linear spectra carry a hierarchy of finite-size 1/N quantum corrections from vacuum-mediated Raman processes, visible when cavity loss is comparable to the single-molecule coupling.

desk verdict A credible 1/N expansion showing vacuum-mediated Raman sidebands as finite-size corrections to polaritonic linear spectra, but the printed continued fraction has a dagger-ordering typo and the central formula lacks a numerical check. read the letter →

arxiv 2411.09039 v2 pith:3UO3I6SY submitted 2024-11-13 quant-ph cond-mat.otherphysics.chem-phphysics.optics

classification quant-phcond-mat.otherphysics.chem-phphysics.optics
keywords molecularpolaritonsstronglight-mattercoupling1/Nexpansionvacuum-mediatedRamanscatteringlinearresponsephotonGreen'sfunctionhigh-Qcavitiesnonlinearsusceptibilities
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Molecular polaritons — N molecules strongly coupled to one cavity photon mode — have linear spectra that are usually captured by classical linear optics. This paper argues that this classical description is only the thermodynamic-limit term of an exact 1/N expansion, and that the first quantum correction contains vacuum-mediated molecular Raman processes that produce weak sidebands around every polariton peak. Those sidebands are hidden in ordinary low-Q cavities because the cavity linewidth κ vastly exceeds the single-molecule coupling λ, but they should become visible in high-Q single-mode cavities with κ ≈ λ. If true, this explains where the expected nonlinear many-body effects of strong coupling live in linear spectroscopy and gives a concrete experimental regime in which to observe them.

What carries the argument

The central object is the block-tridiagonal representation of the single-excitation Hamiltonian, a nearest-neighbor chain of photonic and molecular blocks that the paper calls the CUT-E diagram. This chain separates fast collective couplings $V_n\propto\lambda\sqrt{N}$ — Rayleigh processes that conserve the number of ground-state phonons — from slow single-molecule couplings $v_n\propto\lambda$ — Raman processes that create or destroy one vibrational quantum. Matrix continued-fraction techniques applied to this chain produce the photon Green's function, and the timescale separation organizes the $1/N$ expansion that isolates the hidden sidebands.

What would settle it

Measure the linear absorption or transmission of a high-Q single-mode microcavity containing a known number N of identical molecules at low temperature, with κ ≈ λ and the collective coupling λ√N kept fixed while N is varied. The paper predicts sidebands offset from each zeroth-order polariton by the ground-state vibrational frequency, with peak height ∝1/N and blue shift Δ = 2λ(√N − √(N−1)). If those sidebands are absent, do not scale with N, or persist unchanged in a multimode cavity, the central claim fails.

Watch

Extended reading notes

Core claim

The paper derives the exact photon Green's function $D^R_N(\omega)$ for a single-mode cavity coupled to N vibronic molecules as a continued fraction (Eq. 6) and expands it in powers of $1/N$. The zeroth-order term $d_{N,0}$ reproduces classical linear optics: the self-energy is just the molecular linear susceptibility $\chi^{(1)}$, and in the limit $N\to\infty$ nothing else survives. The first-order term $d_{N,1}$ (Eq. 8) is the new claim: it factorizes into $d_{N,0}\cdot V_0G_{e,0}v_0\bigl(\omega-H_{\mathrm{ph},1}+i\kappa/2-V_1G_{e,1}V_1^\dagger\bigr)^{-1} v_0^\dagger G_{e,0} V_0^\dagger\cdot d_{N,0}$, which describes a Stokes Raman transition, propagation through first-order polaritons formed by the remaining $N-1$ molecules, and an anti-Stokes transition back. These vacuum-mediated processes are each penalized by $1/\sqrt{N}$, so the sidebands have height $\propto 1/N$ and blue shift $\Delta = 2\lambda(\sqrt{N}-\sqrt{N-1})$. The same continued fraction is rewritten (Eq. 9) as the photon self-energy equal to a sum of irreducible Rayleigh and Raman nonlinear susceptibilities of the molecular ensemble, making the hidden nonlinear character of the linear spectrum explicit.

Load-bearing premise

The prediction that the sidebands are visible assumes a perfectly identical set of molecules at zero temperature with only simple Markovian loss and a single cavity mode; if static disorder, thermal vibrations, or extra cavity modes blur the resonances, the 1/N features could be washed out.

Editorial extensions

If this is right

  • In the thermodynamic limit ($N\to\infty$) the polariton linear response becomes exactly the classical linear-optics result built from the molecular linear susceptibility, so no nonlinearity is visible.
  • In low-$Q$ cavities ($\kappa\gg\lambda$) the $1/N$ terms are suppressed both by the timescale separation and by limited spectral resolution, which is why transfer-matrix and effective-medium models work so well.
  • In high-$Q$ single-mode cavities with $\kappa\sim\lambda$, each zeroth-order polariton acquires sidebands shifted by the ground-state vibrational frequency $\omega_v$; their height scales as $1/N$ and their Rabi splitting is $2\lambda\sqrt{N-1}$.
  • The photon self-energy can be written as a sum over irreducible Rayleigh and Raman nonlinear susceptibilities, so a linear measurement in this regime directly carries nonlinear molecular response information.
  • For two molecular species, the $O(N^{-2})$ terms produce sidebands at $2\omega_{v,A}$, $2\omega_{v,B}$, and $\omega_{v,A}+\omega_{v,B}$, revealing collective Raman processes involving different molecules.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the sidebands survive at the predicted $1/N$ height, a careful measurement of sideband intensity versus $N$ would give a direct estimate of the single-molecule coupling $\lambda$, which is otherwise hard to isolate in the collective strong-coupling regime.
  • The two-species sum-frequency sideband at $\omega_{v,A}+\omega_{v,B}$ suggests using linear transmission as a probe of cooperative Raman events between distinct molecular species, a measurement that does not require a pulsed nonlinear experiment.
  • The paper's own remark that multimode cavities blur the sidebands implies the observable doubles as a diagnostic: the presence of sharp sidebands certifies effective single-mode behavior, while their absence could indicate multimode or disordered broadening.
  • Because the first-order sidebands are genuine Raman transitions driven by vacuum fluctuations, one might expect weak nonclassical photon statistics or entanglement in the transmitted field at those frequencies; the paper leaves this as future work.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 7 minor

Summary. This manuscript asks why molecular polaritons' linear spectra are so well described by classical linear optics, despite the expected nonlinear character of strong coupling. Using a bosonic (Holstein–Tavis–Cummings) model restricted to the first excitation manifold, the authors express the photon Green's function as a matrix continued fraction (Eq. 6) over the block-tridiagonal 'CUT-E' chain, in which V-type couplings are collective (∝ λ√N) and v-type couplings are single-molecule Raman processes (∝ λ). They then develop a 1/N expansion of the Green's function. The zeroth order reproduces the classical linear-optics formula in terms of the ensemble linear susceptibility; the first-order term, dN,1 (Eq. 8), is interpreted as a vacuum-mediated Stokes–anti-Stokes process and predicts Raman sidebands at ±ωv around the polaritons, with height ~1/N and reduced first-order Rabi splitting 2λ√(N−1). The authors argue these sidebands are washed out in low-Q and multimode cavities but should be observable in high-Q single-mode cavities with κ ~ λ. They also classify the Dyson diagrams into reducible and irreducible contributions and claim (Eq. 9) that the photon self-energy equals the sum of irreducible odd-order nonlinear susceptibilities of the molecular ensemble, with O(N^{-2}) terms generating two-molecule Raman sidebands.

Significance. The central claim is significant and cleanly falsifiable: linear spectroscopy of high-Q single-mode molecular polaritons should reveal finite-size Raman sidebands with a specific N-scaling, providing a molecule-specific signature of nonlinear response mediated by vacuum fluctuations. If Eq. (8) holds, the paper resolves an active conceptual puzzle (polaritons as 'optical filters') by showing where the nonlinearity goes in the N→∞ limit. Strengths: the derivation is self-contained and parameter-free; the N-counting is consistent throughout; and the reduction to classical optics in the thermodynamic limit is a genuine consistency check rather than a circular assumption. The weaknesses are local but load-bearing: (i) the printed continued fraction (Eq. 6 and SI Eqs. S18, S22) uses a dagger ordering opposite to that of Eq. (5b), so Eq. (8) does not follow from the displayed equations as written; (ii) no numerical comparison of the truncated 1/N expansion with the exact Green's function is provided; and (iii) the proof that the self-energy equals the sum of irreducible susceptibilities (Eq. 9) is delegated to the reader. These issues are repairable within the scope of the manuscript.

major comments (3)
  1. [Eq. (6); SI Eqs. (S18), (S22)] The continued fraction printed in Eq. (6), and repeated in SI Eqs. (S18) and (S22), is not the iteration of the recursion in Eq. (5b). Since v0 is the (He,0; Hph,1) block in Eq. (S3), the self-energy on He,0 from Hph,1 is v0(ω − Hph,1 + ...)^{-1}v†0, and the self-energy on Hph,1 from He,1 is V1(ω − He,1 + ...)^{-1}V†1; Eq. (5b) states exactly this. Eq. (6) instead writes v†0(...)v0 and V†1(...)V1 in those denominators. Because Eq. (8) and the SI's Dyson terms (e.g., V0Ge,0v0Gph,1v†0Ge,0V†0 in the fourth-order table of Sec. S5) use the opposite ordering, Eq. (8) cannot be obtained from Eq. (6) as printed; the conversion in Sec. S6b silently switches to the correct ordering when it introduces X = v0Gph,1(I − V1Ge,1V†1Gph,1)^{-1}v†0Ge,0. The ordering is immaterial in the fully scalar 3-level model used in Fig. 2, which may explain the oversight, but the manuscript claims validity for arbitrary vibronic manifolds in which the blocks are non-commuting matrices. Please correct the dagger placement in Eq. (6), Eq. (S18), and Eq. (S22), or state the alternative convention explicitly and make Eq. (5b), the diagram rules, and Eq. (8) consistent with it.
  2. [Figs. 2a–2b; SI Sec. S6] No comparison is made between the 1/N-truncated Green's function and the exact DR_N. The exact result is readily available either by numerically solving the finite tridiagonal system defined by Eq. (5) (or the continued fraction once its ordering is fixed) or by exact diagonalization of Eq. (2) for small N. Such a comparison would serve two purposes: it would independently confirm that Eq. (8) is indeed the O(1/N) term of the exact photon Green's function — essential given the ordering inconsistency in Eq. (6) — and it would demonstrate that the 1/N expansion actually converges for the N values plotted in Fig. 2. I recommend adding a panel overlaying the exact and truncated absorption spectra for a few values of N (e.g., N=5–50) at fixed λ√N, including an error estimate for the truncated series.
  3. [Eq. (9); SI Sec. S5A, Eqs. (S18)–(S19)] The paper's title-level claim, that linear polaritonic spectra are governed by hidden nonlinear molecular susceptibilities, is carried by Eq. (9). The supporting argument equates the continued fraction (S18) with the nested matrix-geometric series (S19), but the SI only says that 'this result can be explicitly checked by the reader'; no induction, convergence statement, or justification of the order of non-commuting factors in the nested powers of Eq. (S19) is given. Furthermore, with Eq. (S18) exactly as printed, the two expressions are not equal as operator identities in the general multi-level case (see major comment 1 on the dagger ordering), so the equality cannot simply be taken on trust. I ask the authors to replace the sketch with a direct derivation of Eq. (9), e.g., an induction on the truncation level of the continued fraction, together with a statement of the sense in which the infinite sums converge.
minor comments (7)
  1. [SI Sec. S6c] In the opening sentence of Sec. S6c, 'To obtain dN,2(ω) ∝ N^{−k}' should read 'dN,2(ω) ∝ N^{−2}'.
  2. [Fig. 2 caption] The caption to Fig. 2 does not identify which colors in panel (b) correspond to the zeroth- and first-order spectra, and it is unclear whether the blue curve shows the O(N^{-2}) increment alone or the sum of all included orders; please clarify.
  3. [Main text, '1/N expansion'] The statement that 'we cannot provide an expression for dN,k when k ≥ 2' is stronger than what the Supplemental Material achieves, since Sec. S6c gives d^{(1)}_{N,2} in closed form; the wording should be changed to 'no complete closed-form expression for dN,k is obtained for k ≥ 2'.
  4. [Experimental considerations] The observability prediction assumes a homogeneous ensemble at zero temperature with Markovian losses and a single cavity mode; the text acknowledges the multimode-cavity blurring but not the possible suppression of the 1/N sidebands by static disorder or finite-temperature phonon populations. A brief quantitative estimate (or an explicit statement that these effects are beyond the present scope) would make the experimental claim more complete.
  5. [SI Eqs. (S17a)–(S17c)] The index placement in Eqs. (S17) is inconsistent (e.g., µ_{0g}^{me} versus µ_{me 0g}), which makes the susceptibility expressions difficult to parse; a uniform tensor notation consistent with Eq. (S.2) is recommended.
  6. [SI Sec. S5] The equivalence of the CUT-E diagram rules with the ket-only Dyson diagrams is asserted with the parenthetical promise of an induction proof that is not given; please either supply the induction or state the induction hypothesis explicitly.
  7. [Introduction / Ref. [23]] Since the block-tridiagonal (CUT-E) structure and the 1/N expansion were introduced in Ref. [23] by the same group, the authors should state explicitly which elements (the linear-spectroscopy application, the Raman-sideband prediction, and the susceptibility identification) are new in the present work.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the 1/N corrections are derived from the stated Hamiltonian via standard Green's-function algebra, not fitted or defined into existence.

full rationale

The paper's central result, Eq. (8) for dN,1, is obtained in Supplemental Sec. S6b by truncating the continued-fraction self-energy (Eq. S18) and expanding the resulting matrix geometric series; this is a derivation, not a restatement of the model. dN,0 reduces to the input linear susceptibility in the thermodynamic limit, which is a consistency check rather than a circular input. The self-citations to the authors' CUT-E work supply the block-tridiagonal representation and timescale language, but the recurrence relations (Eqs. 5a/5b) and the scaling Vn proportional to lambda*sqrt(N) and vn proportional to lambda are stated and used directly, so the citation is not load-bearing. Eq. (9)'s sum over 'irreducible' susceptibilities is an identification/rewriting of the self-energy expansion, and the Raman-sideband prediction rests on Eq. (8), not on that rewriting. No parameters are fitted to the predicted spectra. The printed dagger ordering in Eq. (6)/S18 is inconsistent with the Schur-complement derivation and with the ordering used in S6b; that is a correctness/reproducibility concern, not circularity, and does not affect the circularity score.

Assumptions & free parameters 0 free parameters · 5 assumptions · 1 invented entities

No free parameters fitted to data; the model inputs are lambda, kappa, gamma, omega_v, omega_ph and N. The central assumptions are the bosonic mapping, the block tridiagonal structure from the CUT-E framework, the single-mode Markovian input-output description, and the zero-temperature initial state. One conceptual entity is introduced: the classification of irreducible nonlinear susceptibilities.

assumptions (5)
  • domain assumption The molecular ensemble can be described by a bosonic Hamiltonian (Eq. 1) derived from the Holstein-Tavis-Cummings model under the assumption of a permutationally invariant initial state at zero temperature.
    Invoked in the Model section and SI S1; this mapping is established in prior work (Refs 17-22) and assumes bosonic statistics for molecular vibronic excitations.
  • domain assumption The Hamiltonian in the first excitation manifold has a block tridiagonal form (Eq. 2) with Vn scaling as lambda*sqrt(N) and vn scaling as lambda.
    Structural assumption from Ref 23 (CUT-E); underlies the timescale separation and all subsequent 1/N counting.
  • domain assumption Linear response observables are expressed via the retarded photon Green's function with input-output relations using Markovian loss rates kappa (cavity) and gamma (molecular dephasing).
    Standard input-output theory; stated in the Model section and SI S2, following Ref 16.
  • domain assumption The initial state at zero temperature has no photons and all molecules in the ground vibronic state.
    Ensures only ket-only diagrams contribute; stated in SI S2.
  • standard math The matrix continued fraction expansion and Schur complement formulas (SI S3) are valid for the non-Hermitian effective Hamiltonians with complex loss terms.
    Standard linear algebra; used without proof.
invented entities (1)
  • Irreducible nonlinear susceptibilities of the N-molecule ensemble
    purpose: To represent the cavity self-energy as a geometric sum over nonlinear response functions (Eq. 9) and to classify Feynman diagrams by whether they revisit the photonic block Hph,0.
    A new classification construct introduced in the Feynman pathways section; no direct experimental observable outside the paper, though the underlying susceptibilities are physical.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Hidden nonlinear optical susceptibilities in linear polaritonic spectra." pith.science (2026). https://pith.science/paper/3UO3I6SY

@misc{pith2026241109039,
  author       = {Pith},
  title        = {Pith review of: Hidden nonlinear optical susceptibilities in linear polaritonic spectra},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3UO3I6SY}},
  note         = {Machine review of arXiv:2411.09039}
}
abstract

Linear spectra of molecular polaritons formed by $N$ molecules coupled to a microcavity photon mode are usually well described by classical linear optics, raising the question of where the expected nonlinear effects in these strongly coupled systems are. In this work, we derive a general expression for the polaritonic linear spectra that reveal previously overlooked finite-size quantum corrections due to vacuum-mediated molecular Raman processes. Using a $1/N$ expansion, we demonstrate that these nonlinearities are suppressed in typical low-Q cavities due to an emergent timescale separation in polariton dynamics yet manifest in high-Q single-mode cavities where the photon loss is comparable to the single-molecule light-matter coupling strength.

Figures

Figures reproduced from arXiv: 2411.09039 by the authors.

Figure 1
Figure 1. FIG. 1. Strong coupling between a cavity photon and molec [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. (a) Linear absorption spectrum of an ensemble of 3-level systems (two ground-state vibrational levels and one excited [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

49 extracted references · 45 canonical work pages

  1. [1]

    D. M. Coles, N. Somaschi, P. Michetti, C. Clark, P. G. Lagoudakis, P. G. Savvidis, and D. G. Lidzey, “Polariton- mediated energy transfer between organic dyes in a 5 FIG. 3. Cavity-mediated nonlinear molecular processes : Ladder diagrams for an ensemble of identical three-level systems with two ground-state vibrational levels and one excited state, showin...

  2. [2]

    Energy transfer and correlations in cavity-embedded donor-acceptor configu- rations,

    M. Reitz, F. Mineo, and C. Genes, “Energy transfer and correlations in cavity-embedded donor-acceptor configu- rations,” Scientific Reports 8, 9050 (2018)

  3. [3]

    A path towards single molecule vi- brational strong coupling in a fabry–p´ erot microcavity,

    A. Koner, M. Du, S. Pannir-Sivajothi, R. H. Goldsmith, and J. Yuen-Zhou, “A path towards single molecule vi- brational strong coupling in a fabry–p´ erot microcavity,” Chemical Science 14, 7753–7761 (2023)

  4. [4]

    Modifying chemical landscapes by cou- pling to vacuum fields,

    J. A. Hutchison, T. Schwartz, C. Genet, E. Devaux, and T. W. Ebbesen, “Modifying chemical landscapes by cou- pling to vacuum fields,” Angewandte Chemie Interna- tional Edition 51, 1592–1596 (2012)

  5. [5]

    Molecular polaritons for controlling chemistry with quantum optics,

    F. Herrera and J. Owrutsky, “Molecular polaritons for controlling chemistry with quantum optics,” The Journal of Chemical Physics 152 (2020)

  6. [6]

    Room-temperature polariton lasing from cdse core-only nanoplatelets,

    F. Freire-Fern´ andez, N. G. Sinai, M. J. Hui Tan, S.-M. Park, E. R. Koessler, T. Krauss, P. Huo, and T. W. Odom, “Room-temperature polariton lasing from cdse core-only nanoplatelets,” ACS Nano 18, 15177–15184 (2024)

  7. [7]

    Polarized room-temperature polariton lasing in elliptical microcav- ities filled with fluorescent proteins,

    C. Bennenhei, M. Struve, S. Stephan, N. Kunte, V. N. Mitryakhin, F. Eilenberger, J. Ohmer, U. Fischer, M. Silies, C. Schneider, and M. Esmann, “Polarized room-temperature polariton lasing in elliptical microcav- ities filled with fluorescent proteins,” Opt. Mater. Ex- press 13, 2633–2643 (2023)

  8. [8]

    Po- lariton condensation in organic semiconductors,

    K. S. Daskalakis, S. A. Maier, and S. K´ ena-Cohen, “Po- lariton condensation in organic semiconductors,” Quan- tum Plasmonics , 151–163 (2017)

Show all 49 references
  1. [9]

    When do molecular polaritons behave like optical filters?

    K. Schwennicke, A. Koner, J. B. P´ erez-S´ anchez, W. Xiong, N. C. Giebink, M. L. Weichman, and J. Yuen- Zhou, “When do molecular polaritons behave like optical filters?” (2024), arXiv:2408.05036 [physics.chem-ph]

  2. [10]

    Vacuum rabi splitting as a feature of linear-dispersion theory: Analysis and experimental observations,

    Y. Zhu, D. J. Gauthier, S. Morin, Q. Wu, H. Carmichael, and T. Mossberg, “Vacuum rabi splitting as a feature of linear-dispersion theory: Analysis and experimental observations,” Physical Review Letters 64, 2499 (1990)

  3. [11]

    Polarization-dependent optical parame- ters of arbitrarily anisotropic homogeneous layered sys- tems,

    M. Schubert, “Polarization-dependent optical parame- ters of arbitrarily anisotropic homogeneous layered sys- tems,” Physical Review B 53, 4265 (1996)

  4. [12]

    Yariv and P

    A. Yariv and P. Yeh, Photonics: optical electronics in modern communications (Oxford University Press, 2007)

  5. [13]

    Transfer-matrix approach to determining the linear re- sponse of all-fiber networks of cavity-qed systems,

    N. N´ emet, D. White, S. Kato, S. Parkins, and T. Aoki, “Transfer-matrix approach to determining the linear re- sponse of all-fiber networks of cavity-qed systems,” Phys- ical Review Applied 13, 064010 (2020)

  6. [14]

    Excitonic spectral features in strongly coupled organic polaritons,

    J. A. ´Cwik, P. Kirton, S. De Liberato, and J. Keeling, “Excitonic spectral features in strongly coupled organic polaritons,” Physical Review A 93, 033840 (2016)

  7. [15]

    Con- tinuum model of strong light-matter coupling for molecu- lar polaritons,

    S. Gunasekaran, R. F. Pinard, and A. J. Musser, “Con- tinuum model of strong light-matter coupling for molecu- lar polaritons,” arXiv preprint arXiv:2308.08744 (2023)

  8. [16]

    Linear response of molec- ular polaritons,

    J. Yuen-Zhou and A. Koner, “Linear response of molec- ular polaritons,” The Journal of Chemical Physics 160, 154107 (2024)

  9. [17]

    Radiative pump- ing vs vibrational relaxation of molecular polaritons: a bosonic mapping approach,

    J. B. P´ erez-S´ anchez and J. Yuen-Zhou, “Radiative pump- ing vs vibrational relaxation of molecular polaritons: a bosonic mapping approach,” (2024), arXiv:2407.20594 [quant-ph]

  10. [18]

    Efficient and exact numerical approach for many multi-level systems in open system cqed,

    M. Gegg and M. Richter, “Efficient and exact numerical approach for many multi-level systems in open system cqed,” New J. Phys. 18, 043037 (2016). 6

  11. [19]

    Open quantum systems with local and col- lective incoherent processes: Efficient numerical simula- tions using permutational invariance,

    N. Shammah, S. Ahmed, N. Lambert, S. De Liberato, and F. Nori, “Open quantum systems with local and col- lective incoherent processes: Efficient numerical simula- tions using permutational invariance,” Phys. Rev. A 98, 063815 (2018)

  12. [20]

    Efficient linear scaling mapping for permu- tation symmetric fock spaces,

    M. A. Zeb, “Efficient linear scaling mapping for permu- tation symmetric fock spaces,” Comp. Phys. Commun. 276, 108347 (2022)

  13. [21]

    Permutational symmetry for identical multilevel systems: A second-quantized ap- proach,

    R. E. F. Silva and J. Feist, “Permutational symmetry for identical multilevel systems: A second-quantized ap- proach,” Phys. Rev. A 105, 043704 (2022)

  14. [22]

    Second quantization of open quantum systems in liouville space,

    V. Sukharnikov, S. Chuchurka, A. Benediktovitch, and N. Rohringer, “Second quantization of open quantum systems in liouville space,” Phys. Rev. A 107, 053707 (2023)

  15. [24]

    Theory of the Coupling of Electronic and Vibrational Excitations in Molecular Crystals and Helical Polymers,

    M. R. Philpott, “Theory of the Coupling of Electronic and Vibrational Excitations in Molecular Crystals and Helical Polymers,” J. Chem. Phys.55, 2039–2054 (1971)

  16. [25]

    Dark vibronic polaritons and the spectroscopy of organic microcavities,

    F. Herrera and F. C. Spano, “Dark vibronic polaritons and the spectroscopy of organic microcavities,” Phys. Rev. Lett. 118, 223601 (2017)

  17. [26]

    Simulating molecular polaritons in the collective regime using few-molecule models,

    J. B. P´ erez-S´ anchez, A. Koner, N. P. Stern, and J. Yuen- Zhou, “Simulating molecular polaritons in the collective regime using few-molecule models,” Proceedings of the National Academy of Sciences 120, e2219223120 (2023)

  18. [27]

    On the equivalence of differ- ent techniques for evaluating the green function for a semi-infinite system using a localized basis,

    J. Velev and W. Butler, “On the equivalence of differ- ent techniques for evaluating the green function for a semi-infinite system using a localized basis,” Journal of Physics: Condensed Matter 16, R637 (2004)

  19. [28]

    Matrix continued-fraction calculation of localization length,

    H. M. Pastawski, J. F. Weisz, and S. Albornoz, “Matrix continued-fraction calculation of localization length,” Phys. Rev. B 28, 6896–6903 (1983)

  20. [29]

    Continued frac- tion expansions in scattering theory and statistical non- equilibrium mechanics,

    P. H¨ anggi, F. R¨ osel, and D. Trautmann, “Continued frac- tion expansions in scattering theory and statistical non- equilibrium mechanics,” Zeitschrift f¨ ur Naturforschung A 33, 402–417 (1978)

  21. [30]

    Mukamel, Principles of nonlinear optical spectroscopy (Oxford University Press, 1995)

    S. Mukamel, Principles of nonlinear optical spectroscopy (Oxford University Press, 1995)

  22. [31]

    Exact states and spectra of vibrationally dressed polaritons,

    M. A. Zeb, P. G. Kirton, and J. Keeling, “Exact states and spectra of vibrationally dressed polaritons,” ACS Photonics 5, 249–257 (2018)

  23. [32]

    Frequency-dependent photoreactivity in disordered molecular polaritons,

    J. B. P´ erez-S´ anchez, F. Mellini, N. C. Giebink, and J. Yuen-Zhou, “Frequency-dependent photoreactivity in disordered molecular polaritons,” arXiv preprint arXiv:2308.03954 (2023)

  24. [33]

    In- ducing disallowed two-atom transitions with temporally entangled photons,

    A. Muthukrishnan, G. S. Agarwal, and M. O. Scully, “In- ducing disallowed two-atom transitions with temporally entangled photons,” Phys. Rev. Lett. 93, 093002 (2004)

  25. [34]

    Two- photon excitation sets limit to entangled photon pair gen- eration from quantum emitters,

    T. Seidelmann, C. Schimpf, T. K. Bracht, M. Cosacchi, A. Vagov, A. Rastelli, D. E. Reiter, and V. M. Axt, “Two- photon excitation sets limit to entangled photon pair gen- eration from quantum emitters,” Phys. Rev. Lett. 129, 193604 (2022)

  26. [35]

    Photon-exchange induces optical nonlinearities in har- monic systems,

    R. Glenn, K. Bennett, K. E. Dorfman, and S. Mukamel, “Photon-exchange induces optical nonlinearities in har- monic systems,” Journal of Physics B: Atomic, Molecular and Optical Physics 48, 065401 (2015)

  27. [36]

    Fetter and J

    A. Fetter and J. Walecka, Quantum Theory of Many- particle Systems , Dover Books on Physics (Dover Publi- cations, 2003)

  28. [37]

    Polariton localization and dispersion properties of disordered quantum emitters in multimode microcavities,

    G. Engelhardt and J. Cao, “Polariton localization and dispersion properties of disordered quantum emitters in multimode microcavities,” Phys. Rev. Lett. 130, 213602 (2023)

  29. [38]

    Dispersionless gaps and cavity modes in photonic crystals containing hyperbolic meta- materials,

    C.-h. Xue, Y. Ding, H.-t. Jiang, Y. Li, Z.-s. Wang, Y.- w. Zhang, and H. Chen, “Dispersionless gaps and cavity modes in photonic crystals containing hyperbolic meta- materials,” Phys. Rev. B 93, 125310 (2016)

  30. [39]

    Dispersion engineering and measurement of whispering gallery mode microresonator for kerr frequency comb generation,

    S. Fujii and T. Tanabe, “Dispersion engineering and measurement of whispering gallery mode microresonator for kerr frequency comb generation,” Nanophotonics 9, 1087–1104 (2020)

  31. [40]

    Intercavity polariton slows down dynamics in strongly coupled cavities,

    Y. A. Garc´ ıa Jomaso, B. Vargas, D. L. Dom´ ınguez, R. J. Armenta-Rico, H. E. Sauceda, C. L. Ordo˜ nez-Romero, H. A. Lara-Garc´ ıa, A. Camacho-Guardian, and G. Pir- ruccio, “Intercavity polariton slows down dynamics in strongly coupled cavities,” Nature Communications 15, 2915 (2024)

  32. [41]

    Bosonic condensation and disorder-induced localization in a flat band,

    F. Baboux, L. Ge, T. Jacqmin, M. Biondi, E. Galopin, A. Lemaˆ ıtre, L. Le Gratiet, I. Sagnes, S. Schmidt, H. E. T¨ ureci, A. Amo, and J. Bloch, “Bosonic condensation and disorder-induced localization in a flat band,” Phys. Rev. Lett. 116, 066402 (2016). Hidden nonlinear optica...

  33. [42]

    The (2m)th order term in DR N(ω) consists of all the 2 m step paths in the CUT-E diagram starting from and returning to Hph,0; the order here is calculated with respect to the light-matter interaction, V ,

  34. [43]

    It can be seen (and proven explicitly using mathematical induction) that the CUT-E diagram with the aforementioned 3 a b c d e Figure S1: The ket-only DS-FDs for the a

    Every time we jump from one box to another in the CUT-E diagram, we encounter a Raman process mediated by a v(†) k which are penalized by a factor of 1 / √ N. It can be seen (and proven explicitly using mathematical induction) that the CUT-E diagram with the aforementioned 3 a...

  35. [44]

    We have DR,(0)(ω) = Gph,0(ω) = 1 ω−ωph+iκ/2

    Zeroth order term: 0 step paths. We have DR,(0)(ω) = Gph,0(ω) = 1 ω−ωph+iκ/2

  36. [45]

    We have just one of them: Hph,0 1 ← → 2 He,0 (the numbers on the top (bottom) represent the counting of the step forward (backward))

    Second order term : We want to count the 2 step paths. We have just one of them: Hph,0 1 ← → 2 He,0 (the numbers on the top (bottom) represent the counting of the step forward (backward)). Thus, second-order photon Green’s function, DR,(2)(ω), is then given as DR,(2)(ω) = Gph,...

  37. [46]

    Fourth order term: All the 4 step paths are presented in the table below . The fourth-order photon Green’s function, DR,(4)(ω), is then given as the sum of all the Dyson series terms, DR,(4) N (ω) Pathways in CUT-E diagram Dyson series term Hph,0 1,3 ←→ 2,4 He,0 = 2· (Hph,0↔ H...

  38. [47]

    We can have five of them presented in the table below

    Sixth order term: We want all the 6-step paths in the diagram. We can have five of them presented in the table below. Note that DR,(6) N (ω) is the sum of all the Dyson series terms. 4 DR,(6) N (ω) Pathways in CUT-E diagram Dyson series term Hph,0 1,3,5 ←−→ 2,4,6 He,0 = 3· (Hp...

  39. [48]

    Proceeding this way, we get a recursive algorithm, where the self-energy at the kth step in the matrix depends on the ( k + 1)th step

    (S.14) 6 Again Σph,1 = v0 ( ω−Hph,1+iκ/2−Σe,1 )−1 v† 0 is obtained in terms of the self-energy of the next step, Σe,1. Proceeding this way, we get a recursive algorithm, where the self-energy at the kth step in the matrix depends on the ( k + 1)th step. Thus we have, Σe,k = Vk...

  40. [49]

    [vN· Gph,N v† N Ge,N−1]nph,N

    (S.18) On the other hand, we have the result that the sum over the irreducible nonlinear susceptibilities computed diagram- matically using the CUT-E diagram and the rules, given as a nested summation, ∞∑ l=0 (ωph 2 )l χ(2l+1) N ({ω,ω−ωph}l,ω ) = − ∞∑ ne,0=0 ∞∑ nph,1=0 ··· ∞∑ ...

  41. [50]

    Cut-e as a 1/n expansion for multiscale molecular polariton dynamics,

    (S.19) This result can be explicitly checked by the reader by isolating the pathways that start and end at Hph,0, and do not involve Hph,0 in between, in other words, they must have V0 and V† 0 only in the left and right, respectively (see Sec. S5). We notice that series const...

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.