""" use the symetry of lineshape function in the exciton indices """ #D = CumulantExpr(expr)._leading_index(a) #expr = D._getExpr() A = Anorm.rewrite(gg) norm = CumulantExpr(A) """ use option large=T to calculate evaluate in T --> oo """ norm = norm.evaluate(large=T) """ use the symetry of lineshape function in the exciton indices """ #D = CumulantExpr(expr)._leading_index(a) #expr = D._getExpr() expr = (expr-norm).simplify() if verbatim: print("Cumulant: ") print(" ") print(expr) print(" ") print("Foerster rate") expr = expr.subs(e,d).simplify() a1 = Symbol('a1') b1 = Symbol('b1') expr = expr.subs({d:b1,c:a1,b:a1,a:b1}) expr = expr.subs({a1:a,b1:b}).simplify() print(expr)
A = A.rewrite(gg) expr = CumulantExpr(A) """ use option large=T to calculate evaluate in T --> oo """ expr = expr.evaluate(large=T) """ use the symetry of lineshape function in the exciton indices """ #D = CumulantExpr(expr)._leading_index(a) #expr = D._getExpr() A = Anorm.rewrite(gg) norm = CumulantExpr(A) """ use option large=T to calculate evaluate in T --> oo """ norm = norm.evaluate(large=T) """ use the symetry of lineshape function in the exciton indices """ #D = CumulantExpr(expr)._leading_index(a) #expr = D._getExpr() expr = (expr-norm).simplify() if verbatim: print("Cumulant: ") print(" ") print(expr) print(" ") print("In simplified form: ") print(" ") expr = expr.subs(c,b) expr = expr.subs(d,b).simplify() print(expr)
expr = expr.evaluate(large=T) """ use the symetry of lineshape function in the exciton indices """ #D = CumulantExpr(expr)._leading_index(a) #expr = D._getExpr() A = Anorm.rewrite(gg) norm = CumulantExpr(A) """ use option large=T to calculate evaluate in T --> oo """ norm = norm.evaluate(large=T) """ use the symetry of lineshape function in the exciton indices """ #D = CumulantExpr(expr)._leading_index(a) #expr = D._getExpr() expr = (expr - norm).simplify() if verbatim: print("Cumulant: ") print(" ") print(expr) print(" ") print("Foerster rate") expr = expr.subs(e, d).simplify() a1 = Symbol('a1') b1 = Symbol('b1') expr = expr.subs({d: b1, c: a1, b: a1, a: b1}) expr = expr.subs({a1: a, b1: b}).simplify() print(expr)