Many-Body effects beyond excitons in second-harmonic generation of monolayer MoS 2

arXiv:2602.22463 · cond-mat.mtrl-sci, cond-mat.mes-hall · Submitted 2026-02-25 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Quantum Radio. Generated commentary on the latest quantum physics and condensed matter papers.

Kai: We're looking at the paper: "Many-Body effects beyond excitons in second-harmonic generation of monolayer MoS 2".

Mira: Here's the gist of the paper. The gist: The inclusion of three-particle trionic correlations through a static long-range quadratic exchange-correlation kernel achieves significantly improved quantitative agreement with experiment for second-harmonic generation in monolayer MoS2, highlighting the role of many-body interactions beyond the excitonic level in accurately describing second-order optical responses.

Kai: Theoretical Framework The study employs a combination of many-body perturbation theory and time-dependent current-density functional theory within an ab initio framework in the optical limit.

Mira: The theoretical description utilizes tensorial many-body description of second-harmonic generation, employing time-dependent current-density functional theory (TDCDFT) to represent the tensorial extension of TDDFT.

Kai: The response functions are formally equivalent to contracting the correlation parts of the two- and three-particle Green’s functions, establishing a correspondence between two- and three-body interactions and the linear fxc and quadratic gxc xc kernels of TDDFT.

Title and authors: Kai: So let's look at the title again, "Many-Body effects beyond excitons in second-harmonic generation of monolayer MoS two <ref:2602.22463#pg1,Many-Body effects beyond excitons in second-harmonic generation of monolayer>." It immediately tells you they're looking past the usual two-particle picture for this material.

Mira: I think that’s the crucial part because it sets up exactly what they are trying to address: those missing many-body interactions that standard models usually ignore.

Lev: From an engineering standpoint, if we were building a device based on this material, knowing the two-particle model is insufficient means our initial design assumptions about the nonlinear response might be wrong.

Kai: Exactly. The authors are showing that ignoring these higher-order effects leads to underestimation of the SHG intensity by almost a factor of two when compared to experimental values.

Mira: They're not just saying "it's complicated"; they’re providing a quantitative way to correct that deficiency using these three-particle correlations.

Lev: So, what’s the immediate implication for us who are trying to model this stuff for real hardware? It means we can't stop at the excitonic level description when dealing with nonlinear optical phenomena in these 2D materials <ref:2602.22463#pg1>.

Kai: That’s right, and it points toward a much more complete picture of how these systems respond to light beyond the simplest approximations.

Mira: We need to look at this paper because it shows that for accurate predictions of second-order responses, those three-body correlations aren't just small corrections; they are significant contributors.

Lev: So we're looking at a way to fix a quantitative gap in modeling nonlinear optical materials using these specific many-body terms.

The paper's summary: Kai: Let’s look at the actual summary of "Many-Body effects beyond excitons in second-harmonic generation of monolayer MoS two <ref:2602.22463#pg1,Many-Body effects beyond excitons in second-harmonic generation of monolayer>." It boils down to this: they use many-body perturbation theory combined with TDCDFT.

Mira: They use this ab initio framework to calculate the macroscopic SHG susceptibility tensor, which is how you figure out the light response of the material.

Lev: And they show that including two-particle excitonic effects alone yields an intensity of about zero point one four nm two/V, but we know the experiment is around zero point zero eight nm two/V.

Kai: That’s almost double the experimental value, which confirms that the two-body description is significantly underestimating the response because it misses those higher-order effects.

Mira: The paper then shows that adding that second term in their equation—the three-particle level effects—which are themselves screened by a linear kernel—that brings it right in line with experiment.

Lev: So, what’s the big number they found? What specific value did they find where things matched up best?

Kai: They found that for values of Bxx yLRC between negative four and zero atomic units, the impact on the maximum SHG response is substantial, reducing it by nearly a factor of two.

Mira: And they pinpointed an optimal parameter value, BLRC equals negative three point five atomic units, which matches the experimental peak intensity exactly.

Lev: That’s a lot of detail packed into this paper. It shows that these three-particle correlations are comparable in size to the two-particle contributions themselves.

The paper's improvements: Kai: The main improvement suggested by the authors is precisely this inclusion of trionic effects via a static long-range quadratic xc kernel to fix the quantitative agreement.

Mira: They argue that neglecting these three-particle trionic contributions is a significant limitation in state-of-the-art TDDFT approaches, which usually start from independent-particle responses.

Lev: So, what’s the practical change this means for us when we think about running these calculations on actual hardware? It means the overhead of including these three-body terms needs to be factored into our simulation costs.

Kai: It suggests that if we want quantitative predictions for nonlinear optical responses like SHG, we have to move beyond just two-particle excitonic states in our basis.

Mira: They propose a framework where you can model the frequency dependence of the macroscopic SHG susceptibility tensor using a Dyson like equation that includes both the two-particle effects and these three-particle effects explicitly.

