银团簇结构与特性的分子动力学研究

用分子动力学方法模拟了银团簇的结构与力能学.计算模拟中使用了一种基于第一性原理的原子间互作用多体势函数.通过分子动力学模拟确定了银微团簇(原子个数3～13)的稳态结构;模拟了原子个数为13～141的银FCC晶体结构理想球形团簇的力能学,发现球形银团簇形成三雏紧密结构;计算了平均结合能,给出了结合能随团簇原子数N的变化图,发现随N增大团簇结合能逐渐接近块材的数值.     

PSO-SA模拟铝团簇的几何结构与基态能量

介绍了采用基于模拟退火的粒子群算法(PSO-SA),应用Gupta多体势函数和Sutton-Chen势函数对铝团簇Aln(n≤13)的基态能量进行了计算,并给出其相应的稳定结构.同时研究了铝团簇的半径分布图、平均结合能和能量二阶差分与团簇尺寸变化的关系.结果表明:铝团簇的结构稳定性随团簇原子数n的增加而增大,并在原子数n为7,9和11处出现极值,表明其具有“幻数”效应和相对较高的结构稳定性,由半径的分布图得出Al6,Al7,Al13具有高度的对称性.     

荷电(±1)对二十面体X_(13)(X=Cr,Mn,Fe,Co)团簇的稳定性与磁性的影响

采用基于密度泛函理论的计算方法,对正二十面体金属X13(X=Cr,Mn,Fe,Co)中性和荷电团簇进行了全面的结构优化计算,研究了荷电对团簇的稳定性和磁性的影响.结果表明:荷负电能够使团簇的稳定性增强;荷电对不同团簇的原子间距离的影响不同;同时荷电对不同团簇磁性的影响也是不一样的,尤其是荷负电能够使Fe13和Co13团簇的磁性大大增强;荷电对不同团簇磁性的影响不是通过原子间距离的变化来实现的,而是受到原子内部电荷的转移和杂化程度的影响.     

金镍二元团簇结构和电子性质的理论研究

在UBP86/LANL2DZ和UBP86/def2-TZVP水平下详细研究了AumNin (m＋n≤6)团簇的几何结构和电子性质. 详细地分析了团簇的结构特征, 平均结合能, 垂直电离势, 垂直电子亲和能, 电荷转移以及成键特征. 所有混合团簇中, 镍原子趋于聚集到一起, 形成最多Ni—Ni键, 金原子分布在镍原子聚集体周围以形成最多Au—Ni键. Ni原子较少团簇的电子性质与纯金团簇类似, 呈现一定奇偶振荡. 混合团簇中存在镍到金原子间的电荷转移. Ni原子较少团簇中, 自旋电子主要定域在Ni原子上, Ni原子较多团簇中, Au原子明显受到自旋极化. 混合团簇的分波态密度表明, AuNi混合团簇对小分子的反应活性要高于纯金团簇.     

Ge_nEu(n=1-13)团簇的结构、电子及磁性质(英文)

利用密度泛函理论在广义梯度近似下研究了GenEu(n=1-13)团簇的生长模式和磁性.结果表明:对于GenEu(n=1-13)团簇的基态结构而言,没有Eu原子陷入笼中.这和SinEu以及其它过渡金属掺杂半导体团簇的生长模式不同.除GeEu团簇外,GenEu(n=2-13)团簇的磁矩均为7μB.团簇的总磁矩与Eu原子的4f轨道磁矩基本相等.Ge、Eu原子间的电荷转移以及Eu原子的5d、6p和6s间的轨道杂化可以增强Eu原子的局域磁矩,却不能增强团簇总磁矩.     

金钯二元小团簇的几何结构与电子性质

