new_LOM_Analysis.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "17802b23-9957-497d-bb25-6d2f922ec8d2",
   "metadata": {},
   "source": [
    "# Standard imports"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "03173c59-9f9e-4f9d-a1ad-b4e54514f27d",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "12:36AM 06s INFO [__init__]: TransmonBuilder with system_type TRANSMON registered to QuantumSystemRegistry\n",
      "12:36AM 06s INFO [__init__]: FluxoniumBuilder with system_type FLUXONIUM registered to QuantumSystemRegistry\n",
      "12:36AM 06s INFO [__init__]: TLResonatorBuilder with system_type TL_RESONATOR registered to QuantumSystemRegistry\n",
      "12:36AM 06s INFO [__init__]: LumpedResonatorBuilder with system_type LUMPED_RESONATOR registered to QuantumSystemRegistry\n"
     ]
    }
   ],
   "source": [
    "%load_ext autoreload\n",
    "%autoreload 2\n",
    "\n",
    "\n",
    "import scqubits as scq\n",
    "\n",
    "from qiskit_metal.analyses.quantization.lumped_capacitive import load_q3d_capacitance_matrix\n",
    "from qiskit_metal.analyses.quantization.lom_core_analysis import CompositeSystem, Cell, Subsystem, QuantumSystemRegistry\n",
    "                        \n",
    "from scipy.constants import speed_of_light as c_light\n",
    "\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "bc734149-fd46-4772-92c7-3f057658d78f",
   "metadata": {},
   "source": [
    "# Quantum System Reegistry"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "497d7464-7d12-4d51-aa8c-e9d0987e5a8e",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'TRANSMON': qiskit_metal.analyses.quantization.lom_core_analysis.TransmonBuilder,\n",
       " 'FLUXONIUM': qiskit_metal.analyses.quantization.lom_core_analysis.FluxoniumBuilder,\n",
       " 'TL_RESONATOR': qiskit_metal.analyses.quantization.lom_core_analysis.TLResonatorBuilder,\n",
       " 'LUMPED_RESONATOR': qiskit_metal.analyses.quantization.lom_core_analysis.LumpedResonatorBuilder}"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "QuantumSystemRegistry.registry()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c1497969-e4ad-4c26-8519-d4cf7801fc23",
   "metadata": {},
   "source": [
    "## Load the capacitance matrices for each qubit"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "6e99011d-55d4-4339-806f-a1a5c7e795c6",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Imported capacitance matrix with UNITS: [fF] now converted to USER UNITS:[fF]                 from file:\n",
      "\t./capacitance_matrices/Q1_cap_matrix.txt\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<style type=\"text/css\">\n",
       "#T_cae0f_row0_col0, #T_cae0f_row1_col1, #T_cae0f_row2_col2, #T_cae0f_row3_col3, #T_cae0f_row4_col4, #T_cae0f_row5_col5, #T_cae0f_row6_col6 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 100.0%, transparent 100.0%);\n",
       "}\n",
       "#T_cae0f_row0_col1 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 39.9%, transparent 39.9%);\n",
       "}\n",
       "#T_cae0f_row0_col2 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 39.7%, transparent 39.7%);\n",
       "}\n",
       "#T_cae0f_row0_col3, #T_cae0f_row1_col3 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 49.2%, transparent 49.2%);\n",
       "}\n",
       "#T_cae0f_row0_col4, #T_cae0f_row1_col4, #T_cae0f_row2_col4 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 1.3%, transparent 1.3%);\n",
       "}\n",
       "#T_cae0f_row0_col5 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 16.0%, transparent 16.0%);\n",
       "}\n",
       "#T_cae0f_row0_col6, #T_cae0f_row2_col6 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 29.8%, transparent 29.8%);\n",
       "}\n",
       "#T_cae0f_row1_col0 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 39.8%, transparent 39.8%);\n",
       "}\n",
       "#T_cae0f_row1_col2 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 40.0%, transparent 40.0%);\n",
       "}\n",
       "#T_cae0f_row1_col5 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 25.6%, transparent 25.6%);\n",
       "}\n",
       "#T_cae0f_row1_col6 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 19.8%, transparent 19.8%);\n",
       "}\n",
       "#T_cae0f_row2_col0 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 39.6%, transparent 39.6%);\n",
       "}\n",
       "#T_cae0f_row2_col1, #T_cae0f_row3_col0 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 40.1%, transparent 40.1%);\n",
       "}\n",
       "#T_cae0f_row2_col3 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 49.1%, transparent 49.1%);\n",
       "}\n",
       "#T_cae0f_row2_col5 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 15.9%, transparent 15.9%);\n",
       "}\n",
       "#T_cae0f_row3_col1, #T_cae0f_row3_col2 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 40.2%, transparent 40.2%);\n",
       "}\n",
       "#T_cae0f_row3_col4 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 7.8%, transparent 7.8%);\n",
       "}\n",
       "#T_cae0f_row3_col5 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 26.8%, transparent 26.8%);\n",
       "}\n",
       "#T_cae0f_row3_col6 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 30.9%, transparent 30.9%);\n",
       "}\n",
       "#T_cae0f_row4_col0, #T_cae0f_row4_col1, #T_cae0f_row4_col2, #T_cae0f_row4_col3, #T_cae0f_row4_col5, #T_cae0f_row4_col6, #T_cae0f_row6_col4 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "}\n",
       "#T_cae0f_row5_col0 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 24.2%, transparent 24.2%);\n",
       "}\n",
       "#T_cae0f_row5_col1 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 38.5%, transparent 38.5%);\n",
       "}\n",
       "#T_cae0f_row5_col2 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 24.3%, transparent 24.3%);\n",
       "}\n",
       "#T_cae0f_row5_col3 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 48.8%, transparent 48.8%);\n",
       "}\n",
       "#T_cae0f_row5_col4 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 1.6%, transparent 1.6%);\n",
       "}\n",
       "#T_cae0f_row5_col6 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 9.5%, transparent 9.5%);\n",
       "}\n",
       "#T_cae0f_row6_col0 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 38.4%, transparent 38.4%);\n",
       "}\n",
       "#T_cae0f_row6_col1 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 24.1%, transparent 24.1%);\n",
       "}\n",
       "#T_cae0f_row6_col2 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 38.3%, transparent 38.3%);\n",
       "}\n",
       "#T_cae0f_row6_col3 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 48.6%, transparent 48.6%);\n",
       "}\n",
       "#T_cae0f_row6_col5 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 5.7%, transparent 5.7%);\n",
       "}\n",
       "</style>\n",
       "<table id=\"T_cae0f_\">\n",
       "  <thead>\n",
       "    <tr>\n",
       "      <th class=\"blank level0\" >&nbsp;</th>\n",
       "      <th class=\"col_heading level0 col0\" >a_connector_pad_Q1</th>\n",
       "      <th class=\"col_heading level0 col1\" >b_connector_pad_Q1</th>\n",
       "      <th class=\"col_heading level0 col2\" >c_connector_pad_Q1</th>\n",
       "      <th class=\"col_heading level0 col3\" >cl_metal_Q1</th>\n",
       "      <th class=\"col_heading level0 col4\" >ground_main_plane</th>\n",
       "      <th class=\"col_heading level0 col5\" >pad_bot_Q1</th>\n",
       "      <th class=\"col_heading level0 col6\" >pad_top_Q1</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th id=\"T_cae0f_level0_row0\" class=\"row_heading level0 row0\" >a_connector_pad_Q1</th>\n",
       "      <td id=\"T_cae0f_row0_col0\" class=\"data row0 col0\" >50.69</td>\n",
       "      <td id=\"T_cae0f_row0_col1\" class=\"data row0 col1\" >-0.34</td>\n",
       "      <td id=\"T_cae0f_row0_col2\" class=\"data row0 col2\" >-0.45</td>\n",
       "      <td id=\"T_cae0f_row0_col3\" class=\"data row0 col3\" >-0.01</td>\n",
       "      <td id=\"T_cae0f_row0_col4\" class=\"data row0 col4\" >-34.01</td>\n",
       "      <td id=\"T_cae0f_row0_col5\" class=\"data row0 col5\" >-13.50</td>\n",
       "      <td id=\"T_cae0f_row0_col6\" class=\"data row0 col6\" >-1.52</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_cae0f_level0_row1\" class=\"row_heading level0 row1\" >b_connector_pad_Q1</th>\n",
       "      <td id=\"T_cae0f_row1_col0\" class=\"data row1 col0\" >-0.34</td>\n",
       "      <td id=\"T_cae0f_row1_col1\" class=\"data row1 col1\" >50.75</td>\n",
       "      <td id=\"T_cae0f_row1_col2\" class=\"data row1 col2\" >-0.14</td>\n",
       "      <td id=\"T_cae0f_row1_col3\" class=\"data row1 col3\" >-0.01</td>\n",
       "      <td id=\"T_cae0f_row1_col4\" class=\"data row1 col4\" >-34.19</td>\n",
       "      <td id=\"T_cae0f_row1_col5\" class=\"data row1 col5\" >-1.50</td>\n",
       "      <td id=\"T_cae0f_row1_col6\" class=\"data row1 col6\" >-13.69</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_cae0f_level0_row2\" class=\"row_heading level0 row2\" >c_connector_pad_Q1</th>\n",
       "      <td id=\"T_cae0f_row2_col0\" class=\"data row2 col0\" >-0.45</td>\n",
       "      <td id=\"T_cae0f_row2_col1\" class=\"data row2 col1\" >-0.14</td>\n",
       "      <td id=\"T_cae0f_row2_col2\" class=\"data row2 col2\" >50.73</td>\n",
       "      <td id=\"T_cae0f_row2_col3\" class=\"data row2 col3\" >-0.05</td>\n",
       "      <td id=\"T_cae0f_row2_col4\" class=\"data row2 col4\" >-34.13</td>\n",
       "      <td id=\"T_cae0f_row2_col5\" class=\"data row2 col5\" >-13.51</td>\n",
       "      <td id=\"T_cae0f_row2_col6\" class=\"data row2 col6\" >-1.59</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_cae0f_level0_row3\" class=\"row_heading level0 row3\" >cl_metal_Q1</th>\n",
       "      <td id=\"T_cae0f_row3_col0\" class=\"data row3 col0\" >-0.01</td>\n",
       "      <td id=\"T_cae0f_row3_col1\" class=\"data row3 col1\" >-0.01</td>\n",
       "      <td id=\"T_cae0f_row3_col2\" class=\"data row3 col2\" >-0.05</td>\n",
       "      <td id=\"T_cae0f_row3_col3\" class=\"data row3 col3\" >16.01</td>\n",
       "      <td id=\"T_cae0f_row3_col4\" class=\"data row3 col4\" >-15.53</td>\n",
       "      <td id=\"T_cae0f_row3_col5\" class=\"data row3 col5\" >-0.12</td>\n",
       "      <td id=\"T_cae0f_row3_col6\" class=\"data row3 col6\" >-0.21</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_cae0f_level0_row4\" class=\"row_heading level0 row4\" >ground_main_plane</th>\n",
       "      <td id=\"T_cae0f_row4_col0\" class=\"data row4 col0\" >-34.01</td>\n",
       "      <td id=\"T_cae0f_row4_col1\" class=\"data row4 col1\" >-34.19</td>\n",
       "      <td id=\"T_cae0f_row4_col2\" class=\"data row4 col2\" >-34.13</td>\n",
       "      <td id=\"T_cae0f_row4_col3\" class=\"data row4 col3\" >-15.53</td>\n",
       "      <td id=\"T_cae0f_row4_col4\" class=\"data row4 col4\" >248.10</td>\n",
       "      <td id=\"T_cae0f_row4_col5\" class=\"data row4 col5\" >-33.24</td>\n",
       "      <td id=\"T_cae0f_row4_col6\" class=\"data row4 col6\" >-37.77</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_cae0f_level0_row5\" class=\"row_heading level0 row5\" >pad_bot_Q1</th>\n",
       "      <td id=\"T_cae0f_row5_col0\" class=\"data row5 col0\" >-13.50</td>\n",
       "      <td id=\"T_cae0f_row5_col1\" class=\"data row5 col1\" >-1.50</td>\n",
       "      <td id=\"T_cae0f_row5_col2\" class=\"data row5 col2\" >-13.51</td>\n",
       "      <td id=\"T_cae0f_row5_col3\" class=\"data row5 col3\" >-0.12</td>\n",
       "      <td id=\"T_cae0f_row5_col4\" class=\"data row5 col4\" >-33.24</td>\n",
       "      <td id=\"T_cae0f_row5_col5\" class=\"data row5 col5\" >90.53</td>\n",
       "      <td id=\"T_cae0f_row5_col6\" class=\"data row5 col6\" >-26.22</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_cae0f_level0_row6\" class=\"row_heading level0 row6\" >pad_top_Q1</th>\n",
       "      <td id=\"T_cae0f_row6_col0\" class=\"data row6 col0\" >-1.52</td>\n",
       "      <td id=\"T_cae0f_row6_col1\" class=\"data row6 col1\" >-13.69</td>\n",
       "      <td id=\"T_cae0f_row6_col2\" class=\"data row6 col2\" >-1.59</td>\n",
       "      <td id=\"T_cae0f_row6_col3\" class=\"data row6 col3\" >-0.21</td>\n",
       "      <td id=\"T_cae0f_row6_col4\" class=\"data row6 col4\" >-37.77</td>\n",
       "      <td id=\"T_cae0f_row6_col5\" class=\"data row6 col5\" >-26.22</td>\n",
       "      <td id=\"T_cae0f_row6_col6\" class=\"data row6 col6\" >83.74</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n"
      ],
      "text/plain": [
       "<pandas.io.formats.style.Styler at 0x7fc77a7d77c0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Q1\n",
    "path1 = \"./capacitance_matrices/Q1_cap_matrix.txt\"\n",
    "t1_mat, _, _, _ = load_q3d_capacitance_matrix(path1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "24275e48-f418-4763-9e7a-1539e3909b14",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Imported capacitance matrix with UNITS: [fF] now converted to USER UNITS:[fF]                 from file:\n",
      "\t./capacitance_matrices/Q2_cap_matrix.txt\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<style type=\"text/css\">\n",
       "#T_472d2_row0_col0, #T_472d2_row1_col1, #T_472d2_row2_col2, #T_472d2_row3_col3, #T_472d2_row4_col4, #T_472d2_row5_col5, #T_472d2_row6_col6, #T_472d2_row7_col7 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 100.0%, transparent 100.0%);\n",
       "}\n",
       "#T_472d2_row0_col1, #T_472d2_row1_col0, #T_472d2_row4_col2 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 39.6%, transparent 39.6%);\n",
       "}\n",
       "#T_472d2_row0_col2, #T_472d2_row2_col0, #T_472d2_row4_col1 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 39.4%, transparent 39.4%);\n",
       "}\n",
       "#T_472d2_row0_col3, #T_472d2_row1_col3 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 49.1%, transparent 49.1%);\n",
       "}\n",
       "#T_472d2_row0_col4, #T_472d2_row3_col4 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 39.8%, transparent 39.8%);\n",
       "}\n",
       "#T_472d2_row0_col5, #T_472d2_row5_col0, #T_472d2_row5_col1, #T_472d2_row5_col2, #T_472d2_row5_col3, #T_472d2_row5_col4, #T_472d2_row5_col6, #T_472d2_row5_col7 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "}\n",
       "#T_472d2_row0_col6, #T_472d2_row2_col6 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 14.3%, transparent 14.3%);\n",
       "}\n",
       "#T_472d2_row0_col7, #T_472d2_row2_col7 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 24.5%, transparent 24.5%);\n",
       "}\n",
       "#T_472d2_row1_col2, #T_472d2_row2_col1, #T_472d2_row4_col0 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 39.9%, transparent 39.9%);\n",
       "}\n",
       "#T_472d2_row1_col4 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 39.3%, transparent 39.3%);\n",
       "}\n",
       "#T_472d2_row1_col5, #T_472d2_row2_col5, #T_472d2_row4_col5 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 0.0%, transparent 0.0%);\n",
       "}\n",
       "#T_472d2_row1_col6, #T_472d2_row4_col6 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 24.6%, transparent 24.6%);\n",
       "}\n",
       "#T_472d2_row1_col7, #T_472d2_row4_col7 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 14.2%, transparent 14.2%);\n",
       "}\n",
       "#T_472d2_row2_col3 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 48.9%, transparent 48.9%);\n",
       "}\n",
       "#T_472d2_row2_col4 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 39.5%, transparent 39.5%);\n",
       "}\n",
       "#T_472d2_row3_col0 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 40.1%, transparent 40.1%);\n",
       "}\n",
       "#T_472d2_row3_col1, #T_472d2_row3_col2 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 40.0%, transparent 40.0%);\n",
       "}\n",
       "#T_472d2_row3_col5 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 6.2%, transparent 6.2%);\n",
       "}\n",
       "#T_472d2_row3_col6 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 25.8%, transparent 25.8%);\n",
       "}\n",
       "#T_472d2_row3_col7 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 25.7%, transparent 25.7%);\n",
       "}\n",
       "#T_472d2_row4_col3, #T_472d2_row7_col3 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 48.6%, transparent 48.6%);\n",
       "}\n",
       "#T_472d2_row6_col0, #T_472d2_row6_col2, #T_472d2_row7_col1 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 24.3%, transparent 24.3%);\n",
       "}\n",
       "#T_472d2_row6_col1, #T_472d2_row7_col0, #T_472d2_row7_col2 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 38.3%, transparent 38.3%);\n",
       "}\n",
       "#T_472d2_row6_col3 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 48.8%, transparent 48.8%);\n",
       "}\n",
       "#T_472d2_row6_col4 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 38.2%, transparent 38.2%);\n",
       "}\n",
       "#T_472d2_row6_col5, #T_472d2_row7_col5 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 1.4%, transparent 1.4%);\n",
       "}\n",
       "#T_472d2_row6_col7, #T_472d2_row7_col6 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 5.4%, transparent 5.4%);\n",
       "}\n",
       "#T_472d2_row7_col4 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 24.2%, transparent 24.2%);\n",
       "}\n",
       "</style>\n",
       "<table id=\"T_472d2_\">\n",
       "  <thead>\n",
       "    <tr>\n",
       "      <th class=\"blank level0\" >&nbsp;</th>\n",
       "      <th class=\"col_heading level0 col0\" >a_connector_pad_Q2</th>\n",
       "      <th class=\"col_heading level0 col1\" >b_connector_pad_Q2</th>\n",
       "      <th class=\"col_heading level0 col2\" >c_connector_pad_Q2</th>\n",
       "      <th class=\"col_heading level0 col3\" >cl_metal_Q2</th>\n",
       "      <th class=\"col_heading level0 col4\" >d_connector_pad_Q2</th>\n",
       "      <th class=\"col_heading level0 col5\" >ground_main_plane</th>\n",
       "      <th class=\"col_heading level0 col6\" >pad_bot_Q2</th>\n",
       "      <th class=\"col_heading level0 col7\" >pad_top_Q2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th id=\"T_472d2_level0_row0\" class=\"row_heading level0 row0\" >a_connector_pad_Q2</th>\n",
       "      <td id=\"T_472d2_row0_col0\" class=\"data row0 col0\" >51.41</td>\n",
       "      <td id=\"T_472d2_row0_col1\" class=\"data row0 col1\" >-0.38</td>\n",
       "      <td id=\"T_472d2_row0_col2\" class=\"data row0 col2\" >-0.59</td>\n",
       "      <td id=\"T_472d2_row0_col3\" class=\"data row0 col3\" >-0.01</td>\n",
       "      <td id=\"T_472d2_row0_col4\" class=\"data row0 col4\" >-0.16</td>\n",
       "      <td id=\"T_472d2_row0_col5\" class=\"data row0 col5\" >-34.37</td>\n",
       "      <td id=\"T_472d2_row0_col6\" class=\"data row0 col6\" >-13.50</td>\n",
       "      <td id=\"T_472d2_row0_col7\" class=\"data row0 col7\" >-1.50</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_472d2_level0_row1\" class=\"row_heading level0 row1\" >b_connector_pad_Q2</th>\n",
       "      <td id=\"T_472d2_row1_col0\" class=\"data row1 col0\" >-0.38</td>\n",
       "      <td id=\"T_472d2_row1_col1\" class=\"data row1 col1\" >51.42</td>\n",
       "      <td id=\"T_472d2_row1_col2\" class=\"data row1 col2\" >-0.16</td>\n",
       "      <td id=\"T_472d2_row1_col3\" class=\"data row1 col3\" >-0.01</td>\n",
       "      <td id=\"T_472d2_row1_col4\" class=\"data row1 col4\" >-0.58</td>\n",
       "      <td id=\"T_472d2_row1_col5\" class=\"data row1 col5\" >-34.36</td>\n",
       "      <td id=\"T_472d2_row1_col6\" class=\"data row1 col6\" >-1.50</td>\n",
       "      <td id=\"T_472d2_row1_col7\" class=\"data row1 col7\" >-13.50</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_472d2_level0_row2\" class=\"row_heading level0 row2\" >c_connector_pad_Q2</th>\n",
       "      <td id=\"T_472d2_row2_col0\" class=\"data row2 col0\" >-0.59</td>\n",
       "      <td id=\"T_472d2_row2_col1\" class=\"data row2 col1\" >-0.16</td>\n",
       "      <td id=\"T_472d2_row2_col2\" class=\"data row2 col2\" >51.38</td>\n",
       "      <td id=\"T_472d2_row2_col3\" class=\"data row2 col3\" >-0.05</td>\n",
       "      <td id=\"T_472d2_row2_col4\" class=\"data row2 col4\" >-0.38</td>\n",
       "      <td id=\"T_472d2_row2_col5\" class=\"data row2 col5\" >-34.31</td>\n",
       "      <td id=\"T_472d2_row2_col6\" class=\"data row2 col6\" >-13.49</td>\n",
       "      <td id=\"T_472d2_row2_col7\" class=\"data row2 col7\" >-1.50</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_472d2_level0_row3\" class=\"row_heading level0 row3\" >cl_metal_Q2</th>\n",
       "      <td id=\"T_472d2_row3_col0\" class=\"data row3 col0\" >-0.01</td>\n",
       "      <td id=\"T_472d2_row3_col1\" class=\"data row3 col1\" >-0.01</td>\n",
       "      <td id=\"T_472d2_row3_col2\" class=\"data row3 col2\" >-0.05</td>\n",
       "      <td id=\"T_472d2_row3_col3\" class=\"data row3 col3\" >15.94</td>\n",
       "      <td id=\"T_472d2_row3_col4\" class=\"data row3 col4\" >-0.15</td>\n",
       "      <td id=\"T_472d2_row3_col5\" class=\"data row3 col5\" >-15.38</td>\n",
       "      <td id=\"T_472d2_row3_col6\" class=\"data row3 col6\" >-0.11</td>\n",
       "      <td id=\"T_472d2_row3_col7\" class=\"data row3 col7\" >-0.15</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_472d2_level0_row4\" class=\"row_heading level0 row4\" >d_connector_pad_Q2</th>\n",
       "      <td id=\"T_472d2_row4_col0\" class=\"data row4 col0\" >-0.16</td>\n",
       "      <td id=\"T_472d2_row4_col1\" class=\"data row4 col1\" >-0.58</td>\n",
       "      <td id=\"T_472d2_row4_col2\" class=\"data row4 col2\" >-0.38</td>\n",
       "      <td id=\"T_472d2_row4_col3\" class=\"data row4 col3\" >-0.15</td>\n",
       "      <td id=\"T_472d2_row4_col4\" class=\"data row4 col4\" >51.48</td>\n",
       "      <td id=\"T_472d2_row4_col5\" class=\"data row4 col5\" >-34.29</td>\n",
       "      <td id=\"T_472d2_row4_col6\" class=\"data row4 col6\" >-1.50</td>\n",
       "      <td id=\"T_472d2_row4_col7\" class=\"data row4 col7\" >-13.50</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_472d2_level0_row5\" class=\"row_heading level0 row5\" >ground_main_plane</th>\n",
       "      <td id=\"T_472d2_row5_col0\" class=\"data row5 col0\" >-34.37</td>\n",
       "      <td id=\"T_472d2_row5_col1\" class=\"data row5 col1\" >-34.36</td>\n",
       "      <td id=\"T_472d2_row5_col2\" class=\"data row5 col2\" >-34.31</td>\n",
       "      <td id=\"T_472d2_row5_col3\" class=\"data row5 col3\" >-15.38</td>\n",
       "      <td id=\"T_472d2_row5_col4\" class=\"data row5 col4\" >-34.29</td>\n",
       "      <td id=\"T_472d2_row5_col5\" class=\"data row5 col5\" >271.68</td>\n",
       "      <td id=\"T_472d2_row5_col6\" class=\"data row5 col6\" >-30.08</td>\n",
       "      <td id=\"T_472d2_row5_col7\" class=\"data row5 col7\" >-30.03</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_472d2_level0_row6\" class=\"row_heading level0 row6\" >pad_bot_Q2</th>\n",
       "      <td id=\"T_472d2_row6_col0\" class=\"data row6 col0\" >-13.50</td>\n",
       "      <td id=\"T_472d2_row6_col1\" class=\"data row6 col1\" >-1.50</td>\n",
       "      <td id=\"T_472d2_row6_col2\" class=\"data row6 col2\" >-13.49</td>\n",
       "      <td id=\"T_472d2_row6_col3\" class=\"data row6 col3\" >-0.11</td>\n",
       "      <td id=\"T_472d2_row6_col4\" class=\"data row6 col4\" >-1.50</td>\n",
       "      <td id=\"T_472d2_row6_col5\" class=\"data row6 col5\" >-30.08</td>\n",
       "      <td id=\"T_472d2_row6_col6\" class=\"data row6 col6\" >86.24</td>\n",
       "      <td id=\"T_472d2_row6_col7\" class=\"data row6 col7\" >-23.77</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_472d2_level0_row7\" class=\"row_heading level0 row7\" >pad_top_Q2</th>\n",
       "      <td id=\"T_472d2_row7_col0\" class=\"data row7 col0\" >-1.50</td>\n",
       "      <td id=\"T_472d2_row7_col1\" class=\"data row7 col1\" >-13.50</td>\n",
       "      <td id=\"T_472d2_row7_col2\" class=\"data row7 col2\" >-1.50</td>\n",
       "      <td id=\"T_472d2_row7_col3\" class=\"data row7 col3\" >-0.15</td>\n",
       "      <td id=\"T_472d2_row7_col4\" class=\"data row7 col4\" >-13.50</td>\n",
       "      <td id=\"T_472d2_row7_col5\" class=\"data row7 col5\" >-30.03</td>\n",
       "      <td id=\"T_472d2_row7_col6\" class=\"data row7 col6\" >-23.77</td>\n",
       "      <td id=\"T_472d2_row7_col7\" class=\"data row7 col7\" >86.26</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n"
      ],
      "text/plain": [
       "<pandas.io.formats.style.Styler at 0x7fc77aae0cd0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Q2\n",
    "path2 = \"./capacitance_matrices/Q2_cap_matrix.txt\"\n",
    "t2_mat, _, _, _ = load_q3d_capacitance_matrix(path2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "88f7a91f-f955-4e77-954a-70d99a9b969c",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Imported capacitance matrix with UNITS: [fF] now converted to USER UNITS:[fF]                 from file:\n",
      "\t./capacitance_matrices/Q3_cap_matrix.txt\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<style type=\"text/css\">\n",
       "#T_9af88_row0_col0, #T_9af88_row1_col1, #T_9af88_row2_col2, #T_9af88_row3_col3, #T_9af88_row4_col4, #T_9af88_row5_col5, #T_9af88_row6_col6, #T_9af88_row7_col7 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 100.0%, transparent 100.0%);\n",
       "}\n",
       "#T_9af88_row0_col1, #T_9af88_row1_col0, #T_9af88_row4_col2 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 39.6%, transparent 39.6%);\n",
       "}\n",
       "#T_9af88_row0_col2, #T_9af88_row2_col0, #T_9af88_row4_col1 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 39.2%, transparent 39.2%);\n",
       "}\n",
       "#T_9af88_row0_col3, #T_9af88_row1_col3 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 49.1%, transparent 49.1%);\n",
       "}\n",
       "#T_9af88_row0_col4, #T_9af88_row1_col2, #T_9af88_row3_col4 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 39.8%, transparent 39.8%);\n",
       "}\n",
       "#T_9af88_row0_col5, #T_9af88_row2_col5, #T_9af88_row4_col5 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 0.0%, transparent 0.0%);\n",
       "}\n",
       "#T_9af88_row0_col6, #T_9af88_row1_col7, #T_9af88_row2_col6, #T_9af88_row4_col7 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 13.1%, transparent 13.1%);\n",
       "}\n",
       "#T_9af88_row0_col7, #T_9af88_row1_col6, #T_9af88_row2_col7, #T_9af88_row4_col6 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 24.0%, transparent 24.0%);\n",
       "}\n",
       "#T_9af88_row1_col4 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 39.1%, transparent 39.1%);\n",
       "}\n",
       "#T_9af88_row1_col5, #T_9af88_row5_col0, #T_9af88_row5_col1, #T_9af88_row5_col2, #T_9af88_row5_col3, #T_9af88_row5_col4, #T_9af88_row5_col6, #T_9af88_row5_col7 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "}\n",
       "#T_9af88_row2_col1, #T_9af88_row4_col0 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 39.9%, transparent 39.9%);\n",
       "}\n",
       "#T_9af88_row2_col3 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 49.0%, transparent 49.0%);\n",
       "}\n",
       "#T_9af88_row2_col4 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 39.5%, transparent 39.5%);\n",
       "}\n",
       "#T_9af88_row3_col0, #T_9af88_row3_col1 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 40.1%, transparent 40.1%);\n",
       "}\n",
       "#T_9af88_row3_col2 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 40.0%, transparent 40.0%);\n",
       "}\n",
       "#T_9af88_row3_col5 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 6.4%, transparent 6.4%);\n",
       "}\n",
       "#T_9af88_row3_col6 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 25.3%, transparent 25.3%);\n",
       "}\n",
       "#T_9af88_row3_col7 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 25.2%, transparent 25.2%);\n",
       "}\n",
       "#T_9af88_row4_col3 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 48.6%, transparent 48.6%);\n",
       "}\n",
       "#T_9af88_row6_col0, #T_9af88_row6_col2, #T_9af88_row7_col1, #T_9af88_row7_col4 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 24.6%, transparent 24.6%);\n",
       "}\n",
       "#T_9af88_row6_col1, #T_9af88_row6_col4, #T_9af88_row7_col0, #T_9af88_row7_col2 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 38.3%, transparent 38.3%);\n",
       "}\n",
       "#T_9af88_row6_col3 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 48.8%, transparent 48.8%);\n",
       "}\n",
       "#T_9af88_row6_col5, #T_9af88_row7_col5 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 2.3%, transparent 2.3%);\n",
       "}\n",
       "#T_9af88_row6_col7 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 5.3%, transparent 5.3%);\n",
       "}\n",
       "#T_9af88_row7_col3 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 48.7%, transparent 48.7%);\n",
       "}\n",
       "#T_9af88_row7_col6 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 5.4%, transparent 5.4%);\n",
       "}\n",
       "</style>\n",
       "<table id=\"T_9af88_\">\n",
       "  <thead>\n",
       "    <tr>\n",
       "      <th class=\"blank level0\" >&nbsp;</th>\n",
       "      <th class=\"col_heading level0 col0\" >a_connector_pad_Q3</th>\n",
       "      <th class=\"col_heading level0 col1\" >b_connector_pad_Q3</th>\n",
       "      <th class=\"col_heading level0 col2\" >c_connector_pad_Q3</th>\n",
       "      <th class=\"col_heading level0 col3\" >cl_metal_Q3</th>\n",
       "      <th class=\"col_heading level0 col4\" >d_connector_pad_Q3</th>\n",
       "      <th class=\"col_heading level0 col5\" >ground_main_plane</th>\n",
       "      <th class=\"col_heading level0 col6\" >pad_bot_Q3</th>\n",
       "      <th class=\"col_heading level0 col7\" >pad_top_Q3</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th id=\"T_9af88_level0_row0\" class=\"row_heading level0 row0\" >a_connector_pad_Q3</th>\n",
       "      <td id=\"T_9af88_row0_col0\" class=\"data row0 col0\" >52.29</td>\n",
       "      <td id=\"T_9af88_row0_col1\" class=\"data row0 col1\" >-0.42</td>\n",
       "      <td id=\"T_9af88_row0_col2\" class=\"data row0 col2\" >-0.76</td>\n",
       "      <td id=\"T_9af88_row0_col3\" class=\"data row0 col3\" >-0.01</td>\n",
       "      <td id=\"T_9af88_row0_col4\" class=\"data row0 col4\" >-0.18</td>\n",
       "      <td id=\"T_9af88_row0_col5\" class=\"data row0 col5\" >-34.97</td>\n",
       "      <td id=\"T_9af88_row0_col6\" class=\"data row0 col6\" >-13.47</td>\n",
       "      <td id=\"T_9af88_row0_col7\" class=\"data row0 col7\" >-1.52</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_9af88_level0_row1\" class=\"row_heading level0 row1\" >b_connector_pad_Q3</th>\n",
       "      <td id=\"T_9af88_row1_col0\" class=\"data row1 col0\" >-0.42</td>\n",
       "      <td id=\"T_9af88_row1_col1\" class=\"data row1 col1\" >52.29</td>\n",
       "      <td id=\"T_9af88_row1_col2\" class=\"data row1 col2\" >-0.18</td>\n",
       "      <td id=\"T_9af88_row1_col3\" class=\"data row1 col3\" >-0.01</td>\n",
       "      <td id=\"T_9af88_row1_col4\" class=\"data row1 col4\" >-0.76</td>\n",
       "      <td id=\"T_9af88_row1_col5\" class=\"data row1 col5\" >-34.98</td>\n",
       "      <td id=\"T_9af88_row1_col6\" class=\"data row1 col6\" >-1.52</td>\n",
       "      <td id=\"T_9af88_row1_col7\" class=\"data row1 col7\" >-13.48</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_9af88_level0_row2\" class=\"row_heading level0 row2\" >c_connector_pad_Q3</th>\n",
       "      <td id=\"T_9af88_row2_col0\" class=\"data row2 col0\" >-0.76</td>\n",
       "      <td id=\"T_9af88_row2_col1\" class=\"data row2 col1\" >-0.18</td>\n",
       "      <td id=\"T_9af88_row2_col2\" class=\"data row2 col2\" >52.31</td>\n",
       "      <td id=\"T_9af88_row2_col3\" class=\"data row2 col3\" >-0.06</td>\n",
       "      <td id=\"T_9af88_row2_col4\" class=\"data row2 col4\" >-0.42</td>\n",
       "      <td id=\"T_9af88_row2_col5\" class=\"data row2 col5\" >-34.94</td>\n",
       "      <td id=\"T_9af88_row2_col6\" class=\"data row2 col6\" >-13.48</td>\n",
       "      <td id=\"T_9af88_row2_col7\" class=\"data row2 col7\" >-1.52</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_9af88_level0_row3\" class=\"row_heading level0 row3\" >cl_metal_Q3</th>\n",
       "      <td id=\"T_9af88_row3_col0\" class=\"data row3 col0\" >-0.01</td>\n",
       "      <td id=\"T_9af88_row3_col1\" class=\"data row3 col1\" >-0.01</td>\n",
       "      <td id=\"T_9af88_row3_col2\" class=\"data row3 col2\" >-0.06</td>\n",
       "      <td id=\"T_9af88_row3_col3\" class=\"data row3 col3\" >15.96</td>\n",
       "      <td id=\"T_9af88_row3_col4\" class=\"data row3 col4\" >-0.16</td>\n",
       "      <td id=\"T_9af88_row3_col5\" class=\"data row3 col5\" >-15.42</td>\n",
       "      <td id=\"T_9af88_row3_col6\" class=\"data row3 col6\" >-0.10</td>\n",
       "      <td id=\"T_9af88_row3_col7\" class=\"data row3 col7\" >-0.14</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_9af88_level0_row4\" class=\"row_heading level0 row4\" >d_connector_pad_Q3</th>\n",
       "      <td id=\"T_9af88_row4_col0\" class=\"data row4 col0\" >-0.18</td>\n",
       "      <td id=\"T_9af88_row4_col1\" class=\"data row4 col1\" >-0.76</td>\n",
       "      <td id=\"T_9af88_row4_col2\" class=\"data row4 col2\" >-0.42</td>\n",
       "      <td id=\"T_9af88_row4_col3\" class=\"data row4 col3\" >-0.16</td>\n",
       "      <td id=\"T_9af88_row4_col4\" class=\"data row4 col4\" >52.37</td>\n",
       "      <td id=\"T_9af88_row4_col5\" class=\"data row4 col5\" >-34.91</td>\n",
       "      <td id=\"T_9af88_row4_col6\" class=\"data row4 col6\" >-1.52</td>\n",
       "      <td id=\"T_9af88_row4_col7\" class=\"data row4 col7\" >-13.48</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_9af88_level0_row5\" class=\"row_heading level0 row5\" >ground_main_plane</th>\n",
       "      <td id=\"T_9af88_row5_col0\" class=\"data row5 col0\" >-34.97</td>\n",
       "      <td id=\"T_9af88_row5_col1\" class=\"data row5 col1\" >-34.98</td>\n",
       "      <td id=\"T_9af88_row5_col2\" class=\"data row5 col2\" >-34.94</td>\n",
       "      <td id=\"T_9af88_row5_col3\" class=\"data row5 col3\" >-15.42</td>\n",
       "      <td id=\"T_9af88_row5_col4\" class=\"data row5 col4\" >-34.91</td>\n",
       "      <td id=\"T_9af88_row5_col5\" class=\"data row5 col5\" >270.02</td>\n",
       "      <td id=\"T_9af88_row5_col6\" class=\"data row5 col6\" >-27.93</td>\n",
       "      <td id=\"T_9af88_row5_col7\" class=\"data row5 col7\" >-27.91</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_9af88_level0_row6\" class=\"row_heading level0 row6\" >pad_bot_Q3</th>\n",
       "      <td id=\"T_9af88_row6_col0\" class=\"data row6 col0\" >-13.47</td>\n",
       "      <td id=\"T_9af88_row6_col1\" class=\"data row6 col1\" >-1.52</td>\n",
       "      <td id=\"T_9af88_row6_col2\" class=\"data row6 col2\" >-13.48</td>\n",
       "      <td id=\"T_9af88_row6_col3\" class=\"data row6 col3\" >-0.10</td>\n",
       "      <td id=\"T_9af88_row6_col4\" class=\"data row6 col4\" >-1.52</td>\n",
       "      <td id=\"T_9af88_row6_col5\" class=\"data row6 col5\" >-27.93</td>\n",
       "      <td id=\"T_9af88_row6_col6\" class=\"data row6 col6\" >82.20</td>\n",
       "      <td id=\"T_9af88_row6_col7\" class=\"data row6 col7\" >-22.03</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_9af88_level0_row7\" class=\"row_heading level0 row7\" >pad_top_Q3</th>\n",
       "      <td id=\"T_9af88_row7_col0\" class=\"data row7 col0\" >-1.52</td>\n",
       "      <td id=\"T_9af88_row7_col1\" class=\"data row7 col1\" >-13.48</td>\n",
       "      <td id=\"T_9af88_row7_col2\" class=\"data row7 col2\" >-1.52</td>\n",
       "      <td id=\"T_9af88_row7_col3\" class=\"data row7 col3\" >-0.14</td>\n",
       "      <td id=\"T_9af88_row7_col4\" class=\"data row7 col4\" >-13.48</td>\n",
       "      <td id=\"T_9af88_row7_col5\" class=\"data row7 col5\" >-27.91</td>\n",
       "      <td id=\"T_9af88_row7_col6\" class=\"data row7 col6\" >-22.03</td>\n",
       "      <td id=\"T_9af88_row7_col7\" class=\"data row7 col7\" >82.21</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n"
      ],
      "text/plain": [
       "<pandas.io.formats.style.Styler at 0x7fc77ac3b250>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Q3\n",
    "path3 = \"./capacitance_matrices/Q3_cap_matrix.txt\"\n",
    "t3_mat, _, _, _ = load_q3d_capacitance_matrix(path3)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "1b86aec8-1a07-40ba-8627-905de7601360",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Imported capacitance matrix with UNITS: [fF] now converted to USER UNITS:[fF]                 from file:\n",
      "\t./capacitance_matrices/Q4_cap_matrix.txt\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<style type=\"text/css\">\n",
       "#T_cebf3_row0_col0, #T_cebf3_row1_col1, #T_cebf3_row2_col2, #T_cebf3_row3_col3, #T_cebf3_row4_col4, #T_cebf3_row5_col5, #T_cebf3_row6_col6 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 100.0%, transparent 100.0%);\n",
       "}\n",
       "#T_cebf3_row0_col1 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 39.8%, transparent 39.8%);\n",
       "}\n",
       "#T_cebf3_row0_col2, #T_cebf3_row2_col0 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 39.2%, transparent 39.2%);\n",
       "}\n",
       "#T_cebf3_row0_col3, #T_cebf3_row1_col3 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 49.3%, transparent 49.3%);\n",
       "}\n",
       "#T_cebf3_row0_col4, #T_cebf3_row2_col4 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 0.1%, transparent 0.1%);\n",
       "}\n",
       "#T_cebf3_row0_col5, #T_cebf3_row2_col5 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 13.3%, transparent 13.3%);\n",
       "}\n",
       "#T_cebf3_row0_col6 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 29.0%, transparent 29.0%);\n",
       "}\n",
       "#T_cebf3_row1_col0 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 39.6%, transparent 39.6%);\n",
       "}\n",
       "#T_cebf3_row1_col2 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 40.0%, transparent 40.0%);\n",
       "}\n",
       "#T_cebf3_row1_col4, #T_cebf3_row4_col0, #T_cebf3_row4_col1, #T_cebf3_row4_col2, #T_cebf3_row4_col3, #T_cebf3_row4_col5, #T_cebf3_row4_col6 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "}\n",
       "#T_cebf3_row1_col5 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 24.3%, transparent 24.3%);\n",
       "}\n",
       "#T_cebf3_row1_col6 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 17.5%, transparent 17.5%);\n",
       "}\n",
       "#T_cebf3_row2_col1, #T_cebf3_row3_col0, #T_cebf3_row3_col2 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 40.1%, transparent 40.1%);\n",
       "}\n",
       "#T_cebf3_row2_col3 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 49.1%, transparent 49.1%);\n",
       "}\n",
       "#T_cebf3_row2_col6 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 28.9%, transparent 28.9%);\n",
       "}\n",
       "#T_cebf3_row3_col1 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 40.3%, transparent 40.3%);\n",
       "}\n",
       "#T_cebf3_row3_col4 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 7.2%, transparent 7.2%);\n",
       "}\n",
       "#T_cebf3_row3_col5 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 25.6%, transparent 25.6%);\n",
       "}\n",
       "#T_cebf3_row3_col6 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 30.3%, transparent 30.3%);\n",
       "}\n",
       "#T_cebf3_row5_col0, #T_cebf3_row5_col2 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 24.8%, transparent 24.8%);\n",
       "}\n",
       "#T_cebf3_row5_col1 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 38.5%, transparent 38.5%);\n",
       "}\n",
       "#T_cebf3_row5_col3 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 49.0%, transparent 49.0%);\n",
       "}\n",
       "#T_cebf3_row5_col4 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 2.7%, transparent 2.7%);\n",
       "}\n",
       "#T_cebf3_row5_col6 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 9.9%, transparent 9.9%);\n",
       "}\n",
       "#T_cebf3_row6_col0, #T_cebf3_row6_col2 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 38.3%, transparent 38.3%);\n",
       "}\n",
       "#T_cebf3_row6_col1 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 24.6%, transparent 24.6%);\n",
       "}\n",
       "#T_cebf3_row6_col3 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 48.8%, transparent 48.8%);\n",
       "}\n",
       "#T_cebf3_row6_col4 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 1.1%, transparent 1.1%);\n",
       "}\n",
       "#T_cebf3_row6_col5 {\n",
       "  width: 10em;\n",
       "  height: 80%;\n",
       "  background: linear-gradient(90deg,#5fba7d 5.5%, transparent 5.5%);\n",
       "}\n",
       "</style>\n",
       "<table id=\"T_cebf3_\">\n",
       "  <thead>\n",
       "    <tr>\n",
       "      <th class=\"blank level0\" >&nbsp;</th>\n",
       "      <th class=\"col_heading level0 col0\" >a_connector_pad_Q4</th>\n",
       "      <th class=\"col_heading level0 col1\" >b_connector_pad_Q4</th>\n",
       "      <th class=\"col_heading level0 col2\" >c_connector_pad_Q4</th>\n",
       "      <th class=\"col_heading level0 col3\" >cl_metal_Q4</th>\n",
       "      <th class=\"col_heading level0 col4\" >ground_main_plane</th>\n",
       "      <th class=\"col_heading level0 col5\" >pad_bot_Q4</th>\n",
       "      <th class=\"col_heading level0 col6\" >pad_top_Q4</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th id=\"T_cebf3_level0_row0\" class=\"row_heading level0 row0\" >a_connector_pad_Q4</th>\n",
       "      <td id=\"T_cebf3_row0_col0\" class=\"data row0 col0\" >52.54</td>\n",
       "      <td id=\"T_cebf3_row0_col1\" class=\"data row0 col1\" >-0.44</td>\n",
       "      <td id=\"T_cebf3_row0_col2\" class=\"data row0 col2\" >-0.84</td>\n",
       "      <td id=\"T_cebf3_row0_col3\" class=\"data row0 col3\" >-0.01</td>\n",
       "      <td id=\"T_cebf3_row0_col4\" class=\"data row0 col4\" >-35.22</td>\n",
       "      <td id=\"T_cebf3_row0_col5\" class=\"data row0 col5\" >-13.46</td>\n",
       "      <td id=\"T_cebf3_row0_col6\" class=\"data row0 col6\" >-1.60</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_cebf3_level0_row1\" class=\"row_heading level0 row1\" >b_connector_pad_Q4</th>\n",
       "      <td id=\"T_cebf3_row1_col0\" class=\"data row1 col0\" >-0.44</td>\n",
       "      <td id=\"T_cebf3_row1_col1\" class=\"data row1 col1\" >52.44</td>\n",
       "      <td id=\"T_cebf3_row1_col2\" class=\"data row1 col2\" >-0.20</td>\n",
       "      <td id=\"T_cebf3_row1_col3\" class=\"data row1 col3\" >-0.01</td>\n",
       "      <td id=\"T_cebf3_row1_col4\" class=\"data row1 col4\" >-35.42</td>\n",
       "      <td id=\"T_cebf3_row1_col5\" class=\"data row1 col5\" >-1.56</td>\n",
       "      <td id=\"T_cebf3_row1_col6\" class=\"data row1 col6\" >-13.82</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_cebf3_level0_row2\" class=\"row_heading level0 row2\" >c_connector_pad_Q4</th>\n",
       "      <td id=\"T_cebf3_row2_col0\" class=\"data row2 col0\" >-0.84</td>\n",
       "      <td id=\"T_cebf3_row2_col1\" class=\"data row2 col1\" >-0.20</td>\n",
       "      <td id=\"T_cebf3_row2_col2\" class=\"data row2 col2\" >52.48</td>\n",
       "      <td id=\"T_cebf3_row2_col3\" class=\"data row2 col3\" >-0.06</td>\n",
       "      <td id=\"T_cebf3_row2_col4\" class=\"data row2 col4\" >-35.26</td>\n",
       "      <td id=\"T_cebf3_row2_col5\" class=\"data row2 col5\" >-13.47</td>\n",
       "      <td id=\"T_cebf3_row2_col6\" class=\"data row2 col6\" >-1.67</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_cebf3_level0_row3\" class=\"row_heading level0 row3\" >cl_metal_Q4</th>\n",
       "      <td id=\"T_cebf3_row3_col0\" class=\"data row3 col0\" >-0.01</td>\n",
       "      <td id=\"T_cebf3_row3_col1\" class=\"data row3 col1\" >-0.01</td>\n",
       "      <td id=\"T_cebf3_row3_col2\" class=\"data row3 col2\" >-0.06</td>\n",
       "      <td id=\"T_cebf3_row3_col3\" class=\"data row3 col3\" >16.01</td>\n",
       "      <td id=\"T_cebf3_row3_col4\" class=\"data row3 col4\" >-15.59</td>\n",
       "      <td id=\"T_cebf3_row3_col5\" class=\"data row3 col5\" >-0.10</td>\n",
       "      <td id=\"T_cebf3_row3_col6\" class=\"data row3 col6\" >-0.16</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_cebf3_level0_row4\" class=\"row_heading level0 row4\" >ground_main_plane</th>\n",
       "      <td id=\"T_cebf3_row4_col0\" class=\"data row4 col0\" >-35.22</td>\n",
       "      <td id=\"T_cebf3_row4_col1\" class=\"data row4 col1\" >-35.42</td>\n",
       "      <td id=\"T_cebf3_row4_col2\" class=\"data row4 col2\" >-35.26</td>\n",
       "      <td id=\"T_cebf3_row4_col3\" class=\"data row4 col3\" >-15.59</td>\n",
       "      <td id=\"T_cebf3_row4_col4\" class=\"data row4 col4\" >241.52</td>\n",
       "      <td id=\"T_cebf3_row4_col5\" class=\"data row4 col5\" >-27.95</td>\n",
       "      <td id=\"T_cebf3_row4_col6\" class=\"data row4 col6\" >-32.50</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_cebf3_level0_row5\" class=\"row_heading level0 row5\" >pad_bot_Q4</th>\n",
       "      <td id=\"T_cebf3_row5_col0\" class=\"data row5 col0\" >-13.46</td>\n",
       "      <td id=\"T_cebf3_row5_col1\" class=\"data row5 col1\" >-1.56</td>\n",
       "      <td id=\"T_cebf3_row5_col2\" class=\"data row5 col2\" >-13.47</td>\n",
       "      <td id=\"T_cebf3_row5_col3\" class=\"data row5 col3\" >-0.10</td>\n",
       "      <td id=\"T_cebf3_row5_col4\" class=\"data row5 col4\" >-27.95</td>\n",
       "      <td id=\"T_cebf3_row5_col5\" class=\"data row5 col5\" >80.67</td>\n",
       "      <td id=\"T_cebf3_row5_col6\" class=\"data row5 col6\" >-21.97</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th id=\"T_cebf3_level0_row6\" class=\"row_heading level0 row6\" >pad_top_Q4</th>\n",
       "      <td id=\"T_cebf3_row6_col0\" class=\"data row6 col0\" >-1.60</td>\n",
       "      <td id=\"T_cebf3_row6_col1\" class=\"data row6 col1\" >-13.82</td>\n",
       "      <td id=\"T_cebf3_row6_col2\" class=\"data row6 col2\" >-1.67</td>\n",
       "      <td id=\"T_cebf3_row6_col3\" class=\"data row6 col3\" >-0.16</td>\n",
       "      <td id=\"T_cebf3_row6_col4\" class=\"data row6 col4\" >-32.50</td>\n",
       "      <td id=\"T_cebf3_row6_col5\" class=\"data row6 col5\" >-21.97</td>\n",
       "      <td id=\"T_cebf3_row6_col6\" class=\"data row6 col6\" >74.13</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n"
      ],
      "text/plain": [
       "<pandas.io.formats.style.Styler at 0x7fc77a71f8e0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Q4\n",
    "path4 = \"./capacitance_matrices/Q4_cap_matrix.txt\"\n",
    "t4_mat, _, _, _ = load_q3d_capacitance_matrix(path4)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "7c1fd689-e715-479b-b0f1-6c952c91f76c",
   "metadata": {},
   "source": [
    "## Create LOM Cells from Capacitance Matrix"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "64d38655-9b20-488f-96f8-b8768a8b6554",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Cell 2: Transmon-2\n",
    "opt2 = dict(\n",
    "    node_rename = {'b_connector_pad_Q2': 'coupler23',\n",
    "                   'c_connector_pad_Q2': 'readout2'}, \n",
    "    cap_mat = t2_mat,\n",
    "    ind_dict = {('pad_top_Q2', 'pad_bot_Q2'): 12},  # junction inductance in nH\n",
    "    jj_dict = {('pad_top_Q2', 'pad_bot_Q2'): 'j2'},\n",
    "    cj_dict = {('pad_top_Q2', 'pad_bot_Q2'): 1}, # junction capacitance in fF\n",
    "\n",
    ")\n",
    "cell_2 = Cell(opt2)\n",
    "\n",
    "# Cell 3: Transmon-3\n",
    "opt3 = dict(\n",
    "    node_rename = {'b_connector_pad_Q3': 'coupler23',\n",
    "                   'c_connector_pad_Q3': 'readout3'}, \n",
    "    cap_mat = t3_mat,\n",
    "    ind_dict = {('pad_top_Q3', 'pad_bot_Q3'): 12},  # junction inductance in nH\n",
    "    jj_dict = {('pad_top_Q3', 'pad_bot_Q3'): 'j3'},\n",
    "    cj_dict = {('pad_top_Q3', 'pad_bot_Q3'): 1}, # junction capacitance in fF\n",
    "\n",
    ")\n",
    "cell_3 = Cell(opt3)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1c6f6e43-32be-4216-9a12-0b8dc298015f",
   "metadata": {},
   "source": [
    "## Create Qubit Subsystems"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "id": "ba9b3e3f-69a3-4ca7-9372-28cd51f7163a",
   "metadata": {},
   "outputs": [],
   "source": [
    "# subsystem 1: Transmon-2\n",
    "transmon2 = Subsystem(name='transmon2', sys_type='TRANSMON', nodes=['j2'])\n",
    "\n",
    "# subsystem 2: Transmon-3\n",
    "transmon3 = Subsystem(name='transmon3', sys_type='TRANSMON', nodes=['j3'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "id": "e20d9fa4-d046-4ba2-916a-100ac6d95d22",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Resonator Subsystems\n",
    "q_opts = dict(\n",
    "    Z0 = 50, # characteristic impedance in Ohm\n",
    "    vp = 0.404314 * c_light # phase velocity \n",
    ")\n",
    "\n",
    "# subsystem 3: Q2 readout resonator\n",
    "ro2 = Subsystem(name='readout2', sys_type='TL_RESONATOR', nodes=['readout2'], q_opts=dict(f_res = 6.96, **q_opts))\n",
    "\n",
    "# subsystem 4: Q3 readout resonator\n",
    "ro3 = Subsystem(name='readout3', sys_type='TL_RESONATOR', nodes=['readout3'], q_opts=dict(f_res = 6.98, **q_opts))\n",
    "\n",
    "# subsystem 15: Q2-Q3 coupling resonator\n",
    "coup23 = Subsystem(name='coupler23', sys_type='TL_RESONATOR', nodes=['coupler23'], q_opts=dict(f_res = 7.36, **q_opts))"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "cfcb2fa2-66de-4719-a281-626dc9df2888",
   "metadata": {},
   "source": [
    "# Create Composite System"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "id": "8d8c7099-7896-4dbf-9fdd-2015984c9f38",
   "metadata": {},
   "outputs": [],
   "source": [
    "composite_sys = CompositeSystem(\n",
    "    subsystems=[transmon2, transmon3, coup23, ro2, ro3], \n",
    "    cells=[cell_2, cell_3],\n",
    "    grd_node='ground_main_plane',\n",
    "    nodes_force_keep=['readout2', 'readout3', 'coupler23']\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "8acc3ad6-9408-4e1e-8c39-6660437306f9",
   "metadata": {},
   "source": [
    "# Print the Circuit Graph\n",
    "This method generates a graph consisting of different circuit nodes and then reduces this to form an equivalent capacitance matrix only consisting of the reduced set of components."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "id": "61b28436-5540-413f-b2a2-e145ce979ebf",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "node_jj_basis:\n",
      "-------------\n",
      "\n",
      "['j2', 'pad_bot_Q2', 'j3', 'pad_bot_Q3', 'coupler23', 'readout2', 'readout3', 'a_connector_pad_Q2', 'a_connector_pad_Q3', 'cl_metal_Q2', 'cl_metal_Q3', 'd_connector_pad_Q2', 'd_connector_pad_Q3']\n",
      "\n",
      "nodes_keep:\n",
      "-------------\n",
      "\n",
      "['j2', 'j3', 'coupler23', 'readout2', 'readout3']\n",
      "\n",
      "\n",
      "L_inv_k (reduced inverse inductance matrix):\n",
      "-------------\n",
      "\n",
      "                 j2        j3  coupler23  readout2  readout3\n",
      "j2         0.083333  0.000000        0.0       0.0       0.0\n",
      "j3         0.000000  0.083333        0.0       0.0       0.0\n",
      "coupler23  0.000000  0.000000        0.0       0.0       0.0\n",
      "readout2   0.000000  0.000000        0.0       0.0       0.0\n",
      "readout3   0.000000  0.000000        0.0       0.0       0.0\n",
      "\n",
      "C_k (reduced capacitance matrix):\n",
      "-------------\n",
      "\n",
      "                  j2         j3  coupler23   readout2   readout3\n",
      "j2         54.614341   0.000000  -6.022198   6.023767   0.000000\n",
      "j3          0.000000  51.756296  -6.018143   0.000000   6.015938\n",
      "coupler23  -6.022198  -6.018143  99.571810  -2.180863  -2.301963\n",
      "readout2    6.023767   0.000000  -2.180863  49.367061   0.000000\n",
      "readout3    0.000000   6.015938  -2.301963   0.000000  50.186987\n",
      "\n",
      "\n"
     ]
    }
   ],
   "source": [
    "cg = composite_sys.circuitGraph()\n",
    "print(cg)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "b353a4a0-97ee-42af-b955-dc789274c4fc",
   "metadata": {},
   "source": [
    "# Get Hilbert space of the Composite System"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "id": "1051b2a3-b43a-40a9-9563-f3663ccff988",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "HilbertSpace:  subsystems\n",
      "-------------------------\n",
      "\n",
      "Transmon------------| [Transmon_5]\n",
      "                    | EJ: 13621.792733898432\n",
      "                    | EC: 361.764776983381\n",
      "                    | ng: 0.001\n",
      "                    | ncut: 22\n",
      "                    | truncated_dim: 10\n",
      "                    |\n",
      "                    | dim: 45\n",
      "\n",
      "\n",
      "Transmon------------| [Transmon_6]\n",
      "                    | EJ: 13621.792733898432\n",
      "                    | EC: 382.04618957580647\n",
      "                    | ng: 0.001\n",
      "                    | ncut: 22\n",
      "                    | truncated_dim: 10\n",
      "                    |\n",
      "                    | dim: 45\n",
      "\n",
      "\n",
      "Oscillator----------| [Oscillator_1]\n",
      "                    | E_osc: 7360.0\n",
      "                    | l_osc: None\n",
      "                    | truncated_dim: 3\n",
      "                    |\n",
      "                    | dim: 3\n",
      "\n",
      "\n",
      "Oscillator----------| [Oscillator_2]\n",
      "                    | E_osc: 6960.0\n",
      "                    | l_osc: None\n",
      "                    | truncated_dim: 3\n",
      "                    |\n",
      "                    | dim: 3\n",
      "\n",
      "\n",
      "Oscillator----------| [Oscillator_3]\n",
      "                    | E_osc: 6980.0\n",
      "                    | l_osc: None\n",
      "                    | truncated_dim: 3\n",
      "                    |\n",
      "                    | dim: 3\n",
      "\n",
      "\n"
     ]
    }
   ],
   "source": [
    "hilbertspace = composite_sys.create_hilbertspace()\n",
    "print(hilbertspace)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "id": "595a9c9b-40fb-4675-9d4a-b6cd03e58cb8",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/latex": [
       "Quantum object: dims = [[10, 10, 3, 3, 3], [10, 10, 3, 3, 3]], shape = (2700, 2700), type = oper, isherm = True\\begin{equation*}\\left(\\begin{array}{*{11}c}-2.107\\times10^{+04} & 0.097j & 0.0 & 0.092j & -0.119 & \\cdots & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\-0.097j & -1.409\\times10^{+04} & 0.138j & 0.119 & 0.092j & \\cdots & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\0.0 & -0.138j & -7.110\\times10^{+03} & 0.0 & 0.168 & \\cdots & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\-0.092j & 0.119 & 0.0 & -1.411\\times10^{+04} & 0.097j & \\cdots & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\-0.119 & -0.092j & 0.168 & -0.097j & -7.130\\times10^{+03} & \\cdots & 0.0 & 0.0 & 0.0 & 0.0 & 0.0\\\\\\vdots & \\vdots & \\vdots & \\vdots & \\vdots & \\ddots & \\vdots & \\vdots & \\vdots & \\vdots & \\vdots\\\\0.0 & 0.0 & 0.0 & 0.0 & 0.0 & \\cdots & 1.047\\times10^{+05} & 676.309j & 0.168 & 638.773j & -0.238\\\\0.0 & 0.0 & 0.0 & 0.0 & 0.0 & \\cdots & -676.309j & 1.117\\times10^{+05} & 0.0 & 0.238 & 638.773j\\\\0.0 & 0.0 & 0.0 & 0.0 & 0.0 & \\cdots & 0.168 & 0.0 & 1.047\\times10^{+05} & 478.223j & 0.0\\\\0.0 & 0.0 & 0.0 & 0.0 & 0.0 & \\cdots & -638.773j & 0.238 & -478.223j & 1.117\\times10^{+05} & 676.309j\\\\0.0 & 0.0 & 0.0 & 0.0 & 0.0 & \\cdots & -0.238 & -638.773j & 0.0 & -676.309j & 1.187\\times10^{+05}\\\\\\end{array}\\right)\\end{equation*}"
      ],
      "text/plain": [
       "Quantum object: dims = [[10, 10, 3, 3, 3], [10, 10, 3, 3, 3]], shape = (2700, 2700), type = oper, isherm = True\n",
       "Qobj data =\n",
       "[[-21069.94377723+0.00000000e+00j      0.        +9.72724112e-02j\n",
       "       0.        +0.00000000e+00j ...      0.        +0.00000000e+00j\n",
       "       0.        +0.00000000e+00j      0.        +0.00000000e+00j]\n",
       " [     0.        -9.72724112e-02j -14089.94377723+0.00000000e+00j\n",
       "       0.        +1.37563963e-01j ...      0.        +0.00000000e+00j\n",
       "       0.        +0.00000000e+00j      0.        +0.00000000e+00j]\n",
       " [     0.        +0.00000000e+00j      0.        -1.37563963e-01j\n",
       "   -7109.94377723+0.00000000e+00j ...      0.        +0.00000000e+00j\n",
       "       0.        +0.00000000e+00j      0.        +0.00000000e+00j]\n",
       " ...\n",
       " [     0.        +0.00000000e+00j      0.        +0.00000000e+00j\n",
       "       0.        +0.00000000e+00j ... 104725.46523701+0.00000000e+00j\n",
       "       0.        +4.78222943e+02j      0.        +0.00000000e+00j]\n",
       " [     0.        +0.00000000e+00j      0.        +0.00000000e+00j\n",
       "       0.        +0.00000000e+00j ...      0.        -4.78222943e+02j\n",
       "  111705.46523701+0.00000000e+00j      0.        +6.76309372e+02j]\n",
       " [     0.        +0.00000000e+00j      0.        +0.00000000e+00j\n",
       "       0.        +0.00000000e+00j ...      0.        +0.00000000e+00j\n",
       "       0.        -6.76309372e+02j 118685.46523701+0.00000000e+00j]]"
      ]
     },
     "execution_count": 25,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Convert the hilbert space into\n",
    "# \"Interaction Picture\"\n",
    "hilbertspace = composite_sys.add_interaction()\n",
    "hilbertspace.hamiltonian()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "id": "16ac830a-8498-40f7-9a9e-55b51c5abab4",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Finished eigensystem.     \n",
      "\n",
      "system frequencies in GHz:\n",
      "--------------------------\n",
      "{'transmon2': 5.875345402531805, 'transmon3': 6.024975402956764, 'coupler23': 7.372511788398417, 'readout2': 6.96758399935126, 'readout3': 6.989434164823506}\n",
      "\n",
      "Chi matrices in MHz\n",
      "--------------------------\n",
      "            transmon2   transmon3  coupler23      readout2      readout3\n",
      "transmon2 -418.646319    1.502700  -3.407336 -5.453703e+00 -4.588201e-02\n",
      "transmon3    1.502700 -447.200698  -4.470614 -3.150036e-02 -7.439598e+00\n",
      "coupler23   -3.407336   -4.470614   4.225685 -3.458875e-02 -6.315097e-02\n",
      "readout2    -5.453703   -0.031500  -0.034589  2.032900e+00 -1.216649e-07\n",
      "readout3    -0.045882   -7.439598  -0.063151 -1.216649e-07  2.134685e+00\n"
     ]
    }
   ],
   "source": [
    "hamiltonian_results = composite_sys.hamiltonian_results(hilbertspace, evals_count=30)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "id": "dda41b37-82b1-4a39-a688-cc56bc5c6412",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "           transmon2   transmon3   coupler23   readout2   readout3\n",
       "transmon2   0.000000   18.985900   95.566620 -93.568702   1.506729\n",
       "transmon3  18.985900    0.000000  100.649134   1.497971 -98.781480\n",
       "coupler23  95.566620  100.649134    0.000000   7.540125   7.987558\n",
       "readout2  -93.568702    1.497971    7.540125   0.000000   0.118880\n",
       "readout3    1.506729  -98.781480    7.987558   0.118880   0.000000"
      ]
     },
     "execution_count": 27,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "composite_sys.compute_gs()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "4efcc1ba-573f-44e0-bef0-8704ab5c043c",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "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.8.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}