一种适用于原子团簇的多体展开新方案:以氮团簇为例（英文）

多体展开方法虽然已经广泛地用于估算弱相互作用体系的能量,但是其并不适用于计算共价团簇和金属团簇的能量.本文提出了一种适用于计算共价体系能量的相互作用多体展开(IMBE)方法.在相互作用多体展开方法中,体系的能量表示为孤立原子的能量及该原子与其他周围原子间相互作用的和.首先将该方法应用于计算氮团簇的能量,且多体展开截断至四体项.结果表明:与传统的多体展开方法相比,相互作用多体展开方法可以显著地降低能量误差.另外,以密度泛函理论计算结果为参考,相互作用多体展开方法估算能量的误差不依赖于体系的大小和结构,说明相互作用多体展开方法比较适合用于估算共价相互作用大体系的能量.     

基于内核构建的Cu-Au-Pd团簇稳定结构优化

具有特殊催化、磁性和化学活性的三元合金团簇已成为基础科学研究的热点问题.确定其稳定结构是研究团簇性质的重要前提.针对大尺寸Cu-Au-Pd团簇结构优化,提出了内核构建的方法改进了自适应免疫优化算法的效率(称为AIOA-IC算法).采用基于紧束缚势二阶矩近似的多体Gupta势函数来描述三元合金团簇原子间相互作用.为测试算法效率优化了原子数为60的Ag-Pd-Pt团簇稳定结构.结果显示新得到的结构比文献报道的团簇结构势能量值更低,由此可知AIOA-IC算法具有更强的势能面搜索能力.运用该算法研究了38及55原子Cu-Au-Pd团簇的稳定结构.所研究的38原子Cu-Au-Pd团簇包含了五折叠、六折叠和截角八面体结构,并且原子成分比例影响了团簇的结构类型.而55原子Cu-Au-Pd团簇均为完整二十面体结构,序列参数显示Cu,Au和Pd原子分层现象明显.对于147原子Cu_(12)Au_(93)Pd_(42)团簇完整二十面体结构,中心原子为Au,内层和次外层分别被12个Cu原子和42个Pd原子占据,最外层则被92个Au原子占满.通过原子半径及表面能分析了Cu,Pd和Au原子分别倾向于分布在内层、次外层和最外层的规律.     

铍团簇Be_n(n=2-55)的基态能量、结构及其拓扑性质

采用一种第一性原理多体势来描述团簇Ben中原子之间的相互作用,用模拟自然界生物进化的方法系统地优化计算了Ben团簇(n=2-55)的基态能量和几何结构,分析了其拓扑性质,总结了能量和结构随原子数目的变化规律,并将得到的Ben的结果和长程Morse势所描述的金属团簇Mn,以及LennardJones势所描述的惰性气体团簇LJn进行了比较     

纳米铜团簇在常温和升温过程中能量特征的分子动力学研究

采用分子动力学方法和F-S多体势函数，模拟研究纳米铜团簇常温下能量特征及其在升温直到熔化过程中的变化，确定了常温下纳米铜团簇的表面原子厚度和表面能，给出在不同温度下纳米铜团簇能量大小分布比例和能量的概率密度，细致描述了团簇升温过程团簇内部原子和表面原子之间不同的变化特征. 关键词： 铜团簇 分子动力学 能量特征 温度     

平面团簇稳定结构的蒙特卡罗树搜索

以平面团簇为例提出了一种结合结构识别和蒙特卡罗树技术搜索稳定结构的新方法.体系原子之间的相互作用由两类模型势能函数来描述:Lennard-Jones二体势函数与基于Lennard-Jones势的三体势函数.考虑可能的三角晶格碎片作为候选结构,引入编号策略对结构进行快速识别,并运用蒙特卡罗树搜索研究稳定结构随着原子数增大的演化过程;对于能量较低的候选结构,进一步采取局域优化来获得对应体系的稳定结构.计算表明,Lennard-Jones二体势函数对应的三角晶格团簇更稳定;在特定的参数下,三体势函数对应的六角晶格团簇更稳定.结合结构识别和蒙特卡罗树搜索可以对候选结构空间进行高效扫描,在较短时间内更容易搜索到稳定的团簇结构,并可以与第一原理计算结合实现材料的结构预测.     

Agn-1Si（n=5-10）团簇结构与性质研究

采用密度泛函理论研究了Agn-1Si团簇（n = 5-10）的几何结构和物理性质。首先得到体系最低能量结构,其中Ag5Si 和 Ag7Si团簇的结构比之前研究中结构能量更低。通过分析相应结构的能隙,平均结合能,二阶能量差可以发现Si原子的加入可以加强团簇结构稳定性,使团簇更加紧凑。在团簇尺寸n = 5-10的范围里,拥有八个价电子的Ag4Si团簇在以上三个方面都显示出非常稳定的特点。通过分析Agn-1Si团簇 (n = 5-10)的差分电荷发现,电荷的转移主要发生在Si原子与其相邻的Ag原子之间, Si原子和附近的Ag原子之间产生了强烈的共价相互作用,是团簇稳定性增强的重要原因。     

第一性原理对SinMg(n=1-12)团簇结构、稳定性与电子性质的研究

运用密度泛函理论(DFT)的杂化密度泛函B3LYP方法,在6-311G基组水平上对 团簇进行了构型优化、频率分析与电子性质计算.同时讨论了团簇的平均结合能、能级间隙、二阶能量差分、自然电子布居、极化率.研究结果表明： 团簇的基态绝大多数为立体结构. 时,体系的基态为自旋三重度, 时,则为单重态.镁原子的掺入使得主团簇的电子性质发生了明显的变化,掺杂使体系的平均结合能降低,能隙减小,化学硬度减小,电子亲和能增大.电子总是从 原子向 原子转移.团簇中原子之间的成键相互作用随n的增大 而增强,团簇的电子结构随n的增大而趋于紧凑.     

第一性原理计算Li_2Be_N(N=1~10)团簇的最低能量结构及其电子性质

本文从第一性原理出发,利用密度泛函理论(DFT)计算了Li2BeN(N=1~10)团簇的最低能量结构及其电子性质.计算结果表明,Li2BeN(N=1~10)团簇最低能量结构的对称性与单一成分的铍团簇相比有所下降,能量二阶有限差分随团簇尺寸的增加出现了奇偶振荡现象.从结构稳定性上来看,N=9是Li2BeN(N=1~10)团簇的一个幻数.我们同时发现了团簇中从Be到Li的电荷转移以及小的Li2BeN(N≤6)团簇中类范德瓦尔斯相互作用和共价相互作用共存的现象.     

密度泛函理论研究Be_nNa(n=2~14)团簇的结构及其电子性质(英文)

基于密度泛函理论研究了BenNa(n=2~14)团簇的最低能量结构及其电子性质,具体采用的计算方法是B3LYP/6-311G(d).对于每一尺寸的团簇,我们从多个可能构型出发来寻找其最低能量结构.计算结果表明:钠原子易于附着在铍团簇的表面;Be4Na和Be9Na均具有较高的平均结合能及能隙,n=4,9是团簇的幻数;在所有的BenNa(n=2~14)团簇中,钠原子均失去电子;随着团簇尺寸n的增大,团簇内部原子间的相互作用是共价键和金属键共存.     

Co-Ag与Co-Au合金团簇核层结构模拟

钴基二元合金团簇具有特殊催化和电磁性能已成为研究的热点.探究Co-Ag与Co-Au团簇稳定结构是研究其性质的重要步骤.运用基于内核构建的自适应免疫优化算法对大尺寸Co_(55)Ag_n、Co_(55)Au_n(n=1-55, 60, 70, 80, 92)团簇结构进行优化.基于紧束缚势二阶矩近似的多体Gupta势函数被用于描述合金团簇原子间相互作用.计算结果显示Co_(55)Ag_n、Co_(55)Au_n(n=1-55, 60)最稳定结构的主要构型为Mackay二十面体.原子分布分析显示Co原子位于内层,而Ag、Au原子均趋于分布在外层,呈现出核层状态.对于大尺寸Co_(55)Ag_n、Co_(55)Au_n(n=80, 92)团簇,Co-Ag团簇采取无定形结构,而Co-Au团簇为规则的二十面体结构.     

