OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 20, Iss. 10 — May. 7, 2012
  • pp: 11159–11165
« Show journal navigation

Optimization of refractive liquid crystal lenses using an efficient multigrid simulation

Harry Milton, Paul Brimicombe, Philip Morgan, Helen Gleeson, and John Clamp  »View Author Affiliations


Optics Express, Vol. 20, Issue 10, pp. 11159-11165 (2012)
http://dx.doi.org/10.1364/OE.20.011159


View Full Text Article

Acrobat PDF (1012 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

Abstract: A multigrid computational model has been developed to assess the performance of refractive liquid crystal lenses, which is up to 40 times faster than previous techniques. Using this model, the optimum geometries producing an ideal parabolic voltage distribution were deduced for refractive liquid crystal lenses with diameters from 1 to 9 mm. The ratio of insulation thickness to lens diameter was determined to be 1:2 for small diameter lenses, tending to 1:3 for larger lenses. The model is used to propose a new method of lens operation with lower operating voltages needed to induce specific optical powers. The operating voltages are calculated for the induction of optical powers between + 1.00 D and + 3.00 D in a 3 mm diameter lens, with the speed of the simulation facilitating the optimization of the refractive index profile. We demonstrate that the relationship between additional applied voltage and optical power is approximately linear for optical powers under + 3.00 D. The versatility of the computational simulation has also been demonstrated by modeling of in-plane electrode liquid crystal devices.

© 2012 OSA

1. Introduction

Liquid crystal lenses have shown potential for a number of display and optical applications, with recent developments in electronic ophthalmic devices [1

1. G. Li, D. L. Mathine, P. Valley, P. Ayräs, J. N. Haddock, M. S. Giridhar, G. Williby, J. Schwiegerling, G. R. Meredith, B. Kippelen, S. Honkanen, and N. Peyghambarian, “Switchable electro-optic diffractive lens with high efficiency for ophthalmic applications,” Proc. Natl. Acad. Sci. U.S.A. 103(16), 6100–6104 (2006). [CrossRef] [PubMed]

]. A unique property of liquid crystal lenses is that they are switchable, and application of a voltage activates the lens. This is perfect for presbyopia correction, non-mechanical optical zooms and holographic applications. There are many types of liquid crystal lens design, including solid lenses [2

2. S. Sato, “Liquid-crystal lens-cells with variable focal length,” Jpn. J. Appl. Phys. 18(9), 1679–1684 (1979). [CrossRef]

], refractive gradient index lenses [3

3. H. W. Ren, D. W. Fox, B. Wu, and S. T. Wu, “Liquid crystal lens with large focal length tunability and low operating voltage,” Opt. Express 15(18), 11328–11335 (2007). [CrossRef] [PubMed]

8

8. S. Sato, “Applications of liquid crystals to variable-focusing lenses,” Opt. Rev. 6(6), 471–485 (1999). [CrossRef]

] and diffractive lenses [1

1. G. Li, D. L. Mathine, P. Valley, P. Ayräs, J. N. Haddock, M. S. Giridhar, G. Williby, J. Schwiegerling, G. R. Meredith, B. Kippelen, S. Honkanen, and N. Peyghambarian, “Switchable electro-optic diffractive lens with high efficiency for ophthalmic applications,” Proc. Natl. Acad. Sci. U.S.A. 103(16), 6100–6104 (2006). [CrossRef] [PubMed]

,9

9. G. Q. Li, P. Valley, P. Ayras, D. L. Mathine, S. Honkanen, and N. Peyghambarian, “High-efficiency switchable flat diffractive ophthalmic lens with three-layer electrode pattern and two-layer via structures,” Appl. Phys. Lett. 90(11), 111105 (2007). [CrossRef]

].

A RLCL design by Ye et al. is shown in Fig. 1
Fig. 1 An illustration of a RLCL device designed by Ye et al. [9]. The voltage profile across the LC layer is determined by the powering of the active and control electrodes, with the control electrode providing a baseline voltage across the lens and the active electrodes providing additional voltage at the peripheries. Typically a thick LC layer is required of between 50 – 100 µm in order to produce sufficient phase retardation for RLCLs with diameters of several mm. A ratio of 2-3:1 between RLCL diameter and insulation thickness is usually required.
[1

1. G. Li, D. L. Mathine, P. Valley, P. Ayräs, J. N. Haddock, M. S. Giridhar, G. Williby, J. Schwiegerling, G. R. Meredith, B. Kippelen, S. Honkanen, and N. Peyghambarian, “Switchable electro-optic diffractive lens with high efficiency for ophthalmic applications,” Proc. Natl. Acad. Sci. U.S.A. 103(16), 6100–6104 (2006). [CrossRef] [PubMed]

,9

9. G. Q. Li, P. Valley, P. Ayras, D. L. Mathine, S. Honkanen, and N. Peyghambarian, “High-efficiency switchable flat diffractive ophthalmic lens with three-layer electrode pattern and two-layer via structures,” Appl. Phys. Lett. 90(11), 111105 (2007). [CrossRef]

], which uses two electrodes to induce a parabolic voltage distribution across a LC layer via a glass insulation layer placed between the LC and the electrodes. Calibration of the thickness of the insulating layer is vital for RLCL operation; insufficient insulation will result in a non-parabolic voltage distribution and excessive thickness results in high operating voltages. The insulation thickness required to produce a parabolic voltage distribution can be calculated by modeling the electric potential within the device, with a ratio of 2-3:1 between lens diameter and insulation thickness reported [6

6. M. Ye and S. Sato, “Optical properties of liquid crystal lens of any size,” Jpn. J. Appl. Phys. 41(Part 2, No. 5B), L571–L573 (2002). [CrossRef]

]. Further optimization of the RLCL dimensions using numerical modeling can be conducted to deduce the most effective insulation thickness within this range. An issue arises from simulation times and computational resources required for high resolution numerical simulations of the electric potential. Typical simulations of 3 mm diameter RLCLs require a simulation box of 300 by 150 data points, resulting in calculation times of 1-2 hours using a standard desktop computer, which is acceptable for testing existing lenses, but is generally considered too slow to use as a design tool.

2. Electric field and potential calculations

The process of solving the electric potential within a LC device requires a numerical solution to Maxwell’s equations. Assuming no enclosed charge in the system, the electric potential in a cylindrically symmetric system is given as
2φx2ε+εφxz+2φz2ε+εφzx=0,
where ϕ denotes the electric potential, ε is the permittivity of the material(s) and x and z are coordinates shown in Fig. 2
Fig. 2 Simplified simulation boxes for the modeling of RLCL. Mirror symmetry is applied to the x-boundary. In the case of RLCLs there are set boundaries in the z direction corresponding to the control and ground electrodes.
. In terms of computational modeling, standard numerical techniques are utilized using MATLAB (R2008A, MathsWorks, Natick, USA), including a successive over-relaxation (SOR) solution and red black ordering. The multigrid algorithm reduces simulation time by initially solving the equation numerically for a coarse grid, with further iterations on progressively finer grids until an adequately accurate solution is realized [11

11. U. Trottenberg, Multigrid, 1st ed. (Academic Press, London, 2001).

].

3. Simulation of RLCLs

A RLCL has relatively simple boundary conditions, as shown in Fig. 2. The system lacks free boundaries resulting in a closed simulation box. The model assumes an isotropic permittivity; hence the average electric permittivity of the LC is used for calculations. This is an approximation; in a LC material a field induced change in the director modifies the electric permittivity in the system. However, this simplification allows quick, approximate calculations of the potential profile across the system to give a good indication of device performance.

A discretized profile of the director across the LC layer is calculated using a simplified 1-dimensional model. The extraordinary refractive index (ne) is calculated for each electric potential simulation point across the top of the LC layer, giving an approximation of the refractive profile. In the following simulations, the dielectric and refractive properties of 4-cyano-4-pentylbiphenyl (5CB) are used for modeling the LC layer. While this 1-dimensional model neglects non-linear director field interactions, the approximations are acceptable for providing quick estimates of operating voltages and geometries rather than in-depth calculations. A typical simulation of a 3 mm diameter RLCL is shown in Fig. 3
Fig. 3 Electric potential calculations for a 3 mm, + 2.00 D RLCL similar to the design shown in Fig. 1 with (a) showing the potential throughout the entire RLCL device, with the bottom 100 µm constituting the LC layer, and (b) showing the potential of the LC layer. The insulating glass layer is ~1.4 mm thick yielding a bell shaped voltage distribution across the LC layer. The calculation used 4 multigrid iteration grids and a 400x150 point final simulation box to ensure high resolution across the LC layer. The average permittivity of the glass and the LC are taken to be 6.7 and 11.5 respectively (the latter measured for 5CB at room temperature).
. The electric potential in the device can be calculated using the multigrid model and the electric potential across the LC layer is deduced. The thickness of the glass and the active and control voltages have been set to produce a parabolic voltage profile.

The optimum operational geometries of RLCLs were calculated for various lens diameters. To select the optimum geometry of a RLCL, a series of comparisons were made between the potential profiles across the liquid crystal layer and an ideal parabola for a range of insulation thicknesses. To allow comparison of shape, the voltage was normalized with the minimum potential across the LC layer defined as zero and the maximum potential across the LC layer defined as one. A comparison of normalized potential curves with an ideal parabola is shown graphically in Fig. 4
Fig. 4 A comparison of insulation thickness voltage profiles for a 1 mm RLCL. Greater insulation thickness results in a better parabolic profile approximation across the optical area, with a 50 µm insulation thickness closest to a parabolic voltage profile. Here, the permittivity of the glass is 6.7 and the average permittivity of 5CB is 11.5.
. High quality parabolic voltage distributions, which correspond to superior optical quality, are shown to be clearly dependent upon the device geometry.

The relationship between the RLCL diameter and the insulation thickness for high quality parabolic voltage distributions is shown in Fig. 5(a)
Fig. 5 The relationship between insulation thickness against RLCL diameter (a) and operating voltages and optical power of a 3 mm diameter RLCL with a 140 µm thick insulation layer (b). (a) shows a comparison of minimum operational RLCL geometries for 5CB. The general ratio of 2-3:1 of the RLCL diameter to insulation thickness is observed, with a ratio of 2:1 observed at 1 mm and 2 mm, which increases to 3:1 as the RLCL diameter increases. (b) shows an approximately linear relationship between the additional active variable voltage applied and the optical power achieved with a slight reduction in control voltage due the effect of the increased variable voltages on the central area of the lens.
. Previously, the ratio of RLCL diameter to insulation thickness was found to be broadly in the range from 2:1 to 3:1 [6

6. M. Ye and S. Sato, “Optical properties of liquid crystal lens of any size,” Jpn. J. Appl. Phys. 41(Part 2, No. 5B), L571–L573 (2002). [CrossRef]

]. Our model allows the exact relationship of this ratio with respect to lens diameter to be specified, thus providing valuable insight into the optimization of the lens geometry required for parabolic voltage profiles. Many different simulations are needed to identify the optimum insulation thickness, so the speed of the multigrid solution is vital. We have discovered using our model that the optimum ratio of insulation thickness to lens diameter is strongly dependent on the lens geometry, taking values close to 2:1 for smaller lenses (~1-2 mm in diameter) and increasing towards 3:1 as the lens diameter increases.

Combining the electric potential simulation with a 1-dimesional director model allows us to assess and minimize the driving voltages needed to induce different positive focal powers. Our model compares the resulting refractive profile to a parabolic profile to assess the optimum driving voltages. In order to reduce the voltages, we propose a new method of operation whereby the refractive profile is controlled by both the active and control electrodes to give deviations in refractive index from an initial value of ne in the centre of the lens, rather than varying from no in the lens peripheries as used previously. Figure 5(b) shows the calculated operating voltages for a 3 mm diameter RLCL with a 50 µm 5CB LC layer giving optical powers between + 0.50 D and + 3.00 D. This new method of driving the lens reduces the operating voltages for low optical powers. Furthermore, the voltages can be tuned to give the most parabolic refractive profile, hence optimizing the optical quality. The calculation speed is again vital in optimizing the refractive profile due to the number of calculations required. In these lenses, an approximately linear relationship between the variable voltage and the additional optical power was observed, with each additional + 0.50 D corresponding to an increase of approximately 4 V. It was noted that at the peripheries of the lens the potential profile deviated from a parabola as the lens became more powerful, as shown in Fig. 6
Fig. 6 The refractive profiles of a + 0.50 D RLCL (a) and a + 3.00 D RLCL (b). The computer model indicates a greater deviation in the peripheries of the optical area (0-30 µm and 270-300 µ) for a stronger lens than in a weaker lens. Lab measurements give a value of ne = 1.698 and ordinary refractive index (no) = 1.535 (at 589 nm) for 5CB. In the case of 5CB, deviations from the idea parabola occur when a change in ne greater than 0.06 is required.
. This is expected when a large change in ne is required for lenses of high optical powers, as the relationship between ne and potential is increasingly non-linear at high voltages.

4. Further applications

Due to the efficient nature of the multigrid computer model, the system was adapted to simulate in-plane electrode systems which are vital in many LC display applications such as in-plane switching devices and proposed blue phase displays. The system contains free boundaries above and below the electrodes requiring boundary conditions approximating infinity, an incredibly time consuming constraint for standard numerical simulations. Analysis was undertaken using the computer model described here to assess the smoothness of the electric field within these devices for different LC materials. Low ratios of LC to glass electric permittivities (ζ) resulted in smoother electric fields as shown in Fig. 7
Fig. 7 A comparison of the electric field between inter-digited in-plane electrodes using the multigrid simulation, with a 5 µm thick LC layer above the electrodes. A smoother electric field was observed in simulations with lower LC to glass electric permittivity ratios, resulting in an improved performance.
, with distortions in the electric field at the interface minimized if the permittivity of the glass and LC are matched.

5. Conclusions

An efficient multigrid algorithm has been developed to model the electric potential in various LC devices with vastly reduced calculation times. This is an important tool for RLCL design; optimum geometries and operating voltages for RLCLs of various optical powers and sizes can be quickly screened using the model. The simulation gives a greater understanding to the previously calculated ratio of diameter to insulation thickness of between 2:1 and 3:1 for RLCLs, with additional insight into how this ratio changes with respect to RLCL diameter. The model predicts a ratio of 2:1 for small lenses of 1-2 mm in diameter, tending towards 3:1 as the lens diameter increases. Using the model we propose a new method of lens operation which minimizes operating voltages by varying the refractive profile from ne in the centre, rather than from no at the lens periphery. The speed of the multigrid simulation is vital in the optimization of the refractive profile, which maximizes the parabolicity of the refractive index profile. The model shows an approximately linear relationship between active electrode operating voltage and the optical power of the lens. Details of the limitations of RLCL performance are also shown using the model, with director profiles of high powered lenses exhibiting deviation from the parabolic distribution at the peripheries due to the increasingly non-linear relationship between ne and applied voltage at higher voltages. Finally we illustrate the versatility of our model with the simulation of in-plane electrodes systems.

Acknowledgments

The authors acknowledge a Collaborative Award in Science and Industry (CASE) from the Medical Research Council (MRC) and Ultravision CLPL.

References and links

1.

G. Li, D. L. Mathine, P. Valley, P. Ayräs, J. N. Haddock, M. S. Giridhar, G. Williby, J. Schwiegerling, G. R. Meredith, B. Kippelen, S. Honkanen, and N. Peyghambarian, “Switchable electro-optic diffractive lens with high efficiency for ophthalmic applications,” Proc. Natl. Acad. Sci. U.S.A. 103(16), 6100–6104 (2006). [CrossRef] [PubMed]

2.

S. Sato, “Liquid-crystal lens-cells with variable focal length,” Jpn. J. Appl. Phys. 18(9), 1679–1684 (1979). [CrossRef]

3.

H. W. Ren, D. W. Fox, B. Wu, and S. T. Wu, “Liquid crystal lens with large focal length tunability and low operating voltage,” Opt. Express 15(18), 11328–11335 (2007). [CrossRef] [PubMed]

4.

H. W. Ren and S. T. Wu, “Adaptive liquid crystal lens with large focal length tunability,” Opt. Express 14(23), 11292–11298 (2006). [CrossRef] [PubMed]

5.

A. F. Naumov, G. D. Love, M. Y. Loktev, and F. L. Vladimirov, “Control optimization of spherical modal liquid crystal lenses,” Opt. Express 4(9), 344–352 (1999). [CrossRef] [PubMed]

6.

M. Ye and S. Sato, “Optical properties of liquid crystal lens of any size,” Jpn. J. Appl. Phys. 41(Part 2, No. 5B), L571–L573 (2002). [CrossRef]

7.

M. Ye, B. Wang, and S. Sato, “Liquid-crystal lens with a focal length that is variable in a wide range,” Appl. Opt. 43(35), 6407–6412 (2004). [CrossRef] [PubMed]

8.

S. Sato, “Applications of liquid crystals to variable-focusing lenses,” Opt. Rev. 6(6), 471–485 (1999). [CrossRef]

9.

G. Q. Li, P. Valley, P. Ayras, D. L. Mathine, S. Honkanen, and N. Peyghambarian, “High-efficiency switchable flat diffractive ophthalmic lens with three-layer electrode pattern and two-layer via structures,” Appl. Phys. Lett. 90(11), 111105 (2007). [CrossRef]

10.

B. Wang, M. Ye, and S. Sato, “Liquid crystal lens with focal length variable from negative to positive values,” IEEE Photon. Technol. Lett. 18(1), 79–81 (2006). [CrossRef]

11.

U. Trottenberg, Multigrid, 1st ed. (Academic Press, London, 2001).

OCIS Codes
(220.3630) Optical design and fabrication : Lenses
(230.3720) Optical devices : Liquid-crystal devices

ToC Category:
Optical Design and Fabrication

History
Original Manuscript: March 22, 2012
Revised Manuscript: April 27, 2012
Manuscript Accepted: April 27, 2012
Published: April 30, 2012

Citation
Harry Milton, Paul Brimicombe, Philip Morgan, Helen Gleeson, and John Clamp, "Optimization of refractive liquid crystal lenses using an efficient multigrid simulation," Opt. Express 20, 11159-11165 (2012)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-20-10-11159


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. G. Li, D. L. Mathine, P. Valley, P. Ayräs, J. N. Haddock, M. S. Giridhar, G. Williby, J. Schwiegerling, G. R. Meredith, B. Kippelen, S. Honkanen, and N. Peyghambarian, “Switchable electro-optic diffractive lens with high efficiency for ophthalmic applications,” Proc. Natl. Acad. Sci. U.S.A.103(16), 6100–6104 (2006). [CrossRef] [PubMed]
  2. S. Sato, “Liquid-crystal lens-cells with variable focal length,” Jpn. J. Appl. Phys.18(9), 1679–1684 (1979). [CrossRef]
  3. H. W. Ren, D. W. Fox, B. Wu, and S. T. Wu, “Liquid crystal lens with large focal length tunability and low operating voltage,” Opt. Express15(18), 11328–11335 (2007). [CrossRef] [PubMed]
  4. H. W. Ren and S. T. Wu, “Adaptive liquid crystal lens with large focal length tunability,” Opt. Express14(23), 11292–11298 (2006). [CrossRef] [PubMed]
  5. A. F. Naumov, G. D. Love, M. Y. Loktev, and F. L. Vladimirov, “Control optimization of spherical modal liquid crystal lenses,” Opt. Express4(9), 344–352 (1999). [CrossRef] [PubMed]
  6. M. Ye and S. Sato, “Optical properties of liquid crystal lens of any size,” Jpn. J. Appl. Phys.41(Part 2, No. 5B), L571–L573 (2002). [CrossRef]
  7. M. Ye, B. Wang, and S. Sato, “Liquid-crystal lens with a focal length that is variable in a wide range,” Appl. Opt.43(35), 6407–6412 (2004). [CrossRef] [PubMed]
  8. S. Sato, “Applications of liquid crystals to variable-focusing lenses,” Opt. Rev.6(6), 471–485 (1999). [CrossRef]
  9. G. Q. Li, P. Valley, P. Ayras, D. L. Mathine, S. Honkanen, and N. Peyghambarian, “High-efficiency switchable flat diffractive ophthalmic lens with three-layer electrode pattern and two-layer via structures,” Appl. Phys. Lett.90(11), 111105 (2007). [CrossRef]
  10. B. Wang, M. Ye, and S. Sato, “Liquid crystal lens with focal length variable from negative to positive values,” IEEE Photon. Technol. Lett.18(1), 79–81 (2006). [CrossRef]
  11. U. Trottenberg, Multigrid, 1st ed. (Academic Press, London, 2001).

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited