An accelerated simulation method and related apparatus for graded periodic polarizer holographic gratings

Abstract

This invention discloses an accelerated simulation method and related apparatus for a graded periodic polarizer holographic grating. The invention performs D(k) calculations on the predicted direction of the incident light after passing through the grating layer. in After calculation, the result is multiplied by the corresponding incident light amplitude to calculate the diffraction efficiency. Compared with the traditional simulation method that uses complete and rigorous coupled wave calculation, this method reduces the number of integrations, matrix multiplications and matrix inversions, which can greatly accelerate the simulation speed.

Description

Technical Field

[0001] This invention relates to an accelerated simulation method and related apparatus for a graded periodic polarizer holographic grating, belonging to the field of holographic optical material simulation. Background Technology

[0002] Augmented reality (AR), as a next-generation smart display technology, offers immersive and practical experiences, enhances interactivity and visual communication, and improves efficiency, quality, and user engagement, thus attracting widespread attention. Among AR display technologies, polarization volume gratings (PVGs) possess unique advantages in the fields of optics and optoelectronics, making them particularly suitable for augmented reality display systems.

[0003] To meet the ever-increasing demand for field of view, graded-period polarimetric holographic gratings have been proposed in recent years. Compared with traditional fixed-period PVGs, graded-period PVGs enhance their diffraction characteristics and wavelength bandwidth. Simulating light passing through graded-period polarimetric holographic gratings has become a research direction of great interest, which involves using computer programming to simulate the behavior of light, i.e., simulating graded-period polarimetric holographic gratings.

[0004] Traditional simulation methods using rigorous coupled-wave analysis, such as those in Zemax software, require multiple integrations, matrix multiplications, and matrix inversions. Furthermore, since gradient periodic PVGs typically need to be simulated in multiple layers, the simulation efficiency is significantly lower compared to fixed-period PVGs, which poses a challenge for parameter optimization in optical waveguide design. Summary of the Invention

[0005] This invention provides an accelerated simulation method and related apparatus for a graded periodic polarizer holographic grating, which solves the problem of slow simulation efficiency in traditional methods.

[0006] According to one aspect of this disclosure, an accelerated simulation method for a graded periodic polarizer holographic grating is provided, comprising:

[0007] The gradient periodic polarization holographic grating to be simulated is segmented, and a mixing matrix of the gradient periodic polarization holographic grating is constructed based on each segmented grating layer; wherein, the mixing matrix includes amplitude phase factors in the electric field forward propagation mode, electric field backward propagation mode, magnetic field forward propagation mode, and magnetic field backward propagation mode.

[0008] The mixing matrix of the gradient periodic polarizer holographic grating and the preset wave vector direction k are used to... in and amplitude E inThe incident light is matched with the electromagnetic field boundary conditions. Based on the matched incident light information and the mixing matrix of the gradient periodic polarizer holographic grating, the coupling relationship between the electromagnetic fields of the two surfaces of the grating layer is calculated.

[0009] Based on the coupling relationship between the electromagnetic fields of the two surfaces of the grating layer and the preset grating parameters, the direction of the incident light passing through the grating layer is deduced.

[0010] The direction of the predicted incident light after passing through the grating layer is calculated using D(k). in Accelerate computation;

[0011] D(k) in Accelerated calculation results and corresponding amplitude E in Multiply the diffraction efficiency of the grating layer for the incident light to be calculated.

[0012] In some embodiments of this disclosure, a mixing matrix of a gradient periodic polarizer holographic grating is constructed based on each segmented grating layer, including:

[0013] Each grating layer after segmentation is regarded as a single-layer polarizing holographic grating with a constant longitudinal period, and a hybrid matrix of the single-layer polarizing holographic grating is constructed.

[0014] Based on the mixing matrix of a single-layer polarizer holographic grating, the mixing matrix of a gradient periodic polarizer holographic grating is obtained.

[0015] In some embodiments of this disclosure, the coupling relationship between the electromagnetic fields of the two surfaces of the grating layer is calculated using the following formula:

[0016] ;

[0017] ;

