Ivar Waller 90 years on June 11, 1988

Linkage properties of the diagrammatic representation of the energies obtained in the multireference many-body perturbation calculations with respect to the incompleteness or completeness of the model space are discussed. The case of not completely degenerate model space is considered for which a comparison with the standard single-reference many-body perturbation expansion is possible. The Hose–Kaldor type of graphical representation of the perturbation expansion for the effective Hamiltonian is used in this comparison. It is shown that for an incomplete model space the perturbation expansion is not size-extensive. In this case, for a truncated expansion of the effective Hamiltonian, the energies obtained by diagonalization of the effective Hamiltonian matrix are represented by both linked and unlinked irreducible contributions. The unlinked ones do not appear when the complete model space is used.     

The viscosity and temperature dependence of X-band ESR lineshapes of Gd(III) aqueous complex

X-band ESR spectra of Gd-aqua complex in various weight concentration of glycerol have been recorded at four temperatures. The interpretation of the ESR linewidth is preformed using both the stochastic Liouville approach (SLA) and a perturbation theory. The SLA uses a one dynamic model of the zero-field splitting whereas the perturbation approach uses a two dynamic model. Both models can reproduce the variation of the linewidth with respect to viscosity. In the SLA model, both the zero-field splitting (ZFS) interaction and the correlation time vary with the glycerol content. In the two dynamic perturbation model, only the correlation times are viscosity dependent. The two models give different NMRD profiles.     

调制BZ-CSTR反应体系小振荡诱导非平衡相变

以流量为扰动参数,对连续搅拌釜式反应器（ＣＳＴＲ）中的Ｂｅｌｏｕｓｏｖ－Ｚｈａｂｏｔｉｎｄｓｋｙ（ＢＺ）反应体系小振荡进行周期或噪声扰动,改变扰动周期或扰动强度,在某些条件下可以诱导非平衡相变。以ＢＺ反应的四变量Ｍｏｎｔａｎａｔｏｒ模型为基础,对相应流率区进行周期或噪声扰动,计算结果与实验基本定性一致。     

A perturbation method for the Ornstein-Zernike equation and the generic van der Waals equation of state for a square well potential model

We calculate the generic van der Waals parameters A and B for a square well model by means of a perturbation theory. To calculate the pair distribution function or the cavity function necessary for the calculation of A and B, we have used the Percus-Yevick integral equation, which is put into an equivalent form by means of the Wiener-Hopf method. This latter method produces a pair of integral equations, which are solved by a perturbation method treating the Mayer function or the well width or the functions in the square well region exterior to the hard core as the perturbation. In the end, the Mayer function times the well width is identified as the perturbation parameter in the present method. In this sense, the present perturbation method is distinct from the existing thermodynamic perturbation theory, which expands the Helmholtz free energy in a perturbation series with the inverse temperature treated as an expansion parameter. The generic van der Waals parameters are explicitly calculated in analytic form as functions of reduced temperature and density. The van der Waals parameters are recovered from them in the limits of vanishing density and high temperature. The equation of state thus obtained is tested against Monte Carlo simulation results and found reliable, provided that the temperature is in the supercritical regime. By scaling the packing fraction with a temperature-dependent hard core, it is suggested to construct an equation of state for fluids with a temperature-dependent hard core that mimicks a soft core repulsive force on the basis of the equation of state derived for the square well model.     

A nondiagonal quasidegenerate fourth-order perturbation theory

A canonical quasidegenerate Rayleigh-Schrödinger perturbation theory, correct through fourth order in the energy, is explored for a block-diagonal unperturbed Hamiltonian. The theory is developed completely within a Lie Algebra in Hilbert space. Explicit equations forn-particle transition elements in terms of solutions of simultaneous linear equations are presented. A two-dimensional anisotropic anharmonic oscillator is used to provide numerical results. The perturbation theory is shown to be stable under small separation of model and complement spaces. An iterative variant of the fourth order perturbation theory is introduced; the iterative variant is related to the non-iterative one in much the same way as nondegenerate coupled cluster theories are related to nondegenerate perturbation theory. The quasidegenerate coupled cluster theory appears to be stable in the presence of multiple intruder states.     

