wf_x.value = wf_x.value * 0.05
Ex = electron_state_2.energy_levels[0]
Ex.convert_to('eV')
electron_state_1.static_potential.convert_to('eV')
n = spsc_data.DensityValue(electron_state_1.sum_density.value,
electron_state_1.sum_density.units)
n.convert_to(n.units_default)
electron_state_1.mass.convert_to(electron_state_1.mass.units_default)
well_start = electron_state_1.static_potential.meta_info['well_start']
g = electron_state_1.mass.value[well_start] / (np.pi *
spsc_constants.h_plank**2)
E1.convert_to(E1.units_default)
E2.convert_to(E2.units_default)
Ef = spsc_data.EnergyValue(0.5 * (n.value / g + E1.value + E2.value))
Ef1 = Ef - E1
Ef2 = Ef - E2
Ef.convert_to('eV')
E1.convert_to('eV')
E2.convert_to('eV')
Ef1.convert_to('eV')
Ef2.convert_to('eV')
x = np.linspace(0, length.value, len(potential_1))
np.savetxt("z_C236.csv", x, delimiter=",")
ax = plt.subplot(111)
# Hide the right and top spines
ax.spines['right'].set_visible(False)
ax.spines['top'].set_visible(False)
ax.set_xlabel('Z, nm', fontsize=16)
def test_j_to_erg(self):
value = spsc_data.EnergyValue(8, "J")
value.convert_to("erg")
self.assertEqual(value.value, 8 * 10**7)
def test_erg_to_ert(self):
value = spsc_data.EnergyValue(5, "erg")
value.convert_to("erg")
self.assertEqual(value.value, 5)
def test_erg_to_ev(self):
value = spsc_data.EnergyValue(7, "erg")
value.convert_to("eV")
assert_value = value.value / 10**11
self.assertAlmostEqual(assert_value, 7 * 6.24150965)
potential_2.convert_to('eV')
wf_x = electron_state_2.wave_functions[0]
wf_x.value = wf_x.value * 0.1
Ex = electron_state_2.energy_levels[0]
Ex.convert_to('eV')
electron_state_1.static_potential.convert_to('eV')
n = spsc_data.DensityValue(6.81 * 10**15, 'm^-2')
n.convert_to(n.units_default)
electron_state_1.mass.convert_to(electron_state_1.mass.units_default)
well_start = electron_state_1.static_potential.meta_info['well_start']
g = electron_state_1.mass.value[well_start] / (np.pi *
spsc_constants.h_plank**2)
E1.convert_to(E1.units_default)
Ef = spsc_data.EnergyValue(n.value / g + E1.value)
Ef1 = Ef - E1
Ef.convert_to('eV')
E1.convert_to('eV')
Ef1.convert_to('eV')
x = np.linspace(0, length.value, len(potential_1))
ax = plt.subplot(111)
# Hide the right and top spines
ax.spines['right'].set_visible(False)
ax.spines['top'].set_visible(False)
ax.set_xlabel('Z, nm', fontsize=16)
ax.set_ylabel('V, eV', fontsize=16)
ax.plot(x, potential_1, label=r'$\Gamma$')
ax.plot(x, potential_2, label='X')
ax.plot(x, wf_1, label=r'$\psi_\Gamma$')