Transformation of Bloch representation into atomic orbital representation and combined matrix element calculation technique for spatially bounded basis functions in crystals

A method is proposed for transforming the Hamiltonian from Bloch to atomic function representation. For spatially bounded functions, this is a rigorous method based on solution of a certain algebraic system of equations. Unlike the conventional procedure based on integration over the Brillouin zone, the new method requires knowledge of the matrix elements of the Bloch representation only at several points of the Brillouin zone. The number of these points is determined by the trimming radius for the spatially bounded functions and by the lattice constant. The method can be used for calculating matrix elements in a basis of atomic functions and for reducing computations in matrix element calculations of the Bloch representation for procedures using numerical integration.     

An optimum minimal basis of exponential atomic orbitals

An analysis of some of the energetic properties of the conventional minimal STO basis is used to suggest a new optimum set of exponential functions for use in molecular calculations.     

Quantum-chemical modeling of the atomic and electronic structure of thin diamond-like nanocrystallites

The structural and electronic characteristics of thin diamond-like nanocrystallites (with cross-sectional areas 5 < S < 280 Ų) were investigated by the electron density functional tight binding (DFTB) method. A new type of extended “hybrid” (sp³ + sp²) nanostructures, in the form of monolithic diamond-like sp³ crystallites encapsulated in a graphite-like sp² shell, was discovered. __________ Translated from Teoreticheskaya i éksperimental’naya Khimiya, Vol. 42, No. 4, pp. 199–203, July–August, 2006.     

Intrinsic local constituents of molecular electronic wave functions. I. Exact representation of the density matrix in terms of chemically deformed and oriented atomic minimal basis set orbitals

A coherent, intrinsic, basis-set-independent analysis is developed for the invariants of the first-order density matrix of an accurate molecular electronic wavefunction. From the hierarchical ordering of the natural orbitals, the zeroth-order orbital space is deduced, which generates the zeroth-order wavefunction, typically an MCSCF function in the full valence space. It is shown that intrinsically embedded in such wavefunctions are elements that are local in bond regions and elements that are local in atomic regions. Basis-set-independent methods are given that extract and exhibit the intrinsic bond orbitals and the intrinsic minimal-basis quasi-atomic orbitals in terms of which the wavefunction can be exactly constructed. The quasi-atomic orbitals are furthermore oriented by a basis-set independent method (viz. maximization of the sum of the fourth powers of all off-diagonal density matrix elements) so as to exhibit clearly the chemical interactions. The unbiased nature of the method allows for the adaptation of the localized and directed orbitals to changing geometries. Contribution of the Mark S. Gordon 65th Birthday Fegtschrift Issue.     

Optimal diffuse augmented atomic basis sets for extrapolation of the correlation energy

We seek correlation-consistent diffuse-augmented double-zeta and triple-zeta basis sets that perform optimally in extrapolating the correlation energy to the one-electron complete basis set limit, denoted oAVXZ and oAV(X + d)Z. The novel basis sets are method-dependent in that they are trained to perform optimally for the correlation energy at each specific level of theory. They are shown to yield accurate results in calculating both the energy and tensorial properties such as polarizabilities while not significantly altering the Hartree-Fock energy. Quantitatively, complete basis set limit (CBS)-/(oAVdZ,oAVtZ)-extrapolated correlation energies typically outperform, by 3- to 5-fold, the ones calculated with traditional ansatzes of similar flexibility. Attaining energies of CBS/(AVtZ,AVqZ) type or better accuracy, they frequently outperform expensive raw explicitly correlated ones. Promisingly, a limited test on CBS-extrapolated energies based on conventional basis sets has shown that they compare well even with extrapolated explicitly correlated ones. Calculated atomization and dissociation energies, molecular geometries, ionization potentials, and electron affinities also tend to outperform the ones obtained with traditional Dunning's ansatzes from which the new basis sets have been determined. The method for basis set generation is simple, and there is no reason of principle why the approach could not be adapted for handling other bases in the literature.     

Implementation of atomic basis set composed of 1s Gaussian and 1s Slater-type orbitals to carry out quantum mechanics molecular calculations

A mixed atomic basis set formed with ls Slater-type orbitals and 1s floating spherical Gaussian orbitals is implemented. Evaluation of multicenter integrals is carried out using a method based on expansion of binary products of atomic basis functions in terms of a complete basis set, and a systematic analysis is performed. The proposed algorithm is very stable and furnishes fairly good results for total energy and geometry. An LCAO-SCF test calculation is carried out on LiH. The trends observed show that there are some combinations of mixed orbitals that are appropriate to describe the system. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 604–609, 1999     

Electron localizability indicators ELI and ELIA: the case of highly correlated wavefunctions for the argon atom

Electron localizability indicators based on the same-spin electron pair density and the opposite-spin electron pair density are studied for correlated wavefunctions of the argon atom. Different basis sets and reference spaces are used for the multireference configuration interaction method following the complete active space calculations aiming at the understanding of the effect of local electron correlation when approaching the exact wavefunction. The populations of the three atomic shells of Ar atom in real space are calculated for each case.     