Bond length alternation in cyclic polyenes. V. Local finite-order perturbation theory approach

The finite-order many-body perturbation theory using the localized Wannier orbital basis is applied to the problem of bond length alternation in the Pariser–Parr–Pople model of cyclic polyenes C_N H_N, N = 4v + 2, which may be regarded as a simplified model of polyacetylene. Both the Møller–Plesset and the Epstein–Nesbet-type partitionings of the model Hamiltonian are employed. The localized orbital basis enables an efficient truncation of the perturbation theory summations over the intermediate states as well as an elimination of energetically unimportant diagrams, thus enabling one to obtain the fourth-order Møller–Plesset-perturbation energies with a relatively small computational effort even for large polyenes. The results obtained with the second-, third-, and fourth-order Møller–Plesset and with the third-order Epstein–Nesbet perturbation theories yield very similar bond length distortions (about 0.05 Å) and stabilization energies per site (about 0.04 eV) as obtained earlier with the RHF , one-parameter AMO , and delocalized orbital perturbation theories. The effects of truncation and diagram elimination in the fourth-order Møller–Plesset perturbation theory and the abnormal behavior of the second-order Epstein–Nesbet perturbation theory results in the localized Wannier basis near the instability threshold of the RHF solutions are discussed.     

A Simple Approximate Perturbation Approach to Quasi-degenerate Systems

In this paper a simple approximate approach for the study of quasi-degenerate systems is presented in the frame of multireference perturbation theory. The formulation can be considered as an approximation of the quasi-degenerate perturbation theory (QDPT) with the simplification that only the state specific (diagonal) perturbation corrections to the energy have to be computed. The new approach is discussed and compared with previous QDPT formulations using the weakly avoided crossing model (for which new properties are here presented) and applied to the case of the neutral/ionic energy crossing in the LiF molecule.     

Perturbational approach to the electron correlation cusp applied to helium‐like atoms

A recently proposed perturbational approach to the electron correlation cusp problem 1 is tested in the context of three spherically symmetrical two‐electron systems: helium atom, hydride anion, and a solvable model system. The interelectronic interaction is partitioned into long‐ and short‐range components. The long‐range interaction, lacking the singularities responsible for the electron correlation cusp, is included in the reference Hamiltonian. Accelerated convergence of orbital‐based methods for this smooth reference Hamiltonian is shown by a detailed partial wave analysis. Contracted orbital basis sets constructed from atomic natural orbitals are shown to be significantly better for the new Hamiltonian than standard basis sets of the same size. The short‐range component becomes the perturbation. The low‐order perturbation equations are solved variationally using basis sets of correlated Gaussian geminals. Variational energies and low‐order perturbation wave functions for the model system are shown to be in excellent agreement with highly accurate numerical solutions for that system. Approximations of the reference wave functions, described by fewer basis functions, are tested for use in the perturbation equations and shown to provide significant computational advantages with tolerable loss of accuracy. Lower bounds for the radius of convergence of the resulting perturbation expansions are estimated. The proposed method is capable of achieving sub‐μHartree accuracy for all systems considered here. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Large-Order perturbation theory

The original motivation for studying the asymptotic behavior of the coefficients of perturbation series came from quantum field theory. An overview is given of some of the attempts to understand quantum field theory beyond finite-order perturbation series. At least in the case of the Thirring model and probably in general, the full content of a relativistic quantum field theory cannot be recovered from its perturbation series. This difficulty, however, does not occur in quantum mechanics, and the anharmonic oscillator is used to illustrate the methods used in large-order perturbation theory. Two completely different methods are discussed, the first one using the WKB approximation, and a second one involving the statistical analysis of Feynman diagrams. The first one is well developed and gives detailed information about the desired asymptotic behavior, while the second one is still in its infancy and gives instead information about the distribution of vertices of the Feynman diagrams.     

