{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "0", "metadata": { "execution": { "iopub.execute_input": "2026-06-30T22:30:58.163942Z", "iopub.status.busy": "2026-06-30T22:30:58.163727Z", "iopub.status.idle": "2026-06-30T22:30:58.168079Z", "shell.execute_reply": "2026-06-30T22:30:58.167422Z" }, "tags": [ "hide-in-docs" ] }, "outputs": [], "source": [ "# Check whether easydiffraction is installed; install it if needed.\n", "# Required for remote environments such as Google Colab.\n", "import importlib.util\n", "\n", "if importlib.util.find_spec('easydiffraction') is None:\n", " %pip install easydiffraction==0.19.1" ] }, { "cell_type": "markdown", "id": "1", "metadata": {}, "source": [ "# PbSO₄ — powder X-ray CW — round robin\n", "\n", "Verifies the anglesite X-ray round-robin case with the Cu Kα doublet\n", "and FullProf empirical asymmetry.\n", "\n", "**Refinement:** atomic ADPs, the overall scale, and Bérar-Baldinozzi\n", "coefficients (round-robin case); see the known difference below." ] }, { "cell_type": "code", "execution_count": 2, "id": "2", "metadata": { "execution": { "iopub.execute_input": "2026-06-30T22:30:58.170037Z", "iopub.status.busy": "2026-06-30T22:30:58.169843Z", "iopub.status.idle": "2026-06-30T22:31:00.982709Z", "shell.execute_reply": "2026-06-30T22:31:00.981768Z" } }, "outputs": [], "source": [ "import easydiffraction as edi\n", "from easydiffraction import ExperimentFactory\n", "from easydiffraction import StructureFactory\n", "from easydiffraction.analysis import verification as verify" ] }, { "cell_type": "markdown", "id": "3", "metadata": {}, "source": [ "## Build the project" ] }, { "cell_type": "code", "execution_count": 3, "id": "4", "metadata": { "execution": { "iopub.execute_input": "2026-06-30T22:31:00.984534Z", "iopub.status.busy": "2026-06-30T22:31:00.984254Z", "iopub.status.idle": "2026-06-30T22:31:01.195548Z", "shell.execute_reply": "2026-06-30T22:31:01.194661Z" } }, "outputs": [], "source": [ "project = edi.Project()" ] }, { "cell_type": "markdown", "id": "5", "metadata": {}, "source": [ "## Define the structure" ] }, { "cell_type": "code", "execution_count": 4, "id": "6", "metadata": { "execution": { "iopub.execute_input": "2026-06-30T22:31:01.197703Z", "iopub.status.busy": "2026-06-30T22:31:01.197543Z", "iopub.status.idle": "2026-06-30T22:31:01.206517Z", "shell.execute_reply": "2026-06-30T22:31:01.205835Z" } }, "outputs": [], "source": [ "structure = StructureFactory.from_scratch(name='pbso4')\n", "\n", "structure.space_group.name_h_m = 'P n m a' # FullProf Space group symbol\n", "\n", "structure.cell.length_a = 8.485900 # FullProf a\n", "structure.cell.length_b = 5.402259 # FullProf b\n", "structure.cell.length_c = 6.964587 # FullProf c\n", "\n", "structure.atom_sites.create(\n", " id='Pb', # FullProf Atom\n", " type_symbol='Pb', # FullProf Typ\n", " fract_x=0.18822, # FullProf X\n", " fract_y=0.25, # FullProf Y\n", " fract_z=0.16711, # FullProf Z\n", " adp_type='Biso', # FullProf Biso\n", " adp_iso=1.86290, # FullProf Biso\n", ")\n", "structure.atom_sites.create(\n", " id='S', # FullProf Atom\n", " type_symbol='S', # FullProf Typ\n", " fract_x=0.06306, # FullProf X\n", " fract_y=0.25, # FullProf Y\n", " fract_z=0.68485, # FullProf Z\n", " adp_type='Biso', # FullProf Biso\n", " adp_iso=1.34971, # FullProf Biso\n", ")\n", "structure.atom_sites.create(\n", " id='O1', # FullProf Atom\n", " type_symbol='O', # FullProf Typ\n", " fract_x=0.90281, # FullProf X\n", " fract_y=0.25, # FullProf Y\n", " fract_z=0.59724, # FullProf Z\n", " adp_type='Biso', # FullProf Biso\n", " adp_iso=1.19713, # FullProf Biso\n", ")\n", "structure.atom_sites.create(\n", " id='O2', # FullProf Atom\n", " type_symbol='O', # FullProf Typ\n", " fract_x=0.18443, # FullProf X\n", " fract_y=0.25, # FullProf Y\n", " fract_z=0.54586, # FullProf Z\n", " adp_type='Biso', # FullProf Biso\n", " adp_iso=2.56174, # FullProf Biso\n", ")\n", "structure.atom_sites.create(\n", " id='O3', # FullProf Atom\n", " type_symbol='O', # FullProf Typ\n", " fract_x=0.08094, # FullProf X\n", " fract_y=0.02239, # FullProf Y\n", " fract_z=0.81289, # FullProf Z\n", " adp_type='Biso', # FullProf Biso\n", " adp_iso=1.37463, # FullProf Biso\n", ")\n", "\n", "project.structures.add(structure)" ] }, { "cell_type": "markdown", "id": "7", "metadata": {}, "source": [ "## Load the FullProf reference" ] }, { "cell_type": "code", "execution_count": 5, "id": "8", "metadata": { "execution": { "iopub.execute_input": "2026-06-30T22:31:01.208370Z", "iopub.status.busy": "2026-06-30T22:31:01.208201Z", "iopub.status.idle": "2026-06-30T22:31:01.217831Z", "shell.execute_reply": "2026-06-30T22:31:01.217039Z" } }, "outputs": [], "source": [ "FULLPROF_PROJECT_DIR = 'pd-xray-cwl_pbso4_round-robin'\n", "FULLPROF_PRF_FILE = 'pbsox.prf'\n", "FULLPROF_SUM_FILE = 'pbsox.sum'\n", "FULLPROF_BAC_FILE = 'pbsox.bac'\n", "FULLPROF_LABEL = verify.fullprof_label(FULLPROF_PROJECT_DIR, FULLPROF_SUM_FILE)\n", "\n", "FULLPROF_ZERO = 0.00363 # FullProf Zero\n", "FULLPROF_SCALE = 0.0004703142 # FullProf Scale\n", "FULLPROF_WAVELENGTH_1 = 1.540560 # FullProf Lambda1\n", "FULLPROF_WAVELENGTH_2 = 1.544400 # FullProf Lambda2\n", "FULLPROF_WAVELENGTH_2_TO_1_RATIO = 0.50000 # FullProf Ratio\n", "FULLPROF_U = 0.048457 # FullProf U\n", "FULLPROF_V = -0.083053 # FullProf V\n", "FULLPROF_W = 0.035188 # FullProf W\n", "FULLPROF_X = 0.0 # FullProf X\n", "FULLPROF_Y = 0.058360 # FullProf Y\n", "FULLPROF_WDT = 48.0 # FullProf Wdt\n", "FULLPROF_ASY_1 = -0.41356 # FullProf Asy1\n", "FULLPROF_ASY_2 = 0.0 # FullProf Asy2\n", "FULLPROF_ASY_3 = 1.26777 # FullProf Asy3\n", "FULLPROF_ASY_4 = 0.0 # FullProf Asy4\n", "\n", "x, calc_fullprof = verify.load_fullprof_calc_profile(\n", " FULLPROF_PROJECT_DIR,\n", " FULLPROF_PRF_FILE,\n", " FULLPROF_BAC_FILE,\n", " FULLPROF_ZERO,\n", ")" ] }, { "cell_type": "markdown", "id": "9", "metadata": {}, "source": [ "## Create the experiment" ] }, { "cell_type": "code", "execution_count": 6, "id": "10", "metadata": { "execution": { "iopub.execute_input": "2026-06-30T22:31:01.219579Z", "iopub.status.busy": "2026-06-30T22:31:01.219410Z", "iopub.status.idle": "2026-06-30T22:31:02.103844Z", "shell.execute_reply": "2026-06-30T22:31:02.103068Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "⚠️ Switching peak profile type adds these settings with defaults: \n", " • asym_beba_a0=0.0 \n", " • asym_beba_a1=0.0 \n", " • asym_beba_b0=0.0 \n", " • asym_beba_b1=0.0 \n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "\u001b[1;36mPeak profile type for experiment \u001b[0m\u001b[32m'pbso4'\u001b[0m\u001b[1;36m changed to\u001b[0m\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "pseudo-voigt + berar-baldinozzi asymmetry\n" ] } ], "source": [ "experiment = ExperimentFactory.from_scratch(\n", " name='pbso4',\n", " sample_form='powder',\n", " beam_mode='constant wavelength',\n", " radiation_probe='xray',\n", " scattering_type='bragg',\n", ")\n", "verify.set_reference_as_measured(experiment, x, calc_fullprof)\n", "\n", "experiment.linked_structures.create(structure_id='pbso4', scale=FULLPROF_SCALE)\n", "\n", "experiment.instrument.setup_wavelength = FULLPROF_WAVELENGTH_1\n", "experiment.instrument.setup_wavelength_2 = FULLPROF_WAVELENGTH_2\n", "experiment.instrument.setup_wavelength_2_to_1_ratio = FULLPROF_WAVELENGTH_2_TO_1_RATIO\n", "experiment.instrument.calib_twotheta_offset = FULLPROF_ZERO\n", "\n", "experiment.peak.type = 'pseudo-voigt + berar-baldinozzi asymmetry'\n", "experiment.peak.broad_gauss_u = FULLPROF_U\n", "experiment.peak.broad_gauss_v = FULLPROF_V\n", "experiment.peak.broad_gauss_w = FULLPROF_W\n", "experiment.peak.broad_lorentz_x = FULLPROF_X\n", "experiment.peak.broad_lorentz_y = FULLPROF_Y\n", "experiment.peak.asym_beba_a0 = FULLPROF_ASY_1\n", "experiment.peak.asym_beba_b0 = FULLPROF_ASY_2\n", "experiment.peak.asym_beba_a1 = FULLPROF_ASY_3\n", "experiment.peak.asym_beba_b1 = FULLPROF_ASY_4\n", "\n", "# FullProf excludes 0-10 deg and 154-180 deg in the PCR.\n", "experiment.excluded_regions.create(id='1', start=0.0, end=10.0)\n", "experiment.excluded_regions.create(id='2', start=154.0, end=180.0)\n", "\n", "experiment.peak.cutoff_fwhm = FULLPROF_WDT\n", "\n", "project.experiments.add(experiment)" ] }, { "cell_type": "markdown", "id": "11", "metadata": {}, "source": [ "## edi-cryspy VS FullProf" ] }, { "cell_type": "code", "execution_count": 7, "id": "12", "metadata": { "execution": { "iopub.execute_input": "2026-06-30T22:31:02.105568Z", "iopub.status.busy": "2026-06-30T22:31:02.105342Z", "iopub.status.idle": "2026-06-30T22:31:02.703737Z", "shell.execute_reply": "2026-06-30T22:31:02.703026Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\u001b[1;36mCalculator for experiment \u001b[0m\u001b[32m'pbso4'\u001b[0m\u001b[1;36m already set to\u001b[0m\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "cryspy\n" ] }, { "data": { "text/html": [ "
Loading plot…
" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "experiment.calculator.type = 'cryspy'\n", "\n", "project.analysis.calculate()\n", "calc_ed_cryspy = experiment.data.intensity_calc\n", "LABEL_ED_CRYSPY = verify.engine_label('cryspy')\n", "\n", "project.display.pattern_comparison(\n", " 'pbso4',\n", " reference=verify.restrict_to_included(experiment, calc_fullprof),\n", " candidate=calc_ed_cryspy,\n", " reference_label=FULLPROF_LABEL,\n", " candidate_label=LABEL_ED_CRYSPY,\n", ")" ] }, { "cell_type": "markdown", "id": "13", "metadata": {}, "source": [ "## Fit edi-cryspy to FullProf" ] }, { "cell_type": "code", "execution_count": 8, "id": "14", "metadata": { "execution": { "iopub.execute_input": "2026-06-30T22:31:02.708111Z", "iopub.status.busy": "2026-06-30T22:31:02.707918Z", "iopub.status.idle": "2026-06-30T22:31:24.749184Z", "shell.execute_reply": "2026-06-30T22:31:24.748480Z" } }, "outputs": [ { "data": { "text/html": [], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "\n", "(function() {\n", " const button = document.getElementById('ed-fit-stop-0cb04812138a41498821eb953180d800-button');\n", " const status = document.getElementById('ed-fit-stop-0cb04812138a41498821eb953180d800-status');\n", " const kernelId = '';\n", " if (!button) {\n", " return;\n", " }\n", "\n", " function setStatus(text) {\n", " if (status) {\n", " status.textContent = text;\n", " }\n", " }\n", "\n", " function pageConfig() {\n", " const element = document.getElementById('jupyter-config-data');\n", " if (!element || !element.textContent) {\n", " return {};\n", " }\n", " try {\n", " return JSON.parse(element.textContent);\n", " } catch (error) {\n", " return {};\n", " }\n", " }\n", "\n", " function baseUrl(config) {\n", " const configured = config.baseUrl || config.base_url ||\n", " (window.Jupyter && Jupyter.notebook && Jupyter.notebook.base_url);\n", " if (configured) {\n", " return configured.endsWith('/') ? configured : configured + '/';\n", " }\n", " const markers = ['/lab/', '/notebooks/', '/tree/'];\n", " for (const marker of markers) {\n", " const index = window.location.pathname.indexOf(marker);\n", " if (index >= 0) {\n", " return window.location.pathname.slice(0, index + 1);\n", " }\n", " }\n", " return '/';\n", " }\n", "\n", " function token(config) {\n", " return config.token || new URLSearchParams(window.location.search).get('token') || '';\n", " }\n", "\n", " function cookie(name) {\n", " const prefix = name + '=';\n", " for (const part of document.cookie.split(';')) {\n", " const trimmed = part.trim();\n", " if (trimmed.startsWith(prefix)) {\n", " return decodeURIComponent(trimmed.slice(prefix.length));\n", " }\n", " }\n", " return '';\n", " }\n", "\n", " function notebookPath() {\n", " const decoded = decodeURIComponent(window.location.pathname);\n", " const markers = ['/lab/tree/', '/notebooks/', '/tree/'];\n", " for (const marker of markers) {\n", " const index = decoded.indexOf(marker);\n", " if (index >= 0) {\n", " return decoded.slice(index + marker.length);\n", " }\n", " }\n", " return '';\n", " }\n", "\n", " async function kernelFromSessions(config) {\n", " const url = new URL(baseUrl(config) + 'api/sessions', window.location.origin);\n", " const authToken = token(config);\n", " if (authToken) {\n", " url.searchParams.set('token', authToken);\n", " }\n", " const response = await fetch(url, {credentials: 'same-origin'});\n", " if (!response.ok) {\n", " return '';\n", " }\n", " const sessions = await response.json();\n", " const path = notebookPath();\n", " const session = sessions.find((item) => item.path === path) || sessions[0];\n", " return session && session.kernel ? session.kernel.id : '';\n", " }\n", "\n", " async function interruptKernel(config, resolvedKernelId) {\n", " const url = new URL(\n", " baseUrl(config) + 'api/kernels/' + resolvedKernelId + '/interrupt',\n", " window.location.origin\n", " );\n", " const authToken = token(config);\n", " if (authToken) {\n", " url.searchParams.set('token', authToken);\n", " }\n", " const xsrfToken = cookie('_xsrf');\n", " const headers = {};\n", " if (xsrfToken) {\n", " headers['X-XSRFToken'] = xsrfToken;\n", " }\n", " const response = await fetch(url, {\n", " method: 'POST',\n", " credentials: 'same-origin',\n", " headers: headers\n", " });\n", " return response.ok;\n", " }\n", "\n", " button.addEventListener('click', async function() {\n", " button.disabled = true;\n", " setStatus('Stopping...');\n", " const config = pageConfig();\n", " try {\n", " const resolvedKernelId = kernelId || await kernelFromSessions(config);\n", " if (!resolvedKernelId) {\n", " throw new Error('Could not resolve the current kernel id.');\n", " }\n", " const interrupted = await interruptKernel(config, resolvedKernelId);\n", " if (!interrupted) {\n", " throw new Error('Jupyter Server rejected the interrupt request.');\n", " }\n", " setStatus('Interrupt sent...');\n", " } catch (error) {\n", " button.disabled = false;\n", " setStatus('Use Kernel > Interrupt to stop this fit.');\n", " }\n", " });\n", "})();\n" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "\u001b[1;36mStandard fitting\u001b[0m\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "📋 Using experiment 🔬 \u001b[32m'pbso4'\u001b[0m for \u001b[32m'single'\u001b[0m fitting\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "🚀 Starting fit process with \u001b[32m'lmfit \u001b[0m\u001b[32m(\u001b[0m\u001b[32mleastsq\u001b[0m\u001b[32m)\u001b[0m\u001b[32m'\u001b[0m\u001b[33m...\u001b[0m\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "📈 Goodness-of-fit progress:\n" ] }, { "data": { "text/html": [ "
iterationtime (s)χ²change / status
110.23236207.78
2123.2262567.9273.5% ↓
3215.821188.6898.1% ↓
4308.191078.969.2% ↓
54313.621078.96
65518.651078.96
75819.811078.96
" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "🏆 Best goodness-of-fit \u001b[1m(\u001b[0mreduced χ²\u001b[1m)\u001b[0m is \u001b[1;36m1078.96\u001b[0m at iteration \u001b[1;36m57\u001b[0m\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "✅ Fitting complete.\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "⚙️ Settings used:\n" ] }, { "data": { "text/html": [ "
NameValueDescription
1max_iterations1000Maximum solver iterations.
" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "📋 Least-squares fit results:\n" ] }, { "data": { "text/html": [ "
MetricValue
1🧪 Minimizerlmfit (leastsq)
2✅ Overall statussuccess
3⏱️ Fitting time (seconds)19.81
4🔁 Iterations55
5📏 Goodness-of-fit (reduced χ²)1078.96
6📏 R-factor (Rf, %)7.81
7📏 R-factor squared (Rf², %)3.66
8📏 Weighted R-factor (wR, %)3.66
" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "📈 Refined parameters:\n" ] }, { "data": { "text/html": [ "
datablockcategoryentryparameterunitsstartvalues.u.change
1pbso4atom_sitePbadp_isoŲ1.86293.37820.009581.34 % ↑
2pbso4atom_siteSadp_isoŲ1.34972.55130.048689.03 % ↑
3pbso4atom_siteO1adp_isoŲ1.19713.46520.1333189.46 % ↑
4pbso4atom_siteO2adp_isoŲ2.56175.36100.1345109.27 % ↑
5pbso4atom_siteO3adp_isoŲ1.37463.26840.0849137.76 % ↑
6pbso4linked_structurepbso4scale0.00050.00080.000069.25 % ↑
7pbso4peakasym_beba_a0-0.41360.46100.0066211.46 % ↓
8pbso4peakasym_beba_a11.2678-1.36880.0140207.97 % ↓
" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "
• start = parameter value before refinement
• value = refined value from least-squares minimization
• s.u. = standard uncertainty (one sigma), from the covariance matrix
• change = relative change from start, in %; ↑ = increase, ↓ = decrease
" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "
Loading plot…
" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# The anomalous-dispersion table discrepancy is partially compensated by\n", "# scale and isotropic ADPs. The asymmetry terms are refined separately\n", "# because cryspy and FullProf use different Berar-Baldinozzi conventions.\n", "for atom_id in ('Pb', 'S', 'O1', 'O2', 'O3'):\n", " atom = structure.atom_sites[atom_id]\n", " atom.adp_iso.free = True\n", "\n", "experiment.linked_structures['pbso4'].scale.free = True\n", "experiment.peak.asym_beba_a0.free = True\n", "experiment.peak.asym_beba_a1.free = True\n", "\n", "project.analysis.fit()\n", "project.display.fit.results()\n", "\n", "project.analysis.calculate()\n", "calc_ed_cryspy_refined = experiment.data.intensity_calc\n", "LABEL_ED_CRYSPY_REFINED = verify.engine_label('cryspy', note='refined')\n", "\n", "project.display.pattern_comparison(\n", " 'pbso4',\n", " reference=verify.restrict_to_included(experiment, calc_fullprof),\n", " candidate=calc_ed_cryspy_refined,\n", " reference_label=FULLPROF_LABEL,\n", " candidate_label=LABEL_ED_CRYSPY_REFINED,\n", ")" ] }, { "cell_type": "code", "execution_count": 9, "id": "15", "metadata": { "execution": { "iopub.execute_input": "2026-06-30T22:31:24.761841Z", "iopub.status.busy": "2026-06-30T22:31:24.761528Z", "iopub.status.idle": "2026-06-30T22:31:24.766666Z", "shell.execute_reply": "2026-06-30T22:31:24.765726Z" } }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "experiment.linked_structures['pbso4'].scale" ] }, { "cell_type": "markdown", "id": "16", "metadata": {}, "source": [ "## Agreement check" ] }, { "cell_type": "code", "execution_count": 10, "id": "17", "metadata": { "execution": { "iopub.execute_input": "2026-06-30T22:31:24.768993Z", "iopub.status.busy": "2026-06-30T22:31:24.768771Z", "iopub.status.idle": "2026-06-30T22:31:24.780664Z", "shell.execute_reply": "2026-06-30T22:31:24.779900Z" } }, "outputs": [ { "data": { "text/html": [ "
ComparisonMetricExpectedActualOK
1edi 0.19.1 (cryspy 0.12.1) vs FullProf 7.95Profile diff (%)< 2.554.11❌
2Max deviation (%)< 671.43❌
3Area ratio0.99 to 1.010.7406❌
4Shape correlation> 0.9990.8523❌
5edi 0.19.1 (cryspy 0.12.1, refined) vs FullProf 7.95Profile diff (%)< 2.53.66❌
6Max deviation (%)< 61.66✅
7Area ratio0.99 to 1.010.9509❌
8Shape correlation> 0.9990.9994✅
" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "
• Known discrepancy = Laboratory X-ray PbSO4: the remaining strict-tolerance difference is partially compensated by refining scale and isotropic ADPs, consistent with a Cu Kalpha anomalous-dispersion table mismatch in f-prime/f-double-prime, especially for Pb. The Berar-Baldinozzi asymmetry terms are refined separately because cryspy and FullProf use different conventions.
" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "True" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "verify.assert_patterns_agree(\n", " [\n", " (\n", " f'{LABEL_ED_CRYSPY} vs {FULLPROF_LABEL}',\n", " verify.restrict_to_included(experiment, calc_fullprof),\n", " calc_ed_cryspy,\n", " ),\n", " (\n", " f'{LABEL_ED_CRYSPY_REFINED} vs {FULLPROF_LABEL}',\n", " verify.restrict_to_included(experiment, calc_fullprof),\n", " calc_ed_cryspy_refined,\n", " ),\n", " ],\n", " known_discrepancy=True,\n", " reason=(\n", " 'Laboratory X-ray PbSO4: the remaining strict-tolerance '\n", " 'difference is partially compensated by refining scale and '\n", " 'isotropic ADPs, consistent with a Cu Kalpha '\n", " 'anomalous-dispersion table mismatch in f-prime/f-double-prime, '\n", " 'especially for Pb. The Berar-Baldinozzi asymmetry terms are '\n", " 'refined separately because cryspy and FullProf use different '\n", " 'conventions.'\n", " ),\n", ")" ] } ], "metadata": { "jupytext": { "cell_metadata_filter": "-all", "main_language": "python", "notebook_metadata_filter": "-all" }, "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.14.5" } }, "nbformat": 4, "nbformat_minor": 5 }