Competition between face-centered cubic and icosahedral cluster structures

The binding energy of atoms in icosahedral and face-centered-cubic clusters is calculated numerically for pairwise Morse-potential interactions between atoms and for clusters containing from 561 to 923 atoms, which corresponds to gradual filling of the sixth layer of the icosahedral cluster. Perturbation theory is used to calculate the cluster binding energy, in which the small parameter is the ratio of the interaction energy between non-nearest neighbor atoms to the interaction energy between nearest neighbors. Values of the Morse interaction potential parameter are found for which the energies of clusters with different structures coincide. Under the conditions used in these computations, the strain energy of a cluster can be neglected. Although the contribution of the interaction energy between non-nearest neighbors to the total cluster energy is small, it turns out to be important in finding the level crossing parameter. Zh. éksp. Teor. Fiz. 112, 1082–1090 (September 1997)     

The structure and energetics of carbon-nitrogen clusters

The structures and energies of clusters of carbon and nitrogen with up to 12 atoms have been investigated by density functional theory using the hybrid B3LYP functional and the cc-pVTZ basis set. This is the first systematic study of these clusters. Geometries are reported for the lowest energy states at this level of theory. Linear structures tend to be the global minima for clusters containing one or two nitrogen atoms, and patterns in the electronic structure of these clusters are reported. More complex branched structures lie close in energy to the linear conformations and, for clusters greater than six atoms and containing three or more nitrogen atoms, these branched structures are the minimum energy conformers. Comparisons are made with pure carbon and silicon-carbon clusters.     

Microscopic theory of the second- and third-order elastic constants of cobalt in a thermodynamic module

The total potential energy of a crystal is presented as an expansion by irreducible interactions in clusters containing pairs of atomic triplets and quadruplets. The potential energy of the clusters in the adiabatic approximation is a function of vectors \(\vec r_{ik}\) connecting the centers of atoms in the clusters. Arguments (basic functions) affecting the potential energy of clusters are found using the model with allowance for the exchange symmetry of atoms and the irreducibility of the considered energies of doublet, triplet, and quadruplet atoms and interactions. This allows us to present arguments of the potential energy in the form of summed integral numbers (latticed sums) multiplied by a fixed value of the crystal unit cell parameter, and to set the numerical values of the potential energy arguments as a function of vectors \(\vec r_{ik}\). A potential corresponding to the thermodynamic additivity concept of crystal energy is selected as the model potential of pairwise interaction. Secondand third-order moduli of the rigidity of Co crystals with A1, A2, and HCP structures are calculated using this model of multiatom interactions.     

Graphene on Ir(111): physisorption with chemical modulation

The nonlocal van der Waals density functional approach is applied to calculate the binding of graphene to Ir(111). The precise agreement of the calculated mean height h = 3.41 ? of the C atoms with their mean height h = (3.38±0.04) ? as measured by the x-ray standing wave technique provides a benchmark for the applicability of the nonlocal functional. We find bonding of graphene to Ir(111) to be due to the van der Waals interaction with an antibonding average contribution from chemical interaction. Despite its globally repulsive character, in certain areas of the large graphene moiré unit cell charge accumulation between Ir substrate and graphene C atoms is observed, signaling a weak covalent bond formation.     

原子团簇P10模型的理论计算

用分子图形学方法设计出２６种Ｐ１２模型,并对其进行了分子力学、ＰＭ３半经验量子化学和Ｈａｒｔｒｅｅ－Ｆｏｃｋ从头算优化。在Ｐ１０原子团簇模型设计中,磷原子采用一、二、三或四配位。大部分Ｐ１０的模型是在Ｐ９＾＋和Ｐ８的模型上分别增加、,２个原子生成的。这些模型包括１５种在势能面上局域极小点和１１种鞍点（或过流态）。从模型优化后的能量比较可知,２个四面体Ｐ４与１个Ｐ２通过４个单键连接的桥式结构最稳定。     

Theoretical study of interaction between atoms of Au, Ag, Cu and clean Si(1 1 1) surface

