Reconstruction of a wavefront of a light beam containing optical vortices

Patent metadata

Field	Value
Publication number	US-12510416-B2
Application number	US-202118019407-A
Country	US
Kind code	B2
Filing date	Jul 23, 2021
Priority date	Aug 4, 2020
Publication date	Dec 30, 2025
Grant date	Dec 30, 2025

Abstract

A method for reconstructing the wavefront of a light beam by analyzing wavefront-gradient data of said light beam, the light beam containing at least one optical vortex, considering the contribution of the optical vortices to the wavefront. The method including providing a phase-gradient map g of the wavefront of the light beam, generating a Laplacian of a vector potential based on the phase gradient map g, the resulting Laplacian of the vector potential map, called “Laplacian map”, exhibiting peaks, the location of each peak corresponding to the location of an optical vortex and the integral of the peak being proportional to a topological charge n of said optical vortex, computing a singular phase map φ s based on the topological charge n and location of each optical vortex, the singular phase φ s map being representative of the contribution of the optical vortex.

First claim

The invention claimed is:

1 . A method for reconstructing the wavefront of a light beam by analyzing wavefront-gradient data of said light beam, the light beam containing at least one optical vortex (V), considering the contribution of the optical vortices (V) to the wavefront, the method comprising:
a) providing a phase-gradient map g of the wavefront of the light beam, b) generating a Laplacian of a vector potential based on the phase gradient map g, the resulting Laplacian of the vector potential map (−ΔA z ), called “Laplacian map”, exhibiting peaks (P), the Laplacian map being obtained by applying a curl operator to the phase-gradient map g, the location of each peak (P) corresponding to the location (u) of an optical vortex (V) and the integral of the peak being proportional to a topological charge n of said optical vortex (V), c) applying an image processing method to the Laplacian map to extract the location (u) and the charge n of optical vortices, the image processing method comprising a segmentation of the Laplacian map into a plurality of segments (Ω) to obtain a segmented Laplacian map, segments (Ω) being determined so that different peaks of the Laplacian map belong to different segments; d) for these different segments (Ω) of the segmented Laplacian map (−ΔA z ):
i) determining the location (u) of the corresponding peak (P) in the segment (Ω) of the Laplacian map (−ΔA z ) to locate the corresponding optical vortex (V),
ii) computing an approximate value of the integral of the corresponding peak (P) over the segment (Ω), and determining the topological charge n of the corresponding vortex (V) based on said approximate value of the integral of said peak (P);
e) computing a singular phase map φ s based on the topological charge n and location (u) of each optical vortex (V) computed in steps i) and ii), the singular phase φs map being representative of the contribution of the optical vortex and being computed by:
iii) generating a vortex-location map, where map values are equal to the corresponding topological charges n of optical vortices at optical vortices locations (u), and
iv) convoluting said vortex-location map with a single +1 spiral phase profile;
f) reconstructing the wavefront of the light beam with:
computing a regular phase φ R map of the light beam, the regular phase corresponding to a vortex-free phase computed from the irrotational component,
adding the regular phase φ R to the singular phase φ s to obtain the said reconstructed wavefront of the light beam.

2 . The method according to claim 1 , wherein in step ii) the topological charge n of an optical vortex is obtained by:
dividing the approximate value of the integral of the corresponding peak by 2kπ, with k an integer rounding the result of the former operation to the closest integer to get the topological charge n.

3 . The method according to claim 2 , wherein the integer k is equal to −1 or +1.

4 . The method according to claim 1 , wherein the segmentation at step c) is realized with a numerical tool among the list: the watershed algorithm, the normalized cut algorithm, the SLIC algorithm, or convolutional neuronal networks, and/or a deep learning algorithm.

5 . The method according to claim 1 , step i) comprising computing the weighted centroid of each peak (P) in each segment (Ω) corresponding to said peak, determining the location (u) of the said peak based on the weighted centroid.

6 . The method according to claim 1 , wherein in step c) the image processing method comprises fitting the peaks (P) with bell-shaped functions, for example with Gaussian functions.

