IP Library Granted Patent US 12,320,989
Granted Patent B2
US 12,320,989 · App. 17/666,717 · Granted Jun 3, 2025

Method of conical anisotropic rigorous coupled wave analysis for grating and computing device

Inventor: Lingshan Li (San Jose, CA)
Assignee: Goertek Inc.
G02B27/44G02B5/1814G02B5/1842G02B27/0101G02B27/0172G02B27/4205G02B27/4266G06F17/11
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 12,320,989
App. No.
17/666,717
Granted
Jun 3, 2025
Kind
B2
Abstract

A method of conical anisotropic rigorous coupled wave analysis for grating and a computing device are disclosed. The method includes: obtaining a target geometric phase δ′ g for the anisotropic-material-based grating; obtaining a slow axis azimuth angle ϕ c (x) of the anisotropic-material-based grating according to the target geometric phase δ′ g ; obtaining a permittivity tensor of the anisotropic-material-based grating, wherein the anisotropic-material-based grating has an ordinary index n o and an extraordinary index n e , the anisotropic-material-based grating has a slow axis polar angle θ c and slow axis azimuth angle ϕ c (x), and the permittivity tensor is based on n o , n e , θ c and ϕ c (x); applying the permittivity tensor into Maxwell equations; obtaining electromagnetic field for the anisotropic-material-based grating by using boundary conditions of at least two layers or sublayers of the anisotropic-material-based grating to obtain a diffraction efficiency for the anisotropic-material-based grating.

Claims (261)

1. A method of obtaining a diffraction efficiency of an anisotropic-material-based grating having at least two layers or sublayers in augmented reality (AR)/virtual reality (VR) devices, comprising:

obtaining, by a processor, a target geometric phase δ′ g for the anisotropic-material-based grating;

obtaining, by the processor, a slow axis azimuth angle ϕ c (x) of the anisotropic-material-based grating according to the target geometric phase δ′ g ;

obtaining, by the processor, a permittivity tensor of the anisotropic-material-based grating, wherein the anisotropic-material-based grating has an ordinary index n, and an extraordinary index ne, the anisotropic-material-based grating has a slow axis polar angle θ c and slow axis azimuth angle ϕ c (x), and the permittivity tensor is based on n o , n e , θ c and ϕ c (x);

applying, by the processor, the permittivity tensor into Maxwell equations;

obtaining, by the processor, electromagnetic field for the anisotropic-material-based grating by using boundary conditions of the at least two layers or sublayers of the anisotropic-material-based grating and Maxwell equations for each layer or sublayer, to obtain the diffraction efficiency for the anisotropic-material-based grating,

determining if the at least two layers or sublayers are suitable to be integrated into one stack for an output coupler grating or an input coupler grating of a waveguide in the augmented reality/virtual reality devices based on the diffraction efficiency; and

analysing the diffraction efficiency by a color dispersion and light control within the augmented reality/virtual reality devices.

2. The method according to claim 1 , wherein obtaining electromagnetic field for the anisotropic-material-based grating by using boundary conditions of the at least two layers or sublayers of the anisotropic-material-based grating and Maxwell equations for each layer or sublayer further includes:

obtaining a matrix F 1 for a first layer or sublayer of the anisotropic-material-based grating by using a first boundary condition between the first layer or sublayer and a region 1 through a first equation as below:

[

sin

ψ

δ

i

0

j

sin

ψ

n

1

cos

θ

δ

i

0

-

j

cos

ψ

n

1

δ

i

0

cos

ψ

cos

θ

δ

i

0

]

+

[

I

0

-

j

Y

I

0

0

I

0

-

j

Z

I

]

[

R

s

R

p

]

=

F

1

·

C

obtaining a matrix F L for a final layer or sublayer of the anisotropic-material-based grating by using a last boundary condition between the final layer or sublayer and a region 3 through a second equation as below:

F

L

·

C

=

[

I

0

-

j

Y

I

I

0

0

I

0

-

jZ

II

]

[

T

s

T

p

]

;

obtaining respective matrix F l of lth layer or sublayer by using respective boundary condition between adjacent layers or sublayers according to an iteration equation as below:

[

sin

ψ

δ

i

0

j

sin

ψ

n

1

cos

θ

δ

i

0

-

j

cos

ψ

n

1

δ

i

0

cos

ψcosθδ

i

0

]

+

[

I

0

-

j

Y

I

0

0

I

0

-

j

Z

I

]

[

R

s

R

p

]

=

l

=

1

L

F

l

·

(

F

l

)

-

1

[

I

0

-

j

Y

I

I

0

0

I

0

-

jZ

II

]

[

T

s

T

p

]

;

and