Real‐space grid representation of momentum and kinetic energy operators for electronic structure calculations

We show that the central finite difference formula for the first and the second derivative of a function can be derived, in the context of quantum mechanics, as matrix elements of the momentum and kinetic energy operators on discrete coordinate eigenkets defined on a uniform grid. Starting from the discretization of integrals involving canonical commutations, simple closed‐form expressions of the matrix elements are obtained. A detailed analysis of the convergence toward the continuum limit with respect to both the grid spacing and the derivative approximation order is presented. It is shown that the convergence from below of the eigenvalues in electronic structure calculations is an intrinsic feature of the finite difference method. © 2018 Wiley Periodicals, Inc.     

Big data analysis of ab Initio molecular integrals in the neglect of diatomic differential overlap approximation

Most modern semiempirical quantum-chemical (SQC) methods are based on the neglect of diatomic differential overlap (NDDO) approximation to ab initio molecular integrals. Here, we check the validity of this approximation by computing all relevant integrals for 32 typical organic molecules using Gaussian-type orbitals and various basis sets (from valence-only minimal to all-electron triple-ζ basis sets) covering in total more than 15.6 million one-electron (1-e) and 10.3 billion two-electron (2-e) integrals. The integrals are calculated in the nonorthogonal atomic basis and then transformed by symmetric orthogonalization to the Löwdin basis. In the case of the 1-e integrals, we find strong orthogonalization effects that need to be included in SQC models, for example, by strategies such as those adopted in the available OMx methods. For the valence-only minimal basis, we confirm that the 2-e Coulomb integrals in the Löwdin basis are quantitatively close to their counterparts in the atomic basis and that the 2-e exchange integrals can be safely neglected in line with the NDDO approximation. For larger all-electron basis sets, there are strong multishell orthogonalization effects that lead to more irregular patterns in the transformed 2-e integrals and thus cast doubt on the validity of the NDDO approximation for extended basis sets. Focusing on the valence-only minimal basis, we find that some of the NDDO-neglected integrals are reduced but remain sizable after the transformation to the Löwdin basis; this is true for the two-center 2-e hybrid integrals, the three-center 1-e nuclear attraction integrals, and the corresponding three-center 2-e hybrid integrals. We consider a scheme with a valence-only minimal basis that includes such terms as a possible strategy to go beyond the NDDO integral approximation in attempts to improve SQC methods. © 2018 Wiley Periodicals, Inc.     

Variational problem for ground and excited rovibronic states on the basis of a molecular Hamiltonian in the principal central axes

A general scheme of realization of nonempirical methods for calculation of the ground and excited rovibronic states of flexible molecules is proposed. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

推、拉电子基团对水杨醛缩胺分子电子结构和非线性光学性质的影响

采用量子化学半经验FF／PM3方法,讨论了水杨醛缩苯胺分子中两苯环的对位被推、拉电子基团取代后,体系电子结构和非线性光学性质的变化,考查了分子电子结构对非线性光学性质影响的微观本质,得到的具有给体－共轭桥键－受体型结构的水杨醛缩苯胺分子应显示良好的非线光学性质。     

An atomic charge study of highly ionic diatomic molecular systems

Four atomic charge formalisms are compared using highly ionic diatomic molecules, such as LiF, NaF, KF, LiCl, NaCl, KCl, BF, AlF, GaF, BeO, and MgO. All calculations were done at the QCISD/6‐311G(2df) level. The only formalism consistent with the characteristics of all these systems is Quantum theory of atoms in molecules (QTAIM). Absolute Mulliken charge values are small. ChelpG charges are not reliable for systems in which the atoms are largely anisotropic. Generalized atomic polar tensor values are contaminated with charge fluxes and atomic dipole fluxes and fail when these contributions are important and do not cancel each other. Finally, the charge–charge flux–dipole flux model was applied to dipole moment derivatives with QTAIM. This analysis shows that charge flux and atomic dipole flux contributions during bond stretching are almost null, except for oxides. There are also evidences that the lone electron pair at Group 13 elements in fluorides becomes less localized as the bond is stretched. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Radial basis neural network for the classification of fresh edible oils using an electronic nose

An electronic nose utilising an array of six-bulk acoustic wave polymer coated Piezoelectric Quartz (PZQ) sensors has been developed. The nose was presented with 346 samples of fresh edible oil headspace volatiles, generated at 45°C. Extra virgin olive (EVO), Non-virgin olive oil (OI) and Sunflower oil (SFO), were used over a period of 30 days. The sensor responses were then analysed producing an architecture for the Radial Basis Function Artificial Neural Network (RBF). It was found that the RBF results were excellent, giving classifications of above 99% for the vegetable oil test samples.This revised version was published online in November 2005 with corrections to the Cover Date.     

Accuracy and efficiency of atomic basis set methods versus plane wave calculations with ultrasoft pseudopotentials for DNA base molecules

