Sum frequency generation from real-time simulations in two-dimensional crystals

Sum frequency generation (SFG) and difference frequency generation (DFG) are second order nonlinear processes where two lasers with frequencies ω1 and ω2 combine to produce a response at frequency ω=ω1±ω2 . Compared with other nonlinear responses such as second-harmonic generation, SFG and DFG allow for tunability over a larger range. Moreover, the optical response can be enhanced by selecting the two laser frequencies in order to match specific electron-hole transitions. Here, we propose a first-principles framework based on the real-time solution of an effective Schr\"odinger equation to calculate the SFG and DFG in various systems, such as bulk materials, 2D materials, and molecules. Within this framework, one can select from various levels of theory for the effective one-particle Hamiltonian to account for local-field effects and electron-hole interactions. To assess the approach, we calculate the SFG and DFG of two-dimensional crystals, h-BN and MoS2 monolayers, both within the independent-particle picture and including many-body effects. Additionally, we demonstrate that our approach can also extract higher-order response functions, such as field-induced second-harmonic generation. We provide an example using bilayer h-BN.