Lev: That kind of explicit inclusion would definitely increase the complexity of whatever error correction or simulation we have to do, but it should yield much better results if implemented correctly.

Kai: They also suggest extracting material-dependent parameters like ALRC and BLRC from the BSE results to parameterize the strength of those long-range contributions in both linear and quadratic kernels.

Mira: So, you can use those calculated tensors to model how material properties dictate the strength of these many-body effects in the optical limit.

Conclusion: Kai: So we wrap up with what this paper really means for us. The main point is that many-body interactions beyond the usual two-particle excitonic level are necessary for a quantitative description of second-harmonic generation in monolayer MoS two <ref:2602.22463#pg1,second-harmonic generation in monolayer>.

Mira: They’ve shown that by incorporating these three-particle effects through that static long-range quadratic xc kernel, they bridge the gap between theory and experiment, getting excellent agreement with experiment on both peak positions and intensities.

Lev: It confirms that three-body effects can play a substantial role in second-order optical processes, which is important for our understanding of how these materials behave under intense light fields.

Kai: And they suggest that we shouldn't expect those three-particle correlations to vanish; theoretical frameworks need to explicitly include them for predicting other quadratic optical responses too.

Mira: They are establishing a framework for improving predictions of other quadratic optical responses in emerging materials by showing that these three-body effects are indispensable.

Lev: For us, it means that when we move toward modeling these complex nonlinear optical systems, we need to be prepared to include these higher-order many-body effects from the start.

Kai: So that’s it for this paper on "Many-Body effects beyond excitons in second-harmonic generation of monolayer MoS two <ref:2602.22463#pg1,Many-Body effects beyond excitons in second-harmonic generation of monolayer>." It closes a quantitative gap between theory and experiment for SHG in low-dimensional semiconductors.

Department of Physics, University of the Basque Country UPV/EHU · Centro de Física de Materiales, Universidad del País Vasco UPV/EHU · Ikerbasque Foundation · Donostia International Physics Center

cond-mat.mtrl-sci, cond-mat.mes-hall

Submitted: 2026-02-25

Updated: 2026-10-22

Comments: 14 pages, 4 figures

DOI: 10.21468/SciPostPhys.21.3.068

License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/

Importance score: 83/100

The gist: The gist: The inclusion of three-particle trionic correlations through a static long-range quadratic exchange-correlation kernel achieves significantly improved quantitative agreement with experiment

Key concepts

Second Harmonic Generation (SHG): SHG is a nonlinear optical process where light of one frequency is converted into light at exactly twice the original frequency. It measures how strongly a material responds to this specific type of intense, non-linear interaction with light.
Excitonic Effects: These effects arise from the interaction between an electron and a hole bound together in a semiconductor, forming an exciton. The paper uses these two-particle interactions to describe linear optical absorption and initial SHG responses, which are important but insufficient for full accuracy.
Three-Particle Trionic Correlations: These refer to complex many-body interactions involving three interacting particles (e.g., electron-hole-hole). The paper shows that including these three-particle correlations through a specific kernel drastically improves the theoretical prediction of SHG intensity, showing they are comparable in importance to two-particle effects.

Terminology

Many-Body Effects — 2× in this paper · explained in 2 episodes
Ab Initio — 1× in this paper · explained in 1 episode

Summary

The gist: The inclusion of three-particle trionic correlations through a static long-range quadratic exchange-correlation kernel achieves significantly improved quantitative agreement with experiment for second-harmonic generation in monolayer MoS2, highlighting the role of many-body interactions beyond the excitonic level in accurately describing second-order optical responses.

Theoretical Framework

The study employs a combination of many-body perturbation theory and time-dependent current-density functional theory within an ab initio framework in the optical limit. The theoretical description utilizes tensorial many-body description of second-harmonic generation, employing time-dependent current-density functional theory (TDCDFT) to represent the tensorial extension of TDDFT.

The response functions are formally equivalent to contracting the correlation parts of the two- and three-particle Green’s functions, establishing a correspondence between two- and three-body interactions and the linear fxc and quadratic gxc xc kernels of TDDFT. The linear contribution is encapsulated by fxc encapsulating electron-hole interactions leading to the formation of two-particle bound states, or excitons. In turn, gxc incorporates electron-electron-hole and electron-hole-hole interactions. The in-plane lattice parameter was set equal to the experimental bulk value, i.e., 3.16 Å. A supercell height of Lsc = 20 Å was used in the out-of-plane direction. The effective monolayer thickness employed here is Leff = 7.3 Å, obtained from the spatial extent of the charge density along the out-of-plane direction.

Electronic Band Structure and Linear Optical Response

The LDA calculation predicts that the system is a semiconductor with a direct energy band gap of about 1.65 eV at high symmetry point K. Including quasiparticle corrections via the one-shot G0W0 approximation predicts a direct band gap of 2.52 eV, much closer to the experimentally reported value of ∼ 2.34 eV. The BSE results exhibit very good agreement with experimental data, capturing the spin-split A and B exciton double-peak structure between 2 and 2.15 eV, followed by the C peak around 2.7 eV.

