首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
    
We present a thorough derivation of the mathematical foundations of the representation of the molecular linear electronic density-density response function in terms of a computationally highly efficient moment expansion. Our new representation avoids the necessities of computing and storing numerous eigenfunctions of the response kernel by means of a considerable dimensionality reduction about from 103 to 101. As the scheme is applicable to any compact, self-adjoint, and positive definite linear operator, we present a general formulation, which can be transferred to other applications with little effort. We also present an explicit application, which illustrates the actual procedure for applying the moment expansion of the linear density-density response function to a water molecule that is subject to a varying external perturbation potential. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.  相似文献   

2.
    
We present a low rank moment expansion of the linear density‐density response function. The general interacting (fully nonlocal) density‐density response function is calculated by means of its spectral decomposition via an iterative Lanczos diagonalization technique within linear density functional perturbation theory. We derive a unitary transformation in the space of the eigenfunctions yielding subspaces with well‐defined moments. This transformation generates the irreducible representations of the density‐density response function with respect to rotations within SO(3). This allows to separate the contributions to the electronic response density from different multipole moments of the perturbation. Our representation maximally condenses the physically relevant information of the density‐density response function required for intermolecular interactions, yielding a considerable reduction in dimensionality. We illustrate the performance and accuracy of our scheme by computing the electronic response density of a water molecule to a complex interaction potential. © 2015 Wiley Periodicals, Inc.  相似文献   

3.
    
We have tested the reference interaction site model (RISM) for the case of the hypernetted chain (HNC) and the partially linearized hypernetted chain (PLHNC) closures improved by a repulsive bridge correction (RBC) for ionic hydrated species. We have analyzed the efficiency of the RISM/HNC+RBC and RISM/PLHNC+RBC techniques for decomposition of the electrostatic and the nonpolar hydration energies on the energetic and the enthalpic parts for polyatomic ions when the repulsive bridge correction is treated as a thermodynamic perturbation, and investigate the repulsive bridge effect on the electrostatic potential induced by solvent on solute atoms. For a number of univalent and bivalent atomic ions, molecular cations, and anions, the method provides hydration energies deviating only by several percents from the experimental data. In most cases, the enthalpic contributions to the free energies are also close to the experimental results. The above models are able to satisfactory predict the hydration energies as well as the electrostatic potential around the ionic species. For univalent atomic ions, they also provide qualitative estimates of the Samoilov activation energies.  相似文献   

4.
提出了一个改进的密度泛函理论模型.该模型同时考虑了流体中两分子和三分子相互作用对体系Helm hlotz自由能的贡献,对Ar,N2,CH4,CO2等四种流体的气液界面张力进行了预测,结果与实验值均吻合良好.通过与仅考虑两分子作用的理论计算值和分子模拟值进行比较,表明流体中的三分子相互作用对描述非均相流体的结构与性质有重要影响.  相似文献   

5.
    
Multiconfigurational ab initio methods predict that the 5Πu state as the ground state instead of the 7Δu state. Although multiconfigurational perturbation theory correctly predicts the ground state, they overestimate the bond dissociation energy (BDE). Only multireference configuration interaction method can reasonably calculate the BDE. The spin‐orbit effect on the spectroscopic constants is not significant. The results calculated by density functional theory (DFT) vary significantly depending on the selection of a DFT functional. No DFT functional gives the same energy ordering as calculated by the second‐order multiconfigurational perturbation theory (CASPT2). The old generalized gradient approximations functionals are well suited for predicting the ground state and calculating the bond length and the vibrational frequency of Os2. According to the CASPT2 calculation, the ground state of Os2 has a quadruple bond. © 2014 Wiley Periodicals, Inc.  相似文献   

6.
张义权  罗成林  余智 《无机化学学报》2002,18(11):1147-1150
Magnetic coupling constants J between Cu ions and spin density map for a binuclear complex Cu2(2,2′-bpy)2(4,4′-bpy)2L2 have been calculated by the combination of the broken symmetry approach with the spin project method under the DFT framework and the effect of nonorthogonality between magnetic α and β orbitals in the broken symmetry solution. The result indicates that using SV/SVP basis sets in the calculation of Gaussian 98 is much better than using LanL2DZ pseudo-potential function in above system. It is interesting that calculated J value is quite near to the experimental value in the studies of magnetic susceptibility.  相似文献   

7.
通过对桥联双核铁 (III) [Cl3FeOFeCl3] 2 - 的磁耦合常数的计算 ,探讨了密度泛函理论计算条件对计算结果的影响 .基于密度泛函理论下的破损态方法 ,着重讨论了双核Fe(III) 2 的d5 d5 电子通过氧桥的超交换作用 .研究发现分子的反铁磁通道主要是Fe(III)的dyz和dz2 与 μ O的p轨道形成的 ,具有π /π 和σ /σ 特征的超交换通道  相似文献   

8.
    
We present a new implementation of a recent open‐ended response theory formulation for time‐ and perturbation‐dependent basis sets (Thorvaldsen et al., J. Chem. Phys. 2008, 129, 214108) at the Hartree–Fock and density functional levels of theory. A novel feature of the new implementation is the use of recursive programming techniques, making it possible to write highly compact code for the analytic calculation of any response property at any valid choice of rule for the order of perturbation at which to include perturbed density matrices. The formalism is expressed in terms of the density matrix in the atomic orbital basis, allowing the recursive scheme presented here to be used in linear‐scaling formulations of response theory as well as with two‐ and four‐component relativistic wave functions. To demonstrate the new code, we present calculations of the third geometrical derivatives of the frequency‐dependent second hyperpolarizability for HSOH at the Hartree–Fock level of theory, a seventh‐order energy derivative involving basis sets that are both time and perturbation dependent. © 2014 Wiley Periodicals, Inc.  相似文献   

9.
    
In this work, a graphical tool for postdeprotonation analysis of molecular systems is presented. This tool was applied to a series of molecules to distinguish the information given by its resulting graphics. The outcome identified the h function sensitivity toward electron density rearrangement, being able to recognize bond cleavage and increments in π electron population through simple visual analysis.  相似文献   

10.
In these years there was considerable interest inunderstanding of intermolecular forces in energetic(explosive) systems[1—3]. The supermolecular approach(SM) is widely adopted for calculating ab initio in-termolecular interactions. Nevertheless, it is unable toprovide physically meaningful interaction contribu-tions such as electrostatic, induction, repulsion anddispersion energies. In contrast, the symmetry-adaptedperturbation theory (SAPT)[4—8] has the ability to de-rive these correlated…  相似文献   

11.
    
There are various experimental studies regarding the toxicity and the time of action of local anesthetics, which contain general insights about their pharmacological and physicochemical properties. Although a detailed microscopic analysis of the local anesthetics would contribute to understanding these properties, there are relatively few theoretical studies about these molecules. In this article, we present the results from calculations performed for three local anesthetics: tetracaine, procaine, and lidocaine, both in their charged and uncharged forms, in aqueous environment. We have used the density functional theory and molecular dynamics simulations to study the structural characteristics of these compounds. The radial distribution function g(r) was used to examine the structure of water molecules surrounding different regions of the local anesthetics. We demonstrated the nonhomogeneous character of the anesthetics with respect to their affinity to water solvent molecules as well as the modifications in their affinity to water caused by changes in their charge state. We also observed that the biological potency of the anesthetics is more related to the behavior of specific groups within the molecule, which are responsible for the interaction with the lipid phase of membranes, rather than the general properties of the molecule as a whole. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007  相似文献   

12.
    
The action of local anesthetics (LA) is controversial. There is experimental evidence that the unprotonated form of LA penetrates the axon, while the charged form acts in the intracellular phase. To obtain some insight on the structure of the local anesthetics tetracaine and its pharmacological action, we made calculations using the density functional theory (DFT) method. After those calculations, we performed molecular dynamics (MD) simulations in a p, N, T ensemble, in an aqueous environment, on both unprotonated and protonated forms of the molecule. The radial distribution function was used to study water solvent effects, through the characterization of the affinity of tetracaine to water. The results indicate that the molecule has regions with different degree of hydrophobicity, and the N‐terminal of the anesthetic was primarily affected by changes in the protonation state of the anesthetic. The pH‐dependent activity of TTC should then be analyzed in view of local changes in different regions of the molecule, rather than in terms of general effects on the hydrophobicity of the molecule as a whole. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006  相似文献   

13.
黄酮类化合物的密度泛函理论研究   总被引:1,自引:0,他引:1  
在混合密度泛函B3LYP理论下,用6-31G*基函数研究了几种典型黄酮类化合物分子的几何结构、电子结构和分子的静电势,讨论了电子结构和分子活性部位的关系.  相似文献   

14.
    
Theoretical methods involving molecular dynamics (MD) simulation and density functional theory were performed to investigate the different molecular ratios, mechanical Properties, structure, trigger bond, and intermolecular interaction of hexaazaisowurtzitane (CL‐20)/nitroguanidine (NQ) cocrystal explosive. Results of MD simulation show that CL‐20 and NQ packed in ratios of 1:1 present the larger binding energy and better mechanical properties than any other molecular ratios, which indicates 1:1 cocrystal can form the stable crystal structure. Shorter length and larger dissociation energy of trigger bond in composite structure than in isolated CL‐20 component suggests that the cocrystal may exhibit less sensitive than CL‐20. Analyses of atoms in molecules, reduced density gradient, and natural bond orbital confirm that intermolecular interactions are mainly derived from a series of weak hydrogen bond and strong vdW forces, involving of NH···O, CH···O, CH···N, O···N, and O···O. Additionally, composite structures of 2 and 3 bringing us more attractive performance will act as a key role in constructing of CL‐20/NQ cocrystal explosive. © 2015 Wiley Periodicals, Inc.  相似文献   

