All quantum chemical calculations in this work were carried out by Gaussian 16(A.03) package. The analysis of electronic structures was performed by Multiwfn 3.8(dev) code based on the output files and checkpoint files of Gaussian. The isosurfaces of various real space functions and the unit sphere representations of the (hyper)polarizability were rendered by means of VMD software based on the files exported from Multiwfn.