[0018] In the formula, E > and E < Electric fields H in forward and reverse propagation modes, respectively. > and H < Magnetic fields in forward and reverse propagation modes, respectively. H For the mixing matrix of a graded periodic polarizer holographic grating, p > (d) and p < (-d) < The amplitude phase factors, e, are the amplitude and phase factors in the forward and reverse propagation modes, respectively. > e < h > and h <Let W be the sub-blocks of the electric field forward propagation mode, the electric field backward propagation mode, the magnetic field forward propagation mode, and the magnetic field backward propagation mode, and let W be the eigenvector matrix derived from a single-layer polarization holographic grating using a rigorous coupled-wave algorithm.

[0019] In some embodiments of this disclosure, the diffraction efficiency of the grating layer for incident light is calculated using the following formula:

[0020] ;

[0021] In the formula, E R E represents the amplitude of the reflected light. T Let V be the amplitude of the transmitted light, V be the linear variation of the magnetic field amplitude with respect to the electric field amplitude in the waveguide medium, and I be the identity matrix. H 11 , H 12 , H 21 and H 22 For H The submatrix obtained by midpoint partitioning H It is the mixing matrix of a gradient periodic polarizer holographic grating.

[0022] According to another aspect of this disclosure, an accelerated simulation apparatus for a graded periodic polarizer holographic grating is provided, comprising:

[0023] The global mixing matrix module divides the gradient periodic polarization holographic grating to be simulated into segments and constructs a mixing matrix for each segmented grating layer. The mixing matrix includes amplitude and phase factors in the forward propagation mode of the electric field, the backward propagation mode of the electric field, the forward propagation mode of the magnetic field, and the backward propagation mode of the magnetic field.

[0024] The coupling module combines the mixing matrix of the graded periodic polarizer holographic grating with a preset wave vector direction k. in and amplitude E in The incident light is matched with the electromagnetic field boundary conditions. Based on the matched incident light information and the mixing matrix of the gradient periodic polarizer holographic grating, the coupling relationship between the electromagnetic fields of the two surfaces of the grating layer is calculated.

[0025] The direction module deduces the direction of the incident light through the grating layer based on the coupling relationship between the electromagnetic fields of the two surfaces of the grating layer and the preset grating parameters.

