Padé approximants to physical properties via inner projections

The [N, M] Padé approximants to functions formally associated to power series expanssions are expressed in terms of expectation values of inverse matrices. These formulae, which can be derived by the inner-projection technique, lead to a simple analysis of the properties of serveral approximation methods and their inter-relationships, in particular Gaussian integration, continued factorization and Padé approximations, which are of current interest in the calculation of physical properties. A relation with Fredholm integral equations and expansions of the resolvent is also discussed. The use of operator inequalities in a systematic fashion is particularly convenient when both the function being approximated and the coefficients of the power series have physically meaningful expressions as moments of operators.     

Generalized Euler transformation in extracting useful information from divergent (asymptotic) perturbation series and the construction of Padé approximants

Euler transformation for accelerating convergence of a series is considered in the context of handling divergent (asymptotically convergent) perturbation series. A generalized (parametrized) version of this transformation is developed, based on the conjecture of Dalgarno and Stewart, which works better. Viewed from this standpoint, the Padé approximants follow as a special case of the parametrized Euler transformation (PET ), as is the case with the μ transformation procedure of Feenberg in a perturbative context. The PET is shown to serve as a more general method of handling a divergent series and is able to appreciate the construction and convergence behavior of specific sequences of Padé approximants. The role of parametrization in the context of the Z^?1 perturbation theory of atoms is also noted and the workability of the adopted strategy is demonstrated by choosing some specific test cases.     

Electrostatic energy of polygonal charge distributions

An asymptotic series for the electrostatic energy E₁(N){\mathcal{E}_1(N)} of an N-gonal charge distribution, i.e., a set of unit charges occupying vertices of a regular N-gon with a unit circumradius, is derived. Application of Padé approximants to truncations of this expansion produces compact approximate formulae capable of estimating E₁(N){\mathcal{E}_1(N)} with great accuracy. A closed-form expression for the energy of electrostatic interaction of two polygonal charge distributions is obtained from the respective Fourier series. The availability of this expression allows for a rapid calculation of the relevant energy with computational effort independent of the numbers of particles involved.     

Fifth-order constant denominator perturbation theory within a localized bond model: Contributions from triple and quadruple excitations

The contributions of the triple and quadruple excitations to the fifth-order perturbation energy for the perturbation configuration interaction using localized orbitals (PCILO ) method are derived. This completes the development of a fifth-order constant denominator perturbation theory initiated in a previous paper [5] with the single and double excitations. This theory is tested on molecules containing strained ring geometries, stretched bonds, strongly polarized bonds, and delocalized pi systems: cases where the starting zero order reference wave function poorly describes the system. Although the perturbation expansions turn out to be slowly convergent, the Padé approximant taken from an energy series which itself is constructed from Padé approximants provides results accurate to within a few kilocalories/mole of benchmark calculations. Computational times as in the original PCILO procedure remain proportional to N³, where N is the number of bonds.     

Critical point,singularities and extrapolations in the helium iso-electronic sequence

The Z-expansion of two-electron systems is analyzed with the Padé technique with emphasis on establishing analytical properties of the function E(Z) formally associated with the power series expansion. The concept of critical point in this connection is stressed. For this sequence it occurs at Z_c = 0.911246 with E(Z_c) = ?0.415184. The structure of E(Z) for Z < Z_c is investigated. The use of Padé approximants to extrapolate values of electron affinities is emphasized.     

Accurate estimates of asymptotic indices via fractional calculus

We devise a three-parameter random search strategy to obtain accurate estimates of the large-coupling amplitude and exponent of an observable from its divergent Taylor expansion, known to some desired order. The endeavor exploits the power of fractional calculus, aided by an auxiliary series and subsequent construction of Padé approximants. Pilot calculations on the ground-state energy perturbation series of the octic anharmonic oscillator reveal the spectacular performance.     

Additive self-similar approximants

A novel type of approximants is introduced, being based on the ideas of self-similar approximation theory, allowing for accurate extrapolation and interpolation of asymptotic series. The method is illustrated by the examples possessing the structure typical of many nonlinear problems in mathematical chemistry. Good numerical convergence is demonstrated for the cases that can be compared with exact solutions, when these are available. The method is shown to be not less and as a rule essentially more accurate than that of Padé approximants. Comparison with other approximation methods is also given.     

