Excitation Energy Transfer in Nanohybrid System of Organic Molecule and Inorganic Transition Metal Dichalcogenides Nanoflake

Technical Summary: Excitation Energy Transfer in an Organic/Inorganic Nanohybrid System

Problem Statement
Two-dimensional transition-metal dichalcogenides (TMDCs), such as MoS$_2$, possess unique electronic and optical properties but suffer from narrow absorption bands and high production costs. To enhance their functionality, they are often integrated with organic materials into heterostructures. While surface charge transfer mechanisms in these systems are relatively well-understood, the role of long-range Coulomb interactions—specifically Excitation Energy Transfer (EET) or Förster Resonance Energy Transfer (FRET)—remains a critical area of investigation, particularly in the absence of wavefunction overlap or energy level alignment. The specific challenge addressed in this work is the theoretical characterization of EET rates in a finite-sized MoS$_2$ nanoflake coupled with a single organic molecule (para-sexiphenyl, 6P), accounting for edge effects, excitonic states, and geometric dependencies without relying on phenomenological fitting parameters.

Methodology
The authors employ a multi-scale theoretical framework combining tight-binding modeling with configuration interaction and Fermi's Golden Rule:

• Electronic Structure of MoS$_2$: The electronic structure of finite MoS$_2$ nanoflakes (side lengths of 60, 70, and 74 Å) is modeled using an 11-band tight-binding Hamiltonian. This model utilizes an orbital basis of five $d$-orbitals from Molybdenum and six $p$-orbitals from Sulfur. To address the unphysical in-gap edge states inherent to finite systems, the dangling bonds at the sulfur edges are passivated with Hydrogen atoms, restoring a well-defined bandgap.
• Exciton Construction: Excitonic states are constructed within a configuration-interaction (CI) scheme. However, given that the 6P molecule's excitation energy (4 eV) is significantly higher than the MoS$_2$ bandgap (1.84 eV), the authors approximate the relevant high-energy states as uncorrelated electron-hole pairs rather than strongly bound Wannier-Mott excitons. This approximation reduces computational cost while retaining physical relevance for the high-energy window.
• EET Rate Calculation: The EET rates ($k_{mol \to sem}$ and $k_{sem \to mol}$) are evaluated using Fermi's Golden Rule. The coupling term ($V_{DA}$) is derived explicitly from Coulomb integrals between atomic-centered transition charges of the 6P molecule and the transition charge densities of the MoS$_2$ electron-hole pairs. The calculation incorporates thermal broadening and spectral overlap between the molecular excitation and the semiconductor's density of states.
• System Configuration: The 6P molecule is modeled as lying flat on a plane above the MoS$_2$ surface. The study varies the vertical separation (2–16 Å) and the lateral position of the molecule relative to the nanoflake center.

Key Contributions and Results

• Edge Passivation Validation: The study demonstrates that H-passivation effectively removes edge-localized states, resulting in a clear bandgap (~1.9 eV) and charge densities that resemble those of an infinite 2D sheet. Without passivation, edge states dominate the low-energy spectrum.
• Dominance of Molecule-to-Nanoflake Transfer: The results indicate that energy transfer from the 6P molecule to the MoS$_2$ nanoflake ($k_{mol \to sem}$) is the dominant process, exceeding the reverse transfer rate ($k_{sem \to mol}$) by approximately five orders of magnitude.
• Distance Dependence: The EET rate is highly sensitive to the vertical separation. At the minimum separation of 2 Å, the transfer rate is approximately $8.89 \times 10^{14}$ s$^{-1}$, corresponding to an excited-state lifetime of ~1 fs. As the distance increases to 16 Å, the rate drops to $0.41 \times 10^{14}$ s$^{-1}$. This decay is driven by the weakening of the Coulomb coupling, which falls from a maximum of ~140 meV at 2 Å to lower values at greater distances.
• Size and Position Effects:
  • Size: Larger nanoflakes (e.g., 74 Å vs. 60 Å) exhibit stronger Coulomb couplings and higher EET rates.
  • Lateral Position: The EET rate is non-uniform across the nanoflake surface. Rates are maximized when the molecule is positioned within the central region (approximately $\pm 35$ Å from the center) and decrease radically as the molecule moves toward the periphery. This spatial variation reflects the symmetry and distribution of the electron-hole states within the nanoflake.
• Coupling Strength: Individual transfer couplings for the 70 Å system are found to be less than 10 meV in the relevant energy window, confirming that the process is driven by weak, long-range Coulomb interactions rather than wavefunction overlap.

Significance and Claims
The paper claims to provide a quantitative theoretical description of non-contact energy transfer in a model organic-inorganic heterostructure. By combining an efficient tight-binding model for large-scale nanoflakes with explicit Coulomb integral calculations, the authors avoid the need for fitted parameters or periodic boundary conditions. The framework successfully highlights the critical roles of Coulomb coupling, spectral overlap, and interfacial geometry in determining EET efficiency. The authors note that their theoretical approach aligns with experimental trends observed in other TMDC-molecule hybrid systems and suggest that this framework could be extended to study few-layer MoS$_2$ systems, where properties like indirect bandgaps and interlayer excitons may introduce new complexities. The work establishes a baseline for understanding how nanoflake size and molecular positioning influence energy transfer in next-generation nanoscale optoelectronic devices.