An iterative method to solve the algebraic eigenvalue problem

An iterative method based on perturbation theory in matrix form is described as a procedure to obtain the eigenvalues and eigenvectors of square matrices. Practical vector notation and elementary schematic algorithm codes are given. The particular programming characteristics of the present computational scheme are based upon eigenvector corrections, obtained through a simple Rayleigh–Schrödinger perturbation theory algorithm. The proposed methodological processes can be used to evaluate the eigensystem of large matrices.     

Combination of perturbative and variational methods for calculating molecular spectra: calculation of the upsilon=3-5 CH stretch overtone spectrum of CHF3

Highly excited states of the CHF3 molecule belonging to the third, fourth, and fifth Fermi polyad are calculated using a combination of the Van Vleck perturbation theory and a variational treatment. The perturbation theory preconditions the Hamiltonian matrix by transforming away all couplings except those between nearly degenerate states. This transformation is implemented so that eigenvalues can be found with significantly smaller matrices than that which would be needed in the original normal mode representation. Even with preconditioning, at the energies as high as 3-5 quanta in the CH stretch, it is not possible to directly diagonalize the Hamiltonian matrix due to the large basis sets required. Iterative methods, particularly the block-Davidson method, are explored for finding the eigenvalues. The methods are compared and the advantages discussed.     

A new approach to the determination of good action-angle variables for coupled oscillator systems

A theory of action-angle variables for coupled oscillator systems is developed which involves solving the Schrödinger equation using a basis of WKB eigenfunctions, then using the logarithm of the resulting wavefunction to define the generator for the canonical transformation which determines the action-angle variables. This theory is based on the marriage between Miller's method for solving the Hamilton-Jacobi equation using the logarithm of the generating function, and the Ratner-Buch-Gerber method for solving the Schrödinger equation using WKB basis functions. A perturbation-theory analysis of this theory indicates that the semiclassical eigenvalues and canonical transformations obtained from it should become identical to their exact classical counterparts in the limit of large actions for each vibrational mode. Two methods for systematically improving the theory for the lower eigenstates are also proposed. Numerical applications of the theory are presented for two systems, the Morse oscillator and the Henon-Heiles two-mode hamiltonian. The resulting semiclassical eigenvalues are in excellent agreement with their exact quantum counterparts, with the magnitude of the error roughly independent of the energy of the eigenstate. Analogous good agreement is found in comparing the approximate and exact classical canonical transformations. In particular, for the Morse oscillator, good results are obtained for certain higher energy states where second-order classical perturbation theory makes serious errors. Other information examined includes surfaces of section for the Henon-Heiles system (comparing the analytical functions obtained from the present theory with results based on exact trajectory calculations) and vibrational distributions chosen to simulate trajectory calculations (using the present theory to determine bin boundaries for a histogram calculation). Again, the comparison in each case with accurate results is excellent, with maximum errors in action calculations of 0.02 h, and in angle calculations of 0.01 rad.     

Accurate eigenvalues of three- through ten-electron atomic isoelectronic sequences from low-order Z-dependent perturbation theory combined with variational constraints and screening

A method is developed, based on Rayleigh-Schrödinger perturbation theory combined with variational constraints and screening, for obtaining accurate atomic eigenvalues from third-order 1/Z expansions. Application of the procedure to the ground states of the 3NV10 electron atomic sequences yields energies of 99.95–100.05% or greater accuracy, a marked improvement over those obtained from other third-order summations including Padé approximants. In the important test cases of the Be and Ne atoms, our results are found to exceed in accuracy all but the most elaborate ab initio calculations.     

Asymptotic Functional Form Preservation Method for the Perturbation Theory: Application to Anharmonic Oscillators

The asymptotic functional form preservation method is developed for the perturbation theory to obtain the energy eigenvalues of anharmonic oscillators. The conventional energy perturbative series expansion for the anharmonic oscillator is strongly divergent even if the anharmonicity is small. Employing a transformation containing an unphysical parameter, we analytically continue this series expansion into a new series expansion applicable to all the range of the perturbation parameter. The unphysical parameter is determined by the principle of minimal sensitivity. This new series expansion is reduced to the conventional energy perturbative series expansion for small anharmonicity, and it preserves the correct asymptotic functional form when the perturbation parameter tends to infinity. Then, we use the full‐range energy series expansion to calculate the energy eigenvalues of the anharmonic oscillator. In addition to excellent energy eigenvalues obtained for the oscillator with small and strong anharmonicity, accurate energy eigenvalues can be obtained using the full‐range energy series expansion when the perturbation parameter tends to infinity.     