7 . The method according to claim 6 , wherein in step ii), the integral of the peak is estimated by computing the area under the corresponding fitting bell-shaped curve.

8 . The method according to claim 1 , wherein in step c) the image processing method comprising solving an inverse problem to estimate the distribution of the peaks (P), by using tools such as a compressed sensing algorithm, a Markov Chain Monte Carlo Algorithm (MCMC) methods or neuronal networks.

9 . The method according to claim 1 , the image processing method comprising identifying the extrema values of the Laplacian, notably using a threshold to ignore noise, wherein in step i) the location of the peaks is determined based on the location of the extrema.

10 . The method according to claim 1 , wherein in step ii) the integral of the peaks is determined from the extrema magnitude only.

11 . The method according to claim 1 , wherein the convolution in step e)iv) to obtain the singular phase φ s is: achieved numerically by adding as many spiral phase profiles as the number of detected optical vortices, each spiral phase profile being centered at the same location and having the same topological charge as the detected optical vortices (V).

12 . The method according to claim 1 , wherein in the step f), the regular phase φ R map is computed by applying a divergence operator to the phase gradient map g, followed by a double spatial integration.

13 . The method according to claim 1 , wherein the phase-gradient g is decomposed according to Helmholtz decomposition as follows: g=∇φ r +∇×A where φr is the regular phase, the gradient ∇φr corresponding to the irrotational component of g and where A is the potential vector corresponding to the solenoidal component of g.

14 . A computer-readable storage medium comprising instructions which, when executed by a computer, causes the computer to carry out the steps of the method according to claim 1 .

15 . A device for reconstructing the wavefront of a light beam by analyzing wavefront-gradient data of the light beam, the light beam containing at least one optical vortex (V), in order to estimate the contribution of the optical vortices (V) to the wavefront, comprising:
a wavefront sensor to measure the phase-gradient of the wavefront of the light beam; coupled to the wavefront sensor, means for carrying out the steps of the method of claim 1 , some calculating means comprising a memory storage unit, a set of instructions and a processor for executing said set of instructions, these instructions being: a) providing a phase-gradient map g of the wavefront of the light beam; b) generating a Laplacian of a vector potential based on the phase gradient map g, the resulting Laplacian of the vector potential map (−ΔA z ), called “Laplacian map”, exhibiting peaks (P), the Laplacian map being obtained by applying a curl operator to the phase-gradient map g,
the location of each peak (P) corresponding to the location (u) of an optical vortex (V) and the integral of the peak being proportional to a topological charge n of said optical vortex (V);
c) applying an image processing method to the Laplacian map to extract the location (u) and the charge n of optical vortices, the image processing method comprising a segmentation of the Laplacian map into a plurality of segments (Ω) to obtain a segmented Laplacian map, segments (Ω) being determined so that different segments contain different peaks of the Laplacian map;
d) for these different segments (Ω) of the segmented Laplacian map (−ΔA z ):
i) determining the location (u) of the corresponding peak (P) in the segment (Ω) of the Laplacian map (−ΔA z ) to locate the corresponding optical vortex (V),
ii) computing an approximate value of the integral of the corresponding peak (P) over the segment (Ω), and determining the topological charge n of the corresponding vortex (V) based on the said approximate value of the integral of the said peak (P);
e) computing a singular ph

Assignees

• Univ Paris Cite
• Centre Nat Rech Scient
• Univ Sorbonne
• Inst Nat Sante Rech Med

Inventors

• Guillon Marc
• Berto Pascal

Classifications

• G01J2009/002

  Wavefront phase distribution · CPC title

• G01J1/4257

  applied to monitoring the characteristics of a beam, e.g. laser beam, headlamp beam (monitoring arrangements for lasers in general H01S3/0014) · CPC title

• G01J9/0215 Primary

  by shearing interferometric methods · CPC title

• G01J9/00 Primary

  Measuring optical phase difference (devices or arrangements for controlling the phase of light beams G02F1/01); Determining degree of coherence; Measuring optical wavelength (spectrometry G01J3/00) · CPC title