Comparison of high-order many-body perturbation theory and configuration interaction for H2O

Diagrammatic many-body perturbation theory, coupled with a recursive computational procedure, is employed to obtain the correlation energy of H₂O within a 39-STO basis set by evaluating all double-excitation diagrams through twelfth order without any approximations. This provides, for the first time, the complete double-excitation diagrams contributions to the correlation energy, which is ?0.28826 hartree, compared with a correlation energy of ?0.27402 hartree obtained from a configuration interaction calculation which includes all double excitations. The difference of 0.0142 hartree includes the “size consistency” correction to the all-double-excitations CI energy, due to the “pathological” unliked-diagram terms remaining in that result, but also involves certain fourth- and higher-order rearrangement diagrams. Contrary to conventional belief, the unshifted, or Møller-Plesset partitioning of the hamiltonian provides a much more rapid convergence of the perturbation series that does the shifted, or Epstein-Nesbet partitioning. In both cases. Padé approximants enhance the convergence of the series considerably. A simple variation-perturbation scheme based on the first-order MBPT wavefunction is sufficient to provide 97.5% of the all-doubles CI correlation energy.     

Equations of state for tetra-dimensional hard-sphere fluids

An equation of state, a kind of generalised Padé approximants, first proposed for the hard-sphere fluid in cases of two, three, four and five dimensions is extended for the tetra-dimensional case with new simulation data available. The corresponding equations of state show good to excellent agreement with the above-mentioned data.     

Electron correlation effects in the 4f14 shell

This paper serves a twofold purpose. First, Löwdin's inner projection in both nonperturbative and perturbative forms is applied to the quartic anharmonic oscillator. Inner projection with perturbation theory yields rational approximations to Brillouin–Wigner-type perturbation expansions. These lower bounds are compared with [N ? 1, N] Padé approximants to the Rayleigh–Schrödinger perturbation series for this problem. These Padés are also expressible as the even convergents, w_2N, of a Stieltjes-type continued fraction. The latter representation has certain advantages with respect to its Padé counterpart. Inner projection without perturbation theory provides significantly better results than the perturbative version. The application of inner projection techniques to a perturbed hydrogen atom is not straightforward. The usual problems associated with the continuum spectrum of hydrogen are present. By means of a nonunitary “tilting” transformation associated with the Lie group SO(4, 2), these problems may be bypassed. In the SO(4, 2)-reformulated eigenvalue problem, a reinterpretation of the basic variables, as developed by Silverstone and Moats, yields a new Hamiltonian that permits direct use of the inner projection method. This method has been applied to the ground state of the hydrogen atom in a magnetic field, using both four- and eight-dimensional basis manifolds. This represents the first application of inner projection to this problem.     

Variational functions of the padé type—application to the ground state of the helium atom

It is shown that trial functions involving Padé approximants yield satisfactory results for the ground state of the helium atom. In particular, the five-parameter form reproduces the best variational function of the type Ψ = e^–zs?(u) obtained numerically, to a remarkable extent.     

Estimation of eigenvalues using padé approximants. Stark effect for a model system

The perturbation series for a one-dimensional model system in a static electric field is obtained and shown to be divergent. Padé approximants to the series are calculated and found to converge, and to provide a good estimation of the Stark eigenvalue.     

Second-order dispersion energy between interacting hydrogen pairs by using the Borel integral

The second-order dispersion energy between interacting hydrogen pairs has been evaluated by the Borel-integral method. The results are in complete harmony with those obtained by the method of Padé approximants.     

Applications of Padé approximants of type II in partitioning technique

Padé approximants of type II have been applied to solve equations of the type F(z) = 0. The method is compared with the first-order iteration procedure, the Aitken-Samuelson formula, and the Newton-Raphson tangential method. As a test example the partitioning technique in its most simple form is applied to the Hamiltonian of a rigid symmetric-top molecule in a static electric field. The proposed algorithm is found to be superior to the first-order procedure and the Aitken-Samuleson formula, and at least as effective as the Newton-Raphson method.     

Perturbation for a rigid rotator in an electric field

