Toward high precision position control using laser interferometry: main sources of error. DUCOURTIEUX SΓ©bastien Nanometrology, Laboratoire National de MΓ©trologie et dβEssais (LNE)
[email protected]
28th of June 2018
INTRODUCTION Position and displacement sensors are one of the key elements when designing a positioning system or a scanning device. In many applications, the performances of the instrument are dependent on the position sensor performances and on the way it is integrated into the machine. Many position sensors are available on the market, with many technologies. This is the reason why the choice is not easy to do and the best compromise not always easy to reach when designing a new setup. But when the development request severe constrains as for example: β’ Accuracy in the nanometer range ( bright fringe
Movable mirror
Photodetector
11
28th June 2018
π
PRINCIPLE OF MICHELSON INTERFEROMETER
reference arm
Mirror supposed to be fixed Displacement 50/50 beam splitter
Laser source measurement arm
Destructive interferences => dark fringe
Movable mirror
Photodetector
The electronic counts the number of fringes passing over the detector
12
28th June 2018
π
PRINCIPLE OF MICHELSON INTERFEROMETER Mirror supposed to be fixed reference arm
Displacement
π
50/50 beam splitter
Laser source measurement arm
Movable mirror
Electromagnetic wave propagating in the reference arm: πΈπ (π‘, πππΏπ ) = πΈ0
2π π ππ‘β πππΏπ ππ£ π
Electromagnetic wave propagating in the measurement arm: πΈπ (π‘, πππΏπ ) = πΈ0 13
2π π ππ‘β πππΏπ ππ£ π
28th June 2018
πΈ0 : electric field amplitude π: angular frequency πππΏ: optical path length ππ£ : vacuum wavelength π‘: time
PRINCIPLE OF MICHELSON INTERFEROMETER Mirror supposed to be fixed reference arm
Displacement 50/50 beam splitter
Laser source measurement arm
Movable mirror
Intensity on the photodetector : πΌ = πΈπ + πΈπ
πΌ = πΈ0
14
2π π ππ‘β πππΏπ ππ£ π
2π π ππ‘β πππΏπ ππ£ +π
Γ
2
2π βπ ππ‘β πππΏπ ππ£ π
28th June 2018
+
2π βπ ππ‘β πππΏπ ππ£ π
π
PRINCIPLE OF MICHELSON INTERFEROMETER
reference arm
Mirror supposed to be fixed Displacement 50/50 beam splitter
Laser source measurement arm
Movable mirror
Intensity on the photodetector : πΌ = πΈπ + πΈπ
πΌ = 2πΈ0 2 + 2πΈ0 2 cos
15
2
2π πππΏπ β πππΏπ ππ£
28th June 2018
π
PRINCIPLE OF MICHELSON INTERFEROMETER
reference arm
Mirror supposed to be fixed Displacement
π
50/50 beam splitter
Laser source measurement arm
Movable mirror
Intensity on the photodetector : πΌ = 2πΈ0 2 cos
2π πππΏπ β πππΏπ ππ£
As the interferometer counts some fringes, only this part as to be taken into account 16
28th June 2018
PRINCIPLE OF MICHELSON INTERFEROMETER Counting : π 2 = πππΏπ β πππΏπ π π 2 π£
ππ£ : vacuum wavelength π : displacement πππΏ: optical path length k: counting i: interpolation factor
1 - Interferometer are sensitive to the evolution of the optical path length difference between the two arms. 2 - It means that everything that will unbalance the difference between the two optical paths will produce a signal that will be interpreted as a displacement !!!
πππΏ =
π π ππ
π : beam path π: refractive index of air
π
Assumption 1: if the refractive index of air is constant along the beam path then: πππΏ = π Γ πΏπ
17
28th June 2018
π : beam path πΏπ : length of the beam path π: refractive index of air
PRINCIPLE OF MICHELSON INTERFEROMETER
reference arm πΏπ
Mirror supposed to be fixed Displacement
π
50/50 beam splitter
Laser source measurement arm πΏπ
Movable mirror
Counting: π 2 = 2 ππ π + 2 ππ πΏπ β 2ππ πΏπ 2π ππ£
ππ£ : vacuum wavelength ππ ππ ππ : refractive index πΏπ ππ πΏπ : physical length π : displacement k: counting i: interpolation factor
Assumption 2: the refractive index of air in the two arms is homogeneous (ππ = ππ = π). 18
28th June 2018
PRINCIPLE OF MICHELSON INTERFEROMETER
reference arm πΏπ
Mirror supposed to be fixed Displacement
π
50/50 beam splitter
Laser source measurement arm πΏπ
Movable mirror
Counting: π 4π = π + πΏπ β πΏπ 2π ππ£
ππ£ : vacuum wavelength π: refractive index πΏπ ππ πΏπ : physical length π : displacement k: counting i: interpolation factor
Interferometers measure a displacement, not an absolute position !! ο the initial positions of the mirrors are never known, only the relative change in position is determined
19
28th June 2018
PRINCIPLE OF MICHELSON INTERFEROMETER At π‘ = 0, interferometer is reseted: π0 4π0 = πΏπ β πΏπ = 0 π π 2 π£ At t and with π = π0 + βπ:
π 4π0 = ππ£ 2π
π 4 π0 + βπ = π + πΏπ β πΏπ ππ£ 2π 4βπ 4 π0 + βπ π πΏπ β πΏπ + πΏπ β πΏπ + ππ£ ππ£ =0
π 4 π0 + βπ π 4βπ = + πΏπ β πΏπ π π π 2 π£ π£ π=
πππ£ βπ πΏπ β πΏπ β π0 + βπ 4 Γ π0 + βπ Γ 2π π=
20
πππ£ βπ β πΏ β πΏπ π π 4π2π 28th June 2018
PRINCIPLE OF MICHELSON INTERFEROMETER
πππ£ βπ π= β πΏ β πΏπ π π 4π2π Displacement term
dead path term
1 β the vacuum wavelength ππ£ needs to be known (calibrated) to determine π 2 β in air, π and βπ also need to be determined and tracked to properly measure π
3 β for dead path term, difference of length between the measurement and reference arms (πΏπ β πΏπ ) need to be known 4 β in vacuum, π = 1 and simplified:
βπ = 0: the measurement using interferometry is greatly π=
21
πππ£ 4 Γ 2π 28th June 2018
WAVELENGTH CALIBRATION
22
28th June 2018
VACUUM WAVELENGTH CALIBRATION AND TRACEABILITY TO SI METER π=
πππ£ βπ β πΏ β πΏπ π π 4π2π The laser source vacuum wavelength ππ£ππ need to be calibrated if measurement need to be comparable and traceable to meter definition (SI)
Laser wavelength Displacement
Calibration of a laser frequency (wavelength) using a frequency beat comparison setup. The operation is generally directly realized by the supplier but can also be provided by a metrology institute (LNE, PTB, NIST, NMIJ, β¦) It gives you ππ£ with an associated uncertainty (most of time the wavelength stability of the laser source). LNEβs beat comparison setup
23
28th June 2018
VACUUM WAVELENGTH UNCERTAINTY Uncertainty in wavelength produces uncertainty in displacement measurement. Uncertainty mainly originates from the laser stability in wavelength (or frequency). π = 632,8 nm for He-Ne sources with an uncertainty typically comprised between 1-60 femtometers (1 fm = 10-15 m). It represents an equivalent stability in part per billion: 1 to 100 ppb. Displacement
Depending on the laser source and the time frame
Laser source stability
100 Β΅m
1 mm
10 mm
100 mm
400 mm
1000 mm
1 ppb
0,0001 nm
0,001 nm
0,01 nm
0,1 nm
0,4 nm
1 nm
10 ppb
0,001 nm
0,01 nm
0,1 nm
1 nm
4 nm
10 nm
100 ppb
0,01 nm
0,1 nm
1 nm
10 nm
40 nm
100 nm
ο Need to check the laser source stability οThe higher the displacements are, the more stable the source must be !
24
28th June 2018
VARIATION OF REFRACTIVE INDEX OF AIR
25
28th June 2018
VARIATION OF REFRACTIVE INDEX OF AIR
πππ£ βπ π= β πΏ β πΏπ π π 4π2π Changes of the environmental condition during the measurement are the largest source of error in interferometry.
π is sensitive to variations of: β’ temperature of air (T) β’ pressure (P) β’ relative humidity (%RH) β’ sometime even exact gas composition (for very precise applications) 26
28th June 2018
VARIATION OF REFRACTIVE INDEX OF AIR Sensitivity of refractive air index to:
β’
Temperature (T) ππ πΎπ = = β9.56 Γ 10β7 πΎ β1 β
1 πππ πΎ β1 ππ
β’
Pressure (P) πΎπ =
β’
ππ = 2.68 Γ 10β7 βππβ1 β
0.27 πππ βππβ1 ππ
Relative humidity (%RH) πΎπ»π
=
27
ππ = β8.5 Γ 10β9 %π
π»β1 β
0,01 πππ % π
π» β1 ππ
π»
28th June 2018
VARIATION OF REFRACTIVE INDEX OF AIR Temperature:
Displacement Temperature variation
100 Β΅m
1 mm
10 mm
100 mm
400 mm
1000 mm
0.001 Β°C
0,0001 nm
0.001 nm
0.01 nm
0.1 nm
0.4 nm
1 nm
0.01 Β°C
0,001 nm
0.01 nm
0.1 nm
1 nm
4 nm
10 nm
0.1 Β°C
0,01 nm
0.1 nm
1 nm
10 nm
40 nm
100 nm
1 Β°C
0,1 nm
1 nm
10 nm
100 nm
400 nm
1000 nm
Pressure:
28
Displacement Pressure variation
100 Β΅m
1 mm
10 mm
100 mm
400 mm
1000 mm
0.01 hPa
0,00027 nm
0.0027 nm
0.027 nm
0.27 nm
1 nm
2,7 nm
0.1 hPa
0,0027 nm
0.027 nm
0.27 nm
2.7 nm
10,8 nm
27 nm
1 hPa
0,027 nm
0.27 nm
2.7 nm nm
27 nm
108 nm
270 nm
10 hPa
0,27 nm
2.7 nm nm
27 nm nm
270 nm
1080 nm
2700 nm
28th June 2018
VARIATION OF REFRACTIVE INDEX OF AIR Hygrometry:
Displacement Hygrometry variation
100 Β΅m
1 mm
10 mm
100 mm
400 mm
1000 mm
0,01 %RH
0,00001 nm
0.0001 nm
0.001 nm
0.01 nm
0.04 nm
0,1 nm
0,1 %RH
0,0001 nm
0.001 nm
0.01 nm
0.1 nm
0.4 nm
1 nm
1 %RH
0,001 nm
0.01 nm
0.1 nm
1 nm
4 nm
10 nm
10 %RH
0,01 nm
0.1 nm
1 nm
10 nm
4 0nm
100 nm
ο Very sensitive to T,P and HR (depending on the displacement range) ο If possible, limit the variation of T, P and RH (enclosure around the instrument, control of the temperature and hygrometry of the room)
29
28th June 2018
VARIATION OF REFRACTIVE INDEX OF AIR
Clean room: β’ temperature 20 Β± 0.05 Β°C β’ hygrometry 50 Β± 5%RH Instruments use large amount of matter to improve thermal inertia and then stability
Thermal stability achieved in the instruments for the interferometer measurement: β’ hp-CMM: about 50 mK per 20 hours (impact due to βπ βΆ 15 ππ) β’ mAFM: about 1 mK per hour (impact due to βπ βΆ negligible)
30
28th June 2018
VARIATION OF REFRACTIVE INDEX OF AIR Hygrometry:
Displacement Hygrometry variation
100 Β΅m
1 mm
10 mm
100 mm
400 mm
1000 mm
0,01 %RH
0,00001 nm
0.0001 nm
0.001 nm
0.01 nm
0.04 nm
0,1 nm
0,1 %RH
0,0001 nm
0.001 nm
0.01 nm
0.1 nm
0.4 nm
1 nm
1 %RH
0,001 nm
0.01 nm
0.1 nm
1 nm
4 nm
10 nm
10 %RH
0,01 nm
0.1 nm
1 nm
10 nm
4 0nm
100 nm
ο Very sensitive to T,P and HR (depending on the displacement range) ο If possible, limit the variation of T, P and RH (enclosure around the instrument, thermal and hygrometry control of the room)
ο Implement some corrections to compensate for refractive air index variations when displacement is larger than 1 mm
31
28th June 2018
HOW TO APPLY CORRECTIONS Weather station : Measurement of temperature (T), pressure (P) and relative humidity (RH) using a calibrated weather station then correction using modified EdlΓ©n equation (empirical equation).
Edlen formulae uncertainty: 2 Γ 10β8 Uncertainty for Edlen correction: π’ π =
2 Γ 10β8
2
+ πΎπ π’ π
formulae uncertainty
32
2
+ πΎπ π’ π
2
+ πΎπ
π» π’ π
π»
T, P and RH measurement uncertainties
28th June 2018
2
HOW TO APPLY CORRECTIONS π’ π = 2 Γ 10β8 Numerical application:
2
+ πΎπ π’ π
2
+ πΎπ π’ π
2
+ πΎπ π’ π
π»
2
with π’ π = 0.1 πΎ, π’ π = 0.01βππ and π’ π
π» = 1 % β π’ π = 10β7 hp-CMM
mAFM
π = 300 ππ
π = 60 Β΅π
π’ π 33
π
= 30 ππ
π’ π
π
= 0,006 ππ
28th June 2018
HOW TO APPLY CORRECTIONS Weather station : Measurement of temperature (T), pressure (P) and relative humidity (RH) using a calibrated weather station then correction using modified EdlΓ©n equation (empirical equation).
β’ Easy to implement β’ Not really real-time (measurement of P, T, %HR are slow due to the low bandwidth of the sensor). β’ P, T, %HR need to be calibrated β’ Measurement are not really representative of the interferometer conditions because not realized at the same place. β’ Does not take in account π changes due gas composition. β’ T, P and RH must be measured with high precaution for high precision measurement 34
28th June 2018
HOW TO APPLY CORRECTIONS Wavelength tracker It tracks changes in the refractive index from the initial value (π‘ = 0 ) by using an interferometer with fixed physical path length. Principle :
βπππΏ = πππΏπ‘ β πππΏπ‘0 = 2ππ‘ Γ πΏππ‘ β πΏππ‘ β 2ππ‘0 Γ πΏππ‘0 β πΏππ‘0 If πΏπ and πΏπ are fixed (constant in time) βπππΏ β 2 Γ βnΓ πΏπ β πΏπ 35
28th June 2018
If πΏπ β πΏπ is known, βn is known
HOW TO APPLY CORRECTIONS Wavelength tracker Example of implementation on the hp-CMM
Fixed mirror
Metrology frame made of Invar
2 superposed interferometers to measurement one displacement, 1 for measurement of displacement, 1 for wavelength tracking.
36
28th June 2018
HOW TO APPLY CORRECTIONS Wavelength tracker Other solution proposed by Zygo : measurement beam passes through air, reference beam passes through a vacuum cell
βπππΏ = πππΏπ‘ β πππΏπ‘0 = 4πΏ Γ ππ‘0 β 1 β 4πΏ Γ π β 1 βπππΏ = 4πΏ Γ βn 37
28th June 2018
HOW TO APPLY CORRECTIONS Wavelength tracker It tracks changes in the refractive index from the initial value (π‘ = 0 ) by using an interferometer with fixed physical path length.
β’ High bandwidth compare to weather station β’ Need an additional interferometer with a good stability. β’ Need to determine the physical length of the optical path. β’ Based on the assumption that interferometer and wavelength tracker detect the same index changes ! Not true in perturbed environment. β’ Need to be placed very close from the interferometer to compensate. β’ Does not measure an absolute index. β’ Need to be completed with a weather station to fix π at π‘ = 0. β’ Uncertainty is limited by the uncertainty in the determination of π at π‘ = 0 by the weather station.
38
28th June 2018
DEAD PATH ERROR
39
28th June 2018
DEAD PATH ERROR π=
πππ£ βπ β πΏ β πΏπ π π 4π2π Zero location at π‘=0
reference arm πΏπ
Mirror supposed to be fixed
Beam splitter
Laser source measurement arm πΏπ reference arm πΏπ
dead path πΏπ·π
The dead path is the difference in optical path length between the reference and measurement beams when the interferometer electronic is zeroed.
40
28th June 2018
DEAD PATH ERROR π=
πππ£ βπ β πΏ β πΏπ π π 4π2π
reference arm πΏπ
Mirror supposed to be fixed
Zero location at π‘
Beam splitter
Laser source measurement arm πΏπ reference arm πΏπ
dead path πΏπ·π = πΏπ β πΏπ
The unbalance between the two arms (non symmetry) will produce a sensitivity to both changes in wavelength βπ or in air index βπ that will shift in space the zero location. The signal produced is then interpreted as a displacement even though the mirrors are immobile. 41
28th June 2018
DEAD PATH ERROR
42
28th June 2018
DEAD PATH ERROR
1. πΏπ·π must be reduced and if possible equal to zero 2. If not possible, the effect must be corrected
43
28th June 2018
REDUCTION OF DEAD PATH ERROR
πΏπ
Mirror supposed to be fixed
Laser source πΏπ
Zero location at π‘=0
The first solution to reduce the dead path error is to balance the optical path lengths of the two arms (symmetry) at π‘ = 0
If πΏπ = πΏπ then πΏπ·π = 0 44
28th June 2018
REDUCTION OF DEAD PATH ERROR One solution consists in using a differential interferometer where the two mirrors can be placed by the user in the configuration he wants: ο Possibility to balance reference and measurement path length reference arm
right angle prism
Reference mirror Measurement mirror
50/50 beam splitter
Laser source measurement arm
Relative displacement
Photodetector
45
28th June 2018
π
CORRECTION OF DEAD PATH ERROR Zero location at π‘
reference arm πΏπ
Mirror supposed to be fixed
Laser source measurement arm πΏπ dead path πΏπ·π
πΏπ·π = πΏπ β πΏπ If not possible due to the design of the instrument, then the solution is to compensate by applying a corrective term to the measurement:
π=
46
πππ£ βπ β πΏ π π·π 4π2π 28th June 2018
CORRECTION OF DEAD PATH ERROR Zero location at π‘
reference arm πΏπ
Mirror supposed to be fixed
Laser source measurement arm πΏπ dead path πΏπ·π
πΏπ·π = πΏπ β πΏπ Uncertainty of the correction: π’ππππ·π =
47
π’πΏπ·π
Ξπ π
2
+ πΏπ·π
π’Ξπ π
28th June 2018
2
CORRECTION OF DEAD PATH ERROR hp-CMM
πΏπ·π = 178 ππ
48
28th June 2018
Measurement mirror position at π‘ = 0
CORRECTION OF DEAD PATH ERROR πΏπ·π = 178 ππ π’πΏπ·π = 5 ππ
hp-CMM
Ξπ = 10β7 with π’Ξπ = 5 Γ 10β9 (equivalent to Ξπ = 0.1 πΎ ππ Ξπ = 0.4 βππ ππ Ξπ
π» = 10%) Without correction, dead path error = 17,8 nm
π’ππππ·π =
π’πΏπ·π
Ξπ π
2
+ πΏπ·π
π’ππππ·π = 1 nm
49
28th June 2018
π’Ξπ π
2
CORRECTION OF DEAD PATH ERROR mAFM Use of differential interferometers. Reference and measurement mirrors are in the same plane => no dead path generated outside the interferometer head.
50
28th June 2018
CORRECTION OF DEAD PATH ERROR mAFM Use of differential interferometers. Reference and measurement mirrors are in the same plane => no dead path generated outside the interferometer head. But unfortunately, instrument still suffers from dead path error because of the unbalance between reference and measurement paths inside the interferometer head.
51
28th June 2018
CORRECTION OF DEAD PATH ERROR mAFM Use of differential interferometers. Reference and measurement mirrors are in the same plane => no dead path generated outside the interferometer head. But unfortunately, instrument still suffers from dead path error because of the unbalance between reference and measurement paths inside the interferometer head (estimated to 15 Β± 1 mm and corrected). Without correction, dead path error = 1,5 nm With correction π’ππππ·π = 0,125 nm
52
28th June 2018
NON-SYMMETRY OF INTERFEROMETER
53
28th June 2018
NON-SYMMETRY OF INTERFEROMETER 2 massive aluminum enclosures
ο special attention should be given to the design of the interferometer because non symmetry of the reference and measurement paths can be a great source or errors.
Differential interferometer Only one plane mirror for both arms frame
Another funny experiment:
54
28th June 2018
NON-SYMMETRY OF INTERFEROMETER 2 massive aluminum enclosures
ο special attention should be given to the design of the interferometer because non symmetry of the reference and measurement paths can be a great source or errors.
Differential interferometer Only one plane mirror for both arms frame
Another funny experiment:
Variation of room pressure
160 nm
Pressure (hPa)
Position (nm)
Displacement detected interferometer
Time (hour)
Time (hour)
55
28th June 2018
NON-SYMMETRY OF INTERFEROMETER 2 massive aluminum enclosures
ο special attention should be given to the design of the interferometer because non symmetry of the reference and measurement paths can be a great source or errors.
Differential interferometer Only one plane mirror for both arms frame
Another funny experiment:
variation of room pressure
Position (nm)
Pressure (hPa)
Displacement detected interferometer
Drifts are presentβ¦
β¦ and correlated to room pressure variations (due to the unbalance between reference an measurement paths) Time (hour)
Time (hour)
56
28th June 2018
NON-SYMMETRY OF INTERFEROMETER 2 massive aluminum enclosures
ο special attention should be given to the design of the interferometer because non symmetry of the reference and measurement paths can be a great source or errors.
Differential interferometer Only one plane mirror for both arms frame
Another funny experiment:
Corrected measurements
ο differential thermal expansion between reference and measurement paths
Position (nm)
After correction of the dead path inside the interferometer head, the measurement exhibit a thermal drift !!!
ο can be assimilated to a drift of the interferometer in nm/hour Time (hour)
57
28th June 2018
AIR TURBULENCES
58
28th June 2018
AIR TURBULENCES Air turbulences can have a strong impact on the measurement and the stability of the positioning if no precaution are taken ! Realized with a differential interferometer and only one mirror
In clean room, no protection
59
28th June 2018
AIR TURBULENCES Air turbulences can have a strong impact on the measurement and the stability of the positioning if no precaution are taken ! Realized with a differential interferometer and only one mirror (only sensitive to refractive index of air variation between the two arms
In clean room, no protection
In clean room, with a metallic tube and aluminum foils to protect the beams. The system is enclosed in an aluminum box
60
28th June 2018
AIR TURBULENCES Air turbulences can have a strong impact on the measurement and the stability of the positioning if no precaution are taken ! Realized with a differential interferometer and only one mirror (only sensitive to refractive index of air variation between the two arms
ο Reduction of the turbulences effect from 15 nm to only 0.1 nm
In clean room, no protection
In clean room, with a metallic tube and aluminum foils to protect the beams, The system is enclosed in an aluminum box
61
ο Protection of the beams can improves measurement quality
28th June 2018
AIR TURBULENCES Origin of the effect : n2(T2,P2)
mirror
Displacement of nonhomogeneous (T, P) air into the beams
n1(T1,P1)
interferometer n3(T3,P3)
T1