{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Calculate and Plot Activity\n",
    "\n",
    "\n",
    "Given an existing database for Al-Zn, we would like to calculate the activity of the liquid.\n",
    "\n",
    "## Experimental activity results\n",
    "\n",
    "In order to make sure we are correct, we'll compare the values with experimental results.\n",
    "Experimental activities are digtized from Fig 18 in A. Yazawa, Y.K. Lee, Thermodynamic Studies of the Liquid Aluminum Alloy Systems, Trans. Japan Inst. Met. 11 (1970) 411–418. https://doi.org/10.2320/matertrans1960.11.411.\n",
    "\n",
    "The measurements at at 1073 K and they used a reference state of the pure Zn at that temperature.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "exp_x_zn =   [0.0482, 0.1990, 0.3550, 0.5045, 0.6549, 0.8070, 0.9569]\n",
    "exp_acr_zn = [0.1154, 0.3765, 0.5411, 0.6433, 0.7352, 0.8384, 0.9531]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Set up the database\n",
    "\n",
    "Al-Zn database is taken from S. Mey, Reevaluation of the Al-Zn system, Zeitschrift F{ü}r Met. 84 (1993) 451–455."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "from pycalphad import Database, equilibrium, variables as v\n",
    "import numpy as np\n",
    "\n",
    "dbf = Database('alzn_mey.tdb') \n",
    "\n",
    "comps = ['AL', 'ZN', 'VA']\n",
    "phases = list(dbf.phases.keys())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Calculate the reference state\n",
    "\n",
    "Because all chemical activities must be specified with a reference state, we're going to choose a reference state as the pure element at the same temperature, consistent with the experimental data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "ref_eq = equilibrium(dbf, comps, phases, {v.P: 101325, v.T: 1023, v.X('ZN'): 1})"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Calculate the equilibria\n",
    "\n",
    "Do the calculation over the composition range"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "eq = equilibrium(dbf, comps, phases, {v.P: 1013325, v.T: 1023, v.X('ZN'): (0, 1, 0.005)})"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Get the chemical potentials and calculate activity\n",
    "\n",
    "We need to select the chemical potentials from the xarray Dataset and calculate the activity."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "chempot_ref = ref_eq.MU.sel(component='ZN').squeeze()\n",
    "chempot = eq.MU.sel(component='ZN').squeeze()\n",
    "\n",
    "acr_zn = np.exp((chempot - chempot_ref)/(8.315*1023))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Plot the result"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.legend.Legend at 0x11352f5c0>"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAARQAAAEWCAYAAACnuGhyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd4VGX2wPHvAQKhoxQFQhNCJxCIVEVQlOKCooiABRtY\n186uuFZ2fzYs69pRESvFQlkFUVgQUanSQXpLQkmA0FJIOb8/7hCHYUICzM1MJufzPPOQufedO2eG\nmTP3vve95xVVxRhjAqFEsAMwxoQPSyjGmICxhGKMCRhLKMaYgLGEYowJGEsoxpiAsYQS4kTkBhH5\noYBt14hIN5dDQhwficgBEVnk9vOZosMSiotEZK7nS1emgO3ri4iKSKnjy1T1c1W9oiCPV9UWqjrX\ns61nROSzMwo8fxcBlwNRqtred6UnsR3xuWWISE6gAxGRbiISn0+b7iIyR0QOisg2P+vre9anisgf\nItLDa927fl7HYa/1n4nILhE5JCIbROSOvGITkdIi8o2I/CIilc76xYcgSyguEZH6wMWAAv2CGkzg\n1QO2qepRfys9ia3C8RtwPrAF+GdhBunlKDAWGJHH+vHAMqAq8A/gKxGpDqCqd/m8lvHAl16PfR6o\nr6qVcP6f/yUi7XyfwPOj8g1QBbhCVQ8F5qWFGFW1mws34CngF+BV4FufdWWBV4DtwEFgvmfZDpwE\ndMRz6wTcAsz3PO5d4GWfbU0FHvb8vQ3oAfQCjgGZnu2sAK4Dlvo89hFgSh7x1wKmAfuBTcAwz/Lb\ngXQg27PtZwvwXkwAfgBKeO53A+I9z78X2AXceorH3wqsAw7jJKY7PcvLA2lAjtd7VusU2+mBkwi9\nlzUGMoCKXst+Bu7y8/jynhguyWP7TTyvZaDP6yznef0zgbLB/my6+rkPdgDhevN8Ce8B2nm+2Od5\nrXsLmAvUBkoCnYEyQH1PQinl1dY7oXQFdgLiuX+O5wtVy3N/G9DD8/czwGde2ynjSQ7NvJYtA67N\nI/6fgLeBSKANkARc5htTAd6H+z0xV/Na1g3IAkYBEUAfIBU4J49tXAk0BAS4xNO2rde24gsYi7+E\n0h9Y57PsTeANP4+/2ZPQxGf5256YFPgdqOAVW5LnvZwGlAn259Ltmx3yuEBELsI5LJikqkuBzcAQ\nz7oSwG3AA6qaoKrZqvqrqmYUYNM/43xoL/bcHwD8pqqJ+T3Qs/2JwI2eOFrgJLBv/cRfB6ef5O+q\nmq6qy4EPgJsKEKP3djoCzwHXqWqyz+pMYJSqZqrqdJy9iyZ5xP6dqm5Wx084v/YX+2t7Birg7CV6\nOwhU9NN2KPCJerKFV3z3eNpfjHNY4/1/WRFnT/PjAv4fF2mWUNwxFPjB60v0hWcZQDWcX/3Np7tR\nzwd5AjDYs2gI8PlpbOJjYIiICE5ymJTHh7wWsF9VD3st246zR1UgIlINp69hpKou8NNkn6pmed1P\nxfly+9tWbxFZICL7RSQFZ4+mWkFjyccRwLeDtBLOoY13DHVw9o4+8bcRzw/DfCAKuNtrVTIwCPhY\nRHoGKOaQZQklwESkLDAQuEREdovIbuAhoLWItMb5gKXj7ML7Ksil3+OBASJSD+gAfJ1Hu5O25fli\nH8P5JR0CfJrHYxOBc0XE+1e6LpBQgPiO74V9Afyiqm8U5DGn2FYZnNf4Ms5hYxVgOs7hDxTsPTuV\nNcAFPq+1tWe5t5uBX1V1Sz7bK4XP/62qfgMMw+ns7X6W8YY0SyiBdzVOh2VznL6HNkAznMOVm1U1\nB+eMw6siUktESopIJ88XJwmng/GCvDauqss87T4AZqpqSh5N9wD1PV9ub5/g9BFkeX5R/T3HTuBX\n4HkRiRSRGJzO2ILuDT0D1AHuyKddQZTG6f9JArJEpDfgfRp9D1BVRCrntQERKSEikTj9NeJ5TaUB\nVHUDsBx42rO8PxDDyYn6ZmCcz3ZriMggEang+X/sibP3+D/fGFR1PHAfMFVEupzG6y9agt2JE243\n4HvgFT/LBwK7cX7BygL/xvnFPwjMw9P7j9NRmQSkAB3x0wEKPInzy3ydz/Jt/NkpWxXn7NEB4Hev\nNnVxktYpz87g7Lp/i9ORuxmvsx7+YvJ5bA5OP8IRP7e6+OlI9Y7dz/buxUkcKTh7VROAf3mtHwvs\n86w/6SyP5/nU5zbXa319nE7yNGC9bxw4fSBH8ToT5FleHafDNQU4BKzCczbM63l9X+cwT/v2wf6s\nunE7frbAFBOeQ7K9OGdJNgY7HhNe7JCn+LkbWGzJxLjBtYQiImNFZK+IrM5jvYjIf0Rkk4isFJG2\nbsViHJ5h5w/gDCgzJuDc3EMZhzNiMy+9gWjPbTjwjouxGEBV66tqPXU6do0JONcSiqrOw+nQy8tV\neAYJqXM6s4qI1HQrHmOM+0rl38Q1tXGGZB8X71m2y7ehiAzH2YuhfPny7Zo2bVooARpTnChwMPUY\nVcqVZunSpcmqWv10txHMhCJ+lvk95aSqY4AxAHFxcbpkyRI34zKm2MnIyua+L5bx49o9vH1nJzpc\nUHX7mWwnmAklHmfw03FROCM0jTGFKD0zm7s/W8qc9Uk8268F7Ruce8bbCuZp42nAzZ6zPR2Bg6p6\n0uGOMcY9aceyGfbJEuasT+L/+rdkaOf6Z7U91/ZQRGQ8zkjBap6qVU/jDH1GVd/FuR6jD85l/qk4\nNS+MMYUk9VgWt49bwoKt+3hpQAwD4+rk/6B8uJZQVHVwPusVZ0j1WcvMzCQ+Pp709PRAbM4EQGRk\nJFFRUURERAQ7FOPHwbRMbh+3mN93HODVga3pHxsVkO0Gsw8lYOLj46lYsSL169fHuTLfBJOqsm/f\nPuLj42nQoEGwwzE+ko9kcPOHi9i49zBvDG7LlTGBG60RFkPv09PTqVq1qiWTECEiVK1a1fYYQ1Bi\nShoD3/2NLclHeP/muIAmEwiTPRTAkkmIsf+P0DJlWQLPTV/H3sMZCHDfpY3o1qRGwJ8nLPZQjDF5\nm7Isgb9/vZK9h53ifAp88PNWpiwrUL2s02IJJYB2797NoEGDaNiwIc2bN6dPnz5s2LAhz/YVKvit\neJivuXPn8pe//OWUbZYvX8706dNPe9vdunXDBg6Gl399t5aMrBOnRErLzGb0zPUBfy5LKAGiqvTv\n359u3bqxefNm1q5dy3PPPceePXuCEs+ZJhQTXuau30vykWN+1yWmpAX8+SyhBMicOXOIiIjgrrvu\nyl3Wpk0bYmNjueyyy2jbti2tWrVi6tSpfh//0ksv0apVK1q3bs1jjz0GnLi3kJycTP369U963KJF\ni+jcuTOxsbF07tyZ9evXc+zYMZ566ikmTpxImzZtmDhxIkePHuW2227jwgsvJDY2NjeOtLQ0Bg0a\nRExMDNdffz1paYH/kJng+HppPHd8vISIEv77s2pVKRvw5wybTtnjnv3vGtYmBnZStua1KvF03xan\nbLN69WratTtpwjgiIyOZPHkylSpVIjk5mY4dO9KvX78TOi1nzJjBlClTWLhwIeXKlWP//lNdpH2i\npk2bMm/ePEqVKsWsWbN4/PHH+frrrxk1ahRLlizhzTffBODxxx/n0ksvZezYsaSkpNC+fXt69OjB\ne++9R7ly5Vi5ciUrV66kbVsrS1PUqSrv/rSFF7//gy6NqtI3phbP/nctaZnZuW3KRpRkRE+/s5ac\nlbBLKKFGVXn88ceZN28eJUqUICEhgT179nD++efntpk1axa33nor5cqVA+Dccwt+LcXBgwcZOnQo\nGzduRETIzMz02+6HH35g2rRpvPzyy4Bzqn3Hjh3MmzeP+++/H4CYmBhiYmLO9KWaEJCdo/zz27WM\n+3Ub/VrX4uXrWlO6VAkiI0oyeuZ6ElPSqFWlLCN6NuHq2ALPilJgYZdQ8tuTcEuLFi346quvTlr+\n+eefk5SUxNKlS4mIiKB+/fonjc9QVb+nWUuVKkVOjtOZlteYjieffJLu3bszefJktm3bRrdu3fy2\nU1W+/vprmjQ5+VfJTvGGh/TMbB6ZtILvVu3ijosa8HifZpTwHO5cHVvblQTiy/pQAuTSSy8lIyOD\n999/P3fZ4sWL2b59OzVq1CAiIoI5c+awffvJV4VfccUVjB07ltTUVIDcQ5769euzdOlSAL/JCpw9\nlNq1nQ/KuHHjcpdXrFiRw4f/nKuqZ8+evPHGG8crr7NsmVO0rWvXrnz+uTM7xurVq1m5cuUZvX4T\nXAfTMhk6dhHfrdrFP/o044m/NM9NJoXJEkqAiAiTJ0/mxx9/pGHDhrRo0YJnnnmGPn36sGTJEuLi\n4vj888/xVxyqV69e9OvXj7i4ONq0aZN7WPLoo4/yzjvv0LlzZ5KTfWfydPztb39j5MiRdOnShezs\nP4+Ru3fvztq1a3M7ZZ988kkyMzOJiYmhZcuWPPnkkwDcfffdHDlyhJiYGF566SXat2/vwrtj3BR/\nIJWB7/7G7zsO8PqgNgzrmue0Tq4rctNo+CuwtG7dOpo1axakiExe7P/FfSt2pnD7x0vIyMrm3Rvb\n0aVRYGZoFZGlqhp3uo8Luz4UY4qL71fv4sGJy6lWoQzjh3Ug+jx/87sXLksoxhQxqsr7P2/h+Rl/\n0KZOFd6/OY5qFcoEOyzAEooxRUpmdg5PTV3D+EU7uDKmJq9c15rIiJLBDiuXJRRjiohD6Znc+/nv\n/LwxmXu7N+SRy5sE5UzOqVhCMaYI2Lk/lds/XsyWpKMBK9foBksoxoS43zbv457Pl5Kdo3xyW3s6\nB+hMjhtsHEoAqCoXXXQRM2bMyF02adIkevU61UyshWvv3r1069aN8uXL8+CDD56w7osvvqBVq1a0\naNGCkSNH5i6///77adOmDW3atCE6Oppq1f78IH/44YdER0cTHR3NZ599Vmivo7j5dMF2bvpwIVUr\nlGHqfReFdDIBnC9DUbq1a9dOfa1du/akZYVt1apV2rRpU01LS9MjR45oo0aNdNOmTcEOK9fhw4d1\n/vz5+sYbb+gDDzyQu3zPnj1at25dTU5O1pycHB0yZIjOnTv3pMe/+uqrOmzYMFVVTUpK0gsuuEAP\nHDigycnJWr9+fU1JSTnpMaHw/1JUZWRm68hvVmq9v3+rt360SA+lHSvU5weW6Bl8P4vlHsqUZQl0\neeF/NHjsO7q88L+AVK5q2bIlffv25cUXX+TZZ5/l5ptvpmHDhvTt25d27drRokULPvjgAwAmT56c\n+8vfuHFjoqOjAXj66ae58MILadmyJXfddReqSmJiIh06dABg6dKliAiJic58aA0aNCA9PZ2pU6fS\noUMHYmNjueKKK9i7d+9J8VWoUIEuXboQGRl5wvLNmzfTrFmz3Jq8PXr04Ouvvz7p8ePHj2fwYGci\ngxkzZtCrVy+qVKlC1apVufTSS/nhhx/O+j00jn1HMrjxg4V8sXAHd3dryPs3x1ExsmjMHlDs+lCm\nLEtg5Derci/lTkhJY+Q3qwDO+uKpp59+mrZt21K6dOncOiYff/wx5557LqmpqcTFxXHttdfSv39/\n+vfvj6oyYMAArrjiCgAeeOABnn32WVSVIUOG8P3339O7d28OHjzI0aNH+fnnn4mLi8v9NyoqisjI\nSLp27ZpbEuHdd9/llVde4cUXXyxQzNHR0axevZodO3ZQs2ZNpk6detLFglu2bCEhIYFLLrnEec8S\nEqhT589OwaioKBISAl9OsDham3iIYZ8sIflIBq8PasNVbdy/oC+Qil1CGT1z/Ql1IeDPcnhnm1DK\nly/P9ddfT4UKFShTxhlo9NprrzFt2jTAme5j8+bNxMU5I5qff/55KleuzJ133gnA7NmzGT16NOnp\n6SQnJ9OuXTt69+5Np06d+PXXX/n55595/PHHmTVrFmlpaVx88cUA7Nixg4EDB7J7924yMjJo3Lhx\ngWOuVq0ab731FgMGDKBUqVJ06NCB+Pj4E9qMHz+egQMHUqKEs0Orfq6OtiuWz970Vbt4ZNIKKpeN\n4Mu7OhETVSXYIZ22YnfIk1fZu0CVwytRokTuF2/WrFnMmzePBQsWsGLFCmJiYnLLEMycOZOpU6fy\n9ttvA5Camsp9993H5MmTWblyJbfddltu24svvph58+aRkJBA3759WbZsGfPnz6dr164A3HvvvTz0\n0EOsWrWKt99++7Snr7jqqqtYtGgRv/76a25Hq7cJEybkHu6As0eyc+fO3Pvx8fHUqlXrNN8pc1xW\ndg7Pz1jHPZ//TtOaFZn21y5FMplAMUwoeZW9c6Mc3sGDBzn33HMpW7Ysa9asYfHixQBs3bqV+++/\nny+//DK3TyMtLY0SJUpQrVo1Dh8+fEI/RteuXfn4449p2rQppUqVomLFivzwww907tw593lq166N\nqvLxxx+fdpzH+1z279/Pu+++y+233567bs2aNaSlpZ1wFXKvXr2YMWMGKSkp7Nu3j9mzZ+cetpnT\nk3wkg5s+XMR7P23hhg51mTC8IzUqRub/wBBV7A55RvRsckIfCrhXDu/KK69kzJgxtG7dmqZNm+Z2\nro4dO5b9+/fTr18/AOrWrcu0adMYOnQoLVu2pF69erltARo1akRWVlbuHkmXLl1ISkqiUqVKADzz\nzDP079+fqKgo2rdvz65d/uecj4qKIjU1lczMTL766itmz55NkyZNuPfee1mzZk3utho2bJj7mPHj\nxzNo0KATtlO9enVGjhyZe+g2atQoKleuHIi3LGxNWZZwUsW0elXLcc/nv7P/6DFGD4jhuhAdrHY6\nimX5An//uYVRzaq4sfIFDt8TAQARJYQcoFaVSN65oR0ta4dWQrbyBaehsMrhGQP+TwRk5ihlSpXg\nv/ddRJVypYMUWeAVuz4UYwpbXh3+GVk5YZVMIIwSSlE7dAt39v/xp7w6/Gu7cCIg2MIioURGRrJv\n3z77EIcIVWXfvn0njcotjjKzc2h6/smV1Nw6ERBsYdGHEhUVRXx8PElJScEOxXhERkYSFRUV7DCC\nKjEljb+OX8bS7Qfo3LAq25KPsutgelifCAiLhBIREUGDBg2CHYYxuf73xx4enrSCzKwc/jM4ln6t\ni8fAv7BIKMaEiszsHF6euZ735m2hWc1KvH1DWxpUKx/ssAqNq30oItJLRNaLyCYReczP+roiMkdE\nlonIShHp42Y8xrgpISWN69/7jffmOaNeJ9/TuVglE3BxD0VESgJvAZcD8cBiEZmmqmu9mj0BTFLV\nd0SkOTAdqO9WTMa4Zfa6PTzy5QqyspU3BsfSt5gc4vhy85CnPbBJVbcAiMgE4CrAO6EoUMnzd2Ug\n0cV4jAk470Oc5jUr8VYxO8Tx5WZCqQ3s9LofD3TwafMM8IOI/BUoD/TwtyERGQ4MB+e6F2NCwfZ9\nR7l/wnJW7Ezhxo51eeLK5iE1pUUwuNmH4q9Ahu9AkcHAOFWNAvoAn4rISTGp6hhVjVPVuOrVq7sQ\nqjGn55vf4+nz+s9sTTrCW0Pa8q+rWxX7ZALu7qHEA96XT0Zx8iHN7UAvAFX9TUQigWrAyTUMjQkB\nh9IzeWrKaqYsT6R9/XN5bVCbsBzxeqbcTCiLgWgRaQAkAIOAIT5tdgCXAeNEpBkQCdjoNBOSft9x\ngAcmLCMxJZ2HL2/Mvd0bUTLEJtoKNtcSiqpmich9wEygJDBWVdeIyCicitrTgEeA90XkIZzDoVvU\nxs+bEJOdo7wzdxOvzdpIzcqRTLqzI+3qnRvssEKSqwPbVHU6zqlg72VPef29FujiZgzGnI3ElDQe\nnLicRVv307d1Lf6vf0sqFZEK9MFgI2WNycP3q3fx969XOaeGr2vNtW1rWzHufFhCMcZH6rEs/vnt\nWsYv2klMVGVeHxRbrMeWnA5LKMZ4WbbjAA9PWsG2fUe565KGPHx5Y0qXCosqH4XCEooxOCNe3/zf\nJt6cs4nzKpbhizs60qlh1WCHVeRYQjHF3uakIzw8cTkr4g9yTWxtnu7XgsplreP1TFhCMWHtVDMc\nqCqfLdzB/323lsiIkrw1pC1XxtQMcsRFmyUUE7ZONY9154ZV+dvXK5m7PomLo6vx8nWtOa+Slaw8\nW5ZQTNjKax7rUf9di6KkHstm1FUtuKljPTsdHCCWUEzYymv6iv2px4iJqsyrA9vQqEaFQo4qvFlC\nMWGrVpWyJPhJKhXLlOLruzsTUdJOBweavaMmbI3o2YRInzEkpUuW4J9Xt7Rk4hJ7V03YalazEudW\n+HNmvpqVI3lpQExYTl8RKuyQx4SdnBxl7C9been79VQqW4oPh8ZxWbPzgh1WsWAJxYSVhJQ0Hp20\ngt+27OPy5ufxwjWtqFqhTLDDKjYsoZiwoKpMXZ7Ik1NXk5OjvHRtDNfFRdnp4EJmCcUUeSmpx/jH\nlNV8t3IXcfXO4dWBbahbtVywwyqWLKGYIu2nDUmM+HIFB1KPMaJnE+66pKGVZQwiSyimSEo7ls3z\nM9bxyW/bia5RgbG3XEjL2pWDHVaxZwnFFDnLd6bw8MTlbEk+yh0XNeDRnk1sCosQYQnFFBkn1SwZ\n1oHODasFOyzjxRKKKRKsZknRYAnFhDRV5dMF23lu+joiI0ry9g1t6dPKapaEKksoJmTtPpjOiK9W\n8PPGZC5pXJ3RA2KoYTVLQpolFBOS/rsikSemrOZYVg7/urolN3Soa4PUigBLKCakHEzN5Klpq5m6\nPJHWdarw2sDWXFDdapYUFZZQjKtOVdPV1/yNyTz65QqSjmTw8OWNuadbQ0pZmYEixRKKcc2parp6\nJ5X0zGxe/P4PPvplGw2rl2fMzZ2JiaoSlJjN2bGEYlyTV03X0TPX5yaUP3Yf4oHxy1m/5zBDO9Xj\nsd7NKFvaBqkVVZZQjGvyqumamJKGqvLRL9t44fs/qBQZwbhbL6RbkxqFHKEJNEsoxjV51XQ9r1Ik\nt3y0mJ82JHFZ0xq8NCDGapaECevxMq4Z0bMJZX2usSldsgRHMrJYsGUf/7y6JR8MjbNkEkZsD8W4\n5ng/yeiZ60lISaN86ZIcPZZNwxoV+M+gNkSfVzHIEZpAs4RiXHV1bG0an1eR+ycsY9PeI9xxUQNG\n9GpCmVLW8RqOLKEY13gXi65cLoJPbmtP18bVgx2WcZGrfSgi0ktE1ovIJhF5LI82A0VkrYisEZEv\n3IzHFJ6kwxkM/WgR//puHV0bV2fmg10tmRQDru2hiEhJ4C3gciAeWCwi01R1rVebaGAk0EVVD4iI\nnTcMA/M3JvPgxOUcTs+063CKGTcPedoDm1R1C4CITACuAtZ6tRkGvKWqBwBUda+L8RiXZWXn8O9Z\nG3lr7iYaVq/AZ3e0p+n5lYIdlilEbiaU2sBOr/vxQAefNo0BROQXoCTwjKp+77shERkODAeoW7eu\nK8Gas7PrYBr3j1/G4m0HGBgXxTP9WlCutHXRFTdu/o/728dVP88fDXQDooCfRaSlqqac8CDVMcAY\ngLi4ON9tmCCbvW4Pj365gmNZOfz7+jY21Wcx5mZCiQfqeN2PAhL9tFmgqpnAVhFZj5NgFrsYlwmQ\nY1k5vPT9H3wwfyvNa1bizSGxVmqgmHMzoSwGokWkAZAADAKG+LSZAgwGxolINZxDoC0uxmQCZOf+\nVO774ndWxB9kaKd6jOzTzCrPm/wTioh0Am4ELgZqAmnAauA74DNVPejvcaqaJSL3ATNx+kfGquoa\nERkFLFHVaZ51V4jIWiAbGKGq+wLwuoyLZq3dw0OTliPAuze2pVdLq/FqHKKad5eEiMzAOUyZCiwB\n9gKROHsS3YG+wKue5FAo4uLidMmSJYX1dMZLVnYOr/64gbfnbqZl7Uq8c0M76pxrU36GIxFZqqpx\np/u4/PZQblLVZJ9lR4DfPbdXPIcqJswlHc7g/vHL+G3LPga3r8vTfZvbIY45ySkTip9kckZtTNG2\nZNt+7v3id1JSMxk9IIbr4urk/yBTLJ0yoYjIVk481Ste91VVG7oVmAk+VWXsL9t4fvo6ap9Tlsn3\ntKd5LRuoZvKW3yGP7zFUCWAg8CiwzJWITEg4kpHF379ayXerdnF58/N4+brWNlOfyVd+hzz7AESk\nBHATMAJYDlzpfU2OCU2nU3He2/Z9Rxn2yRI27T3CY72bcmfXC+xaHFMg+R3yRAC3AQ8B84GrVHVz\nYQRmzk5BK877mr8xmXu/+B2AT27rwEXR1uduCi6/Q56tQBbwb2AH0FpEWh9fqarfuBibOQsFqTjv\nTVX5cP5Wnpu+jugaFRlzczvqVS1fWOGaMJFfQpmF0wnb2nPzpoAllBB1qorzvtIzs3l88iq++T2B\nXi3O55WBrSlfxi7sM6cvvz6UWwopDhNgeVWcr1Wl7An3dx9M585Pl7Ai/iAPX96Y+7o3okQJ6y8x\nZ+aUFdtE5DcRuTSPdbPdCckEgr+K82UjSjKiZ5Pc+0u3H6Dvm/PZtPcI793Ujvsvi7ZkYs5KfiUg\n6wJvisjLng5ab+e6FJMJgKtja/P8Na2oXaUsAtSuUpbnr2mV238yZVkCg8csoFzpkky+tws9W5wf\n3IBNWMjvQHkPcBHwH2ChiAxW1fWedVaXJMRdHVv7pA5YVeXfszby+uyNdLzgXN69sR1VypUOUoQm\n3ORbpFpVU1X1DuCfwI8icpdnle0bFzEZWdk8NHE5r8/eyLVto/jktg6WTExAFbgrX1Uni8ginNol\nfQCrpFOE7D96jDs/XcLibQd49IrG3Nu9kQ1WMwGXX0I5oWi0qiYAl4vICOAK16IyAbU56Qi3jVvM\nroPp/GdwLP1a1wp2SCZMnfKQR1V75bF8tKpGuhOSCaSFW/Zxzdu/cjg9i/HDOlgyMa7K77TxEyKS\n59kcEblURP4S+LBMIHy/ehc3jV1E1QqlmXJPF9rVsxNzxl35HfKsAv4rIuk4BZWScCq2RQOxwI/A\nc65GaM7IZwu289TU1bSuU4WxQy/knPLW+Wrcl99I2anAVM8Mf11wasoeAj4Dhquq//HdJmi8Twtf\n2rQGbw6JtflxTKEp0CdNVTcCG72XiUg9YLsbQZkzk52jPDFlNeMX7WBAuyiev6YVESVdnb7amBMU\ntOp9bWCequ4VkRjgMZwq+FYLMESkZ2bzwIRlzFyzh3u6NWREzyZ2WtgUuvw6ZUcDY4Frge9E5Gmc\nfpOFOP1RCRpHAAASO0lEQVQoJgQcTs/k5rGL+GHtHp7u25y/9WpqycQERX57KFcCsaqaLiLn4Eyp\nEeM5BDIh4MDRYwz9aBFrEw/x+iAbY2KCK7+Ekqaq6QCqekBE1lsyCR17D6dz0weL2LrvKO/d1I7L\nmp0X7JBMMZdfQmkoIt6TeNX3vq+q/dwJy+QnMSWNGz5YyJ5D6Xx0y4V0aWSlGk3w5ZdQrvK5/4pb\ngZiC277vKEPeX8ihtEw+vb29DVgzISO/hLIWqO5b4V5EWuBznY8pHJv2HmbI+wvJzM5h/PCOtKxd\nOdghGZMrv0EKbwDV/SyPAl4PfDjmVDbtPcKgMQvJUZh4ZydLJibk5JdQWqnqT74LVXUmEONOSMaf\nzUlHGPz+AgAmDO9A4/MqBjkiY06WX0I51VRxNo1cIdmSdITBYxagqowf1oFGNSyZmNCUX0LZ6Cmm\ndAIR6Q1scSck421r8lEGv7+A7Bzli2EdibY9ExPC8uuUfQj4VkQGAks9y+KAToCVLXDZ1uSjDBrz\nG1nZTjKxwxwT6vIrsLQBaAX8BNQH6gFzcaYnfcDl2Iq1+AOpDHl/AZnZyufDOtDkfEsmJvTle3Gg\nqmYAH4lILDAYeBpnitKvXY6t2Eo6nMGNHyzkaEYWE4Z3oun5lYIdkjEFkt9k6Y2BQTiJZB8wERBV\n7V4IsRVLB1MzuenDhew5lMFnd3SgeS1LJqboyK9T9g/gMqCvql6kqm8A2fk8JpeI9BKR9SKySUQe\nO0W7ASKiIhJX0G2Ho6MZWdw6bhFbko4y5uZ2tKt3TrBDMua05JdQrgV2A3NE5H0RuYwCzscjIiWB\nt4DeQHNgsIg099OuInA/TkmEYis9M5s7P13K8p0p/GdwLBdH+xtPaExoy69TdrKqXg80xemMfQg4\nT0TeEZH8ptFoD2xS1S2qegyYwMnXBoEzgdhLQPrpBh8usrJzuH/8MuZvSualAa3p1dKmBTVFU4Hq\nA6rqUVX9XFX/gjPsfjlO1bZTqQ3s9Lof71mWy9PRW0dVvz3VhkRkuIgsEZElSUlJBQm5yFBVnpy6\nOrc40oB2UcEOyZgzdtoFR1V1v6q+p6qX5tPU36FR7nzIIlICeA14pADPOUZV41Q1rnr18DoUePN/\nmxi/aCf3dm/IrV0aBDscY86Km+XQ4zmx5mwUTsW34yoCLYG5nnKF5wPTRKSfqi5xMS7XTVmWwOiZ\n60lMSaNWlbKM6NnkpEnLAb5cspNXftzANbG1efSKJkGI1JjAcjOhLAaiRaQBkIBz+nnI8ZWqehDI\nrQokInOBR8MhmYz8ZhVpmc7JsISUNEZ+swrghKQyb0MSI79ZRZdGVXnh2hirAWvCgmtzLKhqFnAf\nMBNYB0xS1TUiMkpEwrbS2+iZ63OTyXFpmdmMnrk+9/6axIPc/dlSGtWowDs3tqN0KZvqwoQHV2eA\nUtXpwHSfZU/l0babm7EUlsQU/3OfHV+ekJLGrR8tplLZCMbd2p5KkXbRtgkf9tMYYLWqlM1z+eH0\nTG79aBFpmdmMu7U951e2+eZNeLGEEmAjejahbETJE5aVjSjJI5c35sEJy9mcdJR3bmhnF/uZsGST\n3gbY8Y5X37M863YfYvYfexl1VQsuirYK9SY8WUJxwdWxtU84o/P10nje+2kLN3asy82d6gcvMGNc\nZoc8Llu6/QAjv1lFpwuq8nTfFsEOxxhXWUJxUUJKGnd+uoSaVSJ5+4a2RJS0t9uENzvkcUnasWyG\nfbyEjMwcJgyP45zypYMdkjGus4TiAlXl8cmrWLf7EGOHXmhV6k2xYfvgLvh0wXYmL0vg4R6N6d60\nRrDDMabQWEIJsKXb9zPqv2u5rGkN7u3eKNjhGFOoLKEEUNLhDO75/HdqVSnLqwPbUKKEXfBnihfr\nQwmQrOwc/jr+d1JSM5l8T3sql7NrdEzxYwklQF7+YQMLtuznletaW6V6U2zZIU8AzF2/l3d/2szg\n9nW41ko4mmLMEspZ2nMonUcmraDJeRVtJKwp9iyhnIXsHOXBCctJPZbNm0NiifS5ytiY4sb6UM7C\nW3M28duWfbw0IIZom8jcGNtDOVMLt+zj37M2cHWbWlxn/SbGAJZQzsjB1EwenLicelXL86/+razA\ntDEedshzBp6YupqkwxlMvqcLFcrYW2jMcbaHcpqmLk/gvysSebBHNK2iKgc7HGNCiiWU05CQksYT\nU1bTrt453HVJw2CHY0zIsYRSQDk5yqOTVpCTo7w2sA2lrFiSMSexb0UBjf1lK79t2cdTfZtTt2q5\nYIdjTEiyhFIAG/Yc5qXv13NF8/MYGFcn/wcYU0xZQslHVnYOI75cQYXIUjx3jZ0iNuZU7JxnPj6Y\nv5UV8Qd5Y3As1SqUCXY4xoQ020M5hU17j/Dqjxvo2eI8/hJTM9jhGBPyLKHkITtH+dtXKygbUZJ/\nXt3SDnWMKQBLKHkY9+s2ft+RwtN9m1Ojok1qbkxBWELxY1vyUUbP/INLm9agv9eUosaYU7OE4uP4\nnDoRJUrwnF34Z8xpsYTiY/KyBH7dvI+/9W7K+ZXtUMeY02EJxUtK6jH+77t1tKlThRva1w12OMYU\nOa4mFBHpJSLrRWSTiDzmZ/3DIrJWRFaKyGwRqedmPPl5YcYfpKRl8lz/VjanjjFnwLWEIiIlgbeA\n3kBzYLCINPdptgyIU9UY4CvgJbfiyc+SbfuZsHgnt3Wpb9NgGHOG3NxDaQ9sUtUtqnoMmABc5d1A\nVeeoaqrn7gIgKLUUM7Nz+Mfk1dSqHMmDPRoHIwRjwoKbCaU2sNPrfrxnWV5uB2b4WyEiw0VkiYgs\nSUpKCmCIjg9+3sr6PYd59qqWlLcKbMacMTcTir9OCPXbUORGIA4Y7W+9qo5R1ThVjatevXoAQ4TE\nlDRen72BK5qfx+XNzwvoto0pbtz8OY4HvK/1jwISfRuJSA/gH8AlqprhYjx+PTd9HarwVF/f7h1j\nzOlycw9lMRAtIg1EpDQwCJjm3UBEYoH3gH6qutfFWPxasGUf367cxd3dGhJ1jhVNMuZsuZZQVDUL\nuA+YCawDJqnqGhEZJSL9PM1GAxWAL0VkuYhMy2NzAZeVncMz09ZQu0pZ7uxq9WGNCQRXeyBVdTow\n3WfZU15/93Dz+U9l/OKd/LH7MG/f0JaypW0KUWMCoViOlD1w9Biv/LCeThdUpXfL84MdjjFho1gm\nlFd/3MDh9Cye7tfcLv4zJoCKXUL5Y/chPl+4nZs61qPp+TYi1phAKnYJ5bnpf1AxMoIHe0QHOxRj\nwk6xSig/b0xi3oYk/nppI6qUKx3scIwJO8UmoeTkKM9P/4Ooc8pyU6egXtRsTNgqNgll8rIE1u46\nxIieTShTyk4TG+OGYpFQ0jOzeeWH9bSqXZm+MbWCHY4xYatYJJSPftlG4sF0Hu/TzAonGeOisE8o\n+48e4+05m7isaQ06Nawa7HCMCWthWfxjyrIERs9cT2JKGuXKlCT1WDaP9W4a7LCMCXtht4cyZVkC\nI79ZRUJKGgoczcimBMKaxEPBDs2YsBd2CWX0zPWkZWafsCxbldEz1wcpImOKj7BLKIkpaae13BgT\nOGGXUGpVKXtay40xgRN2CcUZuHbiyyobUZIRPZsEKSJjio+wSyhXx9amWc1KuRWya1cpy/PXtOJq\nm/TcGNeF3WnjtYmHWL4zhfu6N+JR2ysxplCF3R7Kqz9uoGJkKYZdfEGwQzGm2AmrhLI64SCz1u1h\n2MUXULlcRLDDMabYCauE8ub/NlExshS3dKkf7FCMKZbCJqFs2HOY79fs5pbO9akUaXsnxgRD2CSU\nt+dsolzpktzWpUGwQzGm2AqLhLIt+SjTViRyU8d6nFPeSjsaEyxhkVDembuZiJIluP1i2zsxJpiK\nfEJJSEnj69/jGdy+LjUqRgY7HGOKtSKfUN77aTMiMLyrjTsxJtiKdELZeyidCYt3MqBdlF38Z0wI\nKNIJ5aNft5GVncNdlzQMdijGGIpwQjmSkcVnC7bTu2VN6lUtH+xwjDEU4YQycfFODqdnWd+JMSGk\nSCaUrOwcxs7fSvsG59K6TpVgh2OM8SiSCWX66t0kpKQx3K4oNiakFMmEMmbeZi6oXp5Lm9YIdijG\nGC9FLqEcychidcIhhl18gc0CaEyIcTWhiEgvEVkvIptE5DE/68uIyETP+oUiUj+/bSYfzqBahdL0\nt5KOxoQc1xKKiJQE3gJ6A82BwSLS3KfZ7cABVW0EvAa8mN92D2dkcXOn+kRGlAx0yMaYs+TmHkp7\nYJOqblHVY8AE4CqfNlcBH3v+/gq4TEROeRwjAjd2rBfwYI0xZ8/NItW1gZ1e9+OBDnm1UdUsETkI\nVAWSvRuJyHBguOduRtUKZVa7ErE7quHzekJYUYoVila8RSlWgDOq8O5mQvG3p6Fn0AZVHQOMARCR\nJaoad/bhFY6iFG9RihWKVrxFKVZw4j2Tx7l5yBMP1PG6HwUk5tVGREoBlYH9LsZkjHGRmwllMRAt\nIg1EpDQwCJjm02YaMNTz9wDgf6p60h6KMaZocO2Qx9Mnch8wEygJjFXVNSIyCliiqtOAD4FPRWQT\nzp7JoAJseoxbMbukKMVblGKFohVvUYoVzjBesR0CY0ygFLmRssaY0GUJxRgTMCGbUNwYtu+WAsT6\nsIisFZGVIjJbRII6Mi+/eL3aDRARFZGgne4sSKwiMtDz/q4RkS8KO0afWPL7LNQVkTkisszzeegT\njDg9sYwVkb0i4ndclzj+43ktK0Wkbb4bVdWQu+F04m4GLgBKAyuA5j5t7gHe9fw9CJgYwrF2B8p5\n/r47WLEWNF5Pu4rAPGABEBeqsQLRwDLgHM/9GqH83uJ0dt7t+bs5sC2I8XYF2gKr81jfB5iBM16s\nI7Awv22G6h6KK8P2XZJvrKo6R1VTPXcX4IzJCZaCvLcA/wReAtILMzgfBYl1GPCWqh4AUNW9hRyj\nt4LEq0Alz9+VOXlsVqFR1XmcetzXVcAn6lgAVBGRmqfaZqgmFH/D9n0vLz5h2D5wfNh+YStIrN5u\nx8n6wZJvvCISC9RR1W8LMzA/CvLeNgYai8gvIrJARHoVWnQnK0i8zwA3ikg8MB34a+GEdkZO97Pt\n6tD7sxGwYfuFoMBxiMiNQBxwiasRndop4xWREjhXft9SWAGdQkHe21I4hz3dcPb8fhaRlqqa4nJs\n/hQk3sHAOFV9RUQ64YzDaqmqOe6Hd9pO+zsWqnsoRWnYfkFiRUR6AP8A+qlqRiHF5k9+8VYEWgJz\nRWQbzrHztCB1zBb0czBVVTNVdSuwHifBBENB4r0dmASgqr8BkTgXDoaiAn22TxCsDqF8OotKAVuA\nBvzZudXCp829nNgpOymEY43F6ayLLgrvrU/7uQSvU7Yg720v4GPP39VwdtGrhnC8M4BbPH8383xB\nJYifh/rk3Sl7JSd2yi7Kd3vBeiEFeKF9gA2eL+I/PMtG4fzCg5PZvwQ2AYuAC0I41lnAHmC55zYt\nlN9bn7ZBSygFfG8FeBVYC6wCBoXye4tzZucXT7JZDlwRxFjHA7uATJy9kduBu4C7vN7btzyvZVVB\nPgc29N4YEzCh2odijCmCLKEYYwLGEooxJmAsoRhjAsYSijEmYCyhmAITkToislVEzvXcP8dzv56I\n1BSRbz3L/09ElnvdNohItohUEJFbRCRHRGK8trv6+NXiIjJLRM4JxuszZ88SiikwVd0JvAO84Fn0\nAjBGVbcDDwPve9r9Q1XbHL/h1Bd+XlWPeB4XjzNq2J9Pca4kN0WQJRRzul4DOorIg8BFwCue5dcC\n3/s29ly/1AjnorjjvgVaiIi/uV+m4VzvYoogSyjmtKhqJjACJ7E8qKrHRKQBzpSyJ1yj5DmMeQG4\nQZ0rwo/LwSmN8Lif7R8AyohIMK4cN2fJEoo5E71xhmy39NyvCSR5N/DMbf0Z8KSqbvKzjS9w9nQa\n+Fm3F6gVuHBNYbGEYk6LiLQBLse5WOwhT8GdNJxrq7w9AexS1Y/8bcezx/IK8Hc/qyM92zRFjCUU\nU2Ceinjv4Bzq7ABGAy/jXAxX36tdR5x6KsNP3soJxgE9gOo+z3E+sC1ggZtCYwnFnI5hwA5V/dFz\n/22gKU7RqM0i0siz/FmgHDDH5/RxQ++NqVMm8T9ADa/F7YAFPn0upoiwq41NQIhIf6Cdqj5xltt5\nHae8w+zARGYKU6iWgDRFjKpODtCZmdWWTIou20MxxgSM9aEYYwLGEooxJmAsoRhjAsYSijEmYCyh\nGGMC5v8BIXGhEi4OEEMAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x110d70e48>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "plt.plot(eq.X.sel(component='ZN', vertex=0).squeeze(), acr_zn, label='Calculated')\n",
    "# add experimental data\n",
    "plt.scatter(exp_x_zn, exp_acr_zn, label='Yazawa 1970')\n",
    "\n",
    "plt.xlim((0, 1))\n",
    "plt.ylim((0, 1))\n",
    "plt.gca().set_aspect(1)\n",
    "plt.xlabel('X(ZN)')\n",
    "plt.ylabel('ACR(ZN)')\n",
    "plt.title('Activity of Zn at 1073K')\n",
    "plt.legend(loc=0)"
   ]
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "kernelspec": {
   "display_name": "Python [conda env:espei]",
   "language": "python",
   "name": "conda-env-espei-py"
  },
  "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.6.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}