共查询到20条相似文献,搜索用时 0 毫秒
1.
HONG Gong-Yi LI Le-MinCollege of Chemistry Molecular Engineering State Key Laboratory of Rare Earth Materials Chemistry Applications Peking University Beijing China 《中国化学》1996,14(4):289-296
Several numerical integration schemes for the evaluation of matrix elements in density functional theory calculations have been studied and compared by computational practice. The best scheme was found to be the combination of the atomic partition function proposed by Becke with the scaled generalized Gauss-Laguerre quadrature formula for radial integration suggested by Yang, which achieve the highest convergence rate to the numerical integration. With the same number of integration points, the accuracy of the calculated results by this scheme is higher by 1 to 2 orders of magnitudes than that by other schemes. The reason for achieving higher accuracy by this scheme has been proposed preliminarily. 相似文献
2.
We introduce a Gaussian quadrature, based on the polynomials that are orthogonal with respect to the weight function ln(2)x on the interval [0, 1], which is suitable for the evaluation of radial integrals. The quadrature is exact if the non-Jacobian part of the integrand is a linear combination of a geometric sequence of exponential functions. We find that the new scheme is a useful alternative to existing approaches, particularly for integrands that exhibit multiexponential behavior. 相似文献
3.
Recently, Gill and Chien introduced a new radial quadrature for multiexponential integrands (MultiExp grid) to deal with the radial part of the numerical integration. In this article, the MultiExp grid is studied and used to integrate the charge density. The MultiExp grid, along with an optimal pruning scheme, performed very well both in terms of accuracy and efficiency compared to other radial mappings commonly used in Density Functional Theory. 相似文献
4.
Hideaki Ishikawa Kazuo Yamamoto Kazumi Fujima Misako Iwasawa 《International journal of quantum chemistry》1999,72(5):509-523
A powerful and accurate numerical three‐dimensional integration scheme was developed especially for molecular orbital calculations. A multicenter integral is decomposed into the sum of single‐center integrals using nuclear weight functions and calculated using Gaussian quadrature rules. The decomposed single‐center integrands show strong anisotropy. With a careful selection of the Gaussian quadrature rule according to the anisotropy, it is possible to obtain an accuracy of 13 digits with a small number of integration points for the overlap integrals, normalization integrals, and molecular integrals for the hydrogen molecule. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 72: 509–523, 1999 相似文献
5.
Mirko Franchini Pierre Herman Theodoor Philipsen Lucas Visscher 《Journal of computational chemistry》2013,34(21):1819-1827
In this article, we document a new implementation of the fuzzy cells scheme for numerical integration in polyatomic systems [Becke, J. Chem. Phys. 1998, 88, 2547] and compare its efficiency and accuracy with respect to an integration scheme based on the Voronoi space partitioning. We show that the accuracy of both methods is comparable, but that the fuzzy cells scheme is better suited for geometry optimization. For this method, we also introduce the locally dense grid concept and present a proof‐of‐concept application. © 2013 Wiley Periodicals, Inc. 相似文献
6.
Rodríguez JI Köster AM Ayers PW Santos-Valle A Vela A Merino G 《Journal of computational chemistry》2009,30(7):1082-1092
We introduce a method to compute atomic properties according to the \"quantum theory of atoms in molecules.\" An integration grid in real space is partitioned into subsets, omega(i). The subset, omega(i), is composed of all grid points contained in the atomic basin, Omega(i), so that integration over Omega(i) is reduced to simple quadrature over the points in omega(i). The partition is constructed from deMon2k's atomic center grids by following the steepest ascent path of the density starting from each point in the grid. We also introduce a technique that exploits the cellular nature of the grid to make the algorithm faster. The performance of the method is tested by computing properties of atoms and nonnuclear attractors (energies, charges, dipole, and quadrupole moments) for a set of representative molecules. 相似文献
7.
Geudtner G Janetzko F Köster AM Vela A Calaminici P 《Journal of computational chemistry》2006,27(4):483-490
The parallelization of the LCGTO-KS-DFT code deMon2k is presented. The parallelization of the three-center electron repulsion integrals, the numerical integration using a direct grid algorithm and the matrix multiplication and diagonalization are described. The efficiency of the parallelization is analyzed by selected benchmark calculations. It is shown that geometry optimizations of systems with more than 8,000 basis functions are feasible on cluster architectures. 相似文献
8.
P. Kolandaivel N. Suba K. Senthilkumar 《International journal of quantum chemistry》2000,76(5):662-669
Balint Kurti's Fourier grid Hamiltonian method is employed to obtain the molecular wave function and equilibrium bond length for H2 and HF molecules. The density functional theory parameter, namely, the chemical hardness (η) value, was determined for some diatomic hydride molecules using this wave function and the results are found to be in good agreement with the values obtained from the ab initio HF–SCF method. A new formula for chemical hardness (η=1/2Dr, where D is the proportionality constant and r is the internuclear distance) is introduced in binding energy and change of hardness equations to determine the chemical hardness and chemical potential values for different bond lengths. The binding energy and change of hardness values are calculated for H2, H, H, HF, HF+, and HF− molecules and the bond stability is discussed. Finally, the concept of an atom in a molecule is examined in the context of DFT parameters and comparison is made between an atom in a molecule and the isolated atom. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 76: 662–669, 2000 相似文献
9.
Dr. Scott Collins Dr. Anuj Joshi Prof. Mikko Linnolahti 《Chemistry (Weinheim an der Bergstrasse, Germany)》2021,27(62):15460-15471
Methylaluminoxane (MAO) activators have sheet structures which form ion-pairs on reaction of neutral donors such as octamethyltrisiloxane (OMTS). The ion-pairs can be detected by electrospray ionization mass spectrometry (ESI-MS) in polar media. The growth of these reactive precursors during hydrolysis of Me3Al can be monitored using ESI-MS. Density functional theory, combined with numerical simulation of growth, indicates that this process involves rapid formation of low MW oligomers, followed by assembly of these species into low MW sheets. These can grow through further addition of low MW oligomers or by fusion into larger sheets. The mechanism of these growth processes leads to the prediction that even-numbered sheets should be favored, and this surprising result is confirmed by ESI-MS monitoring experiments of both activator growth and MAO aging. 相似文献
10.
Alves-Santos M Dávila LY Petrilli HM Capaz RB Caldas MJ 《Journal of computational chemistry》2006,27(2):217-227
We present a detailed analysis of the application of density functional theory (DFT) methods to the study of structural properties of molecular and supramolecular systems, using as a paradigmatic example three para-phenylene-based systems: isolated biphenyl, single chain poly-para-phenylene, and crystalline biphenyl. We use different functionals for the exchange correlation potential, the local density (LDA), and generalized gradient approximations (GGA), and also different basis sets expansions, localized, plane waves (PW), and mixed (localized plus PW), within the reciprocal space formulation for the hamiltonian. We find that regardless of the choice of basis functions, the GGA calculations yield larger interring distances and torsion angles than LDA. For the same XC approximation, the agreement between calculations with different basis functions lies within 1% (LDA) or 0.5% (GGA) for distances, and while PW and mixed basis calculations agree within 1 degrees for torsion angles, the localized basis results show larger angles by approximately 8 degrees and a nonmonotonic dependence on basis size, with differences within 6 degrees. The most prominent features, namely the torsion between rings for isolated molecule and infinite chain, and planarity for the molecule in crystalline environment, are well reproduced by all DFT calculations. 相似文献
11.
Bottom-up nanoparticle (NP) formation is assumed to begin with the reduction of the precursor metallic ions to form zero-valent atoms. Studies in which this assumption was made are reviewed. The standard reduction potential for the formation of aqueous metallic atoms—E0(Mn+aq/M0aq)—is significantly lower than the usual standard reduction potential for reducing metallic ions Mn+ in aqueous solution to a metal in solid state. E0(Mn+aq/M0solid). E0(Mn+aq/M0aq) values are negative for many typical metals, including Ag and Au, for which E0(Mn+aq/M0solid) is positive. Therefore, many common moderate reduction agents that do not have significantly high negative reduction standard potentials (e.g., hydrogen, carbon monoxide, citrate, hydroxylamine, formaldehyde, ascorbate, squartic acid, and BH4−), and cannot reduce the metallic cations to zero-valent atoms, indicating that the mechanism of NP production should be reconsidered. Both AgNP and AuNP formations were found to be multi-step processes that begin with the formation of clusters constructed from a skeleton of M+-M+ (M = Ag or Au) bonds that is followed by the reduction of a cation M+ in the cluster to M0, to form Mn0 via the formation of NPs. The plausibility of M+-M+ formation is reviewed. Studies that suggest a revised mechanism for the formation of AgNPs and AuNPs are also reviewed. 相似文献
12.
13.
Pomelli CS 《Journal of computational chemistry》2004,25(12):1532-1541
A new integration grid for the Polarizable Continuum Model reaction field is presented. In this scheme (called TsLess) the molecular surface is not partitioned in tesserae, but a set of sampling points is chosen. These points are weighted according to they position on the surface. The new procedure is tested against the GEPOL tesselation method. 相似文献
14.
《Journal of computational chemistry》2017,38(12):869-882
Density‐functional approximations developed in the past decade necessitate the use of quadrature grids that are far more dense than those required to integrate older generations of functionals. This category of difficult‐to‐integrate functionals includes meta‐generalized gradient approximations, which depend on orbital gradients and/or the Laplacian of the density, as well as functionals based on B97 and the popular “Minnesota” class of functionals, each of which contain complicated and/or oscillatory expressions for the exchange inhomogeneity factor. Following a strategy introduced previously by Gill and co‐workers to develop the relatively sparse “SG‐0” and “SG‐1” standard quadrature grids, we introduce two higher‐quality grids that we designate SG‐2 and SG‐3, obtained by systematically “pruning” medium‐ and high‐quality atom‐centered grids. The pruning procedure affords computational speedups approaching a factor of two for hybrid functionals applied to systems of atoms, without significant loss of accuracy. The grid dependence of several popular density functionals is characterized for various properties. © 2017 Wiley Periodicals, Inc. 相似文献
15.
V. M. Tapilin 《Journal of Structural Chemistry》2003,44(4):531-537
A numerical integration method is suggested for calculating Hamiltonian matrix elements in a basis of functions with polynomial tails with allowance for discontinuities of higher-order derivatives of the basis function within the domain of integration. The method is tested by calculating matrix elements for a copper crystal. The results for the overlap matrix elements are presented demonstrating efficiency of the method. 相似文献
16.
Zbigniew Romanowski 《International journal of quantum chemistry》2008,108(2):249-256
The efficient algorithm calculating the overlap and the kinetic integrals for the numerical atomic orbitals is presented. The described algorithm exploits the properties of the prolate spheroidal coordinates. The overlap and the kinetic integrals in ?3 are reduced to the integrals over the rectangular domain in ?2, what substantially reduces the complexity of the problem. We prove that the integrand over the rectangular domain is continuous and does not have any slope singularities. For calculation of the integral over the rectangle any adaptive algorithm can be applied. The exemplary results were obtained by application of the adaptive Gauss quadrature. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008 相似文献
17.
Erin Day Brice Kauffmann Matthieu Scarpi-Luttenauer Alain Chaumont Marc Henry Pierre Mobian 《Chemistry (Weinheim an der Bergstrasse, Germany)》2022,28(19):e202200047
The formation of a tetranuclear self-assembled species constructed around a TiO4N2 motif is reported. This aggregate is generated from Ti(OiPr)4, 2,2’-bipyrimidine (bpym) and a bis-biphenol strand (L2H4) where two 2,2’-biphenol units are connected with a biphenyl spacer. The solid-state structure of the [Ti4(L2)4(bpym)4] architecture reveals the formation of an unprecedented chiral alternate [2×2] grid. In addition to the structural characterization of the [Ti4(L2)4(bpym)4] architecture, geometry optimisation on various possible isomeric tetrameric assemblies ([2×2] grid, alternate [2×2] grid, circular helicate or cyclic hemihelicate) is performed using DFT calculations. These results confirm the higher stability of the alternate [2×2] grid over the other possible tetranuclear isomers and allow examining the replacement of the bpym ligands by two novel diimine chelates within the tetranuclear assembly (2,2’-bipyridine=bipy and 2,2’-bipyrazine=bipyraz). From this initial theoretical investigation, the competition between these three nitrogen ligands in the course of the self-assembly process is next evaluated. Overall, this investigation shows that the exclusive formation of the alternate [2×2] grid is driven by CH⋅⋅⋅N interactions. 相似文献
18.
We present a novel algorithm to integrate property densities over the volume of a quantum topological atom. Atoms are grown outward, starting from a sphere centered on the nucleus, by means of a finite element meshing algorithm. Bond critical points and ring critical points require special treatment. The overall philosophy as well as intricate features of this meshing algorithm are given, followed by details of the quadrature over the finite elements. An effort has been made to design a streamlined and compact algorithm, focusing on the core of challenges arising in tracing the electron density's gradient vector field. The current algorithm also generates a new type of pictures that can be a Graphical User Interface. Excellent integration errors, L(Omega), are obtained, even for atoms with (narrow) tails or sharp corners. 相似文献
19.
离子选择性电极标准加入法的测量误差研究 总被引:1,自引:0,他引:1
研究了离子选择性电极标准加入法的测量误差与标准溶液加入量的关系。从理论上导出了离子选择性电极一次标准加入法测量结果的相对误差计算公式。当加入标准溶液使试液增加的浓度等于试液原来的浓度,测量结果的相对误差最小。 相似文献