Recent results from Preuss et al. (J Comput Chem 2004, 25, 112) on DNA base molecules, obtained by plane wave density functional calculations using ultrasoft pseudopotentials, are compared with calculations using Gaussian basis sets. Bond lengths and angles agree closely, but dihedral angles and vibrational frequencies show significant differences. The Gaussian basis calculations are at least an order of magnitude more efficient than the plane wave/ultrasoft pseudopotential calculations at a similar level of accuracy; the advantage is even larger if the Fourier Transform Coulomb method is used. To obtain definite benchmark values, the geometries of the four DNA bases were optimized at the MP2 level with large basis sets, up to cc-pVQZ and aug-cc-pVTZ.     

Contracted well-tempered Gaussian basis sets in SCF calculations on the ground and excited electronic states of neutral and ionized diatomic molecules containing first-row atoms

Summary Calculations were done on ground and excited states of C₂, C ₂ ⁺ , C ₂ ^– , N₂, N ₂ ⁺ , O₂, O ₂ ⁺ , O ₂ ^– , CO, CO⁺, CO²⁺, and CO^– using contracted well-tempered basis sets. The (14s 10p) basis sets were augmented with threed, one or twof, and oneg functions. Total energies, orbital energies, and spectroscopic constants were compared with the best available computational data.     

Calculation of wave‐functions with frozen orbitals in mixed quantum mechanics/molecular mechanics methods. II. Application of the local basis equation

The application of the local basis equation (Ferenczy and Adams, J. Chem. Phys. 2009 , 130, 134108) in mixed quantum mechanics/molecular mechanics (QM/MM) and quantum mechanics/quantum mechanics (QM/QM) methods is investigated. This equation is suitable to derive local basis nonorthogonal orbitals that minimize the energy of the system and it exhibits good convergence properties in a self‐consistent field solution. These features make the equation appropriate to be used in mixed QM/MM and QM/QM methods to optimize orbitals in the field of frozen localized orbitals connecting the subsystems. Calculations performed for several properties in divers systems show that the method is robust with various choices of the frozen orbitals and frontier atom properties. With appropriate basis set assignment, it gives results equivalent with those of a related approach [G. G. Ferenczy previous paper in this issue] using the Huzinaga equation. Thus, the local basis equation can be used in mixed QM/MM methods with small size quantum subsystems to calculate properties in good agreement with reference Hartree–Fock–Roothaan results. It is shown that bond charges are not necessary when the local basis equation is applied, although they are required for the self‐consistent field solution of the Huzinaga equation based method. Conversely, the deformation of the wave‐function near to the boundary is observed without bond charges and this has a significant effect on deprotonation energies but a less pronounced effect when the total charge of the system is conserved. The local basis equation can also be used to define a two layer quantum system with nonorthogonal localized orbitals surrounding the central delocalized quantum subsystem. © 2013 Wiley Periodicals, Inc.     

Non-empirical molecular orbital theory of the electronic structure of molecular crystals

A non-empirical molecular orbital treatment of molecular crystals, based on SCF perturbation theory and matrix partitioning methods is presented.     

Calculation of electrostatic interaction energies in molecular dimers from atomic multipole moments obtained by different methods of electron density partitioning

Accurate and fast evaluation of electrostatic interactions in molecular systems is still one of the most challenging tasks in the rapidly advancing field of macromolecular chemistry, including molecular recognition, protein modeling and drug design. One of the most convenient and accurate approaches is based on a Buckingham-type approximation that uses the multipole moment expansion of molecular/atomic charge distributions. In the mid-1980s it was shown that the pseudoatom model commonly used in experimental X-ray charge density studies can be easily combined with the Buckingham-type approach for calculation of electrostatic interactions, plus atom-atom potentials for evaluation of the total interaction energies in molecular systems. While many such studies have been reported, little attention has been paid to the accuracy of evaluation of the purely electrostatic interactions as errors may be absorbed in the semiempirical atom-atom potentials that have to be used to account for exchange repulsion and dispersion forces. This study is aimed at the evaluation of the accuracy of the calculation of electrostatic interaction energies with the Buckingham approach. To eliminate experimental uncertainties, the atomic moments are based on theoretical single-molecule electron densities calculated at various levels of theory. The electrostatic interaction energies for a total of 11 dimers of alpha-glycine, N-acetylglycine and L-(+)-lactic acid structures calculated according to Buckingham with pseudoatom, stockholder and atoms-in-molecules moments are compared with those evaluated with the Morokuma-Ziegler energy decomposition scheme. For alpha-glycine a comparison with direct "pixel-by-pixel" integration method, recently developed Gavezzotti, is also made. It is found that the theoretical pseudoatom moments combined with the Buckingham model do predict the correct relative electrostatic interactions energies, although the absolute interaction energies are underestimated in some cases. The good agreement between electrostatic interaction energies computed with Morokuma-Ziegler partitioning, Gavezzotti's method, and the Buckingham approach with atoms-in-molecules moments demonstrates that reliable and accurate evaluation of electrostatic interactions in molecular systems of considerable complexity is now feasible.