Ab initio molecular orbital calculations on large lattice cluster models: Use of translational symmetry

Summary Translational symmetry has been shown to be useful in the calculation of electronic structures of large lattice models. The number of unique integrals has been derived for cases of different dimensionality. For the unique integrals zero screening and approximation methods are described. The method has been applied to arrays of hydrogen atoms and to a zincblende surface model. When the size of the system is increased the translationally unique integrals are shown to become either zero or they can be calculated by simple coulombic approximations.     

Using ordinary differential equations system to solve isoconversional problems in non-isothermal kinetic analysis

The mathematical evaluation of the activation energy, E, of non-isothermal degradation reactions is usually made using the Ozawa/Flynn–Wall isoconversion principle and involves the numerical resolution of a set of integrals without closed form solution, which are solved by polynomial approximation or by numeric integration. In the present work, the isoconversion principle, originally described and maintained until now as an algebraic problem, was written as a set of ordinary differential equations (ODEs). The individual ODEs obtained are integrated by numeric methods and are used to estimate the activation energy of simulated examples. A least square error (LSE) objective function using the introduced ODEs was written to deal with multiple heating rate CaCO₃ thermal decomposition TG experiments.     

Exploiting non-abelian point group symmetry in direct two-electron integral transformations

Summary A simple and general scheme to exploit any discrete point group symmetry in two-electron integral transformations is introduced. It has been implemented together with integral pre-screening techniques in direct two-electron integral transformations. Its application has also been extended to subsequent MO integral processing steps like MP2 or solution of the coupled-perturbed Hartree-Fock equations (CPHF). Timings for representative applications are presented.     

On the complete solution of ∈y″=y 3

In Ref. 1, the author claimed that the problem y=y ³ is soluble only for a certain range of the parameter . An analytic approach, as adopted in the following contribution, reveals that a unique solution exists for any positive value of . The solution is given in closed form by means of Jacobian elliptic functions, which can be numerically computed very efficiently. In the limit 0+, the solutions exhibit boundary-layer behavior at both endpoints. An easily interpretable approximate solution for small is obtained using a three-variable approach.     

Estimating the complexity of a class of path-following methods for solving linear programs by curvature integrals

In this paper we study a particular class of primal-dual path-following methods which try to follow a trajectory of interior feasible solutions in primal-dual space toward an optimal solution to the primal and dual problem. The methods investigated are so-called first-order methods: each iteration consists of a long step along the tangent of the trajectory, followed by explicit recentering steps to get close to the trajectory again. It is shown that the complexity of these methods, which can be measured by the number of points close to the trajectory which have to be computed in order to achieve a desired gain in accuracy, is bounded by an integral along the trajectory. The integrand is a suitably weighted measure of the second derivative of the trajectory with respect to a distinguished path parameter, so the integral may be loosely called a curvature integral.     

A new algorithm to evaluate bielectronic integrals with 1s Slater‐type orbitals obtained by the integral transforms

This article presents a variation of the integral transform method to evaluate multicenter bielectronic integrals (12|34), with 1s Slater‐type orbitals. It is proved that it is possible to define, out of the expression of (12|34) given by the integral transform method, a function F(q) that has the property of having a unique Q, such that F(Q) = (12|34). Therefore, F(q) may be used to calculate (12|34). It is shown that the evaluation of F(Q) turns out to be simpler than the three‐dimensional integral involved in the calculation of (12|34), and an algorithm is presented to calculate Q. The results show that relative errors on the order of 10^?3 or lower are obtained very efficiently. In addition, it is shown that the proposed algorithm is very stable. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Molecular integrals over Gaussian-type geminal basis functions

We present formulas for the evaluation of molecular integrals over basis functions with an explicit Gaussian dependence on interelectronic coordinates. These formulas use expansions in Hermite Gaussian functions and represent an extension to the work of McMurchie and Davidson to two-electron basis functions. Integrals that depend on the coordinates of up to four electrons are discussed explicitly. A key feature of this approach is that it allows full exploitation of the shell structure of the orbital part of the basis. Received: 24 February 1997 / Accepted: 4 March 1997     

A semiempirical model for the two-center repulsion integrals in the NDDO approximation

A self-consistent formalism is proposed for the two-center electron repulsion integrals in the NDDO approximation, based on their expansion in terms of multipole-multipole interactions and free from adjustable parameters.     

Quasi‐random integration in quantum chemistry: Efficiency,stability, and application to the study of small atoms and molecules constrained in spherical boxes

This project consists of two parts. In the first part, a series of test calculations is performed to verify that the integrals involved in the determination of atomic and molecular properties by standard self‐consistent field (SCF) methods can be obtained through Halton, Korobov, or Hammersley quasi‐random integration procedures. Through these calculations, we confirm that all three methods lead to results that meet the levels of precision required for their use in the calculation of properties of small atoms or molecules at least at a Hartree–Fock level. Moreover, we have ensured that the efficiency of quasi‐random integration methods that we have tested is Halton=Korobov>Hammersley?pseudo‐random. We also find that these results are comparable to those yielded by ordinary Monte Carlo (pseudo‐random) integration, with a calculation effort of two orders of smaller magnitude. The second part, which would not have been possible without the integration method previously analyzed, contains a first study of atoms constrained in spherical boxes through SCF calculations with basis functions adapted to the features of the problem: Slater‐type orbitals (STOs) trimmed by multiplying them by a function that yields 1 for 0 < r < (R‐δ), polynomial values for (R‐δ) < r < R and null for r > R, R being the radius of the box and δ a variationally determined interval. As a result, we obtain a equation of state for electrons of small systems, valid just in the limit of low temperatures, but fairly simple. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Evaluation of Multicenter Electronic Attraction, Electric Field and Electric Field Gradient Integrals with Screened and Nonscreened Coulomb Potentials over Integer and Noninteger n Slater Orbitals

By the use of complete orthonormal sets of -exponential-type orbitals, where ( = 1, 0, –1, –2,...) the multicenter electronic attraction (EA), electric field (EF) and electric field gradient (EFG) integrals of nonscreened and Yukawa-like screened Coulomb potentials are expressed through the two-center overlap integrals with the same screening constants and the auxiliary functions introduced in our previous paper (I.I. Guseinov, J. Phys. B, 3 (1970) 1399). The recurrence relations for auxiliary functions are useful for the calculation of multicenter EA, EF and EFG integrals for arbitrary integer and noninteger values of principal quantum numbers, screening constants, and location of slater-type orbitals. The convergence of the series is tested by calculating concrete cases.