Skip to content
Snippets Groups Projects
Commit 8559a997 authored by Anton Akhmerov's avatar Anton Akhmerov
Browse files

migrate one more figure to code

parent 34101aa4
No related branches found
No related tags found
No related merge requests found
Pipeline #
Hide whitespace changes
Inline Side-by-side
code/band_structures.py
+ 21
− 0
View file @ 8559a997
......@@ -147,6 +147,26 @@ def DOS_2masses():
pyplot.savefig('phonons8.svg')
def consistency_1_2_masses():
k = np.linspace(0, 2*pi, 300)
k_dos = np.linspace(0, pi, 20)
fig, ax = pyplot.subplots()
ax.plot(k, dispersion_2m(k, acoustic=False), label='optical')
ax.plot(k, dispersion_2m(k), label='acoustic')
omega_max = dispersion_2m(0, acoustic=False)
ax.plot(k, omega_max * np.sin(k/4), label='equal masses')
ax.set_xlabel('$ka$')
ax.set_ylabel(r'$\omega$')
ax.set_xticks([0, pi, 2*pi])
ax.set_xticklabels(['$0$', r'$\pi/2a$', r'$\pi/a$'])
ax.set_yticks([], [])
ax.set_xlim(-pi/8, 2*pi+.4)
ax.set_ylim((0, dispersion_2m(0, acoustic=False) + .2))
ax.legend(loc='lower right')
pyplot.vlines([pi, 2*pi], 0, 2.2, linestyles='dashed')
draw_classic_axes(ax, xlabeloffset=.2)
pyplot.savefig('phonons7.svg')
def DOS_finite_phonon_chain(n, output_name):
rhs = 2 * np.eye(n) - np.eye(n, k=1) - np.eye(n, k=-1)
......@@ -173,6 +193,7 @@ def main():
fermi_dirac()
phonons_1d()
phonons_1d_2masses()
consistency_1_2_masses()
DOS_2masses()
meff_1d_tb()
tight_binding_1d()
......
docs/figures/phonons7.svg deleted 100644 → 0
+ 0
− 131
View file @ 34101aa4
<?xml version="1.0" encoding="utf-8"?>
<!-- Generator: Adobe Illustrator 21.0.0, SVG Export Plug-In . SVG Version: 6.00 Build 0) -->
<svg version="1.1" id="Layer_1" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink" x="0px" y="0px"
viewBox="0 0 425.2 105.3" style="enable-background:new 0 0 425.2 105.3;" xml:space="preserve">
<style type="text/css">
.st0{fill:none;stroke:#000000;stroke-width:0.5;stroke-miterlimit:10;}
.st1{fill:none;stroke:#F7941D;stroke-width:0.5;stroke-miterlimit:10;}
.st2{fill:#FFFFFF;}
.st3{fill:none;stroke:#000000;stroke-width:0.5;stroke-miterlimit:10;stroke-dasharray:2;}
.st4{clip-path:url(#SVGID_2_);}
.st5{clip-path:url(#SVGID_4_);}
.st6{clip-path:url(#SVGID_6_);}
.st7{clip-path:url(#SVGID_8_);fill:#F7941D;}
</style>
<g>
<path class="st0" d="M212.3,27.1c6.4-13,16.8-13,23.2,0"/>
<path class="st0" d="M258.7,27.1c-6.4,13-16.8,13-23.2,0"/>
<path class="st0" d="M258.7,27.1c6.4-13,16.8-13,23.2,0"/>
<path class="st0" d="M305.2,27.1c-6.4,13-16.8,13-23.2,0"/>
<path class="st0" d="M165.8,27.1c6.4-13,16.8-13,23.2,0"/>
<path class="st0" d="M212.3,27.1c-6.4,13-16.8,13-23.2,0"/>
<path class="st0" d="M119.3,27.1c6.4-13,16.8-13,23.2,0"/>
<path class="st0" d="M165.8,27.1c-6.4,13-16.8,13-23.2,0"/>
<path class="st0" d="M305.2,27.1c6.4-13,16.8-13,23.2,0"/>
</g>
<path class="st0" d="M176.8,84c0,0,14.8-34,23.2-34c8.4,0,23.2,34,23.2,34"/>
<path class="st1" d="M176.8,84c0,0,29.4-66.4,45.7-66.5c16.3-0.1,44.6,65.7,44.6,65.7"/>
<path class="st0" d="M223.3,84c0,0,14.8-34,23.2-34c8.4,0,23.2,34,23.2,34"/>
<path class="st0" d="M130.3,84c0,0,14.8-34,23.2-34c8.4,0,23.2,34,23.2,34"/>
<rect x="223.3" y="3.3" class="st2" width="114.7" height="86.2"/>
<rect x="61.7" y="3.3" class="st2" width="115.1" height="86.2"/>
<g>
<g>
<line class="st0" x1="176.8" y1="84" x2="176.8" y2="9.7"/>
<g>
<polygon points="178.3,10.7 176.8,5.3 175.4,10.7 "/>
</g>
</g>
</g>
<line class="st3" x1="223.3" y1="84" x2="223.3" y2="5.3"/>
<line class="st3" x1="200" y1="84" x2="200" y2="5.3"/>
<g>
<defs>
<rect id="SVGID_1_" x="168" y="5.7" width="4.8" height="3.7"/>
</defs>
<clipPath id="SVGID_2_">
<use xlink:href="#SVGID_1_" style="overflow:visible;"/>
</clipPath>
<path class="st4" d="M172.9,6.4c0-0.4-0.1-0.6-0.3-0.6c-0.2,0-0.4,0.2-0.4,0.4c0,0.1,0,0.2,0.1,0.2c0.1,0.1,0.3,0.3,0.3,0.7
c0,0.3-0.2,0.8-0.4,1.1c-0.2,0.3-0.5,0.6-0.9,0.6c-0.5,0-0.7-0.3-0.8-0.7c0.1-0.2,0.3-0.7,0.3-0.9c0-0.1,0-0.2-0.1-0.2
c-0.1,0-0.2,0-0.2,0.1c-0.1,0.2-0.2,0.7-0.2,1c-0.3,0.4-0.6,0.7-1.1,0.7c-0.5,0-0.7-0.5-0.7-0.9c0-1,0.8-1.8,0.8-1.9
c0-0.1-0.1-0.1-0.1-0.1c-0.1,0-0.2,0.1-0.2,0.2c-0.4,0.6-0.7,1.5-0.7,2.2c0,0.5,0.2,1.1,0.8,1.1c0.6,0,0.9-0.4,1.2-0.8
c0.1,0.5,0.4,0.8,0.9,0.8c0.6,0,1-0.5,1.3-1.1C172.6,7.9,172.9,6.9,172.9,6.4"/>
</g>
<g>
<defs>
<rect id="SVGID_3_" x="245.2" y="89.5" width="3.7" height="5.7"/>
</defs>
<clipPath id="SVGID_4_">
<use xlink:href="#SVGID_3_" style="overflow:visible;"/>
</clipPath>
<path class="st5" d="M247.2,89.7c0,0,0-0.1-0.1-0.1c-0.2,0-0.8,0.1-1,0.1c-0.1,0-0.2,0-0.2,0.2c0,0.1,0.1,0.1,0.2,0.1
c0.4,0,0.4,0.1,0.4,0.1l0,0.2l-1.1,4.6c0,0.1,0,0.1,0,0.2c0,0.2,0.2,0.2,0.2,0.2c0.1,0,0.2-0.1,0.3-0.2c0-0.1,0.4-1.5,0.4-1.7
c0.3,0,0.9,0.2,0.9,0.7c0,0.1,0,0.1,0,0.2c0,0.1,0,0.2,0,0.3c0,0.5,0.3,0.8,0.7,0.8c0.2,0,0.5-0.1,0.6-0.4c0.2-0.3,0.3-0.8,0.3-0.8
c0-0.1-0.1-0.1-0.1-0.1c-0.1,0-0.1,0-0.1,0.1c-0.2,0.6-0.3,1-0.7,1c-0.1,0-0.3-0.1-0.3-0.4c0-0.1,0-0.3,0.1-0.4c0-0.1,0-0.2,0-0.2
c0-0.5-0.5-0.7-1.2-0.8c0.2-0.1,0.5-0.4,0.7-0.6c0.4-0.4,0.7-0.8,1.1-0.8c0.1,0,0.1,0,0.1,0c0.1,0,0.1,0,0.2,0.1c0,0,0,0,0,0
c-0.4,0-0.5,0.3-0.5,0.4c0,0.1,0.1,0.3,0.3,0.3c0.2,0,0.4-0.2,0.4-0.5c0-0.2-0.2-0.5-0.5-0.5c-0.2,0-0.6,0.1-1.2,0.7
c-0.3,0.3-0.6,0.6-0.9,0.7L247.2,89.7"/>
</g>
<g>
<defs>
<rect id="SVGID_5_" x="197.6" y="87.5" width="4.8" height="14.5"/>
</defs>
<clipPath id="SVGID_6_">
<use xlink:href="#SVGID_5_" style="overflow:visible;"/>
</clipPath>
<path class="st6" d="M199.7,88.1h0.9c-0.3,1.2-0.3,1.5-0.3,2c0,0.1,0,0.3,0.1,0.6c0.1,0.4,0.2,0.4,0.3,0.4c0.2,0,0.3-0.1,0.3-0.3
c0,0,0-0.1-0.1-0.2c-0.2-0.6-0.2-1.1-0.2-1.3c0-0.4,0.1-0.8,0.1-1.3h0.9c0.1,0,0.4,0,0.4-0.3c0-0.2-0.2-0.2-0.3-0.2h-2.7
c-0.2,0-0.5,0-0.8,0.4c-0.3,0.3-0.5,0.7-0.5,0.7c0,0,0,0.1,0.1,0.1c0.1,0,0.1,0,0.1-0.1c0.4-0.6,0.8-0.6,1-0.6h0.5
c-0.3,1-0.7,1.9-1,2.6c-0.1,0.1-0.1,0.1-0.1,0.2c0,0.2,0.1,0.2,0.2,0.2c0.2,0,0.3-0.2,0.4-0.5c0.1-0.4,0.1-0.4,0.2-0.8L199.7,88.1"
/>
<rect x="197.6" y="94.2" class="st6" width="4.8" height="0.3"/>
<path class="st6" d="M200.4,100.8c0,0.1,0,0.2-0.2,0.3c-0.4,0.4-0.7,0.6-0.9,0.6c-0.4,0-0.5-0.4-0.5-0.7c0-0.4,0.3-1.4,0.4-1.7
c0.2-0.5,0.6-0.8,0.9-0.8c0.5,0,0.6,0.6,0.6,0.7c0,0.1,0,0.1,0,0.1L200.4,100.8 M200.9,98.8c-0.1-0.3-0.4-0.5-0.7-0.5
c-0.9,0-1.9,1.2-1.9,2.3c0,0.7,0.4,1.3,1.1,1.3c0.2,0,0.6,0,1-0.6c0.1,0.3,0.3,0.6,0.7,0.6c0.3,0,0.5-0.2,0.6-0.4
c0.1-0.3,0.2-0.8,0.2-0.8c0-0.1-0.1-0.1-0.1-0.1c-0.1,0-0.1,0-0.1,0.1c-0.1,0.5-0.3,1-0.6,1c-0.2,0-0.2-0.2-0.2-0.4
c0-0.2,0-0.2,0.1-0.6c0.1-0.3,0.1-0.4,0.2-0.7l0.3-1.1c0.1-0.2,0.1-0.2,0.1-0.3c0-0.1-0.1-0.2-0.2-0.2
C201,98.4,200.9,98.6,200.9,98.8z"/>
</g>
<g>
<g>
<line class="st0" x1="176.8" y1="84" x2="244.8" y2="84"/>
<g>
<polygon points="243.8,85.4 249.2,84 243.8,82.5 "/>
</g>
</g>
</g>
<g>
<defs>
<rect id="SVGID_7_" x="220.8" y="87.5" width="4.8" height="14.5"/>
</defs>
<clipPath id="SVGID_8_">
<use xlink:href="#SVGID_7_" style="overflow:visible;"/>
</clipPath>
<path class="st7" d="M223,88.1h0.9c-0.3,1.2-0.3,1.5-0.3,2c0,0.1,0,0.3,0.1,0.6c0.1,0.4,0.2,0.4,0.3,0.4c0.2,0,0.3-0.1,0.3-0.3
c0,0,0-0.1-0.1-0.2c-0.2-0.6-0.2-1.1-0.2-1.3c0-0.4,0.1-0.8,0.1-1.3h0.9c0.1,0,0.4,0,0.4-0.3c0-0.2-0.2-0.2-0.3-0.2h-2.7
c-0.2,0-0.5,0-0.8,0.4c-0.3,0.3-0.5,0.7-0.5,0.7c0,0,0,0.1,0.1,0.1c0.1,0,0.1,0,0.1-0.1c0.4-0.6,0.8-0.6,1-0.6h0.5
c-0.3,1-0.7,1.9-1,2.6c-0.1,0.1-0.1,0.1-0.1,0.2c0,0.2,0.1,0.2,0.2,0.2c0.2,0,0.3-0.2,0.4-0.5c0.1-0.4,0.1-0.4,0.2-0.8L223,88.1"/>
<rect x="220.9" y="94.2" class="st7" width="4.8" height="0.3"/>
<path class="st7" d="M223.6,100.8c0,0.1,0,0.2-0.2,0.3c-0.4,0.4-0.7,0.6-0.9,0.6c-0.4,0-0.5-0.4-0.5-0.7c0-0.4,0.3-1.4,0.4-1.7
c0.2-0.5,0.6-0.8,0.9-0.8c0.5,0,0.6,0.6,0.6,0.7c0,0.1,0,0.1,0,0.1L223.6,100.8 M224.1,98.8c-0.1-0.3-0.4-0.5-0.7-0.5
c-0.9,0-1.9,1.2-1.9,2.3c0,0.7,0.4,1.3,1.1,1.3c0.2,0,0.6,0,1-0.6c0.1,0.3,0.3,0.6,0.7,0.6c0.3,0,0.5-0.2,0.6-0.4
c0.1-0.3,0.2-0.8,0.2-0.8c0-0.1-0.1-0.1-0.1-0.1c-0.1,0-0.1,0-0.1,0.1c-0.1,0.5-0.3,1-0.6,1c-0.2,0-0.2-0.2-0.2-0.4
c0-0.2,0-0.2,0.1-0.6c0.1-0.3,0.1-0.4,0.2-0.7l0.3-1.1c0.1-0.2,0.1-0.2,0.1-0.3c0-0.1-0.1-0.2-0.2-0.2
C224.3,98.4,224.1,98.6,224.1,98.8z"/>
</g>
<g>
</g>
<g>
</g>
<g>
</g>
<g>
</g>
<g>
</g>
<g>
</g>
<g>
</g>
</svg>
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment