The invention belongs to the field of optical detection, and in particular relates to a design method of a spherical mirror multi-pass pool with dense light spot atterns.
A long-path multipass cell based on two spherical mirrors with a reflectivity as high as 98% has been widely used in optical absorption spectroscopy and gas-phase optical delay lines. Early spherical mirror-based multipass cells developed by White and Heriott Heriott are still used today in laser-based spectroscopic trace gas sensors due to their simplicity, reliability, robustness, and operability. In order to improve the utilization efficiency of the mirror surface and obtain a longer total optical path length, Heriott designed a multipass cell based on an astigmatic mirror, in which the spherical mirror has different focal lengths in the x-z and y-z planes, resulting in Lissajous For example, the light spot pattern increases the number of light reflections in the multi-pass cell while minimizing the overlap of the light spots. Recently, a number of multipass cell variants based on astigmatic mirrors and with similarly high fill factors have been reported in the literature, in which at least one spherical mirror is replaced by a cylindrical mirror. Therefore, the current multi-pass cells with dense spot patterns are based on aspheric mirrors. In actual manufacturing, spherical surfaces of mirrors with high surface precision are produced by natural grinding and polishing techniques. The manufacturing process of aspheric surfaces is more complicated, and it is difficult to manufacture mirrors with sufficient surface accuracy to match the designed spot pattern. Therefore, spherical mirrors are more popular in the design of multi-pass cells, because their manufacturing process is simple, it is easier to control the surface quality, and more importantly, the use of spherical mirrors to design multi-pass cells can reduce costs. like
In traditional calculations based on the paraxial approximation theory, two assumptions are included: (i) of any ray between two mirrors the optical path length dn is a constant D, where n is the number of passes; (ii) all rays make small angles to the optical axis of the system such that three important angular approximations are valid, namely sin θ˜θ, tan θ˜θ, and cos θ˜1. However, with more curved reflective surfaces, especially marginal rays, and with greater number of passes, the paraxial approximation theory suffers from increasing deviations from actual performance due to spherical mirror aberration. For example, at an angle of 10°, the paraxial approximation of sin θ˜θ has 0.5% error, however these errors would not be accounted for in the conventional ABCD matrix. As the number of passes increases, these errors will be accumulated and amplified, distorting the actual spot pattern. The presence of paraxial approximations and the lack of analytical equations for actual ray trajectories during multipass cell design limit the ability to develop dense spot patterns. In particular, the aberration effect of the spherical mirror is not considered in the traditional ABCD matrix method, so that the traditional ABCD matrix cannot calculate and simulate the actual ray trajectory at all.
The present invention overcomes the deficiencies in the prior art, and the technical problem to be solved is: to provide design method for spherical mirror multi-pass pool with dense spot pattern, using an accurate ABCD matrix without paraxial approximation to describe the relationship between two spherical mirrors. The propagation of off-axis and edge rays between them and the reflection of these rays on the mirror surface, through numerical simulation, multi-pass cell with dense spot pattern can be obtained.
In order to solve the above-mentioned technical problems, the technical solution adopted in the present invention is: the design method of spherical mirror multi-pass pool with dense spot pattern, the multi-pass pool is composed of two identical spherical mirrors, and the design method includes the following steps:
S1. Determine the distance of the spherical mirror, the focal length of the spherical mirror, and the incident position (x0,y0) and the incident direction (x′0,y′0);
S2. Perform iterative calculation through the ABCD matrix without paraxial approximation to obtain the light parameters passing through the multi-pass cell each time, and the light parameters include the three-dimensional coordinates of the light beam on the two mirrors (xn,yn,zn), light inclination angle (x′n,y′n,z′n) and optical path length dn, n represents the number of passes;
S3. Project the light spots of all passing times of the light beam on the two mirrors onto the x-y plane, and observe the light spot pattern;
S4. Change the distance of the spherical mirror, the focal length of the spherical mirror, the initial incident point of the light beam, the incident direction or the number of iterations, and repeat steps S2 and S3 until the desired spot pattern is obtained;
The expression of the ABCD matrix without paraxial approximation is:
Among them, S and L represent operators respectively, Sφ=sin φ and Lφ=−2 arcsin (φ/R).
In the step S2, when the iterative operation is performed through the ABCD matrix without paraxial approximation, the iterative formula is:
In the formula, (xn,yn) and (xn+1,yn+1) represents the spot positions of the nth and n+1th incidents on the mirror surface of the multi-pass pool respectively, (x′n,y′n) and (x′n+1,y′n+1) represent the inclination angles of the nth and n+1th reflected beams respectively; dn+1 Indicates the optical path length of the n+1th transmission of the light beam in the multi-pass cell, and R represents the curvature radius of the spherical mirror that constitutes the multi-pass cell.
The optical path length dn+1 of the n+1th transmission calculation formula is:
In the formula,
The design method of a spherical mirror multi-pass pool with dense spot patterns also includes:
Step S5, when the required spot pattern is obtained, keep the initial position of the incident light on the M1 mirror surface (x0,y0) and tilt angle (x′0,y′0) remains unchanged, multiply the focal length f and the mirror distance D by the expected gain factor to adjust the spot density; or the initial position of the incident light on the M1 mirror surface (x0,y0), the focal length f and the mirror distance D are multiplied by the same scaling factor at the same time to adjust the size of the pattern in the spot.
In the step S2, when performing the iterative operation, the incident light beam is defined as a plurality of edge rays and a center ray, and then the light parameters of each edge ray and center ray passing through the multi-pass pool are calculated respectively, and these The light spot area obtained after the same number of reflections is used for contour fitting to obtain the spot position of the actual beam.
Compared with the prior art, the present invention has the following beneficial effects: the present invention proposes a design method of a spherical mirror multi-pass pool with a dense spot pattern, and describes the distance between two spherical mirrors through an accurate ABCD matrix that does not contain a paraxial approximation. The propagation of off-axis and edge rays and the reflection of these rays on the mirror surface, by changing the number of iterations, the distance of the spherical mirror, and the initial incident point and incident direction of the beam, a rich spot pattern can be generated on the spherical mirror through numerical simulation, The utilization efficiency of the spherical mirror is improved. Long total optical path lengths are usually produced by using multipass cells composed of aspheric mirrors. This design method allows the use of lower cost spherical mirrors to produce a total optical path length comparable to that of aspheric mirrors.
image 3 It is a schematic diagram of six kinds of singular spot patterns produced by a double-spherical mirror multi-pass pool according to an embodiment of the present invention;
Image 6 for the invention through
In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below. Obviously, the described embodiments are part of the embodiments of the present invention, rather than All the embodiments, based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without creative work all belong to the protection scope of the present invention.
The embodiment of the present invention provides a design method of a spherical mirror multi-pass pool with a dense spot pattern, and the design method of the multi-pass pool includes the following steps:
S1. Determine the distance of the spherical mirror, the focal length of the spherical mirror, and the incident point (x0,y0) and the incident direction (x′0,y′0);
S2. Perform iterative calculation through the ABCD matrix without paraxial
approximation to obtain the light parameters passing through the multi-pass cell each time, and the light parameters include the three-dimensional coordinates of the light beam on the two mirrors (xn,yn,zn), light inclination angle (x′,y′n,z′n) and optical path length dn, n represents the number of light passes;
S3. Project the light spots of all passing times of the light beam on the two mirrors onto the x-y plane, and observe the light spot pattern;
S4. Change the distance of the spherical mirror, the focal length of the spherical mirror, the initial incident point of the light beam, the incident direction or the number of iterations, and repeat steps S2 and S3 until the desired spot pattern is obtained.
like
The embodiment of the present invention uses an ABCD matrix that does not include a paraxial approximation for iterative operations. In this matrix, we modify the corresponding reflection matrix in the traditional ABCD. The modified reflection matrix as follows:
Define an operator L such that Lφ=−2 arcsin(φ/R), At this time, errors due to angular approximation can be avoided.
For the space free transfer matrix, it is modified as:
Define an operator S such that Sφ=sin φ, At this time, errors due to the approximation of the angle as well as the approximation of the actual propagation distance can be avoided.
In the absence of paraxial approximation, using the modified ABCD matrix to describe the x-z plane ray after a space-free transmission and a reflection in the multi-pass cell will have the following equations (the same applies to the y-z plane rays):
Therefore, the ABCD matrix given by (2) does not contain paraxial approximation.
Since d in the modified ABCD matrix dn+1 will vary with each pass, an iterative approach is used to compute the ray parameters for each pass through the multipass cell. In the Cartesian coordinate system, the transmitted and reflected ray parameters between the nth and (n+1)th passes through the multipass cell are determined by the coordinates of the light spot on the mirror surface of M1 or M2 (xn,yn) description, and its rays have an oblique angle (x′n,y′n), like
In the formula, ((xn,yn) and ((xn+1,yn+1) represents the spot positions of the nth and n+1th incidents on the mirror surface of the multi-pass pool respectively, (x′n,y′n) and (x′n,y′n) represent the inclination angles of the nth and n+1th reflected beams respectively; dn+1 Indicates the optical path length of the n+1th transmission of the light beam in the multi-pass cell, and R represents the curvature radius of the spherical mirror that constitutes the multi-pass cell.
If the incident ray enters the multipass cell from M1, the initial position on the surface of M1 (x0,y0) and tilt angle (x′n,y′n), then can using the spherical equation of M1 calculated z0 and z′0:
Therefore, zn+1 can be expressed by solving the equations for rays and spherical mirrors:
In addition, in formula (3), the optical path length dn+1 of the n+1th transmission calculation formula is:
In formula (6),
in,
As shown in Table 1, it is the beam parameter table for multiple iterative calculations through formula (3), where the incident coordinates of the beam are (8.56 mm, −5.35 mm), and the incident angle is (6.56°, 6.56°). The focal length of two identical spherical mirrors is 25 mm, and the distance between the mirrors is 61.16 mm.
, y
)(
)
, y
)(°)
)(°)
(
)
)
)
)
)
)
)
)
)
)
)
)
)
indicates data missing or illegible when filed
like
Further, the embodiment of the present invention obtains a numerical solution for incident light beams with different parameters passing through a double-spherical mirror multi-pass cell. A set of rich and singular dense spot patterns is calculated by equation (3) and drawn in image 3 middle. All light spots on the two mirrors are projected onto the x-y plane. The initial incidence angle of the incident light is (6.56°, 6.56°), and the focal length of the spherical mirror is 25 mm. Among all these spot patterns, the only different parameters of each pattern are the initial entrance position of the incident light on the M1 mirror surface, two The distance between the spherical mirrors and the number of passes are shown in Table 2. N in the table represents the total number of passes of the beam between the spherical mirrors. image 3 (a) gives a spot pattern similar to morning glory, where the “petals” face the “heart”, while image 3 (d) shows a sunflower-like pattern of light spots, with the “petals” facing outward. A rosette-like spot pattern is drawn on the image 3 (e). By adjusting the distance between the mirrors and the incident position of the incident light on the M1 mirror, a more compact spot pattern can be obtained. image 3 (b) and (f) show three and four nested ring patterns, respectively, while image 3 (c) shows seven non-intersecting ring patterns.
, y
) (mm)
, y
) (mm)
, y
)(°)
3a
3b
3c
3d
3e
3f
indicates data missing or illegible when filed
The double-spherical mirror multi-pass cell with a dense spot pattern achieves a long total optical path length, and such a multi-pass cell has a small volume. A longer total optical path length can improve detection sensitivity, and a small-volume multi-pass cell can not only improve the gas exchange rate, shorten the response time, but also realize a more compact multi-pass cell. The total optical path length L and volume V of each multipass cell are listed in Table 2. exist image 3, choose the distance of each pattern from the farthest spot to the origin as their mirror radius. The volume is defined as the product of the mirror area and the mirror pitch. The longest total optical path length is at image 3 The rosette-like spot pattern in (e) corresponds to a total optical path length of 15.30 m, while the minimum volume is three nested ring patterns. In fact, the ratio of total path length to volume, RLV, can better reflect the space utilization of ray trajectories in a multipass cell. According to column 7 in Table 2, the most efficient use of space is image 3 Of the three nested ring patterns in (b), the second is image 3 Seven non-intersecting ring patterns in (c).
The inapplicability of the paraxial approximation to these unusual spot patterns is mainly due to the large initial oblique angle of the incident rays, which plays an important role in the evolution of the spot pattern caused by spherical aberration. When the incident ray parameters are kept constant, it is possible to manipulate image 3 Each spot pattern in to increase the spot density, the coordinates of the spot pattern are scaled up or down. In order to increase the spot density, the initial position of the incident light on the M1 mirror surface (x0,y0) and tilt angle (x′0,y′0) remains constant, and the focal length f and the mirror distance D are multiplied by the expected gain factor respectively, then the number of passes n will increase or decrease with the gain factor. For example, for the following parameters: f=50 mm, D=122.79 mm, such as
When the inclination angle of the incident ray (x′0,y′0) remains unchanged, and the initial incident position (x0,y0), the focal length f and the mirror spacing D are simultaneously multiplied by the same scaling factor. the spot pattern can be scaled up or down, as shown in equation (3). For example, if image 3(d) in (x0,y0), f and D are increased by 2 times at the same time, which is (6.28 mm, −11.26 mm), f=50 mm and D=123.52 mm, using the same number of passes n can obtain a similar pattern without any distortion and a larger spot distribution area, like
Due to the nonlinear calculation in equation (3), in order to verify whether the double-spherical mirror multi-pass cell can meet the free-flowing conditions, a numerical simulation method was used to verify. In other words, after a period of N passes, the incident light can continue to be transmitted and reflected again, and the same light spot pattern can be generated again. In this simulation, using a fixed tilt angle (6.56°, 6.56°), the mirror spacing D and the initial position (x0,y0), so as to produce the final (xN,yN) and (x′N,y′N) and the initial (x0,y0) and (x′0,y′0) overlapping. In fact, image 3 The spot pattern in is a manually adjusted and screened spot pattern, and its corresponding multipass cell satisfies the condition of freedom. The last two columns in Table 2 list the spot position on the M1 mirror and for image 3 The inclination angle of the light after the Nth pass through the multi-pass cell of the six spot patterns. Based on these values, the free occurrence conditions for all six patterns can be met. In order to understand free behavior, at Attention should be paid to the periodic function sin φ in the expression of xn+1 and yn+1. if x′n+1 and y′n+1 is linear, free behavior can be achieved for any initial beam parameters. However, for nonlinear functions arcsin φ appears in x′n+1 and y′n+1, The interaction and confinement between arcsin φ and sin φ generate a periodic spot pattern and satisfy the free-flowing condition by some specific initial beam parameters.
Further, in the embodiment of the present invention, the incident light beam I can be defined as being composed of 12 marginal rays and 1 central ray, wherein the initial positions of the 12 marginal rays are placed on a circular contour with a radius of ω to present a uniform distribution, And take the position of the central ray as the center of the circle. In order to simplify the calculation, the incident beam is defined here as an ideal parallel beam. The specific laser used in practice can also match the output beam of the laser by setting the angles of the beam's contour edge rays and the x, y axes and the diameter of the spot at the incident position. The divergence angle of the beam and the spot corresponding to the incident position, for the beam whose spot of the output beam is ellipse, the beam profile model matching the actual beam can also be obtained by changing the radius of the beam on the x, y plane. The actual spot shape is obtained by tracing each individual ray as it bounces back and forth between a pair of spherical mirrors. Due to the existence of the spherical aberration of the coaxial spherical mirror, we can numerically simulate each ray that can represent the actual beam profile by using the exact ABCD matrix without paraxial approximation, and at the same time combine these rays after the same number of reflections to obtain The spot area is profile fitted to obtain the spot position of the actual beam, thus generating a rich spot pattern on the surface of the spherical mirror, and selecting the appropriate incident angle and position, as well as the appropriate beam diameter, so that the spots do not overlap.
like
The present invention uses the ABCD matrix without paraxial approximation to iteratively calculate the light beam entering the multi-pass cell, which can more accurately study the double-spherical mirror multi-pass cell in the case of spherical aberration. It turns out that every transmission and reflection accumulates and amplifies the spherical aberration, producing a massively complex spot pattern that is very different from the standard Heriott ring spot pattern. The excellent ratio of total optical path length to volume allows for more compact multipass cells. And the cost of using a pair of ordinary spherical mirrors is lower. Compact or portable multipass cells have many uses in protection as well as atmospheric detection and medical diagnostics. Also, the fabrication and testing of double spherical mirror multipass cells and the effect on beam interference in this type of multipass cells are very interesting studies.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solutions of the present invention, rather than limiting them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: It is still possible to modify the technical solutions described in the foregoing embodiments, or perform equivalent replacements for some or all of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the technical solutions of the various embodiments of the present invention scope.
Number | Date | Country | Kind |
---|---|---|---|
202210216123.3 | Mar 2022 | CN | national |
202221981020.8 | Jul 2022 | CN | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/CN2022/127415 | 10/25/2022 | WO |