在UBP86/LANL2DZ和UB3LYP/def2-TZVP水平下详细研究了AumPdn(m+n≤6)团簇的几何结构和电子性质.阐明了团簇的结构特征、平均结合能、垂直电离势、垂直电子亲和能、电荷转移以及成键特征.除单取代混合团簇(AunPd和AuPdn,n=5或6)外,五和六原子混合团簇中钯原子趋于聚集到一起形成Pdcore,金原子分布在Pdcore周围形成PdcoreAushell结构.含一个和两个钯原子团簇的电子性质与纯金团簇类似,呈现一定奇偶振荡.混合团簇的电子性质,如最高占据分子轨道(HOMO),最低未占据分子轨道(LUMO),垂直电离势,垂直电子亲和能,Fermi能级和化学硬度等均与团簇空间结构和金、钯原子数之比直接相关.混合团簇中存在钯原子到金原子间的电荷转移,表明团簇中存在明显金钯间成键作用.分析团簇的电荷分布、前线轨道和化学硬度表明,金钯混合团簇对小分子如O2、H2和CO等的反应活性要强于纯金团簇.     

嵌入纳米Fe颗粒的Fe液凝固过程的分子动力学模拟

采用Sutton-Chen 势函数及分子动力学(MD)方法对嵌入了Fe纳米团簇(半径从0.4-1.8 nm)的Fe液凝固过程进行了模拟. 模拟结果表明只有当嵌入的纳米团簇半径超过0.82 nm才能降低凝固时所需要的临界过冷度(ΔT^*), 起到诱导凝固的作用. 同时采用原子键型指数法(CTIM-2)对样本在凝固过程中的原子结构进行了标定, 通过观察微观结构演变发现当嵌入纳米团簇能够作为凝固核心时, 体系按照hcp-fcc 交叉形核的方式长大. 同时还发现嵌入纳米团簇对体系凝固过程晶核的生长方向及凝固的最终构型存在“结构遗传效应”.     

硅、锗、锡、铅/磷二元原子团簇的形成、光解和结构

利用激光溅射产生了第IV主族 (硅、锗、锡、铅 ) /磷二元团簇正负离子 ,用飞行时间质谱研究了团簇离子的组成规律和激光光解产物 .研究表明二元团簇稳定性受团簇电子结构和几何结构的影响 ,但随着第IV主族元素自上而下 ,几何结构对团簇稳定性的作用越来越大 .在二元团簇离子中存在两类幻数团簇 :一类可以用Wade规则解释 ,其中磷原子或者充当给电子配体结合在第IV主族原子构成的团簇骨架外 ,或者直接参与团簇骨架的构成 ;另一类则与稳定的第IV主族中性团簇 (或磷中性团簇 )是等电子体 .利用从头计算和Wade规则对幻数团簇的结构和价键进行了分析 .     

(BN)n(n≤12)团簇的结构及成键性质

利用遗传算法和Gastreich提出的经验势函数研究了(BN)_n(n≤12)团簇的可能稳定结构, 并对能量较低的异构体在HF/6-31G(d)水平进行优化, 得到了(BN)_n(n≤12)团簇的线状、蒲扇形、单环、双环、三环和笼状结构, 讨论了各种结构的特征及相对稳定性. 分析了BN团簇中原子的成键性质, 在单环结构中, N原子以sp2杂化成键, B原子以sp杂化成键, 而在节点处B原子以sp2杂化成键. (BN)6是唯一没有张力的单环结构.     

微波频率下氯化钠溶液电导率的非线性特性

在微波频段研究了电磁场对氯化钠溶液电导率的影响. 运用不对称周期双阱势函数和Langvin方程描述了恒温条件下氯化钠溶液团簇结构中平均氢键的动力学行为, 由Bertolini公式建立团簇的平均氢键形成的概率变化与电导率变化的关系. 计算结果表明, 在恒温条件下, 当微波场的场强达到10 kV·m-1 时, 微波能会部分转化为团簇里分子间势能, 使得氯化钠溶液的电导率与电场强度有关, 这一结果与实验结果符合. 这种微波与氯化钠溶液的相互作用表现出强烈的非线性变化, 应该属于微波的非热效应范畴.     

Linear augmented Slater-type orbital method for free standing clusters