The behavior of a rigid rotor perturbed by an electric field is studied. An hypervirial calculation, to obtain the perturbation series without calculation of the perturbed wavefunctions, is used to determine the Stark-shifted eigenvalues. A numerical estimation of the m equals; 0 → m = 1 transition for the lowest vibrational state (v = 0) and the first rotational state (J = 1) for the HF and DF molecules using Padé approximants is made.     

Diagrammatic perturbation theory: Potential curves for the ground state of the carbon monoxide molecule

Diagrammatic many-body perturbation theory is used to calculate the potential energy function for the X¹ σ+ state of the CO molecule near the equilibrium nuclear configuration. Spectroscopic constants are derived from a number of curves which are obtained from calculations taken through third order in the energy. By forming [2/1] Padé approximants to the constants we obtain: r_e = 1.125 Å (1.128 Å), B_e = 1.943 cm^?1 (1.9312 cm^?1), a^B_e = 0.0156 cm^?1 (0.0175 cm^?1), W_e = 2247 cm^?1 (2170 cm^?1), W_eX_e = 12.16 cm^?1 (13.29 cm^?1), where the experimental values are given in parenthesis.     

Extrapolation of the total energy of polymers to the bulk limit using generalized padé approximants

Generalized Padé approximants are used to extrapolate the total energy of polymers described by alternant Hamiltonians to the bulk limit. The method provides an upper bound to the energy. The origin of (quasi) periodic oscillations of the energy per unit cell as the function of the number of unit cells is enlightened through analysis of the moments of the Hamiltonian.     

Bounding the extrapolated correlation energy using Padé approximants

We present a rigorous strategy, based on Stieltjes series and Padé approximants, to obtain suitable bounds for extrapolation of the quantum chemical correlation energy. Computational tests are performed for the second‐order Møller–Plesset (MP2) correlation energy, and the bounds obtained are tight enough for practical calculational purposes: The associated error in most cases is much less than 1 kcal/mol. The bounds presented here are also shown to be rigorous for functional forms that represent a wide variety of methods in quantum chemistry and hence may be used in extrapolating a wide range of expressions, some of them yielding significant computational advantages compared to traditional techniques. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 79: 222–234, 2000     

Quadratic Padé approximants and the intruder state problem of multireference perturbation methods

Simple and quadratic Padé resummation methods are applied to high‐order series from multireference many‐body perturbation theory (MR‐MBPT) calculations using various partitioning schemes (Møller–Plesset, Epstein–Nesbet, and forced degeneracy) to determine their efficacy in resumming slowly convergent or divergent series. The calculations are performed for the ground and low‐lying excited states of (i) CH₂, (ii) BeH₂ at three geometries, and (iii) Be, for which full configuration interaction (CI) calculations are available for comparison. The 49 perturbation series that are analyzed include those with oscillatory and monotonic divergence and convergence, including divergences that arise from either frontdoor or backdoor intruder states. Both the simple and quadratic Padé approximations are found to speed the convergence of slowly convergent or divergent series. However, the quadratic Padé method generally outperforms the simple Padé resummation. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

The convergence of energy expansions for molecules in electrostatic fields: A linear‐response coupled‐cluster study

The radii of convergence of power series expansions describing energy of a molecule in external electrostatic field are investigated usingD’Alembert ratio test, standard and generalized Cauchy–Hadamardcriteria, and Padé approximants. The corresponding coefficients at various field and field‐gradient components, representing multipole moments and (hyper)polarizabilities and including terms of tenth or even twentieth order, are determined using an ab initio linear responsecoupled‐cluster theory. Most calculations are performed for the HF molecule described by the basis set of double zeta quality, while the role of basis set is discussed by comparing the results with estimates of the radii of convergence obtained with the basis set of [5s3p2d/3s2p] quality. Emphasis is placed on the dependence of the interval of convergence of power series expansion describing energy of a molecule in applied electrostatic field on the nuclear geometry. The results might have important implications for various numerical methods used to calculate electrostatic molecular properties as functions of the internuclear geometry, including the finite‐field andfixed‐point‐charge approaches. This revised version was published online in July 2006 with corrections to the Cover Date.