#!/usr/bin/env python3
from flat import Flat
f = Flat(dumping=True,
m="bigraphene",
Exc=1,
Hc=200,
Ex=0.3,
omega=5e10,
T=300,
n=100,
dt=1e-13,
alltime=1e-8)
f.run(True)
#!/usr/bin/env python3
import numpy as np
from flat import Flat
freqs = np.linspace(5.25e14, 5.35e14, 3)
waves = lambda freq: [[1, 0, 0, 1e12, 0], [10, 0, 0, freq, 0]]
population = []
for freq in freqs:
f = Flat(m="bigraphene_0.01_2",
waves=waves(freq),
T=300,
n=100,
dt=1e-15,
alltime=1e-9)
d = f.run()
population.append(d["population"])
population = np.array(population)
data = np.hstack([freqs.reshape(len(freqs), 1), population[:, :, 0]])
np.savetxt("ps.dat", data)
#!/usr/bin/env python3
import sys
import numpy as np
from matplotlib import pyplot as plt
from flat import Flat
Excl = np.linspace(0, 2, 6)
Hcl = np.linspace(20, 5000, 50)
for Exc in Excl:
line = []
for Hc in Hcl:
f = Flat(m="bigraphene_0.01_1",
Exc=Exc,
Hc=Hc,
waves=[[1, 0, 0, 1e12, 0]],
T=300,
n=100,
dt=1e-13,
alltime=1e-9)
line.append(f.run(True)["power"][1][0])
plt.plot(Hcl, line, label="Exc=%.2f" % Exc)
plt.legend()
plt.gca().set_ylim(bottom=0)
plt.savefig("magnetic.png")
#!/usr/bin/env python3
import sys
import numpy as np
from flat import Flat
result = np.zeros((2, 21))
result[0, :] = np.linspace(-2, 2, 21)
for i, Exc in enumerate(result[0, :]):
f = Flat(m="graphene", Exc=Exc, T=300, n=100, dt=1e-14, alltime=1e-8)
d = f.run(True)
result[1, i] = d["v_y"][0]
np.savetxt("cvc.dat", result.transpose())