[0026] D(k in The acceleration module performs D(k) calculations on the predicted direction of the incident light after passing through the grating layer. in Accelerate computation;

[0027] The diffraction efficiency module will convert D(k) into D(k) in Accelerated calculation results and corresponding amplitude E in Multiply the diffraction efficiency of the grating layer for the incident light to be calculated.

[0028] In some embodiments of this disclosure, the global mixing matrix module constructs a mixing matrix for a gradient periodic polarizer holographic grating based on each segmented grating layer, including:

[0029] Each grating layer after segmentation is regarded as a single-layer polarizing holographic grating with a constant longitudinal period, and a hybrid matrix of the single-layer polarizing holographic grating is constructed.

[0030] Based on the mixing matrix of a single-layer polarizer holographic grating, the mixing matrix of a gradient periodic polarizer holographic grating is obtained.

[0031] In some embodiments of this disclosure, the coupling relationship module calculates the coupling relationship between the electromagnetic fields of the two surfaces of the grating layer using the following formula:

[0032] ;

[0033] ;

[0034] In the formula, E > and E < Electric fields H in forward and reverse propagation modes, respectively. > and H < Magnetic fields in forward and reverse propagation modes, respectively. H For the mixing matrix of a graded periodic polarizer holographic grating, p > (d) and p < (-d) < The amplitude phase factors, e, are the amplitude and phase factors in the forward and reverse propagation modes, respectively. > e < h > and h < Let W be the sub-blocks of the electric field forward propagation mode, the electric field backward propagation mode, the magnetic field forward propagation mode, and the magnetic field backward propagation mode, and let W be the eigenvector matrix derived from a single-layer polarization holographic grating using a rigorous coupled-wave algorithm.

[0035] In some embodiments of this disclosure, the diffraction efficiency module calculates the diffraction efficiency of the grating layer for incident light using the following formula:

[0036] ;

[0037] In the formula, E R E represents the amplitude of the reflected light. TLet V be the amplitude of the transmitted light, V be the linear variation of the magnetic field amplitude with respect to the electric field amplitude in the waveguide medium, and I be the identity matrix. H 11 , H 12 , H 21 and H 22 For H The submatrix obtained by midpoint partitioning H It is the mixing matrix of a gradient periodic polarizer holographic grating.

[0038] According to another aspect of this disclosure, a computer-readable storage medium is provided that stores one or more programs, the one or more programs including instructions that, when executed by a computing device, cause the computing device to perform an accelerated simulation method for a gradient periodic polarizer holographic grating.

[0039] According to another aspect of this disclosure, a computer device is provided, including one or more processors and one or more memories, wherein one or more programs are stored in the one or more memories and configured to be executed by the one or more processors, and the one or more programs include instructions for performing an accelerated simulation method for a gradient periodic polarizer holographic grating.

[0040] The beneficial effects achieved by this invention are: This invention performs D(k) on the predicted direction of the incident light after passing through the grating layer. in After calculation, the result is multiplied by the corresponding incident light amplitude to calculate the diffraction efficiency. Compared with the traditional simulation method that uses complete and rigorous coupled wave calculation, this method reduces the number of integration, matrix multiplication and matrix inversion operations, which can significantly accelerate the simulation speed. Attached Figure Description

[0041] Figure 1 A flowchart of an accelerated simulation method for a graded periodic polarizer holographic grating;

[0042] Figure 2 A schematic diagram of a holographic grating for a commonly shaped gradient periodic polarizer;

[0043] Figure 3 This is an internal cross-sectional view of a holographic grating with a gradient periodic polarizer.

[0044] Figure 4 This is a graph showing the relationship between speedup ratio and the number of grating layers.

[0045] Figure 5 This is a block diagram of an accelerated simulation device for a graded periodic polarizer holographic grating. Detailed Implementation

[0046] The technical solutions of the embodiments of this disclosure will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of this disclosure, and not all of them. The following description of at least one exemplary embodiment is merely illustrative and is in no way intended to limit this disclosure or its application or use. All other embodiments obtained by those skilled in the art based on the embodiments of this disclosure without creative effort are within the scope of protection of this disclosure.

[0047] Unless otherwise stated, the relative arrangement, numerical expressions, and values ​​of the components and steps set forth in these embodiments do not limit the scope of this disclosure.

[0048] At the same time, it should be understood that, for ease of description, the dimensions of the various parts shown in the accompanying drawings are not drawn according to actual scale.

[0049] Techniques, methods, and equipment known to those skilled in the art may not be discussed in detail, but where appropriate, such techniques, methods, and equipment should be considered part of the specification.

[0050] In all examples shown and discussed herein, any specific values ​​should be interpreted as merely exemplary and not as limitations. Therefore, other examples of exemplary embodiments may have different values.

[0051] It should be noted that similar symbols and letters in the following figures represent similar items; therefore, once an item is defined in one figure, it does not need to be discussed further in subsequent figures.

[0052] To address the slow simulation efficiency of traditional methods, this disclosure proposes an accelerated simulation method for gradient periodic polarization holographic gratings.

[0053] Figure 1 This is a schematic diagram of an embodiment of the accelerated simulation method for a graded periodic polarizer holographic grating disclosed herein. Figure 1 The implementation can be executed by a simulation terminal (such as a computer).

[0054] like Figure 1 As shown, in step 1 of the embodiment, the gradient periodic polarization holographic grating to be simulated is segmented, and a hybrid matrix of the gradient periodic polarization holographic grating is constructed according to each segmented grating layer; wherein, the hybrid matrix includes amplitude phase factors in the electric field forward propagation mode, the electric field backward propagation mode, the magnetic field forward propagation mode, and the magnetic field backward propagation mode.

[0055] It should be noted that before the simulation, a series of definitions will be made for the holographic grating of the gradually polarized body to be simulated, such as defining its parameters and shape.

[0056] The defined parameters may include vertex, length, width, thickness above the grating, thickness below the grating, refractive index of the waveguide and the surrounding air, determination of the state of the light on the waveguide, and changes in position and amplitude of the light after exiting the exit surface. Among these, position changes may include changes in the position of the light reaching the upper and lower walls of the waveguide after total internal reflection, and amplitude changes may include calculations of the amplitude changes of the TE and TM components in the light.

[0057] See Figure 2 and 3 The defined shape mainly determines whether the gradient periodic polarizer holographic grating to be simulated is a trapezoidal grating, a circular grating, a triangular grating, etc. Specifically, the shape can be defined by selecting the grating shape identifier from the established grating library, and defining the position of the four vertices of the trapezoidal grating, defining the center position and radius of the circular grating, initializing the parameters of the circular grating, and the result of the light rays exiting the gradient periodic polarizer holographic grating at the exit surface, etc.

[0058] In some embodiments, a mixing matrix of a gradient periodic polarizer holographic grating is constructed based on each segmented grating layer, which may specifically include:

[0059] 1) Treat each grating layer after segmentation as a single-layer polarizing holographic grating with a constant longitudinal period, and construct a hybrid matrix of the single-layer polarizing holographic grating.

