Numerical generation of hyperspherical harmonics for tetra-atomic systems

A numerical generation method of hyperspherical harmonics for tetra-atomic systems, in terms of row-orthonormal hyperspherical coordinates-a hyper-radius and eight angles-is presented. The nine-dimensional coordinate space is split into three three-dimensional spaces, the physical rotation, kinematic rotation, and kinematic invariant spaces. The eight-angle principal-axes-of-inertia hyperspherical harmonics are expanded in Wigner rotation matrices for the physical and kinematic rotation angles. The remaining two-angle harmonics defined in kinematic invariant space are expanded in a basis of trigonometric functions, and the diagonalization of the kinetic energy operator in this basis provides highly accurate harmonics. This trigonometric basis is chosen to provide a mathematically exact and finite expansion for the harmonics. Individually, each basis function does not satisfy appropriate boundary conditions at the poles of the kinetic energy operator; however, the numerically generated linear combination of these functions which constitutes the harmonic does. The size of this basis is minimized using the symmetries of the system, in particular, internal symmetries, involving different sets of coordinates in nine-dimensional space corresponding to the same physical configuration.     

A high dimensional model representation based numerical method for solving ordinary differential equations

A new numerical method for solving ordinary differential equations by using High Dimensional Model Representation (HDMR) has been developed in this work. Higher order ordinary differential equations can be reduced to a set of first order ODEs. Although HDMR is generally used for multivariate functions, univariate functions are taken into account throughout the work because of the ODEs’ natures. Not the numerical solution but its image under an appropriately chosen linear ordinary differential operator is expressed as a linear combination of the positive deviation powers of independent variable from its initial value. The linear combination of these image functions are expected to form a basis set under consideration. The unknown constants in the linear combination are found by maximizing the constancy measurer formed in terms of the HDMR components after they are evaluated. Results are compared with well-known step size based numerical methods. A semi qualitative error analysis of the proposed method is also established.     

Symmetrical overlap transformations of function basis sets: the LCAO MO and quantum similarity practical cases

Quantum chemical computational procedures, like LCAO MO theory and quantum similarity, use non orthogonal function basis sets, which define finite dimensional subspaces of a Hilbert space. Based on the original overlap metric matrices, generated by the chosen finite non orthogonal basis sets, there are several symmetrical overlap and basis set transformations possible. This study tries to find out the general point of view, from where all these procedures can be studied in a clear generalized perspective.     

Symmetry-adapted integrals over many-electron basis functions and operators

We consider integrals over symmetry-adapted basis functions that involve the coordinates of more than one electron. We focus on basis functions that can be written as products of one-electron functions and (say) a two-electron function. We show first that the two-electron parts of the basis functions can be absorbed into the operator, resulting in an integral over only one-electron basis functions, but a more complicated many-electron operator. We then prove a general formula for expressing such integrals in terms of symmetry-distinct integrals only. Received: 16 June 2000 / Accepted: 10 July 2000 / Published online: 19 January 2001     

Constrained active space unrestricted mean-field methods for controlling spin-contamination

We have recently proposed a novel approach for obtaining high-spin restricted open-shell Hartree-Fock wave functions by imposing constraints on the unrestricted Hartree-Fock (UHF) method [T. Tsuchimochi and G. E. Scuseria, J. Chem. Phys. 133, 141102 (2010)]. We here extend these ideas to the case where the constraints are released in an active space but imposed elsewhere. If the active space is properly chosen, our constrained UHF (CUHF) method greatly benefits from a controlled broken-symmetry effect while avoiding the massive spin contamination of traditional UHF. We also revisit and apply Lo?wdin's projection operator to CUHF and obtain multireference wave functions with moderate computational cost. We report singlet-triplet energy splittings showing that our constrained scheme outperforms fully unrestricted methods. This constrained approach can be readily used in spin density functional theory with similar favorable effects.     

Rapidly converging threshold value calculations in screened coulomb potential systems: Critical values of the screening parameter for the Yukawa case

In this work, the threshold values of the screening parameter for Yukawa potential systems are obtained by means of a basis set constructed from Laguerre type functions. The Laguerre basis set is modified by an appropriately chosen extra function in order to imitate the true behaviour of the solutions at the boundary points. The method used is a variational scheme and the numerical results are accurate to thirty decimal points. A scaling parameter is also inserted into the structure of the basis functions, the optimized values of which accelerate the convergence. The main goal of this paper is to develop a method which enables us to calculate the threshold values of the screening parameter for low-lying states. The method is quite general and can be extended to all systems whose potentials decay exponentially when the radial variable goes to infinity.     

The symmetrized random matrix approach, an efficient method to obtain relativistic molecular symmetry adapted functions

