Unitary group tensor operator algebras for many-electron systems: II. One- and two-body matrix elements

Relying on our earlier results in the unitary group Racah-Wigner algebra, specifically designed to facilitate quantum chemical calculations of molecular electronic structure, the tensor operator formalism required for an efficient evaluation of one- and two-body matrix elements of molecular electronic Hamiltonians within the spin-adapted Gel'fand-Tsetlin basis is developed. Introducing the second quantization-like creation and annihilation vector operators at the unitary group [U(n)] level, appropriate two-box symmetric and antisymmetric irreducible tensor operators as well as adjoint tensors are defined and their matrix elements evaluated in the electronic Gel'fand-Tsetlin basis as single products of segment values. Using these tensor operators, the matrix elements of one- and two-body components of a general electronic Hamiltonian are found. Explicit expressions for all relevant quantities pertaining to at most two-column irreducible representations that are required in molecular electronic structure calculations are given. Relationships with other approaches and possible future extensions of the formalism to partitioned bases or spin-dependent Hamiltonians are discussed.On leave from: Department of Chemistry, Xiamen University, Xiamen, Fujian, PR China.     

Study of simulation method of time evolution of atomic and molecular systems by quantum electrodynamics

We discuss a method to follow step‐by‐step time evolution of atomic and molecular systems based on quantum electrodynamics. Our strategy includes expanding the electron field operator by localized wavepackets to define creation and annihilation operators and following the time evolution using the equations of motion of the field operator in the Heisenberg picture. We first derive a time evolution equation for the excitation operator, the product of two creation or annihilation operators, which is necessary for constructing operators of physical quantities such as the electronic charge density operator. We, then, describe our approximation methods to obtain time differential equations of the electronic density matrix, which is defined as the expectation value of the excitation operator. By solving the equations numerically, we show “electron‐positron oscillations,” the fluctuations originated from virtual electron‐positron pair creations and annihilations, appear in the charge density of a hydrogen atom and molecule. We also show that the period of the electron‐positron oscillations becomes shorter by including the self‐energy process, in which the electron emits a photon and then absorbs it again, and it can be interpreted as the increase in the electron mass due to the self‐energy. © 2014 Wiley Periodicals, Inc.     

New realizations of electronic configuration interaction spaces

The general properties of new models of the electronic spaces based on the notion of (p, q)‐sheaves are studied. The interrelation between simple sheaves and density operators is established. Explicit expressions for the transformed reduced Hamiltonians in terms of the standard creation‐annihilation operators are presented. The general scheme of parametrization of p‐electron states by κ‐electron means (κ = 2, 3, …) is described and studied in detail for the case of sheaves induced by κ‐electron wave functions. It is demonstrated that, under certain conditions, a p‐electron problem may be reformulated as the eigenvalue problem in κ‐electron space equipped with certain p‐electron metric. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Transition probabilities for a time-dependent bilinear Hamiltonian with two harmonic limits

A method is developed for calculating transition probabilities for a time-dependent bilinear Hamiltonian that yields two different Harmonic oscillators at initial and final times. The procedure is based on the equations of motion for the creation and annihilation operators and the transition probabilities are obtained from recurrence relations that are most suitable for computational purposes.     

Method of computer algebraic calculation of the matrix elements in the second quantization language

An automated method by the algebraic programming language REDUCE3 for specifying the matrix elements expressed in second quantization language is presented and then applied to the case of the matrix elements in the TDHF theory. This program works in a very straightforward way by commuting the electron creation and annihilation operators (a^? and a) until these operators have completely vanished from the expression of the matrix element under the appropriate elimination conditions. An improved method using singlet generators of unitary transformations in the place of the electron creation and annihilation operators is also presented. This improvement reduces the time and memory required for the calculation. These methods will make programming in the field of quantum chemistry much easier. © 1995 John Wiley & Sons, Inc.     

Solutions to the time-dependent,forced quantum oscillator equation in the coordinate representation

Time-dependent creation and annihilation operators are derived which are used to obtain exact solutions, in the coordinate representation, to the time-dependent, forced quantum oscillator equation. The solutions are used to obtain a general formula for the transition probabilities, valid for any time-dependent force.     

Algebraic modifications to second quantization for non‐Hermitian complex scaled hamiltonians with application to a quadratically convergent multiconfigurational self‐consistent field method

