def get_w0_00_D(H_00_D):
return \
WaveFunction(states=H_00_D.states, init_state=[0, 0, [0, 1, 2], [0, 0]], amplitude=1./sqrt(6)) + \
WaveFunction(states=H_00_D.states, init_state=[0, 0, [1, 2, 0], [0, 0]], amplitude=1./sqrt(6)) + \
WaveFunction(states=H_00_D.states, init_state=[0, 0, [2, 0, 1], [0, 0]], amplitude=1./sqrt(6)) - \
WaveFunction(states=H_00_D.states, init_state=[0, 0, [0, 2, 1], [0, 0]], amplitude=1./sqrt(6)) - \
WaveFunction(states=H_00_D.states, init_state=[0, 0, [1, 0, 2], [0, 0]], amplitude=1./sqrt(6)) - \
WaveFunction(states=H_00_D.states, init_state=[0, 0, [2, 1, 0], [0, 0]], amplitude=1./sqrt(6))
def get_w0_11_000(H_11_000):
return WaveFunction(states=H_11_000.states,
init_state=[1, 1, [0, 0, 0], [0, 0]])
def get_w0_11_2D(H_11_2D):
return \
WaveFunction(states=H_11_2D.states, init_state=[1, 1, [0, 1, 0], [0, 0]], amplitude=1./sqrt(2)) - \
WaveFunction(states=H_11_2D.states, init_state=[
1, 1, [1, 0, 0], [0, 0]], amplitude=1./sqrt(2))
def get_w0_01_1D(H_01_1D):
return \
WaveFunction(states=H_01_1D.states, init_state=[0, 1, [0, 0, 1], [0, 0]], amplitude=1./sqrt(2)) - \
WaveFunction(states=H_01_1D.states, init_state=[
0, 1, [1, 0, 0], [0, 0]], amplitude=1./sqrt(2))