Effective Floquet Hamiltonian for spin I = 1 in magic angle spinning NMR using contact transformation

Contact transformation is an operator transformation method in time-independent perturbation theory which is used successfully in molecular spectroscopy to obtain an effective Hamiltonian. Floquet theory is used to transform the periodic time-dependent Hamiltonian, to a time-independent Floquet Hamiltonian. In this article contact transformation method has been used to get the analytical representation of Floquet Hamiltonian for quadrupolar nuclei with spin I = 1 in the presence of an RF field and first order quadrupolar interaction in magic angle spinning NMR experiments. The eigenvalues of contact transformed Hamiltonian as well as Floquet Hamiltonian have been calculated and a comparison is made between the eigenvalues obtained using the two Hamiltonians.     

Perturbational and variational treatments of the Morse oscillator

The Morse oscillator hamiltonian is expressed as an infinite expansion in powers of a natural perturbation parameter, the square root of the anharmonicity constant, relative to the simple harmonic oscillator as zeroth-order hamiltonian. A transformation of variables leads to a hamiltonian which involves terms no higher than second order in this natural perturbation parameter. In both cases, the exact bound state eigenvalues of the Morse oscillator are given by second-order perturbation theory. The Schrödinger equation corresponding to the transformed Morse hamiltonian is solved variationally, via a complete set expansion in simple harmonic oscillator eigenstates. Accurate approximations to the exact eigenvalues and eigenfunctions of bound states of the Morse oscillator can be obtained for all but the very highest levels.     

Numerical solution of Dalgarno-Lewis equations by a mapped Fourier grid method

Inhomogeneous radial differential equations emerging in applications of standard perturbation theory are numerically solved by a novel approach making use of Fourier grid methods in conjunction with a simple mapping scheme. The proposed algorithm is applied along the lines of the Dalgarno-Lewis method [Proc. R. Soc. London, Ser. A 223, 70 (1955)] to the calculation of the static dipole polarizabilities and hyperpolarizabilities of 1s, 2s, and 2p states of hydrogen atom and their frequency dependent dynamic dipole polarizabilities. The high efficiency and accuracy of the algorithm are demonstrated for the above test cases, where exact values are available. Then, the frequency dependent dipole polarizability of the ground state of lithium atom is computed by a variationally stable method combining an effective local potential approach with a second-order energy correction. The obtained results are in perfect agreement with other elaborate theoretical approaches.     

Convergence property of multireference many-body perturbation theory analyzed by the use of a norm of effective Hamiltonian

Summary It is proposed to use a norm of anth order effective Hamiltonian, for analyzing the convergence property of the multireference many-body perturbation theory (MR-MBPT). The utilization of the norm allows us to employ only (1) asingle number for all the states that we are interested in, and (2) values which decreases from thepositive side to zero as the ordern of the perturbation increases. This characteristic features are in contrast to those in the usually used scheme whereseveral numbers, namely, the eigenvalues of the target states, should be used and they mayoscillate around exact eigenvalues. The present method has been applied to MR-MBPT calculations of the (H₂)₂, CH₂, and LiH molecules based on the multireference versions of Rayleigh-Schrödinger PT, Kirtman-Certain-Hirschfelder PT, and the canonical Van Vleck PT; and following features are found: (1) the above three versions of the perturbation theories have essentially the same convergence property judged from the lowering of the norm; (2) the lower order truncation of the perturbation series gives reasonable solutions; (3) the norm decreases irrespective of the perturbation expansion being convergent or divergent for the first several orders (up to about the sixth order).     

Perturbation theory without wave function for multidimensional systems

A new method to obtain perturbation corrections to the eigenvalues of multidimensional quantummechanical models is developed. It consists of rearranging the Rayleigh-Schrodinger perturbation theory so that any coefficient of the perturbation series is obtained from a simple and compact recursion relationship. The Zeeman effect in hydrogen and the hydrogen molecule-ion are used to illustrate the procedure.     

Variational constraints on perturbed eigenvalues and their perturbation expansions. II. Perturbational inequalities

