{ "cells": [ { "cell_type": "markdown", "id": "c2e0135b-da8a-42cb-bf44-a7891a5fa7e4", "metadata": {}, "source": [ "# Mono Calibration and EXAFS Results\n", "\n", "This notebook explores the effect of energy miscalibration on EXAFS results.\n", "\n", "EXAFS measurements need well-calibrated energies, typically coming from a double-crystal monochromator. While EXAFS can be measured in modes without a double-crystal monochromator, that is sort of unusual, and out of the scope of the discussion here. \n", "\n", "A double-crystal monochromator uses Bragg's Law to relate the monochromator angle $\\theta$ to X-ray energy $E$ with\n", "\n", " $\\lambda = \\frac{hc}{E} = 2d\\sin(\\theta) $\n", "\n", "With $h$ being Planck's constant, $c$ the speed of light, and $d$ the lattice spacing of the monochromator crystal. With $d$ in Angstroms and $E$ in eV, $hc \\approx 12398.419$ eV A. The crystal $d$ values for Si, Ge, and C (the most common crystals used) are typically very well known, but there may be thermally induced changes as synchrotron beamlines. That is, both the values of $d$ and an angular offset, $\\theta_0$, may need to be adjusted to calibrate a monochromator. This is typically done by scanning across known (or at least reproducible) energies, say the $K$-edges of metal foils, and adjusting $d$ and/or $\\theta_0$ until the energy is correct. \n", "\n", "Most commonly, the XANES of some metal foil is measured and the maximum of the first derivative is set to the tabulated edge energy for the metal. This method will give a bit of a variation as the energy resolution changes. These days, most good XAFS beamlines have high enough energy resolution that this is not much of a problem, but it can be a problem for older data sets. It should also be said that those tabulated values are not always accurate to better than 1 or 2 eV. In fact, rather than use the tabulated values such as at https://xraydb.xrayabsorption.org/ the values from Kraft et al Review of Scientific Instruments 67, p681 (1996): (https://doi.org/10.1063/1.1146657) should be used, as they were carefully measured with a single, very high-resolution and consistently calibrated monochromator.\n", "\n", "\n", "Ideally, a well-calibrated monochromator would have a single $d$-spacing (and so good thermal stability) and angular offset that stays calibrated across the edges of many edges. If the energy needs to be recalibrated very frequently or for each edge, then the energy will be drifting between those edges, and the reason for any need to recalibrate should be investigated. If the beamline is well collimated and the angle of the beam incident on the monochromator is stable, it is certainly possible to have a monochromator set up so that stays calibrated for weeks and over an energy range of 10 keV or more. This can be more complicated when a collimating (or any other) mirror is placed before the monochromator, as then very high stability of that mirror is also needed.\n", "\n", "For any XAS analysis, the energy scale is critical. Because of the variations above in calibrating energy, small energy shifts between beamlines or even runs at the same beamline but different months are not uncommon. If the shift is relatively small, say 1 to 5 eV at 10 keV, simply adding a constant energy offset is satisfactory.\n", " \n", "You may hear people say that poor calibration leads to large errors in EXAFS results or even leads to \"the wrong results\". You should know that when people describe scientific measurements and analysis as Wrong, they are almost certainly wrong. Measurements have precision and accuracy, and science always allows for improvements in measurement and analysis.\n", "\n", "Here, we'll look at these effects. From this, we will conclude that claims that poor calibration leads to large errors in EXAFS results are not accurate." ] }, { "cell_type": "code", "execution_count": 1, "id": "7c518804", "metadata": {}, "outputs": [], "source": [ "# We start with some imports of Python functions we will need to read, \n", "# modify, and fit xafs data:\n", "\n", "import time, os, requests\n", "import numpy as np\n", "import scipy.constants as consts\n", "\n", "from larch.io import read_ascii\n", "from larch import Group\n", "from larch.fitting import param_group, param, guess\n", "from larch.xafs import (autobk, feffpath, feffit_transform,\n", " feffit_dataset, feffit, feffit_report)\n", "from larch.plot.plotly_xafsplots import plot, multi_plot\n", "\n", "PLANCK_HC = 1.e10 * consts.Planck * consts.c / consts.e\n", "RAD2DEG = 180.0/np.pi\n", "DEG2RAD = np.pi/180.0" ] }, { "cell_type": "code", "execution_count": 2, "id": "9905f684", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Got data files: ../xafsdata/cu_10k.xmu ../feffit/feff0001.dat\n" ] } ], "source": [ "# then we will some XAFS data (Cu metal) at the Cu K edge, and a Feff.dat file\n", "# so we'll download these to the current directory\n", "cudata_file = '../xafsdata/cu_10k.xmu'\n", "feffdat_file = '../feffit/feff0001.dat'\n", "\n", "root_url = 'https://raw.githubusercontent.com/xraypy/xraylarch/master/examples'\n", "current_dir = os.getcwd()\n", "downloads = {}\n", "if not os.path.exists(cudata_file):\n", " cudata_file = 'cu_10k.xmu'\n", " downloads[cudata_file] = f'{root_url}/xafsdata/cu_10k.xmu'\n", "if not os.path.exists(feffdat_file):\n", " feffdat_file = 'feff0001.dat'\n", " downloads[feffdat_file] = f'{root_url}/feffit/feff0001.dat'\n", "\n", "for dest, source_url in downloads.items():\n", " t0 = time.time()\n", " req = requests.get(source_url, verify=True, timeout=1.0)\n", " print(source_url, req)\n", " if req.status_code == 200:\n", " full_dest = os.path.join(current_dir, dest)\n", " with open(full_dest, 'wb') as fh:\n", " fh.write(req.content)\n", " print(f'Downloaded {full_dest}: {time.time()-t0:.3f} sec')\n", " time.sleep(0.25)\n", "time.sleep(0.1)\n", "print(\"Got data files: \", cudata_file, feffdat_file)" ] }, { "cell_type": "code", "execution_count": 3, "id": "2aedfd2a", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Using a d spacing 3.135253 Ang\n" ] }, { "data": { "text/html": [ " \n", " " ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/vnd.plotly.v1+json": { "config": { "plotlyServerURL": "https://plot.ly" }, "data": [ { "line": { "color": "red", "width": 2 }, "name": "angle", "type": "scatter", "uid": "307c6074-aa98-4923-8276-5b07f2a67c4e", "x": [ 8786.204, 8796.258, 8806.253, 8816.27, 8826.27, 8836.292, 8846.337, 8856.322, 8866.372, 8876.361, 8886.416, 8896.409, 8906.425, 8916.464, 8926.483, 8936.483, 8936.998, 8937.471, 8937.944, 8938.417, 8938.976, 8939.406, 8939.878, 8940.352, 8940.824, 8941.384, 8941.814, 8942.287, 8942.761, 8943.319, 8943.793, 8944.224, 8944.697, 8945.257, 8945.73, 8946.204, 8946.635, 8947.108, 8947.668, 8948.143, 8948.616, 8949.048, 8949.564, 8950.038, 8950.513, 8950.987, 8951.461, 8951.979, 8952.453, 8952.927, 8953.401, 8953.876, 8954.394, 8954.868, 8955.343, 8955.817, 8956.336, 8956.811, 8957.285, 8957.761, 8958.235, 8958.71, 8959.229, 8959.703, 8960.179, 8960.654, 8961.13, 8961.647, 8962.123, 8962.599, 8963.074, 8963.55, 8964.068, 8964.544, 8965.02, 8965.496, 8965.972, 8966.447, 8966.967, 8967.442, 8967.919, 8968.395, 8968.871, 8969.347, 8969.909, 8970.386, 8970.862, 8971.339, 8971.772, 8972.248, 8972.812, 8973.288, 8973.765, 8974.198, 8974.675, 8975.151, 8975.628, 8976.192, 8976.626, 8977.103, 8977.58, 8978.057, 8978.534, 8979.055, 8979.531, 8980.009, 8980.486, 8980.964, 8981.441, 8981.919, 8982.396, 8982.918, 8983.396, 8983.873, 8984.352, 8984.829, 8985.307, 8985.785, 8986.307, 8986.784, 8987.263, 8987.741, 8988.219, 8988.697, 8989.176, 8989.653, 8990.132, 8990.698, 8991.133, 8991.611, 8992.09, 8992.568, 8993.048, 8993.526, 8994.006, 8994.484, 8995.007, 8995.485, 8995.965, 8996.444, 8996.924, 8997.402, 8997.882, 8998.361, 8998.84, 8999.319, 8999.799, 9000.366, 9000.803, 9001.281, 9001.762, 9002.241, 9002.721, 9003.201, 9003.681, 9004.161, 9004.642, 9005.121, 9005.602, 9006.125, 9006.605, 9007.086, 9007.566, 9008.047, 9008.526, 9009.007, 9009.487, 9009.969, 9011.062, 9012.153, 9013.246, 9014.384, 9015.565, 9016.746, 9017.928, 9019.154, 9020.38, 9021.606, 9022.877, 9024.191, 9025.464, 9026.822, 9028.14, 9029.499, 9030.903, 9032.309, 9033.759, 9035.164, 9036.615, 9038.109, 9039.604, 9041.145, 9042.642, 9044.271, 9045.768, 9047.397, 9049.029, 9050.616, 9052.336, 9053.969, 9055.646, 9057.324, 9059.047, 9060.814, 9062.627, 9064.353, 9066.166, 9067.981, 9069.842, 9071.702, 9073.563, 9075.47, 9077.377, 9079.329, 9081.283, 9083.237, 9085.281, 9087.237, 9089.283, 9091.33, 9093.378, 9095.472, 9097.61, 9099.705, 9101.847, 9104.033, 9106.221, 9108.409, 9110.644, 9112.879, 9115.16, 9117.398, 9119.728, 9122.057, 9124.389, 9126.721, 9129.099, 9131.478, 9133.904, 9136.331, 9138.805, 9141.279, 9143.756, 9146.233, 9148.803, 9151.328, 9153.899, 9156.474, 9159.094, 9161.716, 9164.339, 9167.009, 9169.681, 9172.399, 9175.074, 9177.842, 9180.564, 9183.381, 9186.153, 9188.973, 9191.794, 9194.663, 9197.533, 9200.406, 9203.326, 9206.247, 9209.217, 9212.143, 9215.115, 9218.183, 9221.16, 9224.185, 9227.258, 9230.333, 9233.502, 9236.581, 9239.755, 9242.838, 9246.016, 9249.196, 9252.425, 9255.655, 9258.889, 9262.171, 9265.454, 9268.787, 9272.122, 9275.459, 9278.8, 9282.188, 9285.58, 9289.021, 9292.465, 9295.91, 9299.404, 9302.949, 9306.403, 9309.953, 9313.506, 9317.061, 9320.666, 9324.32, 9327.932, 9331.592, 9335.256, 9338.922, 9342.639, 9346.358, 9350.128, 9353.9, 9357.677, 9361.455, 9365.285, 9369.165, 9373.048, 9376.84, 9380.777, 9384.671, 9388.614, 9392.563, 9396.561, 9400.515, 9404.567, 9408.576, 9412.637, 9416.7, 9420.768, 9424.935, 9429.01, 9433.184, 9437.313, 9441.495, 9445.729, 9449.966, 9454.208, 9458.453, 9462.702, 9467.004, 9471.358, 9475.668, 9480.029, 9484.396, 9488.814, 9493.189, 9497.665, 9502.097, 9506.581, 9511.118, 9515.563, 9520.109, 9524.611, 9529.265, 9533.775, 9538.438, 9543.007, 9547.679, 9552.306, 9556.986, 9561.673, 9566.462, 9571.157, 9575.956, 9580.661, 9585.47, 9590.283, 9595.102, 9599.976, 9604.854, 9609.737, 9614.625, 9619.568, 9624.518, 9629.471, 9634.431, 9639.445, 9644.466, 9649.541, 9654.571, 9659.708, 9664.75, 9669.847, 9675.001, 9680.159, 9685.324, 9690.545, 9695.722, 9700.953, 9706.243, 9711.486, 9716.788, 9722.094, 9727.407, 9732.776, 9738.204, 9743.535, 9748.923, 9754.317, 9759.769, 9765.227, 9770.742, 9776.212, 9781.793, 9787.275, 9792.816, 9798.416, 9804.022, 9809.584, 9815.254, 9820.88, 9826.564, 9832.308, 9837.953, 9843.659, 9849.422, 9855.192, 9860.971, 9866.755, 9872.6, 9878.45, 9884.31, 9890.175, 9896.1, 9901.98, 9907.92, 9913.867, 9919.876, 9925.838, 9931.859, 9937.89, 9943.98, 9950.025, 9956.132, 9962.3, 9968.367, 9974.497, 9980.687, 9986.832, 9993.038, 9999.253, 10005.53, 10011.76, 10018.05, 10024.35, 10030.66, 10037.03, 10043.36, 10049.74, 10056.14, 10062.54, 10068.96, 10075.43, 10081.92, 10088.46, 10094.91, 10101.42, 10107.94, 10114.52, 10121.11, 10127.71, 10134.32, 10141.05, 10147.68, 10154.31, 10161.01, 10167.78, 10174.49, 10181.28, 10188.07, 10194.82, 10201.63, 10208.45, 10215.28, 10222.17, 10229.13, 10235.99, 10242.86, 10249.85, 10256.79, 10263.8, 10270.82, 10277.85, 10284.89, 10291.88, 10299.05, 10306.12, 10313.2, 10320.41, 10327.5, 10334.73, 10341.85, 10349.09, 10356.35, 10363.61, 10370.89, 10378.17, 10385.47, 10392.78, 10400.1, 10407.48, 10414.82, 10422.23, 10429.65, 10437.14, 10444.58, 10452.03, 10459.55, 10467.08, 10474.63, 10482.18, 10489.75, 10497.32, 10504.97, 10512.57, 10520.24, 10527.92, 10535.62, 10543.32, 10551.04, 10558.77, 10566.5, 10574.32, 10582.2, 10589.91, 10597.76, 10605.62, 10613.55, 10621.43, 10629.33, 10637.3, 10645.21, 10653.21, 10661.21, 10669.23, 10677.26, 10685.3, 10693.35, 10701.54, 10709.56, 10717.71, 10725.81, 10733.99, 10742.24, 10750.38, 10758.6, 10766.83, 10775.14, 10783.39, 10791.72, 10800, 10808.36, 10816.73, 10825.05, 10833.45, 10841.86, 10850.41, 10858.78, 10867.23, 10875.76, 10884.23, 10892.85, 10901.29, 10909.94, 10918.47, 10927.08, 10935.7, 10944.4, 10952.99, 10961.65, 10970.4, 10979.09, 10987.86, 10996.65, 11005.32, 11014.14, 11022.97, 11031.81, 11040.67, 11049.54, 11058.43, 11067.33, 11076.25, 11085.18, 11094.12, 11103.15, 11112.13, 11121.25, 11130.19, 11139.34, 11148.38, 11157.43, 11166.56, 11175.78, 11184.87, 11194.05, 11203.31, 11212.45, 11221.68, 11230.98, 11240.17, 11249.51, 11258.8, 11268.17, 11277.55, 11286.82, 11296.24, 11305.6, 11315.05, 11324.52, 11333.93, 11343.43, 11352.94, 11362.47 ], "y": [ 13.005314993100693, 12.99018989542178, 12.975188699245399, 12.960189514821664, 12.945250642722556, 12.930313724478918, 12.915377381346755, 12.9005647158221, 12.885690195476501, 12.870940207096115, 12.856127124395691, 12.841439418801889, 12.826751839185006, 12.812064500675222, 12.797440210482673, 12.782877194707783, 12.782128104535529, 12.78144018322484, 12.78075233659776, 12.780064564642021, 12.779251839493499, 12.77862673724659, 12.77794064951247, 12.777251729405881, 12.776565790586398, 12.775752061213693, 12.775127304259154, 12.774440142790285, 12.773751603344179, 12.772941140186147, 12.772252763715159, 12.771626899858974, 12.770940118149888, 12.770127111378503, 12.769440492348897, 12.76875249639954, 12.768126978490766, 12.767440576381365, 12.766628018949788, 12.765938877894609, 12.76525271305077, 12.764626090499938, 12.763877705992412, 12.763190314497775, 12.762501547719905, 12.761814305717149, 12.761127138363147, 12.760376268672317, 12.759689257505544, 12.759002320949591, 12.758315458992199, 12.757627222781439, 12.756876768446714, 12.756190137289483, 12.75550213232777, 12.754815650412302, 12.754064081634928, 12.753376307962222, 12.752690056810886, 12.752000985043116, 12.751314883165694, 12.750627408513205, 12.749876337470374, 12.749190466163826, 12.748501775780797, 12.747814607030602, 12.747126066558781, 12.746378304015446, 12.745689920022826, 12.745001611020292, 12.744314822786702, 12.743626663569788, 12.742877869586623, 12.74218986687526, 12.741501939091144, 12.740814086221935, 12.740126308255283, 12.739440049854723, 12.738688862994097, 12.738002760770659, 12.73731384470511, 12.736626447826438, 12.735939125775213, 12.735251878539104, 12.73444056108046, 12.733752033573541, 12.733065024374989, 12.732376646879159, 12.731751832516034, 12.731065040977239, 12.730251376300776, 12.729564748016804, 12.728876752172892, 12.728252284202377, 12.727564431434699, 12.72687809544987, 12.72619039247091, 12.725377356095436, 12.724751793051453, 12.724064321809799, 12.723376925484347, 12.722689604062749, 12.722002357532654, 12.721251802750587, 12.720566153233491, 12.719877697846343, 12.719190757661885, 12.718502452415388, 12.717815662032761, 12.717127506877375, 12.716440866247153, 12.715689533892585, 12.715001610770143, 12.714315201649251, 12.71362598971456, 12.712939730376966, 12.712252107292148, 12.711564559225032, 12.710813807893533, 12.710127854705732, 12.709439100540225, 12.708751859320333, 12.70806469305512, 12.707377601732214, 12.706689148145003, 12.706003643863882, 12.705315340412884, 12.704502118448046, 12.703877187052136, 12.703190552327749, 12.702502556218723, 12.701816071347798, 12.70112678945542, 12.700440454544987, 12.69975132321567, 12.699065138216167, 12.698314439881834, 12.69762841146912, 12.69693958786909, 12.696252274451904, 12.69556360142048, 12.694877872733903, 12.694189350064372, 12.693502336901085, 12.692815398731137, 12.692128535542128, 12.691440313604227, 12.690627448408813, 12.690001026076652, 12.689315903272421, 12.68862655584346, 12.68794014978583, 12.687252385857938, 12.68656469712373, 12.685877083570764, 12.68518954518661, 12.68450064982222, 12.68381469387501, 12.683125949099141, 12.68237715004626, 12.681689994203895, 12.681001482114802, 12.68031447660267, 12.679626115132068, 12.67894069077081, 12.678252479712201, 12.677565774555008, 12.676876283614446, 12.675313050116976, 12.67375306476755, 12.672190607950471, 12.67056423587987, 12.668876855498825, 12.667189928305124, 12.665502026301294, 12.663751771299186, 12.662002004075603, 12.660252724424911, 12.65843975246154, 12.656565995143357, 12.65475123705771, 12.652815883499645, 12.650938106719481, 12.649002505035053, 12.647003437789822, 12.645002161353615, 12.6429389249843, 12.640940367703097, 12.638877045952023, 12.636753287851272, 12.634628828587187, 12.632439755157243, 12.63031391848581, 12.628001452908274, 12.625877122245488, 12.623566294592766, 12.621252066478062, 12.619002470259682, 12.616565257566148, 12.614252200884325, 12.611877710549594, 12.609502706403719, 12.607064948704457, 12.604565925362984, 12.602002883590181, 12.599563810366375, 12.597002819352484, 12.594440054559225, 12.59181342973909, 12.589189319984118, 12.586564902968727, 12.58387676027102, 12.581189775173673, 12.578440582664898, 12.575689786502064, 12.572940203292593, 12.570065273188218, 12.567315358078938, 12.56444021109142, 12.561564986058615, 12.558689683998873, 12.555751171969165, 12.552752344957396, 12.549815231932108, 12.546813659190414, 12.54375192178685, 12.540688891578622, 12.537627369355553, 12.534501638927209, 12.531377479598419, 12.528190638594955, 12.525065461965527, 12.521813485092622, 12.518564605054678, 12.51531324282342, 12.512063583049352, 12.508751573688674, 12.505439940125768, 12.502064701494485, 12.498689909249862, 12.495251652773351, 12.491815303045211, 12.488376694976536, 12.484939995077516, 12.481376278110298, 12.47787695795754, 12.4743159199373, 12.470751395076995, 12.467126684958385, 12.463501334108358, 12.459876726956509, 12.4561893549187, 12.452501423005268, 12.448752259435826, 12.445064631026776, 12.44125111532392, 12.437503271674425, 12.433627021556203, 12.429815069452962, 12.425939526780661, 12.422065046802254, 12.418127139211558, 12.414190376619135, 12.410252018463959, 12.406251812415908, 12.402252836684234, 12.398189441693564, 12.394188869181011, 12.390128066725318, 12.385938907286844, 12.38187673159242, 12.377751809439488, 12.373564269925179, 12.369376863162543, 12.36506443949462, 12.360877391540027, 12.356564146809752, 12.352377468933648, 12.348064774617246, 12.343752403723576, 12.339376690083162, 12.335002748881667, 12.33062652151394, 12.326188539885763, 12.321752426579419, 12.317252042746073, 12.312752273952393, 12.308253122469518, 12.303751897019287, 12.29919073761235, 12.294627606673153, 12.290002043825123, 12.285375959768272, 12.280752043051255, 12.276065940209978, 12.271315119487998, 12.266689815650603, 12.261939617648512, 12.257189116824653, 12.252439653669137, 12.247627177741698, 12.242753176402024, 12.237939039352087, 12.233064817900118, 12.228189187635042, 12.223314814791577, 12.218376628679934, 12.213439808567173, 12.208439390284221, 12.203440447589692, 12.198439010989613, 12.193440381900196, 12.188377164854504, 12.183252167767368, 12.17812755495401, 12.173127231516448, 12.167940080604533, 12.162813964468732, 12.157627777427633, 12.152438164150432, 12.14718870356654, 12.142001509331067, 12.136690380274311, 12.131440219344704, 12.12612662449287, 12.120815104984507, 12.115501744807139, 12.110063942220238, 12.104750951961725, 12.099313753675139, 12.093940014746588, 12.08850219925739, 12.083001786357794, 12.077502523238767, 12.072001821362177, 12.06650228145395, 12.061002613460825, 12.055439489484307, 12.049814384759184, 12.044251332611166, 12.038627718808282, 12.033001667893112, 12.027315302393376, 12.021689615796248, 12.01593954144686, 12.010251449933126, 12.004502140284913, 11.998690517096778, 11.993002236194021, 11.987190323372829, 11.981440253948977, 11.975501886157502, 11.969752913953279, 11.963814754771246, 11.958002056271926, 11.952064202571139, 11.94618939429882, 11.94025320912262, 11.934314099216035, 11.928251885660044, 11.922314683580325, 11.916252116027248, 11.910314325976938, 11.904251443868812, 11.898189742016484, 11.89212671127019, 11.886000811079127, 11.879876246142187, 11.873751768530012, 11.867627387798134, 11.861440564350948, 11.85525148827542, 11.849065170228037, 11.842876623551113, 11.836627316849489, 11.830375939910386, 11.824064086382284, 11.817814891188187, 11.81143962508341, 11.805188992922801, 11.798876946573863, 11.79250122238885, 11.786127495230422, 11.779752068784907, 11.773314575807284, 11.766938329028216, 11.760502635424476, 11.754001561582179, 11.747565387315404, 11.741064002147445, 11.734564964530911, 11.728064610246713, 11.721503105821844, 11.714877010024903, 11.708376664880168, 11.70181419551368, 11.69525183590304, 11.688626442269692, 11.682001326022867, 11.675314699755424, 11.668690244555302, 11.661939161169563, 11.655315488864819, 11.648628225540294, 11.641877604639046, 11.635127639006853, 11.62843843797001, 11.62162731803138, 11.614876991018741, 11.608065088181775, 11.601189448954512, 11.594440301158938, 11.587626254122943, 11.580752320399439, 11.573878260252195, 11.567001710331311, 11.560127445349178, 11.553189035777592, 11.546253083879941, 11.539313678424891, 11.532376757304716, 11.525377392974962, 11.518439641605354, 11.51143963145188, 11.504439949387407, 11.49737599348888, 11.490375916596188, 11.483315272321256, 11.476251655114455, 11.469127809370976, 11.46206540428259, 11.454939451632232, 11.44775137098873, 11.440689849880346, 11.433563904793083, 11.426377275255538, 11.419251882742765, 11.412064833540711, 11.404876484288472, 11.397625673594517, 11.390438327132513, 11.383191013412045, 11.37594147905347, 11.368689750323485, 11.361378502166518, 11.354122537003466, 11.346818693261556, 11.339501452110898, 11.332193705853928, 11.324872643635842, 11.317504192115633, 11.3101226546699, 11.302694047958477, 11.295377285448804, 11.28800212354198, 11.280625344034853, 11.273190513832125, 11.26575426606794, 11.258316625896091, 11.250877618411502, 11.24331372112887, 11.235872218640504, 11.22844062525282, 11.220940614872049, 11.213372483220246, 11.205881558293472, 11.198311571007986, 11.190751871602181, 11.183246882786973, 11.175685443810355, 11.168123209201383, 11.160560203086629, 11.152941196926678, 11.145255407459437, 11.137690472141383, 11.130124861061766, 11.122437709902995, 11.114816111130487, 11.107128302575124, 11.09944024619948, 11.091751965247209, 11.08406348289816, 11.076440212383993, 11.068631591351524, 11.060942719311077, 11.053253737901873, 11.045434623262063, 11.037756495713742, 11.029937810221728, 11.022248967924323, 11.014441592602397, 11.006623818793884, 10.998817205158629, 10.991000267878064, 10.983194504016074, 10.975378491039123, 10.967562980629749, 10.959747994522592, 10.951880288390003, 10.944066492285357, 10.936189550142963, 10.928313402465601, 10.920374520698271, 10.912500117604942, 10.904626572232518, 10.896690632892057, 10.888755777800746, 10.880811513072102, 10.8728789036451, 10.864936956421648, 10.857006674484186, 10.849004409586087, 10.841066191272299, 10.833066699209802, 10.82506866991672, 10.817061732270224, 10.809066703203458, 10.801062836005844, 10.793060541091336, 10.785070167041257, 10.77699886397654, 10.76887791783464, 10.760944078571372, 10.752878249077925, 10.744814323257327, 10.73669090277593, 10.728630933154967, 10.720562717885755, 10.712435373362185, 10.704381461214611, 10.696248293754332, 10.68812754986115, 10.679998945555672, 10.671872657338106, 10.66374870294245, 10.655627100034911, 10.647377012822927, 10.639310616916525, 10.631126054186389, 10.623004244283376, 10.614814877440116, 10.606568284744819, 10.598444269057111, 10.59025310383537, 10.582064722549852, 10.57380965859254, 10.565627002827775, 10.557377914170022, 10.549191169340713, 10.540938271691973, 10.532688503078763, 10.524500878070823, 10.516247503918416, 10.507997336750122, 10.49962316593525, 10.491438286471281, 10.483188184915655, 10.47487320088849, 10.466629820433134, 10.458253839965234, 10.450065816728491, 10.441687431817027, 10.433438499753056, 10.42512548543335, 10.41681615079163, 10.40844319537712, 10.40018937503701, 10.391881609224544, 10.383501049295452, 10.375191398855259, 10.366818798883426, 10.358440729499895, 10.350190365720255, 10.341810812085095, 10.3334354138662, 10.32506418388654, 10.31668768832435, 10.308315417258036, 10.299937967494074, 10.29156479805361, 10.283186536323601, 10.274812610315767, 10.266443032216959, 10.258003103058957, 10.249623732486596, 10.241127802248615, 10.23281328874568, 10.224317512592439, 10.2159377893319, 10.207562618574123, 10.199127391920639, 10.190623226943446, 10.182252915727865, 10.173813749080184, 10.165315275352437, 10.156940919676977, 10.148498169038554, 10.140005648878276, 10.131627603380904, 10.12312705984157, 10.114686235830591, 10.106187046868802, 10.097693157567827, 10.089312966996953, 10.080811487149731, 10.072378410480974, 10.063878623481507, 10.055375302186315, 10.046940152352295, 10.038438745149529, 10.029942862395018, 10.021443601749656 ] } ], "layout": { "height": 500, "hovermode": "closest", "legend": { "bgcolor": "#F2F2F2", "borderwidth": 0.5 }, "plot_bgcolor": "#FDFDFF", "showlegend": true, "template": { "data": { "bar": [ { "error_x": { "color": "#2a3f5f" }, "error_y": { "color": "#2a3f5f" }, "marker": { "line": { "color": "#E5ECF6", "width": 0.5 }, "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "bar" } ], "barpolar": [ { "marker": { "line": { "color": "#E5ECF6", "width": 0.5 }, "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "barpolar" } ], "carpet": [ { "aaxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "baxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "type": "carpet" } ], "choropleth": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "choropleth" } ], "contour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "contour" } ], "contourcarpet": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "contourcarpet" } ], "heatmap": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmap" } ], "heatmapgl": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmapgl" } ], "histogram": [ { "marker": { "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "histogram" } ], "histogram2d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2d" } ], "histogram2dcontour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2dcontour" } ], "mesh3d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "mesh3d" } ], "parcoords": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "parcoords" } ], "pie": [ { "automargin": true, "type": "pie" } ], "scatter": [ { "fillpattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 }, "type": "scatter" } ], "scatter3d": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatter3d" } ], "scattercarpet": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattercarpet" } ], "scattergeo": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergeo" } ], "scattergl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergl" } ], "scattermapbox": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattermapbox" } ], "scatterpolar": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolar" } ], "scatterpolargl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolargl" } ], "scatterternary": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterternary" } ], "surface": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "surface" } ], "table": [ { "cells": { "fill": { "color": "#EBF0F8" }, "line": { "color": "white" } }, "header": { "fill": { "color": "#C8D4E3" }, "line": { "color": "white" } }, "type": "table" } ] }, "layout": { "annotationdefaults": { "arrowcolor": "#2a3f5f", "arrowhead": 0, "arrowwidth": 1 }, "autotypenumbers": "strict", "coloraxis": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "colorscale": { "diverging": [ [ 0, "#8e0152" ], [ 0.1, "#c51b7d" ], [ 0.2, "#de77ae" ], [ 0.3, "#f1b6da" ], [ 0.4, "#fde0ef" ], [ 0.5, "#f7f7f7" ], [ 0.6, "#e6f5d0" ], [ 0.7, "#b8e186" ], [ 0.8, "#7fbc41" ], [ 0.9, "#4d9221" ], [ 1, "#276419" ] ], "sequential": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "sequentialminus": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ] }, "colorway": [ "#636efa", "#EF553B", "#00cc96", "#ab63fa", "#FFA15A", "#19d3f3", "#FF6692", "#B6E880", "#FF97FF", "#FECB52" ], "font": { "color": "#2a3f5f" }, "geo": { "bgcolor": "white", "lakecolor": "white", "landcolor": "#E5ECF6", "showlakes": true, "showland": true, "subunitcolor": "white" }, "hoverlabel": { "align": "left" }, "hovermode": "closest", "mapbox": { "style": "light" }, "paper_bgcolor": "white", "plot_bgcolor": "#E5ECF6", "polar": { "angularaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "radialaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "scene": { "xaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "yaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "zaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" } }, "shapedefaults": { "line": { "color": "#2a3f5f" } }, "ternary": { "aaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "baxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "caxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "title": { "x": 0.05 }, "xaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 }, "yaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 } } }, "title": { "text": "Si(111) at Cu K edge" }, "width": 650, "xaxis": { "color": "#004", "gridcolor": "#D8D8D8", "showgrid": true, "title": { "text": "Energy (eV)" }, "zerolinecolor": "#DDD" }, "yaxis": { "color": "#004", "gridcolor": "#D8D8D8", "showgrid": true, "title": { "text": "$\\theta$" }, "zerolinecolor": "#DDD" } } }, "text/html": [ "
" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# now we will read that data file we just downloaded\n", "\n", "\n", "cudat = read_ascii(cudata_file, labels='energy, mu')\n", "\n", "# This data for Cu metal at the Cu K edge measured \n", "# with a Si(111) monochromator. It is a little cryptic, \n", "# but this data file set says that hc_2d = 1977.260\n", "hc_2d = 1977.260\n", "# use the current value of HC to get the dspace they would have used\n", "dspace = PLANCK_HC / (2 * hc_2d)\n", "print(f'Using a d spacing {dspace:.6f} Ang')\n", "\n", "\n", "# now let's calculate the mono angles in degrees and plot that\n", "cudat.theta = RAD2DEG * np.arcsin(PLANCK_HC / (2.0*dspace*cudat.energy))\n", "\n", "\n", "plot(cudat.energy, cudat.theta, xlabel='Energy (eV)', ylabel=r\"$\\theta$\", \n", " color='red', title='Si(111) at Cu K edge', label='angle', linewidth=2)\n" ] }, { "cell_type": "markdown", "id": "a363cf36", "metadata": {}, "source": [ "note from the plot above that the angle is around 10 to 12 degrees, and that Energy and Angle are not perfectly linear." ] }, { "cell_type": "code", "execution_count": 4, "id": "e97bd9f5", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "E0 = 8977.580 eV\n" ] }, { "data": { "application/vnd.plotly.v1+json": { "config": { "plotlyServerURL": "https://plot.ly" }, "data": [ { "line": { "color": "#1f77b4", "width": 2 }, "name": "norm", "type": "scatter", "uid": "00526c9a-bd94-47c9-8816-c25d116a5a64", "x": [ 8786.204, 8796.258, 8806.253, 8816.27, 8826.27, 8836.292, 8846.337, 8856.322, 8866.372, 8876.361, 8886.416, 8896.409, 8906.425, 8916.464, 8926.483, 8936.483, 8936.998, 8937.471, 8937.944, 8938.417, 8938.976, 8939.406, 8939.878, 8940.352, 8940.824, 8941.384, 8941.814, 8942.287, 8942.761, 8943.319, 8943.793, 8944.224, 8944.697, 8945.257, 8945.73, 8946.204, 8946.635, 8947.108, 8947.668, 8948.143, 8948.616, 8949.048, 8949.564, 8950.038, 8950.513, 8950.987, 8951.461, 8951.979, 8952.453, 8952.927, 8953.401, 8953.876, 8954.394, 8954.868, 8955.343, 8955.817, 8956.336, 8956.811, 8957.285, 8957.761, 8958.235, 8958.71, 8959.229, 8959.703, 8960.179, 8960.654, 8961.13, 8961.647, 8962.123, 8962.599, 8963.074, 8963.55, 8964.068, 8964.544, 8965.02, 8965.496, 8965.972, 8966.447, 8966.967, 8967.442, 8967.919, 8968.395, 8968.871, 8969.347, 8969.909, 8970.386, 8970.862, 8971.339, 8971.772, 8972.248, 8972.812, 8973.288, 8973.765, 8974.198, 8974.675, 8975.151, 8975.628, 8976.192, 8976.626, 8977.103, 8977.58, 8978.057, 8978.534, 8979.055, 8979.531, 8980.009, 8980.486, 8980.964, 8981.441, 8981.919, 8982.396, 8982.918, 8983.396, 8983.873, 8984.352, 8984.829, 8985.307, 8985.785, 8986.307, 8986.784, 8987.263, 8987.741, 8988.219, 8988.697, 8989.176, 8989.653, 8990.132, 8990.698, 8991.133, 8991.611, 8992.09, 8992.568, 8993.048, 8993.526, 8994.006, 8994.484, 8995.007, 8995.485, 8995.965, 8996.444, 8996.924, 8997.402, 8997.882, 8998.361, 8998.84, 8999.319, 8999.799, 9000.366, 9000.803, 9001.281, 9001.762, 9002.241, 9002.721, 9003.201, 9003.681, 9004.161, 9004.642, 9005.121, 9005.602, 9006.125, 9006.605, 9007.086, 9007.566, 9008.047, 9008.526, 9009.007, 9009.487, 9009.969, 9011.062, 9012.153, 9013.246, 9014.384, 9015.565, 9016.746, 9017.928, 9019.154, 9020.38, 9021.606, 9022.877, 9024.191, 9025.464, 9026.822, 9028.14, 9029.499, 9030.903, 9032.309, 9033.759, 9035.164, 9036.615, 9038.109, 9039.604, 9041.145, 9042.642, 9044.271, 9045.768, 9047.397, 9049.029, 9050.616, 9052.336, 9053.969, 9055.646, 9057.324, 9059.047, 9060.814, 9062.627, 9064.353, 9066.166, 9067.981, 9069.842, 9071.702, 9073.563, 9075.47, 9077.377, 9079.329, 9081.283, 9083.237, 9085.281, 9087.237, 9089.283, 9091.33, 9093.378, 9095.472, 9097.61, 9099.705, 9101.847, 9104.033, 9106.221, 9108.409, 9110.644, 9112.879, 9115.16, 9117.398, 9119.728, 9122.057, 9124.389, 9126.721, 9129.099, 9131.478, 9133.904, 9136.331, 9138.805, 9141.279, 9143.756, 9146.233, 9148.803, 9151.328, 9153.899, 9156.474, 9159.094, 9161.716, 9164.339, 9167.009, 9169.681, 9172.399, 9175.074, 9177.842, 9180.564, 9183.381, 9186.153, 9188.973, 9191.794, 9194.663, 9197.533, 9200.406, 9203.326, 9206.247, 9209.217, 9212.143, 9215.115, 9218.183, 9221.16, 9224.185, 9227.258, 9230.333, 9233.502, 9236.581, 9239.755, 9242.838, 9246.016, 9249.196, 9252.425, 9255.655, 9258.889, 9262.171, 9265.454, 9268.787, 9272.122, 9275.459, 9278.8, 9282.188, 9285.58, 9289.021, 9292.465, 9295.91, 9299.404, 9302.949, 9306.403, 9309.953, 9313.506, 9317.061, 9320.666, 9324.32, 9327.932, 9331.592, 9335.256, 9338.922, 9342.639, 9346.358, 9350.128, 9353.9, 9357.677, 9361.455, 9365.285, 9369.165, 9373.048, 9376.84, 9380.777, 9384.671, 9388.614, 9392.563, 9396.561, 9400.515, 9404.567, 9408.576, 9412.637, 9416.7, 9420.768, 9424.935, 9429.01, 9433.184, 9437.313, 9441.495, 9445.729, 9449.966, 9454.208, 9458.453, 9462.702, 9467.004, 9471.358, 9475.668, 9480.029, 9484.396, 9488.814, 9493.189, 9497.665, 9502.097, 9506.581, 9511.118, 9515.563, 9520.109, 9524.611, 9529.265, 9533.775, 9538.438, 9543.007, 9547.679, 9552.306, 9556.986, 9561.673, 9566.462, 9571.157, 9575.956, 9580.661, 9585.47, 9590.283, 9595.102, 9599.976, 9604.854, 9609.737, 9614.625, 9619.568, 9624.518, 9629.471, 9634.431, 9639.445, 9644.466, 9649.541, 9654.571, 9659.708, 9664.75, 9669.847, 9675.001, 9680.159, 9685.324, 9690.545, 9695.722, 9700.953, 9706.243, 9711.486, 9716.788, 9722.094, 9727.407, 9732.776, 9738.204, 9743.535, 9748.923, 9754.317, 9759.769, 9765.227, 9770.742, 9776.212, 9781.793, 9787.275, 9792.816, 9798.416, 9804.022, 9809.584, 9815.254, 9820.88, 9826.564, 9832.308, 9837.953, 9843.659, 9849.422, 9855.192, 9860.971, 9866.755, 9872.6, 9878.45, 9884.31, 9890.175, 9896.1, 9901.98, 9907.92, 9913.867, 9919.876, 9925.838, 9931.859, 9937.89, 9943.98, 9950.025, 9956.132, 9962.3, 9968.367, 9974.497, 9980.687, 9986.832, 9993.038, 9999.253, 10005.53, 10011.76, 10018.05, 10024.35, 10030.66, 10037.03, 10043.36, 10049.74, 10056.14, 10062.54, 10068.96, 10075.43, 10081.92, 10088.46, 10094.91, 10101.42, 10107.94, 10114.52, 10121.11, 10127.71, 10134.32, 10141.05, 10147.68, 10154.31, 10161.01, 10167.78, 10174.49, 10181.28, 10188.07, 10194.82, 10201.63, 10208.45, 10215.28, 10222.17, 10229.13, 10235.99, 10242.86, 10249.85, 10256.79, 10263.8, 10270.82, 10277.85, 10284.89, 10291.88, 10299.05, 10306.12, 10313.2, 10320.41, 10327.5, 10334.73, 10341.85, 10349.09, 10356.35, 10363.61, 10370.89, 10378.17, 10385.47, 10392.78, 10400.1, 10407.48, 10414.82, 10422.23, 10429.65, 10437.14, 10444.58, 10452.03, 10459.55, 10467.08, 10474.63, 10482.18, 10489.75, 10497.32, 10504.97, 10512.57, 10520.24, 10527.92, 10535.62, 10543.32, 10551.04, 10558.77, 10566.5, 10574.32, 10582.2, 10589.91, 10597.76, 10605.62, 10613.55, 10621.43, 10629.33, 10637.3, 10645.21, 10653.21, 10661.21, 10669.23, 10677.26, 10685.3, 10693.35, 10701.54, 10709.56, 10717.71, 10725.81, 10733.99, 10742.24, 10750.38, 10758.6, 10766.83, 10775.14, 10783.39, 10791.72, 10800, 10808.36, 10816.73, 10825.05, 10833.45, 10841.86, 10850.41, 10858.78, 10867.23, 10875.76, 10884.23, 10892.85, 10901.29, 10909.94, 10918.47, 10927.08, 10935.7, 10944.4, 10952.99, 10961.65, 10970.4, 10979.09, 10987.86, 10996.65, 11005.32, 11014.14, 11022.97, 11031.81, 11040.67, 11049.54, 11058.43, 11067.33, 11076.25, 11085.18, 11094.12, 11103.15, 11112.13, 11121.25, 11130.19, 11139.34, 11148.38, 11157.43, 11166.56, 11175.78, 11184.87, 11194.05, 11203.31, 11212.45, 11221.68, 11230.98, 11240.17, 11249.51, 11258.8, 11268.17, 11277.55, 11286.82, 11296.24, 11305.6, 11315.05, 11324.52, 11333.93, 11343.43, 11352.94, 11362.47 ], "y": [ 0.006099431461700133, 9.117887986916555e-05, 5.4670997746302065e-05, -2.7191318620089223e-05, -6.938989746250113e-05, -0.00010056449982715447, -7.513226728440632e-05, -1.3385541299481825e-05, 5.272897741641901e-05, 8.708466946053892e-05, 0.00017065148104311614, 0.00041148213917538364, 0.0005611003115116054, 0.0009345899709346704, 0.001294474862171137, 0.0021675119184885456, 0.0022607822078094674, 0.0023875422375935563, 0.0023762273633420768, 0.0023721337978712797, 0.0025079414162668275, 0.002596333204568099, 0.0026458845236113394, 0.002719086672601256, 0.002780035478997492, 0.0028462330302125366, 0.0028928195309341903, 0.002940493179011494, 0.0029652777576828714, 0.0030315329927198486, 0.003082027170724749, 0.003230835951758087, 0.003293082722977108, 0.0033462297161551994, 0.003441320391768271, 0.0035178286821276475, 0.003617741372381377, 0.0036910811179321993, 0.0037889915251784114, 0.003887461163920137, 0.0039618884559740304, 0.004083327476008513, 0.0041815063845434056, 0.004229346949080708, 0.004403031303977713, 0.004464922969334836, 0.004629870541327335, 0.004732298450192962, 0.004796756725297867, 0.004937082244459865, 0.005072579057148074, 0.005231211768441568, 0.005365396035198512, 0.005503415955773532, 0.005692760980315324, 0.00584918218772303, 0.006005023289841032, 0.006139512468765868, 0.006347526882281607, 0.0064916228813130285, 0.006707163116630306, 0.006988601579055301, 0.007225965167046599, 0.007393348843206092, 0.007674388697909812, 0.007874130667021008, 0.008181967667561034, 0.008389155146587373, 0.008714218883736713, 0.009033975393950906, 0.009430780789426454, 0.009781902140790373, 0.01024961516912145, 0.010635494506724855, 0.011077012723427917, 0.011505654389533273, 0.012058319858854423, 0.012576538945937448, 0.013250784624178266, 0.013968024721330552, 0.014837420143237271, 0.015396132371115685, 0.016182663840465787, 0.017379504854510168, 0.018681054047579807, 0.019737420464301995, 0.021237686952714403, 0.02338142586510873, 0.02435640779166551, 0.026455259875391095, 0.029764405731245495, 0.03331113022549566, 0.03785164762330877, 0.041459318622740333, 0.05100611708740155, 0.05791562237299598, 0.06834898874699495, 0.09645415137195536, 0.1310924751406965, 0.19162954393169235, 0.29135526341584717, 0.3803821354209815, 0.44203154677536083, 0.4717576761114283, 0.48211339875608555, 0.47868162767092126, 0.4696792329784953, 0.45905921452635834, 0.45370314574956777, 0.460077563973669, 0.4715597601952697, 0.49142640335243826, 0.5145646983941182, 0.533314315368118, 0.5491166671367219, 0.5660078846462031, 0.5931265998892751, 0.6100014971102291, 0.6415096730553159, 0.6720896208901742, 0.7170621449652635, 0.7511228870470033, 0.7992947625159293, 0.8462098692458588, 0.8807676129885406, 0.9149103631675959, 0.957964020192593, 1.0050294543708573, 1.023806775337923, 1.043785095260159, 1.0472594548880798, 1.0409952776902094, 1.0314432157221682, 1.0148108051822236, 0.9982605989759534, 0.9809573897528895, 0.9620462501947948, 0.9475567412845837, 0.939415879658981, 0.9337626789955784, 0.9314536767382668, 0.9352097180014745, 0.9441046198192283, 0.9585810348912375, 0.9864468874909668, 1.0012952434670397, 1.0253875201192892, 1.0474691764091684, 1.0639987420816401, 1.0716903966303117, 1.0704552306991924, 1.0622623914417215, 1.0439051179658794, 1.033929782898826, 1.0160123132288432, 1.006598152157396, 0.982101671580503, 0.9711447241307211, 0.9575593801102559, 0.9439850487522351, 0.935379587260503, 0.926429366436303, 0.9145634405282396, 0.9063584065679724, 0.8975626315291759, 0.8916597160066978, 0.8881291815170896, 0.8873744865400649, 0.8921382376259682, 0.898582479750442, 0.9095865726712046, 0.9219882787729393, 0.9452540520192998, 0.9673239591308375, 0.9847572133900242, 1.0101331602861956, 1.034947063149558, 1.0543436793716165, 1.069669289938395, 1.07162799654185, 1.0586169162601677, 1.0324870298152853, 1.0070408412321845, 0.9768354759305574, 0.9539400616194232, 0.9388914986555852, 0.9286304384139012, 0.926698780915085, 0.9311711410786876, 0.9416993077873551, 0.9514925034782377, 0.9544538408700594, 0.9505117450632707, 0.9462159131614261, 0.9493802812245763, 0.9614241830040963, 0.9738532724820766, 0.9770239626545879, 0.9749541281839764, 0.976383112710594, 0.9860384220022077, 0.9997852916568509, 1.0123174877919885, 1.0205431227902535, 1.0287801771056422, 1.0385941382379584, 1.0534782767803708, 1.0688987790761573, 1.080862732399942, 1.0816425820859978, 1.0709190102990507, 1.0485318286090124, 1.0183369330638294, 0.9934030508972742, 0.9741765985970925, 0.9573353532875375, 0.9483594761374367, 0.9468426169119785, 0.9376883679534058, 0.9146467955870067, 0.8935177754100005, 0.8783586243313712, 0.8752596169505101, 0.8846152795261619, 0.9008402050672836, 0.9219339818089923, 0.9453459584947743, 0.967569029343354, 0.982210045224044, 0.9856441118994594, 0.9789802006176652, 0.9708634146105803, 0.9686652743404159, 0.98666875247515, 1.02993352200194, 1.084369569208702, 1.1122391392709219, 1.1138429316268879, 1.1015133974367877, 1.0893041387989504, 1.0751454460421443, 1.048832768585153, 1.0134927813514836, 0.9702354555751549, 0.9232465333040473, 0.8770964284809694, 0.8375387043746495, 0.8131230793765442, 0.8105940086594493, 0.8321323278742748, 0.8699657842728565, 0.9144644799129241, 0.9578799650539309, 0.9866828491940066, 1.0023647896861978, 0.999979530698181, 0.9836309358139109, 0.9683620359826701, 0.961117155836223, 0.9584139645383036, 0.9467626510310942, 0.9295096976471825, 0.9262762991439139, 0.9469869557081324, 0.9898368981787231, 1.0325974372357776, 1.04757973356158, 1.049494818127731, 1.0432286166937583, 1.025126618460171, 0.9927860191405699, 0.9524872841813219, 0.9147826597928286, 0.8862377728901889, 0.8691971563824952, 0.8618789737636843, 0.8652035001090724, 0.879343016651351, 0.8898809930786616, 0.8901420758881987, 0.8929287467478129, 0.8970219888852865, 0.8990062324800004, 0.9062812669545872, 0.9209622685150891, 0.9344722804706077, 0.9437486772168484, 0.9487027149591579, 0.9498185552187592, 0.9497105723862685, 0.9545448075185998, 0.9617684698788191, 0.9710090747107682, 0.9789383088782843, 0.9702975545756256, 0.9499977226102866, 0.9313857555691152, 0.9110443438724538, 0.8924462479473396, 0.8843812432322101, 0.8877159814322754, 0.8970755973108331, 0.9045270844069161, 0.8993174497652923, 0.882638956356522, 0.8611334823663758, 0.8515325832871878, 0.8626938028076607, 0.8745021858760732, 0.8756438986599964, 0.8812186330450035, 0.8893014600845084, 0.8933365053270356, 0.8976108912698163, 0.9022439287763816, 0.9063732074940425, 0.9103585544089995, 0.9133756373204529, 0.915218040059793, 0.9167389215069297, 0.9178664161897893, 0.916793463772099, 0.9104334703959749, 0.8969610102923625, 0.8794398983545699, 0.8652462295157394, 0.8606248920331012, 0.8617483607193223, 0.8667941744659874, 0.8742751728632395, 0.8788629641284397, 0.8780955180996398, 0.8726582136740815, 0.8608232869163714, 0.8491764702487967, 0.8478998232340591, 0.8529805096939942, 0.8526102547676554, 0.8476306090301311, 0.8457142412671218, 0.8464765739198666, 0.8509578380487717, 0.8583057306291684, 0.8648946031532602, 0.8726357762909063, 0.8819550738984164, 0.8899563314317614, 0.8931957170432046, 0.8906252115680379, 0.8807428497820006, 0.8626093671706966, 0.8440054216595896, 0.8367682282444255, 0.8355007977206106, 0.8358588282246266, 0.8359000873438316, 0.837341337906672, 0.8401344679695073, 0.8415113453391251, 0.8394873940413847, 0.8381775398018678, 0.8377687416607446, 0.8352303239815054, 0.8286650716857552, 0.8259897273584161, 0.8267535975466362, 0.8271092477922529, 0.8300862917650494, 0.8328847251471315, 0.8334195037827323, 0.8341529134244667, 0.8372856185245032, 0.8415756693402106, 0.8416779639541678, 0.8389183232961801, 0.8357094762955424, 0.8325824648975845, 0.830888941375558, 0.8280094509395097, 0.8226977926191602, 0.8190046660194084, 0.8173437596480766, 0.8166998984025333, 0.8176849146724874, 0.8162390849223692, 0.8118250002405033, 0.8109559342458291, 0.811605752348145, 0.8107205055734571, 0.8078578983975526, 0.805905756790131, 0.8048573166899318, 0.8021697644239261, 0.8025410621033174, 0.8043396911106152, 0.80617772013403, 0.8099854708379974, 0.8111178808860322, 0.8109521528820676, 0.810604835444787, 0.8078302157973883, 0.8030065239671748, 0.7998768440493922, 0.8002031860599436, 0.8015459656830464, 0.8022801117371239, 0.8025895566209336, 0.8009976394099932, 0.796909897468257, 0.7932419722184467, 0.7900634524547986, 0.7880912887515856, 0.7876145039589342, 0.7868782614033442, 0.785246483420226, 0.7843909323122046, 0.783347635222842, 0.7825851451093653, 0.7829225640742327, 0.7823746492757602, 0.7821314295525564, 0.7835617359481047, 0.7848651184969174, 0.7849939603821315, 0.7844469512334469, 0.7835397420607317, 0.7819705503437115, 0.7808474868530114, 0.7799414267345256, 0.7783073796375992, 0.7757655503047028, 0.7742460100625337, 0.7730368580611338, 0.7730601774182054, 0.7731459102313756, 0.7732016363208223, 0.7710036326220469, 0.7695471205978753, 0.7695911414198725, 0.7682587587653312, 0.7663118105940208, 0.7646010870829482, 0.7634310207834172, 0.7627968264908127, 0.7616976068732406, 0.7603241741575081, 0.7597978031279683, 0.7597206644939014, 0.7610586413696856, 0.762547862899089, 0.7621283032757984, 0.7612322564645312, 0.7598995314471827, 0.7592987108867857, 0.7585560838221451, 0.7574266543111912, 0.7559113728606949, 0.7548388694364594, 0.7534543867653245, 0.7528541717224281, 0.7522174342365513, 0.7512690624930461, 0.7506617453329167, 0.7510166279686973, 0.7501812383614269, 0.7488438277285572, 0.7474044904831457, 0.7460625974044128, 0.7451586166777585, 0.7443407791939128, 0.7436646868642668, 0.743368294863701, 0.7434273226201144, 0.7437307694803829, 0.7436995260971572, 0.7432288430168671, 0.7420460536058634, 0.7410941786025518, 0.7409979238378449, 0.7403682653582049, 0.739818579391388, 0.7395681153124208, 0.7387985935725994, 0.7380338857166782, 0.7377131709134914, 0.7370039135895862, 0.7360797052863829, 0.7349492790466655, 0.7338078546418284, 0.7328877646756945, 0.7321834630553017, 0.7317587365758971, 0.7314967070533577, 0.7308734702721713, 0.7301762803287742, 0.7293149465530941, 0.7285736070015517, 0.7279178952046477, 0.7277037368477696, 0.7278148210252474, 0.7276056701621647, 0.7266999836086818, 0.7256258835064731, 0.7253177706914882, 0.7246754926809403, 0.7240709999386562, 0.7240128196497297, 0.7234851124704469, 0.7230692050716838, 0.7227023129552484, 0.7221548881192963, 0.7216925946693372, 0.7210405526396667, 0.7202722992934193, 0.7195700978647549, 0.7189303972950564, 0.7183602997015571, 0.7177743996188006, 0.717000500590499, 0.7162466124178547, 0.7156439075487863, 0.7151325661457415, 0.714687247980993, 0.7143844932173643, 0.7141947723802528, 0.7136195637123931, 0.7130842962862802, 0.7124497027994698, 0.7121513173413774, 0.7119320957089497, 0.711502121903862, 0.7110660578383569, 0.7104680250335754, 0.7099604051880466, 0.709472290269767, 0.709124180431282, 0.7088472909433123, 0.7085246249131186, 0.708015284112709, 0.7077625238174652, 0.7073395389892835, 0.7068769061183443, 0.7063645910975325, 0.7058712846353243, 0.7052922085988642, 0.7046420922705505, 0.7043256256496293, 0.7038480698863968, 0.7035966433262432, 0.703508651500307, 0.7033583811045853, 0.7030000913479006, 0.7024968081371503, 0.7023658299392591, 0.7019593622666904, 0.7016389995978827, 0.7013480859341081, 0.7009082636488576, 0.7007172624437563, 0.7002376502961347, 0.699992484495298, 0.6996871020773557, 0.6994656591299294, 0.6991375043653474, 0.6988701812951335, 0.6986666736994509, 0.6983012491547751, 0.6980262471291135, 0.697650552636907, 0.6972539676920789, 0.6970247326488789, 0.6966422292384273, 0.6964924416601179, 0.6961910739986797, 0.6959139964691048, 0.6957529979962559, 0.6955144243927722, 0.6953545182463075, 0.6952958296143197, 0.6948686093172707, 0.6948562549461965, 0.6946215351600799, 0.6944132805988328, 0.6947329354498389, 0.6942020484658425, 0.6939759868212309, 0.6938412504035937, 0.6937141602706455, 0.6935769746805385, 0.6934062217483287, 0.6932797084535949, 0.693152556970603, 0.6930733093240299, 0.6929081598433539, 0.6927189520797588, 0.6928429385895357, 0.6925256130679929, 0.6924422662038753, 0.6923897156311157, 0.6923517620809029, 0.6922710962381989, 0.6921676771683798, 0.6921446656299842, 0.6922438042892081, 0.6921404656889575, 0.6919601192365243, 0.6920292702915485, 0.6919782143930162, 0.6921381329563252, 0.6920137812033639, 0.6922082271753086, 0.692031025442177 ] } ], "layout": { "height": 500, "hovermode": "closest", "legend": { "bgcolor": "#F2F2F2", "borderwidth": 0.5 }, "plot_bgcolor": "#FDFDFF", "showlegend": true, "template": { "data": { "bar": [ { "error_x": { "color": "#2a3f5f" }, "error_y": { "color": "#2a3f5f" }, "marker": { "line": { "color": "#E5ECF6", "width": 0.5 }, "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "bar" } ], "barpolar": [ { "marker": { "line": { "color": "#E5ECF6", "width": 0.5 }, "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "barpolar" } ], "carpet": [ { "aaxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "baxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "type": "carpet" } ], "choropleth": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "choropleth" } ], "contour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "contour" } ], "contourcarpet": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "contourcarpet" } ], "heatmap": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmap" } ], "heatmapgl": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmapgl" } ], "histogram": [ { "marker": { "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "histogram" } ], "histogram2d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2d" } ], "histogram2dcontour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2dcontour" } ], "mesh3d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "mesh3d" } ], "parcoords": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "parcoords" } ], "pie": [ { "automargin": true, "type": "pie" } ], "scatter": [ { "fillpattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 }, "type": "scatter" } ], "scatter3d": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatter3d" } ], "scattercarpet": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattercarpet" } ], "scattergeo": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergeo" } ], "scattergl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergl" } ], "scattermapbox": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattermapbox" } ], "scatterpolar": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolar" } ], "scatterpolargl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolargl" } ], "scatterternary": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterternary" } ], "surface": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "surface" } ], "table": [ { "cells": { "fill": { "color": "#EBF0F8" }, "line": { "color": "white" } }, "header": { "fill": { "color": "#C8D4E3" }, "line": { "color": "white" } }, "type": "table" } ] }, "layout": { "annotationdefaults": { "arrowcolor": "#2a3f5f", "arrowhead": 0, "arrowwidth": 1 }, "autotypenumbers": "strict", "coloraxis": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "colorscale": { "diverging": [ [ 0, "#8e0152" ], [ 0.1, "#c51b7d" ], [ 0.2, "#de77ae" ], [ 0.3, "#f1b6da" ], [ 0.4, "#fde0ef" ], [ 0.5, "#f7f7f7" ], [ 0.6, "#e6f5d0" ], [ 0.7, "#b8e186" ], [ 0.8, "#7fbc41" ], [ 0.9, "#4d9221" ], [ 1, "#276419" ] ], "sequential": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "sequentialminus": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ] }, "colorway": [ "#636efa", "#EF553B", "#00cc96", "#ab63fa", "#FFA15A", "#19d3f3", "#FF6692", "#B6E880", "#FF97FF", "#FECB52" ], "font": { "color": "#2a3f5f" }, "geo": { "bgcolor": "white", "lakecolor": "white", "landcolor": "#E5ECF6", "showlakes": true, "showland": true, "subunitcolor": "white" }, "hoverlabel": { "align": "left" }, "hovermode": "closest", "mapbox": { "style": "light" }, "paper_bgcolor": "white", "plot_bgcolor": "#E5ECF6", "polar": { "angularaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "radialaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "scene": { "xaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "yaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "zaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" } }, "shapedefaults": { "line": { "color": "#2a3f5f" } }, "ternary": { "aaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "baxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "caxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "title": { "x": 0.05 }, "xaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 }, "yaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 } } }, "title": {}, "width": 650, "xaxis": { "color": "#004", "gridcolor": "#D8D8D8", "range": [ 8900, 9500 ], "showgrid": true, "title": { "text": "Energy (eV)" }, "zerolinecolor": "#DDD" }, "yaxis": { "color": "#004", "gridcolor": "#D8D8D8", "showgrid": true, "title": { "text": "$\\mu(E)$" }, "zerolinecolor": "#DDD" } } }, "text/html": [ "
" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# let's plot the normalized XAFS data around the edge:\n", "\n", "autobk(cudat.energy, cudat.mu, group=cudat, rbkg=1.0, kw=2)\n", "print(f'E0 = {cudat.e0:.3f} eV')\n", "plot(cudat.energy, cudat.norm, xmin=8900, xmax=9500, label='norm',\n", " xlabel='Energy (eV)', ylabel='$\\mu(E)$')\n", " " ] }, { "cell_type": "code", "execution_count": 5, "id": "b8f6a375", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Result: E0 = 8977.580 eV, R= 2.5494 +/- 0.0017 Ang\n", "=================== FEFFIT RESULTS ====================\n", "[[Statistics]]\n", " n_function_calls = 31\n", " n_variables = 4\n", " n_data_points = 104\n", " n_independent = 13.7323954\n", " chi_square = 80.4996792\n", " reduced chi_square = 8.27131200\n", " r-factor = 0.00122781\n", " Akaike info crit = 32.2856799\n", " Bayesian info crit = 34.7647106\n", " \n", "[[Variables]]\n", " amp = 0.9216005 +/- 0.0271218 (init= 1.0000000)\n", " del_e0 = 5.6772818 +/- 0.3759213 (init= 0.0000000)\n", " del_r = 0.0015526 +/- 0.0016816 (init= 0.0000000)\n", " sig2 = 0.0035496 +/- 1.8541e-4 (init= 0.0000000)\n", " \n", "[[Correlations]] (unreported correlations are < 0.100)\n", " amp, sig2 = 0.924\n", " del_e0, del_r = 0.914\n", " del_r, sig2 = 0.135\n", " amp, del_r = 0.131\n", " \n", "[[Dataset]]\n", " unique_id = 'dzhvfzs6'\n", " fit space = 'r'\n", " r-range = 1.400, 3.000\n", " k-range = 2.500, 15.000\n", " k window, dk = 'kaiser', 4.000\n", " paths used in fit = ['../feffit/feff0001.dat']\n", " k-weight = 2\n", " epsilon_k = Array(mean=0.0013512, std=9.2103e-4)\n", " epsilon_r = 0.0332145\n", " n_independent = 13.732\n", " \n", "[[Paths]]\n", " = Path 'p3tf6mewg' = Cu K Edge\n", " feffdat file = ../feffit/feff0001.dat, from feff run 'feffit'\n", " geometry atom x y z ipot\n", " Cu 0.0000, 0.0000, 0.0000 0 (absorber)\n", " Cu 0.0000, -1.8016, 1.8016 1\n", " reff = 2.5478000\n", " degen = 12.000000\n", " n*s02 = 0.9216005 +/- 0.0271218 := 'amp'\n", " e0 = 5.6772818 +/- 0.3759213 := 'del_e0'\n", " r = 2.5493526 +/- 0.0016816 := 'reff + del_r'\n", " deltar = 0.0015526 +/- 0.0016816 := 'del_r'\n", " sigma2 = 0.0035496 +/- 1.8541e-4 := 'sig2'\n", "\n", "=======================================================\n" ] }, { "data": { "application/vnd.plotly.v1+json": { "config": { "plotlyServerURL": "https://plot.ly" }, "data": [ { "line": { "color": "#1f77b4", "width": 3 }, "name": "data", "type": "scatter", "uid": "57c08797-0fe0-437d-b16c-a71e32aa735b", "x": [ 0, 0.030679615757712823, 0.06135923151542565, 0.09203884727313846, 0.1227184630308513, 0.15339807878856412, 0.18407769454627693, 0.21475731030398976, 0.2454369260617026, 0.2761165418194154, 0.30679615757712825, 0.3374757733348411, 0.36815538909255385, 0.3988350048502667, 0.4295146206079795, 0.46019423636569234, 0.4908738521234052, 0.521553467881118, 0.5522330836388308, 0.5829126993965437, 0.6135923151542565, 0.6442719309119693, 0.6749515466696822, 0.705631162427395, 0.7363107781851077, 0.7669903939428205, 0.7976700097005334, 0.8283496254582462, 0.859029241215959, 0.8897088569736719, 0.9203884727313847, 0.9510680884890975, 0.9817477042468103, 1.012427320004523, 1.043106935762236, 1.0737865515199487, 1.1044661672776617, 1.1351457830353744, 1.1658253987930873, 1.1965050145508, 1.227184630308513, 1.2578642460662257, 1.2885438618239387, 1.3192234775816514, 1.3499030933393643, 1.380582709097077, 1.41126232485479, 1.4419419406125027, 1.4726215563702154, 1.5033011721279284, 1.533980787885641, 1.564660403643354, 1.5953400194010667, 1.6260196351587797, 1.6566992509164924, 1.6873788666742053, 1.718058482431918, 1.748738098189631, 1.7794177139473437, 1.8100973297050567, 1.8407769454627694, 1.8714565612204823, 1.902136176978195, 1.932815792735908, 1.9634954084936207, 1.9941750242513334, 2.024854640009046, 2.055534255766759, 2.086213871524472, 2.116893487282185, 2.1475731030398975, 2.1782527187976104, 2.2089323345553233, 2.2396119503130363, 2.270291566070749, 2.3009711818284617, 2.3316507975861747, 2.3623304133438876, 2.3930100291016, 2.423689644859313, 2.454369260617026, 2.4850488763747385, 2.5157284921324514, 2.5464081078901644, 2.5770877236478773, 2.60776733940559, 2.6384469551633027, 2.6691265709210157, 2.6998061866787286, 2.730485802436441, 2.761165418194154, 2.791845033951867, 2.82252464970958, 2.8532042654672924, 2.8838838812250054, 2.9145634969827183, 2.945243112740431, 2.9759227284981438, 3.0066023442558567, 3.0372819600135696, 3.067961575771282, 3.098641191528995, 3.129320807286708, 3.160000423044421, 3.1906800388021335, 3.2213596545598464, 3.2520392703175593, 3.2827188860752723, 3.313398501832985, 3.3440781175906977, 3.3747577333484107, 3.4054373491061236, 3.436116964863836, 3.466796580621549, 3.497476196379262, 3.5281558121369745, 3.5588354278946874, 3.5895150436524004, 3.6201946594101133, 3.650874275167826, 3.6815538909255388, 3.7122335066832517, 3.7429131224409646, 3.773592738198677, 3.80427235395639, 3.834951969714103, 3.865631585471816, 3.8963112012295285, 3.9269908169872414, 3.9576704327449543, 3.988350048502667, 4.01902966426038, 4.049709280018092, 4.080388895775806, 4.111068511533518, 4.1417481272912315, 4.172427743048944, 4.2031073588066565, 4.23378697456437, 4.264466590322082, 4.295146206079795, 4.325825821837508, 4.356505437595221, 4.387185053352934, 4.417864669110647, 4.448544284868359, 4.479223900626073, 4.509903516383785, 4.540583132141498, 4.571262747899211, 4.601942363656923, 4.632621979414636, 4.663301595172349, 4.693981210930062, 4.724660826687775, 4.755340442445488, 4.7860200582032, 4.816699673960914, 4.847379289718626, 4.878058905476339, 4.908738521234052, 4.9394181369917645, 4.970097752749477, 5.00077736850719, 5.031456984264903, 5.062136600022616, 5.092816215780329, 5.123495831538041, 5.154175447295755, 5.184855063053467, 5.21553467881118, 5.246214294568893, 5.2768939103266055, 5.307573526084319, 5.338253141842031, 5.368932757599744, 5.399612373357457, 5.43029198911517, 5.460971604872882, 5.491651220630596, 5.522330836388308, 5.553010452146021, 5.583690067903734, 5.6143696836614465, 5.64504929941916, 5.675728915176872, 5.706408530934585, 5.737088146692298, 5.767767762450011, 5.798447378207723, 5.829126993965437, 5.859806609723149, 5.890486225480862, 5.921165841238575, 5.9518454569962875, 5.982525072754001, 6.013204688511713, 6.043884304269426, 6.074563920027139, 6.105243535784852, 6.135923151542564, 6.166602767300278, 6.19728238305799, 6.227961998815704, 6.258641614573416, 6.2893212303311286, 6.320000846088842, 6.350680461846554, 6.381360077604267, 6.41203969336198, 6.442719309119693, 6.473398924877405, 6.504078540635119, 6.534758156392831, 6.565437772150545, 6.596117387908257, 6.62679700366597, 6.657476619423683, 6.6881562351813955, 6.718835850939108, 6.749515466696821, 6.780195082454534, 6.810874698212247, 6.84155431396996, 6.872233929727672, 6.902913545485386, 6.933593161243098, 6.964272777000811, 6.994952392758524, 7.0256320085162365, 7.056311624273949, 7.086991240031662, 7.117670855789375, 7.148350471547088, 7.179030087304801, 7.209709703062513, 7.240389318820227, 7.271068934577939, 7.301748550335652, 7.332428166093365, 7.3631077818510775, 7.39378739760879, 7.424467013366503, 7.455146629124216, 7.485826244881929, 7.516505860639642, 7.547185476397354, 7.577865092155068, 7.60854470791278, 7.639224323670493, 7.669903939428206, 7.7005835551859185, 7.731263170943632, 7.761942786701344, 7.792622402459057, 7.82330201821677, 7.853981633974483, 7.884661249732195, 7.915340865489909, 7.946020481247621, 7.976700097005334, 8.007379712763047, 8.03805932852076, 8.068738944278472, 8.099418560036185, 8.130098175793899, 8.160777791551611, 8.191457407309324, 8.222137023067036, 8.252816638824749, 8.283496254582463, 8.314175870340176, 8.344855486097888, 8.3755351018556, 8.406214717613313, 8.436894333371026, 8.46757394912874, 8.498253564886452, 8.528933180644165, 8.559612796401877, 8.59029241215959, 8.620972027917304, 8.651651643675017, 8.682331259432729, 8.713010875190442, 8.743690490948154, 8.774370106705868, 8.80504972246358, 8.835729338221293, 8.866408953979006, 8.897088569736718, 8.92776818549443, 8.958447801252145, 8.989127417009858, 9.01980703276757, 9.050486648525283, 9.081166264282995, 9.11184588004071, 9.142525495798422, 9.173205111556134, 9.203884727313847, 9.23456434307156, 9.265243958829272, 9.295923574586986, 9.326603190344699, 9.357282806102411, 9.387962421860124, 9.418642037617836, 9.44932165337555, 9.480001269133263, 9.510680884890975, 9.541360500648688, 9.5720401164064, 9.602719732164113, 9.633399347921827, 9.66407896367954, 9.694758579437252, 9.725438195194965, 9.756117810952677, 9.786797426710391, 9.817477042468104, 9.848156658225816, 9.878836273983529, 9.909515889741241, 9.940195505498954, 9.970875121256668 ], "y": [ 0.0010500764195195702, 0.0010266682809029874, 0.0009793994731410354, 0.0009510150764243139, 0.0009477989665394854, 0.0009307792280536157, 0.0008793399999442237, 0.000847483772160443, 0.0009229568276458026, 0.0010831189300459852, 0.0012043858567844345, 0.0011930800746768277, 0.0010608230091035868, 0.0009786996343503485, 0.0011534504169269674, 0.0014689023089074142, 0.0016703854475650266, 0.0016266059694439664, 0.0014314803827807795, 0.0015112451588085024, 0.002078930727349822, 0.0026779349217974932, 0.0029175316230999693, 0.0033235562382948544, 0.006060246682131165, 0.01213564244681699, 0.021528793368822498, 0.034198342832535235, 0.0498509640245669, 0.06777879739323489, 0.08679944339774721, 0.10528219105362474, 0.12126852377137902, 0.1326907328095182, 0.1376810240870295, 0.134958888052281, 0.12431027404132337, 0.10728905812979021, 0.08858742399414933, 0.07815206613344618, 0.08683127365648273, 0.11196336592269487, 0.14327708374162051, 0.17412188732153233, 0.2017067339448505, 0.22574112406214367, 0.2474653597529165, 0.2686093092423925, 0.2903396615900395, 0.3127018045424414, 0.3349372936380533, 0.3565125598804615, 0.3783066367926545, 0.40335892851128957, 0.43667215030086437, 0.483929548957103, 0.5498957343007457, 0.6379296058098175, 0.7511772753359104, 0.8944781464311266, 1.0755419159956716, 1.3045067410940236, 1.5918775562226988, 1.9456678446940345, 2.3688301405631744, 2.857630217749345, 3.4010465241388688, 3.9810277941361747, 4.573462419413825, 5.149771501789923, 5.679030229519041, 6.130462746916682, 6.476086399376074, 6.693233908548924, 6.766672966226011, 6.690076253253228, 6.46666649742484, 6.108959918050503, 5.637642151462975, 5.079716828274695, 4.466152765769695, 3.829308931803635, 3.2004290925717926, 2.6074661984346266, 2.073415157259933, 1.615186605212847, 1.2428149815962508, 0.9585041957707097, 0.7551687763400409, 0.6158471556518877, 0.5173838584728742, 0.4382671702440461, 0.36447672573382195, 0.28989699206373326, 0.21384498776558897, 0.1393329042849531, 0.07870631817659768, 0.08277825381346376, 0.15143123091173274, 0.24078158736248503, 0.34372313632006957, 0.46007260556182084, 0.588984194206225, 0.7275419847731618, 0.870695028957948, 1.0117410440774552, 1.1430757180673803, 1.2570511659931438, 1.346832762341371, 1.4071611528520076, 1.434939087355527, 1.429582545327844, 1.3931062044625082, 1.3299533449035938, 1.2466251999209346, 1.151206480024549, 1.0529068137071977, 0.9617086406356818, 0.8880692703405576, 0.8423369986337814, 0.8333781041605177, 0.8666007756538969, 0.942889207813733, 1.0595493907521107, 1.2121054318929048, 1.3952683154660521, 1.6028409016672367, 1.827214631904138, 2.0590249586673197, 2.2872084231887673, 2.499485517549812, 2.6831822513230446, 2.826250770104734, 2.918330194962036, 2.9516950401135196, 2.9219674505103828, 2.8285183975544346, 2.6745485057201925, 2.4669198730154234, 2.215913057773696, 1.9352373379124859, 1.6428982865974144, 1.3639821670591392, 1.1362514105258539, 1.0127627559570196, 1.036497534091264, 1.1942645746276404, 1.4331065057737395, 1.7072490111044105, 1.988330808603757, 2.258707674286661, 2.505909244379188, 2.7199535755923594, 2.8923450912478703, 3.015934711987843, 3.0851959785499816, 3.096652337031568, 3.0492806751595225, 2.944777132551437, 2.7876272703235907, 2.5849796589951324, 2.3463786229170096, 2.083466033751474, 1.8098118346701337, 1.5410639609903045, 1.295520072890587, 1.094577003547077, 0.9603517317215379, 0.905802049596924, 0.9220515389431886, 0.9816013394002527, 1.0553214038444383, 1.122503758843109, 1.1712260549215174, 1.1958956027208747, 1.1951551363440032, 1.1705013321288313, 1.1253549232801237, 1.0643597055218337, 0.9927910358706103, 0.9160232376541907, 0.8390442881069706, 0.76603012493156, 0.700009739944519, 0.6426670692844786, 0.5943260671192238, 0.5541359268694867, 0.5204167532396715, 0.4910766967643044, 0.46400676401471336, 0.4373946034524892, 0.4099359084201753, 0.38094006785143314, 0.35033851064989413, 0.3186341957132681, 0.28688658296632635, 0.2568807157152263, 0.23158149838991887, 0.2155595344062266, 0.21408112677253166, 0.22973135917132612, 0.2595663201743057, 0.29651428426237564, 0.3323508985988983, 0.3592584785280439, 0.3703823598446884, 0.36022383763028704, 0.32511190218408265, 0.2638351215008961, 0.17920716983789267, 0.09051562530682421, 0.1175036903025314, 0.24835714922391616, 0.39680947956622453, 0.5453911353240208, 0.6846484498492645, 0.8072382983016353, 0.9076087021573619, 0.9821792716861956, 1.0294628025722894, 1.050027509235396, 1.046297789706801, 1.0222307802492925, 0.9829265093378702, 0.9342353299396725, 0.8824119220590887, 0.8338201908136319, 0.7946144295835798, 0.7702482853068296, 0.7647145926640915, 0.7797010448317567, 0.8141482498935529, 0.8645462649144328, 0.9257584320218003, 0.9918867389238565, 1.0569023478448956, 1.115051478911144, 1.1611469623723654, 1.1908245378417512, 1.2007887150552163, 1.1890390588096256, 1.1550559046317896, 1.0999270268944246, 1.0264071446722662, 0.9389146746895276, 0.8434748842580283, 0.747589879793225, 0.6598944967160832, 0.5891806948959527, 0.5421815064653349, 0.5204477579927874, 0.5189240135054122, 0.5283906140396364, 0.5395790781837798, 0.5456843733943463, 0.5428644867743774, 0.5297974893687313, 0.5070326313026419, 0.47636844229022207, 0.44027657075148774, 0.4013505584112461, 0.3617923528539139, 0.32301890708643205, 0.2855428491991873, 0.24929933262599793, 0.2145401519904319, 0.18337942760027146, 0.16182605952420537, 0.15971415205215037, 0.18290860240894036, 0.2265449434629472, 0.2815650501722045, 0.34069528657885395, 0.39886360707757507, 0.45249395568024164, 0.49910037243710986, 0.5371114270499066, 0.5657701653752741, 0.5850383385950925, 0.5954866389733379, 0.5981760983695046, 0.5945406127529805, 0.5862732033961092, 0.5752052688609067, 0.5631575389457308, 0.5517450971214575, 0.5421453243407793, 0.5348811655702548, 0.5297039418623574, 0.5256439853987688, 0.5212306899960217, 0.5148115742169191, 0.5048726534033748, 0.4902856136199302, 0.47045048400259243, 0.44533765831032746, 0.41545212829199757, 0.3817500497358734, 0.3455391985158853, 0.30839336002425044, 0.2721050020200016, 0.23868269450163523, 0.21034842249710095, 0.1893791990680724, 0.17754459085856109, 0.17518395156211494, 0.18072207939821788, 0.1913018020722188, 0.20389125177386663, 0.21595269292166236, 0.22560024707877724, 0.23155666472419054, 0.23310641318837727, 0.23010462835070547, 0.22304780970366572, 0.21319593500493672, 0.20270266647385826, 0.19460178143620163, 0.19231243729964898, 0.19840552854833535, 0.21325695283446064, 0.2348857830534279, 0.260030510129927, 0.28528295088316036, 0.3076723418106493, 0.32485754682265044 ] }, { "line": { "color": "#d62728", "width": 3 }, "name": "fit", "type": "scatter", "uid": "7f8ab199-5c0c-4612-a35d-2459ffeba8a3", "x": [ 0, 0.030679615757712823, 0.06135923151542565, 0.09203884727313846, 0.1227184630308513, 0.15339807878856412, 0.18407769454627693, 0.21475731030398976, 0.2454369260617026, 0.2761165418194154, 0.30679615757712825, 0.3374757733348411, 0.36815538909255385, 0.3988350048502667, 0.4295146206079795, 0.46019423636569234, 0.4908738521234052, 0.521553467881118, 0.5522330836388308, 0.5829126993965437, 0.6135923151542565, 0.6442719309119693, 0.6749515466696822, 0.705631162427395, 0.7363107781851077, 0.7669903939428205, 0.7976700097005334, 0.8283496254582462, 0.859029241215959, 0.8897088569736719, 0.9203884727313847, 0.9510680884890975, 0.9817477042468103, 1.012427320004523, 1.043106935762236, 1.0737865515199487, 1.1044661672776617, 1.1351457830353744, 1.1658253987930873, 1.1965050145508, 1.227184630308513, 1.2578642460662257, 1.2885438618239387, 1.3192234775816514, 1.3499030933393643, 1.380582709097077, 1.41126232485479, 1.4419419406125027, 1.4726215563702154, 1.5033011721279284, 1.533980787885641, 1.564660403643354, 1.5953400194010667, 1.6260196351587797, 1.6566992509164924, 1.6873788666742053, 1.718058482431918, 1.748738098189631, 1.7794177139473437, 1.8100973297050567, 1.8407769454627694, 1.8714565612204823, 1.902136176978195, 1.932815792735908, 1.9634954084936207, 1.9941750242513334, 2.024854640009046, 2.055534255766759, 2.086213871524472, 2.116893487282185, 2.1475731030398975, 2.1782527187976104, 2.2089323345553233, 2.2396119503130363, 2.270291566070749, 2.3009711818284617, 2.3316507975861747, 2.3623304133438876, 2.3930100291016, 2.423689644859313, 2.454369260617026, 2.4850488763747385, 2.5157284921324514, 2.5464081078901644, 2.5770877236478773, 2.60776733940559, 2.6384469551633027, 2.6691265709210157, 2.6998061866787286, 2.730485802436441, 2.761165418194154, 2.791845033951867, 2.82252464970958, 2.8532042654672924, 2.8838838812250054, 2.9145634969827183, 2.945243112740431, 2.9759227284981438, 3.0066023442558567, 3.0372819600135696, 3.067961575771282, 3.098641191528995, 3.129320807286708, 3.160000423044421, 3.1906800388021335, 3.2213596545598464, 3.2520392703175593, 3.2827188860752723, 3.313398501832985, 3.3440781175906977, 3.3747577333484107, 3.4054373491061236, 3.436116964863836, 3.466796580621549, 3.497476196379262, 3.5281558121369745, 3.5588354278946874, 3.5895150436524004, 3.6201946594101133, 3.650874275167826, 3.6815538909255388, 3.7122335066832517, 3.7429131224409646, 3.773592738198677, 3.80427235395639, 3.834951969714103, 3.865631585471816, 3.8963112012295285, 3.9269908169872414, 3.9576704327449543, 3.988350048502667, 4.01902966426038, 4.049709280018092, 4.080388895775806, 4.111068511533518, 4.1417481272912315, 4.172427743048944, 4.2031073588066565, 4.23378697456437, 4.264466590322082, 4.295146206079795, 4.325825821837508, 4.356505437595221, 4.387185053352934, 4.417864669110647, 4.448544284868359, 4.479223900626073, 4.509903516383785, 4.540583132141498, 4.571262747899211, 4.601942363656923, 4.632621979414636, 4.663301595172349, 4.693981210930062, 4.724660826687775, 4.755340442445488, 4.7860200582032, 4.816699673960914, 4.847379289718626, 4.878058905476339, 4.908738521234052, 4.9394181369917645, 4.970097752749477, 5.00077736850719, 5.031456984264903, 5.062136600022616, 5.092816215780329, 5.123495831538041, 5.154175447295755, 5.184855063053467, 5.21553467881118, 5.246214294568893, 5.2768939103266055, 5.307573526084319, 5.338253141842031, 5.368932757599744, 5.399612373357457, 5.43029198911517, 5.460971604872882, 5.491651220630596, 5.522330836388308, 5.553010452146021, 5.583690067903734, 5.6143696836614465, 5.64504929941916, 5.675728915176872, 5.706408530934585, 5.737088146692298, 5.767767762450011, 5.798447378207723, 5.829126993965437, 5.859806609723149, 5.890486225480862, 5.921165841238575, 5.9518454569962875, 5.982525072754001, 6.013204688511713, 6.043884304269426, 6.074563920027139, 6.105243535784852, 6.135923151542564, 6.166602767300278, 6.19728238305799, 6.227961998815704, 6.258641614573416, 6.2893212303311286, 6.320000846088842, 6.350680461846554, 6.381360077604267, 6.41203969336198, 6.442719309119693, 6.473398924877405, 6.504078540635119, 6.534758156392831, 6.565437772150545, 6.596117387908257, 6.62679700366597, 6.657476619423683, 6.6881562351813955, 6.718835850939108, 6.749515466696821, 6.780195082454534, 6.810874698212247, 6.84155431396996, 6.872233929727672, 6.902913545485386, 6.933593161243098, 6.964272777000811, 6.994952392758524, 7.0256320085162365, 7.056311624273949, 7.086991240031662, 7.117670855789375, 7.148350471547088, 7.179030087304801, 7.209709703062513, 7.240389318820227, 7.271068934577939, 7.301748550335652, 7.332428166093365, 7.3631077818510775, 7.39378739760879, 7.424467013366503, 7.455146629124216, 7.485826244881929, 7.516505860639642, 7.547185476397354, 7.577865092155068, 7.60854470791278, 7.639224323670493, 7.669903939428206, 7.7005835551859185, 7.731263170943632, 7.761942786701344, 7.792622402459057, 7.82330201821677, 7.853981633974483, 7.884661249732195, 7.915340865489909, 7.946020481247621, 7.976700097005334, 8.007379712763047, 8.03805932852076, 8.068738944278472, 8.099418560036185, 8.130098175793899, 8.160777791551611, 8.191457407309324, 8.222137023067036, 8.252816638824749, 8.283496254582463, 8.314175870340176, 8.344855486097888, 8.3755351018556, 8.406214717613313, 8.436894333371026, 8.46757394912874, 8.498253564886452, 8.528933180644165, 8.559612796401877, 8.59029241215959, 8.620972027917304, 8.651651643675017, 8.682331259432729, 8.713010875190442, 8.743690490948154, 8.774370106705868, 8.80504972246358, 8.835729338221293, 8.866408953979006, 8.897088569736718, 8.92776818549443, 8.958447801252145, 8.989127417009858, 9.01980703276757, 9.050486648525283, 9.081166264282995, 9.11184588004071, 9.142525495798422, 9.173205111556134, 9.203884727313847, 9.23456434307156, 9.265243958829272, 9.295923574586986, 9.326603190344699, 9.357282806102411, 9.387962421860124, 9.418642037617836, 9.44932165337555, 9.480001269133263, 9.510680884890975, 9.541360500648688, 9.5720401164064, 9.602719732164113, 9.633399347921827, 9.66407896367954, 9.694758579437252, 9.725438195194965, 9.756117810952677, 9.786797426710391, 9.817477042468104, 9.848156658225816, 9.878836273983529, 9.909515889741241, 9.940195505498954, 9.970875121256668 ], "y": [ 0.0008750484537326243, 0.0011360340002200368, 0.0016595124593535358, 0.0021735364097790603, 0.002533647399629476, 0.0026669662553717786, 0.0025982553180954335, 0.002455675614090537, 0.00236482019909277, 0.00226628591123589, 0.0019400231919883188, 0.0013325833386259337, 0.0014831888458286945, 0.0031236184503889404, 0.005202067860233055, 0.0071522439789152, 0.00863267566741532, 0.009562161177007323, 0.010239977790888538, 0.011260899363237648, 0.013020321969257186, 0.015294339365073085, 0.017526113642065848, 0.019404591406862928, 0.02126620964249914, 0.024095386718982863, 0.028826832002760178, 0.035489114212211154, 0.04326938463326748, 0.0511659751407212, 0.058460675021528664, 0.06487664455028769, 0.07046897654215473, 0.07531776755594265, 0.0792214815206071, 0.0816596762928565, 0.08214484242310299, 0.08084139262995986, 0.07916262774677058, 0.07979576496914777, 0.08536480113729247, 0.09635913363416126, 0.11100576008194704, 0.12718679436232047, 0.14401560922327666, 0.16220294716814132, 0.18338133654862712, 0.20882604592235035, 0.23849397482804968, 0.27116205675418725, 0.30550333685084013, 0.34135030018919205, 0.3804748566425326, 0.4264471915440487, 0.4835087429832948, 0.5551003570945215, 0.6431598318456433, 0.7487689475965826, 0.8737133350038732, 1.0220670187399683, 1.201031790777001, 1.4205262702694195, 1.6914274113358976, 2.022949130328954, 2.420024711536505, 2.8813968924643714, 3.3986501370401565, 3.9561125220657316, 4.53151886870787, 5.097378872504928, 5.62299021558978, 6.076963888066189, 6.430032730733529, 6.657838585145212, 6.743364052140001, 6.67870044849556, 6.465920152999401, 6.1169360425484935, 5.652364134313627, 5.099536394036214, 4.489918499931071, 3.856255923880868, 3.2297912511063824, 2.637863873006916, 2.1021244489361353, 1.6374798933354684, 1.251742759510891, 0.9458177639784672, 0.7141998564020298, 0.5457739179165003, 0.42546144612465053, 0.33732001432243686, 0.26819900169742994, 0.2097951466959852, 0.15834029409266065, 0.11306034427027123, 0.07465218981698808, 0.04429028953297575, 0.023338751979668546, 0.013681274260388246, 0.012961125733453998, 0.013258873908515873, 0.013199849137377033, 0.01494118749786803, 0.018650395927394838, 0.02235852304439501, 0.024437678197686136, 0.024202283092862292, 0.02176850112038076, 0.01786201224617372, 0.013618299930848584, 0.01029309379796276, 0.008532020525545747, 0.007556056188546769, 0.006126945404018424, 0.004077816060825335, 0.003870770495925518, 0.0071608824389305685, 0.011144972784431894, 0.014470302407469557, 0.01650334416950208, 0.016971226131610347, 0.01592471732368329, 0.013704681950905871, 0.010884745001501564, 0.008197933533614024, 0.006379013760311803, 0.005665245818274821, 0.005441023553175286, 0.005033404860800883, 0.004353781739149761, 0.00387285936062667, 0.0041477870960556565, 0.0049614803966172485, 0.0056652644690757354, 0.005852679141429397, 0.0054268090972152145, 0.0045598260125352215, 0.003691652530822212, 0.0034128789928146402, 0.003813313195985804, 0.004334392680146898, 0.004548388346102445, 0.004383317998385634, 0.004066632636210979, 0.004005043775000525, 0.00441018931256939, 0.00503253780788854, 0.005474928270250118, 0.005478384745702302, 0.004967167846476999, 0.004028357319268386, 0.0028948353425790295, 0.0019757184747235617, 0.001793261194919291, 0.002174569545743204, 0.0025425495364807525, 0.0027558793257100368, 0.0029641914723262694, 0.003354758809651114, 0.003924783940472646, 0.004490872895289644, 0.0048511444709134635, 0.004881024338575018, 0.004558180330676638, 0.003958140595000579, 0.0032350097670411147, 0.0025843434670336984, 0.0021591264954744283, 0.0019410431257140053, 0.0017566078142462155, 0.001468610187587478, 0.0010873274209250749, 0.0008252711674867783, 0.0009758775611061196, 0.00133432252073266, 0.0016045856664572416, 0.0016699872113989163, 0.0015099566086256427, 0.0011641049921927967, 0.0007144834067291699, 0.0002927553327710803, 0.00029576700304621617, 0.0005299591809355848, 0.0006658752962033732, 0.00075634096514998, 0.0009428614748496741, 0.0012850790442649276, 0.0017069188426993317, 0.0021020237711390525, 0.0023884479074691603, 0.00252268655660062, 0.0025070065684429277, 0.002392371264164414, 0.002268612759445359, 0.0022245980794399337, 0.002283228741552589, 0.0023818669392967964, 0.002429630790749837, 0.002369117362766575, 0.0022041527254084503, 0.0020057593621210707, 0.0018911679835943412, 0.0019410852012167836, 0.0021115342787386467, 0.0022824072897297256, 0.002358722247741991, 0.0023101322630277796, 0.002171728988370778, 0.0020272053145083376, 0.001962179635867542, 0.0019937262675245146, 0.0020544399531374088, 0.0020578062639870676, 0.0019589645641279177, 0.001774913421983477, 0.0015801257265844828, 0.0014716972237871606, 0.0014866813594147474, 0.0015570560944652713, 0.0015835284870767922, 0.0015101114537947454, 0.0013407882048326548, 0.0011328173848291117, 0.0009756826280890534, 0.0009266405081373758, 0.0009395973309850829, 0.0009194503003853684, 0.0008098190956207661, 0.0006117691097840789, 0.0003862138518548717, 0.0003018489118217902, 0.0004347167828217587, 0.0005765562469373891, 0.000644760537081493, 0.0006582898227283857, 0.0006985173720217634, 0.0008394331441237582, 0.0010592896709659175, 0.0012838066560259503, 0.001455411741817875, 0.001553338021337132, 0.0015956279178440328, 0.001629851729975065, 0.0017047894803258394, 0.0018330117684510878, 0.001981819587496369, 0.002101154353036572, 0.0021546959645031105, 0.002135086626168645, 0.0020653933600250085, 0.0019890993933427397, 0.001947753308724039, 0.0019542199602493836, 0.001984505580913656, 0.0019974390461642206, 0.0019610676690979983, 0.0018671986467684253, 0.0017333120906291835, 0.0015950771546645772, 0.0014892583124445303, 0.001430446080331936, 0.0013997336453496882, 0.0013596657663193614, 0.0012791654425923235, 0.0011475870749208291, 0.0009763473136831848, 0.0007935863146575018, 0.000633376896857019, 0.0005169016132683976, 0.00043331785064168536, 0.0003495044575380862, 0.00024720491518262673, 0.00016964399845562547, 0.00023625770252822995, 0.00038460191389528566, 0.0005262984226238868, 0.0006350995152682573, 0.000714917219914814, 0.0007914182266285711, 0.0008943143748307354, 0.0010314038057017056, 0.0011814103467560704, 0.0013110178732011867, 0.0013942875308516226, 0.0014236975051172146, 0.0014131261701766443, 0.0013920625366338317, 0.0013891967011836642, 0.0014120420802047195, 0.0014411003332496653, 0.0014451891416102692, 0.0014015794378921531, 0.0013070898084842919, 0.0011790146542751362, 0.0010475982370960434, 0.0009400321396523084, 0.0008613342131073972, 0.0007901855426176807, 0.0006973683624647614, 0.0005662452613796144, 0.000399644357556042, 0.00021650533059884642, 5.378466510672121e-05, 0.0001187882242740693, 0.00022190360074468294, 0.00029790311646722944, 0.00037615654887456336, 0.0004846852505569301, 0.0006221322173470865, 0.0007639818682103863, 0.0008832664886507676, 0.0009634474155378504, 0.0010041144977438773, 0.0010212285389833184, 0.0010398888577031016, 0.001079020236987495, 0.0011375874212394461, 0.0011963504280390073, 0.0012324019776250043, 0.0012324520079079588, 0.001198664559882683, 0.0011475567801600903, 0.0011022819982601498, 0.0010792322555494554, 0.001076670574648924 ] } ], "layout": { "height": 500, "hovermode": "closest", "legend": { "bgcolor": "#F2F2F2", "borderwidth": 0.5 }, "plot_bgcolor": "#FDFDFF", "showlegend": true, "template": { "data": { "bar": [ { "error_x": { "color": "#2a3f5f" }, "error_y": { "color": "#2a3f5f" }, "marker": { "line": { "color": "#E5ECF6", "width": 0.5 }, "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "bar" } ], "barpolar": [ { "marker": { "line": { "color": "#E5ECF6", "width": 0.5 }, "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "barpolar" } ], "carpet": [ { "aaxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "baxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "type": "carpet" } ], "choropleth": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "choropleth" } ], "contour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "contour" } ], "contourcarpet": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "contourcarpet" } ], "heatmap": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmap" } ], "heatmapgl": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmapgl" } ], "histogram": [ { "marker": { "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "histogram" } ], "histogram2d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2d" } ], "histogram2dcontour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2dcontour" } ], "mesh3d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "mesh3d" } ], "parcoords": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "parcoords" } ], "pie": [ { "automargin": true, "type": "pie" } ], "scatter": [ { "fillpattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 }, "type": "scatter" } ], "scatter3d": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatter3d" } ], "scattercarpet": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattercarpet" } ], "scattergeo": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergeo" } ], "scattergl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergl" } ], "scattermapbox": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattermapbox" } ], "scatterpolar": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolar" } ], "scatterpolargl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolargl" } ], "scatterternary": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterternary" } ], "surface": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "surface" } ], "table": [ { "cells": { "fill": { "color": "#EBF0F8" }, "line": { "color": "white" } }, "header": { "fill": { "color": "#C8D4E3" }, "line": { "color": "white" } }, "type": "table" } ] }, "layout": { "annotationdefaults": { "arrowcolor": "#2a3f5f", "arrowhead": 0, "arrowwidth": 1 }, "autotypenumbers": "strict", "coloraxis": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "colorscale": { "diverging": [ [ 0, "#8e0152" ], [ 0.1, "#c51b7d" ], [ 0.2, "#de77ae" ], [ 0.3, "#f1b6da" ], [ 0.4, "#fde0ef" ], [ 0.5, "#f7f7f7" ], [ 0.6, "#e6f5d0" ], [ 0.7, "#b8e186" ], [ 0.8, "#7fbc41" ], [ 0.9, "#4d9221" ], [ 1, "#276419" ] ], "sequential": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "sequentialminus": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ] }, "colorway": [ "#636efa", "#EF553B", "#00cc96", "#ab63fa", "#FFA15A", "#19d3f3", "#FF6692", "#B6E880", "#FF97FF", "#FECB52" ], "font": { "color": "#2a3f5f" }, "geo": { "bgcolor": "white", "lakecolor": "white", "landcolor": "#E5ECF6", "showlakes": true, "showland": true, "subunitcolor": "white" }, "hoverlabel": { "align": "left" }, "hovermode": "closest", "mapbox": { "style": "light" }, "paper_bgcolor": "white", "plot_bgcolor": "#E5ECF6", "polar": { "angularaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "radialaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "scene": { "xaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "yaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "zaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" } }, "shapedefaults": { "line": { "color": "#2a3f5f" } }, "ternary": { "aaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "baxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "caxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "title": { "x": 0.05 }, "xaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 }, "yaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 } } }, "title": {}, "width": 650, "xaxis": { "color": "#004", "gridcolor": "#D8D8D8", "showgrid": true, "title": {}, "zerolinecolor": "#DDD" }, "yaxis": { "color": "#004", "gridcolor": "#D8D8D8", "showgrid": true, "title": {}, "zerolinecolor": "#DDD" } } }, "text/html": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# let's now write a function to do a first shell fit to the \n", "# XAFS data from a dataset (Group)\n", "def fit_cu_1stshell(dat, full_report=False, with_plot=False):\n", " \n", " # extract the EXAFS chi(k)\n", " autobk(dat.energy, dat.mu, group=dat, rbkg=1.0, kw=2)\n", "\n", " # define fitting parameter group\n", " pars = param_group(amp = param(1.0, vary=True),\n", " del_e0 = param(0.0, vary=True),\n", " sig2 = param(0.0, vary=True),\n", " del_r = guess(0.0, vary=True) )\n", "\n", " # define a Feff Path, give expressions for Path Parameters\n", " path1 = feffpath(feffdat_file,\n", " s02 = 'amp',\n", " e0 = 'del_e0',\n", " sigma2 = 'sig2',\n", " deltar = 'del_r')\n", "\n", " # do the fit\n", " trans = feffit_transform(kmin=2.5, kmax=15, kw=2, dk=4, window='kaiser',\n", " rmin=1.4, rmax=3.0)\n", " dset = feffit_dataset(data=dat, pathlist=[path1], transform=trans)\n", " out = feffit(pars, dset)\n", " \n", " r1_fit = out.params['del_r'].value + path1.reff\n", " r1_err = out.params['del_r'].stderr\n", " print(f'Result: E0 = {dat.e0:.3f} eV, R= {r1_fit:.4f} +/- {r1_err:.4f} Ang')\n", " if full_report:\n", " print(feffit_report(out))\n", " if with_plot:\n", " multi_plot([dict(xdata=dset.data.r, ydata=dset.data.chir_mag, label='data',\n", " xlabel='R (Ang)', ylabel='$|\\chi(R)|$'),\n", " dict(xdata=dset.model.r, ydata=dset.model.chir_mag, label='fit')])\n", " return out\n", " \n", "out = fit_cu_1stshell(cudat, full_report=True, with_plot=True)\n", " " ] }, { "cell_type": "code", "execution_count": 6, "id": "fa9fd5d2", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "theta_off=0.0200, dspace_off=0.0000: E0= 8977.580 to 8991.481\n" ] } ], "source": [ "# OK, so we found E0=8977 eV, and can do a fit with this dataset and get \n", "# R = 2.5494(0017) Ang.\n", "\n", "# If we take the angles for this energy array and change the offset, \n", "# what do we get? \n", "\n", "# we define two functions to shift X-ray energy either by changing the \n", "# angular offset or the d-spacing of the monochromator\n", "\n", "def energy_shift(energy, dspace, theta_off=0, dspace_off=0):\n", " theta = RAD2DEG * np.arcsin(PLANCK_HC / (2.0*dspace*energy))\n", " return PLANCK_HC/(2*(dspace-dspace_off)*np.sin(DEG2RAD*(theta-theta_off)))\n", "\n", "def change_calib(dgroup, dspace, theta_off=0, dspace_off=0):\n", " \"make a new data group with an offset in theta and/or dspacing\"\n", " enew = energy_shift(dgroup.energy, dspace, theta_off=theta_off, \n", " dspace_off=dspace_off)\n", "\n", " dat = Group(energy=enew, mu=dgroup.mu*1.0)\n", " autobk(dat.energy, dat.mu, group=dat, rbkg=1.0, kw=2)\n", " print(f'theta_off={theta_off:.4f}, dspace_off={dspace_off:.4f}: E0= {dgroup.e0:.3f} to {dat.e0:.3f}')\n", " return dat\n", "\n", "# and let's see how big a shift of 0.02 degrees is\n", "# For reference, 0.02 degrees ~= 350 microradians, \n", "# which is about 10x bigger than the Darwin width for Si(111).\n", "new_dat = change_calib(cudat, dspace, theta_off=0.02)" ] }, { "cell_type": "code", "execution_count": 7, "id": "9c7a99c3", "metadata": {}, "outputs": [ { "data": { "application/vnd.plotly.v1+json": { "config": { "plotlyServerURL": "https://plot.ly" }, "data": [ { "line": { "color": "#1f77b4", "width": 3 }, "name": "E0=8977.6", "type": "scatter", "uid": "85995838-d99e-4ac0-b941-017ce447f3f8", "x": [ 8786.204, 8796.258, 8806.253, 8816.27, 8826.27, 8836.292, 8846.337, 8856.322, 8866.372, 8876.361, 8886.416, 8896.409, 8906.425, 8916.464, 8926.483, 8936.483, 8936.998, 8937.471, 8937.944, 8938.417, 8938.976, 8939.406, 8939.878, 8940.352, 8940.824, 8941.384, 8941.814, 8942.287, 8942.761, 8943.319, 8943.793, 8944.224, 8944.697, 8945.257, 8945.73, 8946.204, 8946.635, 8947.108, 8947.668, 8948.143, 8948.616, 8949.048, 8949.564, 8950.038, 8950.513, 8950.987, 8951.461, 8951.979, 8952.453, 8952.927, 8953.401, 8953.876, 8954.394, 8954.868, 8955.343, 8955.817, 8956.336, 8956.811, 8957.285, 8957.761, 8958.235, 8958.71, 8959.229, 8959.703, 8960.179, 8960.654, 8961.13, 8961.647, 8962.123, 8962.599, 8963.074, 8963.55, 8964.068, 8964.544, 8965.02, 8965.496, 8965.972, 8966.447, 8966.967, 8967.442, 8967.919, 8968.395, 8968.871, 8969.347, 8969.909, 8970.386, 8970.862, 8971.339, 8971.772, 8972.248, 8972.812, 8973.288, 8973.765, 8974.198, 8974.675, 8975.151, 8975.628, 8976.192, 8976.626, 8977.103, 8977.58, 8978.057, 8978.534, 8979.055, 8979.531, 8980.009, 8980.486, 8980.964, 8981.441, 8981.919, 8982.396, 8982.918, 8983.396, 8983.873, 8984.352, 8984.829, 8985.307, 8985.785, 8986.307, 8986.784, 8987.263, 8987.741, 8988.219, 8988.697, 8989.176, 8989.653, 8990.132, 8990.698, 8991.133, 8991.611, 8992.09, 8992.568, 8993.048, 8993.526, 8994.006, 8994.484, 8995.007, 8995.485, 8995.965, 8996.444, 8996.924, 8997.402, 8997.882, 8998.361, 8998.84, 8999.319, 8999.799, 9000.366, 9000.803, 9001.281, 9001.762, 9002.241, 9002.721, 9003.201, 9003.681, 9004.161, 9004.642, 9005.121, 9005.602, 9006.125, 9006.605, 9007.086, 9007.566, 9008.047, 9008.526, 9009.007, 9009.487, 9009.969, 9011.062, 9012.153, 9013.246, 9014.384, 9015.565, 9016.746, 9017.928, 9019.154, 9020.38, 9021.606, 9022.877, 9024.191, 9025.464, 9026.822, 9028.14, 9029.499, 9030.903, 9032.309, 9033.759, 9035.164, 9036.615, 9038.109, 9039.604, 9041.145, 9042.642, 9044.271, 9045.768, 9047.397, 9049.029, 9050.616, 9052.336, 9053.969, 9055.646, 9057.324, 9059.047, 9060.814, 9062.627, 9064.353, 9066.166, 9067.981, 9069.842, 9071.702, 9073.563, 9075.47, 9077.377, 9079.329, 9081.283, 9083.237, 9085.281, 9087.237, 9089.283, 9091.33, 9093.378, 9095.472, 9097.61, 9099.705, 9101.847, 9104.033, 9106.221, 9108.409, 9110.644, 9112.879, 9115.16, 9117.398, 9119.728, 9122.057, 9124.389, 9126.721, 9129.099, 9131.478, 9133.904, 9136.331, 9138.805, 9141.279, 9143.756, 9146.233, 9148.803, 9151.328, 9153.899, 9156.474, 9159.094, 9161.716, 9164.339, 9167.009, 9169.681, 9172.399, 9175.074, 9177.842, 9180.564, 9183.381, 9186.153, 9188.973, 9191.794, 9194.663, 9197.533, 9200.406, 9203.326, 9206.247, 9209.217, 9212.143, 9215.115, 9218.183, 9221.16, 9224.185, 9227.258, 9230.333, 9233.502, 9236.581, 9239.755, 9242.838, 9246.016, 9249.196, 9252.425, 9255.655, 9258.889, 9262.171, 9265.454, 9268.787, 9272.122, 9275.459, 9278.8, 9282.188, 9285.58, 9289.021, 9292.465, 9295.91, 9299.404, 9302.949, 9306.403, 9309.953, 9313.506, 9317.061, 9320.666, 9324.32, 9327.932, 9331.592, 9335.256, 9338.922, 9342.639, 9346.358, 9350.128, 9353.9, 9357.677, 9361.455, 9365.285, 9369.165, 9373.048, 9376.84, 9380.777, 9384.671, 9388.614, 9392.563, 9396.561, 9400.515, 9404.567, 9408.576, 9412.637, 9416.7, 9420.768, 9424.935, 9429.01, 9433.184, 9437.313, 9441.495, 9445.729, 9449.966, 9454.208, 9458.453, 9462.702, 9467.004, 9471.358, 9475.668, 9480.029, 9484.396, 9488.814, 9493.189, 9497.665, 9502.097, 9506.581, 9511.118, 9515.563, 9520.109, 9524.611, 9529.265, 9533.775, 9538.438, 9543.007, 9547.679, 9552.306, 9556.986, 9561.673, 9566.462, 9571.157, 9575.956, 9580.661, 9585.47, 9590.283, 9595.102, 9599.976, 9604.854, 9609.737, 9614.625, 9619.568, 9624.518, 9629.471, 9634.431, 9639.445, 9644.466, 9649.541, 9654.571, 9659.708, 9664.75, 9669.847, 9675.001, 9680.159, 9685.324, 9690.545, 9695.722, 9700.953, 9706.243, 9711.486, 9716.788, 9722.094, 9727.407, 9732.776, 9738.204, 9743.535, 9748.923, 9754.317, 9759.769, 9765.227, 9770.742, 9776.212, 9781.793, 9787.275, 9792.816, 9798.416, 9804.022, 9809.584, 9815.254, 9820.88, 9826.564, 9832.308, 9837.953, 9843.659, 9849.422, 9855.192, 9860.971, 9866.755, 9872.6, 9878.45, 9884.31, 9890.175, 9896.1, 9901.98, 9907.92, 9913.867, 9919.876, 9925.838, 9931.859, 9937.89, 9943.98, 9950.025, 9956.132, 9962.3, 9968.367, 9974.497, 9980.687, 9986.832, 9993.038, 9999.253, 10005.53, 10011.76, 10018.05, 10024.35, 10030.66, 10037.03, 10043.36, 10049.74, 10056.14, 10062.54, 10068.96, 10075.43, 10081.92, 10088.46, 10094.91, 10101.42, 10107.94, 10114.52, 10121.11, 10127.71, 10134.32, 10141.05, 10147.68, 10154.31, 10161.01, 10167.78, 10174.49, 10181.28, 10188.07, 10194.82, 10201.63, 10208.45, 10215.28, 10222.17, 10229.13, 10235.99, 10242.86, 10249.85, 10256.79, 10263.8, 10270.82, 10277.85, 10284.89, 10291.88, 10299.05, 10306.12, 10313.2, 10320.41, 10327.5, 10334.73, 10341.85, 10349.09, 10356.35, 10363.61, 10370.89, 10378.17, 10385.47, 10392.78, 10400.1, 10407.48, 10414.82, 10422.23, 10429.65, 10437.14, 10444.58, 10452.03, 10459.55, 10467.08, 10474.63, 10482.18, 10489.75, 10497.32, 10504.97, 10512.57, 10520.24, 10527.92, 10535.62, 10543.32, 10551.04, 10558.77, 10566.5, 10574.32, 10582.2, 10589.91, 10597.76, 10605.62, 10613.55, 10621.43, 10629.33, 10637.3, 10645.21, 10653.21, 10661.21, 10669.23, 10677.26, 10685.3, 10693.35, 10701.54, 10709.56, 10717.71, 10725.81, 10733.99, 10742.24, 10750.38, 10758.6, 10766.83, 10775.14, 10783.39, 10791.72, 10800, 10808.36, 10816.73, 10825.05, 10833.45, 10841.86, 10850.41, 10858.78, 10867.23, 10875.76, 10884.23, 10892.85, 10901.29, 10909.94, 10918.47, 10927.08, 10935.7, 10944.4, 10952.99, 10961.65, 10970.4, 10979.09, 10987.86, 10996.65, 11005.32, 11014.14, 11022.97, 11031.81, 11040.67, 11049.54, 11058.43, 11067.33, 11076.25, 11085.18, 11094.12, 11103.15, 11112.13, 11121.25, 11130.19, 11139.34, 11148.38, 11157.43, 11166.56, 11175.78, 11184.87, 11194.05, 11203.31, 11212.45, 11221.68, 11230.98, 11240.17, 11249.51, 11258.8, 11268.17, 11277.55, 11286.82, 11296.24, 11305.6, 11315.05, 11324.52, 11333.93, 11343.43, 11352.94, 11362.47 ], "y": [ 0.006099431461700133, 9.117887986916555e-05, 5.4670997746302065e-05, -2.7191318620089223e-05, -6.938989746250113e-05, -0.00010056449982715447, -7.513226728440632e-05, -1.3385541299481825e-05, 5.272897741641901e-05, 8.708466946053892e-05, 0.00017065148104311614, 0.00041148213917538364, 0.0005611003115116054, 0.0009345899709346704, 0.001294474862171137, 0.0021675119184885456, 0.0022607822078094674, 0.0023875422375935563, 0.0023762273633420768, 0.0023721337978712797, 0.0025079414162668275, 0.002596333204568099, 0.0026458845236113394, 0.002719086672601256, 0.002780035478997492, 0.0028462330302125366, 0.0028928195309341903, 0.002940493179011494, 0.0029652777576828714, 0.0030315329927198486, 0.003082027170724749, 0.003230835951758087, 0.003293082722977108, 0.0033462297161551994, 0.003441320391768271, 0.0035178286821276475, 0.003617741372381377, 0.0036910811179321993, 0.0037889915251784114, 0.003887461163920137, 0.0039618884559740304, 0.004083327476008513, 0.0041815063845434056, 0.004229346949080708, 0.004403031303977713, 0.004464922969334836, 0.004629870541327335, 0.004732298450192962, 0.004796756725297867, 0.004937082244459865, 0.005072579057148074, 0.005231211768441568, 0.005365396035198512, 0.005503415955773532, 0.005692760980315324, 0.00584918218772303, 0.006005023289841032, 0.006139512468765868, 0.006347526882281607, 0.0064916228813130285, 0.006707163116630306, 0.006988601579055301, 0.007225965167046599, 0.007393348843206092, 0.007674388697909812, 0.007874130667021008, 0.008181967667561034, 0.008389155146587373, 0.008714218883736713, 0.009033975393950906, 0.009430780789426454, 0.009781902140790373, 0.01024961516912145, 0.010635494506724855, 0.011077012723427917, 0.011505654389533273, 0.012058319858854423, 0.012576538945937448, 0.013250784624178266, 0.013968024721330552, 0.014837420143237271, 0.015396132371115685, 0.016182663840465787, 0.017379504854510168, 0.018681054047579807, 0.019737420464301995, 0.021237686952714403, 0.02338142586510873, 0.02435640779166551, 0.026455259875391095, 0.029764405731245495, 0.03331113022549566, 0.03785164762330877, 0.041459318622740333, 0.05100611708740155, 0.05791562237299598, 0.06834898874699495, 0.09645415137195536, 0.1310924751406965, 0.19162954393169235, 0.29135526341584717, 0.3803821354209815, 0.44203154677536083, 0.4717576761114283, 0.48211339875608555, 0.47868162767092126, 0.4696792329784953, 0.45905921452635834, 0.45370314574956777, 0.460077563973669, 0.4715597601952697, 0.49142640335243826, 0.5145646983941182, 0.533314315368118, 0.5491166671367219, 0.5660078846462031, 0.5931265998892751, 0.6100014971102291, 0.6415096730553159, 0.6720896208901742, 0.7170621449652635, 0.7511228870470033, 0.7992947625159293, 0.8462098692458588, 0.8807676129885406, 0.9149103631675959, 0.957964020192593, 1.0050294543708573, 1.023806775337923, 1.043785095260159, 1.0472594548880798, 1.0409952776902094, 1.0314432157221682, 1.0148108051822236, 0.9982605989759534, 0.9809573897528895, 0.9620462501947948, 0.9475567412845837, 0.939415879658981, 0.9337626789955784, 0.9314536767382668, 0.9352097180014745, 0.9441046198192283, 0.9585810348912375, 0.9864468874909668, 1.0012952434670397, 1.0253875201192892, 1.0474691764091684, 1.0639987420816401, 1.0716903966303117, 1.0704552306991924, 1.0622623914417215, 1.0439051179658794, 1.033929782898826, 1.0160123132288432, 1.006598152157396, 0.982101671580503, 0.9711447241307211, 0.9575593801102559, 0.9439850487522351, 0.935379587260503, 0.926429366436303, 0.9145634405282396, 0.9063584065679724, 0.8975626315291759, 0.8916597160066978, 0.8881291815170896, 0.8873744865400649, 0.8921382376259682, 0.898582479750442, 0.9095865726712046, 0.9219882787729393, 0.9452540520192998, 0.9673239591308375, 0.9847572133900242, 1.0101331602861956, 1.034947063149558, 1.0543436793716165, 1.069669289938395, 1.07162799654185, 1.0586169162601677, 1.0324870298152853, 1.0070408412321845, 0.9768354759305574, 0.9539400616194232, 0.9388914986555852, 0.9286304384139012, 0.926698780915085, 0.9311711410786876, 0.9416993077873551, 0.9514925034782377, 0.9544538408700594, 0.9505117450632707, 0.9462159131614261, 0.9493802812245763, 0.9614241830040963, 0.9738532724820766, 0.9770239626545879, 0.9749541281839764, 0.976383112710594, 0.9860384220022077, 0.9997852916568509, 1.0123174877919885, 1.0205431227902535, 1.0287801771056422, 1.0385941382379584, 1.0534782767803708, 1.0688987790761573, 1.080862732399942, 1.0816425820859978, 1.0709190102990507, 1.0485318286090124, 1.0183369330638294, 0.9934030508972742, 0.9741765985970925, 0.9573353532875375, 0.9483594761374367, 0.9468426169119785, 0.9376883679534058, 0.9146467955870067, 0.8935177754100005, 0.8783586243313712, 0.8752596169505101, 0.8846152795261619, 0.9008402050672836, 0.9219339818089923, 0.9453459584947743, 0.967569029343354, 0.982210045224044, 0.9856441118994594, 0.9789802006176652, 0.9708634146105803, 0.9686652743404159, 0.98666875247515, 1.02993352200194, 1.084369569208702, 1.1122391392709219, 1.1138429316268879, 1.1015133974367877, 1.0893041387989504, 1.0751454460421443, 1.048832768585153, 1.0134927813514836, 0.9702354555751549, 0.9232465333040473, 0.8770964284809694, 0.8375387043746495, 0.8131230793765442, 0.8105940086594493, 0.8321323278742748, 0.8699657842728565, 0.9144644799129241, 0.9578799650539309, 0.9866828491940066, 1.0023647896861978, 0.999979530698181, 0.9836309358139109, 0.9683620359826701, 0.961117155836223, 0.9584139645383036, 0.9467626510310942, 0.9295096976471825, 0.9262762991439139, 0.9469869557081324, 0.9898368981787231, 1.0325974372357776, 1.04757973356158, 1.049494818127731, 1.0432286166937583, 1.025126618460171, 0.9927860191405699, 0.9524872841813219, 0.9147826597928286, 0.8862377728901889, 0.8691971563824952, 0.8618789737636843, 0.8652035001090724, 0.879343016651351, 0.8898809930786616, 0.8901420758881987, 0.8929287467478129, 0.8970219888852865, 0.8990062324800004, 0.9062812669545872, 0.9209622685150891, 0.9344722804706077, 0.9437486772168484, 0.9487027149591579, 0.9498185552187592, 0.9497105723862685, 0.9545448075185998, 0.9617684698788191, 0.9710090747107682, 0.9789383088782843, 0.9702975545756256, 0.9499977226102866, 0.9313857555691152, 0.9110443438724538, 0.8924462479473396, 0.8843812432322101, 0.8877159814322754, 0.8970755973108331, 0.9045270844069161, 0.8993174497652923, 0.882638956356522, 0.8611334823663758, 0.8515325832871878, 0.8626938028076607, 0.8745021858760732, 0.8756438986599964, 0.8812186330450035, 0.8893014600845084, 0.8933365053270356, 0.8976108912698163, 0.9022439287763816, 0.9063732074940425, 0.9103585544089995, 0.9133756373204529, 0.915218040059793, 0.9167389215069297, 0.9178664161897893, 0.916793463772099, 0.9104334703959749, 0.8969610102923625, 0.8794398983545699, 0.8652462295157394, 0.8606248920331012, 0.8617483607193223, 0.8667941744659874, 0.8742751728632395, 0.8788629641284397, 0.8780955180996398, 0.8726582136740815, 0.8608232869163714, 0.8491764702487967, 0.8478998232340591, 0.8529805096939942, 0.8526102547676554, 0.8476306090301311, 0.8457142412671218, 0.8464765739198666, 0.8509578380487717, 0.8583057306291684, 0.8648946031532602, 0.8726357762909063, 0.8819550738984164, 0.8899563314317614, 0.8931957170432046, 0.8906252115680379, 0.8807428497820006, 0.8626093671706966, 0.8440054216595896, 0.8367682282444255, 0.8355007977206106, 0.8358588282246266, 0.8359000873438316, 0.837341337906672, 0.8401344679695073, 0.8415113453391251, 0.8394873940413847, 0.8381775398018678, 0.8377687416607446, 0.8352303239815054, 0.8286650716857552, 0.8259897273584161, 0.8267535975466362, 0.8271092477922529, 0.8300862917650494, 0.8328847251471315, 0.8334195037827323, 0.8341529134244667, 0.8372856185245032, 0.8415756693402106, 0.8416779639541678, 0.8389183232961801, 0.8357094762955424, 0.8325824648975845, 0.830888941375558, 0.8280094509395097, 0.8226977926191602, 0.8190046660194084, 0.8173437596480766, 0.8166998984025333, 0.8176849146724874, 0.8162390849223692, 0.8118250002405033, 0.8109559342458291, 0.811605752348145, 0.8107205055734571, 0.8078578983975526, 0.805905756790131, 0.8048573166899318, 0.8021697644239261, 0.8025410621033174, 0.8043396911106152, 0.80617772013403, 0.8099854708379974, 0.8111178808860322, 0.8109521528820676, 0.810604835444787, 0.8078302157973883, 0.8030065239671748, 0.7998768440493922, 0.8002031860599436, 0.8015459656830464, 0.8022801117371239, 0.8025895566209336, 0.8009976394099932, 0.796909897468257, 0.7932419722184467, 0.7900634524547986, 0.7880912887515856, 0.7876145039589342, 0.7868782614033442, 0.785246483420226, 0.7843909323122046, 0.783347635222842, 0.7825851451093653, 0.7829225640742327, 0.7823746492757602, 0.7821314295525564, 0.7835617359481047, 0.7848651184969174, 0.7849939603821315, 0.7844469512334469, 0.7835397420607317, 0.7819705503437115, 0.7808474868530114, 0.7799414267345256, 0.7783073796375992, 0.7757655503047028, 0.7742460100625337, 0.7730368580611338, 0.7730601774182054, 0.7731459102313756, 0.7732016363208223, 0.7710036326220469, 0.7695471205978753, 0.7695911414198725, 0.7682587587653312, 0.7663118105940208, 0.7646010870829482, 0.7634310207834172, 0.7627968264908127, 0.7616976068732406, 0.7603241741575081, 0.7597978031279683, 0.7597206644939014, 0.7610586413696856, 0.762547862899089, 0.7621283032757984, 0.7612322564645312, 0.7598995314471827, 0.7592987108867857, 0.7585560838221451, 0.7574266543111912, 0.7559113728606949, 0.7548388694364594, 0.7534543867653245, 0.7528541717224281, 0.7522174342365513, 0.7512690624930461, 0.7506617453329167, 0.7510166279686973, 0.7501812383614269, 0.7488438277285572, 0.7474044904831457, 0.7460625974044128, 0.7451586166777585, 0.7443407791939128, 0.7436646868642668, 0.743368294863701, 0.7434273226201144, 0.7437307694803829, 0.7436995260971572, 0.7432288430168671, 0.7420460536058634, 0.7410941786025518, 0.7409979238378449, 0.7403682653582049, 0.739818579391388, 0.7395681153124208, 0.7387985935725994, 0.7380338857166782, 0.7377131709134914, 0.7370039135895862, 0.7360797052863829, 0.7349492790466655, 0.7338078546418284, 0.7328877646756945, 0.7321834630553017, 0.7317587365758971, 0.7314967070533577, 0.7308734702721713, 0.7301762803287742, 0.7293149465530941, 0.7285736070015517, 0.7279178952046477, 0.7277037368477696, 0.7278148210252474, 0.7276056701621647, 0.7266999836086818, 0.7256258835064731, 0.7253177706914882, 0.7246754926809403, 0.7240709999386562, 0.7240128196497297, 0.7234851124704469, 0.7230692050716838, 0.7227023129552484, 0.7221548881192963, 0.7216925946693372, 0.7210405526396667, 0.7202722992934193, 0.7195700978647549, 0.7189303972950564, 0.7183602997015571, 0.7177743996188006, 0.717000500590499, 0.7162466124178547, 0.7156439075487863, 0.7151325661457415, 0.714687247980993, 0.7143844932173643, 0.7141947723802528, 0.7136195637123931, 0.7130842962862802, 0.7124497027994698, 0.7121513173413774, 0.7119320957089497, 0.711502121903862, 0.7110660578383569, 0.7104680250335754, 0.7099604051880466, 0.709472290269767, 0.709124180431282, 0.7088472909433123, 0.7085246249131186, 0.708015284112709, 0.7077625238174652, 0.7073395389892835, 0.7068769061183443, 0.7063645910975325, 0.7058712846353243, 0.7052922085988642, 0.7046420922705505, 0.7043256256496293, 0.7038480698863968, 0.7035966433262432, 0.703508651500307, 0.7033583811045853, 0.7030000913479006, 0.7024968081371503, 0.7023658299392591, 0.7019593622666904, 0.7016389995978827, 0.7013480859341081, 0.7009082636488576, 0.7007172624437563, 0.7002376502961347, 0.699992484495298, 0.6996871020773557, 0.6994656591299294, 0.6991375043653474, 0.6988701812951335, 0.6986666736994509, 0.6983012491547751, 0.6980262471291135, 0.697650552636907, 0.6972539676920789, 0.6970247326488789, 0.6966422292384273, 0.6964924416601179, 0.6961910739986797, 0.6959139964691048, 0.6957529979962559, 0.6955144243927722, 0.6953545182463075, 0.6952958296143197, 0.6948686093172707, 0.6948562549461965, 0.6946215351600799, 0.6944132805988328, 0.6947329354498389, 0.6942020484658425, 0.6939759868212309, 0.6938412504035937, 0.6937141602706455, 0.6935769746805385, 0.6934062217483287, 0.6932797084535949, 0.693152556970603, 0.6930733093240299, 0.6929081598433539, 0.6927189520797588, 0.6928429385895357, 0.6925256130679929, 0.6924422662038753, 0.6923897156311157, 0.6923517620809029, 0.6922710962381989, 0.6921676771683798, 0.6921446656299842, 0.6922438042892081, 0.6921404656889575, 0.6919601192365243, 0.6920292702915485, 0.6919782143930162, 0.6921381329563252, 0.6920137812033639, 0.6922082271753086, 0.692031025442177 ] }, { "line": { "color": "#d62728", "width": 3 }, "name": "E0=8991.5", "type": "scatter", "uid": "0b593aa1-7015-4e51-8e3a-df9a7d771b79", "x": [ 8799.503495318811, 8809.588785746506, 8819.61492805748, 8829.663174433894, 8839.694403248652, 8849.747736321095, 8859.82417701776, 8869.840465353249, 8879.921993046302, 8889.942364845994, 8900.02897963994, 8910.053434993935, 8920.100998258371, 8930.171669678528, 8940.222313612707, 8950.253933011923, 8950.77056237059, 8951.24505893233, 8951.719555573367, 8952.194052293695, 8952.754821247221, 8953.18618205606, 8953.659675856557, 8954.135176066182, 8954.608670024927, 8955.170442620767, 8955.601803796582, 8956.076301165665, 8956.55180177998, 8957.11156842799, 8957.587069215671, 8958.019433925068, 8958.49393169816, 8959.055705062681, 8959.530203008942, 8960.005704201656, 8960.43806927934, 8960.912567456608, 8961.474341299649, 8961.950845985772, 8962.42534441584, 8962.85871303013, 8963.376347850546, 8963.851849687342, 8964.328354772455, 8964.803856768676, 8965.27935884453, 8965.799000444642, 8966.274502687143, 8966.750005009273, 8967.225507411034, 8967.702013062299, 8968.221655105779, 8968.697157753988, 8969.173663652218, 8969.64916645985, 8970.169812030423, 8970.64631817579, 8971.121821230041, 8971.599330706427, 8972.074833920271, 8972.551340385335, 8973.071986488054, 8973.547489948516, 8974.024999832825, 8974.50150662516, 8974.979016669908, 8975.49765687251, 8975.975167084782, 8976.452677377358, 8976.929184577099, 8977.40669503011, 8977.926338849633, 8978.403849470335, 8978.881360171343, 8979.358870952654, 8979.836381814268, 8980.312889581855, 8980.83454028227, 8981.311048217365, 8981.789562582135, 8982.26707385252, 8982.74458520321, 8983.222096634207, 8983.785881242242, 8984.264396024084, 8984.741907710666, 8985.220422653618, 8985.654797923233, 8986.132309843642, 8986.698101382693, 8987.175613478557, 8987.654128831651, 8988.088504473575, 8988.56701998051, 8989.044532390675, 8989.523048058729, 8990.088840273425, 8990.524219465751, 8991.002735383168, 8991.48125138123, 8991.959767459932, 8992.438283619278, 8992.960939726134, 8993.438452875233, 8993.917972462912, 8994.396488952267, 8994.876008701736, 8995.354525352543, 8995.834045263808, 8996.312562076067, 8996.83622207624, 8997.31574223773, 8997.794259299695, 8998.274782803659, 8998.753300027249, 8999.232820512498, 8999.712341078726, 9000.2360017059, 9000.714519260006, 9001.195043258174, 9001.674564155783, 9002.154085134373, 9002.633606193947, 9003.114130516886, 9003.592648556034, 9004.073173041275, 9004.640974521142, 9005.077359127885, 9005.55688068114, 9006.037405498797, 9006.516927214185, 9006.998455378078, 9007.477977255765, 9007.959505582634, 9008.439027622628, 9008.96369304368, 9009.443215253264, 9009.924743913414, 9010.40526947026, 9010.886798293566, 9011.366320827921, 9011.847849814207, 9012.328375696512, 9012.808901660135, 9013.289427705082, 9013.770957017503, 9014.339763623, 9014.778156093354, 9015.257679284889, 9015.740212117951, 9016.220738658974, 9016.702268468509, 9017.183798359707, 9017.665328332565, 9018.146858387088, 9018.629391711138, 9019.109918741115, 9019.59245222883, 9020.117119710816, 9020.598650181131, 9021.081183921842, 9021.562714555654, 9022.045248460205, 9022.525776068262, 9023.008310136474, 9023.489841097104, 9023.97337851891, 9025.06986513526, 9026.164345793974, 9027.260833256407, 9028.4024647514, 9029.58723394993, 9030.772003642816, 9031.957777022644, 9033.187691389878, 9034.417606289864, 9035.647521722603, 9036.922581432693, 9038.240779091346, 9039.517846356464, 9040.880185831782, 9042.20239811256, 9043.565742072607, 9044.974230559883, 9046.38472614119, 9047.83936314317, 9049.248856949082, 9050.704498617577, 9052.203278620436, 9053.703062614699, 9055.24899463792, 9056.750786643628, 9058.385002093395, 9059.886795757791, 9061.521013012501, 9063.158240818599, 9064.750325380995, 9066.475837067901, 9068.114070936384, 9069.796446811419, 9071.479826889832, 9073.208352287784, 9074.981019879944, 9076.799836165836, 9078.531374429562, 9080.350192989712, 9082.17101913559, 9084.037994153789, 9085.903967187947, 9087.770944660691, 9089.684071157459, 9091.597198943276, 9093.555472663164, 9095.515754163716, 9097.476037017646, 9099.52661070416, 9101.48890276147, 9103.541485781812, 9105.595073504577, 9107.649665931947, 9109.750407968611, 9111.89529327537, 9113.997041674817, 9116.145943078845, 9118.338987925226, 9120.53404091322, 9122.72909559818, 9124.97130357549, 9127.213513323468, 9129.501873285695, 9131.747096294957, 9134.084618115046, 9136.421138630054, 9138.760670758762, 9141.100204815188, 9143.485889470214, 9145.87257936259, 9148.306423244088, 9150.741272446407, 9153.223275801281, 9155.705281325812, 9158.19029873155, 9160.675318312222, 9163.253641290361, 9165.78682081562, 9168.366151676204, 9170.949497846703, 9173.577992285049, 9176.208495643828, 9178.840004684147, 9181.518668697847, 9184.199341731412, 9186.926166672825, 9189.60985455664, 9192.386847156942, 9195.117692962476, 9197.943850259424, 9200.72486401972, 9203.554036657548, 9206.384215368687, 9209.262553151346, 9212.141897108186, 9215.024253756786, 9217.953766471026, 9220.884285467218, 9223.863967246545, 9226.799508649998, 9229.78120316204, 9232.859214105796, 9235.845931304622, 9238.880808398044, 9241.963845540498, 9245.04889256644, 9248.228250228372, 9251.317317141373, 9254.501698178661, 9257.594785006311, 9260.78318617787, 9263.97359747904, 9267.21317287422, 9270.453755242606, 9273.698354421624, 9276.99111470764, 9280.284882092701, 9283.62881737355, 9286.974763158121, 9290.32271945353, 9293.674692836985, 9297.073824631314, 9300.476973644922, 9303.929287882347, 9307.384616187259, 9310.840951988426, 9314.346453322292, 9317.903126958183, 9321.368505217684, 9324.930204130576, 9328.494917403581, 9332.061641749433, 9335.678535570298, 9339.344595756718, 9342.968521972733, 9346.640611394825, 9350.316718786502, 9353.994837550637, 9357.724129771239, 9361.4554335083, 9365.237910976046, 9369.022400106318, 9372.811910846205, 9376.602429958382, 9380.44512654791, 9384.33799419246, 9388.233877130022, 9392.03846330908, 9395.988536157876, 9399.895471630663, 9403.851575354656, 9407.813704543412, 9411.825002274423, 9415.792159231345, 9419.857648034907, 9423.87999942387, 9427.954529801229, 9432.031072693917, 9436.11263811676, 9440.293539755176, 9444.382140424326, 9448.570077618859, 9452.712870634103, 9456.908846858942, 9461.157003184928, 9465.40817590458, 9469.664371725932, 9473.923583974067, 9478.186816012812, 9482.503232097157, 9486.871829116339, 9491.196285243612, 9495.571919133703, 9499.953579915173, 9504.38641879375, 9508.776120067316, 9513.267167754253, 9517.71407451908, 9522.21316327103, 9526.765437625592, 9531.225409329796, 9535.78672847026, 9540.30390656142, 9544.97360466218, 9549.498824373864, 9554.17756799226, 9558.762001718671, 9563.449791025992, 9568.092435881894, 9572.788267664746, 9577.491130911916, 9582.296347426462, 9587.007253533204, 9591.822520058166, 9596.54347589258, 9601.368792529494, 9606.19813094953, 9611.03349796762, 9615.924060161917, 9620.818644383307, 9625.718254054984, 9630.62288920295, 9635.58272014692, 9640.549583609143, 9645.519465988937, 9650.496380946759, 9655.527488926964, 9660.565629725801, 9665.657963956353, 9670.705153497225, 9675.859717961392, 9680.918966701527, 9686.033412761046, 9691.205063289333, 9696.380736947858, 9701.563444046225, 9706.80235274414, 9711.997120061977, 9717.246082360722, 9722.554257124058, 9727.815280203173, 9733.13551601219, 9738.459775560095, 9743.791069201898, 9749.178565832317, 9754.625276115172, 9759.974662210783, 9765.381255135533, 9770.793879078545, 9776.264713765286, 9781.741579739657, 9787.275653470428, 9792.764582223723, 9798.364906026303, 9803.865897803822, 9809.426104811932, 9815.04552739853, 9820.670981935806, 9826.252294795495, 9831.941993956363, 9837.587551519418, 9843.291322107656, 9849.055313034576, 9854.719970812133, 9860.445852498675, 9866.22894451501, 9872.019072737787, 9877.818244193044, 9883.622444961962, 9889.487870829811, 9895.3583263043, 9901.238828915499, 9907.124361216529, 9913.070115955938, 9918.970725375279, 9924.93155754717, 9930.899426806509, 9936.929526378797, 9942.912473576967, 9948.954640865239, 9955.006856200856, 9961.118292173913, 9967.184582763699, 9973.313104881045, 9979.502855401504, 9985.59126292649, 9991.742906009396, 9997.954774459384, 10004.121497374543, 10010.349449523432, 10016.58644719983, 10022.885678230085, 10029.137756622596, 10035.45006161629, 10041.772416139707, 10048.104820259958, 10054.497451741463, 10060.849955369446, 10067.252651134613, 10073.67543255442, 10080.098228505043, 10086.541110292079, 10093.034185079361, 10099.547346114761, 10106.11070070228, 10112.583748758478, 10119.117026180964, 10125.66035442965, 10132.263912645723, 10138.87752205995, 10145.50118274254, 10152.134894763829, 10158.889053186189, 10165.542868355002, 10172.196699118893, 10178.920797332663, 10185.715163495563, 10192.449329755564, 10199.263800398736, 10206.078287399188, 10212.852646359439, 10219.687238225262, 10226.531882692787, 10233.386579834701, 10240.301510705138, 10247.286711965247, 10254.171567956879, 10261.066476915244, 10268.081838151676, 10275.047035127684, 10282.082503600239, 10289.128025877266, 10296.183602033589, 10303.24923214412, 10310.264697621067, 10317.460836873304, 10324.556629422515, 10331.66247622688, 10338.89881589846, 10346.014735171826, 10353.271184618687, 10360.417249780401, 10367.683772687953, 10374.970387539948, 10382.257021094545, 10389.563746799873, 10396.870491311194, 10404.19732818011, 10411.534220760397, 10418.88116912998, 10426.288357635878, 10433.655417948468, 10441.092755721664, 10448.54014991177, 10456.05782235207, 10463.525329771799, 10471.00289384627, 10478.550736886222, 10486.108637059233, 10493.68663153774, 10501.264646244595, 10508.862755471504, 10516.460885034296, 10524.139332530702, 10531.767614660803, 10539.466178065995, 10547.174799629627, 10554.903516742184, 10562.632254895783, 10570.381088817143, 10578.139981279977, 10585.898894948485, 10593.748166834952, 10601.657685266893, 10609.396588208887, 10617.276037651971, 10625.165546521863, 10633.12534041283, 10641.034968555672, 10648.96469401845, 10656.96470531644, 10664.904512972678, 10672.934682367068, 10680.964874475918, 10689.015164864875, 10697.075515893548, 10705.145927647567, 10713.226400212665, 10721.44742633883, 10729.49783152595, 10737.678752952399, 10745.80950803565, 10754.020590506008, 10762.301963018468, 10770.472940907574, 10778.724247113556, 10786.98561543145, 10795.32731305909, 10803.608806069216, 10811.97062897716, 10820.282285279402, 10828.67427212827, 10837.076322090923, 10845.428205240833, 10853.860419885954, 10862.302697972456, 10870.885539123456, 10879.287714016782, 10887.770221661363, 10896.333062777216, 10904.835698444731, 10913.488938222084, 10921.961509453009, 10930.644917343281, 10939.207887781648, 10947.851193631777, 10956.504564526049, 10965.238271872178, 10973.86157939491, 10982.555184855928, 10991.339166604534, 11000.062942235601, 11008.867055983692, 11017.691274925513, 11026.395053584067, 11035.249444029996, 11044.11390113086, 11052.988424980877, 11061.883054775608, 11070.787751571168, 11079.712554626494, 11088.647424935187, 11097.602401819773, 11106.567446211126, 11115.542558204517, 11124.60805247891, 11133.623378927014, 11142.779285423987, 11151.75451190092, 11160.9405952284, 11170.01627395763, 11179.102021216679, 11188.268113938746, 11197.524592482909, 11206.650586233709, 11215.86696613201, 11225.163693379292, 11234.339974534829, 11243.606643044932, 11252.943620240676, 11262.170190237752, 11271.547387858474, 11280.874417074, 11290.281796203495, 11299.699246448752, 11309.006288204895, 11318.463960598816, 11327.861424248453, 11337.349279683347, 11346.857247045367, 11356.305005395505, 11365.84315677168, 11375.391380403904, 11384.959716637495 ], "y": [ 0.006099561801475914, 9.121497982804132e-05, 5.468058201668235e-05, -2.7201313452935966e-05, -6.941160081994733e-05, -0.00010059036760923766, -7.515424682835094e-05, -1.339585535803716e-05, 5.273792862081695e-05, 8.711989360137464e-05, 0.00017072113569008163, 0.0004115950458938721, 0.000561263095883639, 0.0009348124798592847, 0.0012947643171494793, 0.0021678808032216273, 0.0022611558704725364, 0.0023879207193445245, 0.0023766092917275033, 0.0023725192621727622, 0.0025083324957700715, 0.002596728457538076, 0.0026462838959072166, 0.002719490434170256, 0.00278044350775007, 0.0028466460816932077, 0.002893236412282186, 0.0029409142531597265, 0.0029657028188941493, 0.003031963141991475, 0.0030824616020743065, 0.0032312753263934704, 0.0032935265217702552, 0.0033466785676240784, 0.0034417740341423485, 0.003518286953253211, 0.003618204172131826, 0.003691548538221832, 0.0037894645488575577, 0.0038879391138853074, 0.003962371090448365, 0.0040838149413271515, 0.004181999166657266, 0.004229844206642737, 0.004403534328208882, 0.004465430643815907, 0.0046303839195448275, 0.0047328172976952255, 0.004797280300531415, 0.004937611327382278, 0.005073113615884189, 0.005231752061180497, 0.005365942209705911, 0.005503967683174689, 0.005693318802501169, 0.005849745781367775, 0.006005593066791347, 0.0061400878404700016, 0.006348108596224188, 0.006492210329189507, 0.0067077570159345, 0.006989202618542135, 0.007226573321612011, 0.007393963016457601, 0.007675010068371107, 0.007874758424415677, 0.008182602925409916, 0.008389797290692846, 0.008714868736744399, 0.009034632919313338, 0.009431446769716297, 0.009782576142631028, 0.010250298780710742, 0.01063618652471933, 0.011077713724421324, 0.011506364260828185, 0.01205903986502245, 0.012577268747615767, 0.013251526244374812, 0.013968778174404, 0.014838186996250623, 0.015396909505580298, 0.016183453565177804, 0.017380311314188986, 0.01868187916544051, 0.019738260950590166, 0.02123854728016106, 0.023382312534169315, 0.024357308618397116, 0.02645618661433456, 0.029765371467061156, 0.03331213647648775, 0.03785270441447502, 0.04146041613555418, 0.051007315537800836, 0.05791689523538759, 0.06835037149998752, 0.09645582272114846, 0.13109449948055954, 0.19163218227725748, 0.291358910037144, 0.38038668270034653, 0.44203671930177574, 0.4717631532075253, 0.4821189850894681, 0.47868718457221315, 0.46968470441132887, 0.4590645842444664, 0.4537084667165455, 0.4600829542305232, 0.4715652711343757, 0.4914321198314303, 0.5145706528666153, 0.533320463688335, 0.5491229796929266, 0.5660143723878066, 0.5931333657368766, 0.6100084380240094, 0.6415169367588665, 0.6720971975624215, 0.7170701794404782, 0.7511312695542127, 0.7993036350356275, 0.8462192191514495, 0.8807773159884708, 0.9149204150811487, 0.9579745107067331, 1.005040424860033, 1.0238179397294473, 1.0437964661465482, 1.0472708662616124, 1.0410066315828579, 1.0314544790967999, 1.0148219068015438, 0.9982715397073829, 0.9809681620145565, 0.9620568383576359, 0.9475671893201589, 0.9394262514778627, 0.9337729996298422, 0.9314639798608413, 0.9352200646222621, 0.9441150616781604, 0.9585916281456077, 0.986457766860798, 1.0013062780082969, 1.0253988028586218, 1.0474806881959617, 1.0640104254945961, 1.0717021632731376, 1.0704669908184485, 1.0622740750330495, 1.043916622800525, 1.03394119331949, 1.0160235493516814, 1.0066092995461415, 0.9821125785307698, 0.971155526847618, 0.9575700521932953, 0.9439955908591313, 0.9353900488563577, 0.9264397440825182, 0.9145737048961645, 0.9063685945174959, 0.8975727371088275, 0.8916697683624875, 0.8881392045204346, 0.8873845081599219, 0.8921483213148446, 0.8985926424774291, 0.9095968604249631, 0.921998706298006, 0.9452647292812246, 0.9673348742024391, 0.9847683197419573, 1.010144538526483, 1.0349587077355704, 1.054355535915089, 1.0696813182187748, 1.0716400628030618, 1.0586288694525021, 1.0324987392816334, 1.0070523134178382, 0.9768466636738559, 0.9539510392490925, 0.9389023452896911, 0.9286412030089095, 0.9267095467594987, 0.9311819734306317, 0.9417102683717263, 0.9515035850819976, 0.9544649756281556, 0.9505228630264851, 0.9462270129764181, 0.9493914360864185, 0.961435484480424, 0.9738647246927404, 0.9770354719379496, 0.9749656441395267, 0.9763946693506075, 0.9860501030019376, 0.9997971384045121, 1.0123294890214583, 1.0205552361304266, 1.028792403676955, 1.0386064926747622, 1.0534908117976334, 1.068911500348727, 1.0808756061966025, 1.0816554961303209, 1.0709318491369872, 1.0485444759749338, 1.0183493106794848, 0.9934152127733833, 0.9741886024723375, 0.9573472234419211, 0.9483712936481533, 0.946854455489282, 0.9377001527380614, 0.9146583871973263, 0.8935291934176659, 0.8783699299932003, 0.8752709327765168, 0.8846267303134088, 0.9008518611768543, 0.9219458934398264, 0.9453581493642926, 0.967581487846662, 0.982222696396693, 0.9856568433665914, 0.9789929121141866, 0.9708760910086515, 0.9686779774936013, 0.9866816860020484, 1.0299469405085457, 1.0843835854867931, 1.1122534874565009, 1.1138573479183962, 1.1015277431196635, 1.0893184155431332, 1.0751596357299897, 1.0488467493979499, 1.0135064629946342, 0.970248758853368, 0.9232594232916005, 0.8771089130542854, 0.8375508514153162, 0.8131350418066695, 0.8106060082387234, 0.8321446072888227, 0.8699785079937019, 0.9144777166615767, 0.9578937045125206, 0.9866969460338278, 1.0023791113721372, 0.9999938983594414, 0.9836452083666659, 0.9683762272900938, 0.9611313460431031, 0.9584282011651588, 0.9467768446740905, 0.929523793821685, 0.9262904395310392, 0.9470013818881013, 0.9898518347771084, 1.0326128840146773, 1.04759541305345, 1.0495105982981647, 1.0432444171674593, 1.0251423235455273, 0.9928014836884685, 0.9525024301725136, 0.9147975154001728, 0.8862524310250259, 0.8692117362600537, 0.8618935712644662, 0.8652182258561854, 0.8793579774471044, 0.8898961562895291, 0.890157338937056, 0.8929441373104512, 0.8970375209128544, 0.8990218856566305, 0.9062970968047181, 0.9209783503785676, 0.9344886049754909, 0.9437652026622269, 0.9487193987528164, 0.9498353597066062, 0.9497274876438314, 0.9545618842537196, 0.9617857346344483, 0.9710265487603995, 0.9789559799452655, 0.9703152585285838, 0.9500153448310947, 0.9314033107929001, 0.9110618189670057, 0.8924636614947413, 0.8843987022097612, 0.8877336032657709, 0.8970934453567719, 0.9045451389268625, 0.8993355861207529, 0.8826570603108908, 0.8611515064314851, 0.8515506501466099, 0.8627121224397541, 0.8745207677679817, 0.8756626366348415, 0.8812375729526506, 0.8893206282573574, 0.8933558641489665, 0.897630446220054, 0.9022636847024359, 0.9063931567925319, 0.9103787025692195, 0.9133959740011978, 0.9152385565827422, 0.9167596160418671, 0.9178872879795624, 0.9168144900226616, 0.9104546032198596, 0.8969821775271918, 0.8794610626924344, 0.8652674257448882, 0.860646217899958, 0.861769879581353, 0.8668159230210724, 0.8742971812377489, 0.8788852025007982, 0.8781179362464381, 0.8726807682712453, 0.8608459150715372, 0.8491991754193493, 0.8479227112251112, 0.8530036459893736, 0.8526335883564957, 0.8476540973715642, 0.8457379144975141, 0.846500462810212, 0.850981981716637, 0.8583301618062988, 0.864919313555489, 0.8726607839000418, 0.8819803939032665, 0.889981954624462, 0.8932215996718025, 0.8906512920676875, 0.8807690611963551, 0.8626356257758401, 0.8440317320599182, 0.8367946989373836, 0.8355274982630537, 0.8358857716476137, 0.8359272774603597, 0.8373687880293355, 0.8401621961308771, 0.841539339297854, 0.8395156267545522, 0.8382060150236541, 0.837797475814464, 0.8352592921531311, 0.8286942408027068, 0.8260191384992162, 0.82678328737279, 0.8271392170470482, 0.8301165688164375, 0.8329153102557939, 0.8334503762384469, 0.8341840802056896, 0.8373171059968688, 0.8416074911467063, 0.841710080204657, 0.8389507102530187, 0.8357421317303341, 0.8326153947740326, 0.8309221592932652, 0.8280429533995793, 0.8227315512641629, 0.8190387024321057, 0.8173780997603892, 0.8167345546734396, 0.8177199059102327, 0.8162743922130236, 0.8118605929384592, 0.8109918534807277, 0.8116420192523878, 0.8107571037270381, 0.807894813814943, 0.8059430009857508, 0.80489490131367, 0.8022076788183812, 0.802579343142716, 0.8043783488027101, 0.8062167608573139, 0.8100249169779784, 0.8111577116747961, 0.8109923578856816, 0.8106454190130106, 0.8078711526284005, 0.803047803475827, 0.7999184785946497, 0.800245216842778, 0.8015884094342363, 0.8023229650985042, 0.8026328144736328, 0.801041292788436, 0.7969539204339958, 0.7932863756206243, 0.7901082480643665, 0.7881364837856621, 0.7876601204016794, 0.7869243033609375, 0.7852929448846679, 0.7844378242554926, 0.7833949586241505, 0.7826329099833945, 0.7829707844159259, 0.7824233195039376, 0.7821805556794483, 0.7836113420978178, 0.7849152024384723, 0.7850445176528742, 0.7844979782202915, 0.7835912427671168, 0.7820225169884647, 0.7808999313757504, 0.7799943548350325, 0.7783607915582468, 0.7758194359455133, 0.7743003875901061, 0.7730917384808685, 0.7731155674569825, 0.7732018185883006, 0.7732580705961343, 0.7710605691432192, 0.7696045748816095, 0.7696491322063824, 0.7683172804725253, 0.7663708558624616, 0.7646606664501368, 0.763491143516168, 0.7628575018218774, 0.7617588384447076, 0.7603859586394801, 0.7598601565358646, 0.7597835961905103, 0.7611221686396653, 0.762611992145486, 0.7621930229702158, 0.7612975667811481, 0.759965435783296, 0.7593652113640834, 0.7586231877550085, 0.7574943619077826, 0.7559796891062929, 0.7549078029631522, 0.7535239386267765, 0.7529243540157277, 0.7522882614895788, 0.7513405251880454, 0.7507338501723203, 0.7510893948002223, 0.7502546654828078, 0.7489179075286233, 0.7474792332893988, 0.7461380076283517, 0.7452346981145872, 0.7444175421936909, 0.7437421373249887, 0.743446441092531, 0.7435061778487643, 0.743810346952081, 0.7437798156395342, 0.7433098447575428, 0.7421277765088636, 0.7411766233448657, 0.7410811100438048, 0.740452192385843, 0.7399032527852356, 0.7396535428384583, 0.7388847682845283, 0.738120830881383, 0.737800883891544, 0.7370923952909161, 0.736168971628773, 0.7350393185146421, 0.733898686468569, 0.732979382660327, 0.7322758866437781, 0.7318519747435628, 0.7315907653601448, 0.7309683512901775, 0.7302719872588259, 0.7294114839932159, 0.7286709812832131, 0.728016112329924, 0.7278028122590368, 0.727914757352091, 0.7277064764504001, 0.7268016581231662, 0.725728436986032, 0.7254212090367016, 0.724779817852142, 0.724176224899155, 0.7241189553181419, 0.7235921607811241, 0.7231771713925597, 0.7228112044824179, 0.7222647073022707, 0.7218033565598939, 0.7211522534507865, 0.7203849509236191, 0.71968370661613, 0.7190449707151934, 0.7184758429043467, 0.7178909194099696, 0.7171180007966405, 0.7163650976884778, 0.7157633954884185, 0.7152530700181238, 0.7148087508960973, 0.7145070193782038, 0.7143183288276393, 0.7137441604215904, 0.7132099317009458, 0.7125763832106421, 0.7122790601556386, 0.7120608983777964, 0.7116319991290199, 0.7111970143232377, 0.7106000666404171, 0.7100935389845101, 0.709606522628113, 0.7092595189540329, 0.7089837605966733, 0.7086622065392877, 0.7081539988165712, 0.707902372097329, 0.7074805353109509, 0.707019064989593, 0.7065079014223715, 0.7060157629783462, 0.7054378605640794, 0.7047889337248936, 0.7044736564001453, 0.7039973050174498, 0.7037470829747255, 0.7036603141515169, 0.7035112728090193, 0.7031542078313169, 0.7026521649809185, 0.7025224376346736, 0.7021172442854001, 0.7017981351646088, 0.7015084926486531, 0.701069957467709, 0.7008802421462732, 0.7004019412698302, 0.7001580669771666, 0.6998540131613835, 0.6996338866424846, 0.69930706509448, 0.6990410829347962, 0.6988389349714403, 0.6984748567431905, 0.6982012185066311, 0.6978269066644477, 0.697431700295223, 0.6972038639403115, 0.696822766618147, 0.6966743738573224, 0.6963744293899397, 0.6960987827186378, 0.6959392237114743, 0.6957020978927078, 0.6955436478147908, 0.6954864254351926, 0.6950606752283395, 0.6950498043530402, 0.6948165734241484, 0.6946098156048943, 0.6949309936399996, 0.6944016188892277, 0.6941771023503961, 0.6940438874289263, 0.69391836080425, 0.6937827259323042, 0.6936135311936881, 0.6934885965076228, 0.693363045500804, 0.6932853824256678, 0.6931218385944066, 0.692934256603805, 0.6930598572070068, 0.6927441635345006, 0.6926624696907372, 0.6926115592141331, 0.6925752791303368, 0.6924962838088218, 0.6923945559349418, 0.6923732447569881, 0.6924740715008222, 0.6923724528038473, 0.6921938210647098, 0.6922647124678556, 0.692215406062479, 0.6923770717700464, 0.6922544876865279, 0.6924507131299074, 0.6922752976250847 ] } ], "layout": { "height": 500, "hovermode": "closest", "legend": { "bgcolor": "#F2F2F2", "borderwidth": 0.5 }, "plot_bgcolor": "#FDFDFF", "showlegend": true, "template": { "data": { "bar": [ { "error_x": { "color": "#2a3f5f" }, "error_y": { "color": "#2a3f5f" }, "marker": { "line": { "color": "#E5ECF6", "width": 0.5 }, "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "bar" } ], "barpolar": [ { "marker": { "line": { "color": "#E5ECF6", "width": 0.5 }, "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "barpolar" } ], "carpet": [ { "aaxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "baxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "type": "carpet" } ], "choropleth": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "choropleth" } ], "contour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "contour" } ], "contourcarpet": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "contourcarpet" } ], "heatmap": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmap" } ], "heatmapgl": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmapgl" } ], "histogram": [ { "marker": { "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "histogram" } ], "histogram2d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2d" } ], "histogram2dcontour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2dcontour" } ], "mesh3d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "mesh3d" } ], "parcoords": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "parcoords" } ], "pie": [ { "automargin": true, "type": "pie" } ], "scatter": [ { "fillpattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 }, "type": "scatter" } ], "scatter3d": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatter3d" } ], "scattercarpet": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattercarpet" } ], "scattergeo": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergeo" } ], "scattergl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergl" } ], "scattermapbox": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattermapbox" } ], "scatterpolar": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolar" } ], "scatterpolargl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolargl" } ], "scatterternary": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterternary" } ], "surface": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "surface" } ], "table": [ { "cells": { "fill": { "color": "#EBF0F8" }, "line": { "color": "white" } }, "header": { "fill": { "color": "#C8D4E3" }, "line": { "color": "white" } }, "type": "table" } ] }, "layout": { "annotationdefaults": { "arrowcolor": "#2a3f5f", "arrowhead": 0, "arrowwidth": 1 }, "autotypenumbers": "strict", "coloraxis": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "colorscale": { "diverging": [ [ 0, "#8e0152" ], [ 0.1, "#c51b7d" ], [ 0.2, "#de77ae" ], [ 0.3, "#f1b6da" ], [ 0.4, "#fde0ef" ], [ 0.5, "#f7f7f7" ], [ 0.6, "#e6f5d0" ], [ 0.7, "#b8e186" ], [ 0.8, "#7fbc41" ], [ 0.9, "#4d9221" ], [ 1, "#276419" ] ], "sequential": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "sequentialminus": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ] }, "colorway": [ "#636efa", "#EF553B", "#00cc96", "#ab63fa", "#FFA15A", "#19d3f3", "#FF6692", "#B6E880", "#FF97FF", "#FECB52" ], "font": { "color": "#2a3f5f" }, "geo": { "bgcolor": "white", "lakecolor": "white", "landcolor": "#E5ECF6", "showlakes": true, "showland": true, "subunitcolor": "white" }, "hoverlabel": { "align": "left" }, "hovermode": "closest", "mapbox": { "style": "light" }, "paper_bgcolor": "white", "plot_bgcolor": "#E5ECF6", "polar": { "angularaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "radialaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "scene": { "xaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "yaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "zaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" } }, "shapedefaults": { "line": { "color": "#2a3f5f" } }, "ternary": { "aaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "baxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "caxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "title": { "x": 0.05 }, "xaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 }, "yaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 } } }, "title": {}, "width": 650, "xaxis": { "color": "#004", "gridcolor": "#D8D8D8", "range": [ 8950, 9100 ], "showgrid": true, "title": {}, "zerolinecolor": "#DDD" }, "yaxis": { "color": "#004", "gridcolor": "#D8D8D8", "showgrid": true, "title": {}, "zerolinecolor": "#DDD" } } }, "text/html": [ "
" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# That is, an angular error of 0.02 degrees gives an energy shift of 14 eV.\n", "# Let's plot those XANES spectra together\n", "\n", "multi_plot([dict(xdata=cudat.energy, ydata=cudat.norm, label='E0=8977.6'),\n", " dict(xdata=new_dat.energy, ydata=new_dat.norm, label='E0=8991.5',\n", " xmin=8950, xmax=9100, xlabel='Energy (eV)', ylabel='$\\mu(E)$')])\n", "\n", "# that is a very big energy shift for XANES work." ] }, { "cell_type": "code", "execution_count": 8, "id": "638b0a19", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Result: E0 = 8991.481 eV, R= 2.5455 +/- 0.0017 Ang\n" ] } ], "source": [ "# and now we can fit that shifted data:\n", "\n", "out = fit_cu_1stshell(new_dat, full_report=False, with_plot=False)\n", " \n", "# and we see that the 1st shell R changed from R= 2.5494 to 2.5455 (+/- 0.0017 Ang).\n", "# a shift in R of 0.004 Ang. \n", "#\n", "# Yes, that shift is outside of the error bars for this very good dataset. \n", "# Typical uncertainties and expected accuracies of EXAFS results are 0.01 Ang. " ] }, { "cell_type": "code", "execution_count": 9, "id": "6e00298d", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "3000.0 -> 3001.2 eV\n", "5000.0 -> 5004.1 eV\n", "7000.0 -> 7008.3 eV\n", "9000.0 -> 9014.0 eV\n", "11000.0 -> 11021.1 eV\n", "20000.0 -> 20070.5 eV\n" ] } ], "source": [ "# If we had an angular offset that big, how would we know?\n", "\n", "# We would change the monochromator energy to another edge. \n", "# Let's run our `energy_shift` function with different energies:\n", "\n", "etest = np.array([3000, 5000, 7000, 9000, 11000, 20000])\n", "enew = energy_shift(etest, dspace, theta_off=0.02)\n", "for e_in, e_out in zip(etest, enew):\n", " print(f'{e_in:.1f} -> {e_out:.1f} eV')\n", " " ] }, { "cell_type": "code", "execution_count": 10, "id": "aa39c4f4", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "3000.0 -> 3009.6 eV\n", "5000.0 -> 5016.0 eV\n", "7000.0 -> 7022.4 eV\n", "9000.0 -> 9028.8 eV\n", "11000.0 -> 11035.2 eV\n", "20000.0 -> 20064.0 eV\n" ] } ], "source": [ "# if we thought a shift of 14 eV at Cu K was big (it is!), then a shift \n", "# of 70 eV at Mo K edge is enormous.\n", "# For completeness, what if we had a bad value for $d$?\n", "enew = energy_shift(etest, dspace, dspace_off=0.01)\n", "for e_in, e_out in zip(etest, enew):\n", " print(f'{e_in:.1f} -> {e_out:.1f} eV')" ] }, { "cell_type": "code", "execution_count": 11, "id": "2f1fc8b8", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "theta_off=0.0000, dspace_off=0.0000: E0= 8977.580 to 8977.580\n", "Result: E0 = 8977.580 eV, R= 2.5494 +/- 0.0017 Ang\n", "theta_off=0.0010, dspace_off=0.0000: E0= 8977.580 to 8978.274\n", "Result: E0 = 8978.274 eV, R= 2.5492 +/- 0.0017 Ang\n", "theta_off=0.0020, dspace_off=0.0000: E0= 8977.580 to 8978.968\n", "Result: E0 = 8978.968 eV, R= 2.5490 +/- 0.0017 Ang\n", "theta_off=0.0050, dspace_off=0.0000: E0= 8977.580 to 8981.051\n", "Result: E0 = 8981.051 eV, R= 2.5484 +/- 0.0017 Ang\n", "theta_off=0.0100, dspace_off=0.0000: E0= 8977.580 to 8984.525\n", "Result: E0 = 8984.525 eV, R= 2.5474 +/- 0.0017 Ang\n", "theta_off=0.0200, dspace_off=0.0000: E0= 8977.580 to 8991.481\n", "Result: E0 = 8991.481 eV, R= 2.5455 +/- 0.0017 Ang\n", "theta_off=0.0300, dspace_off=0.0000: E0= 8977.580 to 8998.448\n", "Result: E0 = 8998.448 eV, R= 2.5436 +/- 0.0017 Ang\n", "theta_off=0.0400, dspace_off=0.0000: E0= 8977.580 to 9005.427\n", "Result: E0 = 9005.427 eV, R= 2.5417 +/- 0.0016 Ang\n", "theta_off=0.0500, dspace_off=0.0000: E0= 8977.580 to 9012.416\n", "Result: E0 = 9012.416 eV, R= 2.5397 +/- 0.0016 Ang\n", "theta_off=0.0600, dspace_off=0.0000: E0= 8977.580 to 9019.417\n", "Result: E0 = 9019.417 eV, R= 2.5378 +/- 0.0016 Ang\n", "theta_off=0.0700, dspace_off=0.0000: E0= 8977.580 to 9026.428\n", "Result: E0 = 9026.428 eV, R= 2.5359 +/- 0.0016 Ang\n", "theta_off=0.0800, dspace_off=0.0000: E0= 8977.580 to 9033.451\n", "Result: E0 = 9033.451 eV, R= 2.5340 +/- 0.0016 Ang\n", "theta_off=0.0900, dspace_off=0.0000: E0= 8977.580 to 9040.485\n", "Result: E0 = 9040.485 eV, R= 2.5320 +/- 0.0016 Ang\n", "theta_off=0.1000, dspace_off=0.0000: E0= 8977.580 to 9047.531\n", "Result: E0 = 9047.531 eV, R= 2.5301 +/- 0.0016 Ang\n" ] }, { "data": { "application/vnd.plotly.v1+json": { "config": { "plotlyServerURL": "https://plot.ly" }, "data": [ { "line": { "color": "#1f77b4", "width": 0.0002 }, "name": "deltaR", "type": "scatter", "uid": "886f9140-3c3e-4e07-b522-35e61c23380d", "x": [ 0, 0.6940156509335793, 1.3881413478793547, 3.471178974290524, 6.945111695282321, 13.901251381230395, 20.868445116762814, 27.846719042994664, 34.83609938354857, 41.83661244487848, 48.84828461660982, 55.871142371859605, 62.90521226757119, 69.95052094483435 ], "y": [ 2.5636612778514017e-06, -0.00018962915238568824, -0.00038182027436379985, -0.0009583867707440507, -0.001919327288730371, -0.003841173834707231, -0.005762967591438366, -0.0076847157702999895, -0.009606478155146337, -0.011549873432441118, -0.013471706614229284, -0.015393507787763513, -0.0173154616055411, -0.01923714343375289 ] } ], "layout": { "height": 500, "hovermode": "closest", "legend": { "bgcolor": "#F2F2F2", "borderwidth": 0.5 }, "plot_bgcolor": "#FDFDFF", "showlegend": true, "template": { "data": { "bar": [ { "error_x": { "color": "#2a3f5f" }, "error_y": { "color": "#2a3f5f" }, "marker": { "line": { "color": "#E5ECF6", "width": 0.5 }, "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "bar" } ], "barpolar": [ { "marker": { "line": { "color": "#E5ECF6", "width": 0.5 }, "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "barpolar" } ], "carpet": [ { "aaxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "baxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "type": "carpet" } ], "choropleth": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "choropleth" } ], "contour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "contour" } ], "contourcarpet": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "contourcarpet" } ], "heatmap": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmap" } ], "heatmapgl": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmapgl" } ], "histogram": [ { "marker": { "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "histogram" } ], "histogram2d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2d" } ], "histogram2dcontour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2dcontour" } ], "mesh3d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "mesh3d" } ], "parcoords": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "parcoords" } ], "pie": [ { "automargin": true, "type": "pie" } ], "scatter": [ { "fillpattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 }, "type": "scatter" } ], "scatter3d": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatter3d" } ], "scattercarpet": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattercarpet" } ], "scattergeo": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergeo" } ], "scattergl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergl" } ], "scattermapbox": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattermapbox" } ], "scatterpolar": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolar" } ], "scatterpolargl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolargl" } ], "scatterternary": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterternary" } ], "surface": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "surface" } ], "table": [ { "cells": { "fill": { "color": "#EBF0F8" }, "line": { "color": "white" } }, "header": { "fill": { "color": "#C8D4E3" }, "line": { "color": "white" } }, "type": "table" } ] }, "layout": { "annotationdefaults": { "arrowcolor": "#2a3f5f", "arrowhead": 0, "arrowwidth": 1 }, "autotypenumbers": "strict", "coloraxis": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "colorscale": { "diverging": [ [ 0, "#8e0152" ], [ 0.1, "#c51b7d" ], [ 0.2, "#de77ae" ], [ 0.3, "#f1b6da" ], [ 0.4, "#fde0ef" ], [ 0.5, "#f7f7f7" ], [ 0.6, "#e6f5d0" ], [ 0.7, "#b8e186" ], [ 0.8, "#7fbc41" ], [ 0.9, "#4d9221" ], [ 1, "#276419" ] ], "sequential": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "sequentialminus": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ] }, "colorway": [ "#636efa", "#EF553B", "#00cc96", "#ab63fa", "#FFA15A", "#19d3f3", "#FF6692", "#B6E880", "#FF97FF", "#FECB52" ], "font": { "color": "#2a3f5f" }, "geo": { "bgcolor": "white", "lakecolor": "white", "landcolor": "#E5ECF6", "showlakes": true, "showland": true, "subunitcolor": "white" }, "hoverlabel": { "align": "left" }, "hovermode": "closest", "mapbox": { "style": "light" }, "paper_bgcolor": "white", "plot_bgcolor": "#E5ECF6", "polar": { "angularaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "radialaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "scene": { "xaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "yaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "zaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" } }, "shapedefaults": { "line": { "color": "#2a3f5f" } }, "ternary": { "aaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "baxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "caxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "title": { "x": 0.05 }, "xaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 }, "yaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 } } }, "title": {}, "width": 650, "xaxis": { "color": "#004", "gridcolor": "#D8D8D8", "showgrid": true, "title": { "text": "E0 (eV)" }, "zerolinecolor": "#DDD" }, "yaxis": { "color": "#004", "gridcolor": "#D8D8D8", "showgrid": true, "title": { "text": "$\\Delta R1$ (Ang)" }, "zerolinecolor": "#DDD" } } }, "text/html": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# We can see that an error in angular offset is more pronounced at high energy, \n", "# or lower angle. On the other hand, a bad value for d-spacing gives a slightly less \n", "# energy-dependent error (less of a \"chromatic aberration\"). \n", "#\n", "# When calibrating a monochromator, it is important to use a wide range of \n", "# energies to set the d-spacing and angular offset.\n", "\n", "# But: how far off do you have to be to see real errors in EXAFS results?\n", "\n", "e0val = []\n", "delr1 = []\n", "\n", "for theta_off in (0.000, 0.001, 0.002, 0.005, 0.010, 0.020, 0.03, \n", " 0.04, 0.050, 0.06, 0.07, 0.08, 0.09, 0.10):\n", " new_dat = change_calib(cudat, dspace, theta_off=theta_off)\n", "\n", " out = fit_cu_1stshell(new_dat, full_report=False, with_plot=False)\n", " \n", " e0val.append(new_dat.e0-cudat.e0) \n", " delr1.append(out.params['del_r'].value -0.00155)\n", " \n", "f = plot(e0val, delr1, marker='o', linewidth=0.0002, xlabel='E0 (eV)', ylabel='$\\Delta R1$ (Ang)', label='deltaR')" ] }, { "cell_type": "code", "execution_count": 12, "id": "1353430a", "metadata": {}, "outputs": [], "source": [ "# To get an error in R of 0.01 Ang, E0 needs to shift by about 35 eV at the Cu K edge. \n", "# That is more than double the one shown in the Cu XANES plot above. That is an \n", "# absurdly large energy shift. \n", "\n", "# If you are expecting to get absolute measures that are accurate to 0.005 Ang,\n", "# then making sure your energy calibration is good to 10 eV is important. \n", "# Claiming such absolute accuracies is rare in the EXAFS literature, and do \n", "# require a very careful analysis.\n", "\n", "# Energy errors or drifts of 5 eV in the 5 to 15 keV range will cause significant \n", "# problems for XANES work and would need to be addressed. XANES Analysis is the \n", "# reason to work toward good energy calibration and reliability.\n", "#\n", "# Errors of 5 eV in the 5 to 15keV range will not impact EXAFS results at all. \n", "#\n", "# If someone tells you otherwise, point them to this notebook.\n" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.12.4" } }, "nbformat": 4, "nbformat_minor": 5 }