obtaining electromagnetic field for the anisotropic-material-based grating by using boundary conditions obtained through the first, second equations,

wherein region 1 is an input region for the anisotropic-material-based grating and region 3 is an output region for the anisotropic-material-based grating,

wherein a total layer or sublayer number of the anisotropic-material-based grating is L,

wherein matrix F′ l is obtained through the boundary condition between layer/sublayer (l−1) and layer/sublayer l,

wherein matrix F l represents a matrix function of lth layer or sublayer of the anisotropic-material-based grating, where l=1 . . . L, which can include or is derived from part or all of the following vectors: diagonal matrix K x , convolution permittivity matrix E, eigenvalue matrix Q l , and eigen matrix W l ,

wherein the diagonal matrix K x has i, i element being equal to k xi /k 0 , k 0 is the incident wavevector k 0 ={circumflex over (x)}k x0 +ŷk y0 +{circumflex over (z)}k z0 , and k xi is the ith order wavevector in x axis in wavevector k i ={circumflex over (x)}k xi +ŷk yi +{circumflex over (z)}k zi ,

wherein the convolution permittivity matrix E is formed by permittivity harmonic components, with i, p element being equals to a permittivity ∈ (i-p) , i represents ith diffraction order, and p represents pth layer/sublayer,

wherein the eigenvalue matrix Q l is formed by eigenvalues q l,i,m , wherein q l,i,m represents an mth eigenvalue of ith diffraction order in lth layer/sublayer,

wherein the eigen matrix W l is formed by eigenvector ω l,i,m , wherein w l,i,m represents an mth eigenvector of ith diffraction order in lth layer/sublayer,

wherein matrix I represents an identity matrix,

wherein matrix C represents constant vector used to reconstruct the linear combination of the eigenvectors,

wherein matrices Y I , Y II , Z I , and Z II represent diagonal matrices with diagonal element (k 1,zi /k 0 ), (k 2,zi /k 0 ), (k 1,zi /k 0 n 1 2 ), and (k 2,zi /k 0 n 1 2 ), respectively,

wherein ψ is an angle between an electric field vector and a plane of incidence in region 1, δ i0 is a geometric phase in region 1, n 1 is an index of the region 1, and θ is a polar angle of the electric field vector in region 1.

3. The method according to claim 2 , wherein the anisotropic-material-based grating is a liquid-crystal-based grating, and the sublayers of the liquid-crystal-based grating is formed through hologram in each layer of the liquid-crystal-based grating.

4. The method according to claim 2 , wherein the ordinary index n, and the extraordinary index n e vary with a wave length of an incident light.

5. The method according to claim 1 , wherein the anisotropic-material-based grating is a transmissive and/or reflective grating.

6. A method of obtaining parameters of an anisotropic-material-based grating, comprising:

obtaining diffraction efficiencies of different anisotropic-material-based gratings by using the method according to claim 1 ;

obtaining a merit function involving a diffraction efficiency for an anisotropic-material-based grating;

obtaining a desired anisotropic-material-based grating based on the merit function using the diffraction efficiencies of different anisotropic-material-based gratings; and

obtaining parameters of the desired anisotropic-material-based grating.

7. A computing device, comprising a non-transitory machine-readable memory and a processing unit, wherein the memory includes executable instructions, when the executable instructions are executed by the processing unit, the processing unit is made to implement the processes of the method according to claim 1 .