Simulation of droplet formation and coalescence using lattice Boltzmann-based single-phase model

A lattice Boltzmann method-based single-phase free surface model is developed to study the interfacial dynamics of coalescence, droplet formation and detachment phenomena related to surface tension and wetting effects. Compared with the conventional multiphase models, the lattice Boltzmann-based single-phase model has a higher computational efficiency since it is not necessary to simulate the motion of the gas phase. A perturbation, which is given in the same fashion as the perturbation step in Gunstensen's color model, is added to the distribution functions of the interface cells for incorporating the surface tension into the single-phase model. The assignment of different mass gradients along the fluid-wall interface is used to model the wetting properties of the solid surface. Implementations of the model are demonstrated for simulating the processes of the droplet coalescence, the droplet formation and detachment from ceiling and from nozzles with different shapes and different wall wetting properties.     

SCF perturbation calculations on metalloporphyrins

The influence of metal ion on the oxidation and ionisation potentials of metalloporphyrins is investigated by the simple electrostatic model using SCF perturbation theory. The zero order wavefunctions are obtained from PPP and CNDO/2 methods. The wide variations in redox potentials with metal and the relative insensitivity of the optical transitions with metal are very well accounted for by the perturbation approach.     

Thomas–Fermi theory from a perturbative treatment of the hydrogenic limit energy density functional: A possible link to local density functional theory

A perturbation expansion which connects the hydrogenic limit energy density functional to the Thomas–Fermi functional is discussed. This perturbation series, where the Coulomb energy density functional is treated as the perturbation to the hydrogenic limit functional, is, in fact, the q = (N/Z) expansion of Thomas–Fermi theory. A truncated form of the first-order correction to the functional provides further insight into the model which treats the ground state energy as a local functional of the electron density.     

Second-order perturbational treatment of normal coordinates in the string model for the hydration reaction of formaldehyde

A perturbation theory for normal coordinates of nonadiabatic solvation is presented by means of the “string model” of chemical reactions. The dynamic normal coordinate is introduced for the perturbational treatment of the “intrinsic” normal coordinates that are orthogonal to the reaction path. The reaction is defined as the intrinsic reaction coordinate (IRC ) that is treated as a string. The string is thrown in the external force field that acts as a nonadiabatic source of perturbation. As an application of the present treatment, the effect of a water molecule for hydration reaction of formaldehyde is calculated. A second-order perturbation effect for the enhancement of the reaction rate is found.     

Test of projected perturbation theory on a simplified model of H+2

The perturbation theory for projected states is shown to converge for a schematic one-dimensional model of H⁺₂. Ordinary perturbation theory diverges for the odd state. It converges for the even state, but very marginally for large interactive distance.     

Dynamic behaviors of Monod type chemostat model with impulsive perturbation on the nutrient concentration

In this paper, the dynamic behaviors of a Monod type chemostat model with impulsive perturbation are investigated. Using Floquet theory and small amplitude perturbation method, we prove that the microorganism-eradication periodic solution is asymptotically stable if the impulsive period satisfies some conditions. Moreover, the permanence of the system is discussed in detail. Finally, we verify the main results by numerical simulation.      

On the nature of the Møller-Plesset critical point

