{ "cells": [ { "cell_type": "markdown", "id": "coordinate-balance", "metadata": {}, "source": [ "# Particle motion in electromagnetic fields\n", "\n", "[A. Petrenko](https://www.inp.nsk.su/~petrenko/) (Novosibirsk, 2021)" ] }, { "cell_type": "markdown", "id": "advisory-relations", "metadata": {}, "source": [ "The equations of particle motion under the Lorentz force (in SI units)\n", "\n", "\\begin{equation}\n", "\\frac{d \\mathbf{p}}{dt} = q \\mathbf{E} + q[\\mathbf{v} \\times \\mathbf{B}],\n", "\\end{equation}\n", "\n", "\\begin{equation}\n", "\\frac{d \\mathbf{r}}{dt} = \\mathbf{v},\n", "\\end{equation}\n", "\n", "\\begin{equation}\n", "\\mathbf{p} = \\gamma m \\mathbf{v}.\n", "\\end{equation}\n", "\n", "Particle velocity $\\mathbf{v}$ can be derived from the momentum $\\mathbf{p}$ using the definition of relativistic factor $\\gamma = \\left(1 - \\frac{v^2}{c^2}\\right)^{-1/2} = \\sqrt{1 + \\left(\\frac{p}{mc}\\right)^2}$." ] }, { "cell_type": "markdown", "id": "english-algeria", "metadata": {}, "source": [ "Therefore the equations of motion\n", "\n", "\\begin{equation}\n", "\\frac{ d\\mathbf{r} }{dt} = \\frac{\\mathbf{p}}{\\gamma m},\\\\\n", "\\frac{d\\mathbf{p}}{dt} = q \\mathbf{E} + q\\left[\\frac{\\mathbf{p}}{\\gamma m} \\times \\mathbf{B}\\right].\n", "\\end{equation}\n", "\n", "And the final 6d vector equation to solve $d\\mathbf{Y}/dt = \\mathbf{F}(t,\\mathbf{Y})$:" ] }, { "cell_type": "markdown", "id": "finnish-coverage", "metadata": {}, "source": [ "\\begin{equation}\n", "\\frac{d}{dt}\n", "\\left(\n", "\\begin{array}{c}\n", "x \\\\\n", "y \\\\\n", "z \\\\\n", "p_x \\\\\n", "p_y \\\\\n", "p_z\n", "\\end{array}\n", "\\right) =\n", "\\left(\n", "\\begin{array}{c}\n", "p_x/\\gamma m \\\\\n", "p_y/\\gamma m \\\\\n", "p_z/\\gamma m \\\\\n", "qE_x + q(p_y B_z - p_z B_y)/\\gamma m \\\\\n", "qE_y - q(p_x B_z - p_z B_x)/\\gamma m \\\\\n", "qE_z + q(p_x B_y - p_y B_x)/\\gamma m\n", "\\end{array}\n", "\\right).\n", "\\end{equation}" ] }, { "cell_type": "code", "execution_count": 1, "id": "constitutional-medium", "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "from scipy.integrate import solve_ivp" ] }, { "cell_type": "markdown", "id": "happy-overall", "metadata": {}, "source": [ "Define external fields:" ] }, { "cell_type": "code", "execution_count": 2, "id": "premium-swedish", "metadata": {}, "outputs": [], "source": [ "def E(x,y,z,t): # V/m\n", " Ex = 0; Ey = 0; Ez = 0\n", " return Ex, Ey, Ez" ] }, { "cell_type": "code", "execution_count": 3, "id": "genuine-tanzania", "metadata": {}, "outputs": [], "source": [ "def B(x,y,z,t): # T\n", " Bx = 0; By = 0; Bz = 1.0 # T\n", " return Bx, By, Bz" ] }, { "cell_type": "markdown", "id": "three-politics", "metadata": {}, "source": [ "Define some constants:" ] }, { "cell_type": "code", "execution_count": 4, "id": "junior-minneapolis", "metadata": {}, "outputs": [], "source": [ "c = 299792458 # m/sec\n", "e = 1.602176634e-19 # C\n", "eV = e # J\n", "keV = 1e3*eV; MeV = 1e6*eV; GeV = 1e9*eV\n", "\n", "sec = 1; ns = 1e-9*sec\n", "\n", "# electron:\n", "q = -e\n", "mc2 = 511*keV # J\n", "\n", "mc = mc2/c # kg*m/sec\n", "m = mc/c # kg" ] }, { "cell_type": "code", "execution_count": 5, "id": "indian-output", "metadata": {}, "outputs": [], "source": [ "def F(t,Y):\n", " x = Y[0]; y = Y[1]; z = Y[2]\n", " px = Y[3]; py = Y[4]; pz = Y[5]\n", " gamma = np.sqrt(1 + (px*px + py*py + pz*pz)/(mc*mc))\n", " Ex, Ey, Ez = E(x,y,z,t) # V/m\n", " Bx, By, Bz = B(x,y,z,t) # T\n", " \n", " return np.array([\n", " px/(gamma*m),\n", " py/(gamma*m),\n", " pz/(gamma*m),\n", " q*Ex + q*(py*Bz-pz*By)/(gamma*m),\n", " q*Ey - q*(px*Bz-pz*Bx)/(gamma*m),\n", " q*Ez + q*(px*By-py*Bx)/(gamma*m)\n", " ])" ] }, { "cell_type": "markdown", "id": "interpreted-crash", "metadata": {}, "source": [ "Set the initial particle coordinates" ] }, { "cell_type": "code", "execution_count": 6, "id": "maritime-right", "metadata": {}, "outputs": [], "source": [ "x = 0; y = 0; z = 0\n", "px = 0.5*GeV/c; py = 0; pz = 0.1*GeV/c\n", "\n", "Y0 = np.array([x,y,z,px,py,pz])" ] }, { "cell_type": "code", "execution_count": 7, "id": "bibliographic-radiation", "metadata": {}, "outputs": [], "source": [ "t_span = (0, 200*ns) # sec\n", "\n", "t_eval = np.linspace(t_span[0], t_span[1], 500) # points of output" ] }, { "cell_type": "code", "execution_count": 8, "id": "combined-badge", "metadata": {}, "outputs": [], "source": [ "#help(solve_ivp)" ] }, { "cell_type": "code", "execution_count": 9, "id": "spread-renaissance", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "CPU times: user 20.7 ms, sys: 2.45 ms, total: 23.2 ms\n", "Wall time: 24.1 ms\n" ] } ], "source": [ "%%time\n", "sol = solve_ivp(F, t_span, Y0, t_eval=t_eval, max_step=5*ns) # rtol=1e-4" ] }, { "cell_type": "markdown", "id": "detected-juvenile", "metadata": {}, "source": [ "## Plot the results" ] }, { "cell_type": "code", "execution_count": 10, "id": "separate-sector", "metadata": {}, "outputs": [], "source": [ "import plotly.io as pio\n", "import plotly.graph_objects as go" ] }, { "cell_type": "code", "execution_count": 11, "id": "19628171-50a8-4593-9e73-117b308e5787", "metadata": {}, "outputs": [], "source": [ "#pio.renderers" ] }, { "cell_type": "code", "execution_count": 12, "id": "e304e5d7-07f7-4509-94ef-485f0b713cd4", "metadata": {}, "outputs": [], "source": [ "pio.renderers.default = \"plotly_mimetype+notebook\" # both notebook & html from nbconvert are interactive\n", "#pio.renderers.default = \"svg\" # svg figures (for publications)" ] }, { "cell_type": "code", "execution_count": 13, "id": "5713334f-285e-45d4-8e5a-d778aace4674", "metadata": {}, "outputs": [], "source": [ "pio.templates.default = pio.templates[\"simple_white\"]" ] }, { "cell_type": "code", "execution_count": 14, "id": "protecting-arnold", "metadata": {}, "outputs": [], "source": [ "x, y, z, px, py, pz = sol.y" ] }, { "cell_type": "code", "execution_count": 15, "id": "blind-louisiana", "metadata": {}, "outputs": [ { "data": { "text/html": [ " \n", " " ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/vnd.plotly.v1+json": { "config": { "plotlyServerURL": "https://plot.ly" }, "data": [ { "line": { "color": "blue", "width": 3 }, "mode": "lines", "type": "scatter", "x": [ 0, 0.40080160320641284, 0.8016032064128257, 1.2024048096192386, 1.6032064128256514, 2.0040080160320644, 2.404809619238477, 2.80561122244489, 3.2064128256513027, 3.6072144288577155, 4.008016032064129, 4.408817635270542, 4.809619238476954, 5.2104208416833675, 5.61122244488978, 6.012024048096193, 6.412825651302605, 6.813627254509019, 7.214428857715431, 7.615230460921845, 8.016032064128257, 8.41683366733467, 8.817635270541084, 9.218436873747496, 9.619238476953909, 10.02004008016032, 10.420841683366735, 10.821643286573147, 11.22244488977956, 11.623246492985972, 12.024048096192386, 12.424849699398798, 12.82565130260521, 13.226452905811623, 13.627254509018037, 14.02805611222445, 14.428857715430862, 14.829659318637274, 15.23046092184369, 15.6312625250501, 16.032064128256515, 16.432865731462925, 16.83366733466934, 17.23446893787575, 17.635270541082168, 18.03607214428858, 18.436873747494992, 18.837675350701403, 19.238476953907817, 19.639278557114228, 20.04008016032064, 20.440881763527052, 20.84168336673347, 21.24248496993988, 21.643286573146295, 22.044088176352705, 22.44488977955912, 22.84569138276553, 23.246492985971944, 23.647294589178358, 24.048096192384772, 24.448897795591183, 24.849699398797597, 25.250501002004007, 25.65130260521042, 26.052104208416836, 26.452905811623246, 26.85370741482966, 27.254509018036075, 27.65531062124249, 28.0561122244489, 28.456913827655313, 28.857715430861724, 29.258517034068138, 29.65931863727455, 30.060120240480963, 30.46092184368738, 30.86172344689379, 31.2625250501002, 31.663326653306612, 32.06412825651303, 32.46492985971944, 32.86573146292585, 33.26653306613226, 33.66733466933868, 34.06813627254509, 34.4689378757515, 34.86973947895792, 35.270541082164335, 35.671342685370746, 36.07214428857716, 36.47294589178357, 36.873747494989985, 37.274549098196395, 37.675350701402806, 38.07615230460922, 38.476953907815634, 38.877755511022045, 39.278557114228455, 39.67935871743487, 40.08016032064128, 40.480961923847694, 40.881763527054105, 41.28256513026053, 41.68336673346694, 42.08416833667335, 42.48496993987976, 42.88577154308618, 43.28657314629259, 43.687374749499, 44.08817635270541, 44.48897795591183, 44.88977955911824, 45.29058116232465, 45.69138276553106, 46.09218436873748, 46.49298597194389, 46.8937875751503, 47.294589178356716, 47.695390781563134, 48.096192384769545, 48.496993987975955, 48.897795591182366, 49.29859719438878, 49.699398797595194, 50.100200400801604, 50.501002004008015, 50.90180360721443, 51.30260521042084, 51.703406813627254, 52.10420841683367, 52.50501002004008, 52.90581162324649, 53.3066132264529, 53.70741482965932, 54.10821643286574, 54.50901803607215, 54.90981963927856, 55.31062124248498, 55.71142284569139, 56.1122244488978, 56.51302605210421, 56.91382765531063, 57.31462925851704, 57.71543086172345, 58.11623246492986, 58.517034068136276, 58.91783567134269, 59.3186372745491, 59.71943887775551, 60.120240480961925, 60.521042084168336, 60.92184368737476, 61.32264529058117, 61.72344689378758, 62.12424849699399, 62.5250501002004, 62.92585170340681, 63.326653306613224, 63.727454909819635, 64.12825651302606, 64.52905811623248, 64.92985971943888, 65.3306613226453, 65.7314629258517, 66.13226452905812, 66.53306613226452, 66.93386773547095, 67.33466933867736, 67.73547094188378, 68.13627254509018, 68.5370741482966, 68.937875751503, 69.33867735470942, 69.73947895791584, 70.14028056112225, 70.54108216432867, 70.94188376753507, 71.34268537074149, 71.7434869739479, 72.14428857715431, 72.54509018036072, 72.94589178356713, 73.34669338677355, 73.74749498997997, 74.14829659318637, 74.54909819639279, 74.9498997995992, 75.35070140280561, 75.75150300601203, 76.15230460921843, 76.55310621242486, 76.95390781563127, 77.35470941883769, 77.75551102204409, 78.15631262525051, 78.55711422845691, 78.95791583166333, 79.35871743486975, 79.75951903807616, 80.16032064128257, 80.56112224448898, 80.96192384769539, 81.3627254509018, 81.76352705410821, 82.16432865731463, 82.56513026052106, 82.96593186372746, 83.36673346693388, 83.76753507014028, 84.1683366733467, 84.5691382765531, 84.96993987975952, 85.37074148296593, 85.77154308617236, 86.17234468937876, 86.57314629258518, 86.97394789579158, 87.374749498998, 87.7755511022044, 88.17635270541082, 88.57715430861724, 88.97795591182366, 89.37875751503007, 89.77955911823648, 90.1803607214429, 90.5811623246493, 90.98196392785572, 91.38276553106212, 91.78356713426855, 92.18436873747495, 92.58517034068137, 92.98597194388778, 93.3867735470942, 93.7875751503006, 94.18837675350701, 94.58917835671343, 94.98997995991985, 95.39078156312627, 95.79158316633267, 96.19238476953909, 96.59318637274549, 96.99398797595191, 97.39478957915831, 97.79559118236473, 98.19639278557115, 98.59719438877757, 98.99799599198397, 99.39879759519039, 99.79959919839679, 100.20040080160321, 100.60120240480961, 101.00200400801603, 101.40280561122245, 101.80360721442887, 102.20440881763528, 102.60521042084169, 103.0060120240481, 103.40681362725451, 103.80761523046093, 104.20841683366734, 104.60921843687376, 105.01002004008016, 105.41082164328658, 105.81162324649299, 106.2124248496994, 106.6132264529058, 107.01402805611222, 107.41482965931864, 107.81563126252506, 108.21643286573148, 108.61723446893788, 109.0180360721443, 109.4188376753507, 109.81963927855712, 110.22044088176352, 110.62124248496995, 111.02204408817636, 111.42284569138278, 111.82364729458918, 112.2244488977956, 112.625250501002, 113.02605210420842, 113.42685370741484, 113.82765531062125, 114.22845691382767, 114.62925851703407, 115.03006012024049, 115.4308617234469, 115.83166332665331, 116.23246492985972, 116.63326653306615, 117.03406813627255, 117.43486973947897, 117.83567134268537, 118.23647294589179, 118.6372745490982, 119.03807615230461, 119.43887775551102, 119.83967935871743, 120.24048096192385, 120.64128256513025, 121.04208416833667, 121.4428857715431, 121.84368737474952, 122.24448897795592, 122.64529058116234, 123.04609218436875, 123.44689378757516, 123.84769539078157, 124.24849699398798, 124.64929859719439, 125.0501002004008, 125.45090180360721, 125.85170340681363, 126.25250501002004, 126.65330661322645, 127.05410821643287, 127.45490981963927, 127.85571142284572, 128.25651302605212, 128.65731462925854, 129.05811623246495, 129.45891783567134, 129.85971943887776, 130.26052104208418, 130.6613226452906, 131.06212424849699, 131.4629258517034, 131.86372745490982, 132.26452905811624, 132.66533066132266, 133.06613226452905, 133.46693386773546, 133.8677354709419, 134.26853707414833, 134.66933867735472, 135.07014028056113, 135.47094188376755, 135.87174348697397, 136.27254509018036, 136.67334669338678, 137.0741482965932, 137.4749498997996, 137.875751503006, 138.27655310621242, 138.67735470941884, 139.07815631262525, 139.47895791583167, 139.87975951903806, 140.2805611222445, 140.68136272545092, 141.08216432865734, 141.48296593186373, 141.88376753507015, 142.28456913827657, 142.68537074148298, 143.08617234468937, 143.4869739478958, 143.8877755511022, 144.28857715430863, 144.68937875751502, 145.09018036072143, 145.49098196392785, 145.89178356713427, 146.29258517034071, 146.6933867735471, 147.09418837675352, 147.49498997995994, 147.89579158316636, 148.29659318637275, 148.69739478957916, 149.09819639278558, 149.498997995992, 149.8997995991984, 150.3006012024048, 150.70140280561122, 151.10220440881764, 151.50300601202406, 151.90380761523045, 152.30460921843687, 152.7054108216433, 153.10621242484973, 153.50701402805612, 153.90781563126254, 154.30861723446895, 154.70941883767537, 155.11022044088176, 155.51102204408818, 155.9118236472946, 156.31262525050101, 156.7134268537074, 157.11422845691382, 157.51503006012024, 157.91583166332666, 158.31663326653307, 158.7174348697395, 159.1182364729459, 159.51903807615233, 159.91983967935874, 160.32064128256513, 160.72144288577155, 161.12224448897797, 161.5230460921844, 161.92384769539078, 162.3246492985972, 162.7254509018036, 163.12625250501003, 163.52705410821642, 163.92785571142284, 164.32865731462925, 164.72945891783567, 165.13026052104212, 165.5310621242485, 165.93186372745492, 166.33266533066134, 166.73346693386776, 167.13426853707415, 167.53507014028057, 167.93587174348698, 168.3366733466934, 168.7374749498998, 169.1382765531062, 169.53907815631263, 169.93987975951904, 170.34068136272546, 170.74148296593185, 171.1422845691383, 171.54308617234472, 171.94388777555113, 172.34468937875752, 172.74549098196394, 173.14629258517036, 173.54709418837678, 173.94789579158316, 174.34869739478958, 174.749498997996, 175.15030060120242, 175.5511022044088, 175.95190380761522, 176.35270541082164, 176.75350701402806, 177.15430861723448, 177.5551102204409, 177.9559118236473, 178.35671342685373, 178.75751503006015, 179.15831663326654, 179.55911823647295, 179.95991983967937, 180.3607214428858, 180.76152304609218, 181.1623246492986, 181.56312625250501, 181.96392785571143, 182.36472945891782, 182.76553106212424, 183.16633266533066, 183.5671342685371, 183.96793587174352, 184.3687374749499, 184.76953907815633, 185.17034068136275, 185.57114228456916, 185.97194388777555, 186.37274549098197, 186.7735470941884, 187.1743486973948, 187.5751503006012, 187.9759519038076, 188.37675350701403, 188.77755511022045, 189.17835671342687, 189.57915831663325, 189.9799599198397, 190.38076152304612, 190.78156312625254, 191.18236472945893, 191.58316633266534, 191.98396793587176, 192.38476953907818, 192.78557114228457, 193.18637274549098, 193.5871743486974, 193.98797595190382, 194.3887775551102, 194.78957915831663, 195.19038076152304, 195.59118236472946, 195.99198396793588, 196.3927855711423, 196.79358717434872, 197.19438877755513, 197.59519038076155, 197.99599198396794, 198.39679358717436, 198.79759519038078, 199.1983967935872, 199.59919839679358, 200 ], "y": [ 0, 0.004174654302551986, 0.01666410054864591, 0.03739008744840845, 0.06623906198192593, 0.10306216939924438, 0.1476752532203693, 0.19985885523526614, 0.2593582155038596, 0.3258832723560345, 0.39910866239163495, 0.4786737204804647, 0.5641824797622873, 0.6552118302013306, 0.7513286376886313, 0.8520249845761372, 0.9567765195621194, 1.0650507061962688, 1.1763068228796951, 1.2899959628649276, 1.4055610342559133, 1.5224367600080204, 1.640049677928034, 1.7578181406741602, 1.8751523157560233, 1.9914542759043967, 2.1061570851427938, 2.2186891001391755, 2.3284681853148363, 2.4349370807265394, 2.537563402066516, 2.635839640662466, 2.7292831634775565, 2.8174362131104242, 2.899865907795174, 2.9761642414013787, 3.045948083434079, 3.108859179033785, 3.164569458116879, 3.212809351409044, 3.2533416494184504, 3.2859669413439128, 3.310525646362785, 3.3268980136309567, 3.335004122282855, 3.3348038814314447, 3.3262970301682264, 3.3095231375632395, 3.284561602665058, 3.251531654500795, 3.210592164781084, 3.1619285143510982, 3.1057959552687064, 3.0424974193575083, 2.9723616488710944, 2.8957431964930525, 2.813022425336961, 2.724605508946396, 2.6309244312949245, 2.5324369867861094, 2.429626780253506, 2.323003226960667, 2.2131015526011333, 2.1004690332444964, 1.9856464221954664, 1.8692351224813981, 1.7518352018827388, 1.6340397545275862, 1.5164349008916975, 1.399599787798478, 1.2841065884189928, 1.1705205022719558, 1.0593997552237397, 0.9512955994883667, 0.8467523136275175, 0.7463065502218116, 0.6504487277832627, 0.5596565103994249, 0.47439357531311926, 0.3950889251030581, 0.3221368876838416, 0.25589711630596224, 0.19669458955579844, 0.14481961135561872, 0.1005278109635801, 0.06404014297373051, 0.035542887316004634, 0.015187649256228286, 0.003089870593059014, -0.0006923498448887528, 0.003845870048969742, 0.0166718919531046, 0.0377159226107929, 0.0668710138301172, 0.10399306248396763, 0.1489008105100405, 0.2013758449108399, 0.2611625977536734, 0.32796834617065806, 0.4014632123587177, 0.4812813236040767, 0.5670482510450845, 0.6583193859921922, 0.7546138398415716, 0.8554381003572653, 0.9602860316711959, 1.0686388742831596, 1.1799652450608307, 1.2937211372397524, 1.4093499204233486, 1.5262823405829185, 1.6439365200576366, 1.7617179575545485, 1.8790385751134975, 1.9953316164736221, 2.1099928936714707, 2.2224363990935156, 2.3320971996436235, 2.4384314367430475, 2.540916326330434, 2.639050158861809, 2.7323522993105964, 2.8203631871676045, 2.9026443364410337, 2.9787783356564685, 3.0483699890397578, 3.1110771091916987, 3.1665912601696293, 3.2146351860960833, 3.2549710739755384, 3.2874005536944124, 3.3117646980210655, 3.3279440226058004, 3.33585848598086, 3.335467489560431, 3.3267698776406407, 3.3098039373995585, 3.2846473988971954, 3.2514115552584024, 3.2102606074328395, 3.161422762671411, 3.105156241382888, 3.0417484040375555, 2.9715157511672143, 2.8948039233651746, 2.8119877012862617, 2.7234710056468123, 2.6296868972246807, 2.531097576859231, 2.4281943854513406, 2.3214946266353684, 2.211504117414163, 2.098791074258819, 1.9839417985894614, 1.8675402583180634, 1.7501680878484622, 1.632404588076354, 1.5148267263892952, 1.3980091366667022, 1.2825241192798496, 1.1689416410918696, 1.0578293354577644, 0.9497525022243851, 0.8452519287478123, 0.7448344689694432, 0.6490165090088866, 0.5582833006737489, 0.47308799588921374, 0.39385164669805284, 0.3209632052606103, 0.2547795238548122, 0.1956253548761654, 0.14379335083775668, 0.09954406437025287, 0.06310594822190063, 0.03467469242129988, 0.01439462405548211, 0.002355762096240638, -0.00139069099129447, 0.0031679322520952805, 0.016005764279716672, 0.03705840253896495, 0.06622290947132377, 0.10335781251236642, 0.148283104091751, 0.20078024163322716, 0.2605921475546324, 0.327423209267892, 0.40095571508237193, 0.4808251547822081, 0.5666048718156655, 0.6578507905114354, 0.7541018595367972, 0.854880051897638, 0.9596903649384397, 1.0680208203422823, 1.1793424641308436, 1.2931093666644, 1.408758622641829, 1.5257103511005956, 1.6433736082923789, 1.7611873632537156, 1.8785449284052818, 1.9948418536169494, 2.109490594265964, 2.221920511236943, 2.3315778709218797, 2.437925845220128, 2.540444511538427, 2.638630852790883, 2.7319987573989755, 2.8200790192915552, 2.9024193379048473, 2.978605537927733, 3.048268424301476, 3.1110560026919103, 3.166654657626676, 3.214789270159698, 3.2552232178711864, 3.287758374867639, 3.312235111781839, 3.3285322957728574, 3.3365672905260504, 3.33629595625306, 3.3277126496918163, 3.3108486285834675, 3.285778116098424, 3.2526450947150987, 3.211630323860932, 3.162946687680995, 3.106839195037998, 3.0435849795122785, 2.9734932994018095, 2.896905537722196, 2.814195202206677, 2.725767925306125, 2.6320614641890447, 2.5335457007415703, 2.4306938170150056, 2.32401887609444, 2.214081719544179, 2.1014456411996076, 1.9866763649499994, 1.8703420447385173, 1.7530132645622136, 1.6352630384720264, 1.5176668105727937, 1.4008024550232319, 1.2852502760359505, 1.1715930078774481, 1.0604073475824023, 0.9522257267070585, 0.8476021833724303, 0.7470703760058144, 0.6511349452363722, 0.5602715138951097, 0.4749266870148914, 0.3955180518304364, 0.3224341777783192, 0.25603461649696824, 0.1966499018266653, 0.14458154980955318, 0.10010205373894096, 0.06344030387244501, 0.034769492482053285, 0.014226089221641772, 0.001907174445384907, -0.0021295607922445503, 0.0021341852635248237, 0.014677325067762315, 0.03543938177583576, 0.0643204892434116, 0.10118139202645443, 0.14584344538122712, 0.19808861526429145, 0.2576655100756263, 0.32429240760715805, 0.397611498749232, 0.47723822190502135, 0.5627668489296183, 0.6537704851300363, 0.7498010692652111, 0.8503893735459986, 0.955045003635177, 1.0632563986474495, 1.174490831149429, 1.288194407159662, 1.4037922797864613, 1.5207283483408585, 1.6384158595470684, 1.7562432680017714, 1.8736114396067882, 1.9899336515690629, 2.1046355924006876, 2.217155361918887, 2.326943471246021, 2.433462842809586, 2.5361888103422134, 2.634609118881672, 2.7282239247708655, 2.8165553499516878, 2.899179107548262, 2.975676801581944, 3.045663148288935, 3.1087898339601248, 3.1647455149411194, 3.213255817632221, 3.2540833384884342, 3.287027644019466, 3.3119252707897235, 3.3286497254183174, 3.3371114845790593, 3.3372579051988813, 3.32907528737749, 3.312616519616163, 3.2879791066924846, 3.255295247584236, 3.214731835469394, 3.166490457726134, 3.1108073959328215, 3.0479536258680318, 2.9782348175105255, 2.901991335039263, 2.819598236833402, 2.7314652754722952, 2.638024518780837, 2.539724033647553, 2.4370829528938223, 2.330632046506959, 2.2209093323096334, 2.1084600759598735, 1.9938367909510661, 1.8775992386119553, 1.7603144281066427, 1.6425566164345806, 1.5249073084306024, 1.4079552567648754, 1.2922953738017666, 1.1784848186388226, 1.0670958654226048, 0.9587012326626384, 0.8538481926310293, 0.7530585713624646, 0.656828748654212, 0.5656296580661205, 0.47990678692061994, 0.4000801763027204, 0.32654442106001336, 0.25966866980267156, 0.19979662490344796, 0.14723872046150704, 0.10224710845078995, 0.06504228537503318, 0.03580606035610909, 0.014680571993028996, 0.0017682883619434542, -0.0028679929838581752, 0.0007948549860529175, 0.012740288779244185, 0.0329120943804225, 0.06121438725143355, 0.09751161233126222, 0.14162947351649868, 0.19336754138175327, 0.2524512330719712, 0.3185649712871895, 0.3913680614198731, 0.47049469155494106, 0.5555539324697502, 0.6461297376340981, 0.7417809432102243, 0.842041268052809, 0.9464193137089748, 1.0543985644182845, 1.1654373871127435, 1.278988636645467, 1.3945033110721337, 1.5113774445931525, 1.6290122647810499, 1.7468166448069686, 1.8642071034406686, 1.9806078050505334, 2.0954505596035404, 2.208174822665307, 2.3182276954000605, 2.425063924570643, 2.5281459025385145, 2.626945834415076, 2.7209864267517743, 2.809796440180246, 2.89292483061272, 2.969955406382759, 3.040506828245257, 3.1042326093764396, 3.160821115373866, 3.209995564256426, 3.2515140264643456, 3.285169424859174, 3.3107895347238006, 3.328236983762445, 3.3374114938133896, 3.338272385828168, 3.3308293844188777, 3.315129123569018, 3.2912550170555774, 3.259327258449034, 3.2195028211133545, 3.171975458205996, 3.116975702677904, 3.054770867273513, 2.985665044530747, 2.90999910678102, 2.8281477078552144, 2.740494182622777, 2.647496436944159, 2.5496408206517165, 2.4474257261459806, 2.341361588395656, 2.231970884937622, 2.1197881358769317, 2.0053599038868146, 1.8892447942086719, 1.7720134546520812, 1.6542485755947927, 1.5365448899827328, 1.4194831753756132, 1.3036329433056073, 1.18959435455113, 1.0779470000306164, 0.9692490557108262, 0.8640372826068762, 0.7628270267822306, 0.6661122193487021, 0.5743653764664521, 0.4880375993439898, 0.40755857423817354, 0.3333365724542088, 0.265755872344769, 0.2051449705100583, 0.1518025814179583, 0.10599375549153658, 0.06794420458889394, 0.037840302003164294, 0.015829082462513522, 0.002018242130144454, -0.0035238613957111475, -0.0007682090817868192, 0.010274879540216886, 0.029555746373632508, 0.05698539475682607, 0.0924430095172839, 0.13574875482660742, 0.1866647768598768, 0.24492430999388784, 0.3102318425108874, 0.3822631165985735, 0.4606651283500953, 0.5450561277640535, 0.6350256187445054, 0.7301343591009419, 0.8299143605483227, 0.9338688887070529, 1.0414787241280148, 1.1522344885496225, 1.2655612033831898, 1.3808731368885967, 1.4975872304033382, 1.6151230983425233, 1.732903028198875, 1.8503519805427306, 1.9668975890220424, 2.081970160362376, 2.1950026743669113, 2.3054307839164427, 2.412692814969386, 2.516259198393973, 2.615624014818845, 2.710278247421631, 2.799744945248929, 2.883579301220656, 2.9613686521300444, 3.0327324786436454, 3.0973224053013264, 3.1548222005162723, 3.2049477765749854, 3.2474471896372847, 3.282100639736308, 3.3087218056872567, 3.3271784774925703, 3.3373900484959376, 3.3393142844279518, 3.332947267187351, 3.318323394841015, 3.2955153816239666, 3.2646342579393717, 3.225829370358538, 3.1792883816209168, 3.1252372706341016, 3.0639403324738286, 2.9957001783839776 ] } ], "layout": { "template": { "data": { "bar": [ { "error_x": { "color": "rgb(36,36,36)" }, "error_y": { "color": "rgb(36,36,36)" }, "marker": { "line": { "color": "white", "width": 0.5 }, "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "bar" } ], "barpolar": [ { "marker": { "line": { "color": "white", "width": 0.5 }, "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "barpolar" } ], "carpet": [ { "aaxis": { "endlinecolor": "rgb(36,36,36)", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "rgb(36,36,36)" }, "baxis": { "endlinecolor": "rgb(36,36,36)", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "rgb(36,36,36)" }, "type": "carpet" } ], "choropleth": [ { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" }, "type": "choropleth" } ], "contour": [ { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" }, "colorscale": [ [ 0, "#440154" ], [ 0.1111111111111111, "#482878" ], [ 0.2222222222222222, "#3e4989" ], [ 0.3333333333333333, "#31688e" ], [ 0.4444444444444444, "#26828e" ], [ 0.5555555555555556, "#1f9e89" ], [ 0.6666666666666666, "#35b779" ], [ 0.7777777777777778, "#6ece58" ], [ 0.8888888888888888, "#b5de2b" ], [ 1, "#fde725" ] ], "type": "contour" } ], "contourcarpet": [ { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" }, "type": "contourcarpet" } ], "heatmap": [ { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" }, "colorscale": [ [ 0, "#440154" ], [ 0.1111111111111111, "#482878" ], [ 0.2222222222222222, "#3e4989" ], [ 0.3333333333333333, "#31688e" ], [ 0.4444444444444444, "#26828e" ], [ 0.5555555555555556, "#1f9e89" ], [ 0.6666666666666666, "#35b779" ], [ 0.7777777777777778, "#6ece58" ], [ 0.8888888888888888, "#b5de2b" ], [ 1, "#fde725" ] ], "type": "heatmap" } ], "heatmapgl": [ { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" }, "colorscale": [ [ 0, "#440154" ], [ 0.1111111111111111, "#482878" ], [ 0.2222222222222222, "#3e4989" ], [ 0.3333333333333333, "#31688e" ], [ 0.4444444444444444, "#26828e" ], [ 0.5555555555555556, "#1f9e89" ], [ 0.6666666666666666, "#35b779" ], [ 0.7777777777777778, "#6ece58" ], [ 0.8888888888888888, "#b5de2b" ], [ 1, "#fde725" ] ], "type": "heatmapgl" } ], "histogram": [ { "marker": { "line": { "color": "white", "width": 0.6 } }, "type": "histogram" } ], "histogram2d": [ { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" }, "colorscale": [ [ 0, "#440154" ], [ 0.1111111111111111, "#482878" ], [ 0.2222222222222222, "#3e4989" ], [ 0.3333333333333333, "#31688e" ], [ 0.4444444444444444, "#26828e" ], [ 0.5555555555555556, "#1f9e89" ], [ 0.6666666666666666, "#35b779" ], [ 0.7777777777777778, "#6ece58" ], [ 0.8888888888888888, "#b5de2b" ], [ 1, "#fde725" ] ], "type": "histogram2d" } ], "histogram2dcontour": [ { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" }, "colorscale": [ [ 0, "#440154" ], [ 0.1111111111111111, "#482878" ], [ 0.2222222222222222, "#3e4989" ], [ 0.3333333333333333, "#31688e" ], [ 0.4444444444444444, "#26828e" ], [ 0.5555555555555556, "#1f9e89" ], [ 0.6666666666666666, "#35b779" ], [ 0.7777777777777778, "#6ece58" ], [ 0.8888888888888888, "#b5de2b" ], [ 1, "#fde725" ] ], "type": "histogram2dcontour" } ], "mesh3d": [ { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" }, "type": "mesh3d" } ], "parcoords": [ { "line": { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" } }, "type": "parcoords" } ], "pie": [ { "automargin": true, "type": "pie" } ], "scatter": [ { "fillpattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 }, "type": "scatter" } ], "scatter3d": [ { "line": { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" } }, "marker": { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" } }, "type": "scatter3d" } ], "scattercarpet": [ { "marker": { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" } }, "type": "scattercarpet" } ], "scattergeo": [ { "marker": { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" } }, "type": "scattergeo" } ], "scattergl": [ { "marker": { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" } }, "type": "scattergl" } ], "scattermapbox": [ { "marker": { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" } }, "type": "scattermapbox" } ], "scatterpolar": [ { "marker": { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" } }, "type": "scatterpolar" } ], "scatterpolargl": [ { "marker": { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" } }, "type": "scatterpolargl" } ], "scatterternary": [ { "marker": { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" } }, "type": "scatterternary" } ], "surface": [ { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" }, "colorscale": [ [ 0, "#440154" ], [ 0.1111111111111111, "#482878" ], [ 0.2222222222222222, "#3e4989" ], [ 0.3333333333333333, "#31688e" ], [ 0.4444444444444444, "#26828e" ], [ 0.5555555555555556, "#1f9e89" ], [ 0.6666666666666666, "#35b779" ], [ 0.7777777777777778, "#6ece58" ], [ 0.8888888888888888, "#b5de2b" ], [ 1, "#fde725" ] ], "type": "surface" } ], "table": [ { "cells": { "fill": { "color": "rgb(237,237,237)" }, "line": { "color": "white" } }, "header": { "fill": { "color": "rgb(217,217,217)" }, "line": { "color": "white" } }, "type": "table" } ] }, "layout": { "annotationdefaults": { "arrowhead": 0, "arrowwidth": 1 }, "autotypenumbers": "strict", "coloraxis": { "colorbar": { "outlinewidth": 1, "tickcolor": "rgb(36,36,36)", "ticks": "outside" } }, "colorscale": { "diverging": [ [ 0, "rgb(103,0,31)" ], [ 0.1, "rgb(178,24,43)" ], [ 0.2, "rgb(214,96,77)" ], [ 0.3, "rgb(244,165,130)" ], [ 0.4, "rgb(253,219,199)" ], [ 0.5, "rgb(247,247,247)" ], [ 0.6, "rgb(209,229,240)" ], [ 0.7, "rgb(146,197,222)" ], [ 0.8, "rgb(67,147,195)" ], [ 0.9, "rgb(33,102,172)" ], [ 1, "rgb(5,48,97)" ] ], "sequential": [ [ 0, "#440154" ], [ 0.1111111111111111, "#482878" ], [ 0.2222222222222222, "#3e4989" ], [ 0.3333333333333333, "#31688e" ], [ 0.4444444444444444, "#26828e" ], [ 0.5555555555555556, "#1f9e89" ], [ 0.6666666666666666, "#35b779" ], [ 0.7777777777777778, "#6ece58" ], [ 0.8888888888888888, "#b5de2b" ], [ 1, "#fde725" ] ], "sequentialminus": [ [ 0, "#440154" ], [ 0.1111111111111111, "#482878" ], [ 0.2222222222222222, "#3e4989" ], [ 0.3333333333333333, "#31688e" ], [ 0.4444444444444444, "#26828e" ], [ 0.5555555555555556, "#1f9e89" ], [ 0.6666666666666666, "#35b779" ], [ 0.7777777777777778, "#6ece58" ], [ 0.8888888888888888, "#b5de2b" ], [ 1, "#fde725" ] ] }, "colorway": [ "#1F77B4", "#FF7F0E", "#2CA02C", "#D62728", "#9467BD", "#8C564B", "#E377C2", "#7F7F7F", "#BCBD22", "#17BECF" ], "font": { "color": "rgb(36,36,36)" }, "geo": { "bgcolor": "white", "lakecolor": "white", "landcolor": "white", "showlakes": true, "showland": true, "subunitcolor": "white" }, "hoverlabel": { "align": "left" }, "hovermode": "closest", "mapbox": { "style": "light" }, "paper_bgcolor": "white", "plot_bgcolor": "white", "polar": { "angularaxis": { "gridcolor": "rgb(232,232,232)", "linecolor": "rgb(36,36,36)", "showgrid": false, "showline": true, "ticks": "outside" }, "bgcolor": "white", "radialaxis": { "gridcolor": "rgb(232,232,232)", "linecolor": "rgb(36,36,36)", "showgrid": false, "showline": true, "ticks": "outside" } }, "scene": { "xaxis": { "backgroundcolor": "white", "gridcolor": "rgb(232,232,232)", "gridwidth": 2, "linecolor": "rgb(36,36,36)", "showbackground": true, "showgrid": false, "showline": true, "ticks": "outside", "zeroline": false, "zerolinecolor": "rgb(36,36,36)" }, "yaxis": { "backgroundcolor": "white", "gridcolor": "rgb(232,232,232)", "gridwidth": 2, "linecolor": "rgb(36,36,36)", "showbackground": true, "showgrid": false, "showline": true, "ticks": "outside", "zeroline": false, "zerolinecolor": "rgb(36,36,36)" }, "zaxis": { "backgroundcolor": "white", "gridcolor": "rgb(232,232,232)", "gridwidth": 2, "linecolor": "rgb(36,36,36)", "showbackground": true, "showgrid": false, "showline": true, "ticks": "outside", "zeroline": false, "zerolinecolor": "rgb(36,36,36)" } }, "shapedefaults": { "fillcolor": "black", "line": { "width": 0 }, "opacity": 0.3 }, "ternary": { "aaxis": { "gridcolor": "rgb(232,232,232)", "linecolor": "rgb(36,36,36)", "showgrid": false, "showline": true, "ticks": "outside" }, "baxis": { "gridcolor": "rgb(232,232,232)", "linecolor": "rgb(36,36,36)", "showgrid": false, "showline": true, "ticks": "outside" }, "bgcolor": "white", "caxis": { "gridcolor": "rgb(232,232,232)", "linecolor": "rgb(36,36,36)", "showgrid": false, "showline": true, "ticks": "outside" } }, "title": { "x": 0.05 }, "xaxis": { "automargin": true, "gridcolor": "rgb(232,232,232)", "linecolor": "rgb(36,36,36)", "showgrid": false, "showline": true, "ticks": "outside", "title": { "standoff": 15 }, "zeroline": false, "zerolinecolor": "rgb(36,36,36)" }, "yaxis": { "automargin": true, "gridcolor": "rgb(232,232,232)", "linecolor": "rgb(36,36,36)", "showgrid": false, "showline": true, "ticks": "outside", "title": { "standoff": 15 }, "zeroline": false, "zerolinecolor": "rgb(36,36,36)" } } }, "xaxis": { "title": { "text": "t (ns)" } }, "yaxis": { "title": { "text": "y (m)" } } } }, "text/html": [ "