Assignments (2)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Dec 15, 2022
From: GOERTEK INC.
To: GOER OPTICAL TECHNOLOGY CO., LTD.
Reel/Frame 062111/0787 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Feb 23, 2022
From: LI, LINGSHAN
To: GOERTEK INC.
Reel/Frame 059078/0600 →
Continuity (1)
Related Publication 20230251500A1 · Aug 10, 2023
References Cited (31)
US 6950173B1 · Sutherland · 2005 [cited by examiner]
US 7454103B2 · Parriaux · 2008 [cited by examiner]
US 9465215B2 · Richards et al. · 2016 [cited by applicant]
US 9551880B2 · Amitai · 2017 [cited by applicant]
US 11048087B2 · Wheelwright et al. · 2021 [cited by applicant]
US 20210055554A1 · Chi et al. · 2021 [cited by applicant]
CN 112099228A · 2020 [cited by examiner]
EP 3180647B1 · 2019 [cited by applicant]
The Finite Difference Time Domain Method for Electromagnetics by Kunz KS, Luebbers RJ in CRC press; 1993; pp. 30-49. [cited by applicant]
“Computational electromagnetics: the finite-difference time-domain method” by Taflove A, Hagness SC, Piket-May M in Electr Eng Handb. 2005; pp. 51-61. [cited by applicant]
Modal theory of spatially periodic media by Chu R-S, Kong JA in IEEE Trans Microw Theory Tech. 1977;25(1):18-24. [cited by applicant]
Interpretation of MIMO channel characteristics in rectangular tunnels from modal theory by Molina-Garcia-Pardo JM, Lienard M, Degauque P, Dudley DG, Juan-Llacer L in IEEE, Trans. Veh. Technol., 2008:57 (3) :1974-1979. [cited by applicant]
“The dielectric lamellar diffraction grating” by Botten IC, Craig MS, McPhedran RC, Adams JL, Andrewartha J in Opt Acta Int J Opt. 1981;28(3):413-428. [cited by applicant]
“The finitely conducting lamellar diffraction grating” by Botten LC, Craig MS, McPhedran RC, Adams JL, Andrewartha JR in Opt Acta Int J Opt. 1981;28(8):1087-1102. [cited by applicant]
“Modal theory of diffraction by multilayered gratings containing dielectric and metallic components” by Jiang M, Tamir T, Zhang S in JOSA A. 2001;18(4):807-820. [cited by applicant]
“Coupled Wave Theory for Thick Hologram Gratings” by Kogelnik H . . . Bell Syst Tech J in 1969;48(9):2909-2947. doi:10.1002/j.1538-7305.1969.tb01198.x. [cited by applicant]
“Rigorous coupled-wave analysis of planar-grating diffraction” by Moharam MG, Gaylord TK in J Opt Soc Am. 1981;71(7):811. doi:10.1364/JOSA.71.000811. [cited by applicant]
Coupled-wave theory for multiply exposed thick holographic gratings by Case SK in JOSA. 1975;65(6):724-729. [cited by applicant]
“Rigorous coupled-wave analysis of grating diffraction—E-mode polarization and losses” by Moharam MG, Gaylord TK . . . JOSA. 1983;73(4):451-455. [cited by applicant]
“Two-dimensional coupled wave theory for square-lattice photonic-crystal lasers with TM-polarization” by Sakai K, Miyai E, Noda S in Opt Express. 2007;15(7):3981-3990. [cited by applicant]
“Rigorous coupled-wave analysis and applications of grating diffraction” by Glytsis EN, Gaylord TK, Brundrett DL in Diffractive Miniaturized Opt A Crit Rev. 1993;10271(1992). [cited by applicant]
“Formulation for stable and efficient implementation of the rigorous coupled-wave analysis of binary gratings” by Grann EB, Pommet D a., Moharam MG, Gaylord TK in J Opt Soc Am, A, 1995; 12 (5) :1068-1076. [cited by applicant]
“Three-dimensional (vector) rigorous coupled-wave analysis of anisotropic grating diffraction” by Glytsis EN, Gaylord TK in J Opt Soc Am A. 1990;7(8):1399. doi:10.1364/josaa.7. [cited by applicant]
“Stable implementation of the rigorous coupled-wave analysis for surface-relief gratings: enhanced transmittance matrix approach” by Moharam MG, Gaylord TK, Pommet DA, Grann EB, in J. Opt. Soc. Am, A, 1995; 12 (5) :1077. [cited by applicant]
“Reflective polarization vol. gratings for high efficiency waveguide-coupling augmented reality displays” by Lee Y-H, Yin K, Wu S-T in Opt Express. 2017;25(22):27008-27014. [cited by applicant]
“Bragg polarization gratings for wide angular bandwidth and high efficiency at steep deflection angles” by Xiang X, Kim J, Escuti MJ in Sci Rep. 2018;8(1):1-6. [cited by applicant]
“Computational Methods for Electromagnetic and Optical Systems” by Jarem JM, Banerjee PP in CRC press; 2016. [cited by applicant]
“Stable implementation of the rigorous coupled-wave analysis for surface-relief gratings enhanced transmittance matrix approach” by Moharam MG, Pommet DA, Grann EB in 1995. [cited by applicant]
“Numerical analysis of Bragg regime polarization gratings by rigorous coupled-wave analysis” by Xiang X, Escuti MJ in Pract Hologr XXXI Mater Appl. 2017;10127(Apr. 2017):1012. [cited by applicant]
“Color-selective geometric phase lens for apochromatic lens system” by Li L, Kim J, Shi S, Escuti MJ in In: Khoo IC, ed. Liquid Crystals XXIV. vol 11472. SPIE; 2020:38-51. [cited by applicant]
Color-selective geometric-phase lenses for focusing and imaging based on liquid crystal polymer films by Lingshan Li, Shuojia Shi, Jihwan Kim, and Michael. J. Escuti, in Opt., Express 30, 2478-2502 (2022). [cited by applicant]