Energy gradients with respect to atomic positions and cell parameters for the Kohn-Sham density-functional theory at the Gamma point

The application of theoretical methods based on density-functional theory is known to provide atomic and cell parameters in very good agreement with experimental values. Recently, construction of the exact Hartree-Fock exchange gradients with respect to atomic positions and cell parameters within the Gamma-point approximation has been introduced. In this article, the formalism is extended to the evaluation of analytical Gamma-point density-functional atomic and cell gradients. The infinite Coulomb summation is solved with an effective periodic summation of multipole tensors. While the evaluation of Coulomb and exchange-correlation gradients with respect to atomic positions are similar to those in the gas phase limit, the gradients with respect to cell parameters needs to be treated with some care. The derivative of the periodic multipole interaction tensor needs to be carefully handled in both direct and reciprocal space and the exchange-correlation energy derivative leads to a surface term that has its origin in derivatives of the integration limits that depend on the cell. As an illustration, the analytical gradients have been used in conjunction with the QUICCA algorithm to optimize one-dimensional and three-dimensional periodic systems at the density-functional theory and hybrid Hartree-Fock/density-functional theory levels. We also report the full relaxation of forsterite supercells at the B3LYP level of theory.     

A new energy decomposition scheme for molecular interactions within the Hartree-Fock approximation

A new method is proposed for the analysis of components of molecular interaction energy within the Hartree-Fock approximation. The Hartree-Fock molecular orbitals of the isolated molecules are used as the basis for the construction of Fock matrix of the supermolecule. Then certain blocks of this matrix are set to zero subject to specify boundary conditions of the supermolecule molecular orbitals, and the resultant matrix is diagonalized iteratively to obtain the desired energy components. This method can be considered as an extension of our previous method, but has an advantage in the explicit definition of the charge transfer energy, placing it on an equal footing with the exchange and polarization terms. The new method is compared with existing perturbation methods, and is also applied to the energy and electron density decomposition of (H₂O)₂.     

Hartree-Fock exchange computed using the atomic resolution of the identity approximation

In this work, we apply the atomic resolution of the identity (ARI) fitting approximation to the computation of Hartree-Fock exchange. The ARI approximation is a local modification of the RI approximation that produces an energy which is differentiable with respect to nuclear motion, unlike other local applications of RI. We justify empirically the use of locality and present timing comparisons of ARI, RI, and exact computation for one-, two-, and three-dimensional carbon systems. ARI is found to reduce significantly the cost of RI for large systems, while retaining accuracy.     

Minimization of the energy of a molecule with respect to nonlinear parameters of the basis set

We study characteristic features of minimization of the Hartree-Fock-Roothaan energy with respect to nonlinear parameters of the Gaussian basis set. We describe and apply regularization of the discrete Newton-Raphson method based on the analysis of eigenvalues of the Hessian matrix. We discuss results of groundstate energy calculations for the molecules LiH, CH⁺, CH^–, He₃ ²⁺, BH₂ ⁺, and H₂O in optimal ls-Gaussian basis sets. We find that, for molecules with four to six electrons, good accuracy is obtained with small basis sets consisting of ls-functions only.Translated from Teoreticheskaya i Éxperimental'naya Khimiya, Vol. 24, No. 2, pp. 215–218, March–April, 1988.     

Density approximation to the average Hartree-Fock exchange potential for atoms

A non-local density-based approximation to the average Hartree-Fock (HF) exchange potential is developed. The new potential is formulated within the spin-dependent version of the weighted density approximation, and is based on a novel form of the (exchangeonly) pair-correlation function for electrons in finite systems. The results for total energies and one-electron orbital energies of atoms are reasonably accurate in comparison with those obtained using the exact average HF potential, or the exact orbital-dependent HF potential.     

Analytical Hartree–Fock gradients with respect to the cell parameter for systems periodic in three dimensions

Analytical Hartree–Fock gradients with respect to the cell parameter have been implemented in the electronic structure code CRYSTAL, for the case of three-dimensional periodicity. The code is based on Gaussian-type orbitals, and the summation of the Coulomb energy is performed with the Ewald method. It is shown that a high accuracy of the cell gradient can be achieved.     

Analytical Hartree–Fock gradients with respect to the cell parameter: systems periodic in one and two dimensions

Analytical Hartree–Fock gradients with respect to the cell parameter have been implemented in the electronic structure code CRYSTAL, for the case of one- and two-dimensional periodicity. As in most molecular codes, Gaussian type orbitals are used to express the wavefunction. Examples demonstrate that the gradients have a good.     

Basis set convergence studies of Hartree-Fock calculations of molecular properties within the resolution of the identity approximation

Within the resolution of the identity (RI) method, the convergence of the Hartree-Fock (HF) total molecular energy and the multipole moments in the course of the combined regular expansion of the molecular and auxiliary (RI) basis sets is studied. Dunning's cc-pVXZ series is used for both the molecular and the RI basis sets. The results show the calculated quantities converge to the HF limit when both the molecular and the RI basis sets are expanded from correlation-consistent polarized valence double zeta to correlation-consistent polarized valence sextuple zeta. Combinations of molecular/RI basis sets sufficient for convergence of the total energy and of the multipole moments at various accuracy levels have been determined. A measure of the RI basis set incompleteness is suggested and discussed. As it is significantly faster than the standard HF algorithm for small and midsize molecules, the RI-HF method, together with appropriate expanding series of both molecular and RI basis sets, provide an efficient tool to estimate and control the error of the Hartree-Fock calculations due to the finite basis set.     

