1397 lines
60 KiB
HTML
1397 lines
60 KiB
HTML
<?xml version="1.0" encoding="utf-8"?>
|
||
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN"
|
||
"http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd">
|
||
<html xmlns="http://www.w3.org/1999/xhtml" lang="en" xml:lang="en">
|
||
<head>
|
||
<!-- 2022-06-22 Wed 13:29 -->
|
||
<meta http-equiv="Content-Type" content="text/html;charset=utf-8" />
|
||
<title>ESRF Double Crystal Monochromator - Metrology</title>
|
||
<meta name="author" content="Dehaeze Thomas" />
|
||
<meta name="generator" content="Org Mode" />
|
||
<link rel="stylesheet" type="text/css" href="https://research.tdehaeze.xyz/css/style.css"/>
|
||
<script type="text/javascript" src="https://research.tdehaeze.xyz/js/script.js"></script>
|
||
<script>
|
||
MathJax = {
|
||
svg: {
|
||
scale: 1,
|
||
fontCache: "global"
|
||
},
|
||
tex: {
|
||
tags: "ams",
|
||
multlineWidth: "%MULTLINEWIDTH",
|
||
tagSide: "right",
|
||
macros: {bm: ["\\boldsymbol{#1}",1],},
|
||
tagIndent: ".8em"
|
||
}
|
||
};
|
||
</script>
|
||
<script id="MathJax-script" async
|
||
src="https://cdn.jsdelivr.net/npm/mathjax@3/es5/tex-svg.js"></script>
|
||
</head>
|
||
<body>
|
||
<div id="org-div-home-and-up">
|
||
<a accesskey="h" href="../index.html"> UP </a>
|
||
|
|
||
<a accesskey="H" href="../index.html"> HOME </a>
|
||
</div><div id="content" class="content">
|
||
<h1 class="title">ESRF Double Crystal Monochromator - Metrology</h1>
|
||
<div id="table-of-contents" role="doc-toc">
|
||
<h2>Table of Contents</h2>
|
||
<div id="text-table-of-contents" role="doc-toc">
|
||
<ul>
|
||
<li><a href="#org12e0a2f">1. Metrology Concept</a>
|
||
<ul>
|
||
<li><a href="#orgd6943c0">1.1. Sensor Topology</a></li>
|
||
<li><a href="#orgbea7abe">1.2. Computation of the relative pose between first and second crystals</a></li>
|
||
</ul>
|
||
</li>
|
||
<li><a href="#org4ff5766">2. Relation Between Crystal position and X-ray measured displacement</a>
|
||
<ul>
|
||
<li><a href="#org7c7d2db">2.1. Definition of frame</a></li>
|
||
<li><a href="#org67fb689">2.2. Effect of an error in crystal’s distance</a></li>
|
||
<li><a href="#orgf515721">2.3. Effect of an error in crystal’s x parallelism</a></li>
|
||
<li><a href="#org2fd3524">2.4. Effect of an error in crystal’s y parallelism</a></li>
|
||
<li><a href="#org4ac158b">2.5. Summary</a></li>
|
||
<li><a href="#org86a50e2">2.6. “Channel cut” Scan</a></li>
|
||
</ul>
|
||
</li>
|
||
<li><a href="#orgc33f680">3. Determining relative pose between the crystals using the X-ray</a>
|
||
<ul>
|
||
<li><a href="#org97fbcd2">3.1. Determine the \(y\) parallelism - “Rocking Curve”</a></li>
|
||
<li><a href="#orga0df616">3.2. Determine the \(x\) parallelism - Bragg Scan</a></li>
|
||
<li><a href="#org01554e4">3.3. Determine the \(z\) distance - Bragg Scan</a></li>
|
||
<li><a href="#org9347512">3.4. Use Channel cut scan to determine crystal <code>dry</code> parallelism</a></li>
|
||
<li><a href="#org7a140bf">3.5. Effect of an error on Bragg angle</a></li>
|
||
</ul>
|
||
</li>
|
||
<li><a href="#orgf968de4">4. Deformations of the Metrology Frame</a>
|
||
<ul>
|
||
<li><a href="#org573b7e7">4.1. Measurement Setup</a></li>
|
||
<li><a href="#orgbabe313">4.2. Simulations</a></li>
|
||
<li><a href="#orgaf7afd6">4.3. Comparison</a></li>
|
||
<li><a href="#org99a64a4">4.4. Test</a></li>
|
||
<li><a href="#org907c706">4.5. Measured frame deformation</a></li>
|
||
<li><a href="#orgaa71cc1">4.6. Test</a></li>
|
||
<li><a href="#org73400be">4.7. Repeatability of frame deformation</a></li>
|
||
</ul>
|
||
</li>
|
||
<li><a href="#org0b40d1b">5. Attocube - Periodic Non-Linearity</a>
|
||
<ul>
|
||
<li><a href="#orge0a38ec">5.1. Measurement Setup</a></li>
|
||
<li><a href="#org821469c">5.2. Choice of the reference signal</a></li>
|
||
<li><a href="#org3263c1f">5.3. Repeatability of the non-linearity</a></li>
|
||
<li><a href="#org39fabc8">5.4. Simulation</a></li>
|
||
<li><a href="#org373814c">5.5. Measurements</a></li>
|
||
</ul>
|
||
</li>
|
||
<li><a href="#orga081f50">Bibliography</a></li>
|
||
</ul>
|
||
</div>
|
||
</div>
|
||
<hr>
|
||
<p>This report is also available as a <a href="./dcm-metrology.pdf">pdf</a>.</p>
|
||
<hr>
|
||
|
||
<p>
|
||
In this document, the metrology system is studied.
|
||
First, in Section <a href="#org8b7cd3e">1</a> the goal of the metrology system is stated and the proposed concept is described.
|
||
</p>
|
||
|
||
<p>
|
||
How the relative crystal pose is affecting the pose of the output beam is studied in Section <a href="#orga416efe">2</a>.
|
||
</p>
|
||
|
||
<p>
|
||
In order to increase the accuracy of the metrology system, two problems are to be dealt with:
|
||
</p>
|
||
<ul class="org-ul">
|
||
<li>The deformation of the metrology frame under the action of gravity (Section <a href="#org6019b13">4</a>)</li>
|
||
<li>The periodic non-linearity of the interferometers (Section <a href="#org2a9ae01">5</a>)</li>
|
||
</ul>
|
||
|
||
<div id="outline-container-org12e0a2f" class="outline-2">
|
||
<h2 id="org12e0a2f"><span class="section-number-2">1.</span> Metrology Concept</h2>
|
||
<div class="outline-text-2" id="text-1">
|
||
<p>
|
||
<a id="org8b7cd3e"></a>
|
||
</p>
|
||
<p>
|
||
The goal of the metrology system is to measure the distance and default of parallelism between the first and second crystals.
|
||
</p>
|
||
|
||
<p>
|
||
Only 3 degrees of freedom are of interest:
|
||
</p>
|
||
<ul class="org-ul">
|
||
<li>\(d_z\)</li>
|
||
<li>\(r_y\)</li>
|
||
<li>\(r_x\)</li>
|
||
</ul>
|
||
</div>
|
||
<div id="outline-container-orgd6943c0" class="outline-3">
|
||
<h3 id="orgd6943c0"><span class="section-number-3">1.1.</span> Sensor Topology</h3>
|
||
<div class="outline-text-3" id="text-1-1">
|
||
<p>
|
||
In order to measure the relative pose of the two crystals, instead of performing a direct measurement which is complicated, the pose of the two crystals are measured from a metrology frame.
|
||
Three interferometers are used to measured the 3dof of interest for each crystals.
|
||
Three additional interferometers are used to measured the relative motion of the metrology frame.
|
||
</p>
|
||
|
||
|
||
<p>
|
||
In total, there are 15 interferometers represented in Figure <a href="#org8523680">1</a>.
|
||
The measurements are summarized in Table <a href="#orga3424d2">2</a>.
|
||
</p>
|
||
|
||
<table id="orga1dc5ef" border="2" cellspacing="0" cellpadding="6" rules="groups" frame="hsides">
|
||
<caption class="t-above"><span class="table-number">Table 1:</span> Notations for the metrology frame</caption>
|
||
|
||
<colgroup>
|
||
<col class="org-left" />
|
||
|
||
<col class="org-left" />
|
||
</colgroup>
|
||
<thead>
|
||
<tr>
|
||
<th scope="col" class="org-left"><b>Notation</b></th>
|
||
<th scope="col" class="org-left"><b>Meaning</b></th>
|
||
</tr>
|
||
</thead>
|
||
<tbody>
|
||
<tr>
|
||
<td class="org-left"><code>d</code></td>
|
||
<td class="org-left">“Downstream”: Positive X</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-left"><code>u</code></td>
|
||
<td class="org-left">“Upstream”: Negative X</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-left"><code>h</code></td>
|
||
<td class="org-left">“Hall”: Positive Y</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-left"><code>r</code></td>
|
||
<td class="org-left">“Ring”: Negative Y</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-left"><code>f</code></td>
|
||
<td class="org-left">“Frame”</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-left"><code>1</code></td>
|
||
<td class="org-left">“First Crystals”</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-left"><code>2</code></td>
|
||
<td class="org-left">“Second Crystals”</td>
|
||
</tr>
|
||
</tbody>
|
||
</table>
|
||
|
||
<table id="orga3424d2" border="2" cellspacing="0" cellpadding="6" rules="groups" frame="hsides">
|
||
<caption class="t-above"><span class="table-number">Table 2:</span> List of Interferometer measurements</caption>
|
||
|
||
<colgroup>
|
||
<col class="org-right" />
|
||
|
||
<col class="org-left" />
|
||
|
||
<col class="org-left" />
|
||
</colgroup>
|
||
<thead>
|
||
<tr>
|
||
<th scope="col" class="org-right"><b>Number</b></th>
|
||
<th scope="col" class="org-left"><b>Measurement</b></th>
|
||
<th scope="col" class="org-left"><b>Description</b></th>
|
||
</tr>
|
||
</thead>
|
||
<tbody>
|
||
<tr>
|
||
<td class="org-right">1</td>
|
||
<td class="org-left">\(z_{1r,u}\)</td>
|
||
<td class="org-left">First “Ring” Crystal, “upstream”</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-right">2</td>
|
||
<td class="org-left">\(z_{1r,c}\)</td>
|
||
<td class="org-left">First “Ring” Crystal, “center”</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-right">3</td>
|
||
<td class="org-left">\(z_{1r,d}\)</td>
|
||
<td class="org-left">First “Ring” Crystal, “downstream”</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-right">4</td>
|
||
<td class="org-left">\(z_{1h,u}\)</td>
|
||
<td class="org-left">First “Hall” Crystal, “upstream”</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-right">5</td>
|
||
<td class="org-left">\(z_{1h,c}\)</td>
|
||
<td class="org-left">First “Hall” Crystal, “center”</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-right">6</td>
|
||
<td class="org-left">\(z_{1h,d}\)</td>
|
||
<td class="org-left">First “Hall” Crystal, “downstream”</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-right">7</td>
|
||
<td class="org-left">\(z_{2h,u}\)</td>
|
||
<td class="org-left">Second “Hall” Crystal, “upstream”</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-right">8</td>
|
||
<td class="org-left">\(z_{2h,c}\)</td>
|
||
<td class="org-left">Second “Hall” Crystal, “center”</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-right">9</td>
|
||
<td class="org-left">\(z_{2h,d}\)</td>
|
||
<td class="org-left">Second “Hall” Crystal, “downstream”</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-right">10</td>
|
||
<td class="org-left">\(z_{2r,u}\)</td>
|
||
<td class="org-left">Second “Ring” Crystal, “upstream”</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-right">11</td>
|
||
<td class="org-left">\(z_{2r,c}\)</td>
|
||
<td class="org-left">Second “Ring” Crystal, “center”</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-right">12</td>
|
||
<td class="org-left">\(z_{2r,d}\)</td>
|
||
<td class="org-left">Second “Ring” Crystal, “downstream”</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-right">13</td>
|
||
<td class="org-left">\(z_{mf,u}\)</td>
|
||
<td class="org-left">Metrology Frame, “upstream”</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-right">14</td>
|
||
<td class="org-left">\(z_{mf,dr}\)</td>
|
||
<td class="org-left">Metrology Frame, “downstream-ring”</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-right">15</td>
|
||
<td class="org-left">\(z_{mf,dh}\)</td>
|
||
<td class="org-left">Metrology Frame, “downstream-hall”</td>
|
||
</tr>
|
||
</tbody>
|
||
</table>
|
||
|
||
|
||
<div id="org8523680" class="figure">
|
||
<p><img src="figs/metrology_schematic.png" alt="metrology_schematic.png" />
|
||
</p>
|
||
<p><span class="figure-number">Figure 1: </span>Schematic of the Metrology System</p>
|
||
</div>
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-orgbea7abe" class="outline-3">
|
||
<h3 id="orgbea7abe"><span class="section-number-3">1.2.</span> Computation of the relative pose between first and second crystals</h3>
|
||
<div class="outline-text-3" id="text-1-2">
|
||
<p>
|
||
To understand how the relative pose between the crystals is computed from the interferometer signals, have a look at <a href="https://gitlab.esrf.fr/dehaeze/dcm-kinematics">this repository</a> (<code>https://gitlab.esrf.fr/dehaeze/dcm-kinematics</code>).
|
||
</p>
|
||
|
||
<p>
|
||
Basically, Jacobian matrices are derived from the geometry and are used to convert the 15 interferometer signals to the <b>relative pose</b> of the primary and secondary crystals \([d_{h,z},\ r_{h,y},\ r_{h,x}]\) or \([d_{r,z},\ r_{r,y},\ r_{r,x}]\).
|
||
</p>
|
||
|
||
<div class="note" id="org186878e">
|
||
<p>
|
||
The sign conventions for the relative crystal pose are:
|
||
</p>
|
||
<ul class="org-ul">
|
||
<li>An increase of \(d_{h,z}\) means the two crystals are further apart</li>
|
||
<li>An increase of \(r_{h,}\) means that the second crystals experiences a rotation around \(y\) with respect to the primary crystal</li>
|
||
<li>An increase of \(r_{h,x}\) means that the second crystals experiences a rotation around \(x\) with respect to the primary crystal</li>
|
||
</ul>
|
||
|
||
</div>
|
||
|
||
<p>
|
||
The relative pose can be expressed as a function of the interferometers using the Jacobian matrices for the “hall” crystals:
|
||
</p>
|
||
\begin{equation}
|
||
\boxed{
|
||
\begin{bmatrix} d_{h,z} \\ r_{h,y} \\ r_{h,x} \end{bmatrix}
|
||
=
|
||
\bm{J}_{2h,s}^{-1}
|
||
\begin{bmatrix} z_{2h,u} \\ z_{2h,c} \\ z_{2h,d} \end{bmatrix}
|
||
-
|
||
\bm{J}_{1h,s}^{-1}
|
||
\begin{bmatrix} z_{1h,u} \\ z_{1h,c} \\ z_{1h,d} \end{bmatrix}
|
||
-
|
||
\bm{J}_{mf,s}^{-1}
|
||
\begin{bmatrix} z_{mf,u} \\ z_{mf,dh} \\ z_{mf,dr} \end{bmatrix}
|
||
}
|
||
\end{equation}
|
||
|
||
<p>
|
||
As well as for the “ring” crystals:
|
||
</p>
|
||
\begin{equation}
|
||
\boxed{
|
||
\begin{bmatrix} d_{r,z} \\ r_{r,y} \\ r_{r,x} \end{bmatrix}
|
||
=
|
||
\bm{J}_{2r,s}^{-1}
|
||
\begin{bmatrix} z_{2r,u} \\ z_{2r,c} \\ z_{2r,d} \end{bmatrix}
|
||
-
|
||
\bm{J}_{1r,s}^{-1}
|
||
\begin{bmatrix} z_{1r,u} \\ z_{1r,c} \\ z_{1r,d} \end{bmatrix}
|
||
-
|
||
\bm{J}_{mf,s}^{-1}
|
||
\begin{bmatrix} z_{mf,u} \\ z_{mf,dr} \\ z_{mf,dr} \end{bmatrix}
|
||
}
|
||
\end{equation}
|
||
|
||
<p>
|
||
Values of the matrices can be found in the document describing the kinematics of the DCM (see <code>https://gitlab.esrf.fr/dehaeze/dcm-kinematics</code>).
|
||
</p>
|
||
</div>
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-org4ff5766" class="outline-2">
|
||
<h2 id="org4ff5766"><span class="section-number-2">2.</span> Relation Between Crystal position and X-ray measured displacement</h2>
|
||
<div class="outline-text-2" id="text-2">
|
||
<p>
|
||
<a id="orga416efe"></a>
|
||
</p>
|
||
<p>
|
||
In this section, the impact of an error in the relative pose between the first and second crystals on the output X-ray beam is studied.
|
||
</p>
|
||
|
||
<p>
|
||
This is very important in order to:
|
||
</p>
|
||
<ul class="org-ul">
|
||
<li>link a measurement of the x-ray beam position to a default in the crystal position</li>
|
||
<li>understand which pose default will have large impact on the output beam position/orientation</li>
|
||
<li>calibrate the deformations of the metrology frame using an external metrology measuring the x-ray beam position/orientation</li>
|
||
</ul>
|
||
|
||
<p>
|
||
In order to simplify the problem, the first crystal is supposed to be fixed (i.e. ideally positioned), and only the motion of the second crystal is studied.
|
||
</p>
|
||
|
||
<p>
|
||
In order to easily study that, “ray tracing” techniques are used.
|
||
</p>
|
||
</div>
|
||
<div id="outline-container-org7c7d2db" class="outline-3">
|
||
<h3 id="org7c7d2db"><span class="section-number-3">2.1.</span> Definition of frame</h3>
|
||
<div class="outline-text-3" id="text-2-1">
|
||
<table border="2" cellspacing="0" cellpadding="6" rules="groups" frame="hsides">
|
||
|
||
|
||
<colgroup>
|
||
<col class="org-left" />
|
||
|
||
<col class="org-left" />
|
||
|
||
<col class="org-left" />
|
||
</colgroup>
|
||
<thead>
|
||
<tr>
|
||
<th scope="col" class="org-left"> </th>
|
||
<th scope="col" class="org-left">Position</th>
|
||
<th scope="col" class="org-left">Orientation</th>
|
||
</tr>
|
||
</thead>
|
||
<tbody>
|
||
<tr>
|
||
<td class="org-left">input beam</td>
|
||
<td class="org-left"><code>x1,y1,z1</code></td>
|
||
<td class="org-left"><code>s1</code></td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-left">primary mirror</td>
|
||
<td class="org-left"><code>xp,yp,zp</code></td>
|
||
<td class="org-left"><code>np</code></td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-left">reflected beam</td>
|
||
<td class="org-left"><code>x2,y2,z2</code></td>
|
||
<td class="org-left"><code>s2</code></td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-left">secondary mirror</td>
|
||
<td class="org-left"><code>xs,ys,zz</code></td>
|
||
<td class="org-left"><code>ns</code></td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-left">output beam</td>
|
||
<td class="org-left"><code>x3,y3,z3</code></td>
|
||
<td class="org-left"><code>s3</code></td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-left">Dectector</td>
|
||
<td class="org-left"><code>xd,yd,zd</code></td>
|
||
<td class="org-left"><code>nd</code></td>
|
||
</tr>
|
||
</tbody>
|
||
</table>
|
||
|
||
<ul class="org-ul">
|
||
<li class="off"><code>[ ]</code> Add schematic</li>
|
||
</ul>
|
||
|
||
<p>
|
||
Notations for the pose of the output beam: <code>d_b_y</code>, <code>d_b_z</code>, <code>r_b_y</code>, <code>r_b_z</code>.
|
||
</p>
|
||
|
||
<p>
|
||
\(d_{b,y}, d_{b,z}, r_{b,y}, r_{b,z}\)
|
||
</p>
|
||
|
||
<p>
|
||
The xy position of the beam is taken in the \(x=0\) plane.
|
||
</p>
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-org67fb689" class="outline-3">
|
||
<h3 id="org67fb689"><span class="section-number-3">2.2.</span> Effect of an error in crystal’s distance</h3>
|
||
<div class="outline-text-3" id="text-2-2">
|
||
<p>
|
||
<a id="org491e37f"></a>
|
||
</p>
|
||
|
||
<p>
|
||
In Figure <a href="#org24baa60">2</a> is shown the light path for three bragg angles (5, 55 and 85 degrees) when there is an error in the <code>dz</code> position of 1mm.
|
||
</p>
|
||
|
||
<p>
|
||
Visually, it is clear that this induce a <code>z</code> offset of the output beam.
|
||
</p>
|
||
|
||
|
||
<div id="org24baa60" class="figure">
|
||
<p><img src="figs/ray_tracing_error_dz_overview.png" alt="ray_tracing_error_dz_overview.png" />
|
||
</p>
|
||
<p><span class="figure-number">Figure 2: </span>Visual Effect of an error in <code>dz</code> (1mm). Side view.</p>
|
||
</div>
|
||
|
||
<p>
|
||
The motion of the output beam is displayed as a function of the Bragg angle in Figure <a href="#org692379d">3</a>.
|
||
It is clear that an error in the distance <code>dz</code> between the crystals only induce a <code>z</code> offset of the output beam.
|
||
This offset decreases with the Bragg angle.
|
||
</p>
|
||
|
||
<p>
|
||
This is indeed equal to:
|
||
</p>
|
||
\begin{equation}
|
||
\boxed{d_{b,z} = 2 d_z \cos \theta}
|
||
\end{equation}
|
||
|
||
|
||
<div id="org692379d" class="figure">
|
||
<p><img src="figs/motion_beam_dz_error.png" alt="motion_beam_dz_error.png" />
|
||
</p>
|
||
<p><span class="figure-number">Figure 3: </span>Motion of the output beam with dZ error</p>
|
||
</div>
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-orgf515721" class="outline-3">
|
||
<h3 id="orgf515721"><span class="section-number-3">2.3.</span> Effect of an error in crystal’s x parallelism</h3>
|
||
<div class="outline-text-3" id="text-2-3">
|
||
<p>
|
||
<a id="orgb8a03c5"></a>
|
||
</p>
|
||
|
||
<p>
|
||
The effect of an error in <code>rx</code> crystal parallelism on the output beam is visually shown in Figure <a href="#org8f13c50">4</a> for three bragg angles (5, 55 and 85 degrees).
|
||
The error is set to one degree, and the top view is shown.
|
||
It is clear that the output beam experiences some rotation around a vertical axis.
|
||
The amount of rotation depends on the bragg angle.
|
||
</p>
|
||
|
||
|
||
<div id="org8f13c50" class="figure">
|
||
<p><img src="figs/ray_tracing_error_drx_overview.png" alt="ray_tracing_error_drx_overview.png" />
|
||
</p>
|
||
<p><span class="figure-number">Figure 4: </span>Visual Effect of an error in <code>drx</code> (1 degree). Top View.</p>
|
||
</div>
|
||
|
||
<p>
|
||
The effect of <code>drx</code> as a function of the Bragg angle on the output beam pose is computed and shown in Figure <a href="#orgd270718">5</a>.
|
||
</p>
|
||
|
||
<p>
|
||
It induces a rotation of the output beam in the <code>z</code> direction that depends on the Bragg Angle.
|
||
It is starting at zero for small bragg angles, and it increases with the bragg angle up to 2 times <code>drx</code>.
|
||
This is indeed equal to:
|
||
</p>
|
||
\begin{equation}
|
||
\boxed{r_{b,z} = -2 r_x \sin \theta}
|
||
\end{equation}
|
||
|
||
<p>
|
||
If also induces a small \(y\) shift of the beam.
|
||
This shift is due to the fact that the rotation point (around which the second crystal is moved) is changing as a function of bragg.
|
||
We can note that the \(y\) shift is equal to zero for a bragg angle of 45 degrees, at which point the center of rotation of the second crystal is at \(x = 0\);
|
||
</p>
|
||
|
||
|
||
<div id="orgd270718" class="figure">
|
||
<p><img src="figs/motion_beam_drx_error.png" alt="motion_beam_drx_error.png" />
|
||
</p>
|
||
<p><span class="figure-number">Figure 5: </span>Motion of the output beam with drx error</p>
|
||
</div>
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-org2fd3524" class="outline-3">
|
||
<h3 id="org2fd3524"><span class="section-number-3">2.4.</span> Effect of an error in crystal’s y parallelism</h3>
|
||
<div class="outline-text-3" id="text-2-4">
|
||
<p>
|
||
<a id="orgc411ea5"></a>
|
||
</p>
|
||
|
||
<p>
|
||
The effect of an error in <code>ry</code> crystal parallelism on the output beam is visually shown in Figure <a href="#org5651b67">6</a> for three bragg angles (5, 55 and 85 degrees).
|
||
</p>
|
||
|
||
|
||
<div id="org5651b67" class="figure">
|
||
<p><img src="figs/ray_tracing_error_dry_overview.png" alt="ray_tracing_error_dry_overview.png" />
|
||
</p>
|
||
<p><span class="figure-number">Figure 6: </span>Visual Effect of an error in <code>dry</code> (1 degree). Side view.</p>
|
||
</div>
|
||
|
||
<p>
|
||
The effect of <code>dry</code> as a function of the Bragg angle on the output beam pose is computed and shown in Figure <a href="#org94ce447">7</a>.
|
||
It is clear that this induces a rotation of the output beam in the <code>y</code> direction equals to 2 times <code>dry</code>:
|
||
</p>
|
||
\begin{equation}
|
||
\boxed{r_{b,z} = 2 r_z}
|
||
\end{equation}
|
||
|
||
<p>
|
||
It also induces a small vertical motion of the beam (at the \(x=0\) location) which is simply due to the fact that the \(x\) coordinate of the impact point on the second crystal changes with the Bragg angle.
|
||
</p>
|
||
|
||
|
||
<div id="org94ce447" class="figure">
|
||
<p><img src="figs/motion_beam_dry_error.png" alt="motion_beam_dry_error.png" />
|
||
</p>
|
||
<p><span class="figure-number">Figure 7: </span>Motion of the output beam with dry error</p>
|
||
</div>
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-org4ac158b" class="outline-3">
|
||
<h3 id="org4ac158b"><span class="section-number-3">2.5.</span> Summary</h3>
|
||
<div class="outline-text-3" id="text-2-5">
|
||
<p>
|
||
Effects of crystal’s pose errors on the output beam are summarized in Table <a href="#orgcdf2fe9">3</a>.
|
||
Note that the three pose errors are well decoupled regarding their effects on the output beam.
|
||
Also note that the effect of an error in crystal’s distance does not depend on the Bragg angle.
|
||
</p>
|
||
|
||
<table id="orgcdf2fe9" border="2" cellspacing="0" cellpadding="6" rules="groups" frame="hsides">
|
||
<caption class="t-above"><span class="table-number">Table 3:</span> Summary of the effects of the errors in second crystal’s pose on the output beam</caption>
|
||
|
||
<colgroup>
|
||
<col class="org-left" />
|
||
|
||
<col class="org-right" />
|
||
|
||
<col class="org-left" />
|
||
|
||
<col class="org-left" />
|
||
</colgroup>
|
||
<thead>
|
||
<tr>
|
||
<th scope="col" class="org-left"><b>Beam Motion</b></th>
|
||
<th scope="col" class="org-right"><b>Crystal</b> \(d_z\)</th>
|
||
<th scope="col" class="org-left"><b>Crystal</b> \(r_x\)</th>
|
||
<th scope="col" class="org-left"><b>Crystal</b> \(r_y\)</th>
|
||
</tr>
|
||
</thead>
|
||
<tbody>
|
||
<tr>
|
||
<td class="org-left">\(d_{b,y}\)</td>
|
||
<td class="org-right">0</td>
|
||
<td class="org-left">\(\approx 0\)</td>
|
||
<td class="org-left">0</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-left">\(d_{b,z}\)</td>
|
||
<td class="org-right">\(\boxed{d_{z} 2 \cos(\theta)}\)</td>
|
||
<td class="org-left">0</td>
|
||
<td class="org-left">\(\approx 0\)</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-left">\(r_{b,y}\)</td>
|
||
<td class="org-right">0</td>
|
||
<td class="org-left">0</td>
|
||
<td class="org-left">\(\boxed{2 r_{y}}\)</td>
|
||
</tr>
|
||
|
||
<tr>
|
||
<td class="org-left">\(r_{b,z}\)</td>
|
||
<td class="org-right">0</td>
|
||
<td class="org-left">\(\boxed{- r_x 2 \sin(\theta)}\)</td>
|
||
<td class="org-left">0</td>
|
||
</tr>
|
||
</tbody>
|
||
</table>
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-org86a50e2" class="outline-3">
|
||
<h3 id="org86a50e2"><span class="section-number-3">2.6.</span> “Channel cut” Scan</h3>
|
||
<div class="outline-text-3" id="text-2-6">
|
||
<p>
|
||
A “channel cut” scan is a Bragg scan where the distance between the crystals is fixed.
|
||
</p>
|
||
|
||
<p>
|
||
This is visually shown in Figure <a href="#org02efdb2">8</a> where it is clear that the output beam experiences some vertical motion.
|
||
</p>
|
||
|
||
|
||
<div id="org02efdb2" class="figure">
|
||
<p><img src="figs/ray_tracing_channel_cut.png" alt="ray_tracing_channel_cut.png" />
|
||
</p>
|
||
<p><span class="figure-number">Figure 8: </span>Visual Effect of a channel cut scan</p>
|
||
</div>
|
||
|
||
<p>
|
||
The \(z\) offset of the beam for several channel cut scans are shown in Figure <a href="#org8f490dd">9</a>.
|
||
</p>
|
||
|
||
|
||
<div id="org8f490dd" class="figure">
|
||
<p><img src="figs/channel_cut_scan.png" alt="channel_cut_scan.png" />
|
||
</p>
|
||
<p><span class="figure-number">Figure 9: </span>Z motion of the beam during “channel cut” scans</p>
|
||
</div>
|
||
</div>
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-orgc33f680" class="outline-2">
|
||
<h2 id="orgc33f680"><span class="section-number-2">3.</span> Determining relative pose between the crystals using the X-ray</h2>
|
||
<div class="outline-text-2" id="text-3">
|
||
<p>
|
||
As Interferometers are only measuring <i>relative</i> displacement, it is mandatory to initialize them correctly.
|
||
</p>
|
||
|
||
<p>
|
||
They should be initialize in such a way that:
|
||
</p>
|
||
<ul class="org-ul">
|
||
<li>\(r_x = 0\) and \(r_y = 0\) when the two crystallographic planes are parallel</li>
|
||
<li>measured \(d_z\) is equal to the distance between the crystallographic planes</li>
|
||
</ul>
|
||
|
||
<p>
|
||
In order to do that, an external metrology using the x-ray is used.
|
||
</p>
|
||
</div>
|
||
|
||
<div id="outline-container-org97fbcd2" class="outline-3">
|
||
<h3 id="org97fbcd2"><span class="section-number-3">3.1.</span> Determine the \(y\) parallelism - “Rocking Curve”</h3>
|
||
</div>
|
||
|
||
<div id="outline-container-orga0df616" class="outline-3">
|
||
<h3 id="orga0df616"><span class="section-number-3">3.2.</span> Determine the \(x\) parallelism - Bragg Scan</h3>
|
||
</div>
|
||
|
||
|
||
<div id="outline-container-org01554e4" class="outline-3">
|
||
<h3 id="org01554e4"><span class="section-number-3">3.3.</span> Determine the \(z\) distance - Bragg Scan</h3>
|
||
</div>
|
||
|
||
<div id="outline-container-org9347512" class="outline-3">
|
||
<h3 id="org9347512"><span class="section-number-3">3.4.</span> Use Channel cut scan to determine crystal <code>dry</code> parallelism</h3>
|
||
<div class="outline-text-3" id="text-3-4">
|
||
<p>
|
||
Now, let’s suppose we want to determine the <code>dry</code> angle between the crystals.
|
||
</p>
|
||
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab">dry = 1e<span class="org-builtin">-</span>6; <span class="org-comment-delimiter">% </span><span class="org-comment">[rad]</span>
|
||
</pre>
|
||
</div>
|
||
|
||
<p>
|
||
The error is
|
||
</p>
|
||
\begin{equation}
|
||
\boxed{d_{b,z} = -2 d_z \cos()}
|
||
\end{equation}
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-org7a140bf" class="outline-3">
|
||
<h3 id="org7a140bf"><span class="section-number-3">3.5.</span> Effect of an error on Bragg angle</h3>
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-orgf968de4" class="outline-2">
|
||
<h2 id="orgf968de4"><span class="section-number-2">4.</span> Deformations of the Metrology Frame</h2>
|
||
<div class="outline-text-2" id="text-4">
|
||
<p>
|
||
<a id="org6019b13"></a>
|
||
</p>
|
||
<p>
|
||
The transformation matrices are valid only if the metrology frames are solid bodies.
|
||
</p>
|
||
|
||
<p>
|
||
The metrology frame itself is experiencing some deformations due to the gravity.
|
||
When the bragg axis is scanned, the effect of gravity on the metrology frame is changing and this introduce some measurement errors.
|
||
</p>
|
||
|
||
<p>
|
||
This can be calibrated.
|
||
</p>
|
||
</div>
|
||
<div id="outline-container-org573b7e7" class="outline-3">
|
||
<h3 id="org573b7e7"><span class="section-number-3">4.1.</span> Measurement Setup</h3>
|
||
<div class="outline-text-3" id="text-4-1">
|
||
<p>
|
||
Two beam viewers:
|
||
</p>
|
||
<ul class="org-ul">
|
||
<li>one close to the DCM to measure position of the beam</li>
|
||
<li>one far away to the DCM to measure orientation of the beam</li>
|
||
</ul>
|
||
|
||
<p>
|
||
For each Bragg angle, the Fast Jacks are actuated to that the beam is at the center of the beam viewer.
|
||
Then, then position of the crystals as measured by the interferometers is recorded.
|
||
This position is the wanted position for a given Bragg angle.
|
||
</p>
|
||
|
||
|
||
<div id="org6080aeb" class="figure">
|
||
<p><img src="figs/calibration_setup.png" alt="calibration_setup.png" />
|
||
</p>
|
||
<p><span class="figure-number">Figure 10: </span>Schematic of the setup</p>
|
||
</div>
|
||
|
||
<p>
|
||
Detector:
|
||
</p>
|
||
|
||
<p>
|
||
<a href="https://www.baslerweb.com/en/products/cameras/area-scan-cameras/ace/aca1920-40gc/">https://www.baslerweb.com/en/products/cameras/area-scan-cameras/ace/aca1920-40gc/</a>
|
||
</p>
|
||
|
||
<p>
|
||
Pixel size depends on the magnification used (1x, 6x, 12x).
|
||
</p>
|
||
|
||
<p>
|
||
Pixel size of camera is 5.86 um x 5.86 um.
|
||
With typical magnification of 6x, pixel size is ~1.44um x 1.44um
|
||
</p>
|
||
|
||
<p>
|
||
Frame rate is: 42 fps
|
||
</p>
|
||
</div>
|
||
</div>
|
||
|
||
|
||
<div id="outline-container-orgbabe313" class="outline-3">
|
||
<h3 id="orgbabe313"><span class="section-number-3">4.2.</span> Simulations</h3>
|
||
<div class="outline-text-3" id="text-4-2">
|
||
<p>
|
||
The deformations of the metrology frame and therefore the expected interferometric measurements can be computed as a function of the Bragg angle.
|
||
This may be done using FE software.
|
||
</p>
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-orgaf7afd6" class="outline-3">
|
||
<h3 id="orgaf7afd6"><span class="section-number-3">4.3.</span> Comparison</h3>
|
||
</div>
|
||
|
||
<div id="outline-container-org99a64a4" class="outline-3">
|
||
<h3 id="org99a64a4"><span class="section-number-3">4.4.</span> Test</h3>
|
||
<div class="outline-text-3" id="text-4-4">
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab">aa = importdata(<span class="org-string">"correctInterf-vlm-220201.dat"</span>);
|
||
</pre>
|
||
</div>
|
||
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab"><span class="org-builtin">figure</span>;
|
||
<span class="org-builtin">plot</span>(aa.data(<span class="org-builtin">:</span>,1), aa.data(<span class="org-builtin">:</span>,24))
|
||
</pre>
|
||
</div>
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-org907c706" class="outline-3">
|
||
<h3 id="org907c706"><span class="section-number-3">4.5.</span> Measured frame deformation</h3>
|
||
<div class="outline-text-3" id="text-4-5">
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab">data = table2array(readtable(<span class="org-string">'itf_polynom.csv'</span>,<span class="org-string">'NumHeaderLines'</span>,1));
|
||
th = <span class="org-matlab-math">pi</span><span class="org-builtin">/</span>180<span class="org-builtin">*</span>data(<span class="org-builtin">:</span>,1); <span class="org-comment-delimiter">% </span><span class="org-comment">[rad]</span>
|
||
fj = 0.030427 <span class="org-builtin">-</span> 10.51e<span class="org-builtin">-</span>3<span class="org-builtin">./</span>(2<span class="org-builtin">*</span>cos(th)); <span class="org-comment-delimiter">% </span><span class="org-comment">[m]</span>
|
||
rx2 = 1e<span class="org-builtin">-</span>9<span class="org-builtin">*</span>data(<span class="org-builtin">:</span>,2); <span class="org-comment-delimiter">% </span><span class="org-comment">[rad]</span>
|
||
ry2 = 1e<span class="org-builtin">-</span>9<span class="org-builtin">*</span>data(<span class="org-builtin">:</span>,3); <span class="org-comment-delimiter">% </span><span class="org-comment">[rad]</span>
|
||
rx1 = 1e<span class="org-builtin">-</span>9<span class="org-builtin">*</span>data(<span class="org-builtin">:</span>,4); <span class="org-comment-delimiter">% </span><span class="org-comment">[rad]</span>
|
||
ry1 = 1e<span class="org-builtin">-</span>9<span class="org-builtin">*</span>data(<span class="org-builtin">:</span>,5); <span class="org-comment-delimiter">% </span><span class="org-comment">[rad]</span>
|
||
</pre>
|
||
</div>
|
||
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab"><span class="org-builtin">figure</span>;
|
||
<span class="org-builtin">hold</span> <span class="org-matlab-commanddual-string">on;</span>
|
||
<span class="org-builtin">plot</span>(180<span class="org-builtin">/</span><span class="org-matlab-math">pi</span><span class="org-builtin">*</span>th, 1e6<span class="org-builtin">*</span>detrend(<span class="org-builtin">-</span>rx2, 1), <span class="org-string">'.'</span>)
|
||
<span class="org-comment-delimiter">% </span><span class="org-comment">plot(180/pi*th, detrend(ry2, 1))</span>
|
||
<span class="org-comment-delimiter">% </span><span class="org-comment">plot(180/pi*th, detrend(rx1, 1))</span>
|
||
<span class="org-comment-delimiter">% </span><span class="org-comment">plot(180/pi*th, detrend(ry1, 1))</span>
|
||
<span class="org-builtin">hold</span> <span class="org-matlab-commanddual-string">off;</span>
|
||
<span class="org-builtin">xlabel</span>(<span class="org-string">'Bragg Angle [deg]'</span>);
|
||
<span class="org-builtin">ylabel</span>(<span class="org-string">'Measured $R_x$ [$\mu$rad]'</span>)
|
||
<span class="org-builtin">xlim</span>([10, 75]);
|
||
</pre>
|
||
</div>
|
||
|
||
|
||
<div id="org85f32b7" class="figure">
|
||
<p><img src="figs/calibration_drx_pres.png" alt="calibration_drx_pres.png" />
|
||
</p>
|
||
<p><span class="figure-number">Figure 11: </span>description</p>
|
||
</div>
|
||
|
||
|
||
<p>
|
||
Strange that there is correlation between Rx and Ry.
|
||
</p>
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab"><span class="org-builtin">figure</span>;
|
||
<span class="org-builtin">hold</span> <span class="org-matlab-commanddual-string">on;</span>
|
||
<span class="org-builtin">plot</span>(108<span class="org-builtin">/</span><span class="org-matlab-math">pi</span><span class="org-builtin">*</span>th, 1e9<span class="org-builtin">*</span>detrend(rx1, 0), <span class="org-string">'-'</span>, <span class="org-string">'DisplayName'</span>, <span class="org-string">'$Rx_1$'</span>)
|
||
<span class="org-builtin">plot</span>(108<span class="org-builtin">/</span><span class="org-matlab-math">pi</span><span class="org-builtin">*</span>th, 1e9<span class="org-builtin">*</span>detrend(ry1, 0), <span class="org-string">'-'</span>, <span class="org-string">'DisplayName'</span>, <span class="org-string">'$Ry_1$'</span>)
|
||
<span class="org-builtin">hold</span> <span class="org-matlab-commanddual-string">off;</span>
|
||
<span class="org-builtin">xlabel</span>(<span class="org-string">'Bragg Angle [deg]'</span>); <span class="org-builtin">ylabel</span>(<span class="org-string">'Angle Offset [nrad]'</span>);
|
||
<span class="org-builtin">legend</span>()
|
||
</pre>
|
||
</div>
|
||
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab"><span class="org-builtin">figure</span>;
|
||
<span class="org-builtin">hold</span> <span class="org-matlab-commanddual-string">on;</span>
|
||
<span class="org-builtin">plot</span>(1e3<span class="org-builtin">*</span>fj, detrend(rx2, 1), <span class="org-string">'DisplayName'</span>, <span class="org-string">'$Rx_1$'</span>)
|
||
<span class="org-builtin">plot</span>(1e3<span class="org-builtin">*</span>fj, detrend(ry2, 1), <span class="org-string">'DisplayName'</span>, <span class="org-string">'$Ry_1$'</span>)
|
||
<span class="org-builtin">hold</span> <span class="org-matlab-commanddual-string">off;</span>
|
||
<span class="org-builtin">xlabel</span>(<span class="org-string">'Fast Jack Displacement [mm]'</span>); <span class="org-builtin">ylabel</span>(<span class="org-string">'Angle Offset [nrad]'</span>);
|
||
<span class="org-builtin">legend</span>()
|
||
</pre>
|
||
</div>
|
||
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab"><span class="org-matlab-cellbreak">%% Compute best polynomial fit</span>
|
||
f_rx2 = fit(180<span class="org-builtin">/</span><span class="org-matlab-math">pi</span><span class="org-builtin">*</span>th, 1e9<span class="org-builtin">*</span>rx2, <span class="org-string">'poly4'</span>);
|
||
f_ry2 = fit(180<span class="org-builtin">/</span><span class="org-matlab-math">pi</span><span class="org-builtin">*</span>th, 1e9<span class="org-builtin">*</span>ry2, <span class="org-string">'poly4'</span>);
|
||
f_rx1 = fit(180<span class="org-builtin">/</span><span class="org-matlab-math">pi</span><span class="org-builtin">*</span>th, 1e9<span class="org-builtin">*</span>rx1, <span class="org-string">'poly4'</span>);
|
||
f_ry1 = fit(180<span class="org-builtin">/</span><span class="org-matlab-math">pi</span><span class="org-builtin">*</span>th, 1e9<span class="org-builtin">*</span>ry1, <span class="org-string">'poly4'</span>);
|
||
</pre>
|
||
</div>
|
||
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab"><span class="org-builtin">figure</span>;
|
||
<span class="org-builtin">hold</span> <span class="org-matlab-commanddual-string">on;</span>
|
||
<span class="org-builtin">plot</span>(180<span class="org-builtin">/</span><span class="org-matlab-math">pi</span><span class="org-builtin">*</span>th, f_rx2(180<span class="org-builtin">/</span><span class="org-matlab-math">pi</span><span class="org-builtin">*</span>th))
|
||
<span class="org-builtin">plot</span>(180<span class="org-builtin">/</span><span class="org-matlab-math">pi</span><span class="org-builtin">*</span>th, f_ry2(180<span class="org-builtin">/</span><span class="org-matlab-math">pi</span><span class="org-builtin">*</span>th))
|
||
<span class="org-builtin">plot</span>(180<span class="org-builtin">/</span><span class="org-matlab-math">pi</span><span class="org-builtin">*</span>th, f_rx1(180<span class="org-builtin">/</span><span class="org-matlab-math">pi</span><span class="org-builtin">*</span>th))
|
||
<span class="org-builtin">plot</span>(180<span class="org-builtin">/</span><span class="org-matlab-math">pi</span><span class="org-builtin">*</span>th, f_ry1(180<span class="org-builtin">/</span><span class="org-matlab-math">pi</span><span class="org-builtin">*</span>th))
|
||
<span class="org-builtin">hold</span> <span class="org-matlab-commanddual-string">off;</span>
|
||
</pre>
|
||
</div>
|
||
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab"><span class="org-builtin">figure</span>;
|
||
<span class="org-builtin">hold</span> <span class="org-matlab-commanddual-string">on;</span>
|
||
<span class="org-builtin">plot</span>(th, f_rx2(th) <span class="org-builtin">-</span> rx2)
|
||
<span class="org-builtin">plot</span>(th, f_ry2(th) <span class="org-builtin">-</span> ry2)
|
||
<span class="org-builtin">plot</span>(th, f_rx1(th) <span class="org-builtin">-</span> rx1)
|
||
<span class="org-builtin">plot</span>(th, f_ry1(th) <span class="org-builtin">-</span> ry1)
|
||
<span class="org-builtin">hold</span> <span class="org-matlab-commanddual-string">off;</span>
|
||
</pre>
|
||
</div>
|
||
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab"><span class="org-builtin">figure</span>;
|
||
<span class="org-builtin">hold</span> <span class="org-matlab-commanddual-string">on;</span>
|
||
<span class="org-builtin">plot</span>(th, f(th))
|
||
<span class="org-builtin">plot</span>(th, rx2, <span class="org-string">'.'</span>)
|
||
</pre>
|
||
</div>
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-orgaa71cc1" class="outline-3">
|
||
<h3 id="orgaa71cc1"><span class="section-number-3">4.6.</span> Test</h3>
|
||
<div class="outline-text-3" id="text-4-6">
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab">filename = <span class="org-string">"/home/thomas/mnt/data_id21/22Jan/blc13550/id21/test_xtal1_interf/test_xtal1_interf_0001/test_xtal1_interf_0001.h5"</span>;
|
||
</pre>
|
||
</div>
|
||
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab">data = struct();
|
||
|
||
data.xtal1_111_u = double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA1_SSIM4/data'</span>));
|
||
data.xtal1_111_m = double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA1_SSIM5/data'</span>));
|
||
data.xtal1_111_d = double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA1_SSIM6/data'</span>));
|
||
data.mframe_u = double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA2_SSIM3/data'</span>));
|
||
data.mframe_dh = double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA2_SSIM4/data'</span>));
|
||
data.mframe_dr = double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA2_SSIM5/data'</span>));
|
||
data.bragg = (<span class="org-matlab-math">pi</span><span class="org-builtin">/</span>180)<span class="org-builtin">*</span>double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA2_SSIM6/data'</span>));
|
||
data.fj_pos = 0.030427 <span class="org-builtin">-</span> 10.5e<span class="org-builtin">-</span>3<span class="org-builtin">./</span>(2<span class="org-builtin">*</span>cos(data.bragg));
|
||
data.time = double(h5read(filename, <span class="org-string">'/7.1/instrument/time/data'</span>));
|
||
data.rx = 1e<span class="org-builtin">-</span>9<span class="org-builtin">*</span>double(h5read(filename, <span class="org-string">'/7.1/instrument/xtal1_111_rx/data'</span>));
|
||
data.ry = 1e<span class="org-builtin">-</span>9<span class="org-builtin">*</span>double(h5read(filename, <span class="org-string">'/7.1/instrument/xtal1_111_ry/data'</span>));
|
||
data.z = 1e<span class="org-builtin">-</span>9<span class="org-builtin">*</span>double(h5read(filename, <span class="org-string">'/7.1/instrument/xtal1_111_z/data'</span>));
|
||
data.drx = 1e<span class="org-builtin">-</span>9<span class="org-builtin">*</span>double(h5read(filename, <span class="org-string">'/7.1/instrument/xtal_111_drx_filter/data'</span>));
|
||
data.dry = 1e<span class="org-builtin">-</span>9<span class="org-builtin">*</span>double(h5read(filename, <span class="org-string">'/7.1/instrument/xtal_111_dry_filter/data'</span>));
|
||
data.dz = 1e<span class="org-builtin">-</span>9<span class="org-builtin">*</span>double(h5read(filename, <span class="org-string">'/7.1/instrument/xtal_111_dz_filter/data'</span>));
|
||
data.xtal2_111_u = double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA1_SSIM10/data'</span>))<span class="org-builtin">+</span>10.5e6<span class="org-builtin">./</span>(2<span class="org-builtin">*</span>cos(data.bragg));
|
||
data.xtal2_111_m = double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA2_SSIM1/data'</span>))<span class="org-builtin">+</span>10.5e6<span class="org-builtin">./</span>(2<span class="org-builtin">*</span>cos(data.bragg));
|
||
data.xtal2_111_d = double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA2_SSIM2/data'</span>))<span class="org-builtin">+</span>10.5e6<span class="org-builtin">./</span>(2<span class="org-builtin">*</span>cos(data.bragg));
|
||
</pre>
|
||
</div>
|
||
|
||
|
||
<div id="orgef4ecac" class="figure">
|
||
<p><img src="figs/drifts_xtal2_detrend.png" alt="drifts_xtal2_detrend.png" />
|
||
</p>
|
||
<p><span class="figure-number">Figure 12: </span>Drifts of the second crystal as a function of Bragg Angle</p>
|
||
</div>
|
||
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab"><span class="org-builtin">figure</span>;
|
||
<span class="org-builtin">hold</span> <span class="org-matlab-commanddual-string">on;</span>
|
||
<span class="org-builtin">plot</span>(108<span class="org-builtin">/</span><span class="org-matlab-math">pi</span><span class="org-builtin">*</span>data.bragg, data.drx)
|
||
<span class="org-builtin">plot</span>(108<span class="org-builtin">/</span><span class="org-matlab-math">pi</span><span class="org-builtin">*</span>data.bragg, data.dry)
|
||
<span class="org-builtin">hold</span> <span class="org-matlab-commanddual-string">off;</span>
|
||
</pre>
|
||
</div>
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-org73400be" class="outline-3">
|
||
<h3 id="org73400be"><span class="section-number-3">4.7.</span> Repeatability of frame deformation</h3>
|
||
<div class="outline-text-3" id="text-4-7">
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab">filename = <span class="org-string">"/home/thomas/mnt/data_id21/22Jan/blc13550/id21/test_xtal1_interf/test_xtal1_interf_0001/test_xtal1_interf_0001.h5"</span>;
|
||
</pre>
|
||
</div>
|
||
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab">data_1 = struct();
|
||
data_1.time = double(h5read(filename, <span class="org-string">'/7.1/instrument/time/data'</span>));
|
||
data_1.bragg = (<span class="org-matlab-math">pi</span><span class="org-builtin">/</span>180)<span class="org-builtin">*</span>double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA2_SSIM6/data'</span>));
|
||
data_1.fj_pos = 0.030427 <span class="org-builtin">-</span> 10.5e<span class="org-builtin">-</span>3<span class="org-builtin">./</span>(2<span class="org-builtin">*</span>cos(data_1.bragg));
|
||
data_1.xtal1_111_u = double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA1_SSIM4/data'</span>));
|
||
data_1.xtal1_111_m = double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA1_SSIM5/data'</span>));
|
||
data_1.xtal1_111_d = double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA1_SSIM6/data'</span>));
|
||
data_1.xtal2_111_u = double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA1_SSIM10/data'</span>))<span class="org-builtin">+</span>10.5e6<span class="org-builtin">./</span>(2<span class="org-builtin">*</span>cos(data_1.bragg));
|
||
data_1.xtal2_111_m = double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA2_SSIM1/data'</span>))<span class="org-builtin">+</span>10.5e6<span class="org-builtin">./</span>(2<span class="org-builtin">*</span>cos(data_1.bragg));
|
||
data_1.xtal2_111_d = double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA2_SSIM2/data'</span>))<span class="org-builtin">+</span>10.5e6<span class="org-builtin">./</span>(2<span class="org-builtin">*</span>cos(data_1.bragg));
|
||
data_1.mframe_u = double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA2_SSIM3/data'</span>));
|
||
data_1.mframe_dh = double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA2_SSIM4/data'</span>));
|
||
data_1.mframe_dr = double(h5read(filename, <span class="org-string">'/7.1/instrument/FPGA2_SSIM5/data'</span>));
|
||
data_1.drx = 1e<span class="org-builtin">-</span>9<span class="org-builtin">*</span>double(h5read(filename, <span class="org-string">'/7.1/instrument/xtal_111_drx_filter/data'</span>));
|
||
data_1.dry = 1e<span class="org-builtin">-</span>9<span class="org-builtin">*</span>double(h5read(filename, <span class="org-string">'/7.1/instrument/xtal_111_dry_filter/data'</span>));
|
||
data_1.dz = 1e<span class="org-builtin">-</span>9<span class="org-builtin">*</span>double(h5read(filename, <span class="org-string">'/7.1/instrument/xtal_111_dz_filter/data'</span>));
|
||
</pre>
|
||
</div>
|
||
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab">data_2 = struct();
|
||
data_2.time = double(h5read(filename, <span class="org-string">'/6.1/instrument/time/data'</span>));
|
||
data_2.bragg = (<span class="org-matlab-math">pi</span><span class="org-builtin">/</span>180)<span class="org-builtin">*</span>double(h5read(filename, <span class="org-string">'/6.1/instrument/FPGA2_SSIM6/data'</span>));
|
||
data_2.fj_pos = 0.030427 <span class="org-builtin">-</span> 10.5e<span class="org-builtin">-</span>3<span class="org-builtin">./</span>(2<span class="org-builtin">*</span>cos(data_2.bragg));
|
||
data_2.xtal1_111_u = double(h5read(filename, <span class="org-string">'/6.1/instrument/FPGA1_SSIM4/data'</span>));
|
||
data_2.xtal1_111_m = double(h5read(filename, <span class="org-string">'/6.1/instrument/FPGA1_SSIM5/data'</span>));
|
||
data_2.xtal1_111_d = double(h5read(filename, <span class="org-string">'/6.1/instrument/FPGA1_SSIM6/data'</span>));
|
||
data_2.xtal2_111_u = double(h5read(filename, <span class="org-string">'/6.1/instrument/FPGA1_SSIM10/data'</span>))<span class="org-builtin">+</span>10.5e6<span class="org-builtin">./</span>(2<span class="org-builtin">*</span>cos(data_2.bragg));
|
||
data_2.xtal2_111_m = double(h5read(filename, <span class="org-string">'/6.1/instrument/FPGA2_SSIM1/data'</span>))<span class="org-builtin">+</span>10.5e6<span class="org-builtin">./</span>(2<span class="org-builtin">*</span>cos(data_2.bragg));
|
||
data_2.xtal2_111_d = double(h5read(filename, <span class="org-string">'/6.1/instrument/FPGA2_SSIM2/data'</span>))<span class="org-builtin">+</span>10.5e6<span class="org-builtin">./</span>(2<span class="org-builtin">*</span>cos(data_2.bragg));
|
||
data_2.mframe_u = double(h5read(filename, <span class="org-string">'/6.1/instrument/FPGA2_SSIM3/data'</span>));
|
||
data_2.mframe_dh = double(h5read(filename, <span class="org-string">'/6.1/instrument/FPGA2_SSIM4/data'</span>));
|
||
data_2.mframe_dr = double(h5read(filename, <span class="org-string">'/6.1/instrument/FPGA2_SSIM5/data'</span>));
|
||
data_2.drx = 1e<span class="org-builtin">-</span>9<span class="org-builtin">*</span>double(h5read(filename, <span class="org-string">'/6.1/instrument/xtal_111_drx_filter/data'</span>));
|
||
data_2.dry = 1e<span class="org-builtin">-</span>9<span class="org-builtin">*</span>double(h5read(filename, <span class="org-string">'/6.1/instrument/xtal_111_dry_filter/data'</span>));
|
||
data_2.dz = 1e<span class="org-builtin">-</span>9<span class="org-builtin">*</span>double(h5read(filename, <span class="org-string">'/6.1/instrument/xtal_111_dz_filter/data'</span>));
|
||
</pre>
|
||
</div>
|
||
</div>
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-org0b40d1b" class="outline-2">
|
||
<h2 id="org0b40d1b"><span class="section-number-2">5.</span> Attocube - Periodic Non-Linearity</h2>
|
||
<div class="outline-text-2" id="text-5">
|
||
<p>
|
||
<a id="org2a9ae01"></a>
|
||
</p>
|
||
<p>
|
||
Interferometers have some periodic nonlinearity (NO_ITEM_DATA:thurner15_fiber_based_distan_sensin_inter).
|
||
The period is a fraction of the wavelength (usually \(\lambda/2\)) and can be due to polarization mixing, non perfect alignment of the optical components and unwanted reflected beams (<a href="#citeproc_bib_item_1">Ducourtieux 2018</a>; <a href="#citeproc_bib_item_2">Thurner et al. 2015</a>).
|
||
The amplitude of the nonlinearity can vary from a fraction of a nanometer to tens of nanometers.
|
||
</p>
|
||
|
||
<p>
|
||
In the DCM case, when using Attocube interferometers, the period non-linearity are in the order of several nanometers with a period of \(765\,nm\).
|
||
This is inducing some positioning errors which are too high.
|
||
</p>
|
||
|
||
<p>
|
||
In order to overcome this issue, the periodic non-linearity of the interferometers have to be calibrated.
|
||
To do so, a displacement is imposed and measured both by the interferometers and by another metrology system which does not have this nonlinearity.
|
||
By comparing the two measured displacements, the nonlinearity can be calibration.
|
||
This process is performed over several periods in order to characterize the error over the full stroke.
|
||
</p>
|
||
</div>
|
||
|
||
<div id="outline-container-orge0a38ec" class="outline-3">
|
||
<h3 id="orge0a38ec"><span class="section-number-3">5.1.</span> Measurement Setup</h3>
|
||
<div class="outline-text-3" id="text-5-1">
|
||
<p>
|
||
The metrology that will be compared with the interferometers are the strain gauges incorporated in the PI piezoelectric stacks.
|
||
</p>
|
||
|
||
<p>
|
||
It is here supposed that the measured displacement by the strain gauges are converted to the displacement at the interferometer locations.
|
||
It is also supposed that we are at a certain Bragg angle, and that the stepper motors are not moving: only the piezoelectric actuators are used.
|
||
</p>
|
||
|
||
<div class="note" id="orgfc7f1df">
|
||
<p>
|
||
Note that the strain gauges are measuring the relative displacement of the piezoelectric stacks while the interferometers are measuring the relative motion between the second crystals and the metrology frame.
|
||
</p>
|
||
|
||
<p>
|
||
Only the interferometers measuring the second crystal motion can be calibrated here.
|
||
</p>
|
||
|
||
<p>
|
||
As any deformations of the metrology frame of deformation of the crystal’s support can degrade the quality of the calibration, it is better to perform this calibration without any bragg angle motion.
|
||
</p>
|
||
|
||
</div>
|
||
|
||
<p>
|
||
The setup is schematically with the block diagram in Figure <a href="#org5485724">13</a>.
|
||
</p>
|
||
|
||
<p>
|
||
The signals are:
|
||
</p>
|
||
<ul class="org-ul">
|
||
<li>\(u\): Reference Signal sent to the PI controller (position where we wish to three stacks to be).
|
||
The PI controller takes care or controlling to position as measured by the strain gauges such that it is close to the reference position.</li>
|
||
<li>\(d\): Disturbances affecting the position of the crystals</li>
|
||
<li>\(y\): Displacement of the crystal as measured by one interferometer</li>
|
||
<li>\(y_g\): Measurement of the motion in the frame of the interferometer by the strain gauge with some noise \(n_g\)</li>
|
||
<li>\(y_a\): Measurement of the crystal motion by the interferometer with some noise \(n_a\)</li>
|
||
</ul>
|
||
|
||
|
||
<div id="org5485724" class="figure">
|
||
<p><img src="figs/block_diagram_lut_attocube.png" alt="block_diagram_lut_attocube.png" />
|
||
</p>
|
||
<p><span class="figure-number">Figure 13: </span>Block Diagram schematic of the setup used to measure the periodic non-linearity of the Attocube</p>
|
||
</div>
|
||
|
||
<p>
|
||
The problem is to estimate the periodic non-linearity of the Attocube from the imperfect measurements \(y_a\) and \(y_g\).
|
||
</p>
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-org821469c" class="outline-3">
|
||
<h3 id="org821469c"><span class="section-number-3">5.2.</span> Choice of the reference signal</h3>
|
||
<div class="outline-text-3" id="text-5-2">
|
||
<p>
|
||
The main specifications for the reference signal are;
|
||
</p>
|
||
<ul class="org-ul">
|
||
<li>sweep several periods (i.e. several micrometers)</li>
|
||
<li>stay in the linear region of the strain gauge</li>
|
||
<li>no excitation of mechanical modes (i.e. the frequency content of the signal should be at low frequency)</li>
|
||
<li>no phase shift due to limited bandwidth of both the interferometers and the strain gauge</li>
|
||
<li>the full process should be quite fast</li>
|
||
</ul>
|
||
|
||
<p>
|
||
The travel range of the piezoelectric stacks is 15 micrometers, the resolution of the strain gauges is 0.3nm and the maximum non-linearity is 0.15%.
|
||
If one non-linear period is swept (765nm), the maximum estimation error of the strain gauge is around 1nm.
|
||
</p>
|
||
|
||
<p>
|
||
Based on the above discussion, one suitable excitation signal is a sinusoidal sweep with a frequency of 10Hz.
|
||
</p>
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-org3263c1f" class="outline-3">
|
||
<h3 id="org3263c1f"><span class="section-number-3">5.3.</span> Repeatability of the non-linearity</h3>
|
||
<div class="outline-text-3" id="text-5-3">
|
||
<p>
|
||
Instead of calibrating the non-linear errors of the interferometers over the full fast jack stroke (25mm), one can only calibrate the errors of one period.
|
||
</p>
|
||
|
||
<p>
|
||
For that, we need to make sure that the errors are repeatable from one period to the other and also the period should be very precisely estimated (i.e. the wavelength of the laser).
|
||
</p>
|
||
|
||
<p>
|
||
Also, the laser wavelength should be very stable (specified at 50ppb).
|
||
</p>
|
||
|
||
<p>
|
||
One way to precisely estimate the laser wavelength is to estimate the non linear errors of the interferometer at an initial position, and then to estimate the non linear errors at a large offset, say 10mm.
|
||
</p>
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-org39fabc8" class="outline-3">
|
||
<h3 id="org39fabc8"><span class="section-number-3">5.4.</span> Simulation</h3>
|
||
<div class="outline-text-3" id="text-5-4">
|
||
<p>
|
||
Suppose we have a first approximation of the non-linear period.
|
||
</p>
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab">period_est = 765e<span class="org-builtin">-</span>9; <span class="org-comment-delimiter">% </span><span class="org-comment">Estimated period [m]</span>
|
||
</pre>
|
||
</div>
|
||
|
||
<p>
|
||
And suppose the real period of the non-linear errors is a little bit above (by 0.02nm):
|
||
</p>
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab">period_err = 0.02e<span class="org-builtin">-</span>9; <span class="org-comment-delimiter">% </span><span class="org-comment">Error on the period estimation [m]</span>
|
||
period_nl = period_est <span class="org-builtin">+</span> period_err; <span class="org-comment-delimiter">% </span><span class="org-comment">Period of the non-linear errors [m]</span>
|
||
</pre>
|
||
</div>
|
||
|
||
<p>
|
||
The non-linear errors are first estimated at the beginning of the stroke (Figure <a href="#org49c2d00">14</a>).
|
||
</p>
|
||
|
||
<div id="org49c2d00" class="figure">
|
||
<p><img src="figs/non_linear_errors_start_stroke.png" alt="non_linear_errors_start_stroke.png" />
|
||
</p>
|
||
<p><span class="figure-number">Figure 14: </span>Estimation of the non-linear errors at the beginning of the stroke</p>
|
||
</div>
|
||
|
||
<p>
|
||
From this only measurement, it is not possible to estimate with great accuracy the period of the error.
|
||
To do so, the same measurement is performed with a stroke of several millimeters (Figure <a href="#org8c5db79">15</a>).
|
||
</p>
|
||
|
||
<p>
|
||
It can be seen that there is an offset between the estimated and the measured errors.
|
||
This is due to a mismatch between the estimated period and the true period of the error.
|
||
</p>
|
||
|
||
|
||
<div id="org8c5db79" class="figure">
|
||
<p><img src="figs/non_linear_errors_middle_stroke.png" alt="non_linear_errors_middle_stroke.png" />
|
||
</p>
|
||
<p><span class="figure-number">Figure 15: </span>Estimated non-linear errors at a latter position</p>
|
||
</div>
|
||
|
||
<p>
|
||
Suppose the non-linear error is characterized by a periodic function \(\mathcal{E}\), to simplify let’s take a sinusoidal function (this can be generalized by taking the fourier transform of the function):
|
||
</p>
|
||
\begin{equation}
|
||
\mathcal{E}(x) = \sin\left(\frac{x}{\lambda}\right)
|
||
\end{equation}
|
||
<p>
|
||
with \(x\) the displacement and \(\lambda\) the period of the error.
|
||
</p>
|
||
|
||
<p>
|
||
The measured error at \(x_0\) is then:
|
||
</p>
|
||
\begin{equation}
|
||
\mathcal{E}_m(x_0) = \sin\left( \frac{x_0}{\lambda} \right)
|
||
\end{equation}
|
||
<p>
|
||
And the estimated one is:
|
||
</p>
|
||
\begin{equation}
|
||
\mathcal{E}_e(x_0) = \sin \left( \frac{x_0}{\lambda_{\text{est}}} \right)
|
||
\end{equation}
|
||
<p>
|
||
with \(\lambda_{\text{est}}\) the estimated error’s period.
|
||
</p>
|
||
|
||
<p>
|
||
From Figure <a href="#org8c5db79">15</a>, we can see that there is an offset between the two curves.
|
||
Let’s call this offset \(\epsilon_x\), we then have:
|
||
</p>
|
||
\begin{equation}
|
||
\mathcal{E}_m(x_0) = \mathcal{E}_e(x_0 + \epsilon_x)
|
||
\end{equation}
|
||
|
||
<p>
|
||
Which gives us:
|
||
</p>
|
||
\begin{equation}
|
||
\sin\left( \frac{x_0}{\lambda} \right) = \sin \left( \frac{x_0 + \epsilon_x}{\lambda_{\text{est}}} \right)
|
||
\end{equation}
|
||
|
||
<p>
|
||
Finally:
|
||
</p>
|
||
\begin{equation}
|
||
\boxed{\lambda = \lambda_{\text{est}} \frac{x_0}{x_0 + \epsilon_x}}
|
||
\end{equation}
|
||
|
||
<p>
|
||
The estimated delay is computed:
|
||
</p>
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab"><span class="org-matlab-cellbreak">%% Estimation of the offset between the estimated and measured errors</span>
|
||
i_period = stroke <span class="org-builtin">></span> 5e<span class="org-builtin">-</span>3<span class="org-builtin">-</span>period_nl<span class="org-builtin">/</span>2 <span class="org-builtin">&</span> stroke <span class="org-builtin"><</span> 5e<span class="org-builtin">-</span>3<span class="org-builtin">+</span>period_nl<span class="org-builtin">/</span>2;
|
||
epsilon_x = finddelay(nl_errors(i_period), est_errors(i_period)) <span class="org-comment-delimiter">% </span><span class="org-comment">[m]</span>
|
||
</pre>
|
||
</div>
|
||
|
||
<pre class="example">
|
||
Estimated delay x0 is -120 [nm]
|
||
</pre>
|
||
|
||
|
||
<p>
|
||
And the period \(\lambda\) can be estimated:
|
||
</p>
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab"><span class="org-matlab-cellbreak">%% Computation of the period [m]</span>
|
||
period_fin = period_est <span class="org-builtin">*</span> (5e<span class="org-builtin">-</span>3)<span class="org-builtin">/</span>(5e<span class="org-builtin">-</span>3 <span class="org-builtin">+</span> d_offset); <span class="org-comment-delimiter">% </span><span class="org-comment">Estimated period after measurement [m]</span>
|
||
</pre>
|
||
</div>
|
||
|
||
<pre class="example">
|
||
The estimated period is 765.020 [nm]
|
||
</pre>
|
||
|
||
|
||
<p>
|
||
And the results confirms that this method is working on paper.
|
||
</p>
|
||
|
||
<p>
|
||
When doing this computation, we suppose that there are <b>at most</b> one half period of offset between the estimated and the measured non-linear (to not have any ambiguity whether the estimated period is too large or too small).
|
||
Mathematically this means that the displacement \(x_0\) should be smaller than:
|
||
</p>
|
||
\begin{equation}
|
||
x_0 < \frac{1}{2} \cdot \lambda \cdot \frac{\lambda}{\epsilon_\lambda}
|
||
\end{equation}
|
||
<p>
|
||
With \(\epsilon_\lambda\) the absolute estimation error of the period in meters.
|
||
</p>
|
||
|
||
<p>
|
||
For instance, if we estimate the error on the period to be less than 0.1nm, the maximum displacement is:
|
||
</p>
|
||
<div class="org-src-container">
|
||
<pre class="src src-matlab"><span class="org-matlab-cellbreak">%% Estimated maximum stroke [m]</span>
|
||
max_x0 = 0.5 <span class="org-builtin">*</span> 765e<span class="org-builtin">-</span>9 <span class="org-builtin">*</span> (765e<span class="org-builtin">-</span>9)<span class="org-builtin">/</span>(0.1e<span class="org-builtin">-</span>9);
|
||
</pre>
|
||
</div>
|
||
|
||
<pre class="example">
|
||
The maximum stroke is 2.9 [mm]
|
||
</pre>
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-org373814c" class="outline-3">
|
||
<h3 id="org373814c"><span class="section-number-3">5.5.</span> Measurements</h3>
|
||
<div class="outline-text-3" id="text-5-5">
|
||
<p>
|
||
We have some constrains on the way the motion is imposed and measured:
|
||
</p>
|
||
<ul class="org-ul">
|
||
<li>We want the frequency content of the imposed motion to be at low frequency in order not to induce vibrations of the structure.
|
||
We have to make sure the forces applied by the piezoelectric actuator only moves the crystal and not the fast jack below.
|
||
Therefore, we have to move much slower than the first resonance frequency in the system.</li>
|
||
<li>As both \(y_a\) and \(y_g\) should have rather small noise, we have to filter them with low pass filters.
|
||
The cut-off frequency of the low pass filter should be high as compared to the motion (to not induce any distortion) but still reducing sufficiently the noise.
|
||
Let’s say we want the noise to be less than 1nm (\(6 \sigma\)).</li>
|
||
</ul>
|
||
|
||
<p>
|
||
Suppose we have the power spectral density (PSD) of both \(n_a\) and \(n_g\).
|
||
</p>
|
||
|
||
<ul class="org-ul">
|
||
<li class="off"><code>[ ]</code> Take the PSD of the Attocube</li>
|
||
<li class="off"><code>[ ]</code> Take the PSD of the strain gauge</li>
|
||
<li class="off"><code>[ ]</code> Using 2nd order low pass filter, estimate the required low pass filter cut-off frequency to have sufficiently low noise</li>
|
||
</ul>
|
||
</div>
|
||
</div>
|
||
</div>
|
||
|
||
<div id="outline-container-orga081f50" class="outline-2">
|
||
<h2 id="orga081f50">Bibliography</h2>
|
||
<div class="outline-text-2" id="text-orga081f50">
|
||
<style>.csl-entry{text-indent: -1.5em; margin-left: 1.5em;}</style><div class="csl-bib-body">
|
||
<div class="csl-entry"><a id="citeproc_bib_item_1"></a>Ducourtieux, Sebastien. 2018. “Toward High Precision Position Control Using Laser Interferometry: Main Sources of Error.” doi:<a href="https://doi.org/10.13140/rg.2.2.21044.35205">10.13140/rg.2.2.21044.35205</a>.</div>
|
||
<div class="csl-entry"><a id="citeproc_bib_item_2"></a>Thurner, Klaus, Francesca Paola Quacquarelli, Pierre-François Braun, Claudio Dal Savio, and Khaled Karrai. 2015. “Fiber-Based Distance Sensing Interferometry.” <i>Applied Optics</i> 54 (10). Optical Society of America: 3051–63.</div>
|
||
</div>
|
||
</div>
|
||
</div>
|
||
</div>
|
||
<div id="postamble" class="status">
|
||
<p class="author">Author: Dehaeze Thomas</p>
|
||
<p class="date">Created: 2022-06-22 Wed 13:29</p>
|
||
</div>
|
||
</body>
|
||
</html>
|