[0060] It should be noted that a gradient periodic polarizing holographic grating can be divided into a single-layer polarizing holographic grating with a fixed longitudinal period, and the hybrid matrix of the single-layer polarizing holographic grating can be constructed by the formula of the coupling relationship module.

[0061] The hybrid matrix includes amplitude phase factors and the Fourier series of the dielectric tensor under the forward propagation mode of the electric field, the backward propagation mode of the electric field, the forward propagation mode of the magnetic field, and the backward propagation mode of the magnetic field. To describe the relative transformation between the forward and backward propagation modes, the W matrix can be divided into four sub-blocks of a two-row, two-column matrix: sub-blocks under the forward propagation mode of the electric field, the backward propagation mode of the electric field, the forward propagation mode of the magnetic field, and the backward propagation mode of the magnetic field. These can be represented by the letter e. > e < h > h < W is the eigenvector matrix derived from a single-layer polarization holographic grating using a rigorous coupled-wave algorithm, used to represent the amplitude corresponding to various propagation modes in the electromagnetic field.

[0062] 2) Based on the mixing matrix of a single-layer polarizer holographic grating, obtain the mixing matrix of a graded-periodic polarizer holographic grating. The mixing matrix of the graded-periodic polarizer holographic grating is similar to that of the single-layer grating, also including the amplitude phase factor and the dielectric tensor Fourier series matrix under the forward propagation mode of the electric field, the backward propagation mode of the electric field, the forward propagation mode of the magnetic field, and the backward propagation mode of the magnetic field.

[0063] It should be noted that the mixing matrix of the gradient periodic polarizer holographic grating is first initialized, and then the elements in the initial mixing matrix are replaced by the elements in the mixing matrix of each single-layer polarizer holographic grating, and finally the final mixing matrix of the gradient periodic polarizer holographic grating is obtained.

[0064] return Figure 1 In step 2 of the embodiment, the mixing matrix of the gradient periodic polarizer holographic grating and the preset wave vector direction k are combined. in and amplitude E in The incident light is matched with the electromagnetic field boundary conditions, that is, the mixing matrix of the graded periodic polarizing holographic grating and the preset boundary conditions with wave vector direction k. in and amplitude E in The incident light is substituted into the electromagnetic field boundary conditions for calculation. Based on the matched incident light information and the mixing matrix of the graded periodic polarizer holographic grating, the coupling relationship between the electromagnetic fields of the two surfaces of the grating layer is calculated, that is, the mixing matrix is ​​calculated according to the formula of the coupling relationship module. H .

[0065] It should be noted that the incident light information is predefined, specifically the light structure, amplitude, direction (projection in each direction), position, energy, wavelength, and other independent variables.

[0066] The coupling relationship between the electromagnetic fields on the two surfaces of the grating layer can be expressed as:

[0067] ;

[0068] ;

[0069] In the formula, E > and E < Electric fields H in forward and reverse propagation modes, respectively. > and H < Magnetic fields in forward and reverse propagation modes, respectively. H For the mixing matrix of a graded periodic polarizer holographic grating, p > (d) and p < (-d) < These are the amplitude phase factors in the forward and reverse propagation modes, respectively.

[0070] return Figure 1 In step 3 of the embodiment, the direction of the incident light after passing through the grating layer is predicted based on the coupling relationship between the electromagnetic fields of the two surfaces of the grating layer.

[0071] It should be noted that in a waveguide system, different rays emitted from the same pixel in the image source have the same incident direction on the input coupling grating after passing through the collimating lens. Therefore, the incident direction of such rays is finite and predictable.

