The anisotropy of the current-induced density (ACID) analysis was realized by the ACID code based on the output file of Gaussian program, and the maps were generated by POV-Ray render . The gauge-including magnetically induced currents (GIMIC) analysis was finished via the GIMIC code based on the formatted check point file of Gaussian program, and the maps and animations were rendered by ParaView visualization program . Other electronic structure and wavefunction analyses were performed with the Multiwfn 3.8(dev) code based on the wavefunctions produced by DFT calculations. The isosurface maps of various orbitals and real space functions were rendered by means of Visual Molecular Dynamics (VMD) software based on the cube files exported from Multiwfn. The colored contour maps of various real space functions were plotted directly via the Multiwfn code.