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))