We have developed a Scalable Linear Augmented Slater-Type Orbital (LASTO) method for electronic-structure calculations on free-standing atomic clusters. As with other linear methods we solve the Schr?dinger equation using a mixed basis set consisting of numerical functions inside atom-centered spheres and matched onto tail functions outside. The tail functions are Slater-type orbitals, which are localized, exponentially decaying functions. To solve the Poisson equation between spheres, we use a finite difference method replacing the rapidly varying charge density inside the spheres with a smoothed density with the same multipole moments. We use multigrid techniques on the mesh, which yields the Coulomb potential on the spheres and in turn defines the potential inside via a Dirichlet problem. To solve the linear eigen-problem, we use ScaLAPACK, a well-developed package to solve large eigensystems with dense matrices. We have tested the method on small clusters of palladium.     

Equivariant Morse theory of theN-body problem: Application to potential surfaces in chemistry

Summary The Morse inequalities linking the critical points of a potential function on the whole configuration space and its restrictions to either planar or linear configurations are derived from the Morse theory in its equivariant form. Brute potential functions arising from standard models of quantum chemistry need eventually morsification which can be achieved without altering the main chemical significances of the potential. Illustrative applications follow in the case of magnesium clusters.     

Size-specific interaction of alkali metal ions in the solvation of M+-benzene clusters by Ar atoms

The size-specific influence of the M+ alkali ion (M = Li, Na, K, Rb, and Cs) in the solvation process of the M+-benzene clusters by Ar atoms is investigated by means of molecular dynamic simulations. To fully understand the behavior observed in M+-bz-Ar(n) clusters, solvation is also studied in clusters containing either M+ or benzene only. The potential energy surfaces employed are based on a semiempirical bond-atom decomposition, which has been developed previously by some of the authors. The outcome of the dynamics is analyzed by employing radial distribution functions, studying the evolution of the distances between the Ar atoms and the alkali ion M+ or the benzene molecule for all M+-bz-Ar(n) clusters. For all members, in the M+-bz series, the benzene molecule (bz) is found to remain strongly bound to M+ even in the presence of solvent atoms. The radial distribution functions for the heavier clusters (K+-bz, Rb+-bz, and Cs+-bz), are found to be different than for the lighter (Na+-bz and Li+-bz) ones.     

Theoretical Study of Structural Symmetry in Ternary Clusters

The geometrical symmetry presents an intriguing theoretical problem in many kinds of clusters. The diversity of geometrical structures is associated with cluster sizes, different model functions and potential parameters, and ternary clusters are investigated to study the relationship between geometrical symmetry and homotopic symmetry. Ternary Lennard-Jones model potential is studied with different parameters, and the putative global minimum structures of A₁₃B₃₀C₁₂ clusters are optimized using an adaptive immune optimization algorithm. The results show that there mainly exist five geometrical symmetry structures, i.e., Mackay icosahedral, fivefold partial Mackay icosahedral, sixfold pancake, partial double Mackay icosahedral, and amorphous structures. Furthermore, the number of bonds is used to distinguish the geometrical symmetry. The importance of geometrical symmetry and homotopic symmetry determined by potential parameters is discussed. It was found that in the optimization it is more important to generate geometrical symmetry than to optimize homotopic symmetry.     

The mathematical modelling of cluster geometry

Calculations of minimum energy configurations for aggregates of up to forty atoms, commonly referred to as clusters, are presented. In contrast to previous studies, random initial configurations have been optimised to find the lowest energy structure for a given number of atoms. Three different two-body, bireciprocal potential functions were used in these calculations and in the case of the Lennard-Jones potential, previously calculated results have been confirmed. New structures obtained using softer potentials are also presented. Minimum energy structures of small clusters containing two different types of atoms have also been calculated, and the relationship between the geometry of a cluster and the relative sizes of its constituent atoms examined.     

Importance of Madelung potential in quantum chemical modeling of ionic surfaces