A large family of interlocking perturbational inequalities is derived by variational considerations for the stationary states of all systems described by a Hamiltonian linear in a real perturbing parameter λ. These inequalities constrain in many different ways the perturbation expansions of both exact and variational eigenvalues for these systems; analogous inequalities are derived for the components of the eigenvalues. A special feature of the analysis consists of obtaining inequalities applicable to the separate sums of even- and odd-order perturbation energies. For lowest states of each symmetry and for positive definite perturbing operator, the interlocking effect of the inequalities becomes extremely restrictive. The inequalities are illustrated with several numerical calculations for different systems and states of the helium isoelectronic sequence. The direction of the inequalities is found to be unaffected by low-order truncation, thus rendering them applicable to low-order perturbation expansions. The inequalities are used to study the efficacy of low-order perturbation theory for two- to ten-electron atomic isoelectronic sequences, and to determine the functional dependence upon λ of the eigenvalues and their components for arbitrary atomic isoelectronic sequences.     

Excited-state energy eigenvalue and wave-function evaluation of the Gaussian symmetric double-well potential problem via numerical shooting method 1

This work aims at computing excited-state energy eigenvalues and wave-function of a particle under Gaussian symmetric double-wells potential using numerical shooting method and perturbation theory a method to deal with discrete-eigenvalue problems. We also compare the energy eigenvalue and wave-function with those obtained from other typical means popular among physics students, namely the numerical shooting method and perturbation theory. Show that the idea of program of the numerical shooting method and perturbation theory of this problem (see Sects. 2.2 and 3). The numerical shooting method is generally regarded as one of the most efficient methods that give very accurate results because it integrates the Schrödinger equation directly, though in the numerical sense. The n = even case is shown in Fig. 5. In this case, the wave-function has split up on symmetric nodes under Gaussian symmetric double-well potential.     

Determination of molecular vibrational state energies using the ab initio semiclassical initial value representation: application to formaldehyde

We have demonstrated the use of ab initio molecular dynamics (AIMD) trajectories to compute the vibrational energy levels of molecular systems in the context of the semiclassical initial value representation (SC-IVR). A relatively low level of electronic structure theory (HF/3-21G) was used in this proof-of-principle study. Formaldehyde was used as a test case for the determination of accurate excited vibrational states. The AIMD-SC-IVR vibrational energies have been compared to those from curvilinear and rectilinear vibrational self-consistent field/vibrational configuration interaction with perturbation selected interactions-second-order perturbation theory (VSCF/VCIPSI-PT2) and correlation-corrected vibrational self-consistent field (cc-VSCF) methods. The survival amplitudes were obtained from selecting different reference wavefunctions using only a single set of molecular dynamics trajectories. We conclude that our approach is a further step in making the SC-IVR method a practical tool for first-principles quantum dynamics simulations.     

Renormalized perturbation theory by the moment method for degenerate states: Anharmonic oscillators

We apply renormalized perturbation theory by the moment method to an anharmonic oscillator in two dimensions with a perturbation that couples unperturbed degenerate states. The method leads to simple recurrence relations for the perturbation corrections to the energy and moments of the eigenfunction. We calculate accurate energy eigenvalues, illustrate the general features of the method, and comment on the application of the approach to other quantum mechanical models. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 66 : 261–272, 1998     

Vibrational quasi-degenerate perturbation theory: applications to fermi resonance in CO2, H2CO, and C6H6

A quasi-degenerate perturbation method with vibrational self-consistent field (VSCF) reference wavefunction is developed. It simultaneously accounts for strong anharmonic mode-mode coupling among a few states (static correlation) by a configuration interaction theory and for weak coupling with a vast number of the other states (dynamic correlation) by a perturbation theory. A general formula is derived based on the van Vleck perturbation theory. An algorithm that selects a compact set of the most important VSCF configurations which contribute to the static correlation is proposed and a scheme to limit the number of configurations considered for dynamic correlation is also implemented. This method reproduces the vibrational frequencies of CO2 and H2CO that are subject to the strongest anharmonic mode-mode coupling within 10 cm(-1) of vibrational configuration interaction results in a computational expense reduced by a factor of one to two orders of magnitude. The method also reproduces the infrared absorption of C6H6 in the CH stretching (nu12) frequency region, in which combination tones nu13nu16 and nu2nu13nu18 appear on account of an intensity borrowing from nu12via the anharmonic coupling.     

Excited-state energy eigenvalue and wave-function evaluation of the Gaussian asymmetric double-well potential problem via numerical shooting method 2