In relativistic quantum chemical calculation of molecules, where the spin-orbit interaction is included, the electron orbitals possess both the double point group symmetry and the time-reversal symmetry. If symmetry adapted functions are employed as the basis functions of electron orbitals, it would allow a significant reduction of the computational expense. The point group symmetry adapted functions can be obtained by the group projection operators via its actions on the atomic orbital functions. We have proposed an efficient and simple method to obtain all irreducible representation matrices, which are the basis of the group projection operators, of any finite double point group. Both double point group symmetry and time-reversal symmetry are automatically imposed on the representation matrices. This is achieved by the symmetrized random matrix (SRM) approach, where the SRM is constructed in the regular representation space of a finite group and the eigenfunctions of SRM provide all irreducible representation matrices of the given point group.     

A multidimensional discrete variable representation basis obtained by simultaneous diagonalization

Direct product basis functions are frequently used in quantum dynamics calculations, but they are poor in the sense that many such functions are required to converge a spectrum, compute a rate constant, etc. Much better, contracted, basis functions, that account for coupling between coordinates, can be obtained by diagonalizing reduced dimension Hamiltonians. If a direct product basis is used, it is advantageous to use discrete variable representation (DVR) basis functions because matrix representations of functions of coordinates are diagonal in the DVR. By diagonalizing matrices representing coordinates it is straightforward to obtain the DVR that corresponds to any direct product basis. Because contracted basis functions are eigenfunctions of reduced dimension Hamiltonians that include coupling terms they are not direct product functions. The advantages of contracted basis functions and the advantages of the DVR therefore appear to be mutually exclusive. A DVR that corresponds to contracted functions is unknown. In this paper we propose such a DVR. It spans the same space as a contracted basis, but in it matrix representations of coordinates are diagonal. The DVR basis functions are chosen to achieve maximal diagonality of coordinate matrices. We assess the accuracy of this DVR by applying it to model four-dimensional problems.     

Simple method to obtain symmetry harmonics of point groups

We propose a simple, self-consistent method to obtain basis functions of irreducible representations of a finite point group. Our method is based on eigenproblem formulation of a projection operator represented as a nonhomogeneous polynomial of angular momentum L. The method is shown to be more efficient than the usual numerical methods when applied to the analysis of high-order symmetry harmonics in cubic and icosahedral groups. For low-order symmetry harmonics the method provides rational coefficients of expansion in the Y(L,M) basis.     

The non-additive ligand field

The semi-empirical ligand field is a perturbation operator whose consequences are taken to first order using a basis set ofl functions. Since the basis spans an irreducible representation of the 3-dimensional rotation-inversion groupR _3i it is useful to express the operator as a sum of components of irreducible tensor operators with respect to this group. IfR _3i is reduced with respect to the molecular subgroup the electronic factor of each term in the sum must be totally symmetrical within this group. This choice of operator leads to thecrystal field parameterization without implying an electrostatic model. Alternatively a shift operator withinl space may be chosen as the essential part of the perturbation operator. This leads to theligand field parameterization. Between the two parameterizations there exists a one to one relationship, whose coefficients are proportional to 3l symbols. This relationship is given together with a brief discussion of the reasons for the proposed parameterizations.     

Basis set limit Hartree-Fock and density functional theory response property evaluation by multiresolution multiwavelet basis

We describe the evaluation of response properties using multiresolution multiwavelet (MRMW) basis sets. The algorithm uses direct projection of the perturbed density operator onto the zeroth order density operator on the real space spanned by the MRMW basis set and is applied for evaluating the polarizability of small molecules using Hartree-Fock and Kohn-Sham density functional theory. The computed polarizabilities can be considered to be converged to effectively complete space within the requested precision. The efficiency of the method against the ordinary Gaussian basis computation is discussed.     

A simplified representation of the potential produced by Gaussian charge distributions

A procedure is presented which allows a more economical representation of the potential produced by orbital charge distributions in which the orbitals are expanded in terms of a finite set of polynomial Gaussian functions. The basic idea is that the products of pairs of Gaussian basis functions, on which the charge distributions are expanded, are expressed in terms of a new basis set of optimally chosen single Gaussian functions. Such a procedure has been tested in a particular case and a few possible applications have been suggested.     

The kinetic energy operator in the subspaces of wavelet analysis

At any resolution level of wavelet expansions the physical observable of the kinetic energy is represented by an infinite matrix which is “canonically” chosen as the projection of the operator  − Δ/2 onto the subspace of the given resolution. It is shown, that this canonical choice is not optimal, as the regular grid of the basis set introduces an artificial consequence of its periodicity, and it is only a particular member of possible operator representations. We present an explicit method of preparing a near optimal kinetic energy matrix which leads to more appropriate results in numerical wavelet based calculations. This construction works even in those cases, where the usual definition is unusable (i.e., the derivative of the basis functions does not exist). It is also shown, that building an effective kinetic energy matrix is equivalent to the renormalization of the kinetic energy by a momentum dependent effective mass compensating for artificial periodicity effects.     

Biased sampling of nonequilibrium trajectories: can fast switching simulations outperform conventional free energy calculation methods?