[0072] return Figure 1 In step 4 of the embodiment, the predicted direction of the incident light after passing through the grating layer is calculated using D(k). in The calculation involves multiplication and inverse matrix operations in the diffraction efficiency module. The computationally intensive matrix D depends only on k. in .

[0073] return Figure 1 In step 5 of the embodiment, D(k) in ) Calculation results and corresponding amplitude E in Multiply the diffraction efficiency of the grating layer for the incident light to be calculated.

[0074] It should be noted that, specifically, D(k) in ) Calculation results and corresponding amplitude E in Multiplication can be expressed by the formula:

[0075] ;

[0076] In the formula, E R E represents the amplitude of the reflected light. T Let V be the amplitude of the transmitted light, V be the linear variation of the magnetic field amplitude with respect to the electric field amplitude in the waveguide medium, and I be the identity matrix. H 11 , H 12 , H 21 and H 22 For H The submatrix obtained by midpoint partitioning is specifically the submatrix obtained by... H Divide it into a matrix with two rows and two columns.

[0077] This invention first performs D(k) on the predicted direction of the incident light after passing through the grating layer. in After calculation, the result is multiplied by the corresponding incident light amplitude to calculate the diffraction efficiency. Compared with the traditional simulation method that uses complete and rigorous coupled wave calculation, this method reduces the number of integrations, matrix multiplications and matrix inversions. Therefore, the acceleration method proposed in this invention can greatly speed up the simulation.

[0078] To verify the above method, it was compared with the traditional rigorous coupled-wave method. The verification experiment involved simulation calculations on 10,000 beams of light with different amplitudes in the same direction on graded-periodic polarizer holographic gratings with varying numbers of grating layers. Simulations were performed for grating layers of 3, 5, 10, 15, 20, 25, and 30 layers. When the number of grating layers was 3, the above method took 12.15 seconds, while the traditional method took 649.49 seconds. The speedup ratio obtained by dividing the time taken by the traditional method by the time taken by the above method, and its relationship with the number of grating layers, is shown in [reference needed]. Figure 4 As can be seen from the above comparison, the simulation efficiency of the above method far exceeds that of the traditional method.

[0079] Figure 5 This is a schematic diagram of one embodiment of the accelerated simulation device for the graded periodic polarizer holographic grating disclosed herein. Figure 5 The embodiment is a virtual device that can be loaded and executed by a simulation terminal (such as a computer), including a global hybrid matrix module, a coupling relationship module, a direction module, and D(k) in ) module and diffraction efficiency module.

[0080] The global mixing matrix module of the embodiment is configured to divide the gradient periodic polarization holographic grating to be simulated into segments, and construct a mixing matrix of the gradient periodic polarization holographic grating according to each segmented grating layer; wherein, the mixing matrix includes amplitude phase factors in the electric field forward propagation mode, the electric field backward propagation mode, the magnetic field forward propagation mode, and the magnetic field backward propagation mode.

[0081] It should be noted that in the global mixing matrix module, the mixing matrix of the gradient periodic polarizer holographic grating is constructed based on each segmented grating layer. This includes: treating each segmented grating layer as a single-layer polarizer holographic grating with a constant longitudinal period, and constructing the mixing matrix of the single-layer polarizer holographic grating; and obtaining the mixing matrix of the gradient periodic polarizer holographic grating based on the mixing matrix of the single-layer polarizer holographic grating.

[0082] The coupling module in this embodiment is configured to combine the mixing matrix of the gradient periodic polarizer holographic grating with a preset wave vector direction k. in and amplitude E in The incident light is matched with the electromagnetic field boundary conditions. Based on the matched incident light information and the mixing matrix of the gradient periodic polarizer holographic grating, the coupling relationship between the electromagnetic fields of the two surfaces of the grating layer is calculated.

[0083] It should be noted that in the coupling module, the formula for calculating the coupling relationship between the electromagnetic fields of the two surfaces of the grating layer is as follows:

[0084] ;

[0085] ;

[0086] In the formula, E > and E < Electric fields H in forward and reverse propagation modes, respectively. > and H < Magnetic fields in forward and reverse propagation modes, respectively. H For the mixing matrix of a graded periodic polarizer holographic grating, p > (d) and p < (-d) < The amplitude phase factors, e, are the amplitude and phase factors in the forward and reverse propagation modes, respectively. > e < h > and h < Let W be the sub-blocks of the electric field forward propagation mode, the electric field backward propagation mode, the magnetic field forward propagation mode, and the magnetic field backward propagation mode, and let W be the eigenvector matrix derived from a single-layer polarization holographic grating using a rigorous coupled-wave algorithm.