15.
16.
用量子化学的密度泛函理论计算了12种有机二硫化物和铁原子簇的分子轨道指数及其与铁原子簇的化学吸附作用能, 探讨了这种作用能与抗磨性能的关系; 运用轨道能量近似原则讨论了有机二硫化物与铁原子的作用方式; 以前线电子密度、超离域性指数和原子净电荷作为判据分析了12种有机二硫化物与铁原子间键合的强弱、反应性的大小等表征有机二硫化物与金属作用强弱的参数。结果表明: 有机二硫化物与铁接触时, 在较缓和条件下, SS键优先断裂与金属发生化学吸附形成配位键, 起到抗磨作用; 在高负荷下, 与金属发生常规条件下不能发生的化学反应, 即CS键断裂生成无机膜, 起到极压作用; 且随着碳链的增长, 有机二硫化物的抗磨性能愈来愈好, 但极压性能愈来愈差; 运用量子化学计算得到的预测结果与摩擦学试验结果具有良好的一致性, 可为同类极压添加剂化合物的分子设计提供较为可靠的参考依据和理论方法。  相似文献   

17.
胡宗超  卫海燕  王凡  赵琦华  陈志达 《化学学报》2004,62(20):1973-1980,F005
用密度泛函理论结合对称性破损方法(DFF-BS)研究了混合桥联三核镍配合物的磁交换耦合作用.这类化合物是由三唑和异硫氰酸根桥联形成的混合桥配合物.计算表明,在标题化合物中,三唑桥传递反铁磁耦合作用,而异硫氰酸根桥传递铁磁耦合作用;并且,随着异硫氰酸根取代三唑桥的数目增加,配合物的铁磁作用增强,在一定意义上说明了混合桥磁耦合作用的加合性.Mulliken自旋布居分析表明,无论是三唑桥还是异硫氰酸根桥,它们的磁交换作用机理都是磁中心的自旋离域.分子磁轨道分析显示,对于三唑桥,在局域磁轨道之间存在着强的轨道作用,导致了反铁磁耦合;对于异硫氰酸根桥,局域磁轨道之间弱的相互作用,表现了铁磁耦合作用.对标题化合物的研究说明了DFF-BS方法可用于三核体系磁交换作用的研究.  相似文献   

18.
<正>Intermolecular geometries and energies of the dimers formed between simple aromatic hydrocarbons(benzene,toluene and naphthalene)and sulfate species(sulfuric acid, bisulfate and sulfate)were studied by density functional theory.The AH:H_2SO_4 complexes were determined by the OH-πH-bond,with H_2SO_4 acting as the H-bond donor and AHs as the acceptor. However,the AHs:HSO_4~-and AHs:SO_4~(2-)complexes were established by the H-bond,with AHs serving as acid and HSO_4~-or SO_4~(2-)as alkali.The atmospheric implications of those complexes were strongly supported by their considerable binding energies.  相似文献   

19.
20.
    
A recently developed empirical dispersion correction (Grimme et al., J. Chem. Phys. 2010, 132, 154104) to standard density functional theory (DFT‐D3) is implemented in the plane‐wave program package VASP. The DFT‐D3 implementation is compared with an implementation of the earlier DFT‐D2 version (Grimme, J. Comput. Chem. 2004, 25, 1463; Grimme, J. Comput. Chem. 2006, 27, 1787). Summation of empirical pair potential terms is performed over all atom pairs in the reference cell and over atoms in shells of neighboring cells until convergence of the dispersion energy is obtained. For DFT‐D3, the definition of coordination numbers has to be modified with respect to the molecular version to ensure convergence. The effect of three‐center terms as implemented in the original molecular DFT‐D3 version is investigated. The empirical parameters are taken from the original DFT‐D3 version where they had been optimized for a reference set of small molecules. As the coordination numbers of atoms in bulk and surfaces are much larger than in the reference compounds, this effect has to be discussed. The results of test calculations for bulk properties of metals, metal oxides, benzene, and graphite indicate that the original parameters are also suitable for solid‐state systems. In particular, the interlayer distance in bulk graphite and lattice constants of molecular crystals is considerably improved over standard functionals. With the molecular standard parameters (Grimme et al., J. Chem. Phys. 2010, 132, 154104; Grimme, J. Comput. Chem. 2006, 27, 1787) a slight overbinding is observed for ionic oxides where dispersion should not contribute to the bond. For simple adsorbate systems, such as Xe atoms and benzene on Ag(111), the DFT‐D implementations reproduce experimental results with a similar accuracy as more sophisticated approaches based on perturbation theory (Rohlfing and Bredow, Phys. Rev. Lett. 2008, 101, 266106). © 2012 Wiley Periodicals, Inc.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号