pith. sign in

arxiv: 2512.02623 · v1 · submitted 2025-12-02 · 🪐 quant-ph · physics.atom-ph· physics.optics

High-harmonic generation from two weakly coupled molecules: a simple tight-binding model

Lina Bielke , Samuel Sch\"opa , Falk-Erik Wiechmann , Franziska Fennel , Dieter Bauer This is my paper

Pith reviewed 2026-05-17 02:41 UTC · model grok-4.3

classification 🪐 quant-ph physics.atom-phphysics.optics
keywords high-harmonic generationtight-binding modelmolecular dimerpolarization dependenceintermolecular couplingadiabatic statesharmonic yield
0
0 comments X

The pith

In a model of two weakly coupled molecules, lower harmonic orders maximize when laser polarization aligns with molecular axes while higher orders maximize along the intermolecular axis, with the switch depending on coupling strength.

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

This paper examines the effect of weak intermolecular coupling on high-harmonic generation using a simple two-dimensional tight-binding model of a molecular dimer. It shows that the laser polarization direction that produces the strongest yield changes with harmonic order: lower orders are strongest for polarization along the individual molecular axes, while higher orders are strongest for polarization along the line connecting the two molecules. The order at which this preference reverses depends on the strength of the coupling between the molecules. An adiabatic analysis demonstrates that this directional flip is already present in the states that follow the driving field. A reader would care because the result offers a concrete way to connect measurable harmonic spectra to the weak interactions that matter in molecular crystals and dimers.

Core claim

Using a two-dimensional tight-binding model for a molecular dimer with weak intermolecular coupling, the intensities of lower harmonic orders tend to maximize for a laser polarization direction aligning with the molecular axes, whereas higher harmonic orders rather show the strongest yield for a polarization direction along the intermolecular axis. The harmonic order at which the maximum flips from the molecular to the intermolecular direction strongly depends on the intermolecular coupling strength. A detailed adiabatic analysis shows that the flipping of the maximum yield towards the intermolecular direction is already contained qualitatively in the adiabatically following states.

What carries the argument

The two-dimensional tight-binding model of the molecular dimer with tunable weak intermolecular coupling, solved by propagating the time-dependent Schrödinger equation for varying laser polarization angles.

Load-bearing premise

The two-dimensional tight-binding model with weak intermolecular coupling accurately captures the essential physics of real molecular dimers or crystals for the polarization-dependent harmonic yields under study.

What would settle it

An experiment on an actual molecular dimer or crystal that measures no change in the optimal polarization direction with harmonic order, or finds that the switchover order does not vary with intermolecular distance or coupling strength, would falsify the model's predictive power for these yields.

Figures

Figures reproduced from arXiv: 2512.02623 by Dieter Bauer, Falk-Erik Wiechmann, Franziska Fennel, Lina Bielke, Samuel Sch\"opa.

