“…With few isolated exceptions [13][14][15], a large number of calculations of nonlinear optical properties, and specifically of SHG in 2D crystals [16][17][18], employ the independent-particle approximation (IPA), which is inadequate for low-dimensional systems, where it is expected that the strongly bound excitons significantly modify the SHG.…”