The importance of the inclusion of the Madelung potential in cluster models of ionic surfaces is the subject of this work. We have determined Hartree-Fock all electron wave functions for a series of MgO clusters with and without a large array of surrounding point charges (PC) chosen to reproduce the Madelung potential. The phenomena investigated include: the reactivity of surface oxygens toward CO₂, atomic hydrogen, and H⁺; the geometry and adsorption energy of water and the vibrational shift of CO adsorbed at Mg²⁺ sites; the electronic structure and the hyperfine coupling constants of oxygen vacancies, the paramagnetic F_s⁺ centers. While some clusters give results which are virtually independent of the presence of the PCs, other clusters, typically of small size, give physically incorrect results when the PCs are not included. The embedding of the cluster in PCs guarantees a similar response for clusters of different size, at variance with the bare clusters, where the long range coulombic interactions are not included. © 1997 by John Wiley & Sons, Inc.     

Competitive role of CH4-CH4 and CH-π interactions in C6H6-(CH4)n aggregates: the transition from dimer to cluster features

The intermolecular methane-methane and benzene (Bz)-methane interactions formulated in this paper are suitable to investigate systems of increasing complexity. The proposed CH(4)-CH(4) and Bz-CH(4) potential energy functions are indeed applied to study some macroscopic properties of methane and important features of both small Bz-(CH(4))(n) (n > 1-10) clusters and Bz surrounded by several CH(4) molecules. Relevant parameters of the interaction, derived from molecular polarizability components, have been proved to be useful to describe in a consistent way both size repulsion and dispersion attraction forces. The proposed potential model also allows one to isolate the role of the different intermolecular energy contributions. The spatial distribution of the CH(4) molecules in the clusters is investigated by means of molecular dynamics simulations under various conditions, even when methane phase transition from liquid to gas is likely to occur. In addition, several properties, such as radial distribution functions, density values, and mean diffusion coefficients, are analyzed in detail.     

Simulation of Benzenesulfonyl Chloride Hydrolysis. Influence of the Size and Structural Ordering of Aqueous Clusters on Thermodynamic and Activation Parameters of the Process

The potential energy surfaces for hydrolysis of benzenesulfonyl chloride in clusters incorporating one to three additional water molecules (along with the molecule acting as nucleophile) were calculated by the PM3 method in the supermolecular approximation. As in the gas phase, hydrolysis of benzenesulfonyl chloride in clusters is a two-step exothermic process occurring via relatively unstable five-coordinate inter- mediate. The structural and thermodynamic parameters of the intermediate and transition states of the reaction, and also the thermodynamic functions of the process were calculated for seven clusters of different structures.     

Quantum vibrational analysis and infrared spectra of microhydrated sodium ions using an ab initio potential

We present a full-dimensional potential energy surface and a dipole moment surface (DMS) for hydrated sodium ion. These surfaces are based on an n-body expansion for both the potential energy and the dipole moment, truncated at the two-body level for the H(2)O-Na(+) interaction and also for the DMS. The water-water interaction is truncated at the three-body level. The new full-dimensional two-body H(2)O-Na(+) potential is a fit to roughly 20,000 coupled-cluster single double (triple)/aug-cc-pVTZ energies. Properties of this two-body potential and the potential describing (H(2)O)(n)Na(+) clusters, with n up to 4 are given. We then report anharmonic, coupled vibrational calculations with the "local-monomer model" to obtain infrared spectra and also 0 K radial distribution functions for these clusters. Some comparisons are made with the recent infrared predissociation spectroscopy experiments of Miller and Lisy [J. Am. Chem. Soc. 130, 15381 (2008).].     

The properties of ion-water clusters. II. Solvation structures of Na+, Cl-, and H+ clusters as a function of temperature

Ion-water-cluster properties are investigated both through the multistate empirical valence bond potential and a polarizable model. Equilibrium properties of the ion-water clusters H+(H2O)100, Na+(H2O)100, Na+(H2O)20, and Cl-(H2O)17 in the temperature region 100-450 K are explored using a hybrid parallel basin-hopping and tempering algorithm. The effect of the solid-liquid phase transition in both caloric curves and structural distribution functions is investigated. It is found that sodium and chloride ions largely reside on the surface of water clusters below the cluster melting temperature but are solvated into the interior of the cluster above the melting temperature, while the solvated proton was found to have significant propensity to reside on or near the surface in both the liquid- and solid-state clusters.