Figure 1
Figure 1. Figure 1: FIG. 1. (a) Sketch of the model system. The blue circles [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Normalized harmonic intensities as a function of the [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. High-harmonic spectra for laser polarization direc [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (a) Normalized harmonic intensities as a function of the laser polarization angle for the reduced intermolecular coupling [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Harmonic intensities as a function of the intermolecu [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
Figure 5
Figure 5. Figure 5: These different dependencies explain the behav [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Exact harmonic spectrum along with the results [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Normalized harmonic intensities as a function of the [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Intensities of the different harmonic orders for the [PITH_FULL_IMAGE:figures/full_fig_p008_9.png] view at source ↗
read the original abstract

The generation of high harmonics is a strongly nonlinear effect that allows to probe properties of the target and to study electron dynamics in matter. It has been investigated in many different kinds of targets, including molecular gases, liquids and solids. Recently, high-harmonic generation was studied in organic molecular crystals by Wiechmann et al. [Nat. Commun. 16, 9890 (2025)]. It was found that the laser-polarization-dependent harmonic yield is sensitive to the weak couplings between nearest- and next-nearest-neighbor molecules. In this paper, the impact of the laser polarization angle and the intermolecular interaction on the harmonic yield is examined in detail using a simple but insightful two-dimensional tight-binding system that models a molecular dimer, i.e. two weakly coupled molecules. We find that the intensities of lower harmonic orders tend to maximize for a laser polarization direction aligning with the molecular axes, whereas higher harmonic orders rather show the strongest yield for a polarization direction along the intermolecular axis. We further demonstrate that the harmonic order at which the maximum flips from the molecular to the intermolecular direction strongly depends on the intermolecular coupling strength. To gain a deeper insight into the origins of the findings, we include a detailed adiabatic analysis, showing that the flipping of the maximum yield towards the intermolecular direction is already contained qualitatively in the adiabatically following states.

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 introduces a minimal two-dimensional tight-binding model of a molecular dimer with tunable weak intermolecular coupling to study laser-polarization-dependent high-harmonic generation. Numerical propagation of the time-dependent Schrödinger equation shows that lower-order harmonics maximize for polarization aligned with the molecular axes while higher-order harmonics maximize along the intermolecular axis; the crossover harmonic order depends on the coupling strength. An adiabatic analysis of the instantaneous eigenstates is used to argue that the directional flip is already encoded in the field-dressed states.

Significance. If the reported trends are robust, the work supplies a transparent, reproducible minimal model that isolates the effect of weak intermolecular coupling on HHG polarization dependence, directly motivated by recent experiments on organic molecular crystals. The adiabatic analysis is a clear strength, providing a parameter-free qualitative explanation without requiring full time-dependent simulations for every case. The approach is well-suited for gaining mechanistic insight and could guide interpretation of more complex solid-state HHG data.

major comments (2)
  1. [§3] §3 (Numerical results and figures): the central claim that the crossover harmonic order 'strongly depends' on intermolecular coupling strength requires explicit quantification. The manuscript should report the specific coupling values examined, the precise definition of the crossover (e.g., angle of maximum yield crossing 45°), and an additional panel or table showing crossover order versus coupling strength to substantiate the dependence.
  2. [Adiabatic analysis] Adiabatic analysis section: while the instantaneous eigenstates qualitatively reproduce the directional flip, the manuscript does not provide a quantitative metric (e.g., overlap of yield maxima or relative error in harmonic intensities) comparing the adiabatic prediction to the full TDSE results. This comparison is load-bearing for the assertion that the effect is 'already contained qualitatively in the adiabatically following states'.
minor comments (3)
  1. [Figures] Figure captions (e.g., Figs. 2–4): explicitly state the zero of the polarization angle (molecular axis or intermolecular axis) and the normalization convention for the harmonic yields to improve immediate readability.
  2. [Introduction] Introduction: expand the one-sentence reference to Wiechmann et al. with a brief statement of their key experimental observation on nearest- and next-nearest-neighbor couplings to strengthen the motivation.
  3. [Methods] Methods or appendix: include a short statement on the numerical integrator, time step, grid size, and convergence tests used for the TDSE propagation and HHG spectrum extraction.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment of our work and for the constructive suggestions that will help strengthen the manuscript. We address each major comment below and outline the revisions we will make.

read point-by-point responses

  1. Referee: [§3] §3 (Numerical results and figures): the central claim that the crossover harmonic order 'strongly depends' on intermolecular coupling strength requires explicit quantification. The manuscript should report the specific coupling values examined, the precise definition of the crossover (e.g., angle of maximum yield crossing 45°), and an additional panel or table showing crossover order versus coupling strength to substantiate the dependence.

    Authors: We agree that explicit quantification will improve the clarity of the central claim. In the revised manuscript we will list the specific intermolecular coupling strengths examined (e.g., t⊥ = 0.01, 0.05, 0.10, 0.20 in units of the intramolecular hopping), define the crossover order as the harmonic number at which the polarization angle of maximum yield crosses 45° relative to the molecular axes, and add a new panel (or compact table) to Figure 3 that plots the crossover order versus coupling strength. These additions will directly substantiate the reported dependence. revision: yes

  2. Referee: [Adiabatic analysis] Adiabatic analysis section: while the instantaneous eigenstates qualitatively reproduce the directional flip, the manuscript does not provide a quantitative metric (e.g., overlap of yield maxima or relative error in harmonic intensities) comparing the adiabatic prediction to the full TDSE results. This comparison is load-bearing for the assertion that the effect is 'already contained qualitatively in the adiabatically following states'.

    Authors: We acknowledge that a quantitative metric would strengthen the support for our assertion. In the revised version we will include a direct comparison: for representative coupling strengths and harmonic orders we will report the angular difference (in degrees) between the yield-maximizing polarization direction obtained from the adiabatic eigenstates and from the full TDSE propagation, together with the relative error in the extracted harmonic intensities. This will provide a concrete measure of the agreement while preserving the primarily qualitative character of the adiabatic analysis. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper introduces an explicit two-dimensional tight-binding Hamiltonian for a molecular dimer with tunable intermolecular coupling, then obtains the reported polarization-dependent harmonic yields and crossover order directly from numerical propagation and an adiabatic analysis of the instantaneous eigenstates. These steps constitute independent computation within the model's stated assumptions rather than any reduction of outputs to prior fitted quantities, self-definitions, or load-bearing self-citations. The cited Wiechmann et al. work supplies experimental motivation but is not invoked to justify uniqueness or to close any derivation loop. The central claims therefore remain self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The model introduces a tunable intermolecular hopping (coupling) parameter whose value controls the flip order; the tight-binding approximation itself is an ad-hoc simplification for the dimer geometry.

free parameters (1)
  • intermolecular coupling strength
    Varied parametrically to demonstrate dependence of the flip harmonic order; no specific fitted value is given in the abstract.
axioms (1)
  • domain assumption Electrons in the dimer can be described by a two-dimensional tight-binding Hamiltonian with nearest-neighbor hoppings along molecular axes and a weaker inter-molecular hopping.
    Invoked to reduce the many-electron molecular problem to a simple lattice model.

pith-pipeline@v0.9.0 · 5557 in / 1109 out tokens · 47134 ms · 2026-05-17T02:41:26.776508+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

41 extracted references · 41 canonical work pages

  1. [1]

    The higher harmonic orders (n≥7), in contrast, show the highest yield for polarization angles close toα inter

    maximize when the laser is polarized parallel to the molecular axes, as marked by the solid black line. The higher harmonic orders (n≥7), in contrast, show the highest yield for polarization angles close toα inter. The precise angles of the maxima vary slightly aroundα inter, as they depend on the details of the probability density distribution and the tr...

    work page

  2. [2]

    Also the narrow- ing of the lobes with increasing harmonic order, which was seen in the measurements and theoretical calcula- tions of [25], is reproduced by our dimer model

    for a periodic tight-binding model of an organic molecular crystal, where the same effect of the lobes flip- ping from the molecular to the intermolecular direction at some harmonic order was observed. Also the narrow- ing of the lobes with increasing harmonic order, which was seen in the measurements and theoretical calcula- tions of [25], is reproduced ...

    work page

  3. [3]

    they follow the laser field adiabatically

    The electrons remain in the instantaneous eigen- states, i.e. they follow the laser field adiabatically

    work page

  4. [4]

    Only terms describing transitions between initially populated and initially unpopulated instantaneous eigenstates are considered. In condensed matter HHG, one typically distinguishes between intraband contributions to the harmonic yield (due to electronic movements within the energy bands) and interband contributions (due to transitions between energy ban...

    work page

  5. [5]

    Antoine, A

    P. Antoine, A. L’Huillier, and M. Lewenstein, Attosecond pulse trains using high–order harmonics, Phys. Rev. Lett. 77, 1234 (1996)

    work page 1996

  6. [6]

    P. M. Paul, E. S. Toma, P. Breger, G. Mullot, F. Aug´ e, P. Balcou, H. G. Muller, and P. Agostini, Observation of a train of attosecond pulses from high harmonic generation, Science292, 1689 (2001)

    work page 2001

  7. [7]

    Goulielmakis, M

    E. Goulielmakis, M. Schultze, M. Hofstetter, V. S. Yakovlev, J. Gagnon, M. Uiberacker, A. L. Aquila, E. M. Gullikson, D. T. Attwood, R. Kienberger, F. Krausz, and U. Kleineberg, Single-Cycle nonlinear optics, Science 320, 1614 (2008)

    work page 2008

  8. [8]

    Ghimire, A

    S. Ghimire, A. D. DiChiara, E. Sistrunk, P. Agostini, L. F. DiMauro, and D. A. Reis, Observation of high-order harmonic generation in a bulk crystal, Nat. Phys.7, 138 (2011)

    work page 2011

  9. [9]

    T. T. Luu, M. Garg, S. Y. Kruchinin, A. Moulet, M. T. Hassan, and E. Goulielmakis, Extreme ultraviolet high- harmonic spectroscopy of solids, Nature521, 498 (2015)

    work page 2015

  10. [10]

    Vampa, T

    G. Vampa, T. J. Hammond, N. Thir´ e, B. E. Schmidt, F. L´ egar´ e, C. R. McDonald, T. Brabec, D. D. Klug, and P. B. Corkum, All-optical reconstruction of crystal band structure, Phys. Rev. Lett.115, 193603 (2015)

    work page 2015

  11. [11]

    Lakhotia, H

    H. Lakhotia, H. Y. Kim, M. Zhan, S. Hu, S. Meng, and E. Goulielmakis, Laser picoscopy of valence electrons in solids, Nature583, 55 (2020)

    work page 2020

  12. [12]

    McPherson, G

    A. McPherson, G. Gibson, H. Jara, U. Johann, T. S. Luk, I. A. McIntyre, K. Boyer, and C. K. Rhodes, Studies of multiphoton production of vacuum-ultraviolet radiation in the rare gases, J. Opt. Soc. Am. B4, 595 (1987)

    work page 1987

  13. [13]

    Ferray, A

    M. Ferray, A. L’Huillier, X. F. Li, L. A. Lompre, G. Main- fray, and C. Manus, Multiple-harmonic conversion of 1064 nm radiation in rare gases, J. Phys. B21, L31 (1988)

    work page 1988

  14. [14]

    J. L. Krause, K. J. Schafer, and K. C. Kulander, High- order harmonic generation from atoms and ions in the high intensity regime, Phys. Rev. Lett.68, 3535 (1992)

    work page 1992

  15. [15]

    J. J. Macklin, J. D. Kmetec, and C. L. Gordon, High- order harmonic generation using intense femtosecond pulses, Phys. Rev. Lett.70, 766 (1993)

    work page 1993

  16. [16]

    Lyng˚ a, A

    C. Lyng˚ a, A. L’Huillier, and C.-G. Wahlstr¨ om, High- order harmonic generation in molecular gases, J. Phys. B29, 3293 (1996)

    work page 1996

  17. [17]

    M. Lein, N. Hay, R. Velotta, J. P. Marangos, and P. L. Knight, Interference effects in high-order harmonic gen- eration with molecules, Phys. Rev. A66, 023805 (2002)

    work page 2002

  18. [18]

    Smirnova, Y

    O. Smirnova, Y. Mairesse, S. Patchkovskii, N. Dudovich, D. Villeneuve, P. Corkum, and M. Y. Ivanov, High harmonic interferometry of multi-electron dynamics in molecules, Nature460, 972 (2009)

    work page 2009

  19. [19]

    Tisch, T

    J. Tisch, T. Ditmire, D. Fraser, N. Hay, M. Mason, E. Springate, J. Marangos, and M. Hutchinson, Inves- tigation of high-harmonic generation from xenon atom clusters, J. Phys. B30, L709 (1997)

    work page 1997

  20. [20]

    Peschel, M

    U. Peschel, M. Th¨ ummler, T. Lettau, S. Gr¨ afe, and K. Busch, Two-particle tight-binding description of higher-harmonic generation in semiconductor nanostruc- tures, Phys. Rev. B106, 245307 (2022)

    work page 2022

  21. [21]

    T. T. Luu, Z. Yin, A. Jain, T. Gaumnitz, Y. Pertot, J. Ma, and H. J. W¨ orner, Extreme–ultraviolet high– harmonic generation in liquids, Nat. Commun.9, 3723 (2018)

    work page 2018

  22. [22]

    Neufeld, Z

    O. Neufeld, Z. Nourbakhsh, N. Tancogne-Dejean, and A. Rubio, Ab initio cluster approach for high harmonic generation in liquids, J. Chem. Theory and Comput.18, 4117 (2022). 10

    work page 2022

  23. [23]

    Y. S. You, D. A. Reis, and S. Ghimire, Anisotropic high- harmonic generation in bulk crystals, Nat. Phys.13, 345 (2017)

    work page 2017

  24. [24]

    Tancogne-Dejean and A

    N. Tancogne-Dejean and A. Rubio, Atomic-like high- harmonic generation from two-dimensional materials, Sci. Adv.4, eaao5207 (2018)

    work page 2018

  25. [25]

    J¨ urß and D

    H. J¨ urß and D. Bauer, Edge-state influence on high- order harmonic generation in topological nanoribbons, Eur. Phys. J. D75, 190 (2021)

    work page 2021

  26. [26]

    Zhang, L

    Y. Zhang, L. Li, J. Li, T. Huang, P. Lan, and P. Lu, Ori- entation dependence of high-order harmonic generation in graphene, Phys. Rev. A104, 033110 (2021)

    work page 2021

  27. [27]

    C. P. Schmid, L. Weigl, P. Gr¨ ossing, V. Junk, C. Gorini, S. Schlauderer, S. Ito, M. Meierhofer, N. Hofmann, D. Afanasiev, J. Crewse, K. A. Kokh, O. E. Tereshchenko, J. G¨ udde, F. Evers, J. Wilhelm, K. Richter, U. H¨ ofer, and R. Huber, Tunable non-integer high-harmonic generation in a topological insulator, Na- ture593, 385 (2021)

    work page 2021

  28. [28]

    Korobenko, S

    A. Korobenko, S. Saha, A. T. K. Godfrey, M. Gertsvolf, A. Y. Naumov, D. M. Villeneuve, A. Boltasseva, V. M. Shalaev, and P. B. Corkum, High-harmonic generation in metallic titanium nitride, Nat. Commun.12, 4981 (2021)

    work page 2021

  29. [29]

    Wiechmann, S

    F.-E. Wiechmann, S. Sch¨ opa, L. Bielke, S. Rindelhardt, S. Patchkovskii, F. Morales, M. Richter, D. Bauer, and F. Fennel, High-order harmonic generation in an organic molecular crystal, Nat. Commun.16, 9890 (2025)

    work page 2025

  30. [30]

    Reese and Z

    C. Reese and Z. Bao, Organic single-crystal field-effect transistors, Mater. Today10, 20 (2007)

    work page 2007

  31. [31]

    Hasegawa and J

    T. Hasegawa and J. Takeya, Organic field-effect transis- tors using single crystals, Sci. Technol. Adv. Mater.10, 024314 (2009)

    work page 2009

  32. [32]

    Y.-Y. Lin, D. Gundlach, S. Nelson, and T. Jackson, Stacked pentacene layer organic thin-film transistors with improved characteristics, IEEE Electron Device Lett.18, 606 (1997)

    work page 1997

  33. [33]

    M. W. B. Wilson, A. Rao, B. Ehrler, and R. H. Friend, Singlet exciton fission in polycrystalline pentacene: From photophysics toward devices, Acc. Chem. Res.46, 1330 (2013)

    work page 2013

  34. [34]

    D. N. Congreve, J. Lee, N. J. Thompson, E. Hontz, S. R. Yost, P. D. Reusswig, M. E. Bahlke, S. Reineke, T. Van Voorhis, and M. A. Baldo, External quantum efficiency above 100% in a singlet-exciton-fission-based organic photovoltaic cell, Science340, 334 (2013)

    work page 2013

  35. [35]

    D. V. Skobel’tsyn, The oriented gas model and its ap- plication to molecular crystals, inVolume 25: Optical Methods of Investigating Solid Bodies, edited by D. V. Skobel’tsyn (Springer US, Boston, MA, 1965) pp. 44–66

    work page 1965

  36. [36]

    Schwoerer and H

    M. Schwoerer and H. C. Wolf,Organic Molecular Solids (Wiley, 2006)

    work page 2006

  37. [37]

    Sir W. H. Bragg, The structure of organic crystals, Proc. Phys. Soc. London34, 33 (1921)

    work page 1921

  38. [38]

    B. E. Douglas and S.-M. Ho,Structure and Chemistry of Crystalline Solids(Springer New York, NY, 2006)

    work page 2006

  39. [39]

    R. B. Campbell, J. M. Robertson, and J. Trotter, The crystal structure of hexacene, and a revision of the crys- tallographic data for tetracene, Acta Crystallogr.15, 289 (1962)

    work page 1962

  40. [40]

    J. E. Anthony, The larger acenes: Versatile organic semi- conductors, Angew. Chem. Int. Ed.47, 452 (2008)

    work page 2008

  41. [41]

    Born and V

    M. Born and V. Fock, Beweis des Adiabatensatzes, Zeitschrift f¨ ur Physik51, 165 (1928)

    work page 1928