{ "cells": [ { "cell_type": "markdown", "id": "2980d153", "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 a well-calibrated energy, typically coming from a double-cystal monochromator. While EXAFS can be measured in modes without a double-crystal monochromator, that is sort of unusual, and we are not discussing that situation here. A double-cystal 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 used0 are typically very well known, but there may be thermally-induced changes. The values of $d$ and an angular offset ($\\theta_0$ may need to be adjusted to calibrate a monochromator.\n", "This is typically done by scanning across a known (or at least reproducible) energies, say the $K$-edges of metal foils, and adjusting $d$ or $\\theta_0$ until the energy is correct. \n", "\n", "\n", "Most commonly, the XANES of some metal foil is measured and the maximumn of the first derivative is set to the tabulated edge energy for the metal. This method will give a bit of a variation as \n", "the energy resolution changes. These days, most good XAFS beamlines have high enough energy \n", "resolution that this is not much of a problem, but it can be a problem for older data. It should \n", "also be said that those tabulated values are not always accurate to better than 1 or 2 eV. In fact, \n", "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 angular offset that \n", "stays calibrated across the edges of many edges. That is, if the energy needs to be recalibrated\n", "at every edge, the energy will be drifting between those edges, and the reason for any need to\n", "recalibrate should be investigated. If the beamline is well collimated and the angle of the beam \n", "incident on the monochomator is stable, it is certainly possible to have a monochromator set up t\n", "that stays calibrated for weeks and over an energy range of 10 keV or more. \n", "\n", "\n", "For any XAS analysis, the energy scale is critical. Because of the variations above in calibrating \n", "energy, small energy shifts between beamlines or even runs at the same beamline but different months \n", "are not uncommon. If the shift is relatively small, say 1 to 5 eV at 10 keV, simply adding a constant \n", "energy offset is satisfactory.\n", " \n", "You may hear people say that poor calibration leads to large errors in XAFS results or even leads \n", "\"the wrong results\". You should know that when people say things are Wrong is describing scientific\n", "measurements and analysis, they are almost certainly wrong. \n", "\n", "Here, we'll look at these effects\n" ] }, { "cell_type": "code", "execution_count": 2, "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", "\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": 3, "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": 4, "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": "43c3cde9-6aaa-49b4-87ed-a0764967908c", "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": 4, "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 linear." ] }, { "cell_type": "code", "execution_count": 5, "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": "9d1eb895-98c7-46a7-a191-7bee6bd67d01", "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.006259474761857858, 0.00020884413456537628, 0.00013367051997887315, 1.3030948957645212e-05, -6.785601444994545e-05, -0.0001377976297546854, -0.00015118830971820596, -0.00012801132607241026, -0.00010071531578492878, -0.00010496081280303482, -6.022165904272004e-05, 0.00014211281173055691, 0.0002530926523955572, 0.0005879856706472706, 0.0009093433105783736, 0.0017442259738558153, 0.001835559497609198, 0.001960564689123814, 0.001947414417208297, 0.001941489667325108, 0.0020752152178703845, 0.002161996036250319, 0.0022097513362068725, 0.002281163532364061, 0.002340322969503525, 0.0024043939727043514, 0.0024493451123328335, 0.0024952177794254496, 0.0025181841561758346, 0.0025823206096061294, 0.002631011585879763, 0.002778240780731134, 0.0028386950736110812, 0.0028897079044650628, 0.0029830252644179817, 0.0030577455309030387, 0.003156050106694802, 0.00322760384606963, 0.003323406208073335, 0.003420096769848103, 0.0034927386902514077, 0.003612578288921742, 0.003708819429431651, 0.0037548552439947533, 0.003926804405728195, 0.003986899519184352, 0.004150110666939558, 0.0042505955540944156, 0.004313258774779225, 0.004451833504040383, 0.004585576709532848, 0.00474246544582659, 0.004874725219025508, 0.0050109930045050075, 0.00519861197305503, 0.0053532917815432325, 0.005507217159343734, 0.005639948276795846, 0.005846251393295659, 0.0059885910695491075, 0.0062024243987689395, 0.006482190536846239, 0.006717685964573413, 0.0068833346378333864, 0.007162698069246487, 0.007360720048432485, 0.007666896260397409, 0.007872205705750526, 0.008195618705213311, 0.008513721381255083, 0.00890892176287896, 0.009258407579793916, 0.009724390710937272, 0.010108654793536417, 0.010548590217574743, 0.010975641578229983, 0.011526789103422348, 0.012043474015112746, 0.012716102564277389, 0.013431924604257003, 0.014299983010813377, 0.014857180822526022, 0.015642330796173286, 0.016838029708406442, 0.01813816538297235, 0.019193303872017592, 0.020692605291314385, 0.02283575070030013, 0.023809627336201407, 0.02590786359380735, 0.029216759524742952, 0.03276371294975062, 0.0373050352374768, 0.04091313698090973, 0.05046366123537228, 0.05737535745623545, 0.06781296688546185, 0.09593234674818188, 0.13058920213716907, 0.19115974684044465, 0.2909418067392709, 0.3800187769481447, 0.44170231320778347, 0.4714437717840313, 0.4818036960492877, 0.4783680745795133, 0.46935858320497825, 0.4587305204065562, 0.4533694823857326, 0.459745771617386, 0.47123282283088763, 0.4911090388737969, 0.5142589857573987, 0.5330176978784525, 0.5488274174881866, 0.5657266458670483, 0.5928593353164412, 0.6097422300183754, 0.6412667710670917, 0.6718627164205705, 0.7168596275875857, 0.7509383941806824, 0.7991365272600233, 0.8460771583423621, 0.8806532127041746, 0.9148140391083454, 0.9578909636679608, 1.0049816696656169, 1.0237682643313581, 1.0437563924173747, 1.0472309271558566, 1.0409612470144032, 1.031401756063377, 1.014757793264758, 0.9981960750289359, 0.9808809221716888, 0.961956726850273, 0.9474569159510084, 0.9393094487030882, 0.9336510977021794, 0.9313388924036582, 0.935095276990121, 0.9439935126513754, 0.9584765219506691, 0.9863567808054249, 1.0012119480166035, 1.0253164253757725, 1.0474087729073185, 1.0639462930995314, 1.0716405871933943, 1.0704028408813193, 1.0622033695417121, 1.043833529708262, 1.0338505187007858, 1.015920739275425, 1.0064992296788633, 0.9819865969757522, 0.9710214047330807, 0.9574262746567648, 0.9438420012803899, 0.9352296630969221, 0.9262723605657974, 0.9143976556576733, 0.906185974767554, 0.897383215865522, 0.8914749965773882, 0.8879405463472012, 0.8871835474522016, 0.8919458519851953, 0.8983896357670406, 0.9093959230438889, 0.9218004649480496, 0.9450752463505093, 0.9671534638935854, 0.9845923194765351, 1.0099783317164455, 1.0348019720006179, 1.0542051649436106, 1.0695348030281715, 1.0714895720146236, 1.0584659785594672, 1.0323155961704455, 1.0068494651997573, 0.9766212222615239, 0.9537070212961404, 0.9386442421675207, 0.9283715888948045, 0.9264333721172989, 0.9309027315540709, 0.9414312645241489, 0.9512243937933982, 0.9541815008827447, 0.9502313171063516, 0.9459266804869394, 0.9490871068258206, 0.9611317372570843, 0.973561768524303, 0.9767281726872864, 0.9746504804104761, 0.9760739848812415, 0.9857284436213061, 0.9994768460668217, 1.0120096923082202, 1.0202332946418209, 1.0284681451017978, 1.0382811587950693, 1.0531669717037115, 1.0685894535728213, 1.0805531919033358, 1.0813263051302948, 1.0705892813806233, 1.048181664775831, 1.0179617789283264, 0.9930058019952148, 0.9737605771788065, 0.956901950994992, 0.947912934022556, 0.9463876271578429, 0.9372201265736783, 0.9141571961975239, 0.893007930017307, 0.8778318381846293, 0.8747227563843458, 0.8840757774351256, 0.9003018875828541, 0.9213995195274476, 0.9448166962302658, 0.9670442734322858, 0.9816851902191693, 0.9851126191393962, 0.978436000613713, 0.9703058259414604, 0.9680973945158586, 0.9861023719360671, 1.0293833677749844, 1.08384215914281, 1.111718795380229, 1.1133143253546647, 1.100968217705162, 1.088742451915961, 1.0745659328965276, 1.0482283379704669, 1.012858154711397, 0.9695660135619018, 0.9225397391556798, 0.8763529455899852, 0.8367622011927748, 0.8123223750512327, 0.8097816988491554, 0.8313224468937503, 0.8691678356053407, 0.9136821706528581, 0.9571126554725979, 0.9859218357592383, 1.0016025832763975, 0.9992052304862202, 0.9828365727860549, 0.9675478727780485, 0.960288048032337, 0.9575723763840182, 0.9459033579045835, 0.928629245702154, 0.92538286418184, 0.9460944962180263, 0.9889581495418057, 1.031732343423515, 1.0467118979879293, 1.048616786887183, 1.0423354385858685, 1.0242110167521536, 0.991840038202967, 0.951506095277355, 0.9137675908625736, 0.8851941604226505, 0.8681313490708437, 0.8607969921071046, 0.8641111862588893, 0.8782470324218741, 0.8887788698562009, 0.8890278099168015, 0.8918036221219175, 0.895886764053666, 0.8978596614881845, 0.9051262511128065, 0.919803124946569, 0.9333081326461052, 0.9425770473133332, 0.9475210733762942, 0.9486246471060673, 0.948503501988093, 0.9533274428918959, 0.9605420156753953, 0.9697746961106184, 0.977695236899694, 0.9690359320553655, 0.9487105499087658, 0.9300743693063865, 0.9097073638370113, 0.8910846796724877, 0.8830012244861232, 0.8863239700424359, 0.8956749190218679, 0.9031167883145368, 0.8978899631893845, 0.8811875723598379, 0.859655376928613, 0.8500345049121885, 0.8611878573593885, 0.8729885537549915, 0.874116348696109, 0.8796797323397759, 0.8877526681108092, 0.8917752593669692, 0.8960371376471874, 0.900657865146255, 0.904774891780094, 0.9087473420100676, 0.9117511295599522, 0.9135793621216065, 0.9150858626021513, 0.9161985573456314, 0.9151096912647669, 0.9087303206121661, 0.8952344997385756, 0.8776874638052895, 0.8634698046250509, 0.8588300424175542, 0.8599380553978295, 0.864971057625807, 0.8724402825343934, 0.8770147862647181, 0.8762307232857258, 0.8707738762433272, 0.8589156625958114, 0.8472456494588244, 0.8459518448203727, 0.8510190673585086, 0.850631963256785, 0.845632577956862, 0.8436984280124931, 0.844444344178521, 0.8489113384220258, 0.8562464364577218, 0.8628222378432897, 0.8705506216447824, 0.8798582207887572, 0.8878468098004099, 0.8910705436721164, 0.8884813524000628, 0.8785756482321746, 0.8604141792675287, 0.8417813852236796, 0.8345225319318809, 0.8332363330163536, 0.8335769069384694, 0.8336001264196775, 0.8350243281523274, 0.8378009932171215, 0.8391605522156163, 0.8371169039729289, 0.8357881328656634, 0.8353605414721634, 0.8328024514597049, 0.8262147752828384, 0.8235192611641813, 0.8242649449623668, 0.8246019579925282, 0.8275618787287851, 0.8303430653326241, 0.8308592571354301, 0.8315739831843411, 0.8346893774887623, 0.8389627811596606, 0.8390459582301991, 0.8362653214567936, 0.8330351891879905, 0.8298867314842828, 0.8281727720077072, 0.8252717399659462, 0.8199374883044795, 0.8162225000404564, 0.814540697314283, 0.8138765176403647, 0.81484213878108, 0.8133752791147207, 0.8089386028176101, 0.808048804754039, 0.8086785488596872, 0.8077725141789572, 0.8048877372957235, 0.8029139417239112, 0.8018443478566591, 0.7991344689713755, 0.799484996592239, 0.801264063351763, 0.8030823327317724, 0.8068714498223677, 0.8079834410876029, 0.8077965137077715, 0.8074276705633521, 0.8046302829934532, 0.7997821985478527, 0.7966294971521881, 0.7969346059744276, 0.7982565173305309, 0.7989694168138948, 0.7992575374569093, 0.7976427690924305, 0.7935308899361562, 0.7898388481627849, 0.7866362654351766, 0.7846411253841429, 0.7841420008691155, 0.7833830468278303, 0.7817280077968847, 0.7808496137309235, 0.7797833448119017, 0.7789978108558054, 0.7793128083848722, 0.7787419169414673, 0.7784758790148905, 0.7798841116383484, 0.7811655203486657, 0.7812714711305224, 0.780701149490558, 0.7797701779535879, 0.7781770193602449, 0.777030021169209, 0.7761001142965999, 0.7744415675874925, 0.7718748829864913, 0.7703308442124279, 0.7690971389252758, 0.7690970145823814, 0.7691590965288193, 0.7691909222616156, 0.7669678772611456, 0.7654865207074447, 0.7655065376839782, 0.7641491084100974, 0.7621769367727041, 0.7604408956375445, 0.7592457884950148, 0.7585868273485861, 0.757462337576954, 0.7560636293733969, 0.7555122837528488, 0.7554103553035594, 0.7567243680090808, 0.7581896362839402, 0.7577448164145653, 0.7568231540246451, 0.7554643653312152, 0.754838256096104, 0.7540700256375796, 0.7529147283902344, 0.7513731220965899, 0.7502745135029728, 0.7488637049753226, 0.748237582990777, 0.7475744532891718, 0.7465998941454912, 0.7459665885724747, 0.746295773541647, 0.7454337211605997, 0.7440695868338536, 0.7426031571156732, 0.7412342284176822, 0.7403036222254931, 0.7394589775535491, 0.7387561220728324, 0.738433149792466, 0.738465572655102, 0.7387422865808491, 0.7386845016267082, 0.7381869819219466, 0.736976476444971, 0.7359972134208362, 0.7358737991980561, 0.7352166313837596, 0.7346394440782477, 0.7343616145769476, 0.7335646178909557, 0.7327717419497592, 0.7324235047453705, 0.7316864596625139, 0.7307338355589155, 0.7295753371667244, 0.7284052928998715, 0.727457137510049, 0.7267244324018733, 0.7262713882244468, 0.7259811359289443, 0.72532938830196, 0.7246036443650106, 0.7237135834990172, 0.7229435482036725, 0.7222591519584626, 0.722016334793744, 0.7220991045802327, 0.7218611818396256, 0.7209262783531648, 0.7198225924112435, 0.7194855339883839, 0.7188140767385023, 0.7181801561562398, 0.718092828107957, 0.7175356218977131, 0.7170902806972413, 0.7166939060493905, 0.7161168933507416, 0.7156247523976297, 0.7149429454753782, 0.7141445887870532, 0.7134122838903431, 0.7127424389910128, 0.7121422376775324, 0.7115261473273506, 0.7107219096804316, 0.7099376945644442, 0.7093044029856088, 0.7087622961959579, 0.7082869084496529, 0.7079536259899677, 0.7077334047423458, 0.7071272001058934, 0.7065611533336065, 0.705895645221598, 0.7055662706515405, 0.7053163380776313, 0.7048551823944541, 0.7043879328975156, 0.7037585428335907, 0.7032195798163282, 0.7027000944426434, 0.7023206571705388, 0.7020119405095577, 0.7016580778821513, 0.7011169289407545, 0.7008327035244808, 0.7003778449481578, 0.6998830445501956, 0.6993389583168607, 0.6988135824225057, 0.69820234824813, 0.697519723023541, 0.6971711741558635, 0.6966611328487795, 0.6963775459980092, 0.6962571799268565, 0.6960744602859082, 0.6956837932345346, 0.695147738822982, 0.6949841679807836, 0.6945444056376935, 0.694191494483825, 0.6938677402068127, 0.6933946811179179, 0.6931708202656802, 0.6926576001252616, 0.6923796590664063, 0.6920406543181964, 0.691786101979905, 0.6914244662536525, 0.6911236590500764, 0.6908863952438701, 0.6904875453201085, 0.6901789000257652, 0.6897691555422796, 0.6893387404002828, 0.6890754634896084, 0.6886587514601786, 0.6884753550054226, 0.6881397100721274, 0.6878283307789372, 0.6876330596046916, 0.6873600907113029, 0.6871657965094644, 0.6870727015496341, 0.6866108212422687, 0.6865639715856483, 0.6862945881117769, 0.6860516466400686, 0.6863365746139243, 0.6857706578254781, 0.6855092029331639, 0.6853398224990094, 0.6851772808717037, 0.6850050631984768, 0.6847992199346179, 0.6846373328096086, 0.6844744591501166, 0.6843600199053266, 0.6841592807333867, 0.6839341599316826, 0.6840228801007908, 0.6836696827825032, 0.6835503299894499, 0.6834622167571626, 0.6833881291075529, 0.6832714975639109, 0.6831317902077394, 0.6830724986327636, 0.6831358538249827, 0.6829960336665097, 0.6827793927084341, 0.6828120468525234, 0.6827243465811366, 0.6828479758466008, 0.6826868209652439, 0.6828446111466192, 0.6826304594586352 ] } ], "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": 5, "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)$', delay_draw=True)\n", " " ] }, { "cell_type": "code", "execution_count": 6, "id": "b8f6a375", "metadata": { "scrolled": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Result: E0 = 8977.580 eV, R= 2.5494 +/- 0.0017 Ang\n", "=================== FEFFIT RESULTS ====================\n", "[[Statistics]]\n", " nvarys, npts = 4, 104\n", " n_independent = 13.732\n", " chi_square = 80.5162009\n", " reduced chi_square = 8.27300959\n", " r-factor = 0.00122781\n", " Akaike info crit = 32.2884981\n", " Bayesian info crit = 34.7675287\n", " \n", "[[Data]]\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.00135189, std=9.20931e-4)\n", " epsilon_r = 0.03323044\n", " n_independent = 13.732\n", " \n", "[[Variables]]\n", " amp = 0.92213823 +/- 0.02713764 (init= 1.00000000)\n", " del_e0 = 5.67728180 +/- 0.37592135 (init= 0.00000000)\n", " del_r = 0.00155256 +/- 0.00168163 (init= 0.00000000)\n", " sig2 = 0.00354962 +/- 1.85415e-4 (init= 0.00000000)\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", "[[Paths]]\n", " = Path 'p3tf6mewgr' = 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.54780000\n", " degen = 12.0000000\n", " n*s02 = 0.92213823 +/- 0.02713764 := 'amp'\n", " e0 = 5.67728180 +/- 0.37592135 := 'del_e0'\n", " r = 2.54935256 +/- 0.00168163 := 'reff + del_r'\n", " deltar = 0.00155256 +/- 0.00168163 := 'del_r'\n", " sigma2 = 0.00354962 +/- 1.85415e-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": "a3beac0b-5567-47ee-adfb-17b7998feacd", "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.0010506890834946086, 0.0010272672874688726, 0.0009799709008622966, 0.0009515699433551839, 0.0009483519570436642, 0.0009313222884447521, 0.0008798530482730621, 0.0008479782340696277, 0.0009234953240868927, 0.0010837508725907755, 0.0012050885521563664, 0.001193776173721225, 0.001061441943150244, 0.0009792706537519527, 0.0011541233941536745, 0.0014697593356080394, 0.0016713600290083014, 0.0016275550078825752, 0.0014323155757771234, 0.0015121268903346527, 0.0020801436733455933, 0.002679497355019709, 0.0029192338483079853, 0.0033254953573596106, 0.006063782515140518, 0.012142722951524392, 0.021541354279667933, 0.034218295754562374, 0.04988004942215747, 0.0678183427322074, 0.08685008627640337, 0.10534361763676299, 0.12133927753308348, 0.13276815082532034, 0.13776135367352554, 0.1350376294165554, 0.12438280250323294, 0.10735165561359465, 0.08863911006480207, 0.07819766372539987, 0.08688193510640155, 0.11202869061754515, 0.14336067833255872, 0.17422347822193318, 0.20182441914246574, 0.2258728320536203, 0.24760974268548355, 0.2687660285497804, 0.2905090594072613, 0.3128842494858568, 0.33513271181184345, 0.3567205660795549, 0.3785273586814896, 0.4035942670856644, 0.4369269253760075, 0.4842118962217851, 0.5502165693412508, 0.6383018039524001, 0.7516155474964833, 0.8950000270612962, 1.0761694377469098, 1.3052678516955705, 1.5928063324767692, 1.946803038845992, 2.3702122275056503, 2.8592974936518014, 3.403030854678771, 3.9833505130334435, 4.57613079254178, 5.152776120748234, 5.682343642917385, 6.134039547282718, 6.479864852840364, 6.697139056085482, 6.770620961563661, 6.693979558465139, 6.470439454875185, 6.112524172654608, 5.640931417109535, 5.082680574040849, 4.468758529397476, 3.831543130779206, 3.202296373464338, 2.608987516879579, 2.0746248852037406, 1.6161289810626356, 1.2435400983230926, 0.9590634322100768, 0.7556093773300679, 0.6162064698542051, 0.5176857245555437, 0.43852287592886474, 0.36468937860194073, 0.29006613160672484, 0.21396975499150822, 0.13941419765604793, 0.07875223914525942, 0.08282655054093233, 0.151519583015727, 0.2409220706678012, 0.34392368056779626, 0.4603410335634093, 0.5893278352931033, 0.7279664669933241, 0.8712030333985333, 1.0123313413985833, 1.1437426422157801, 1.257784588779687, 1.3476185679348798, 1.4079821568666102, 1.4357762983238949, 1.4304166310376014, 1.3939190081590238, 1.3307293022509001, 1.2473525397083964, 1.1518781480419766, 1.053521129074918, 0.9622697466989368, 0.8885874117309159, 0.842828457665412, 0.8338643361516537, 0.867106391314499, 0.9434393337345448, 1.060167581711771, 1.2128126312237624, 1.3960823806398406, 1.6037760744527916, 1.8282807148792772, 2.0602262906923023, 2.2885428881815404, 2.5009438349857294, 2.684747746075216, 2.827899737761232, 2.920032886086751, 2.953417197858503, 2.9236722637810972, 2.8301686882514683, 2.6761089631390216, 2.468359190122319, 2.2172059256962964, 1.9363664464161003, 1.643856830745788, 1.364777978422133, 1.1369143530521932, 1.0133536493929227, 1.0371022755133021, 1.19496136476307, 1.433942647527803, 1.7082451004927557, 1.9894908943358842, 2.2600255106017495, 2.5073713097195376, 2.7215405244649706, 2.894032621440644, 3.0176943501794655, 3.0869960270228214, 3.0984590696828938, 3.051059768954774, 2.9464952540637013, 2.7892537032133227, 2.5864878577351287, 2.3477476105878883, 2.0846816258494596, 1.810867764034677, 1.5419630902024564, 1.2962759402471355, 1.0952156312638708, 0.9609120461002989, 0.9063305371249158, 0.9225895071877009, 0.9821740517997933, 1.0559371279976382, 1.123158680342721, 1.1719094033007371, 1.1965933444748271, 1.195852446075037, 1.1711842576706812, 1.1260115082831987, 1.0649807030454288, 0.9933702768654777, 0.9165576887041388, 0.8395338260163796, 0.7664770629433181, 0.700418158557834, 0.6430420314861195, 0.5946728249060625, 0.5544592358378116, 0.520720388855827, 0.4913632140112482, 0.46427748738942165, 0.4376498000407916, 0.4101750843140376, 0.3811623262078369, 0.3505429146182163, 0.3188201018927686, 0.28705396609502565, 0.2570305920093848, 0.23171661393052181, 0.21568530197930938, 0.21420603177305197, 0.22986539525367217, 0.25971776337644764, 0.2966872846449353, 0.33254480774814793, 0.35946808682525594, 0.37059845833756916, 0.3604340091634132, 0.32530158776227536, 0.26398905532255357, 0.17931172773134083, 0.09056843638081911, 0.11757224748326404, 0.24850205246840967, 0.3970409968839822, 0.5457093421695743, 0.6850479059632874, 0.8077092791003218, 0.9081382437714499, 0.9827523212786201, 1.0300634396010238, 1.0506401446814368, 1.046908249059009, 1.0228271977759826, 0.9834999949038933, 0.9347804067810226, 0.8829267626863182, 0.8343066807389197, 0.7950780450233831, 0.7706976843917894, 0.7651607631324384, 0.7801559591013164, 0.8146232622830692, 0.86505068187682, 0.9262985630421231, 0.992465452309345, 1.057518994395201, 1.1157020524000623, 1.161824430134835, 1.1915193208978911, 1.2014893116810932, 1.1897328001333005, 1.155729818584624, 1.1005687760665206, 1.027005998886428, 0.9394624816807337, 0.8439670071856583, 0.7480260588984678, 0.6602795102093263, 0.5895244506305779, 0.5424978406624167, 0.520751411698576, 0.5192267781869847, 0.5286989019813041, 0.5398938940017757, 0.5460027513287832, 0.5431812194543019, 0.5301065981476492, 0.507328457992421, 0.47664637804957527, 0.4405334488150909, 0.40158472520802757, 0.3620034395326099, 0.323207371512822, 0.28570944833039397, 0.24944478558460725, 0.21466532480785858, 0.18348641978517383, 0.16192047646022228, 0.15980733680199646, 0.1830153198922775, 0.22667712044046925, 0.2817293284683284, 0.34089406423660057, 0.3990963228699514, 0.4527579619409608, 0.49939157107390136, 0.5374248031241383, 0.5661002623056566, 0.5853796774841762, 0.5958340738923884, 0.5985251024456889, 0.5948874957158293, 0.58661526276341, 0.5755408706745238, 0.5634861115469465, 0.5520670111672303, 0.5424616374275041, 0.5351932404050211, 0.5300129960612578, 0.5259506708280738, 0.5215348005011717, 0.5151119395079454, 0.5051672198604354, 0.49057166931175356, 0.47072496694667965, 0.4455974892505059, 0.4156945226080384, 0.38197278067381457, 0.3457408023397329, 0.30857329121844523, 0.2722637608789989, 0.23882195321411817, 0.21047114966234814, 0.18948969180190164, 0.17764817872521751, 0.1752861621207527, 0.18082752116116943, 0.19141341654308733, 0.2040102115219455, 0.21607868988192147, 0.22573187287598875, 0.23169176577566855, 0.23324241843601579, 0.23023888221585545, 0.2231779462888388, 0.21332032354294286, 0.20282093274536142, 0.19471532127027832, 0.19242464142265683, 0.19852128766741883, 0.21338137697313753, 0.2350228264691539, 0.26018222416232983, 0.285449398377466, 0.3078518523289017, 0.32504708399796206 ] }, { "line": { "color": "#d62728", "width": 3 }, "name": "fit", "type": "scatter", "uid": "262359b8-3480-4cc4-81c9-91afe6cdfb96", "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.0008755589985193371, 0.0011366968162749925, 0.0016604806974026593, 0.002174804553565905, 0.002535125649030716, 0.0026685222891520803, 0.0025997712625759994, 0.0024571083708010985, 0.0023661999464671883, 0.0022676081690722323, 0.001941155092799849, 0.0013333608305178106, 0.001484054208324042, 0.0031254409167742227, 0.005205102991810471, 0.007156416934963238, 0.008637712376882772, 0.009567740191940389, 0.010245952275822727, 0.011267469501806135, 0.013027918637765996, 0.015303262802333153, 0.01753633920140025, 0.01941591295839824, 0.021278617349559412, 0.0241094451009511, 0.02884365093247615, 0.03550982023080331, 0.04329463002760492, 0.051195827777319285, 0.058494783728984034, 0.06491449663591412, 0.07051009145706205, 0.07536171148338332, 0.07926770305821114, 0.08170732038778579, 0.08219276958652615, 0.080888559299152, 0.07920881494527073, 0.07984232156957424, 0.08541460697555008, 0.09641535408346436, 0.11107052606226414, 0.1272610011190926, 0.14409963470191728, 0.1622975839948058, 0.1834883298434719, 0.20894788485728388, 0.23863312342808818, 0.27132026544966, 0.3056815818640541, 0.3415494600063189, 0.3806968435654068, 0.4266960008887748, 0.48379084472069084, 0.555424228746983, 0.6435350815379611, 0.7492058146148899, 0.8742231004550474, 1.0226633407085004, 1.2017325292118028, 1.4213550721088257, 1.692414269645475, 2.0241294140069543, 2.42143666779152, 2.8830780349588214, 3.4006330694164744, 3.958420704228557, 4.534162770008444, 5.10035292315697, 5.626270932629834, 6.080509474705314, 6.433784314345, 6.661723081406701, 6.747298447978095, 6.68259711653934, 6.4696926750017045, 6.120504950795534, 5.6556619894464735, 5.102511703451101, 4.492538129581615, 3.8585058449697813, 3.2316756632305585, 2.639402926850397, 2.1033509271961224, 1.6384352761042227, 1.2524730850916739, 0.9463695985285621, 0.7146165541731428, 0.5460923480766691, 0.42570968033361, 0.33751682267752875, 0.26835548160482753, 0.20991755104853577, 0.15843267726449425, 0.11312630898972183, 0.0746957454125684, 0.04431613056787182, 0.02335236890231468, 0.013689256557244239, 0.012968687862943942, 0.013266609758483533, 0.013207550549434906, 0.014949904888812255, 0.018661277445311493, 0.022371568058251454, 0.024451936288438135, 0.02421640384289092, 0.02178120188731607, 0.01787243378353725, 0.013626245486695355, 0.010299099272172568, 0.008536998506990476, 0.007560464746693494, 0.006130520151750005, 0.004080195250025942, 0.0038730288852541336, 0.0071650604346966814, 0.011151475286031649, 0.01447874506278477, 0.016512972996921386, 0.016981127943224185, 0.015934008552649342, 0.01371267790662519, 0.0108910956732105, 0.008202716593058736, 0.006382735576463861, 0.005668551188795669, 0.005444198101986995, 0.005036341585700989, 0.004356321939916254, 0.003875118968672429, 0.00415020710998767, 0.004964375157587173, 0.005668569850786967, 0.005856093869619596, 0.005429975352699371, 0.004562486429219832, 0.0036938064142341912, 0.003414870226727423, 0.003815538062172682, 0.004336921568654774, 0.004551042089714263, 0.0043858754321293705, 0.004069005300775588, 0.0040073805057469665, 0.0044127624243420595, 0.005035474027081855, 0.0054781226008309, 0.005481581092896932, 0.0049700659256106915, 0.004030707652107421, 0.002896524325327064, 0.001976871201475572, 0.0017943074676653658, 0.0021758382918279453, 0.0025440329794099307, 0.0027574872352479434, 0.002965920920958545, 0.0033567161336562636, 0.003927073843971562, 0.004493493081710762, 0.004853974856691634, 0.004883872157629551, 0.004560839787301827, 0.003960449960176725, 0.00323689722365187, 0.0025858512943207697, 0.0021603862312118043, 0.0019421756213973394, 0.001757632701716571, 0.0014694670436956557, 0.0010879618187205782, 0.0008257526693515502, 0.0009764469338983548, 0.0013351010272118328, 0.0016055218571684533, 0.001670961560433169, 0.0015108375882609913, 0.0011647841857115742, 0.0007149002698549017, 0.00029292613989382346, 0.00029593956741873645, 0.0005302683840599748, 0.0006662637991697918, 0.0007567822500332043, 0.0009434115845812758, 0.001285828819844227, 0.0017079147394752393, 0.0021032501907206186, 0.002389841440358398, 0.0025241584106155636, 0.002508469273992597, 0.002393767086064457, 0.002269936374802577, 0.002225896014569558, 0.0022845608845848517, 0.0023832566325043213, 0.0024310483516314933, 0.00237049961725181, 0.002205438731742753, 0.0020069296164223816, 0.0018922713799048372, 0.001942217721607351, 0.002112766247160499, 0.0022837389535043133, 0.00236009843723454, 0.0023114801028093616, 0.0021729960771467703, 0.0020283880813652625, 0.0019633244636943414, 0.001994889501155764, 0.0020556386099926114, 0.002059006884893707, 0.001960107516114221, 0.0017759489898534397, 0.0015810476461415603, 0.0014725558810680128, 0.0014875487591668095, 0.0015579645541574161, 0.0015844523919984058, 0.0015109925237451346, 0.001341570483613784, 0.0011334783236413906, 0.0009762518871033192, 0.0009271811536967529, 0.0009401455361783917, 0.0009199867508473765, 0.0008102915820632639, 0.0006121260445302826, 0.00038643918705194927, 0.00030202502462173454, 0.00043497041705982354, 0.0005768926369963447, 0.0006451367207190774, 0.0006586738999811494, 0.0006989249199313814, 0.0008399229089437836, 0.001059907710432825, 0.0012845556892533661, 0.001456260897503171, 0.00155424431179483, 0.0015965588821949773, 0.001630802662096842, 0.0017057841346575333, 0.001834081233702482, 0.0019829758742305158, 0.0021023802652798034, 0.002155953115428788, 0.0021363323360700085, 0.00206659840758653, 0.001990259927413143, 0.0019488897195504417, 0.001955360144035735, 0.001985663434764942, 0.0019986044460113986, 0.0019622118481709655, 0.0018682880582329875, 0.0017343233863739028, 0.0015960077976446836, 0.001490127215740728, 0.0014312806698156936, 0.0014005503157629787, 0.0013604590592495812, 0.0012799117678423885, 0.0011482566311524698, 0.0009769169605724735, 0.0007940493301750085, 0.0006337464386625012, 0.0005172031979178893, 0.00043357066859501153, 0.000349708374812996, 0.00024734914609154555, 0.00016974297672037015, 0.00023639554644471034, 0.0003848263088279218, 0.0005266054899657991, 0.000635470062276795, 0.0007153343363161229, 0.0007918799773091769, 0.000894836159964041, 0.0010320055752590107, 0.0011820996371825557, 0.0013117827827444034, 0.0013951010238119948, 0.0014245281572187335, 0.0014139506544449045, 0.0013928747313748488, 0.0013900072238560793, 0.0014128659319428145, 0.0014419411389324172, 0.0014460323328827655, 0.0014023971851940497, 0.0013078524260732869, 0.0011797025467809442, 0.001048209455076796, 0.0009405805985143123, 0.0008618367559020382, 0.0007906465739326789, 0.0006977752398669309, 0.0005665756353910262, 0.00039987752875998727, 0.00021663164988037714, 5.3816045507530044e-05, 0.00011885753096615988, 0.00022203306979437276, 0.00029807692720234874, 0.0003763760163435654, 0.00048496803878029496, 0.0006224951985971883, 0.0007644276112129765, 0.0008837818278986253, 0.000964009536092805, 0.0010047003453704415, 0.0010218243717340098, 0.0010404955777550701, 0.0010796497881258537, 0.0011382511432272238, 0.0011970484351253932, 0.0012331210188736915, 0.001233171078337426, 0.0011993639171178781, 0.0011482263187087198, 0.0011029251213730143, 0.0010798619303641012, 0.00107729875486592 ] } ], "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": "R (Ang)" }, "zerolinecolor": "#DDD" }, "yaxis": { "color": "#004", "gridcolor": "#D8D8D8", "showgrid": true, "title": { "text": "$|\\chi(R)|$" }, "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": 7, "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", "new_dat = change_calib(cudat, dspace, theta_off=0.02)" ] }, { "cell_type": "code", "execution_count": 8, "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": "2f178dbd-ad9d-4bcf-a513-5a44d1408dcc", "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.006259474761857858, 0.00020884413456537628, 0.00013367051997887315, 1.3030948957645212e-05, -6.785601444994545e-05, -0.0001377976297546854, -0.00015118830971820596, -0.00012801132607241026, -0.00010071531578492878, -0.00010496081280303482, -6.022165904272004e-05, 0.00014211281173055691, 0.0002530926523955572, 0.0005879856706472706, 0.0009093433105783736, 0.0017442259738558153, 0.001835559497609198, 0.001960564689123814, 0.001947414417208297, 0.001941489667325108, 0.0020752152178703845, 0.002161996036250319, 0.0022097513362068725, 0.002281163532364061, 0.002340322969503525, 0.0024043939727043514, 0.0024493451123328335, 0.0024952177794254496, 0.0025181841561758346, 0.0025823206096061294, 0.002631011585879763, 0.002778240780731134, 0.0028386950736110812, 0.0028897079044650628, 0.0029830252644179817, 0.0030577455309030387, 0.003156050106694802, 0.00322760384606963, 0.003323406208073335, 0.003420096769848103, 0.0034927386902514077, 0.003612578288921742, 0.003708819429431651, 0.0037548552439947533, 0.003926804405728195, 0.003986899519184352, 0.004150110666939558, 0.0042505955540944156, 0.004313258774779225, 0.004451833504040383, 0.004585576709532848, 0.00474246544582659, 0.004874725219025508, 0.0050109930045050075, 0.00519861197305503, 0.0053532917815432325, 0.005507217159343734, 0.005639948276795846, 0.005846251393295659, 0.0059885910695491075, 0.0062024243987689395, 0.006482190536846239, 0.006717685964573413, 0.0068833346378333864, 0.007162698069246487, 0.007360720048432485, 0.007666896260397409, 0.007872205705750526, 0.008195618705213311, 0.008513721381255083, 0.00890892176287896, 0.009258407579793916, 0.009724390710937272, 0.010108654793536417, 0.010548590217574743, 0.010975641578229983, 0.011526789103422348, 0.012043474015112746, 0.012716102564277389, 0.013431924604257003, 0.014299983010813377, 0.014857180822526022, 0.015642330796173286, 0.016838029708406442, 0.01813816538297235, 0.019193303872017592, 0.020692605291314385, 0.02283575070030013, 0.023809627336201407, 0.02590786359380735, 0.029216759524742952, 0.03276371294975062, 0.0373050352374768, 0.04091313698090973, 0.05046366123537228, 0.05737535745623545, 0.06781296688546185, 0.09593234674818188, 0.13058920213716907, 0.19115974684044465, 0.2909418067392709, 0.3800187769481447, 0.44170231320778347, 0.4714437717840313, 0.4818036960492877, 0.4783680745795133, 0.46935858320497825, 0.4587305204065562, 0.4533694823857326, 0.459745771617386, 0.47123282283088763, 0.4911090388737969, 0.5142589857573987, 0.5330176978784525, 0.5488274174881866, 0.5657266458670483, 0.5928593353164412, 0.6097422300183754, 0.6412667710670917, 0.6718627164205705, 0.7168596275875857, 0.7509383941806824, 0.7991365272600233, 0.8460771583423621, 0.8806532127041746, 0.9148140391083454, 0.9578909636679608, 1.0049816696656169, 1.0237682643313581, 1.0437563924173747, 1.0472309271558566, 1.0409612470144032, 1.031401756063377, 1.014757793264758, 0.9981960750289359, 0.9808809221716888, 0.961956726850273, 0.9474569159510084, 0.9393094487030882, 0.9336510977021794, 0.9313388924036582, 0.935095276990121, 0.9439935126513754, 0.9584765219506691, 0.9863567808054249, 1.0012119480166035, 1.0253164253757725, 1.0474087729073185, 1.0639462930995314, 1.0716405871933943, 1.0704028408813193, 1.0622033695417121, 1.043833529708262, 1.0338505187007858, 1.015920739275425, 1.0064992296788633, 0.9819865969757522, 0.9710214047330807, 0.9574262746567648, 0.9438420012803899, 0.9352296630969221, 0.9262723605657974, 0.9143976556576733, 0.906185974767554, 0.897383215865522, 0.8914749965773882, 0.8879405463472012, 0.8871835474522016, 0.8919458519851953, 0.8983896357670406, 0.9093959230438889, 0.9218004649480496, 0.9450752463505093, 0.9671534638935854, 0.9845923194765351, 1.0099783317164455, 1.0348019720006179, 1.0542051649436106, 1.0695348030281715, 1.0714895720146236, 1.0584659785594672, 1.0323155961704455, 1.0068494651997573, 0.9766212222615239, 0.9537070212961404, 0.9386442421675207, 0.9283715888948045, 0.9264333721172989, 0.9309027315540709, 0.9414312645241489, 0.9512243937933982, 0.9541815008827447, 0.9502313171063516, 0.9459266804869394, 0.9490871068258206, 0.9611317372570843, 0.973561768524303, 0.9767281726872864, 0.9746504804104761, 0.9760739848812415, 0.9857284436213061, 0.9994768460668217, 1.0120096923082202, 1.0202332946418209, 1.0284681451017978, 1.0382811587950693, 1.0531669717037115, 1.0685894535728213, 1.0805531919033358, 1.0813263051302948, 1.0705892813806233, 1.048181664775831, 1.0179617789283264, 0.9930058019952148, 0.9737605771788065, 0.956901950994992, 0.947912934022556, 0.9463876271578429, 0.9372201265736783, 0.9141571961975239, 0.893007930017307, 0.8778318381846293, 0.8747227563843458, 0.8840757774351256, 0.9003018875828541, 0.9213995195274476, 0.9448166962302658, 0.9670442734322858, 0.9816851902191693, 0.9851126191393962, 0.978436000613713, 0.9703058259414604, 0.9680973945158586, 0.9861023719360671, 1.0293833677749844, 1.08384215914281, 1.111718795380229, 1.1133143253546647, 1.100968217705162, 1.088742451915961, 1.0745659328965276, 1.0482283379704669, 1.012858154711397, 0.9695660135619018, 0.9225397391556798, 0.8763529455899852, 0.8367622011927748, 0.8123223750512327, 0.8097816988491554, 0.8313224468937503, 0.8691678356053407, 0.9136821706528581, 0.9571126554725979, 0.9859218357592383, 1.0016025832763975, 0.9992052304862202, 0.9828365727860549, 0.9675478727780485, 0.960288048032337, 0.9575723763840182, 0.9459033579045835, 0.928629245702154, 0.92538286418184, 0.9460944962180263, 0.9889581495418057, 1.031732343423515, 1.0467118979879293, 1.048616786887183, 1.0423354385858685, 1.0242110167521536, 0.991840038202967, 0.951506095277355, 0.9137675908625736, 0.8851941604226505, 0.8681313490708437, 0.8607969921071046, 0.8641111862588893, 0.8782470324218741, 0.8887788698562009, 0.8890278099168015, 0.8918036221219175, 0.895886764053666, 0.8978596614881845, 0.9051262511128065, 0.919803124946569, 0.9333081326461052, 0.9425770473133332, 0.9475210733762942, 0.9486246471060673, 0.948503501988093, 0.9533274428918959, 0.9605420156753953, 0.9697746961106184, 0.977695236899694, 0.9690359320553655, 0.9487105499087658, 0.9300743693063865, 0.9097073638370113, 0.8910846796724877, 0.8830012244861232, 0.8863239700424359, 0.8956749190218679, 0.9031167883145368, 0.8978899631893845, 0.8811875723598379, 0.859655376928613, 0.8500345049121885, 0.8611878573593885, 0.8729885537549915, 0.874116348696109, 0.8796797323397759, 0.8877526681108092, 0.8917752593669692, 0.8960371376471874, 0.900657865146255, 0.904774891780094, 0.9087473420100676, 0.9117511295599522, 0.9135793621216065, 0.9150858626021513, 0.9161985573456314, 0.9151096912647669, 0.9087303206121661, 0.8952344997385756, 0.8776874638052895, 0.8634698046250509, 0.8588300424175542, 0.8599380553978295, 0.864971057625807, 0.8724402825343934, 0.8770147862647181, 0.8762307232857258, 0.8707738762433272, 0.8589156625958114, 0.8472456494588244, 0.8459518448203727, 0.8510190673585086, 0.850631963256785, 0.845632577956862, 0.8436984280124931, 0.844444344178521, 0.8489113384220258, 0.8562464364577218, 0.8628222378432897, 0.8705506216447824, 0.8798582207887572, 0.8878468098004099, 0.8910705436721164, 0.8884813524000628, 0.8785756482321746, 0.8604141792675287, 0.8417813852236796, 0.8345225319318809, 0.8332363330163536, 0.8335769069384694, 0.8336001264196775, 0.8350243281523274, 0.8378009932171215, 0.8391605522156163, 0.8371169039729289, 0.8357881328656634, 0.8353605414721634, 0.8328024514597049, 0.8262147752828384, 0.8235192611641813, 0.8242649449623668, 0.8246019579925282, 0.8275618787287851, 0.8303430653326241, 0.8308592571354301, 0.8315739831843411, 0.8346893774887623, 0.8389627811596606, 0.8390459582301991, 0.8362653214567936, 0.8330351891879905, 0.8298867314842828, 0.8281727720077072, 0.8252717399659462, 0.8199374883044795, 0.8162225000404564, 0.814540697314283, 0.8138765176403647, 0.81484213878108, 0.8133752791147207, 0.8089386028176101, 0.808048804754039, 0.8086785488596872, 0.8077725141789572, 0.8048877372957235, 0.8029139417239112, 0.8018443478566591, 0.7991344689713755, 0.799484996592239, 0.801264063351763, 0.8030823327317724, 0.8068714498223677, 0.8079834410876029, 0.8077965137077715, 0.8074276705633521, 0.8046302829934532, 0.7997821985478527, 0.7966294971521881, 0.7969346059744276, 0.7982565173305309, 0.7989694168138948, 0.7992575374569093, 0.7976427690924305, 0.7935308899361562, 0.7898388481627849, 0.7866362654351766, 0.7846411253841429, 0.7841420008691155, 0.7833830468278303, 0.7817280077968847, 0.7808496137309235, 0.7797833448119017, 0.7789978108558054, 0.7793128083848722, 0.7787419169414673, 0.7784758790148905, 0.7798841116383484, 0.7811655203486657, 0.7812714711305224, 0.780701149490558, 0.7797701779535879, 0.7781770193602449, 0.777030021169209, 0.7761001142965999, 0.7744415675874925, 0.7718748829864913, 0.7703308442124279, 0.7690971389252758, 0.7690970145823814, 0.7691590965288193, 0.7691909222616156, 0.7669678772611456, 0.7654865207074447, 0.7655065376839782, 0.7641491084100974, 0.7621769367727041, 0.7604408956375445, 0.7592457884950148, 0.7585868273485861, 0.757462337576954, 0.7560636293733969, 0.7555122837528488, 0.7554103553035594, 0.7567243680090808, 0.7581896362839402, 0.7577448164145653, 0.7568231540246451, 0.7554643653312152, 0.754838256096104, 0.7540700256375796, 0.7529147283902344, 0.7513731220965899, 0.7502745135029728, 0.7488637049753226, 0.748237582990777, 0.7475744532891718, 0.7465998941454912, 0.7459665885724747, 0.746295773541647, 0.7454337211605997, 0.7440695868338536, 0.7426031571156732, 0.7412342284176822, 0.7403036222254931, 0.7394589775535491, 0.7387561220728324, 0.738433149792466, 0.738465572655102, 0.7387422865808491, 0.7386845016267082, 0.7381869819219466, 0.736976476444971, 0.7359972134208362, 0.7358737991980561, 0.7352166313837596, 0.7346394440782477, 0.7343616145769476, 0.7335646178909557, 0.7327717419497592, 0.7324235047453705, 0.7316864596625139, 0.7307338355589155, 0.7295753371667244, 0.7284052928998715, 0.727457137510049, 0.7267244324018733, 0.7262713882244468, 0.7259811359289443, 0.72532938830196, 0.7246036443650106, 0.7237135834990172, 0.7229435482036725, 0.7222591519584626, 0.722016334793744, 0.7220991045802327, 0.7218611818396256, 0.7209262783531648, 0.7198225924112435, 0.7194855339883839, 0.7188140767385023, 0.7181801561562398, 0.718092828107957, 0.7175356218977131, 0.7170902806972413, 0.7166939060493905, 0.7161168933507416, 0.7156247523976297, 0.7149429454753782, 0.7141445887870532, 0.7134122838903431, 0.7127424389910128, 0.7121422376775324, 0.7115261473273506, 0.7107219096804316, 0.7099376945644442, 0.7093044029856088, 0.7087622961959579, 0.7082869084496529, 0.7079536259899677, 0.7077334047423458, 0.7071272001058934, 0.7065611533336065, 0.705895645221598, 0.7055662706515405, 0.7053163380776313, 0.7048551823944541, 0.7043879328975156, 0.7037585428335907, 0.7032195798163282, 0.7027000944426434, 0.7023206571705388, 0.7020119405095577, 0.7016580778821513, 0.7011169289407545, 0.7008327035244808, 0.7003778449481578, 0.6998830445501956, 0.6993389583168607, 0.6988135824225057, 0.69820234824813, 0.697519723023541, 0.6971711741558635, 0.6966611328487795, 0.6963775459980092, 0.6962571799268565, 0.6960744602859082, 0.6956837932345346, 0.695147738822982, 0.6949841679807836, 0.6945444056376935, 0.694191494483825, 0.6938677402068127, 0.6933946811179179, 0.6931708202656802, 0.6926576001252616, 0.6923796590664063, 0.6920406543181964, 0.691786101979905, 0.6914244662536525, 0.6911236590500764, 0.6908863952438701, 0.6904875453201085, 0.6901789000257652, 0.6897691555422796, 0.6893387404002828, 0.6890754634896084, 0.6886587514601786, 0.6884753550054226, 0.6881397100721274, 0.6878283307789372, 0.6876330596046916, 0.6873600907113029, 0.6871657965094644, 0.6870727015496341, 0.6866108212422687, 0.6865639715856483, 0.6862945881117769, 0.6860516466400686, 0.6863365746139243, 0.6857706578254781, 0.6855092029331639, 0.6853398224990094, 0.6851772808717037, 0.6850050631984768, 0.6847992199346179, 0.6846373328096086, 0.6844744591501166, 0.6843600199053266, 0.6841592807333867, 0.6839341599316826, 0.6840228801007908, 0.6836696827825032, 0.6835503299894499, 0.6834622167571626, 0.6833881291075529, 0.6832714975639109, 0.6831317902077394, 0.6830724986327636, 0.6831358538249827, 0.6829960336665097, 0.6827793927084341, 0.6828120468525234, 0.6827243465811366, 0.6828479758466008, 0.6826868209652439, 0.6828446111466192, 0.6826304594586352 ] }, { "line": { "color": "#d62728", "width": 3 }, "name": "E0=8991.5", "type": "scatter", "uid": "303a8a70-63d0-452c-8de8-4b3be12c4dc3", "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.006259651514562021, 0.00020891455129295528, 0.0001337034046604888, 1.3033072635282131e-05, -6.787688702161794e-05, -0.00013783410989637407, -0.00015123249346252256, -0.00012805549221203854, -0.00010075207279176699, -0.00010498321947641967, -6.0221760485826464e-05, 0.00014214380966110694, 0.00025316119709842445, 0.0005881014905019659, 0.0009095135164788676, 0.0017444630396850791, 0.0018358006997278705, 0.0019608101286103705, 0.0019476626947666476, 0.0019417408733291375, 0.00207547134781591, 0.0021622558045513267, 0.0022100146270470538, 0.002281430617372045, 0.0023405937265089764, 0.0024046690447266828, 0.0024496234697047882, 0.0024954997296578437, 0.0025184694875466963, 0.002582610322481835, 0.0026313049792042, 0.002778538589814956, 0.002838996708206949, 0.0028900138786474992, 0.002983335436527682, 0.003058059733624645, 0.0031563682997398593, 0.0032279260616319924, 0.003323733320755675, 0.003420428212388289, 0.0034930742185442947, 0.0036129181112235396, 0.0037091639179728493, 0.0037552036019285247, 0.003927157946874797, 0.003987257106631571, 0.004150473373690431, 0.004250963076141894, 0.004313630420183202, 0.004452210067046839, 0.0045859581570418025, 0.0047428520388735785, 0.004875117044303831, 0.005011389790930111, 0.005199014270783654, 0.005353699261927236, 0.005507630175022482, 0.005640366292278217, 0.005846675171211517, 0.005989019986390021, 0.006202859188083989, 0.006482631898452367, 0.006718133806934797, 0.006883787909747761, 0.007163157968168527, 0.0073611857494677874, 0.00766736889623795, 0.007872684588688282, 0.008196104734091233, 0.008514214518205105, 0.008909422806302808, 0.009258916085917656, 0.009724908234767328, 0.010109180170977117, 0.010549124035617238, 0.01097618372140134, 0.011527340859839787, 0.012044035039811081, 0.012716674850532659, 0.013432508233798443, 0.01430057957701442, 0.014857787148188103, 0.01564294923359407, 0.01683866448089126, 0.01813881831965657, 0.01919397174824073, 0.02069329266617665, 0.022836464196346218, 0.023810354602872372, 0.025908616544383513, 0.029217551359379212, 0.03276454534897278, 0.03730591841572141, 0.04091406099853154, 0.0504646873913707, 0.05737645871899703, 0.06781417940946732, 0.09593385252038601, 0.13059106698106393, 0.19116223669279023, 0.29094532340428453, 0.3800232107508796, 0.4417073834732756, 0.4714491516381645, 0.48180918649529497, 0.47837353429423235, 0.4693639550820587, 0.4587357878870818, 0.45337469944639486, 0.459751058552876, 0.4712382320166518, 0.49111465671879767, 0.5142648454029781, 0.5330237543360901, 0.5488336405765148, 0.5657330467472111, 0.5928660188737503, 0.6097490912432127, 0.6412739604369333, 0.6718702239961332, 0.7168676009672889, 0.7509467214966776, 0.799145353205536, 0.8460864700488654, 0.8806628835026118, 0.9148240647406097, 0.9579014355312423, 1.0049926297877438, 1.0237794213777036, 1.0437677591529204, 1.0472423343981405, 1.0409725949240982, 1.031413010968203, 1.0147688825685446, 0.9982069996118936, 0.9808916743095815, 0.9619672905434906, 0.9474673360819917, 0.9393197903999938, 0.933661386476419, 0.9313491625694594, 0.9351055907258804, 0.9440039226820063, 0.9584870855076356, 0.986367635182726, 1.0012229597671307, 1.025327689301762, 1.047420269352543, 1.0639579637523362, 1.0716523419020298, 1.0704145881716547, 1.062215038074122, 1.0438450152980858, 1.033861907301948, 1.0159319494768062, 1.0065103486787543, 0.9819974701683356, 0.9710321709290384, 0.9574369069472214, 0.9438525002672254, 0.9352400792592491, 0.9262826903970542, 0.9144078692696996, 0.9061961097224083, 0.8973932661006019, 0.8914849917918348, 0.8879505108696201, 0.8871935097821791, 0.8919558757935615, 0.8983997383471616, 0.9094061512567722, 0.9218108337431061, 0.9450858677207674, 0.9671643256800989, 0.9846033742552727, 1.009989661559511, 1.0348135712553337, 1.054216978178185, 1.069546789174603, 1.0715015946847073, 1.0584778838804123, 1.0323272508417631, 1.0068608758518154, 0.9766323407578147, 0.9537179233051343, 0.9386550083105121, 0.9283822689840988, 0.9264440511035396, 0.9309134758590292, 0.9414421369733997, 0.9512353870282262, 0.9541925456563676, 0.9502423422016595, 0.9459376842917929, 0.9490981641542022, 0.9611429411924907, 0.9735731232528424, 0.9767395828281723, 0.9746618943616467, 0.9760854374486693, 0.985740019995052, 0.9994885884165909, 1.012021589065715, 1.0202453025404468, 1.0284802651934815, 1.0382934061391096, 1.0531793998596197, 1.06860206831131, 1.080565958757255, 1.0813391096698142, 1.0706020059339838, 1.0481941907629968, 1.0179740265651085, 0.9930178262328928, 0.9737724368453983, 0.9569136708240327, 0.947924596465451, 0.946399307483635, 0.9372317483162659, 0.9141686172981898, 0.8930191704084117, 0.8778429601966469, 0.8747338847744783, 0.8840870394532911, 0.9003133548345671, 0.9214112430755805, 0.9448287002286081, 0.9670565460410266, 0.9816976549387777, 0.9851251614337332, 0.9784485181998466, 0.9703183034706852, 0.9681098948329442, 0.9861151025379588, 1.0293965881126972, 1.083855984148738, 1.1117329540048346, 1.1133285487374582, 1.100982364372789, 1.0887565235521535, 1.0745799109282848, 1.0482420982342289, 1.0128716051641948, 0.9695790734781216, 0.9225523727388395, 0.8763651609394508, 0.8367740673770317, 0.8123340478886059, 0.8097934042276422, 0.8313344321144791, 0.8691802682615916, 0.9136951206639944, 0.9571261122998782, 0.9859356511818057, 1.0016166222953369, 0.9992193105988585, 0.982850550287564, 0.9675617615767151, 0.9603019298776464, 0.9575862995820209, 0.9459172313189198, 0.9286430136828209, 0.925396671092863, 0.946108588227368, 0.988972755443887, 1.0317474629503247, 1.0467272482395107, 1.0486322333775158, 1.0423508992795205, 1.024226373511555, 0.9918551432798062, 0.9515208690308764, 0.9137820618723507, 0.8852084233452567, 0.8681455251635339, 0.8608111792573239, 0.8641254969636336, 0.8782615757028409, 0.8887936122143049, 0.8890426467823086, 0.8918185815595879, 0.8959018602503506, 0.8978748737019009, 0.9051416357892262, 0.9198187588334537, 0.9333240060461079, 0.9425931176934293, 0.9475372972950664, 0.9486409861498136, 0.9485199459005962, 0.953344043309037, 0.9605587994987003, 0.9697916849804353, 0.9777124182697914, 0.9690531385066148, 0.9487276644752567, 0.9300914071921902, 0.9097243113912687, 0.8911015557904373, 0.8830181381499019, 0.8863410407634135, 0.895692211209813, 0.9031342819220146, 0.897907531039829, 0.881205097986656, 0.8596728118941978, 0.8500519740784134, 0.8612055746744468, 0.8730065287189318, 0.874134473058781, 0.8796980527753265, 0.8877712114131245, 0.8917939870267187, 0.8960560550775505, 0.9006769772381978, 0.9047941909829731, 0.9087668334955809, 0.9117708028607784, 0.9135992082214551, 0.9151058795771078, 0.9162187442497, 0.915130024943039, 0.9087507519538081, 0.8952549552739434, 0.8777079053274544, 0.8634902675363435, 0.8588506262733355, 0.8599588244474258, 0.8649920494759978, 0.8724615275659378, 0.8770362541641691, 0.8762523626742211, 0.8707956427741853, 0.8589374921170402, 0.8472675454303857, 0.8459739150114354, 0.8510413784434157, 0.8506544630476945, 0.8456552228988117, 0.8437212488934422, 0.8444673721579035, 0.8489346132952108, 0.8562699913797784, 0.8628460644450207, 0.8705747378998827, 0.8798826422522013, 0.8878715268791122, 0.8910955115638195, 0.8885065085839569, 0.878600924112645, 0.8604394895975108, 0.8418067341764426, 0.8345480304598061, 0.8332620514793623, 0.833602858847165, 0.8336263152833341, 0.8350507676181034, 0.8378277014377796, 0.8391875166337701, 0.8371440966470945, 0.8358155578562704, 0.8353882151254173, 0.8328303486221449, 0.8262428618558944, 0.8235475789350597, 0.8242935312102251, 0.8246308132327118, 0.8275910317931305, 0.8303725164051157, 0.8308889850281765, 0.8316039947622169, 0.834719699531012, 0.838993427492755, 0.8390768881109494, 0.8362965104453679, 0.8330666349208231, 0.8299184397873355, 0.8282047568224762, 0.825303997292612, 0.8199699894962712, 0.8162552668338598, 0.8145737558688507, 0.8139098805663871, 0.8148758251488513, 0.8134092693800455, 0.8089728658286615, 0.808083382162772, 0.8087134619186508, 0.8078077462443506, 0.8049232738289799, 0.8029497943680002, 0.801880528381652, 0.7991709662366406, 0.7995218478837939, 0.8013012791405869, 0.8031199192320357, 0.8069094297456672, 0.8080217929756764, 0.8078352268234854, 0.807466749046193, 0.8046697011422846, 0.7998219450776499, 0.7966695849446145, 0.7969750767127151, 0.7982973877573808, 0.7990106833980295, 0.7992991950633784, 0.7976848080852477, 0.7935732839622973, 0.7898816079676093, 0.7866794026929892, 0.7846846477452006, 0.7841859304115447, 0.7834273874229597, 0.7817727532030866, 0.7808947750028149, 0.7798289228357806, 0.7790438155023807, 0.7793592538026829, 0.7787887973084548, 0.7785232003436816, 0.7799318981762275, 0.7812137699088911, 0.7813201788248698, 0.7807503115112617, 0.7798197980279181, 0.7782270896536694, 0.7770805535067772, 0.7761511143457509, 0.7744930352482502, 0.7719268080532227, 0.7703832448368306, 0.7691500259621104, 0.7691503952123837, 0.7692129792471293, 0.7692453144531921, 0.767022754889598, 0.7655418992186195, 0.765562436031398, 0.7642055205287734, 0.7622338552599398, 0.7604983308621399, 0.7593037497481343, 0.7586453238861606, 0.7575213727347101, 0.7561231998248032, 0.7555724054854208, 0.7554710376573007, 0.7567856284864799, 0.7582514812072878, 0.7578072336944952, 0.7568861436323562, 0.755527930384004, 0.7549023990529802, 0.7541347535612196, 0.7529800413299966, 0.7514390247368825, 0.7503410144751881, 0.7489308051870556, 0.7483052946035492, 0.7476427904209955, 0.7466688474470857, 0.7460361647536271, 0.7463659923737599, 0.7455045804493902, 0.7441410789860835, 0.7426752921453642, 0.7413070106940058, 0.7403770556846307, 0.7395330723244699, 0.7388308839897171, 0.7385085871441447, 0.7385416985008945, 0.7388191139386929, 0.7387620204957985, 0.7382651922529876, 0.7370553866402915, 0.7360768242885176, 0.7359541302505191, 0.7352976817809268, 0.734721219284192, 0.7344441223190638, 0.7336478512112394, 0.7328557234961695, 0.7325082322051107, 0.7317719337735735, 0.7308200717267419, 0.7296623241419736, 0.7284930494086318, 0.7275456576702668, 0.7268137352557904, 0.7263614826685263, 0.726072027484581, 0.7254210793384004, 0.7246961379986592, 0.7238068841621632, 0.7230376621481818, 0.7223540850999084, 0.7221121023947511, 0.7221957093902738, 0.7219586325442666, 0.7210245729780387, 0.7199217413611687, 0.7195855434007448, 0.7189149484031898, 0.7182819027395299, 0.7181954605255196, 0.7176391418521801, 0.7171946935025276, 0.7167992187687652, 0.7162231083230918, 0.7157318843466322, 0.7150509906683566, 0.7142535587872626, 0.7135221849007694, 0.712853278425168, 0.7122540205841427, 0.7116388803621979, 0.7108355965207243, 0.7100523397843773, 0.7094200238292245, 0.7088789056559905, 0.7084044901563381, 0.7080722036078011, 0.7078529852173275, 0.7072477930122689, 0.706682757214963, 0.7060182661971757, 0.7056899258626186, 0.7054410251193368, 0.7049809155253006, 0.7045147167349487, 0.7038863830524222, 0.7033484833805304, 0.7028300676068358, 0.7024517074424507, 0.7021440922768046, 0.7017913125060402, 0.7012512669320233, 0.7009681455252401, 0.7005144050117228, 0.700020736831075, 0.6994777720243689, 0.6989535337416606, 0.6983434426398842, 0.6976619759542354, 0.6973145856287433, 0.6968057175281062, 0.6965233041596028, 0.6964041296958855, 0.6962226075705865, 0.6958331338168516, 0.6952982878613359, 0.6951359359065189, 0.6946974152417336, 0.6943457255992339, 0.694023210039272, 0.6935514051928146, 0.6933287975453658, 0.6928168552660768, 0.6925401729507527, 0.692202463116474, 0.6919491939040613, 0.6915888576615055, 0.6912893575097607, 0.6910534190758545, 0.6906558815116315, 0.6903485656917634, 0.6899401690673079, 0.6895110978435576, 0.6892491845761757, 0.6888338433784749, 0.6886518069093968, 0.6883175496048627, 0.6880075654657946, 0.6878136979315227, 0.6875421407973263, 0.6873492665088176, 0.6872576014795454, 0.6867971547239288, 0.6867517519171447, 0.6864838204651011, 0.6862423387580147, 0.6865287525607133, 0.6859643105867176, 0.6857043627488543, 0.6855364664202589, 0.6853754499247047, 0.6852047449432731, 0.6850004216896612, 0.6848400745512434, 0.6846787622531247, 0.6845658689150506, 0.684366696217474, 0.6841431615548095, 0.6842334566452962, 0.6838818514208315, 0.6837641115226589, 0.6836775986287774, 0.6836051439271783, 0.6834901424802771, 0.6833520854273702, 0.6832944531754417, 0.6833594558049797, 0.6832213140888492, 0.6830063465157885, 0.6830406992128267, 0.6829547064090815, 0.6830800409823086, 0.6829206113597522, 0.6830801384987769, 0.6828677302441154 ] } ], "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": { "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": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# So, that gave an energy shift of 14 eV. 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 enormous energy shift" ] }, { "cell_type": "code", "execution_count": 9, "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", "# 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": 10, "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, would we notice? Or, how would we know?\n", "\n", "# you 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": 11, "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": 12, "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": "77ec89b8-0c7e-4929-9d6f-95a625a7c43b", "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.5636628523662258e-06, -0.00018962915370151394, -0.00038182028435596125, -0.000958386796344108, -0.001919327257184547, -0.003841173846376017, -0.005762967601727115, -0.007684715787462065, -0.009606478242487195, -0.01154987352668411, -0.01347170670542877, -0.015393507881299283, -0.017315461680400983, -0.019237143481264072 ] } ], "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": [ "# you can see that an error is angular offset is more pronounced at high energy, \n", "# which is lower angle. On the other hand, a bad value for d-spacing gives a \n", "# slightly less energy-dependent error (\"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": 13, "id": "1353430a", "metadata": {}, "outputs": [], "source": [ "# To get an error in R of 0.01 Ang, E0 needs to shift by about 35 eV at Cu K edge. \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 requires\n", "# a very careful analysis.\n", "\n", "# Errors or drifts in energy of 5 eV in the 5 to 15 keV range will cause very\n", "# significant problems for XANES work and would need to be addressed.\n", "#\n", "# Errors of that magnitude will not impact EXAFS results at all.\n", "\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.10.10" } }, "nbformat": 4, "nbformat_minor": 5 }