[0087] The orientation module in this embodiment is configured to predict the direction of incident light after passing through the grating layer based on the coupling relationship between the electromagnetic fields of the two surfaces of the grating layer.

[0088] D(k) in the example in The module is configured to perform D(k) calculations on the predicted direction of the incident light after passing through the grating layer. in )calculate.

[0089] The diffraction efficiency module in this embodiment is configured to reduce D(k) in Calculation results and corresponding amplitude E in Multiply the diffraction efficiency of the grating layer for the incident light to be calculated.

[0090] It should be noted that in the diffraction efficiency module, the formula for calculating the diffraction efficiency of the grating layer for incident light is:

[0091] ;

[0092] In the formula, E R E represents the amplitude of the reflected light. T Let V be the amplitude of the transmitted light, V be the linear variation of the magnetic field amplitude with respect to the electric field amplitude in the waveguide medium, and I be the identity matrix. H 11 , H 12 , H 21 and H 22 For H The submatrix obtained by midpoint partitioning HIt is the mixing matrix of a gradient periodic polarizer holographic grating.

[0093] Similar to the method described above, the device first performs D(k) on the predicted direction of the incident light after passing through the grating layer. in After calculation, the result is multiplied by the corresponding incident light amplitude to calculate the diffraction efficiency. Compared with the traditional simulation method that uses complete and rigorous coupled wave calculation, this method reduces the number of integrations, matrix multiplications and matrix inversions. Therefore, the acceleration method proposed in this invention can greatly speed up the simulation.

[0094] Based on the same technical solution, this disclosure also relates to a computer-readable storage medium that stores one or more programs, the one or more programs including instructions that, when executed by a computing device, cause the computing device to perform an accelerated simulation method for a gradient periodic polarizer holographic grating.

[0095] Based on the same technical solution, this disclosure also relates to a computer device, including one or more processors and one or more memories, wherein one or more programs are stored in one or more memories and configured to be executed by one or more processors, and the one or more programs include instructions for performing an accelerated simulation method for a gradient periodic polarizer holographic grating.

[0096] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0097] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0098] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0099] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0100] The above are merely embodiments of the present invention and are not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention are included within the scope of the claims of the present invention pending approval.

Claims

1. An accelerated simulation method for a graded periodic polarizing holographic grating, characterized in that, include: The gradient periodic polarization holographic grating to be simulated is segmented, and a mixing matrix of the gradient periodic polarization holographic grating is constructed based on each segmented grating layer; wherein, the mixing matrix includes amplitude phase factors in the electric field forward propagation mode, electric field backward propagation mode, magnetic field forward propagation mode, and magnetic field backward propagation mode. The mixing matrix of the gradient periodic polarizer holographic grating and the preset wave vector direction k are used to... in and amplitude E in The incident light is matched using electromagnetic field boundary conditions. Based on the information matrix of the matched incident light and the mixing matrix of the gradient periodic polarizer holographic grating, the coupling relationship between the electromagnetic fields of the two surfaces of the grating layer is calculated. Based on the coupling relationship between the electromagnetic fields of the two surfaces of the grating layer and the preset grating parameters, the direction of the incident light after passing through the grating layer is deduced. The direction of the predicted incident light after passing through the grating layer is calculated using D(k). in )calculate; D(k) in ) Calculation results and corresponding amplitude E in Multiply the diffraction efficiency of the grating layer for the incident light; The above-mentioned construction of the mixing matrix of the gradient periodic polarization holographic grating based on each segmented grating layer includes: treating each segmented grating layer as a single-layer polarization holographic grating with a constant longitudinal period, and constructing the mixing matrix of the single-layer polarization holographic grating; obtaining the mixing matrix of the gradient periodic polarization holographic grating based on the mixing matrix of the single-layer polarization holographic grating. The above formula is used to calculate the coupling relationship between the electromagnetic fields on the two surfaces of the grating layer: In the formula, E > and E < Electric fields H in forward and reverse propagation modes, respectively. > and H < Magnetic fields in forward and reverse propagation modes, respectively, where H is the mixing matrix of the graded periodic polarizer holographic grating, and p > (d) and p < (-d) < The amplitude phase factors, e, are the amplitude and phase factors in the forward and reverse propagation modes, respectively. > e < h > and h < The W matrix is ​​divided into sub-blocks in the forward propagation mode of electric field, the backward propagation mode of electric field, the forward propagation mode of magnetic field, and the backward propagation mode of magnetic field. W is the eigenvector matrix derived from a single-layer polarization holographic grating using a rigorous coupled-wave algorithm. The formula for calculating the diffraction efficiency of the grating layer for incident light is as follows: In the formula, E R E represents the amplitude of the reflected light. T V is the amplitude of the transmitted light, V is the linear change of the magnetic field amplitude relative to the electric field amplitude in the waveguide medium, I is the identity matrix, and H is the amplitude of the transmitted light. 11 H 12 H 21 and H 22 This is the submatrix obtained by dividing H by its midpoint.

