Insulating moir\'e homobilayers lack a threefold symmetric second harmonic generation

Atoms within moir\'e bilayers relax in-plane to minimize elastic energy [e.g., Cazeaux et al., J. Elast. 154, 443 (2023)]; such relaxation brings their space group symmetries down to P1. Here, the ab initio second harmonic generation (SHG) of twisted and atomistically optimized hBN bilayers was determined at four twist angles ($\theta=38.21^{\circ}$, $60.00^{\circ}$, $73.17^{\circ}$, and $98.21^{\circ}$) and for three displacements $\boldsymbol{\tau}$ measured away from the ground state $AA^{\prime}$ configuration. All moir\'e bilayers have a P1 space symmetry after structural optimization. This situation is quite different to monolayers with hexagonal lattices, which retain a three-fold symmetry. We point out that the actual symmetries of the SHG reported for hBN bilayers on two experimental works do not coincide with the sixfold symmetric theoretical profiles they provide [either $\sin^2(3\phi)$ or $\cos^2(3\phi)$], and show that the intrinsic low structural symmetry of (atomically optimized) hBN bilayer moir\'es can in fact be read out from experimental SHG intensity profiles--which are tunable by $\theta$ and by the frequency of light $\omega$: The SHG is most definitely not sixfold-symmetric because moir\'es do not retain a three-fold symmetry. Furthermore, an extrinsic twofold symmetry of the SHG emission is realized by tilting the pump by an angle $\alpha$ away from the 2D material's normal, regardless of $\theta$ and $\omega$. The design of in-plane and ultrathin sources of SHG with low symmetry could be useful for the eventual creation of entanglement sources from 2D materials.