This project aims at computation excited-state energy eigenvalues and wave-function of a particle under Gaussian asymmetric double-well potential using numerical shooting method and perturbation theory a method to deal with discrete-eigenvalue problems. We also compare the energy eigenvalue and wave-function with those obtained from other typical means popular among physics students, namely the numerical shooting method and perturbation theory. Show that the idea of program of the numerical shooting method and perturbation theory of this problem (see Sects. 2.1 and 4) The numerical shooting method is generally regarded as one of the most efficient methods that give very accurate results because it integrates the Schr?dinger equation directly, though in the numerical sense. The n = even case is shown in Figs. 4, 5 and 6. In this case, the wave-function has split up on asymmetric nodes under Gaussian asymmetric double-well potential. The n = odd case is shown in Fig. 7. In this case, the wave-function has not split up on asymmetric nodes under Gaussian asymmetric double-well potential.     

The shifted scheme in the general-model-space diagrammatic perturbation theory

The shifted scheme of many-body perturbation theory is applied to open-shell states within the framework of the general-model-space theory. Rules for shifting the denominators of folded diagrams. which appear in open-shell perturbation expansions, are given. The finite-order energies in the shifted scheme obtained in two equivalent representations may differ. This happens, for instance, in the case of triplet states. For ³Σ_u⁺ states of the He₂, differences up to 0.07 mhartree have been found in third order. A similar phenomenon is the size inconsistency of the shifted scheme observed by Silver in the ground state of He₂. A possible advantage of the shifted scheme is its faster convergence for excited states.     

Bennett's acceptance ratio and histogram analysis methods enhanced by umbrella sampling along a reaction coordinate in configurational space

Free energy perturbation, a method for computing the free energy difference between two states, is often combined with non-Boltzmann biased sampling techniques in order to accelerate the convergence of free energy calculations. Here we present a new extension of the Bennett acceptance ratio (BAR) method by combining it with umbrella sampling (US) along a reaction coordinate in configurational space. In this approach, which we call Bennett acceptance ratio with umbrella sampling (BAR-US), the conditional histogram of energy difference (a mapping of the 3N-dimensional configurational space via a reaction coordinate onto 1D energy difference space) is weighted for marginalization with the associated population density along a reaction coordinate computed by US. This procedure produces marginal histograms of energy difference, from forward and backward simulations, with higher overlap in energy difference space, rendering free energy difference estimations using BAR statistically more reliable. In addition to BAR-US, two histogram analysis methods, termed Bennett overlapping histograms with US (BOH-US) and Bennett-Hummer (linear) least square with US (BHLS-US), are employed as consistency and convergence checks for free energy difference estimation by BAR-US. The proposed methods (BAR-US, BOH-US, and BHLS-US) are applied to a 1-dimensional asymmetric model potential, as has been used previously to test free energy calculations from non-equilibrium processes. We then consider the more stringent test of a 1-dimensional strongly (but linearly) shifted harmonic oscillator, which exhibits no overlap between two states when sampled using unbiased Brownian dynamics. We find that the efficiency of the proposed methods is enhanced over the original Bennett's methods (BAR, BOH, and BHLS) through fast uniform sampling of energy difference space via US in configurational space. We apply the proposed methods to the calculation of the electrostatic contribution to the absolute solvation free energy (excess chemical potential) of water. We then address the controversial issue of ion selectivity in the K(+) ion channel, KcsA. We have calculated the relative binding affinity of K(+) over Na(+) within a binding site of the KcsA channel for which different, though adjacent, K(+) and Na(+) configurations exist, ideally suited to these US-enhanced methods. Our studies demonstrate that the significant improvements in free energy calculations obtained using the proposed methods can have serious consequences for elucidating biological mechanisms and for the interpretation of experimental data.     

Theoretical development of a Gaussian potential function in the description of the radial portion of isotropic bending vibrations

Analysis of a Gaussian potential function suitable for modeling degenerate bending vibrations in weakly bound molecular complexes is presented. Approximate eigenvalues and eigenvectors are obtained by application of perturbation theory. Comparison to the “exact” eigenvalues obtained via a numerical solution shows that the first- and higher-order perturbation corrections are consistent with variational principles.     

A theoretical study of the excited states of AmO2n+, n=1,2,3

The ground and excited states of the AmO(2) (+), AmO(2) (2+), and AmO(2) (3+) ions have been studied using the four-component configuration interaction singles doubles, spin-orbit complete active space self-consistent field, and spin-orbit complete active space-order perturbation theory methods. The roles of scalar relativistic effects and spin-orbit coupling are analyzed; results with different methods are carefully compared by a precise analysis of the wave functions. A molecular spinor diagram is used in relation to the four-component calculations while a ligand field model is used for the two-step method. States with the same number of electrons in the four nonbonding orbitals are in very good agreement with the two methods while ligand field and charge transfer states do not have the same excitation energies.