Elementary parameter-differentiation techniques are developed to systematically derive a wide variety of operator identities, expansions, and solutions to differential equations of interest to quantum physics. The treatment is largely centered around a general closed formula for the derivative of an exponential operator with respect to a parameter. Derivations are given of the Baker-Campbell-Hausdorff formula and its dual, the Zassenhaus formula. The continuous analogs of these formulas which solve the differential equation dY(t)/dt = A(t) Y(t), the solutions of Magnus and Fer, respectively, are similarly derived in a recursive manner which manifestly displays the general repeated-commutator nature of these expansions and which is quite suitable for computer programming. An expansion recently obtained by Kumar and another new expansion are shown to be derivable from the Fer and Magnus solutions, respectively, in the same way. Useful similarity transformations involving linear combinations of elements of a Lie algebra are obtained. Some cases where the product eAeB can be written as a closed-form single exponential are considered which generalize results of Sack and of Weiss and Maradudin. Closed-form single-exponential solutions to the differential equation dY(t)/dt = A(t) Y(t) are obtained for two cases and compared with the corresponding multiple-exponential solutions of Wei and Norman. Normal ordering of operators is also treated and derivations, corollaries, or generalization of a number of known results are efficiently obtained. Higher derivatives of exponential and general operators are discussed by means of a formula due to Poincaré which is the operator analog of the Cauchy integral formula of complex variable theory. It is shown how results obtained by Aizu for matrix elements and traces of derivatives may be readily derived from the Poincaré formula. Some applications of the results of this paper to quantum statistics and to the Weyl prescription for converting a classical function to a quantum operator are given. A corollary to a theorem of Bloch is obtained which permits one to obtain harmonic-oscillator canonical-ensemble averages of general operators defined by the Weyl prescription. Solutions of the density-matrix equation are also discussed. It is shown that an initially canonical ensemble behaves as though its temperature remains constant with a ``canonical distribution'' determined by a certain fictitious Hamiltonian.
Elementary parameter-differentiation techniques are developed to systematically derive a wide variety of operator identities, expansions, and solutions to differential equations of interest to quantum physics. The treatment is largely centered around a general closed formula for the derivative of an exponential operator with respect to a parameter. Derivations are given of the Baker-Campbell-Hausdorff formula and its dual, the Zassenhaus formula. The continuous analogs of these formulas which solve the differential equation dY(t)/dt = A(t) Y(t), the solutions of Magnus and Fer, respectively, are similarly derived in a recursive manner which manifestly displays the general repeated-commutator nature of these expansions and which is quite suitable for computer programming. An expansion recently obtained by Kumar and another new expansion are shown to be derivable from the Fer and Magnus solutions, respectively, in the same way. Useful similarity transformations involving linear combinations of elements of a Lie algebra are obtained. Some cases where the product eAeB can be written as a closed-form single exponential are considered which generalize results of Sack and of Weiss and Maradudin. Closed-form single-exponential solutions to the differential equation dY(t)/dt = A(t) Y(t) are obtained for two cases and compared with the corresponding multiple-exponential solutions of Wei and Norman. Normal ordering of operators is also treated and derivations, corollaries, or generalization of a number of known results are efficiently obtained. Higher derivatives of exponential and general operators are discussed by means of a formula due to Poincaré which is the operator analog of the Cauchy integral formula of complex variable theory. It is shown how results obtained by Aizu for matrix elements and traces of derivatives may be readily derived from the Poincaré formula. Some applications of the results of this paper to quantum statistics and to the Weyl prescription for converting a classical function to a quantum operator are given. A corollary to a theorem of Bloch is obtained which permits one to obtain harmonic-oscillator canonical-ensemble averages of general operators defined by the Weyl prescription. Solutions of the density-matrix equation are also discussed. It is shown that an initially canonical ensemble behaves as though its temperature remains constant with a ``canonical distribution'' determined by a certain fictitious Hamiltonian.
This paper studies the exponential of the sum of two non-commuting operators as an infinite product of exponential operators involving repeated commutators of increasing order. It will be shown how to determine the two coefficients in front of the nested commutators in the Zassenhaus formula. The knowledge of one coefficient is enough to generate closed formulas that have several applications in solving problems ranging from linear differential equations and quantum mechanics to non-linear differential equations.
Lie algebra, Zassenhaus, and parameter differentiation techniques are utilized to break up the exponential of a bilinear Hamiltonian operator into a product of noncommuting exponential operators by the virtue of the theory of Wei and Norman [J. Math. Phys. 4, 575 (1963); Proc. Am. Math. Soc., 15, 327 (1964)]. There are about three different ways to find the Zassenhaus exponents, namely, binomial expansion, Suzuki formula, and q-exponential transformation. A fourth, and most reliable method, is provided. Since linearly displaced and distorted (curvature change upon excitation/emission) Hamiltonian and spin-boson Hamiltonian may be classified as bilinear Hamiltonians, the presented algebraic algorithm (exponential operator disentanglement exploiting six-dimensional Lie algebra case) should be useful in spin-boson problems. The linearly displaced and distorted Hamiltonian exponential is only treated here. While the spin-boson model is used here only as a demonstration of the idea, the herein approach is more general and powerful than the specific example treated. The optical linear dipole moment correlation function is algebraically derived using the above mentioned methods and coherent states. Coherent states are eigenvectors of the bosonic lowering operator a and not of the raising operator a(+). While exp(a(+)) translates coherent states, exp(a(+)a(+)) operation on coherent states has always been a challenge, as a(+) has no eigenvectors. Three approaches, and the results, of that operation are provided. Linear absorption spectra are derived, calculated, and discussed. The linear dipole moment correlation function for the pure quadratic coupling case is expressed in terms of Legendre polynomials to better show the even vibronic transitions in the absorption spectrum. Comparison of the present line shapes to those calculated by other methods is provided. Franck-Condon factors for both linear and quadratic couplings are exactly accounted for by the herein calculated linear absorption spectra. This new methodology should easily pave the way to calculating the four-point correlation function, F(tau(1),tau(2),tau(3),tau(4)), of which the optical nonlinear response function may be procured, as evaluating F(tau(1),tau(2),tau(3),tau(4)) is only evaluating the optical linear dipole moment correlation function iteratively over different time intervals, which should allow calculating various optical nonlinear temporal/spectral signals.
Abstract The exponential of an operator or matrix is widely used in quantum theory, but it sometimes can be a challenge to evaluate. For non-commutative operators X and Y, according to the Campbell-Baker-Hausdorff-Dynkin theorem, e X+Y is not equivalent to e X e Y , but is instead given by the well-known infinite series formula. For a Lie algebra of a basis of three operators {X, Y, Z}, such that [X, Y] = κ Z for scalar κ and cyclic permutations, here it is proven that e a X+b Y is equivalent to e p Z e q X e−p Z for scalar p and q. Extensions for e a X+b Y+c Z are also provided. This method is useful for the dynamics of atomic and molecular nuclear and electronic spins in constant and oscillatory transverse magnetic and electric fields.
We derive a useful expression for the matrix elements [∂f[A(t)]/∂t]ij of the derivative of a function f[A(t)] of a diagonalizable linear operator A(t) with respect to the parameter t at t0. The function f[A(t)] is supposed to be an operator acting on the same space as the operator A(t) which is assumed to have a nondegenerate, pure point spectrum. We use the basis which diagonalizes A(t0), i.e., [A(t0]ij=λiδij, and obtain [∂f[A(t)]/∂t|t=t0]ij=[∂A/∂t|t=t0]ij{[f(λj)−f(λi)]/(λj−λi)}. In addition to this, we show that further elaboration on the (not necessarily simple) integral expressions given by Wilcox (who basically considered f[A(t)] of the exponential type) and generalized by Rajagopal [who extended Wilcox results by considering f[A(t)] of the q-exponential type where expq(x)≡[1+(1−q)x]1/(1−q) with q∈R; hence, exp1(x)=exp(x)] yields these same expressions. Some of the lemmas first established by the above authors are easily recovered.
Ordered operator expansions for operators forming physically important low-dimensional Lie algebras are derived in a simple unified way. Starting with the Zassenhaus formula for the disentangling of exponential operators, series expansions of both undisentangled and disentangled exponentials and comparison of the operator coefficients of equal powers of an ordering parameter alpha leads to ordered operator expansions. This 'comparison method' gives an alternative simple derivation of some already known formulae and a number of new formulae in the physical and chemical applications of the harmonic oscillator and for master equation problems with nearest-neighbour transition probabilities. The 'comparison method' cannot be applied to the angular momentum algebra directly. By a slight modification it can be used to derive from one matrix element or trace of J x k J y l J z n all possible combinations k, l, n by simply comparing powers of ordering parameters.
The Baker--Campbell--Hausdorff (BCH) formula plays a critical role in many branches of mathematics and physics. It expresses the logarithm of the product of exponentials of non-commuting operators as an infinite series of nested commutators of the operators involved. The Zassenhaus formula is the dual of the BCH formula: the exponential of a sum of operators is written as an infinite product of exponentials involving the operators and their commutators. In practical computations, however, one typically has to truncate the expansions, and so understanding the error committed by the resulting approximations and eventually providing suitable bounds for this error is of paramount interest. In this work we present a general strategy to derive rigorous error bounds and explicit error constants for the BCH and Zassenhaus formulas when the operators involved are skew-adjoint, as is the case for quantum evolution problems.