def getPhGS(eta): gs = arbOrder.findGS(eta, 3) gsJ = eF.J(gs) gsT = eF.T(gs) gsJ = 0. gsT = -2. / np.pi * prms.chainLength return phState.findPhotonGS([gsT, gsJ], eta, 3)
def expectAnnihil(eta): gs = arbOrder.findGS(eta, 3) gsJ = eF.J(gs) gsT = eF.T(gs) ptGS = phState.findPhotonGS([gsT, gsJ], eta, 3) a = setUpAnnihil(eta) return np.dot(np.conj(ptGS), np.dot(a, ptGS))
def currentGS(eta): gs = arbOrder.findGS(eta, 3) gsJ = eF.J(gs) gsT = eF.T(gs) ptGS = phState.findPhotonGS([gsT, gsJ], eta, 3) jOp = setupCurrentOperator(gsT, gsJ, eta) current = np.dot(np.conj(ptGS), np.dot(jOp, ptGS)) return current
def expectSinA(eta): gs = arbOrder.findGS(eta, 3) gsJ = eF.J(gs) gsT = eF.T(gs) ptGS = phState.findPhotonGS([gsT, gsJ], eta, 3) x = setUpA(eta) sinX = sciLin.sinm(x) return np.dot(np.conj(ptGS), np.dot(sinX, ptGS))
def getPhGS(eta): gs = arbOrder.findGS(eta, 3) gsJ = eF.J(gs) gsT = eF.T(gs) return phState.findPhotonGS([gsT, gsJ], eta, 3)