Aladdin passive
Passive harmonic cavity case w/ standard Aladdin lattice:
Fig. 1 (a) & 1 (b) of [1]
Fig. 3 (a) of [2]
[1]: Bosch, R. A., Kleman, K. J., & Bisognano, J. J. (2001). Robinson instabilities with a higher-harmonic cavity. Physical Review Special Topics-Accelerators and Beams, 4(7), 074401.
[2]: Bosch, R. A. (2005). Instabilities driven by higher-order modes in a radio frequency system with a passive higher-harmonic cavity. Physical Review Special Topics-Accelerators and Beams, 8(8), 084401.
Imports
import numpy as np
import matplotlib.pyplot as plt
from mbtrack2 import CavityResonator
from aladdin import aladdin
from albums.scan import scan_psi_I0
from albums import SystemState
from albums import EquilibriumOptions, TheoryOptions
from albums import DEFAULT_FLOWS
Ring and cavity setup
ring = aladdin(mode="standard", V_RF=90e3)
MC = CavityResonator(ring, m=1, Rs=0.5e6, Q=8000, QL=8000, detune=100e3)
MC.beta = 11
MC.Vc = 90e3
MC.theta = np.arccos(ring.U0/MC.Vc)
HC = CavityResonator(ring, m=4, Rs=1.24e6, Q=20250, QL=20250, detune=100e3)
Benchmark
points = 50
currents = np.linspace(0, 300e-3, points)
psi_HC_vals = np.linspace(80, 90, points)
method = "Bosch"
flow = DEFAULT_FLOWS[method]
flow_no_coupling = DEFAULT_FLOWS[f"{method}_no_coupling"]
N_sigma = 20
tau_boundary = ring.sigma_0*N_sigma
eq_opts = EquilibriumOptions(tau_boundary=tau_boundary,
passive_harmonic_cavity=True,
auto_set_MC_theta=True,
optimal_tunning=True)
Passive harmonic cavity - Robinson instabilities without coupling
Fig. 1 (a) in [1]
result = scan_psi_I0(
SystemState(ring, [MC, HC]),
psi_HC_vals,
currents,
flow=flow_no_coupling,
eq_opts=eq_opts,
theory_opts=TheoryOptions(mode_coupling=False))
result.plot("xi")
scan: 0%| | 0/2500 [00:00<?, ?it/s]
scan: 1%| | 22/2500 [00:00<00:11, 212.06it/s]
scan: 2%|▏ | 44/2500 [00:00<00:12, 195.30it/s]
scan: 3%|▎ | 64/2500 [00:00<00:12, 194.97it/s]
scan: 3%|▎ | 84/2500 [00:00<00:12, 194.78it/s]
scan: 4%|▍ | 104/2500 [00:00<00:12, 193.32it/s]
scan: 5%|▍ | 124/2500 [00:00<00:12, 195.31it/s]
scan: 6%|▌ | 144/2500 [00:00<00:12, 191.85it/s]
scan: 7%|▋ | 164/2500 [00:00<00:12, 191.96it/s]
scan: 7%|▋ | 187/2500 [00:00<00:11, 202.76it/s]
scan: 8%|▊ | 208/2500 [00:01<00:11, 199.18it/s]
scan: 9%|▉ | 230/2500 [00:01<00:11, 204.71it/s]
scan: 10%|█ | 251/2500 [00:01<00:11, 197.15it/s]
scan: 11%|█ | 272/2500 [00:01<00:11, 199.93it/s]
scan: 12%|█▏ | 295/2500 [00:01<00:10, 204.44it/s]
scan: 13%|█▎ | 316/2500 [00:01<00:10, 202.98it/s]
scan: 14%|█▎ | 339/2500 [00:01<00:10, 210.06it/s]
scan: 14%|█▍ | 361/2500 [00:01<00:10, 206.62it/s]
scan: 15%|█▌ | 382/2500 [00:01<00:10, 198.33it/s]
scan: 16%|█▌ | 403/2500 [00:02<00:10, 198.83it/s]
scan: 17%|█▋ | 424/2500 [00:02<00:10, 201.46it/s]
scan: 18%|█▊ | 445/2500 [00:02<00:10, 191.65it/s]
scan: 19%|█▊ | 465/2500 [00:02<00:10, 193.40it/s]
scan: 20%|█▉ | 489/2500 [00:02<00:09, 204.30it/s]
scan: 20%|██ | 510/2500 [00:02<00:10, 194.84it/s]
scan: 21%|██▏ | 533/2500 [00:02<00:09, 203.29it/s]
scan: 22%|██▏ | 554/2500 [00:02<00:09, 199.25it/s]
scan: 23%|██▎ | 576/2500 [00:02<00:09, 202.95it/s]
scan: 24%|██▍ | 597/2500 [00:03<00:09, 194.27it/s]
scan: 25%|██▍ | 617/2500 [00:03<00:09, 195.02it/s]
scan: 26%|██▌ | 640/2500 [00:03<00:09, 204.72it/s]
scan: 26%|██▋ | 661/2500 [00:03<00:09, 196.15it/s]
scan: 27%|██▋ | 684/2500 [00:03<00:08, 204.82it/s]
scan: 28%|██▊ | 706/2500 [00:03<00:08, 207.47it/s]
scan: 29%|██▉ | 727/2500 [00:03<00:08, 207.85it/s]
scan: 30%|██▉ | 748/2500 [00:03<00:08, 207.07it/s]
scan: 31%|███ | 769/2500 [00:03<00:08, 205.07it/s]
scan: 32%|███▏ | 793/2500 [00:03<00:07, 213.58it/s]
scan: 33%|███▎ | 815/2500 [00:04<00:08, 202.51it/s]
scan: 34%|███▎ | 839/2500 [00:04<00:07, 210.76it/s]
scan: 34%|███▍ | 861/2500 [00:04<00:07, 208.04it/s]
scan: 35%|███▌ | 884/2500 [00:04<00:07, 212.33it/s]
scan: 36%|███▌ | 906/2500 [00:04<00:07, 212.26it/s]
scan: 37%|███▋ | 928/2500 [00:04<00:07, 210.97it/s]
scan: 38%|███▊ | 952/2500 [00:04<00:07, 216.68it/s]
scan: 39%|███▉ | 974/2500 [00:04<00:07, 212.03it/s]
scan: 40%|███▉ | 998/2500 [00:04<00:06, 218.94it/s]
scan: 41%|████ | 1020/2500 [00:05<00:06, 213.67it/s]
scan: 42%|████▏ | 1044/2500 [00:05<00:06, 219.52it/s]
scan: 43%|████▎ | 1066/2500 [00:05<00:06, 214.16it/s]
scan: 44%|████▎ | 1089/2500 [00:05<00:06, 217.99it/s]
scan: 44%|████▍ | 1111/2500 [00:05<00:06, 213.94it/s]
scan: 45%|████▌ | 1133/2500 [00:05<00:06, 213.90it/s]
scan: 46%|████▌ | 1156/2500 [00:05<00:06, 217.47it/s]
scan: 47%|████▋ | 1178/2500 [00:05<00:06, 213.11it/s]
scan: 48%|████▊ | 1202/2500 [00:05<00:05, 220.15it/s]
scan: 49%|████▉ | 1225/2500 [00:05<00:05, 213.57it/s]
scan: 50%|████▉ | 1249/2500 [00:06<00:05, 220.04it/s]
scan: 51%|█████ | 1272/2500 [00:06<00:05, 212.37it/s]
scan: 52%|█████▏ | 1295/2500 [00:06<00:05, 216.14it/s]
scan: 53%|█████▎ | 1317/2500 [00:06<00:05, 210.52it/s]
scan: 54%|█████▎ | 1339/2500 [00:06<00:05, 213.07it/s]
scan: 54%|█████▍ | 1361/2500 [00:06<00:05, 212.45it/s]
scan: 55%|█████▌ | 1383/2500 [00:06<00:05, 210.27it/s]
scan: 56%|█████▌ | 1406/2500 [00:06<00:05, 214.28it/s]
scan: 57%|█████▋ | 1428/2500 [00:06<00:05, 209.04it/s]
scan: 58%|█████▊ | 1451/2500 [00:07<00:04, 214.72it/s]
scan: 59%|█████▉ | 1473/2500 [00:07<00:04, 207.94it/s]
scan: 60%|█████▉ | 1495/2500 [00:07<00:04, 210.44it/s]
scan: 61%|██████ | 1517/2500 [00:07<00:04, 206.87it/s]
scan: 62%|██████▏ | 1538/2500 [00:07<00:04, 205.57it/s]
scan: 62%|██████▏ | 1560/2500 [00:07<00:04, 206.02it/s]
scan: 63%|██████▎ | 1581/2500 [00:07<00:04, 204.18it/s]
scan: 64%|██████▍ | 1604/2500 [00:07<00:04, 208.70it/s]
scan: 65%|██████▌ | 1625/2500 [00:07<00:04, 202.79it/s]
scan: 66%|██████▌ | 1646/2500 [00:07<00:04, 204.62it/s]
scan: 67%|██████▋ | 1667/2500 [00:08<00:04, 202.57it/s]
scan: 68%|██████▊ | 1688/2500 [00:08<00:04, 201.34it/s]
scan: 68%|██████▊ | 1709/2500 [00:08<00:03, 203.69it/s]
scan: 69%|██████▉ | 1730/2500 [00:08<00:03, 199.63it/s]
scan: 70%|███████ | 1752/2500 [00:08<00:03, 202.52it/s]
scan: 71%|███████ | 1773/2500 [00:08<00:03, 194.97it/s]
scan: 72%|███████▏ | 1794/2500 [00:08<00:03, 197.56it/s]
scan: 73%|███████▎ | 1815/2500 [00:08<00:03, 198.22it/s]
scan: 73%|███████▎ | 1835/2500 [00:08<00:03, 196.60it/s]
scan: 74%|███████▍ | 1856/2500 [00:09<00:03, 199.49it/s]
scan: 75%|███████▌ | 1876/2500 [00:09<00:03, 196.57it/s]
scan: 76%|███████▌ | 1896/2500 [00:09<00:03, 197.18it/s]
scan: 77%|███████▋ | 1916/2500 [00:09<00:02, 197.10it/s]
scan: 77%|███████▋ | 1936/2500 [00:09<00:02, 194.71it/s]
scan: 78%|███████▊ | 1956/2500 [00:09<00:02, 192.83it/s]
scan: 79%|███████▉ | 1976/2500 [00:09<00:02, 192.35it/s]
scan: 80%|███████▉ | 1996/2500 [00:09<00:02, 192.23it/s]
scan: 81%|████████ | 2016/2500 [00:09<00:02, 193.01it/s]
scan: 81%|████████▏ | 2036/2500 [00:09<00:02, 191.62it/s]
scan: 82%|████████▏ | 2056/2500 [00:10<00:02, 193.49it/s]
scan: 83%|████████▎ | 2076/2500 [00:10<00:02, 191.68it/s]
scan: 84%|████████▍ | 2096/2500 [00:10<00:02, 190.31it/s]
scan: 85%|████████▍ | 2116/2500 [00:10<00:02, 191.37it/s]
scan: 85%|████████▌ | 2136/2500 [00:10<00:01, 189.64it/s]
scan: 86%|████████▌ | 2156/2500 [00:10<00:01, 191.24it/s]
scan: 87%|████████▋ | 2176/2500 [00:10<00:01, 190.49it/s]
scan: 88%|████████▊ | 2196/2500 [00:10<00:01, 189.56it/s]
scan: 89%|████████▊ | 2215/2500 [00:10<00:01, 189.67it/s]
scan: 89%|████████▉ | 2234/2500 [00:11<00:01, 189.00it/s]
scan: 90%|█████████ | 2253/2500 [00:11<00:01, 189.03it/s]
scan: 91%|█████████ | 2272/2500 [00:11<00:01, 188.79it/s]
scan: 92%|█████████▏| 2291/2500 [00:11<00:01, 187.88it/s]
scan: 92%|█████████▏| 2310/2500 [00:11<00:01, 188.48it/s]
scan: 93%|█████████▎| 2329/2500 [00:11<00:00, 188.07it/s]
scan: 94%|█████████▍| 2348/2500 [00:11<00:00, 185.52it/s]
scan: 95%|█████████▍| 2368/2500 [00:11<00:00, 187.21it/s]
scan: 95%|█████████▌| 2387/2500 [00:11<00:00, 187.16it/s]
scan: 96%|█████████▋| 2407/2500 [00:11<00:00, 188.57it/s]
scan: 97%|█████████▋| 2427/2500 [00:12<00:00, 189.38it/s]
scan: 98%|█████████▊| 2446/2500 [00:12<00:00, 189.18it/s]
scan: 99%|█████████▊| 2467/2500 [00:12<00:00, 192.95it/s]
scan: 100%|█████████▉| 2488/2500 [00:12<00:00, 194.62it/s]
scan: 100%|██████████| 2500/2500 [00:12<00:00, 201.45it/s]
Passive harmonic cavity - Robinson instabilities with coupling
Fig. 1 (b) in [1]
result = scan_psi_I0(
SystemState(ring, [MC, HC]),
psi_HC_vals,
currents,
flow=flow,
eq_opts=eq_opts,
theory_opts=TheoryOptions(mode_coupling=True))
result.plot("xi")
scan: 0%| | 0/2500 [00:00<?, ?it/s]
scan: 1%| | 19/2500 [00:00<00:13, 184.62it/s]
scan: 2%|▏ | 38/2500 [00:00<00:14, 164.68it/s]
scan: 2%|▏ | 55/2500 [00:00<00:18, 133.47it/s]
scan: 3%|▎ | 74/2500 [00:00<00:16, 150.22it/s]
scan: 4%|▎ | 90/2500 [00:00<00:16, 146.22it/s]
scan: 4%|▍ | 106/2500 [00:00<00:18, 131.76it/s]
scan: 5%|▌ | 126/2500 [00:00<00:15, 149.25it/s]
scan: 6%|▌ | 142/2500 [00:00<00:16, 144.32it/s]
scan: 6%|▋ | 157/2500 [00:01<00:17, 130.30it/s]
scan: 7%|▋ | 178/2500 [00:01<00:15, 150.45it/s]
scan: 8%|▊ | 194/2500 [00:01<00:15, 149.16it/s]
scan: 8%|▊ | 210/2500 [00:01<00:18, 123.74it/s]
scan: 9%|▉ | 232/2500 [00:01<00:15, 145.53it/s]
scan: 10%|▉ | 248/2500 [00:01<00:16, 135.79it/s]
scan: 11%|█ | 263/2500 [00:01<00:16, 135.40it/s]
scan: 11%|█▏ | 285/2500 [00:01<00:14, 156.85it/s]
scan: 12%|█▏ | 302/2500 [00:02<00:15, 146.03it/s]
scan: 13%|█▎ | 320/2500 [00:02<00:14, 154.76it/s]
scan: 14%|█▎ | 341/2500 [00:02<00:12, 168.89it/s]
scan: 14%|█▍ | 359/2500 [00:02<00:13, 160.23it/s]
scan: 15%|█▌ | 380/2500 [00:02<00:12, 173.60it/s]
scan: 16%|█▌ | 399/2500 [00:02<00:12, 172.26it/s]
scan: 17%|█▋ | 417/2500 [00:02<00:12, 173.28it/s]
scan: 18%|█▊ | 439/2500 [00:02<00:11, 185.95it/s]
scan: 18%|█▊ | 458/2500 [00:02<00:12, 168.87it/s]
scan: 19%|█▉ | 479/2500 [00:03<00:11, 179.77it/s]
scan: 20%|█▉ | 498/2500 [00:03<00:11, 181.35it/s]
scan: 21%|██ | 517/2500 [00:03<00:11, 177.44it/s]
scan: 22%|██▏ | 539/2500 [00:03<00:10, 188.33it/s]
scan: 22%|██▏ | 559/2500 [00:03<00:10, 182.55it/s]
scan: 23%|██▎ | 580/2500 [00:03<00:10, 189.89it/s]
scan: 24%|██▍ | 600/2500 [00:03<00:10, 181.43it/s]
scan: 25%|██▍ | 619/2500 [00:03<00:10, 183.55it/s]
scan: 26%|██▌ | 641/2500 [00:03<00:09, 193.77it/s]
scan: 26%|██▋ | 661/2500 [00:04<00:09, 186.77it/s]
scan: 27%|██▋ | 681/2500 [00:04<00:09, 190.15it/s]
scan: 28%|██▊ | 703/2500 [00:04<00:09, 196.66it/s]
scan: 29%|██▉ | 723/2500 [00:04<00:09, 194.15it/s]
scan: 30%|██▉ | 746/2500 [00:04<00:08, 202.58it/s]
scan: 31%|███ | 767/2500 [00:04<00:09, 191.22it/s]
scan: 32%|███▏ | 790/2500 [00:04<00:08, 199.70it/s]
scan: 32%|███▏ | 811/2500 [00:04<00:08, 191.86it/s]
scan: 33%|███▎ | 831/2500 [00:04<00:08, 191.53it/s]
scan: 34%|███▍ | 853/2500 [00:05<00:08, 197.45it/s]
scan: 35%|███▍ | 873/2500 [00:05<00:08, 194.15it/s]
scan: 36%|███▌ | 896/2500 [00:05<00:07, 202.93it/s]
scan: 37%|███▋ | 917/2500 [00:05<00:07, 198.50it/s]
scan: 38%|███▊ | 939/2500 [00:05<00:07, 204.41it/s]
scan: 38%|███▊ | 960/2500 [00:05<00:07, 203.08it/s]
scan: 39%|███▉ | 981/2500 [00:05<00:07, 200.36it/s]
scan: 40%|████ | 1003/2500 [00:05<00:07, 205.38it/s]
scan: 41%|████ | 1024/2500 [00:05<00:07, 200.65it/s]
scan: 42%|████▏ | 1047/2500 [00:05<00:06, 208.52it/s]
scan: 43%|████▎ | 1068/2500 [00:06<00:07, 196.31it/s]
scan: 44%|████▎ | 1090/2500 [00:06<00:06, 202.84it/s]
scan: 44%|████▍ | 1111/2500 [00:06<00:06, 202.11it/s]
scan: 45%|████▌ | 1132/2500 [00:06<00:06, 201.68it/s]
scan: 46%|████▌ | 1155/2500 [00:06<00:06, 207.28it/s]
scan: 47%|████▋ | 1176/2500 [00:06<00:06, 202.06it/s]
scan: 48%|████▊ | 1199/2500 [00:06<00:06, 209.06it/s]
scan: 49%|████▉ | 1220/2500 [00:06<00:06, 199.48it/s]
scan: 50%|████▉ | 1242/2500 [00:06<00:06, 203.61it/s]
scan: 51%|█████ | 1263/2500 [00:07<00:06, 202.06it/s]
scan: 51%|█████▏ | 1284/2500 [00:07<00:06, 199.60it/s]
scan: 52%|█████▏ | 1306/2500 [00:07<00:05, 205.03it/s]
scan: 53%|█████▎ | 1327/2500 [00:07<00:05, 200.03it/s]
scan: 54%|█████▍ | 1349/2500 [00:07<00:05, 205.52it/s]
scan: 55%|█████▍ | 1370/2500 [00:07<00:05, 193.72it/s]
scan: 56%|█████▌ | 1390/2500 [00:07<00:05, 192.99it/s]
scan: 56%|█████▋ | 1411/2500 [00:07<00:05, 195.38it/s]
scan: 57%|█████▋ | 1431/2500 [00:07<00:05, 190.76it/s]
scan: 58%|█████▊ | 1452/2500 [00:08<00:05, 194.84it/s]
scan: 59%|█████▉ | 1472/2500 [00:08<00:05, 187.11it/s]
scan: 60%|█████▉ | 1493/2500 [00:08<00:05, 191.74it/s]
scan: 61%|██████ | 1513/2500 [00:08<00:05, 193.27it/s]
scan: 61%|██████▏ | 1533/2500 [00:08<00:05, 187.12it/s]
scan: 62%|██████▏ | 1555/2500 [00:08<00:04, 193.92it/s]
scan: 63%|██████▎ | 1575/2500 [00:08<00:04, 188.87it/s]
scan: 64%|██████▍ | 1596/2500 [00:08<00:04, 194.70it/s]
scan: 65%|██████▍ | 1616/2500 [00:08<00:04, 192.85it/s]
scan: 65%|██████▌ | 1636/2500 [00:08<00:04, 191.75it/s]
scan: 66%|██████▋ | 1657/2500 [00:09<00:04, 193.99it/s]
scan: 67%|██████▋ | 1677/2500 [00:09<00:04, 184.64it/s]
scan: 68%|██████▊ | 1697/2500 [00:09<00:04, 188.80it/s]
scan: 69%|██████▊ | 1716/2500 [00:09<00:04, 188.01it/s]
scan: 69%|██████▉ | 1735/2500 [00:09<00:04, 158.29it/s]
scan: 70%|███████ | 1755/2500 [00:09<00:04, 168.44it/s]
scan: 71%|███████ | 1774/2500 [00:09<00:04, 172.75it/s]
scan: 72%|███████▏ | 1792/2500 [00:09<00:04, 172.55it/s]
scan: 72%|███████▏ | 1812/2500 [00:09<00:03, 177.81it/s]
scan: 73%|███████▎ | 1831/2500 [00:10<00:03, 180.06it/s]
scan: 74%|███████▍ | 1850/2500 [00:10<00:03, 178.22it/s]
scan: 75%|███████▍ | 1869/2500 [00:10<00:03, 179.53it/s]
scan: 76%|███████▌ | 1888/2500 [00:10<00:03, 176.20it/s]
scan: 76%|███████▋ | 1907/2500 [00:10<00:03, 180.09it/s]
scan: 77%|███████▋ | 1926/2500 [00:10<00:03, 181.70it/s]
scan: 78%|███████▊ | 1945/2500 [00:10<00:03, 176.53it/s]
scan: 79%|███████▊ | 1965/2500 [00:10<00:02, 180.87it/s]
scan: 79%|███████▉ | 1984/2500 [00:10<00:02, 182.19it/s]
scan: 80%|████████ | 2003/2500 [00:11<00:02, 177.93it/s]
scan: 81%|████████ | 2022/2500 [00:11<00:02, 179.42it/s]
scan: 82%|████████▏ | 2041/2500 [00:11<00:02, 179.96it/s]
scan: 82%|████████▏ | 2060/2500 [00:11<00:02, 173.68it/s]
scan: 83%|████████▎ | 2079/2500 [00:11<00:02, 176.46it/s]
scan: 84%|████████▍ | 2098/2500 [00:11<00:02, 177.59it/s]
scan: 85%|████████▍ | 2116/2500 [00:11<00:02, 176.58it/s]
scan: 85%|████████▌ | 2134/2500 [00:11<00:02, 177.54it/s]
scan: 86%|████████▌ | 2153/2500 [00:11<00:01, 179.09it/s]
scan: 87%|████████▋ | 2172/2500 [00:12<00:01, 180.49it/s]
scan: 88%|████████▊ | 2191/2500 [00:12<00:01, 180.59it/s]
scan: 88%|████████▊ | 2210/2500 [00:12<00:01, 180.94it/s]
scan: 89%|████████▉ | 2229/2500 [00:12<00:01, 180.20it/s]
scan: 90%|████████▉ | 2248/2500 [00:12<00:01, 179.02it/s]
scan: 91%|█████████ | 2267/2500 [00:12<00:01, 180.45it/s]
scan: 91%|█████████▏| 2286/2500 [00:12<00:01, 181.37it/s]
scan: 92%|█████████▏| 2305/2500 [00:12<00:01, 181.40it/s]
scan: 93%|█████████▎| 2324/2500 [00:12<00:00, 182.10it/s]
scan: 94%|█████████▎| 2343/2500 [00:12<00:00, 180.01it/s]
scan: 94%|█████████▍| 2362/2500 [00:13<00:00, 181.45it/s]
scan: 95%|█████████▌| 2381/2500 [00:13<00:00, 182.03it/s]
scan: 96%|█████████▌| 2400/2500 [00:13<00:00, 180.11it/s]
scan: 97%|█████████▋| 2419/2500 [00:13<00:00, 180.70it/s]
scan: 98%|█████████▊| 2438/2500 [00:13<00:00, 179.45it/s]
scan: 98%|█████████▊| 2457/2500 [00:13<00:00, 180.01it/s]
scan: 99%|█████████▉| 2476/2500 [00:13<00:00, 180.29it/s]
scan: 100%|█████████▉| 2495/2500 [00:13<00:00, 173.18it/s]
scan: 100%|██████████| 2500/2500 [00:13<00:00, 180.67it/s]
Passive harmonic cavity - Robinson instabilities with coupling and HOM
Fig. 3 (a) in [2]
result = scan_psi_I0(
SystemState(ring, [MC, HC]),
psi_HC_vals,
currents,
flow=flow,
eq_opts=eq_opts,
theory_opts=TheoryOptions(mode_coupling=True, f_HOM=1e9, Z_HOM=10e3))
result.plot("xi")
scan: 0%| | 0/2500 [00:00<?, ?it/s]
scan: 1%| | 18/2500 [00:00<00:13, 179.89it/s]
scan: 1%|▏ | 36/2500 [00:00<00:14, 170.52it/s]
scan: 2%|▏ | 54/2500 [00:00<00:18, 130.60it/s]
scan: 3%|▎ | 73/2500 [00:00<00:16, 149.03it/s]
scan: 4%|▎ | 89/2500 [00:00<00:16, 149.09it/s]
scan: 4%|▍ | 105/2500 [00:00<00:18, 132.37it/s]
scan: 5%|▌ | 125/2500 [00:00<00:15, 150.28it/s]
scan: 6%|▌ | 141/2500 [00:00<00:15, 147.95it/s]
scan: 6%|▋ | 157/2500 [00:01<00:17, 130.71it/s]
scan: 7%|▋ | 178/2500 [00:01<00:15, 149.90it/s]
scan: 8%|▊ | 194/2500 [00:01<00:15, 146.67it/s]
scan: 8%|▊ | 210/2500 [00:01<00:19, 120.01it/s]
scan: 9%|▉ | 232/2500 [00:01<00:15, 142.03it/s]
scan: 10%|▉ | 248/2500 [00:01<00:16, 133.23it/s]
scan: 11%|█ | 263/2500 [00:01<00:16, 134.86it/s]
scan: 11%|█▏ | 285/2500 [00:01<00:14, 156.47it/s]
scan: 12%|█▏ | 302/2500 [00:02<00:15, 145.74it/s]
scan: 13%|█▎ | 320/2500 [00:02<00:14, 153.98it/s]
scan: 14%|█▎ | 341/2500 [00:02<00:12, 167.38it/s]
scan: 14%|█▍ | 359/2500 [00:02<00:13, 158.61it/s]
scan: 15%|█▌ | 380/2500 [00:02<00:12, 172.27it/s]
scan: 16%|█▌ | 399/2500 [00:02<00:12, 171.94it/s]
scan: 17%|█▋ | 417/2500 [00:02<00:12, 172.67it/s]
scan: 18%|█▊ | 440/2500 [00:02<00:11, 179.67it/s]
scan: 18%|█▊ | 459/2500 [00:03<00:11, 170.97it/s]
scan: 19%|█▉ | 481/2500 [00:03<00:11, 182.05it/s]
scan: 20%|██ | 500/2500 [00:03<00:11, 177.89it/s]
scan: 21%|██ | 518/2500 [00:03<00:11, 176.51it/s]
scan: 22%|██▏ | 541/2500 [00:03<00:10, 189.70it/s]
scan: 22%|██▏ | 561/2500 [00:03<00:10, 182.94it/s]
scan: 23%|██▎ | 583/2500 [00:03<00:09, 192.06it/s]
scan: 24%|██▍ | 603/2500 [00:03<00:10, 183.33it/s]
scan: 25%|██▍ | 623/2500 [00:03<00:10, 187.65it/s]
scan: 26%|██▌ | 645/2500 [00:03<00:09, 190.07it/s]
scan: 27%|██▋ | 665/2500 [00:04<00:09, 186.98it/s]
scan: 27%|██▋ | 687/2500 [00:04<00:09, 195.92it/s]
scan: 28%|██▊ | 707/2500 [00:04<00:09, 196.77it/s]
scan: 29%|██▉ | 727/2500 [00:04<00:09, 196.06it/s]
scan: 30%|██▉ | 748/2500 [00:04<00:08, 194.98it/s]
scan: 31%|███ | 768/2500 [00:04<00:08, 192.86it/s]
scan: 32%|███▏ | 791/2500 [00:04<00:08, 201.08it/s]
scan: 32%|███▏ | 812/2500 [00:04<00:08, 191.22it/s]
scan: 33%|███▎ | 832/2500 [00:04<00:08, 191.76it/s]
scan: 34%|███▍ | 854/2500 [00:05<00:08, 197.63it/s]
scan: 35%|███▍ | 874/2500 [00:05<00:08, 195.76it/s]
scan: 36%|███▌ | 897/2500 [00:05<00:07, 203.71it/s]
scan: 37%|███▋ | 918/2500 [00:05<00:07, 198.70it/s]
scan: 38%|███▊ | 941/2500 [00:05<00:07, 205.11it/s]
scan: 38%|███▊ | 962/2500 [00:05<00:07, 203.31it/s]
scan: 39%|███▉ | 983/2500 [00:05<00:07, 201.04it/s]
scan: 40%|████ | 1005/2500 [00:05<00:07, 205.30it/s]
scan: 41%|████ | 1026/2500 [00:05<00:07, 201.92it/s]
scan: 42%|████▏ | 1049/2500 [00:05<00:06, 208.49it/s]
scan: 43%|████▎ | 1070/2500 [00:06<00:07, 197.20it/s]
scan: 44%|████▎ | 1092/2500 [00:06<00:06, 203.59it/s]
scan: 45%|████▍ | 1113/2500 [00:06<00:06, 201.68it/s]
scan: 45%|████▌ | 1134/2500 [00:06<00:06, 202.43it/s]
scan: 46%|████▌ | 1156/2500 [00:06<00:06, 204.48it/s]
scan: 47%|████▋ | 1177/2500 [00:06<00:06, 200.57it/s]
scan: 48%|████▊ | 1200/2500 [00:06<00:06, 209.00it/s]
scan: 49%|████▉ | 1221/2500 [00:06<00:06, 198.53it/s]
scan: 50%|████▉ | 1243/2500 [00:06<00:06, 204.20it/s]
scan: 51%|█████ | 1264/2500 [00:07<00:06, 202.27it/s]
scan: 51%|█████▏ | 1285/2500 [00:07<00:06, 200.24it/s]
scan: 52%|█████▏ | 1307/2500 [00:07<00:05, 204.21it/s]
scan: 53%|█████▎ | 1328/2500 [00:07<00:05, 199.66it/s]
scan: 54%|█████▍ | 1351/2500 [00:07<00:05, 206.58it/s]
scan: 55%|█████▍ | 1372/2500 [00:07<00:05, 194.67it/s]
scan: 56%|█████▌ | 1394/2500 [00:07<00:05, 199.83it/s]
scan: 57%|█████▋ | 1415/2500 [00:07<00:05, 199.15it/s]
scan: 57%|█████▋ | 1436/2500 [00:07<00:05, 195.48it/s]
scan: 58%|█████▊ | 1458/2500 [00:08<00:05, 199.72it/s]
scan: 59%|█████▉ | 1479/2500 [00:08<00:05, 192.63it/s]
scan: 60%|██████ | 1501/2500 [00:08<00:05, 198.53it/s]
scan: 61%|██████ | 1521/2500 [00:08<00:05, 188.77it/s]
scan: 62%|██████▏ | 1541/2500 [00:08<00:05, 187.54it/s]
scan: 62%|██████▏ | 1561/2500 [00:08<00:04, 190.43it/s]
scan: 63%|██████▎ | 1581/2500 [00:08<00:04, 187.47it/s]
scan: 64%|██████▍ | 1602/2500 [00:08<00:04, 193.56it/s]
scan: 65%|██████▍ | 1622/2500 [00:08<00:04, 191.10it/s]
scan: 66%|██████▌ | 1642/2500 [00:09<00:04, 192.01it/s]
scan: 66%|██████▋ | 1662/2500 [00:09<00:04, 193.08it/s]
scan: 67%|██████▋ | 1682/2500 [00:09<00:04, 184.50it/s]
scan: 68%|██████▊ | 1703/2500 [00:09<00:04, 189.81it/s]
scan: 69%|██████▉ | 1723/2500 [00:09<00:04, 186.79it/s]
scan: 70%|██████▉ | 1742/2500 [00:09<00:04, 180.98it/s]
scan: 70%|███████ | 1762/2500 [00:09<00:04, 183.52it/s]
scan: 71%|███████ | 1781/2500 [00:09<00:04, 178.89it/s]
scan: 72%|███████▏ | 1801/2500 [00:09<00:03, 183.97it/s]
scan: 73%|███████▎ | 1820/2500 [00:09<00:03, 184.07it/s]
scan: 74%|███████▎ | 1839/2500 [00:10<00:03, 178.62it/s]
scan: 74%|███████▍ | 1859/2500 [00:10<00:03, 182.09it/s]
scan: 75%|███████▌ | 1878/2500 [00:10<00:03, 182.13it/s]
scan: 76%|███████▌ | 1897/2500 [00:10<00:03, 179.18it/s]
scan: 77%|███████▋ | 1916/2500 [00:10<00:03, 181.26it/s]
scan: 77%|███████▋ | 1935/2500 [00:10<00:03, 181.66it/s]
scan: 78%|███████▊ | 1954/2500 [00:10<00:03, 175.79it/s]
scan: 79%|███████▉ | 1973/2500 [00:10<00:02, 178.59it/s]
scan: 80%|███████▉ | 1991/2500 [00:10<00:02, 178.17it/s]
scan: 80%|████████ | 2009/2500 [00:11<00:02, 176.18it/s]
scan: 81%|████████ | 2028/2500 [00:11<00:02, 178.21it/s]
scan: 82%|████████▏ | 2046/2500 [00:11<00:02, 177.68it/s]
scan: 83%|████████▎ | 2064/2500 [00:11<00:02, 173.63it/s]
scan: 83%|████████▎ | 2083/2500 [00:11<00:02, 176.02it/s]
scan: 84%|████████▍ | 2102/2500 [00:11<00:02, 177.84it/s]
scan: 85%|████████▍ | 2121/2500 [00:11<00:02, 179.15it/s]
scan: 86%|████████▌ | 2139/2500 [00:11<00:02, 178.55it/s]
scan: 86%|████████▋ | 2158/2500 [00:11<00:01, 179.09it/s]
scan: 87%|████████▋ | 2177/2500 [00:11<00:01, 179.93it/s]
scan: 88%|████████▊ | 2195/2500 [00:12<00:01, 179.45it/s]
scan: 89%|████████▊ | 2214/2500 [00:12<00:01, 180.39it/s]
scan: 89%|████████▉ | 2233/2500 [00:12<00:01, 179.70it/s]
scan: 90%|█████████ | 2252/2500 [00:12<00:01, 180.12it/s]
scan: 91%|█████████ | 2271/2500 [00:12<00:01, 180.28it/s]
scan: 92%|█████████▏| 2290/2500 [00:12<00:01, 179.40it/s]
scan: 92%|█████████▏| 2308/2500 [00:12<00:01, 178.81it/s]
scan: 93%|█████████▎| 2327/2500 [00:12<00:00, 180.10it/s]
scan: 94%|█████████▍| 2346/2500 [00:12<00:00, 178.32it/s]
scan: 95%|█████████▍| 2365/2500 [00:13<00:00, 179.60it/s]
scan: 95%|█████████▌| 2383/2500 [00:13<00:00, 179.31it/s]
scan: 96%|█████████▌| 2401/2500 [00:13<00:00, 179.22it/s]
scan: 97%|█████████▋| 2420/2500 [00:13<00:00, 180.33it/s]
scan: 98%|█████████▊| 2439/2500 [00:13<00:00, 179.50it/s]
scan: 98%|█████████▊| 2457/2500 [00:13<00:00, 179.40it/s]
scan: 99%|█████████▉| 2476/2500 [00:13<00:00, 181.17it/s]
scan: 100%|█████████▉| 2495/2500 [00:13<00:00, 181.92it/s]
scan: 100%|██████████| 2500/2500 [00:13<00:00, 181.35it/s]
