Layer-Resolved Nonlinear Optics in Finite-Thickness Two-Dimensional Systems

Nonlinear optical (NLO) responses in two-dimensional quantum-confined systems are typically described within bulk-based frameworks as macroscopic spatial averages. In finite-thickness van der Waals multilayers directly relevant to nanoscale devices, this picture substantially breaks down. Here, we establish a general symmetry-based framework for classifying second-order NLO responses in multilayers. We reveal a layer-resolved organization into skin, weak-skin, and hidden effects governed by local symmetry and stacking order. First-principles calculations for both nonmagnetic and spin-polarized systems confirm our predictions, demonstrating that stacking alone suffices to dramatically reshape both the spatial pattern and magnitude of the NLO response, a phenomenon not explainable within standard bulk theory. Our results establish stacking geometry as an effective knob for engineering surface-selective NLO responses in layered materials.