The W algorithm and the D transformation for the numerical evaluation of three‐center nuclear attraction integrals

Three‐center nuclear attraction integrals over exponential‐type functions are required for ab initio molecular structure calculations and density functional theory (DFT). These integrals occur in many millions of terms, even for small molecules, and they require rapid and accurate numerical evaluation. The use of a basis set of B functions to represent atomic orbitals, combined with the Fourier transform method, led to the development of analytic expressions for these molecular integrals. Unfortunately, the numerical evaluation of the analytic expressions obtained turned out to be extremely difficult due to the presence of two‐dimensional integral representations, involving spherical Bessel integral functions. % The present work concerns the development of an extremely accurate and rapid algorithm for the numerical evaluation of these spherical Bessel integrals. This algorithm, which is based on the nonlinear D transformation and the W algorithm of Sidi, can be computed recursively, allowing the control of the degree of accuracy. Numerical analysis tests were performed to further improve the efficiency of our algorithm. The numerical results section demonstrates the efficiency of this new algorithm for the numerical evaluation of three‐center nuclear attraction integrals. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Thermodynamic excess properties of alkanol + amine mixtures and application of the ERAS model

New experimental data of the molar excess volume V ^E of the mixtures ethanol/n-butylamine, heptanol/n-butylamine, n-propanol/dibutylamine have been obtained using the technique of the vibrating tube densitometer. Together with the data for the molar excess enthalpy H ^E from the literature, the V ^E data have been used for testing the applicability of the socalled ERAS model which accounts for hydrogen bonding effects as well as for free volume effects in associating mixtures. The results obtained by adjusting the model parameters reveal a strong cross association between the unlike molecules in the mixture resulting from strong negative values for the hydrogen bonding energy and the hydrogen bonding volume.Communicated at the Festsymposium celebrating Dr. Henry V. Kehiaian's 60^th birthday, Clermont-Ferrand, France, 17–18 May 1990.     

Application of coset representations to the construction of symmetry adapted functions

Summary A coset representation (G(/G _i )), which is defined algebraically by a coset decomposition of a finite groupG by its subgroupG _i , is shown to be a method for the decomposition of a regular body into its point group orbits. This proof also shows that each member of theG(/G _i ) orbit belongs to theG _i site-symmetry. In addition, a general equation concerning the multiplicities of such coset representations is derived and shown to involve Brester's equations and thek-value equations of framework groups as special cases. The relationship of the coset representation and the site-symmetry affords a general procedure for obtaining symmetry adapted functions.     

A modified steepest descent method with applications to maximizing likelihood functions

In maximizing a non-linear function G(), it is well known that the steepest descent method has a slow convergence rate. Here we propose a systematic procedure to obtain a 1–1 transformation on the variables , so that in the space of the transformed variables, the steepest descent method produces the solution faster. The final solution in the original space is obtained by taking the inverse transformation. We apply the procedure in maximizing the likelihood functions of some generalized distributions which are widely used in modeling count data. It was shown that for these distributions, the steepest descent method via transformations produced the solutions very fast. It is also observed that the proposed procedure can be used to expedite the convergence rate of the first derivative based algorithms, such as Polak-Ribiere, Fletcher and Reeves conjugate gradient methods as well.     

Lipschitzian inverse functions,directional derivatives,and applications inC 1,1 optimization

The paper shows that Thibault's limit sets allow an iff-characterization of local Lipschitzian invertibility in finite dimension. We consider these sets as directional derivatives and extend the calculus in a way that can be used to clarify whether critical points are strongly stable inC ^1,1 optimization problems.Many fruitful discussions with colleagues D. Klatte and K. Tammer as well as with H. Th. Jongen and F. Nozicka have influenced the present investigations in a very constructive manner. For the original papers concerning the sets f(x; u), the author is indebted to Prof. L. Thibault.     

Metastable systems driven by colored noise: The stationary state

The influence of the Bardeen-Herring back-jump correlations on the Fermi-Dirac statistics of the one-dimensional nonhomogeneous fermionic lattice gas is studied by the Monte Carlo simulation technique and semianalytically. The resulting distribution is obtained, exhibiting increased population of the lower levels in comparison to the Fermi-Dirac statistics.     

ε-optimality criteria for convex programming problems via exact penalty functions

In this paper, we present-optimality criteria for convex programming problems associated with exact penalty functions. Several authors have given various criteria under the assumption that such convex problems and the associated dual problems can be solved. We assume the solvability of neither the convex problem nor the dual problem. To derive our criteria, we estimate the size of the penalty parameter in terms of an-solution for the dual problem.     

Nonequilibrium fluctuations in μ space

The properties of fluctuations in space in or outside thermal equilibrium are obtained by solving hierarchies of equations derived either from the Liouville or the Master equation. In particular we study the one-, two-, etc., time correlation functions that describe the spatial and temporal behavior of the fluctuations in space. Explicit solutions are obtained for a dilute gas. The Langevin approach is briefly discussed. Our results are compared with those obtained in the extensive literature, which is reviewed in some detail.     

Must thermodynamic functions be piecewise analytic?

There seems to be no general theoretical argument to support the idea that thermodynamic functions are piecewise analytic. We suggest that nonanalyticity may be associated with Gibbs states which are quasiperiodic under space translations, or have a more general nonperiodic (turbulent) behavior.     

Variational approximations for renormalization group transformations

Approximate recursion relations which give upper and lower bounds on the free energy are described. Optimal calculations of the free energy can then be obtained by treating parameters within the renormalization equations variationally. As an example, a particularly simple lower bound approximation which preserves the symmetry of the Hamiltonian (the one-hypercube approximation) is described. The approximation is applied to both the Ising model and the Wilson-Fisher model. At the fixed point a parameter is set variationally and critical indices are calculated. For the Ising model the agreement with the exact results atd = 2 is surprisingly good, 0.1%, and is good atd=3 and evend=4. For the Wilson-Fisher model the recursion relation is reduced to a one-dimensional integral equation which can be solved numerically givingv=0.652 atd=3, or by expansion in agreement with the results of Wilson and Fisher to leading order in . The method is also used to calculate thermodynamic functions for thed = 2 Ising model; excellent agreement with the Onsager solution is found.Supported in part by the National Science Foundation under Grants Nos. MPS73-04886A01 and GH-41512 and by the Brown University Materials Research Laboratory supported by the National Science Foundation. M.C.Y. was supported by a grant from the Scientific and Technical Research Council of Turkey.