It has been suggested [F. H. Stillinger, J. Chem. Phys. 112, 9711 (2000)] that the convergence or divergence of M?ller-Plesset perturbation theory is determined by a critical point at a negative value of the perturbation parameter z at which an electron cluster dissociates from the nuclei. This conjecture is examined using configuration-interaction computations as a function of z and using a quadratic approximant analysis of the high-order perturbation series. Results are presented for the He, Ne, and Ar atoms and the hydrogen fluoride molecule. The original theoretical analysis used the true Hamiltonian without the approximation of a finite basis set. In practice, the singularity structure depends strongly on the choice of basis set. Standard basis sets cannot model dissociation to an electron cluster, but if the basis includes diffuse functions then it can model another critical point corresponding to complete dissociation of all the valence electrons. This point is farther from the origin of the z plane than is the critical point for the electron cluster, but it is still close enough to cause divergence of the perturbation series. For the hydrogen fluoride molecule a critical point is present even without diffuse functions. The basis functions centered on the H atom are far enough from the F atom to model the escape of electrons away from the fluorine end of the molecule. For the Ar atom a critical point for a one-electron ionization, which was not previously predicted, seems to be present at a positive value of the perturbation parameter. Implications of the existence of critical points for quantum-chemical applications are discussed.     

A variational principle for weak stimulated absorptions with energy minimizing fields

Optimal time-dependent perturbations in agreement with the Bohr condition for spectroscopic transitions are obtained variationally for two-level systems via the first-order perturbation model. The criterion of optimality is that specified weak absorptions are achieved with minimal applied energy. Results obtain for fixed perturbation intervals of arbitrary length.     

Heat capacities of Stockmayer fluids from Monte Carlo simulations and perturbation theory

The heat capacities of dipolar fluids are investigated using a thermodynamic perturbation theory approach and the NVT and NpT Monte Carlo simulation methods. The theoretical results are compared to corresponding simulation data. The comparison shows that the applied perturbation theory is appropriate for the heat capacity calculations. As an application, the isobaric heat capacity of ammonia is also studied by the Stockmayer fluid model.     

Hard-sphere perturbation theory for a model of liquid Ga

Investigating thermodynamic properties of a model for liquid Ga, we have extended the application of the hard-sphere (HS) perturbation theory to an interatomic pair potential that possesses a soft repulsive core and a long-range oscillatory part. The model is interesting for displaying a discontinuous jump on the main-peak position of the radial distribution function at some critical density. At densities less than this critical value, the effective HS diameter of the model, estimated by the variational HS perturbation theory, has a substantial reduction with increasing density. Thus, the density dependence of the packing fraction of the HS reference fluid has an anomalous behavior, with a negative slope, within a density region below the critical density. By adding a correction term originally proposed by Mon to remedy the inherent deficiency of the HS perturbation theory, the extended Mansoori-Canfield/Rasaiah-Stell theory [J. Chem. Phys. 120, 4844 (2004)] very accurately predicts the Helmholtz free energy and entropy of the model, including an excess entropy anomaly. Almost occurring in the same density region, the excess entropy anomaly is found to be associated with the anomalous packing faction of the HS fluid.     

Anharmonic analysis of the vibrational spectrum of ketene by density functional theory using second-order perturbative approach

The paper reports main results of a comprehensive study of the vibrational spectrum of ketene computed using second-order perturbation theory treatment based on quartic, cubic and semidiagonal quartic force constants. Two different models--a homogeneous model using the same density functionals and basis functions for the harmonic calculations and anharmonic corrections, and a hybrid model in which the two parts of the calculation are conducted using different density functionals and basis sets--have been employed in the present calculations. Different DFT and CCSD methods and DZ and TZ extended basis sets involving diffuse and polarization functions have been used to calculate optimized and vibrationally averaged geometrical parameters, the harmonic and anharmonic vibrational frequencies and the spectroscopic constants such as anharmonicity constants, rotational constants, rotation-vibration coupling constants, Nielsen's centrifugal distortion constants and Coriolis coupling constants. Homogeneous model is found to be superior to the hybrid model in several respects. Difficulties in the hybrid model may arise due to one of the following reasons: (a) the basic requirement that the geometry optimization and frequency calculations must be done at the same level of theory to have valid frequencies is not met in the hybrid model; (b) the molecular structure gets reoptimized at the low level for anharmonic corrections; (c) in addition, the perturbation could also diverge for the above reasons, particularly for the very low, very anharmonic terms where the harmonic approximation is not close enough to make the perturbation work.