We have investigated the maximum computational efficiency of reversible work calculations that change control parameters in a finite amount of time. Because relevant nonequilibrium averages are slow to converge, a bias on the sampling of trajectories can be beneficial. Such a bias, however, can also be employed in conventional methods for computing reversible work, such as thermodynamic integration or umbrella sampling. We present numerical results for a simple one-dimensional model and for a Widom insertion in a soft sphere liquid, indicating that, with an appropriately chosen bias, conventional methods are in fact more efficient. We describe an analogy between nonequilibrium dynamics and mappings between equilibrium ensembles, which suggests that the practical inferiority of fast switching is quite general. Finally, we discuss the relevance of adiabatic invariants in slowly driven Hamiltonian systems for the application of Jarzynski's theorem.     

A discrete look at localization

This article introduces a set of localized orthonormal functions to serve as basis functions for quantum calculations. They are defined to be eigenfunctions of the position operator in a function space. Their properties, including their variances, for a one-dimensional system are developed. The application to simple harmonic motion is considered as an example and, in particular, the time evolution of an initially localized function is calculated and shown to be periodic. The theory can be interpreted as producing a discrete quantization of space with Hamiltonian interactions that are predominantly between nearest neighbors. These functions can also be used in approximate calculations. To illustrate their accuracy, the example of a Morse oscillator treated as a perturbation of a harmonic oscillator is reconsidered. It is shown that the localized functions in a variational calculation lead to a result that is a good approximation for the lowest states. Furthermore, the use of a wave function that is defined only at discrete points can be justified as the first approximation to this, so that its accuracy can also be discussed. © 1995 John Wiley & Sons, Inc.     

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 probabilistic evolution approach trilogy, part 1: quantum expectation value evolutions, block triangularity and conicality, truncation approximants and their convergence

This is the first one of three companion papers focusing on the “probabilistic evolution approach (PEA)” which has been developed for the solution of the explicit ODE involving problems under certain consistent impositions. The main purpose here is the determination of the expectation value of a given operator in quantum mechanics by solving only ODEs, not directly using the wave function. To this end we first define a basis operator set over the Kronecker powers of an appropriately defined “system operator vector”. We assume that the target operator’s commutator with the system’s Hamiltonian can be expressed in terms of the above-mentioned basis operators. This assumption leads us to an infinite set of linear homogeneous ODEs over the expectation values of the basis operators. Its coefficient matrix is in block Hessenberg form when the target operator has no singularity, and beyond that, it may become block triangular when certain conditions over the system’s potential function are satisfied. The initial conditions are the basic determining agents giving the probabilistic nature to the solutions of the obtained infinite set of ODEs. They may or may not have fluctuations depending on the nature of the probability density. All these issues are investigated in a phenomenological and constructive theoretical manner in this paper. The remaining two papers are devoted to further details of PEA in quantum mechanics, and, the application of PEA to systems defined by Liouville equation.     

Toward Exact Analytical Wave Function of Helium Atom: Two Techniques for Constructing Homogeneous Functions of Kinetic Energy Operator

We discuss two general techniques for constructing homogeneous functions c of the kinetic energy operator T for the helium atom in a state of symmetry S. The first technique is based on algebraic identification of the kernel of T in a space spanned by some predetermined set of basis functions. The second technique, analytic in nature, constructs the homogeneous functions of T as formal power series with coefficients deduced from recurrence relations stemming from the requirement Tc=0. Both approaches are capable of producing a great variety of homogeneous functions c with arbitrary homogeneity that can prove useful for constructing the exact ground state wave function for the helium atom.     

Density of non-Bloch electron states in perfect cubic crystals

This paper gives an abbreviated method for the calculation of the density of states of a crystal on the basis of that band theory in which the crystal electron states are represented by the standinglike wave functions classified according to the point-group symmetry species. The crystal is a large but finite sphere filled regularly with atoms, and the wave functions are quantized at the boundary of the sphere. The Bloch theorem is not satisfied in this theory since the wave functions are not basis functions of the irreducible representations of the translation subgroup. On the other hand, a theorem is established that the density of states can be made up of contributions given by all irreducible representations of the crystal point group, any contribution being proportional to the square of the dimension of the irreducible representation. In distinction to a former approach, the band structure is calculated solely from the energy eigenvalues obtained with the aid of the diagonalization process of the Wannier–Slater differential operator. A simple cubic lattice with an s atomic orbital on each lattice site is taken as an example, and the results are compared with Bloch's theory.     

The coupling algebra of trigonally subduced crystal field eigenvectors

The general theory of subduction of eigenvectors between infinite groups is used to derive a finite group subduction operator and define the corresponding subduction coefficients. The coupling behaviour of these subduced eigenvectors can then be described in terms of 3 Γ symbols. These symbols, defined only in relation to complex basis sets are all fully real and have all phases fixed by the subduction operator. They differ from V coefficients in two phase relationships and have the advantage, unlike V coefficients, of retaining all the symmetry properties and selection rules of Wigner 3-j symbols. Appropriate label systems which render these properties in terms of simple algebras are given for all quantizing axes available in O _h . The specific set of 3Γ symbols for each quantization is determined by the orientation of the coordinate axes in the Hamiltonian. The four possible orientations for trigonal quantization are examined and the operator chosen which produces eigenvectors with conventional conjugate phases and a fully real set of 3Γ symbols.