The algebraic structure for creation and annihilation operators defined on orthogonal orbitals is generalized to permit easy development of bound‐state techniques involving the use of non‐Hermitian Hamiltonians arising from the use of complex‐scaling or complex‐absorbing potentials in the treatment of electron scattering resonances. These extensions are made possible by an orthogonal transformation of complex biorthogonal orbitals and states as opposed to the customary unitary transformation of real orthogonal orbitals and states and preserve all other formal and numerical simplicities of existing bound‐state methods. The ease of application is demonstrated by deriving the modified equations for implementation of a quadratically convergent multiconfigurational self‐consistent field (MCSCF) method for complex‐scaled Hamiltonians but the generalizations are equally applicable for the extension of other techniques such as single and multireference coupled cluster (CC) and many‐body perturbation theory (MBPT) methods for their use in the treatment of resonances. This extends the domain of applicability of MCSCF, CC, MBPT, and methods based on MCSCF states to an accurate treatment of resonances while still using L² real basis sets. Modification of all other bound‐state methods and codes should be similarly straightforward. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

On a field-theoretic formulation of many-electron spin-adaptation

We report the development of a set of spin-adapted creation and annihilation operators which provides a suitable generalisation of the usual many-body apparatus.     

Harmonic oscillator tensors. II. angular momentum expressions of matrix elements of vibrational operators

The spatial symmetries of the harmonic oscillator and the recently found irreducible tensors constructed from the associated annihilation and creation operators are exploited to obtain new expressions for the elements of the matrix representatives of several examples of vibrational operators. Since all vibrational operators are expressible in terms of the irreducible tensors, their matrix elements reflect the angular momentum symmetry inherent in them, for the results derived here are in terms of the Clebsch–Gordan coefficients and the isoscalar factors that arise from the couplinig rule of the irreducible tensors. Familiarity with the mathematical properties of these quantities derived from the elementary theory of angular momentum facilitates the evaluation of many vibrational operators that may be of importance in the study of potentials in this basis. In particular, it is shown that the nonvanishing of matrix elements is governed by a law of conservation of angular momentum along the axis of quantization of the nondegenerate harmonic oscillator. © 1993 John Wiley & Sons, Inc.     

Method of recurrent construction of Löuwdin spin-adapted wave functions. II. Local representation of creation–annihilation operators

Some compositions of the addition and subtraction operators and recurrence relations for the Sanibel-type coefficients c_{u, v} (n, s, M) generated by these compositions are studied. A local representation of the fermion creation–annihilation operators via the addition and subtraction operators is obtained. Operators of single excitations, coupling, and decoupling operators, in terms of which the unitary group generators can be expressed are defined. The resulting representation of the nonelementary unitary group generators is much more simple than in the Gelfand–Tzetlin basis and in the most general case contains only six logically different terms, each of them possessing quite transparent physical significance.     

Unitary group tensor operator algebras for many-electron systems. III. Matrix elements in U(n 1 +n 2) ⊃ U(n 1) × U(n 2) partitioned basis

Exploiting our earlier results [J. Math. Chem. 4 (1990) 295–353 and 13 (1993) 273–316] on the unitary group U(n) Racah-Wigner algebra, specifically designed for quantum chemical calculations of molecular electronic structure, and the related tensor operator formalism that enabled us to introduce spin-free orbital equivalents of the second quantization-like creation and annihilation operators as well as higher rank symmetric, antisymmetric and adjoint tensors, we consider the problem of U(n) basis partitioning that is required for group-function type approaches to the many-electron problem. Using the U(n) U(n ₁) × U(n ₂),n =n ₁ +n ₂ adapted basis, we evaluate all required matrix elements of U(n) generators and their products that arise in one- and two-body components of non-relativistic electronic Hamiltonians. The formalism employed naturally leads to a segmented form of these matrix elements, with many of the segments being identical to those of the standard unitary group approach. Relationship with similar approaches described earlier is briefly pointed out.     

Kratzer potential algebraic representation and matrix elements recurrence formulae

Generalized recurrence relations for the calculation of multipole matrix elements for Kratzer potential wave functions are obtained operationally. These formulas have been determined by using a non-analytical procedure based on the algebraic representation of the Kratzer eigenfunctions along with the usual ladder properties and commutation relations. For that, the creation and annihilation operators are adequately derived by means of an alternative approach to the factorization method and the exact expressions for matrix elements are achieved with the aid of a relationship between the ladder operators associated with the bra and theket. The proposed algebraic approach as well as the formulas for the calculation of matrix elements thus derived are quite simple and direct when compared with other alternative expressions already obtained analytically or pseudo-algebraically by means of the hypervirial theorem commutator algebra.     