The linear LRC screening tensor ALRC is derived from the macroscopic optical dielectric tensor using Eq. (4). This dynamical extension of the quantity is essential in order to reproduce subtle excitonic features within TDDFT starting from very good evaluations of the macroscopic optical dielectric function that can be, for instance, the result of BSE.

Two- and Three-Body Effects on Second-Harmonic Generation

The calculation proceeds by gradually turning on many-body interactions to calculate the macroscopic SHG susceptibility tensor from first principles. Including only excitonic effects at the two-particle level (BLRC = 0 a.u.) yields an intensity of approximately 0.14 nm 2/V, significantly above the experimental value of 0.08 nm 2/V.

The second term in Eq. (3) accounts for many-body effects at the three-particle level, which are themselves linearly screened, as the quadratic LRC screening tensor BLRC correlates three interacting linear optical responses. The study finds that for values Bxx yLRC ∈ [−4,0] a.u., the impact on the maximum SHG response is significant, reducing it by nearly a factor of two. The theoretical results match the experimental peak intensity for BLRC = −3.5 a.u., which is identified as optimal for the quadratic LRC screening parameter. This result demonstrates that three-particle correlations can be comparable in magnitude to twoparticle contributions and are therefore indispensable for quantitative theoretical predictions of nonlinear optical responses, such as SHG.

Conclusions and Outlook

The ab initio results demonstrate that many-body interactions beyond the usual two-particle excitonic level are necessary for a quantitative description of the SHG spectrum in monolayer MoS2. By incorporating trionic effects through a static long-range quadratic xc kernel, the study bridges the gap and obtains excellent agreement with experiment in both peak positions and intensities. These findings reveal that three-body effects can play an important role in second-order optical processes.

Building on this analysis, the work suggests that three-particle correlations to quadratic responses cannot be expected to vanish, necessitating theoretical frameworks that explicitly include them. This establishes a framework for improved predictions of other quadratic optical responses in emerging materials.

The paper concludes by comparing the maximum SHG peak intensity with the experimental value as a function of BLRC, showing that three-body effects are indispensable for quantitative theoretical predictions of nonlinear optical responses, such as SHG. The study closes a quantitative gap between theory and experiment for SHG in low-dimensional semiconductors.

Computational Details

Ground-state electronic structure calculations were performed using density functional theory (DFT) in the local density approximation (LDA). Quasiparticle corrections due to electron-electron interactions were calculated within the one-shot G0W0 approximation. Linear optical absorption spectra including excitonic effects were obtained through BSE calculations.

The calculation of independent-particle SHG susceptibility tensors was performed using length-gauge expressions in the optical and long-wavelength limit. The use of a dense k-point interpolated mesh enabled by the Wannier interpolation technique significantly enhanced numerical convergence and accuracy.

The TDCDFT equations were formulated in the optical and long-wavelength limit, yielding current–density response equations that are independent of the wave vector. These equations are free of any kind of divergence as can be seen.

The ultimate goal is to calculate the macroscopic optical susceptibility tensor from first principles by performing a macroscopic average of the interacting optical polarizability tensor. This task is done in atomic units using expressions for first and second order <ref:2602.

Improvements for AI systems

Improved AI system can accurately predict second-harmonic generation (SHG) magnitudes in monolayer MoS2 by incorporating three-particle trionic correlations through a static long-range quadratic exchange-correlation kernel, which the paper shows leads to significantly improved quantitative agreement with experiment.
The AI system will be able to distinguish between two different theoretical approaches for calculating many-body effects, specifically contrasting the limitation of neglecting three-particle contributions in state-of-the-art TDDFT with the necessity of including them, as the paper states: "This constitutes a significant limitation: just as the two-body electron-hole interaction is essential for describing linear optical absorption, nonlinear optical responses are likewise expected to be increasingly governed by many-body interactions beyond the two-particle level."
The improved system can model and predict the frequency dependence of macroscopic SHG susceptibility tensors using a Dyson like equation that explicitly includes both two-particle excitonic effects and three-particle trionic effects, as shown in Eq. B.24, leading to excellent agreement with experiment in both peak positions and intensities.
The AI system can perform accurate material modeling by extracting the material-dependent parameters ALRC and BLRC from BSE results, allowing for the calculation of the linear LRC screening tensor and its quadratic counterpart which are material-dependent tensors that parameterize the strength of the long-range contributions (LRC) of the linear and quadratic xc kernels in the optical limit.
The system will be able to predict shifts and injection current contributions to bulk photovoltaic effects (BPVE) by adapting Eq. B.18 for processes where the first term on the r.h.s. is non-zero, thereby correctly accounting for many-body effects beyond the two-particle level in these specific quadratic optical responses.