Linear-scaling method for calculating nuclear magnetic resonance chemical shifts using gauge-including atomic orbitals within Hartree-Fock and density-functional theory

Details of a new density matrix-based formulation for calculating nuclear magnetic resonance chemical shifts at both Hartree-Fock and density functional theory levels are presented. For systems with a nonvanishing highest occupied molecular orbital-lowest unoccupied molecular orbital gap, the method allows us to reduce the asymptotic scaling order of the computational effort from cubic to linear, so that molecular systems with 1000 and more atoms can be tackled with today's computers. The key feature is a reformulation of the coupled-perturbed self-consistent field (CPSCF) theory in terms of the one-particle density matrix (D-CPSCF), which avoids entirely the use of canonical MOs. By means of a direct solution for the required perturbed density matrices and the adaptation of linear-scaling integral contraction schemes, the overall scaling of the computational effort is reduced to linear. A particular focus of our formulation is to ensure numerical stability when sparse-algebra routines are used to obtain an overall linear-scaling behavior.     

Some statistic parameters of the Schulz distribution with respect to the molecular heterogeneity

Summary The concepts of absolute and relative heterogeneity related to the number average degree of polymerization were formulated.In case ofSchulz distribution it was proved that the relative heterogeneity depends only on the curve parameterk.Different statistic parameters for theSchulz distribution were calculated. The dependence of those parameters on the curve parameterk is illustrated by diagrams.

---

Zusammenfassung Die Begriffe von absoluter und relativer Uneinheitlichkeit in bezug auf das Zahlenmittel des Polymerisationsgrades wurden formuliert.Es wurde festgestellt, daß die relative Uneinheitlichkeit im Falle derSchulzschen Verteilung nur von dem Kurvenparameterk abhängig ist.Verschiedene statistische Parameter für dieSchulzsche Verteilung werden berechnet und die Abhängigkeit dieser Parameter von dem Kurvenparameterk durch Diagramme illustriert.

     

Basis set convergence of molecular correlation energy differences within the random phase approximation

The basis set convergence of energy differences obtained from the random phase approximation (RPA) to the correlation energy is investigated for a wide range of molecular interactions. For dispersion bound systems the basis set incompleteness error is most pronounced, as shown for the S22 benchmark [P. Jurecka et al., Phys. Chem. Chem. Phys. 8, 1985 (2006)]. The use of very large basis sets (> quintuple-zeta) or extrapolation to the complete basis set (CBS) limit is necessary to obtain a reliable estimate of the binding energy for these systems. Counterpoise corrected results converge to the same CBS limit, but counterpoise correction without extrapolation is insufficient. Core-valence correlations do not play a significant role. For medium- and short-range correlation, quadruple-zeta results are essentially converged, as demonstrated for relative alkane conformer energies, reaction energies dominated by intramolecular dispersion, isomerization energies, and reaction energies of small organic molecules. Except for weakly bound systems, diffuse augmentation almost universally slows down basis set convergence. For most RPA applications, quadruple-zeta valence basis sets offer a good balance between accuracy and efficiency.     

On the invariance of the configuration interaction energy with respect to orbital rotations

Summary The invariance of the configuration interaction (CI) energy with respect to orbital rotation is considered. The inclusion of all spin couplings versus only those from the first-order interacting space is considered. A definition for the analog of a second-order CI calculation when inactive electrons are present is proposed.     

Resolution of identity approach for the Kohn-Sham correlation energy within the exact-exchange random-phase approximation

Two related methods to calculate the Kohn-Sham correlation energy within the framework of the adiabatic-connection fluctuation-dissipation theorem are presented. The required coupling-strength-dependent density-density response functions are calculated within exact-exchange time-dependent density-functional theory, i.e., within time-dependent density-functional response theory using the full frequency-dependent exchange kernel in addition to the Coulomb kernel. The resulting resolution-of-identity exact-exchange random-phase approximation (RI-EXXRPA) methods in contrast to previous EXXRPA methods employ an auxiliary basis set (RI basis set) to improve the computational efficiency, in particular, to reduce the formal scaling of the computational effort with respect to the system size N from N(6) to N(5). Moreover, the presented RI-EXXRPA methods, in contrast to previous ones, do not treat products of occupied times unoccupied orbitals as if they were linearly independent. Finally, terms neglected in previous EXXRPA methods can be included, which leads to a method designated RI-EXXRPA+, while the method without these extra terms is simply referred to as RI-EXXRPA. Both EXXRPA methods are shown to yield total energies, reaction energies of small molecules, and binding energies of noncovalently bonded dimers of a quality that is similar and in some cases even better than that obtained with quantum chemistry methods such as Mo?ller-Plesset perturbation theory of second order (MP2) or with the coupled cluster singles doubles method. In contrast to MP2 and to conventional density-functional methods, the presented RI-EXXRPA methods are able to treat static correlation.     

A multireference energy decomposition scheme with respect to fragment valence states

In this paper we describe a multi-reference energy decomposition scheme defined with respect to specific valence states of the interacting fragments. The orbitals of these fragments are computed with a new iterative procedure involving MC-SCF computations with the fragments at infinite separation and CI computations at the internuclear distance of interest. This procedure is applied here for the analysis of the rotational barrier in methylamine.