Vibronic energies and spectra of molecular dimers

We consider three distinct methods of calculating the vibronic levels and absorption spectra of molecular dimers coupled by dipole-dipole interactions. The first method is direct diagonalization of the vibronic Hamiltonian in a basis of monomer eigenstates. The second method is to use creation and annihilation operators leading in harmonic approximation to the Jaynes-Cummings Hamiltonian. The adiabatic approximation to this problem provides insight into spectral behavior in the weak and strong coupling limits. The third method, which serves as a check on the accuracy of the previous methods, is a numerically exact solution of the time-dependent Schrodinger equation. Using these methods, dimer spectra are calculated for three separate dye molecules and show good agreement with measured spectra.     

Improving an algebraic approach to calculate approximate Franck–Condon factors of diatomic molecules

Using the BCH theorem, we express the Hamiltonian of a Morse oscillator as a complete series of powers of the creation and annihilation operators for the harmonic oscillator. In this way, we improve the results of a previous work that uses a Bogoliubov–Tyablikov tranformation to calculate the Franck–Condon factors by means of equivalent harmonic oscillators potentials. © 1992 John Wiley & Sons, Inc.     

Indole‐Fused Epoxy‐1,5‐Diazocines as Cancer‐Selective Cytotoxic Agents

Dimerization of 2‐(diformylmethylene)‐3,3‐dimethylindole by the action of tosyl chloride (TsCl) led to the creation of an epoxy‐[1,5]‐diazocine bicycle. The structure of the molecule resembles that of the naturally occurring C‐curarine‐I. The molecule shows significant cytotoxicity against cancer lines MCF‐7 and MDA‐MB‐231, but not toward the normal cell line CCD‐841. A series of substituted 2‐(diformylmethylene)‐3,3‐dimethylindoles was accordingly dimerized, and the products were studied for their cytotoxic activities.     

Geometry of density matrices. VI. Superoperators and unitary invariance

We examine and compare ways of dividing into subspaces the space whose elements are density matrices or other operators for the class of model problems defined by a finite one-particle basis set. One method of decomposition makes the significance of the subspaces apparent. We show that this decomposition is also complete, in the group-theoretic sense, for the group of unitary transformations of the set of one-electron basis functions. The irreducible subspaces are labeled by particle number and by an additional integer we call the reduction index. For spaces of particle-number-conserving operators, all subspaces with the same reduction index are isomorphic, and an analogous isomorphism exists for non-particle-number-conserving cases. The general linear group also plays a key role, and we introduce the term “canonical superoperators” to characterize those superoperators which commute with this group. When an appropriate basis set is chosen for the matrix spaces, the supermatrices corresponding to these superoperators have a particularly simple form: a block structure with the only nonzero blocks being multiples of unit matrices. The superoperators of interest can be constructed in terms of two operators, , and these two have been expressed simply in terms of creation and annihilation operators. When only real orthogonal transformations of the basis are considered, a further decomposition is possible. We have introduced superoperators associated with this decomposition.     

A general expression for matrix elements of reciprocal powers of the displacement coordinateq in the harmonic oscillator basis

A formula is derived that allows one to determine the matrix elements of an arbitrary integral reciprocal power of the dimensionless displacement coordinate q of the harmonic oscillator from those ofq ^–1 in an exact manner. This relation is obtained from the use of the chain rule and irreducible tensors expressed in terms of the creation and annihilation operators of the harmonic oscillator.     

Direct configuration interaction with a reference state composed of many reference configurations

The configuration interaction method where a single reference state is composed of a linear combination of reference configurations is analyzed in detail. In this method single and double replacements are constructed by applying annihilation and creation operators on the reference state. The analysis is based on the recently derived factorization of the direct CI coupling coefficients into internal and external parts. Using the internal coupling coefficients the integrals are transformed to new entities which are used in the diagonalization step. This two-step procedure differs significantly from the usual straightforward one-step direct CI procedure. A number of operations analysis shows that calculations with the present method should be feasible even with a large number of reference configurations in the reference state. Based on first-order perturbation theory the accuracy of the method is predicted to be close to the accuracy obtained with the usual CI method with many reference configurations.