Method and device for beam analysis
11054305 ยท 2021-07-06
Assignee
Inventors
- Matthias Manger (Aalen-Unterkochen, DE)
- Christoph Husemann (Jena, DE)
- Matus Kalisky (Aalen, DE)
- Lars Stoppe (Jena, DE)
Cpc classification
G01J1/0411
PHYSICS
G02B27/106
PHYSICS
B23K26/04
PERFORMING OPERATIONS; TRANSPORTING
G01J1/4257
PHYSICS
G01J1/0407
PHYSICS
International classification
B23K26/04
PERFORMING OPERATIONS; TRANSPORTING
B23K26/70
PERFORMING OPERATIONS; TRANSPORTING
B23K26/03
PERFORMING OPERATIONS; TRANSPORTING
Abstract
A method and an apparatus for beam analysis in an optical system are disclosed, wherein a plurality of beam parameters of a beam propagating along an optical axis are ascertained. The method includes: splitting the beam into a plurality of partial beams which have a focus offset in the longitudinal direction in relation to the optical axis; recording a measurement image produced by these partial beams; carrying out a forward simulation of the beam in the optical system on the basis of estimated initial values for the beam parameters in order to obtain a simulated image; and calculating a set of values for the beam parameters on the basis of the comparison between the simulated image and the measurement image.
Claims
1. An apparatus, comprising: a beam-splitting optical arrangement configured to split a beam incident on the beam-splitting optical arrangement along an optical axis into a plurality of partial beams having a focus offset in a longitudinal direction relative the optical axis; a sensor arrangement configured to capture an image of the partial beams; and a controller configured to: a) perform a forward simulation of the beam in the optical system based on estimated initial values for beam parameters to obtain a simulated image; and b) calculate a set of values for the beam parameters based on a comparison between the simulated image and the measurement image, wherein the controller is configured to: c) iteratively perform a) and b), wherein, in each case, the calculated set of values for the beam parameters provide the basis of the forward simulation that following; and d) output the sets of values for the beam parameters ascertained in c).
2. The apparatus of claim 1, wherein c) reduces a number of varied beam parameters.
3. The apparatus of claim 1, wherein the controller is configured to vary an algorithm used during c).
4. The apparatus of claim 1, wherein the apparatus is configured to record a near-field image produced by the beam.
5. The apparatus of claim 1, wherein the apparatus is configured to simultaneously record: i) a near-field image produced by the beam; and ii) a far-field image corresponding to the measurement image produced by the partial beams.
6. The apparatus of claim 1, wherein the plurality of beam parameters comprises at least one parameter selected from the group consisting of beam size, beam decentration, beam inclination, beam divergence, astigmatism, coma, and spherical aberration.
7. The apparatus of claim 1, wherein the apparatus is configured so that, during use of the apparatus, the beam is manipulated on the basis of the output values for the beam parameters while adapting at least one of the beam parameters.
8. The apparatus of claim 7, wherein the apparatus is configured so that, during use of the apparatus, outputting the output values and manipulating the beam occur in real time.
9. The apparatus of claim 1, wherein the beam-splitting optical arrangement is configured provide spherical wave-front deformations of the beam.
10. The apparatus of claim 9, wherein the beam-splitting optical arrangement comprises a diffractive structure.
11. The apparatus of claim 9, wherein the beam-splitting optical arrangement is configured so that points of incidence of the partial beams define a two-dimensional, grid-like distribution on a plane extending transversely to a light propagation direction of the beam.
12. The apparatus of claim 10, wherein the beam-splitting optical arrangement comprises two diffractive structures extending in mutually different directions.
13. The apparatus of claim 12, wherein the diffractive structures differ by at least a factor of three in terms of their focal length related to the first positive order of diffraction in each case.
14. The apparatus of claim 1, wherein the beam comprises a laser beam.
15. The apparatus of claim 1, further comprising a laser plasma source.
16. A beam-splitting arrangement, comprising: a first periodic diffractive structure extending in a first direction; and a second period diffractive structure extending in a second direction different from the first direction, wherein the beam-splitting optical arrangement is configured to split a beam incident thereon along an optical axis into a plurality of partial beams having a focus offset in a longitudinal direction relative to the optical axis so that points of incidence of the partial beams define a two-dimensional, grid-like distribution on a plane extending transversely to a propagation direction of the beam.
17. The beam-splitting arrangement of claim 16, wherein the beam-splitting optical arrangement brings about spherical wave-front deformations of the beam.
18. The beam-splitting arrangement of claim 16, wherein the diffractive structures differ by at least a factor of three in terms of their focal length related to the first positive order of diffraction in each case.
19. The beam-splitting arrangement of claim 16, wherein the beam is a laser beam.
20. The beam-splitting arrangement of claim 16, further comprising a laser plasma source.
21. The beam-splitting arrangement of claim 18, wherein the second direction is perpendicular to the first direction.
22. The beam-splitting arrangement of claim 16, wherein the second direction is perpendicular to the first direction.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1) The disclosure is explained in greater detail below on the basis of exemplary embodiments illustrated in the accompanying figures, in which:
(2)
(3)
(4)
(5)
(6)
(7)
(8)
DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS
(9) According to
(10) In the exemplary embodiment, as indicated in
(11)
(12) Here, f.sub.1 denotes the focal length of the first positive order of diffraction and k denotes the beam index or the order of diffraction. Here, the intensity of the respective focus depends directly on the embodiment and approximation form of the underlying (approximately parabolic) phase profile. In combination with a refractive lens element with a focal length of f.sub.0, a multi-focal optical system emerges with a plurality of used focal lengths f.sub.k, k=0, 1, . . . , k.sub.max, wherein the following applies approximately if the distance between the diffractive structure and the refractive lens element is neglected:
(13)
(14) This relation is elucidated in
(15) The disclosure is not restricted to the configuration of the beam-splitting optical arrangement 110 with such a diffractive structure. Rather, what is important in the configuration of the beam-splitting optical arrangement is that it causes, where possible, spherical wavefront deformations of the beam that is incident on the beam-splitting optical arrangement. In further embodiments, use may also be made of a different beam-spitting optical arrangement suitable to this end, for example in the form of an etalon.
(16) The partial beams emanating from the beam splitting optical arrangement are thereupon incident onreference once again being made to
(17)
(18) In principle, the recording of these individual spot images assigned to different focus positions in each case by the application of known, so-called phase retrieval methods (e.g. Gerchberg-Saxton algorithm) would allow a back calculation to the phase of the wavefront if the individual spot images were independent of one another (i.e. if there were no mutual influencing by way of interference). However, unavoidable interferences between the individual spot images are present in this case, the interferences leading to a pronounced mutual disturbance as indicated in
(19) Mathematically, these circumstances mean that no unique back transformation is possible for directly calculating the beam parameters. In order to take account of this problem, iterative comparisons between respectively calculated or simulated images and the recorded measurement image are performed according to the disclosure in a model-based approach, as described below with reference to
(20) As indicated in the schematic diagram of
(21) These parameters for describing the beam may be, for example, the beam size, beam decentration in the x-direction, beam decentration the y-direction, beam inclination in the x-direction, beam inclination in the y-direction, beam divergence, astigmatism in the x-direction, astigmatism in the y-direction, coma in the x-direction, coma in the y-direction and spherical aberration. Here, a Zernike parameterization may also be effectuated when desired in order to describe and ascertain corresponding wavefront aberrations of higher order too.
(22) Thereupon there is a forward simulation (step S420) for ascertaining a calculated image. According to
(23) If the assumption of a beam propagation in the positive z-direction is made, the beam amplitude to be determined (in the range of scalar diffraction) is denoted by u(x, y|z.sub.0) at the location z.sub.0 in the reference plane (ideally near-field plane). After passing over the free space path between the reference plane and the plane of the effective optical element (beam-splitting optical arrangement 110 or diffractive optical structure), the amplitude present at the entrance of the optical element forming the beam-splitting optical arrangement 110 or 310 at the position z.sub.1 is given by
u(x,y|z.sub.1)={circumflex over (P)}.sub.1u(x,y|z.sub.0)=IFT.sub.xy[(z.sub.1z.sub.0).Math.FT.sub.xy[u(x,y|z.sub.0)]](3).
(24) The optical element forming the beam-splitting optical arrangement 110 or 310, in the approximation of the infinitely thin element, multiplicatively impresses the amplitude function T(x,y)=t(x,y).Math.exp(i(x,y))u(x, y|z.sub.1) according to
u(x,y|z.sub.1+)=T(x,y)u(x,y|z.sub.1)(4).
(25) By way of a further free space propagation from the optical element forming the beam-splitting optical arrangement 110 or 310 to the sensor arrangement 120 or 320 (the plane of which lies perpendicular to the z-axis at the position z.sub.2), the amplitude on the plane of the sensor arrangement 120 or 320 is finally arrived at according to
u(x,y|z.sub.2)={circumflex over (P)}.sub.2u(x,y|z.sub.1+)=IFT.sub.xy[(z.sub.2z.sub.1).Math.FT.sub.xy[u(x,y|z.sub.1+)]](5)
(26) The intensity profile detected at the spatially-resolving sensor arrangement 120 or 320 is obtained by forming the square of the absolute value according to
I.sub.Sensor(x,y)=|u(x,y|z.sub.2)|.sup.2(6)
(27) The propagator of the free space propagation is known from the formalism of Fourier optics. During the propagation from a plane perpendicular to the z-axis at the position z to a parallel plane at the position z, the amplitude is initially transformed into the frequency space according to
(f.sub.x,f.sub.y|z)=FT.sub.xy[u(x,y|z)]=dxdy u(x,y|z)exp(2i(f.sub.xx+f.sub.yy))(7)
(28) by way of the 2D Fourier transform and there it is multiplied by the free space propagation function
(29)
(30) over the distance d=zz. Here, the phase in the propagation function is given by
(f.sub.x,f.sub.y)={square root over (1.sup.2(f.sub.x.sup.2+f.sub.y.sup.2))}(9),
(31) where f.sub.x, f.sub.y denote the spatial frequencies and denotes the wavelength of the radiation. The amplitude in the plane at z in the spatial domain is finally obtained by a back transformation by way of the inverse Fourier transform according to
(32)
(33) The correspondingly calculated image (containing the calculated intensity values I.sub.calc) is subtracted from the recorded measurement image (containing the measured intensity values I.sub.meas), whereupon appropriately modified model parameters for describing the beam are ascertained and these form the basis of a new forward simulation (step S460 in
(34) As already explained above, it is particularly advantageous if the absolute value of the amplitude is available from the near-field measurement and if it need not be described and fitted by a model. As a result, firstly, the number of parameters to be described is reduced, possibly significantly reduced, and, secondly, the quality of the information obtained about the beam to be measured is improved.
(35) According to the disclosure, the fact that enabling a large number of parameters leads to high numerical complexity is preferably further taken into account. Consequently, there preferably is initially a start with a comparatively small scope of the parameters set, which is then successively expanded with respect to the simultaneously varied parameters of the parameters set, i.e. adaptive fitting of the model is undertaken. Thus, for instance, if twenty parameters are sought after in principle, only ten dominating parameters may initially be enabled.
(36) Furthermore, it is also possible to adapt the respective evaluation method or the algorithm in order to obtain the fastest possible speed for the beam analysis, for example after reaching a quasi-stationary operation of the respective system (e.g. a largely stably operated plasma light source), in which, typically, only small changes in the beam parameters still occur. Here, in particular, use can be made of the already collected information in order then to be able to determine and correct, in real time, the small changes in the beam parameters that still occur. In this phase, the originally nonlinear optimization problem may also be approximable in linear form. As a result, what may be achieved thus is that, for instance in a plasma light source, the laser beam can be guided, accurately and quickly at the same time, with respect to the beam parameters.
(37)
(38) In order to obtain this two-dimensional beam split, it is possiblewithout the disclosure being restricted heretoto provide, for example in the embodiment schematically illustrated in
(39) The mode of operation of the beam-splitting optical arrangement according to
(40) Proceeding from the basic construction of the beam-splitting optical arrangement according to
(41) Analogous to the embodiment of
(42)
(43) Here, f.sub.1* and f.sub.2* denote (in relation to the respective first positive order of diffraction) the respective focal lengths of the first diffractive optical element 911 and the second diffractive optical element 912, respectively, and f.sub.0 denotes the focal length of the refractive optical element 913, while m and n denote the orders of diffraction of the respective diffraction at the first optical element 911 and second diffractive optical element 912, respectively.
(44) The focal lengths of the first diffractive optical element 911 and the second diffractive optical element 912 are selected to be different from one another, with the consequence that the element with the relatively shorter focal length produces the relatively greater longitudinal focus offset, and vice versa. In a specific exemplary embodiment, it is possible, for instance, for the focal length of the first diffractive optical element 911 to be greater by a factor of five than the focal length of the second diffractive optical element 912.
(45) In the case of a suitable selection of the aforementioned parameters (i.e. the focal lengths f.sub.1*, f.sub.2* and f.sub.0) and of the used value ranges of the orders of diffraction (n, m), it is now possible to obtain and use a measurement range increase in different ways, as will be explained below with reference to
(46)
(47) The above-described measurement range increase can be used to take account of the large focus variations of the respective beam to be characterized, as occur, for example in applications of material processing, in particular at high laser powers, as a consequence of heating and deformation of the individual optical components, namely by virtue of the capture region of the respective focus values being significantly increased (for example, by approximately a factor of seven according to
(48) In further embodiments (some of which are schematically illustrated in
(49) According to
(50) Even though the disclosure has been described on the basis of specific embodiments, numerous variations and alternative embodiments are apparent to a person skilled in the art, for example by combination and/or exchange of features of individual embodiments. Accordingly, it goes without saying for a person skilled in the art that such variations and alternative embodiments are concomitantly encompassed by the present disclosure, and the scope of the disclosure is restricted only within the meaning of the accompanying patent claims and the equivalents thereof.