#!/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())