Fig. 1. (a) Schematic diagram of electrically stimulated LC photonic device. (b), (c) Grids of patterned electrode meshing in similar number of model grids. (b) Uniform grids and (c) nonuniform grids. The scale bars in (b), (c) are 10 μm.

Fig. 2. (a) Schematic diagram of calculation relation. (b) Calculation sequence of SOR. (c), (d) Schematic diagram of (c) periodic and (d) symmetric boundary conditions. (e) Boundary conditions at different p values.

Fig. 3. (a) Schematic illustration of the unit cell of the LC photonic device. (b)–(j) Slice diagrams of the calculated tilt angle of the directors. (b)–(d) FEM, (e)–(g) uniform FDM, and (h)–(j) nonuniform FDM. (b), (e), (h) Along x direction where x is sampled at 50, 55, 60, 65, 70, 75, 80, and 85 μm. (c), (f), (i) Along z direction where z is sampled at 30, 35, 40, 45, and 49.5 μm. (d), (g), (j) Zoomed-in image of tilt angle diagrams of the directors.

Fig. 4. Comparison of mesh grid distribution between nonuniform and uniform methods. (a)–(c) Grid of nonuniform method (a) in xz plane, (b) in xy plane, and (c) inside the LC layer in xz plane. (d), (e) Grid size comparison between the uniform and nonuniform methods. (d) Along x or y direction. (e) Along z direction.

Fig. 5. (a) Schematic illustration of the experimental setup to characterize the LC devices. (b), (f) Experimentally recorded light field intensity distributions at the LC layer surface of (b) array and (f) one unit. (c), (d), (e) Light field simulation results of the LC layer in the transverse plane with (c) FEM, (d) uniform FDM, and (e) nonuniform FDM. (g) Light field intensity distribution of the outermost ring of the LC layer surface in polar coordinate. (h) Simulated light field intensity distribution along the longitudinal plane. (i) Experimentally recorded light field intensity distribution along the longitudinal plane. The scale bar in (b) is 50 μm; the scale bars in (c)–(f) are 10 μm.

Table 1. Calculation Information of Three Methods