To evaluate the interactions between the atoms of Au, Ag and Cu and clean Si(1 1 1) surface, two types of silicon clusters Si₄H₇ and Si₁₆H₂₀ together with their metal complexes were studied by using hybrid (U)B3LYP density functional theory method. Optimized geometries and energies on different adsorption sites indicate that: (1) the binding energies at different adsorption sites are large (ranging from 1.2 to 2.6 eV depend on the metal atoms and adsorption sites), suggesting a strong interaction between metal atom and silicon surface; (2) the most favorable adsorption site is the on top (T) site. Mulliken population analysis indicated that in the system of on top (T) site, a covalent bond is formed between metal atom and dangling bond of surface Si atom.     

GaAs Electron Spectrum with Antisite AsGa Defect Clusters in the Semiempirical Tight-Binding Method

A computational procedure for calculating the band structure of a crystal with antisite-defect clusters is developed using a large unit-cell method. GaAs electron spectra containing clusters of more than 100 As atoms characterized by crystal point-group symmetry are calculated. The spread of energy levels in the band gap is examined. The origin of these energy levels is interpreted as being the result of splitting of energy levels of single defects.     

Efimov Physics in Small Bosonic Clusters

We study small clusters of bosons, A = 2, 3, 4, 5, 6, characterized by a resonant interaction. Firstly, we use a soft-gaussian interaction that reproduces the values of the dimer binding energy and the atom-atom scattering length obtained with LM2M2 potential, a widely used ⁴He-⁴He interaction. We change the intensity of the potential to explore the clusters’ spectra in different regions with large positive and large negative values of the two-body scattering length and we report the clusters’ energies on Efimov plot, which makes the scale invariance explicit. Secondly, we repeat our calculation adding a repulsive three-body force to reproduce the trimer binding energy. In all the region explored, we have found that these systems present two states, one deep and one shallow close to the A ? 1 threshold, and scale invariance has been investigated for these states. The calculations are performed by means of Hyperspherical Harmonics basis set.     

Interaction Energy in Polymer Blends Containing N-1-Alkylitaconamic Acids Moiety

Abstract

Analysis of the interaction energy of N-1-alkylitaconamic acids to explain the miscibility behavior of polymers containing this moiety is performed. In order to calculate the interaction energy, many conformations of the N-1-alkylitaconamic acid fragments containing ethyl, propyl, butyl, hexyl, octyl, decyl, and dodecyl were selected randomly from molecular dynamic simulations. It was assumed that miscibility and other properties are determined by the enthalpy of mixing, ΔH_mix, that ΔH_mix is dominated by the local interactions between segments of the polymer chains, and, furthermore, that PΔV_mix contribution can be ignored. The polymer miscibility is discussed by using interaction energies calculated from molecular mechanics studies on pairs of small polymer fragments. By this method and by using the Blend package (Accelrys), good correlations between the interaction energy and the Gordon Taylor constant (K_GT) values for blends of poly(N-1-alkylitaconamic acids) and poly(4-vinylpyridine) previously reported are obtained. The analysis of the energies shows that Coulombic energies decay rapidly being almost “zero” as the side-chain lengths reach near eight carbon atoms. On the contrary, the van der Waals interaction (VDW) energies decay in a linear way as the length of the side chain increases. The slope of the VDW energies is an estimation of the energetic contribution per methylene unit to the blending process. The slope in this case is 570 cal/mol, which is a value very close to that experimentally reported in the literature for related systems. The blend energy was calculated, and good correlation with the experimental results is obtained.     

Band structure and quantum conductance of nanostructures from maximally localized Wannier functions: the case of functionalized carbon nanotubes

We have combined large-scale, tau-point electronic-structure calculations with the maximally localized Wannier functions approach to calculate efficiently the band structure and the quantum conductance of complex systems containing thousands of atoms while maintaining full first-principles accuracy. We have applied this approach to study covalent functionalizations in metallic single-walled carbon nanotubes. We find that the band structure around the Fermi energy is much less dependent on the chemical nature of the ligands than on the sp(3) functionalization pattern disrupting the conjugation network. Common aryl functionalizations are more stable when paired with saturating hydrogens; even when paired, they still act as strong scattering centers that degrade the ballistic conductance of the nanotubes already at low degrees of coverage.