In [35]:
import math
import meep as mp
from meep import mpb
import numpy as np
import matplotlib
import matplotlib.pyplot as plt
from matplotlib.patches import Polygon, Patch
from matplotlib.ticker import FormatStrFormatter
In [8]:
num_bands = 2
resolution = 32
geometry_lattice = mp.Lattice(size=mp.Vector3(1, 1),
basis1=mp.Vector3(math.sqrt(3) / 2, 0.5),
basis2=mp.Vector3(math.sqrt(3) / 2, -0.5))
geometry = [mp.Cylinder(0.2, material=mp.Medium(epsilon=12.96))]
k_points = [
mp.Vector3(), # Gamma
mp.Vector3(y=0.5), # M
mp.Vector3(1 / -3, 1 / 3), # K
mp.Vector3(), # Gamma
]
k_points = mp.interpolate(4, k_points)
default_material = mp.Medium(epsilon=2.4)
ms = mpb.ModeSolver(
geometry=geometry,
geometry_lattice=geometry_lattice,
k_points=k_points,
resolution=resolution,
num_bands=num_bands,
default_material=default_material,
)
num_x = 30
results = []
radii = np.linspace(0., 0.5, num=num_x)
for r in radii:
geometry = [mp.Cylinder(r, material=mp.Medium(epsilon=12.96))]
ms.geometry = geometry
ms.run_tm()
if len(ms.gap_list) > 0:
results.append(ms.gap_list[0])
else: results.append((0, 0, 0))
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.0020432472229003906
Outputting epsilon.h5...
elapsed time for k point: 0.015569686889648438
elapsed time for k point: 0.014332056045532227
elapsed time for k point: 0.003877878189086914
elapsed time for k point: 0.0050661563873291016
elapsed time for k point: 0.002601146697998047
elapsed time for k point: 0.0053179264068603516
elapsed time for k point: 0.0023088455200195312
elapsed time for k point: 0.006564140319824219
elapsed time for k point: 0.00342559814453125
elapsed time for k point: 0.00859522819519043
elapsed time for k point: 0.0067310333251953125
elapsed time for k point: 0.002608060836791992
elapsed time for k point: 0.0051517486572265625
elapsed time for k point: 0.006944179534912109
elapsed time for k point: 0.0044879913330078125
elapsed time for k point: 0.003942728042602539
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.4303314829119348 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.37267799625183073 at Vector3<0.0, 0.5, 0.0> to 0.7453559925839891 at Vector3<0.0, 0.0, 0.0>
total elapsed time for run: 0.12585997581481934
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.004997968673706055
Outputting epsilon.h5...
elapsed time for k point: 0.015418767929077148
elapsed time for k point: 0.021904945373535156
elapsed time for k point: 0.013065814971923828
elapsed time for k point: 0.006699323654174805
elapsed time for k point: 0.006916999816894531
elapsed time for k point: 0.008975744247436523
elapsed time for k point: 0.005651950836181641
elapsed time for k point: 0.008779048919677734
elapsed time for k point: 0.005898714065551758
elapsed time for k point: 0.02965521812438965
elapsed time for k point: 0.011137008666992188
elapsed time for k point: 0.011982202529907227
elapsed time for k point: 0.00639796257019043
elapsed time for k point: 0.0065724849700927734
elapsed time for k point: 0.011699914932250977
elapsed time for k point: 0.009176492691040039
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.427448599944701 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.37267602613260614 at Vector3<0.0, 0.5, 0.0> to 0.735111350292789 at Vector3<0.0, 0.0, 0.0>
total elapsed time for run: 0.1995527744293213
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.0038139820098876953
Outputting epsilon.h5...
elapsed time for k point: 0.01292562484741211
elapsed time for k point: 0.018057584762573242
elapsed time for k point: 0.009884119033813477
elapsed time for k point: 0.008152246475219727
elapsed time for k point: 0.008969783782958984
elapsed time for k point: 0.006822347640991211
elapsed time for k point: 0.005202531814575195
elapsed time for k point: 0.005193948745727539
elapsed time for k point: 0.00649261474609375
elapsed time for k point: 0.011062145233154297
elapsed time for k point: 0.015573501586914062
elapsed time for k point: 0.013052940368652344
elapsed time for k point: 0.007429361343383789
elapsed time for k point: 0.010134220123291016
elapsed time for k point: 0.027596473693847656
elapsed time for k point: 0.010506391525268555
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.4176806886314358 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.37264271944916555 at Vector3<0.0, 0.5, 0.0> to 0.7001210807167477 at Vector3<0.0, 0.0, 0.0>
total elapsed time for run: 0.20031118392944336
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.0035233497619628906
Outputting epsilon.h5...
elapsed time for k point: 0.010747194290161133
elapsed time for k point: 0.011946439743041992
elapsed time for k point: 0.01031041145324707
elapsed time for k point: 0.010624408721923828
elapsed time for k point: 0.006540536880493164
elapsed time for k point: 0.0073740482330322266
elapsed time for k point: 0.005744457244873047
elapsed time for k point: 0.009068489074707031
elapsed time for k point: 0.00745391845703125
elapsed time for k point: 0.0068280696868896484
elapsed time for k point: 0.02046036720275879
elapsed time for k point: 0.011612415313720703
elapsed time for k point: 0.008424520492553711
elapsed time for k point: 0.01041865348815918
elapsed time for k point: 0.010540485382080078
elapsed time for k point: 0.006359100341796875
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.40190834881427157 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.37251944683613625 at Vector3<0.0, 0.5, 0.0> to 0.6491545075188152 at Vector3<0.0, 0.0, 0.0>
total elapsed time for run: 0.17689776420593262
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.007841825485229492
Outputting epsilon.h5...
elapsed time for k point: 0.016054630279541016
elapsed time for k point: 0.013176918029785156
elapsed time for k point: 0.010593414306640625
elapsed time for k point: 0.008335590362548828
elapsed time for k point: 0.009932279586791992
elapsed time for k point: 0.0067577362060546875
elapsed time for k point: 0.009006261825561523
elapsed time for k point: 0.006498098373413086
elapsed time for k point: 0.008929252624511719
elapsed time for k point: 0.007028818130493164
elapsed time for k point: 0.012192010879516602
elapsed time for k point: 0.010452032089233398
elapsed time for k point: 0.007930517196655273
elapsed time for k point: 0.010270118713378906
elapsed time for k point: 0.010852575302124023
elapsed time for k point: 0.005564451217651367
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.3815122467391655 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.372210800110519 at Vector3<0.0, 0.5, 0.0> to 0.5971987244881549 at Vector3<0.0, 0.0, 0.0>
total elapsed time for run: 0.18337011337280273
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.04052591323852539
Outputting epsilon.h5...
elapsed time for k point: 0.007176876068115234
elapsed time for k point: 0.02064990997314453
elapsed time for k point: 0.012527704238891602
elapsed time for k point: 0.010232686996459961
elapsed time for k point: 0.007762908935546875
elapsed time for k point: 0.008159399032592773
elapsed time for k point: 0.006475687026977539
elapsed time for k point: 0.010475873947143555
elapsed time for k point: 0.007549762725830078
elapsed time for k point: 0.006558418273925781
elapsed time for k point: 0.021432876586914062
elapsed time for k point: 0.015729427337646484
elapsed time for k point: 0.007999897003173828
elapsed time for k point: 0.009196996688842773
elapsed time for k point: 0.007917404174804688
elapsed time for k point: 0.007232189178466797
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.3596300570268307 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.37149461379122356 at Vector3<0.0, 0.5, 0.0> to 0.5561172373107958 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.3596300570268307) to band 2 (0.37149461379122356), 3.2455632364635205%
total elapsed time for run: 0.22519445419311523
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.005449771881103516
Outputting epsilon.h5...
elapsed time for k point: 0.01604175567626953
elapsed time for k point: 0.012882232666015625
elapsed time for k point: 0.009740829467773438
elapsed time for k point: 0.009322881698608398
elapsed time for k point: 0.008375883102416992
elapsed time for k point: 0.00927734375
elapsed time for k point: 0.006365776062011719
elapsed time for k point: 0.008934259414672852
elapsed time for k point: 0.007727384567260742
elapsed time for k point: 0.008385896682739258
elapsed time for k point: 0.006421804428100586
elapsed time for k point: 0.009249448776245117
elapsed time for k point: 0.00833582878112793
elapsed time for k point: 0.009177684783935547
elapsed time for k point: 0.007880210876464844
elapsed time for k point: 0.004593372344970703
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.3379096666906223 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.37023016361300026 at Vector3<0.0, 0.5, 0.0> to 0.5238236774850537 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.3379096666906223) to band 2 (0.37023016361300026), 9.12828103695854%
total elapsed time for run: 0.17769670486450195
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.007081270217895508
Outputting epsilon.h5...
elapsed time for k point: 0.009630918502807617
elapsed time for k point: 0.013149499893188477
elapsed time for k point: 0.010691642761230469
elapsed time for k point: 0.007991790771484375
elapsed time for k point: 0.009687423706054688
elapsed time for k point: 0.030165672302246094
elapsed time for k point: 0.006558895111083984
elapsed time for k point: 0.008785486221313477
elapsed time for k point: 0.006561756134033203
elapsed time for k point: 0.008621931076049805
elapsed time for k point: 0.01914191246032715
elapsed time for k point: 0.018759727478027344
elapsed time for k point: 0.006554365158081055
elapsed time for k point: 0.007618427276611328
elapsed time for k point: 0.011317253112792969
elapsed time for k point: 0.004635810852050781
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.31798968116492854 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.3680958277292817 at Vector3<0.0, 0.5, 0.0> to 0.500701505386161 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.31798968116492854) to band 2 (0.3680958277292817), 14.606385330922118%
total elapsed time for run: 0.2051105499267578
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.005789756774902344
Outputting epsilon.h5...
elapsed time for k point: 0.008813619613647461
elapsed time for k point: 0.009111642837524414
elapsed time for k point: 0.011197090148925781
elapsed time for k point: 0.010294437408447266
elapsed time for k point: 0.0076558589935302734
elapsed time for k point: 0.012696504592895508
elapsed time for k point: 0.024486303329467773
elapsed time for k point: 0.008550405502319336
elapsed time for k point: 0.008898496627807617
elapsed time for k point: 0.0065495967864990234
elapsed time for k point: 0.014751195907592773
elapsed time for k point: 0.10295915603637695
elapsed time for k point: 0.007704019546508789
elapsed time for k point: 0.010135650634765625
elapsed time for k point: 0.011961936950683594
elapsed time for k point: 0.02894306182861328
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.3003724295839907 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.3647495353383683 at Vector3<0.0, 0.5, 0.0> to 0.4847635625715469 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.3003724295839907) to band 2 (0.3647495353383683), 19.3579851965624%
total elapsed time for run: 0.31070709228515625
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.00590968132019043
Outputting epsilon.h5...
elapsed time for k point: 0.009729146957397461
elapsed time for k point: 0.00916147232055664
elapsed time for k point: 0.00933384895324707
elapsed time for k point: 0.008167505264282227
elapsed time for k point: 0.00969076156616211
elapsed time for k point: 0.006464719772338867
elapsed time for k point: 0.006502866744995117
elapsed time for k point: 0.008850812911987305
elapsed time for k point: 0.008008956909179688
elapsed time for k point: 0.006469249725341797
elapsed time for k point: 0.016292810440063477
elapsed time for k point: 0.07082390785217285
elapsed time for k point: 0.011766433715820312
elapsed time for k point: 0.029861927032470703
elapsed time for k point: 0.009109020233154297
elapsed time for k point: 0.0075724124908447266
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.2847289911277431 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.3599829833929977 at Vector3<0.0, 0.5, 0.0> to 0.47360651450710417 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.2847289911277431) to band 2 (0.3599829833929977), 23.344995979389452%
total elapsed time for run: 0.25422096252441406
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.005955219268798828
Outputting epsilon.h5...
elapsed time for k point: 0.010213375091552734
elapsed time for k point: 0.00901341438293457
elapsed time for k point: 0.009730815887451172
elapsed time for k point: 0.009781360626220703
elapsed time for k point: 0.00766444206237793
elapsed time for k point: 0.0077669620513916016
elapsed time for k point: 0.00676274299621582
elapsed time for k point: 0.00953984260559082
elapsed time for k point: 0.006562948226928711
elapsed time for k point: 0.0065190792083740234
elapsed time for k point: 0.015818357467651367
elapsed time for k point: 0.05258297920227051
elapsed time for k point: 0.012832880020141602
elapsed time for k point: 0.031458139419555664
elapsed time for k point: 0.011996030807495117
elapsed time for k point: 0.006952524185180664
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.2712273414430768 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.3534730556941499 at Vector3<0.0, 0.5, 0.0> to 0.46691129158187933 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.2712273414430768) to band 2 (0.3534730556941499), 26.33125082934958%
total elapsed time for run: 0.2398388385772705
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.00758671760559082
Outputting epsilon.h5...
elapsed time for k point: 0.008674383163452148
elapsed time for k point: 0.013330936431884766
elapsed time for k point: 0.011070489883422852
elapsed time for k point: 0.010954141616821289
elapsed time for k point: 0.009789466857910156
elapsed time for k point: 0.00755000114440918
elapsed time for k point: 0.00683903694152832
elapsed time for k point: 0.010393142700195312
elapsed time for k point: 0.00658726692199707
elapsed time for k point: 0.008565664291381836
elapsed time for k point: 0.019547700881958008
elapsed time for k point: 0.11232757568359375
elapsed time for k point: 0.007628202438354492
elapsed time for k point: 0.011955022811889648
elapsed time for k point: 0.011322259902954102
elapsed time for k point: 0.00756072998046875
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.25944201761549357 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.34524243602972143 at Vector3<0.0, 0.5, 0.0> to 0.4630650177553006 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.25944201761549357) to band 2 (0.34524243602972143), 28.378575932289255%
total elapsed time for run: 0.29148244857788086
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.006531476974487305
Outputting epsilon.h5...
elapsed time for k point: 0.011752843856811523
elapsed time for k point: 0.011745691299438477
elapsed time for k point: 0.011689186096191406
elapsed time for k point: 0.010757684707641602
elapsed time for k point: 0.03177452087402344
elapsed time for k point: 0.0064220428466796875
elapsed time for k point: 0.011904716491699219
elapsed time for k point: 0.006444692611694336
elapsed time for k point: 0.007991552352905273
elapsed time for k point: 0.006451845169067383
elapsed time for k point: 0.017339229583740234
elapsed time for k point: 0.020989418029785156
elapsed time for k point: 0.0077440738677978516
elapsed time for k point: 0.009333610534667969
elapsed time for k point: 0.011401176452636719
elapsed time for k point: 0.01043248176574707
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.24909349130006206 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.33546105344525623 at Vector3<0.0, 0.5, 0.0> to 0.46125678325535685 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.24909349130006206) to band 2 (0.33546105344525623), 29.549872778022188%
total elapsed time for run: 0.21841979026794434
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.006104946136474609
Outputting epsilon.h5...
elapsed time for k point: 0.012239456176757812
elapsed time for k point: 0.015768051147460938
elapsed time for k point: 0.03360247611999512
elapsed time for k point: 0.007421970367431641
elapsed time for k point: 0.009049177169799805
elapsed time for k point: 0.00762939453125
elapsed time for k point: 0.008574962615966797
elapsed time for k point: 0.006119966506958008
elapsed time for k point: 0.007734775543212891
elapsed time for k point: 0.006211996078491211
elapsed time for k point: 0.014813899993896484
elapsed time for k point: 0.07193541526794434
elapsed time for k point: 0.0073888301849365234
elapsed time for k point: 0.012048959732055664
elapsed time for k point: 0.011528491973876953
elapsed time for k point: 0.012755632400512695
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.2401807280848055 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.3244438339920495 at Vector3<0.0, 0.5, 0.0> to 0.460803673259308 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.2401807280848055) to band 2 (0.3244438339920495), 29.847481518444585%
total elapsed time for run: 0.2685680389404297
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.005708456039428711
Outputting epsilon.h5...
elapsed time for k point: 0.018053054809570312
elapsed time for k point: 0.02498340606689453
elapsed time for k point: 0.009145259857177734
elapsed time for k point: 0.007386684417724609
elapsed time for k point: 0.009934663772583008
elapsed time for k point: 0.0074160099029541016
elapsed time for k point: 0.008053302764892578
elapsed time for k point: 0.007385969161987305
elapsed time for k point: 0.010335445404052734
elapsed time for k point: 0.006270885467529297
elapsed time for k point: 0.01377415657043457
elapsed time for k point: 0.04055523872375488
elapsed time for k point: 0.012190580368041992
elapsed time for k point: 0.0313563346862793
elapsed time for k point: 0.009708166122436523
elapsed time for k point: 0.012166023254394531
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.23232031812564818 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.3126702545361094 at Vector3<0.0, 0.5, 0.0> to 0.44727605648060653 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.23232031812564818) to band 2 (0.3126702545361094), 29.486725253990596%
total elapsed time for run: 0.27124643325805664
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.007752418518066406
Outputting epsilon.h5...
elapsed time for k point: 0.019285917282104492
elapsed time for k point: 0.022607088088989258
elapsed time for k point: 0.010709285736083984
elapsed time for k point: 0.007485389709472656
elapsed time for k point: 0.008809566497802734
elapsed time for k point: 0.006140947341918945
elapsed time for k point: 0.00793147087097168
elapsed time for k point: 0.00617527961730957
elapsed time for k point: 0.010277271270751953
elapsed time for k point: 0.006182193756103516
elapsed time for k point: 0.007628917694091797
elapsed time for k point: 0.03479456901550293
elapsed time for k point: 0.009055137634277344
elapsed time for k point: 0.009550094604492188
elapsed time for k point: 0.009613990783691406
elapsed time for k point: 0.009310245513916016
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.22541712304426112 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.30052497720124605 at Vector3<0.0, 0.5, 0.0> to 0.42675856637976556 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.22541712304426112) to band 2 (0.30052497720124605), 28.56126334892185%
total elapsed time for run: 0.20967435836791992
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.009538888931274414
Outputting epsilon.h5...
elapsed time for k point: 0.011218547821044922
elapsed time for k point: 0.023782014846801758
elapsed time for k point: 0.010711431503295898
elapsed time for k point: 0.009190559387207031
elapsed time for k point: 0.007518291473388672
elapsed time for k point: 0.008018970489501953
elapsed time for k point: 0.006293296813964844
elapsed time for k point: 0.006314516067504883
elapsed time for k point: 0.008066654205322266
elapsed time for k point: 0.012271404266357422
elapsed time for k point: 0.01322793960571289
elapsed time for k point: 0.019582271575927734
elapsed time for k point: 0.007770538330078125
elapsed time for k point: 0.032518625259399414
elapsed time for k point: 0.01071023941040039
elapsed time for k point: 0.010648250579833984
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.21944963816319563 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.28838962180689576 at Vector3<0.0, 0.5, 0.0> to 0.40936101489946863 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.21944963816319563) to band 2 (0.28838962180689576), 27.150316676091673%
total elapsed time for run: 0.2211611270904541
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.00680851936340332
Outputting epsilon.h5...
elapsed time for k point: 0.016564369201660156
elapsed time for k point: 0.020645618438720703
elapsed time for k point: 0.008875370025634766
elapsed time for k point: 0.00921773910522461
elapsed time for k point: 0.00978541374206543
elapsed time for k point: 0.006398916244506836
elapsed time for k point: 0.006354093551635742
elapsed time for k point: 0.007939815521240234
elapsed time for k point: 0.008458852767944336
elapsed time for k point: 0.009562492370605469
elapsed time for k point: 0.014513731002807617
elapsed time for k point: 0.019338607788085938
elapsed time for k point: 0.006505489349365234
elapsed time for k point: 0.029655933380126953
elapsed time for k point: 0.012897491455078125
elapsed time for k point: 0.009743213653564453
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.21415770433065823 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.27646295382370867 at Vector3<0.0, 0.5, 0.0> to 0.3938464210817432 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.21415770433065823) to band 2 (0.27646295382370867), 25.39854303217985%
total elapsed time for run: 0.2198042869567871
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.005304813385009766
Outputting epsilon.h5...
elapsed time for k point: 0.021590709686279297
elapsed time for k point: 0.02105259895324707
elapsed time for k point: 0.012128114700317383
elapsed time for k point: 0.009975910186767578
elapsed time for k point: 0.007814645767211914
elapsed time for k point: 0.008364677429199219
elapsed time for k point: 0.006476879119873047
elapsed time for k point: 0.010046720504760742
elapsed time for k point: 0.006679534912109375
elapsed time for k point: 0.009398698806762695
elapsed time for k point: 0.008862733840942383
elapsed time for k point: 0.013497352600097656
elapsed time for k point: 0.006445646286010742
elapsed time for k point: 0.007595539093017578
elapsed time for k point: 0.03910565376281738
elapsed time for k point: 0.01199650764465332
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.20954095816741916 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.26491697657189617 at Vector3<0.0, 0.5, 0.0> to 0.3804368593094669 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.20954095816741916) to band 2 (0.26491697657189617), 23.342856910972575%
total elapsed time for run: 0.2225337028503418
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.0050318241119384766
Outputting epsilon.h5...
elapsed time for k point: 0.011129140853881836
elapsed time for k point: 0.024438142776489258
elapsed time for k point: 0.009493827819824219
elapsed time for k point: 0.007880687713623047
elapsed time for k point: 0.009636878967285156
elapsed time for k point: 0.00914144515991211
elapsed time for k point: 0.006769895553588867
elapsed time for k point: 0.009268522262573242
elapsed time for k point: 0.0077953338623046875
elapsed time for k point: 0.00928807258605957
elapsed time for k point: 0.01857304573059082
elapsed time for k point: 0.013285160064697266
elapsed time for k point: 0.0064203739166259766
elapsed time for k point: 0.012028694152832031
elapsed time for k point: 0.030123472213745117
elapsed time for k point: 0.008208990097045898
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.2055368482140373 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.25386742097353143 at Vector3<0.0, 0.5, 0.0> to 0.36899646222556076 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.2055368482140373) to band 2 (0.25386742097353143), 21.040541414629903%
total elapsed time for run: 0.21728515625
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.00667262077331543
Outputting epsilon.h5...
elapsed time for k point: 0.014522790908813477
elapsed time for k point: 0.024562835693359375
elapsed time for k point: 0.008951425552368164
elapsed time for k point: 0.009439229965209961
elapsed time for k point: 0.0075719356536865234
elapsed time for k point: 0.007691621780395508
elapsed time for k point: 0.006344795227050781
elapsed time for k point: 0.008528709411621094
elapsed time for k point: 0.0074138641357421875
elapsed time for k point: 0.008636236190795898
elapsed time for k point: 0.007858514785766602
elapsed time for k point: 0.014900922775268555
elapsed time for k point: 0.006283760070800781
elapsed time for k point: 0.007615566253662109
elapsed time for k point: 0.01648712158203125
elapsed time for k point: 0.02849555015563965
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.20200907494544726 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.2433029553594587 at Vector3<0.0, 0.5, 0.0> to 0.35891010374895965 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.20200907494544726) to band 2 (0.2433029553594587), 18.54604304569874%
total elapsed time for run: 0.21025347709655762
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.008257627487182617
Outputting epsilon.h5...
elapsed time for k point: 0.0175173282623291
elapsed time for k point: 0.02387547492980957
elapsed time for k point: 0.009181499481201172
elapsed time for k point: 0.007776975631713867
elapsed time for k point: 0.007715463638305664
elapsed time for k point: 0.006155490875244141
elapsed time for k point: 0.007084369659423828
elapsed time for k point: 0.006151676177978516
elapsed time for k point: 0.008712291717529297
elapsed time for k point: 0.0075342655181884766
elapsed time for k point: 0.008493423461914062
elapsed time for k point: 0.015182018280029297
elapsed time for k point: 0.006308555603027344
elapsed time for k point: 0.03551983833312988
elapsed time for k point: 0.011893510818481445
elapsed time for k point: 0.006131649017333984
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.19897955569052891 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.23328198085300075 at Vector3<0.0, 0.5, 0.0> to 0.3504655621018816 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.19897955569052891) to band 2 (0.23328198085300075), 15.871143862006566%
total elapsed time for run: 0.2124171257019043
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.009406328201293945
Outputting epsilon.h5...
elapsed time for k point: 0.011652946472167969
elapsed time for k point: 0.02480792999267578
elapsed time for k point: 0.00901174545288086
elapsed time for k point: 0.007868289947509766
elapsed time for k point: 0.007724285125732422
elapsed time for k point: 0.008044958114624023
elapsed time for k point: 0.006327629089355469
elapsed time for k point: 0.007294416427612305
elapsed time for k point: 0.006405830383300781
elapsed time for k point: 0.011241674423217773
elapsed time for k point: 0.016658306121826172
elapsed time for k point: 0.020248889923095703
elapsed time for k point: 0.010139942169189453
elapsed time for k point: 0.03492093086242676
elapsed time for k point: 0.013141393661499023
elapsed time for k point: 0.0060901641845703125
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.1963590677066445 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.22377934628140955 at Vector3<0.0, 0.5, 0.0> to 0.34335905571039826 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.1963590677066445) to band 2 (0.22377934628140955), 13.052973811409071%
total elapsed time for run: 0.22693848609924316
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.008067131042480469
Outputting epsilon.h5...
elapsed time for k point: 0.01746678352355957
elapsed time for k point: 0.021679162979125977
elapsed time for k point: 0.010233879089355469
elapsed time for k point: 0.00881338119506836
elapsed time for k point: 0.008985042572021484
elapsed time for k point: 0.006472349166870117
elapsed time for k point: 0.006815910339355469
elapsed time for k point: 0.0073087215423583984
elapsed time for k point: 0.006671428680419922
elapsed time for k point: 0.00834965705871582
elapsed time for k point: 0.007554769515991211
elapsed time for k point: 0.04231142997741699
elapsed time for k point: 0.006207942962646484
elapsed time for k point: 0.010961294174194336
elapsed time for k point: 0.012956857681274414
elapsed time for k point: 0.007831811904907227
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.19408930325830304 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.21475388766510534 at Vector3<0.0, 0.5, 0.0> to 0.33727998828840955 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.19408930325830304) to band 2 (0.21475388766510534), 10.108806929194303%
total elapsed time for run: 0.22093582153320312
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.009124040603637695
Outputting epsilon.h5...
elapsed time for k point: 0.010072708129882812
elapsed time for k point: 0.027571678161621094
elapsed time for k point: 0.053467512130737305
elapsed time for k point: 0.00766444206237793
elapsed time for k point: 0.014438867568969727
elapsed time for k point: 0.032671451568603516
elapsed time for k point: 0.008462190628051758
elapsed time for k point: 0.0072329044342041016
elapsed time for k point: 0.00627589225769043
elapsed time for k point: 0.007422208786010742
elapsed time for k point: 0.008940458297729492
elapsed time for k point: 0.014473915100097656
elapsed time for k point: 0.0064144134521484375
elapsed time for k point: 0.009017467498779297
elapsed time for k point: 0.009199857711791992
elapsed time for k point: 0.008310794830322266
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.1921797142329697 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.206225454123271 at Vector3<0.0, 0.5, 0.0> to 0.33257400243480123 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.1921797142329697) to band 2 (0.206225454123271), 7.050982771258665%
total elapsed time for run: 0.25934576988220215
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.007510185241699219
Outputting epsilon.h5...
elapsed time for k point: 0.017840862274169922
elapsed time for k point: 0.021982908248901367
elapsed time for k point: 0.011595487594604492
elapsed time for k point: 0.03897833824157715
elapsed time for k point: 0.00899195671081543
elapsed time for k point: 0.009599685668945312
elapsed time for k point: 0.0066111087799072266
elapsed time for k point: 0.007999658584594727
elapsed time for k point: 0.008045434951782227
elapsed time for k point: 0.009558916091918945
elapsed time for k point: 0.008394956588745117
elapsed time for k point: 0.014480352401733398
elapsed time for k point: 0.008033514022827148
elapsed time for k point: 0.008238792419433594
elapsed time for k point: 0.013323545455932617
elapsed time for k point: 0.010603666305541992
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.19053604717673342 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.19812591186529888 at Vector3<0.0, 0.5, 0.0> to 0.32874004371180704 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.19053604717673342) to band 2 (0.19812591186529888), 3.905638054865381%
total elapsed time for run: 0.22994041442871094
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.008875131607055664
Outputting epsilon.h5...
elapsed time for k point: 0.02390766143798828
elapsed time for k point: 0.05195450782775879
elapsed time for k point: 0.009643077850341797
elapsed time for k point: 0.01063847541809082
elapsed time for k point: 0.008606672286987305
elapsed time for k point: 0.00990915298461914
elapsed time for k point: 0.007128000259399414
elapsed time for k point: 0.023891210556030273
elapsed time for k point: 0.016405582427978516
elapsed time for k point: 0.014245033264160156
elapsed time for k point: 0.008666753768920898
elapsed time for k point: 0.023525714874267578
elapsed time for k point: 0.006838321685791016
elapsed time for k point: 0.011743307113647461
elapsed time for k point: 0.012550115585327148
elapsed time for k point: 0.008023738861083984
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.1891282454768944 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.19043026441380492 at Vector3<0.0, 0.5, 0.0> to 0.32585286662626767 at Vector3<0.0, 0.0, 0.0>
Gap from band 1 (0.1891282454768944) to band 2 (0.19043026441380492), 0.6860702120921748%
total elapsed time for run: 0.2761070728302002
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.008478403091430664
Outputting epsilon.h5...
elapsed time for k point: 0.01838207244873047
elapsed time for k point: 0.028977394104003906
elapsed time for k point: 0.009171485900878906
elapsed time for k point: 0.007512807846069336
elapsed time for k point: 0.008861541748046875
elapsed time for k point: 0.007325410842895508
elapsed time for k point: 0.006292581558227539
elapsed time for k point: 0.007428169250488281
elapsed time for k point: 0.008507966995239258
elapsed time for k point: 0.00979161262512207
elapsed time for k point: 0.025058269500732422
elapsed time for k point: 0.01784515380859375
elapsed time for k point: 0.008170127868652344
elapsed time for k point: 0.007470130920410156
elapsed time for k point: 0.014076709747314453
elapsed time for k point: 0.011316537857055664
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.18793424646219037 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.1831404452139093 at Vector3<0.0, 0.5, 0.0> to 0.32390175111875436 at Vector3<0.0, 0.0, 0.0>
total elapsed time for run: 0.22354531288146973
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.01015782356262207
Outputting epsilon.h5...
elapsed time for k point: 0.03270673751831055
elapsed time for k point: 0.03740096092224121
elapsed time for k point: 0.00975346565246582
elapsed time for k point: 0.007510662078857422
elapsed time for k point: 0.010356664657592773
elapsed time for k point: 0.00794219970703125
elapsed time for k point: 0.01813507080078125
elapsed time for k point: 0.006257772445678711
elapsed time for k point: 0.007811307907104492
elapsed time for k point: 0.006304025650024414
elapsed time for k point: 0.011698246002197266
elapsed time for k point: 0.010148048400878906
elapsed time for k point: 0.007124900817871094
elapsed time for k point: 0.012431144714355469
elapsed time for k point: 0.04539036750793457
elapsed time for k point: 0.009196043014526367
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.18686628936045396 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.17620837910120452 at Vector3<0.0, 0.5, 0.0> to 0.32253158215588007 at Vector3<0.0, 0.0, 0.0>
total elapsed time for run: 0.26797938346862793
done
Initializing eigensolver data
Computing 2 bands with 1e-07 tolerance
16 k-points
Vector3<0.0, 0.0, 0.0>
Vector3<0.0, 0.1, 0.0>
Vector3<0.0, 0.2, 0.0>
Vector3<0.0, 0.30000000000000004, 0.0>
Vector3<0.0, 0.4, 0.0>
Vector3<0.0, 0.5, 0.0>
Vector3<-0.06666666666666667, 0.4666666666666667, 0.0>
Vector3<-0.13333333333333333, 0.43333333333333335, 0.0>
Vector3<-0.2, 0.39999999999999997, 0.0>
Vector3<-0.26666666666666666, 0.36666666666666664, 0.0>
Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Vector3<-0.26666666666666666, 0.26666666666666666, 0.0>
Vector3<-0.19999999999999998, 0.19999999999999998, 0.0>
Vector3<-0.1333333333333333, 0.1333333333333333, 0.0>
Vector3<-0.06666666666666665, 0.06666666666666665, 0.0>
Vector3<0.0, 0.0, 0.0>
elapsed time for initialization: 0.006999492645263672
Outputting epsilon.h5...
elapsed time for k point: 0.018413305282592773
elapsed time for k point: 0.0267486572265625
elapsed time for k point: 0.01043558120727539
elapsed time for k point: 0.007714271545410156
elapsed time for k point: 0.006364583969116211
elapsed time for k point: 0.008165836334228516
elapsed time for k point: 0.005194664001464844
elapsed time for k point: 0.0092926025390625
elapsed time for k point: 0.0062291622161865234
elapsed time for k point: 0.008574724197387695
elapsed time for k point: 0.021870136260986328
elapsed time for k point: 0.022403717041015625
elapsed time for k point: 0.007881641387939453
elapsed time for k point: 0.009210586547851562
elapsed time for k point: 0.010993719100952148
elapsed time for k point: 0.010324954986572266
Band 1 range: 0.0 at Vector3<0.0, 0.0, 0.0> to 0.18587838633269957 at Vector3<-0.3333333333333333, 0.3333333333333333, 0.0>
Band 2 range: 0.16963078210290733 at Vector3<0.0, 0.5, 0.0> to 0.3216684288657836 at Vector3<0.0, 0.0, 0.0>
total elapsed time for run: 0.21781206130981445
done
In [9]:
radii
Out[9]:
array([0. , 0.01724138, 0.03448276, 0.05172414, 0.06896552,
0.0862069 , 0.10344828, 0.12068966, 0.13793103, 0.15517241,
0.17241379, 0.18965517, 0.20689655, 0.22413793, 0.24137931,
0.25862069, 0.27586207, 0.29310345, 0.31034483, 0.32758621,
0.34482759, 0.36206897, 0.37931034, 0.39655172, 0.4137931 ,
0.43103448, 0.44827586, 0.46551724, 0.48275862, 0.5 ])
In [10]:
polypoints = [(r, b[2]) for r, b in zip(radii, results) if b != (0, 0, 0)]
polypoints.extend([(r, b[1]) for r, b in zip(radii, results) if b != (0, 0, 0)][::-1])
In [11]:
max_idx = np.argmax([b[0] for b in results])
max_gap_centre = np.mean(results[max_idx][1:])
max_gap_half = results[max_idx][2] - max_gap_centre
In [12]:
print(radii[np.argmax([b[0] for b in results])])
0.22413793103448276
In [37]:
matplotlib.rcParams['axes.labelsize'] = 18
matplotlib.rcParams['xtick.labelsize'] = 14
matplotlib.rcParams['ytick.labelsize'] = 14
matplotlib.rcParams['legend.fontsize'] = 14
fig, ax = plt.subplots(constrained_layout=True)
poly = Polygon(polypoints)
ax.errorbar(radii[max_idx], max_gap_centre, yerr=max_gap_half, color='black')
ax.plot([radii[max_idx], 0.5], [results[max_idx][2]]*2, color='black', ls=':', alpha=0.5)
ax.plot([radii[max_idx], 0.5], [results[max_idx][1]]*2, color='black', ls=':', alpha=0.5)
ax.axvline(radii[max_idx], ls='--', color='black', alpha=0.5, zorder=3)
ax.add_patch(poly)
ax.legend(handles=[Patch(facecolor='C0', label='first photonic band gap')], loc='upper right')
ax.set_xlim([0, 0.5])
ax.set_ylim([0.15, 0.4])
axy = ax.twinx()
axy.set_xlim([0, 0.5])
axy.set_ylim([0.15, 0.4])
axy.set_yticks([results[max_idx][2], results[max_idx][1]])
axy.yaxis.set_major_formatter(FormatStrFormatter('%.3f'))
ax.set_xlabel(r'pillar radius / $r/a$')
ax.set_ylabel(r'frequency / $\omega a/2\pi c$')
fig.savefig('bandgap_vs_ratio.pdf')
In [2]:
print('TM gaps')
for gap in tm_gaps:
print(gap)
print('TE gaps')
for gap in te_gaps:
print(gap)
TM gaps
(47.39883483253789, 0.2746902260167482, 0.4453310810690987)
(5.6740639480231065, 0.5605607990535272, 0.5932960878698785)
(0.014223870207005853, 0.7907195975678281, 0.7908320764965159)
(0.017161966966868752, 0.8374511179880573, 0.8375948534062025)
(0.22277612770986174, 0.8672009266363844, 0.8691349976044138)
TE gaps
(0.007916179285617626, 0.5631059454933676, 0.5631505237340275)
(5.190482251642612, 0.8214230165372104, 0.8651948161669857)
Content source: alkamid/lab-scripts
Similar notebooks: