共查询到20条相似文献,搜索用时 0 毫秒
1.
Patrick Merlot Thomas Kjærgaard Trygve Helgaker Roland Lindh Francesco Aquilante Simen Reine Thomas Bondo Pedersen 《Journal of computational chemistry》2013,34(17):1486-1496
An analysis of Dunlap's robust fitting approach reveals that the resulting two‐electron integral matrix is not manifestly positive semidefinite when local fitting domains or non‐Coulomb fitting metrics are used. We present a highly local approximate method for evaluating four‐center two‐electron integrals based on the resolution‐of‐the‐identity (RI) approximation and apply it to the construction of the Coulomb and exchange contributions to the Fock matrix. In this pair‐atomic resolution‐of‐the‐identity (PARI) approach, atomic‐orbital (AO) products are expanded in auxiliary functions centered on the two atoms associated with each product. Numerical tests indicate that in 1% or less of all Hartree–Fock and Kohn–Sham calculations, the indefinite integral matrix causes nonconvergence in the self‐consistent‐field iterations. In these cases, the two‐electron contribution to the total energy becomes negative, meaning that the electronic interaction is effectively attractive, and the total energy is dramatically lower than that obtained with exact integrals. In the vast majority of our test cases, however, the indefiniteness does not interfere with convergence. The total energy accuracy is comparable to that of the standard Coulomb‐metric RI method. The speed‐up compared with conventional algorithms is similar to the RI method for Coulomb contributions; exchange contributions are accelerated by a factor of up to eight with a triple‐zeta quality basis set. A positive semidefinite integral matrix is recovered within PARI by introducing local auxiliary basis functions spanning the full AO product space, as may be achieved by using Cholesky‐decomposition techniques. Local completion, however, slows down the algorithm to a level comparable with or below conventional calculations. © 2013 Wiley Periodicals, Inc. 相似文献
2.
Philippe Blanchard José M. Gracia‐Bondiacute;a Joseph C. Várilly 《International journal of quantum chemistry》2012,112(4):1134-1164
Forty‐five years after the point de départ [Hohenberg and Kohn, Phys Rev, 1964, 136, B864] of density functional theory, its applications in chemistry and the study of electronic structures keep steadily growing. However, the precise form of the energy functional in terms of the electron density still eludes us—and possibly will do so forever [Schuch and Verstraete, Nat Phys, 2009, 5, 732]. In what follows we examine a formulation in the same spirit with phase space variables. The validity of Hohenberg–Kohn–Levy‐type theorems on phase space is recalled. We study the representability problem for reduced Wigner functions, and proceed to analyze properties of the new functional. Along the way, new results on states in the phase space formalism of quantum mechanics are established. Natural Wigner orbital theory is developed in depth, with the final aim of constructing accurate correlation‐exchange functionals on phase space. A new proof of the overbinding property of the Müller functional is given. This exact theory supplies its home at long last to that illustrious ancestor, the Thomas–Fermi model. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012 相似文献
3.
Biao Du Changchun Yuan Tianzi Yu Li Yang Yang Yang Prof. Dr. Bo Liu Prof. Dr. Song Qin 《Chemistry (Weinheim an der Bergstrasse, Germany)》2014,20(9):2613-2622
In this article, we describe our efforts on the total synthesis of bolivianine ( 1 ) and isobolivianine ( 2 ), involving the synthesis of onoseriolide ( 3 ). The first generation synthesis of bolivianine was completed in 21 steps by following a chiral resolution strategy. Based on the potential biogenetic relationship between bolivianine ( 1 ), onoseriolide ( 3 ), and β‐(E)‐ocimene ( 8 ), the second generation synthesis of bolivianine was biomimetically achieved from commercially available (+)‐verbenone in 14 steps. The improved total synthesis features an unprecedented palladium‐catalyzed intramolecular cyclopropanation through an allylic metal carbene, for the construction of the ABC tricyclic system, and a Diels–Alder/intramolecular hetero‐Diels–Alder (DA/IMHDA) cascade for installation of the EFG tricyclic skeleton with the correct stereochemistry. Transformation from bolivianine to isobolivianine was facilitated in the presence of acid. The biosynthetic mechanism and the excellent regio‐ and endo selectivities in the cascade are well supported by theoretical chemistry based on the DFT calculations. 相似文献
4.
Roman I. Zubatyuk Leonid Gorb Oleg V. Shishkin Mo Qasim Jerzy Leszczynski 《Journal of computational chemistry》2010,31(1):144-150
Performance of the set of density functional approaches for calculation of one‐electron reduction potentials of nitroaromatic compounds was investigated. To select the most precise and affordable method, we selected a set of model molecules and investigated effects of basis set, density functional, and solvation model on the calculation of reduction potentials. It was found that the mPWB1K/TZVP method provides the most accurate gas phase electron affinity values (RMS error is 0.1 eV). This method in conjunction with the PCM (Bondi) method yields also the most accurate difference in solvation energies of neutral oxidized form and anion‐radical reduced form. The final E0 values were calculated with RMS error of 0.10 V, compared with experimental values. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010 相似文献
5.
Paul L. A. Popelier James Burke Nathaniel O. J. Malcolm 《International journal of quantum chemistry》2003,92(3):326-336
The valence charge concentration shell, as determined by the Laplacian of the electron density, is used as a source of quantum topological graphs, called L‐graphs. A considerable number of such graphs are extracted from the ab initio wave functions of 31 molecules calculated at the B3LYP/6‐311+G(2d,p)//B3LYP/6‐311+G(2d,p) level, covering common functional groups in organic chemistry. We show how L‐graphs can be constructed from a largely transferable subgraph called atomic L‐graph. We investigate the topological stability of the L‐graphs as a function of the basis set. Reliable and consistent atomic L‐graphs are only obtained with basis sets of triple‐zeta quality or higher. The recurrence of invariant motifs or subgraphs in the L‐graphs enables the isolation of 16 atomic L‐graphs. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003 相似文献
6.
Younes Valadbeigi Hossein Farrokhpour 《International journal of quantum chemistry》2013,113(12):1717-1721
The proton affinities, gas phase basicities and ionization energies of formic acid, acetic acid, propanoic acid, 2‐propenoic acid, propiolic acid, butanoic acid, 2‐butenoic acid, 3‐botenoic acid, 2‐methyl‐propanoic acid and 2‐methyl‐2‐propenoic acid were calculated using the computational methods including B3LYP/6‐311++G(2df,p), CBS‐Q and G4MP2. Also, the considered properties were calculated using W1BD method only for formic and acetic acids. In addition, the electron affinities of the acids were calculated using B3LYP, CBS‐Q, G4MP2 and G2MP2 methods, separately. The calculations showed that the PA and gas phase basicity increase with the increase in the number of carbon atoms. The calculated Ionization energies of the unsaturated carboxylic acids are less than the corresponding saturated acids, which are in good agreement with the experimental results. © 2013 Wiley Periodicals, Inc. 相似文献
7.
Shubin Liu 《International journal of quantum chemistry》2006,106(8):1762-1768
Based on the Kato's cusp condition of the electron density and our recent relations for local strongly decaying properties in an electronic system, necessary conditions for trial electron densities of atomic and molecular systems are derived. These conditions take the form of integral‐differential equations, and their validity is verified numerically. The relevance of these conditions to the Thomas–Fermi problem in the orbital‐less density functional approach is discussed. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006 相似文献
8.
The distribution of valence and difference densities in crystalline CaMg(CO3)2 was calculated within the framework of the local density functional theory. It is shown that the terminal maxima of difference density located beyond the oxygen atom nuclei have different values due to the polarizing influence of the cations. 相似文献
9.
D. B. Chesnut 《Journal of computational chemistry》2001,22(14):1702-1711
Electron localization function (ELF) theory is used to characterize changes that occur upon excitation from ground singlet to first excited triplet states in a series of isoelectronic 16‐electron molecules including H2CCH2, HNCH2, H2CO, HNNH, HNO, and O2 (ground triplet to excited singlet). ELF allows one to visualize lone pair or nonbonding electrons, and in these cases the π→π* or n→π excitation processes involved lead to an effective 90° rotation of the electronic structure about one heavy atom center and consequent distortion towards pyramidal symmetry about both heavy atom centers. The heavy atom bond lengths change very little in those cases where effectively two‐center three‐electron bonds can be formed (HNNH, HNO, and O2) while a significant lengthening occurs in those cases where hydrogen atoms prevent such interactions (H2CCH2, HNCH2, and H2CO). It is shown that both ELF basin populations and atoms‐in‐molecules (AIM) delocalization indices reflect expected bond orders for conventional single and double bonds provided one compares the ratio of the molecular quantities rather than their absolute magnitudes. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1702–1711, 2001 相似文献
10.
11.
Thermodynamic stability of metal–aminoacid complexes in water is discussed in terms of the Gibbs free energy of water–ligand exchange processes, and the electronic stabilizing factors thoroughly investigated by means of 1‐electron and 2‐electron density properties. Hexacoordinated complexes formed between iron cations and glycine molecules acting as monodentate or bidentate ligands have been chosen as targets for the current study. Results agree with experimental findings, and complexes formed with bidentate ligands are found to be more stable than those formed with monodentate ones. The larger the number of the coordinated glycine molecules the more stable is the complex. Fe(III) complexes are more stable than Fe(II) ones, but differences are small and the Fe3+/Fe2+ exchange process appears to be energetically feasible for these complexes. Formation of the second glycine–iron interaction involving the amino nitrogen in the bidentate ligands is enthalpycally unfavorable but takes place due to the large entropy rise of the process. The larger stability of Fe(III) complexes is due however to the balance between energetic and solvation terms, which is favorable to these complexes. Electron density properties account satisfactorily for the electronic energy changes along the complex formation in terms of ligand–metal electron transfer and covalent bond orders. © 2010 Wiley Periodicals, Inc. J Comput Chem 2010 相似文献
12.
13.
Yu. M. Basalaev Yu. N. Zhuravlev V. S. Permina A. S. Poplavnoi 《Journal of Structural Chemistry》2007,48(6):996-1000
The band structure spectra, densities of states, and valence and difference densities of LiPN2 and NaPN2 crystals were obtained by DFT self-consistent calculations using the nonlocal pseudopotentials and the localized pseudoorbital basis. Crystal-chemical analysis of these compounds shows that they occupy an intermediate position between the ideal structures of β-cristobalite and chalcopyrite, which manifests itself in the peculiarities of the electronic structure and chemical bonding. The valence band consists of three allowed subbands and differs radically from the typical valence band of chalcopyrite crystals in both subband structure and contributions of the s, p, and d atomic orbitals to the crystal orbitals. 相似文献
14.
Juan I. Rodríguez 《Journal of computational chemistry》2013,34(8):681-686
An efficient method for computing the quantum theory of atoms in molecules (QTAIM) topology of the electron density (or other scalar field) is presented. A modified Newton–Raphson algorithm was implemented for finding the critical points (CP) of the electron density. Bond paths were constructed with the second‐order Runge–Kutta method. Vectorization of the present algorithm makes it to scale linearly with the system size. The parallel efficiency decreases with the number of processors (from 70% to 50%) with an average of 54%. The accuracy and performance of the method are demonstrated by computing the QTAIM topology of the electron density of a series of representative molecules. Our results show that our algorithm might allow to apply QTAIM analysis to large systems (carbon nanotubes, polymers, fullerenes) considered unreachable until now. © 2012 Wiley Periodicals, Inc. 相似文献
15.
《Journal of computational chemistry》2017,38(28):2378-2387
In this contribution, we extend our framework for analyzing and visualizing correlated many‐electron dynamics to non‐variational, highly scalable electronic structure method. Specifically, an explicitly time‐dependent electronic wave packet is written as a linear combination of N‐electron wave functions at the configuration interaction singles (CIS) level, which are obtained from a reference time‐dependent density functional theory (TDDFT) calculation. The procedure is implemented in the open‐source Python program det CI@ORBKIT, which extends the capabilities of our recently published post‐processing toolbox (Hermann et al., J. Comput. Chem. 2016, 37, 1511). From the output of standard quantum chemistry packages using atom‐centered Gaussian‐type basis functions, the framework exploits the multideterminental structure of the hybrid TDDFT/CIS wave packet to compute fundamental one‐electron quantities such as difference electronic densities, transient electronic flux densities, and transition dipole moments. The hybrid scheme is benchmarked against wave function data for the laser‐driven state selective excitation in LiH. It is shown that all features of the electron dynamics are in good quantitative agreement with the higher‐level method provided a judicious choice of functional is made. Broadband excitation of a medium‐sized organic chromophore further demonstrates the scalability of the method. In addition, the time‐dependent flux densities unravel the mechanistic details of the simulated charge migration process at a glance. © 2017 Wiley Periodicals, Inc. 相似文献
16.
Various isomers of macropolyhedral borane ions [B20H18]n (n = 0, -2, -4) are investigated by using the density functional theory methods at RB3LYP/6-31+G* and RB3LYP/6-31G* levels to obtain the optimized geometries, harmonic vibrational frequencies, electron structures, and the stability order. The calculated results show that optimized bond lengths are consistent with the available experimental values and the natural populations, taking [a2 -B20H18]4- (4) as an example, are also in agreement with NMR spectra. The calculated vibrational frequencies are all real, so all of these isomers could be stable, among which [a2 -B20H18]2- (3) and [a2 -B20H18]0 (7) are considered for the first time in this paper. On the basis of the contour maps of molecular orbitals, the delocalized characteristic of molecular orbitals and the possible redox mechanism of these ions are also discussed. Moreover, the analysis on counting of skeletal bonding electrons shows that the isomers (1)-(6) obey the electronic requirement predicted by the mno rule, whereas the newly predicted isomer (7) does not match the mno rule. 相似文献
17.
The geometry optimization in delocalized internal coordinates is discussed within the framework of the density functional theory program deMon. A new algorithm for the selection of primitive coordinates according to their contribution to the nonredundant coordinate space is presented. With this new selection algorithm the excessive increase in computational time and the deterioration of the performance of the geometry optimization for floppy molecules and systems with high average coordination numbers is avoided. A new step selection based on the Cartesian geometry change is introduced. It combines the trust radius and line search method. The structure of the new geometry optimizer is described. The influence of the SCF convergence criteria and the grid accuracy on the geometry optimization are discussed. A performance analysis of the new geometry optimizer using different start Hessian matrices, basis sets and grid accuracies is given. 相似文献
18.
Wei‐Xing Xu 《International journal of quantum chemistry》1999,74(1):1-6
Based on the charge density functional theory, a simple method is proposed to calculate the antiparallel electron correlation coefficient, correlation charge, and correlation energy. In contrast to the result stated in the literature, our analysis reveals that the Coulomb hole has a fine structure. To our knowledge, this is the first time it has been proven that the electrons obey the curve correlation law in an intermediate separation range but the electrons obey the linear correlation law in small and large separation ranges. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 74: 1–6, 1999 相似文献
19.
K. Senthilkumar M. Ramaswamy P. Kolandaivel 《International journal of quantum chemistry》2001,81(1):4-10
Neal's procedure has been applied to determine the electron density ρ(x) for the H2 molecule. The chemical hardness has been calculated employing the ab initio and density functional theory methods and the values are found to be reasonably good. The principle of maximum hardness (PMH) was tested. Fukui functions and the distribution of electron density around the internuclear distance were studied employing the electron density of the H2 molecule. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 4–10, 2001 相似文献
20.
Based on the analysis of the general properties for the one‐ and two‐particle reduced density matrices, a new natural orbital functional is obtained. It is shown that by partitioning the two‐particle reduced density matrix in an antisymmeterized product of one‐particle reduced density matrices and a correction Γc we can derive a corrected Hartree–Fock theory. The spin structure of the correction term from the improved Bardeen–Cooper–Schrieffer theory is considered to take into account the correlation between pairs of electrons with antiparallel spins. The analysis affords a nonidempotent condition for the one‐particle reduced density matrix. Test calculations of the correlation energy and the dipole moment of several molecules in the ground state demonstrate the reliability of the formalism. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 94: 317–323, 2003 相似文献