2. An accelerated simulation device for a graded periodic polarizer holographic grating, characterized in that, include: The global mixing matrix module divides the gradient periodic polarization holographic grating to be simulated into segments and constructs a mixing matrix for each segmented grating layer. The mixing matrix includes amplitude and phase factors in the forward propagation mode of the electric field, the backward propagation mode of the electric field, the forward propagation mode of the magnetic field, and the backward propagation mode of the magnetic field. The coupling module combines the mixing matrix of the graded periodic polarizer holographic grating with a preset wave vector direction k. in and amplitude E in The incident light is matched with the electromagnetic field boundary conditions. Based on the matched incident light information matrix and the mixing matrix of the gradient periodic polarizer holographic grating, the coupling relationship between the electromagnetic fields of the two surfaces of the grating layer is calculated. The direction module deduces the direction of the incident light after passing through the grating layer based on the coupling relationship between the electromagnetic fields of the two surfaces of the grating layer and the preset grating parameters. D(k in The acceleration module performs D(k) calculations on the predicted direction of the incident light after passing through the grating layer. in Accelerate computation; The diffraction efficiency module will convert D(k) into D(k) in ) Calculation results and corresponding amplitude E in Multiply the diffraction efficiency of the grating layer for the incident light; In the aforementioned global mixing matrix module, the mixing matrix of the gradient periodic polarizer holographic grating is constructed based on each segmented grating layer, including: treating each segmented grating layer as a single-layer polarizer holographic grating with a constant longitudinal period, and constructing the mixing matrix of the single-layer polarizer holographic grating; and obtaining the mixing matrix of the gradient periodic polarizer holographic grating based on the mixing matrix of the single-layer polarizer holographic grating. In the coupling relationship module described above, the coupling relationship between the electromagnetic fields of the two surfaces of the grating layer is calculated using the following formula: In the formula, E > and E < Electric fields H in forward and reverse propagation modes, respectively. > and H < Magnetic fields in forward and reverse propagation modes, respectively, where H is the mixing matrix of the graded periodic polarizer holographic grating, and p > (d) and p < (-d) < The amplitude phase factors, e, are the amplitude and phase factors in the forward and reverse propagation modes, respectively. > e < h > and h < The W matrix is ​​divided into sub-blocks in the forward propagation mode of electric field, the backward propagation mode of electric field, the forward propagation mode of magnetic field, and the backward propagation mode of magnetic field. W is the eigenvector matrix derived from a single-layer polarization holographic grating using a rigorous coupled-wave algorithm. In the diffraction efficiency module described above, the diffraction efficiency of the grating layer for incident light is calculated using the following formula: In the formula, E R E represents the amplitude of the reflected light. T V is the amplitude of the transmitted light, V is the linear change of the magnetic field amplitude relative to the electric field amplitude in the waveguide medium, I is the identity matrix, and H is the amplitude of the transmitted light. 11 H 12 H 21 and H 22 This is the submatrix obtained by dividing H by its midpoint.

3. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores one or more programs, the one or more programs including instructions that, when executed by a computing device, cause the computing device to perform the method of claim 1.

4. A computer device, characterized in that, include: One or more processors and one or more memories, one or more programs stored in one or more memories and configured to be executed by one or more processors, the one or more programs including instructions for performing the method of claim 1.

Images