A density functional study of YnAl (n=1-14) clusters

The geometries, stabilities, and electronic and magnetic properties of Y(n)Al (n=1-14) clusters have been systematically investigated by using density functional theory with generalized gradient approximation. The growth pattern for different sized Y(n)Al (n=1-14) clusters is Al-substituted Y(n+1) clusters and it keeps the similar frameworks of the most stable Y(n+1) clusters except for Y(9)Al cluster. The Al atom substituted the surface atom of the Y(n+1) clusters for n<9. Starting from n=9, the Al atom completely falls into the center of the Y-frame. The Al atom substituted the center atom of the Y(n+1) clusters to form the Al-encapsulated Y(n) geometries for n>9. The calculated results manifest that doping of the Al atom contributes to strengthen the stabilities of the yttrium framework. In addition, the relative stability of Y(12)Al is the strongest among all different sized Y(n)Al clusters, which might stem from its highly symmetric geometry. Mulliken population analysis shows that the charges always transfer from Y atoms to Al atom in all different sized clusters. Doping of the Al atom decreases the average magnetic moments of most Y(n) clusters. Especially, the magnetic moment is completely quenched after doping Al in the Y(13), which is ascribed to the disappearance of the ininerant 4d electron spin exchange effect. Finally, the frontier orbitals properties of Y(n)Al are also discussed.     

Microhydration of cytosine and its radical anion: cytosine.(H2O)n (n=1-5)

Microhydration effects on cytosine and its radical anion have been investigated theoretically, by explicitly considering various structures of cytosine complexes with up to five water molecules. Each successive water molecule (through n=5) is bound by 7-10 kcal mol(-1) to the relevant cytosine complex. The hydration energies are uniformly higher for the analogous anion systems. While the predicted vertical detachment energy (VDE) of the isolated cytosine is only 0.48 eV, it is predicted to increase to 1.27 eV for the lowest-lying pentahydrate of cytosine. The adiabatic electron affinity (AEA) of cytosine was also found to increase from 0.03 to 0.61 eV for the pentahydrate, implying that the cytosine anion, while questionable in the gas phase, is bound in aqueous solution. Both the VDE and AEA values for cytosine are smaller than those of uracil and thymine for a given hydration number. These results are in qualitative agreement with available experimental results from photodetachment-photoelectron spectroscopy studies of Schiedt et al. [Chem. Phys. 239, 511 (1998)].     

Density functional study of HO(H2O)n (n = 1–3) clusters

The hydrogen bonding complexes HO(H₂O)_n (n = 1–3) were completely investigated in the present study using DFT and MP2 methods at varied basis set levels from 6‐31++G(d,p) to 6‐311++G(2d,2p). For n = 1 two, for n = 2 two, and for n = 3 five reasonable geometries are considered. The optimized geometric parameters and interaction energies for various complexes at different levels are estimated. The infrared spectrum frequencies and IR intensities of the most stable structures are reported. Finally, thermochemistry studies are also carried out. The results indicate that the formation and the number of hydrogen bonding have played an important role in the structures and relative stabilities of different complexes. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

The mechanism of proton exchange: guided ion beam studies of the reactions, H(H2O)n+ (n=1-4)+D2O and D(D2O)n+ (n=1-4)+H2O

Reactions of protonated water clusters, H(H(2)O)(n) (+) (n=1-4) with D(2)O and their "mirror" reactions, D(D(2)O)(n) (+) (n=1-4) with H(2)O, are studied using guided-ion beam mass spectrometry. Absolute reaction cross sections are determined as a function of collision energy from thermal energy to over 10 eV. At low collision energies, we observe reactions in which H(2)O and D(2)O molecules are interchanged and reactions where H-D exchange has occurred. As the collision energy is increased, the H-D exchange products decrease and the water exchange products become dominant. At high collision energies, processes in which one or more water molecules are lost from the reactant ions become important, with simple collision-induced dissociation processes, i.e., those without H-D exchange, being dominant. Threshold energies of endothermic channels are measured and used to determine binding energies of the proton bound complexes, which are consistent with those determined by thermal equilibrium measurements and previous collision-induced dissociation studies. A kinetic scheme that relies only on the ratio of isomerization and dissociation rate constants successfully accounts for the kinetic energy dependence observed in the branching ratios for H-D and water exchange products in all systems. Rice-Ramsperger-Kassel-Marcus theory and ab initio calculations confirm the feasibility and establish the details of this kinetic model.     

Theoretical study of Al(n) and Al(n)O (n = 2-10) clusters

The stable structures, energies, and electronic properties of neutral, cationic, and anionic clusters of Al(n) (n = 2-10) are studied systematically at the B3LYP/6-311G(2d) level. We find that our optimized structures of Al5(+), Al9(+), Al9(-), Al10, Al10(+), and Al10(-) clusters are more stable than the corresponding ones proposed in previous literature reports. For the studied neutral aluminum clusters, our results show that the stability has an odd/even alternation phenomenon. We also find that the Al3, Al7, Al7(+), and Al7(-) structures are more stable than their neighbors according to their binding energies. For Al7(+) with a special stability, the nucleus-independent chemical shifts and resonance energies are calculated to evaluate its aromaticity. In addition, we present results on hardness, ionization potential, and electron detachment energy. On the basis of the stable structures of the neutral Al(n) (n = 2-10) clusters, the Al(n)O (n = 2-10) clusters are further investigated at the B3LYP/6-311G(2d), and the lowest-energy structures are searched. The structures show that oxygen tends to either be absorbed at the surface of the aluminum clusters or be inserted between Al atoms to form an Al(n-1)OAl motif, of which the Al(n-1) part retains the stable structure of pure aluminum clusters.     

Photoionization and ab initio study of Ba(H2O)n (n = 1-4) clusters

An experimental and theoretical study of the photoionization energies (IE's) of Ba(H(2)O)(n) clusters containing up to n = 4 water molecules has been performed. The clusters were generated by a pick-up source combining laser vaporization with pulsed supersonic expansion, and then photoionized by radiation of 272.5-340 nm. The experimentally determined IE(e)'s for n = 1 to 4 are 4.56 ± 0.05, 4.26 ± 0.05, 3.90 ± 0.05 and 3.71 ± 0.05 eV. This cluster size dependence of IE is reproduced within ±0.06 eV employing the mPW1PW91 density-functional and CCSD(T, Full) quantum-chemical methods combined with the 6-311++G(d,p) basis set for the H and O atoms and three different relativistic effective core potentials for Ba atoms. The calculations indicate that the lowest energy hydration structures represent the most relevant contributions to both the vertical and adiabatic ionization energies. Experimental and theoretical evidence correlates with the progressive surface-delocalization of the electron from the hydration cavity around the Ba atom and suggests that the intra-cluster electron transfer is possible even for small aggregates.     

Microhydration of X2 gas (X = Cl, Br, and I): a theoretical study on X2.nH2O clusters (n = 1-8)

Structure and properties of hydrated clusters of halogen gas, X2.nH2O (X = Cl, Br, and I; n = 1-8) are presented following first principle based electronic structure theory, namely, BHHLYP density functional and second-order Moller-Plesset perturbation (MP2) methods. Several geometrical arrangements are considered as initial guess structures to look for the minimum energy equilibrium structures by applying the 6-311++G(d,p) set of the basis function. Results on X2-water clusters (X = Br and I) suggest that X2 exists as a charge separated ion pair, X+delta-X-delta in the hydrated clusters, X2.nH2O (n > or = 2). Though the optimized structures of Cl2.nH2O clusters look like X2.nH2O (X = Br and I) clusters, Cl2 does not exist as a charge separated ion pair in the presence of solvent water molecules. The calculated interaction energy between X2 and solvent water cluster increases from Cl2.nH2O to I2.nH2O clusters, suggesting solubility of gas-phase I2 in water to be a maximum among these three systems. Static and dynamic polarizabilities of hydrated X2 clusters, X2.nH2O, are calculated and observed to vary linearly with the size (n) of these water clusters with correlation coefficient >0.999. This suggests that the polarizability of the larger size hydrated clusters can be reliably predicted. Static and dynamic polarizabilities of these hydrated clusters grow exponentially with the frequency of an external applied field for a particular size (n) of hydrated cluster.     

A density functional study on nitrogen-doped carbon clusters CnN3- (n=1-8)

Using molecular graphics software, we designed numerous models of CnN3- (n=1-8). Geometry optimization and calculation on vibration frequency were carried out by the B3LYP density functional method. After comparison of structure stability, we found that the structures of ground-state CN3- and C2N3- are bent chains with a nitrogen atom at either end, whereas when n=3-8, the ground-state clusters show three branches, each with a nitrogen atom located at the end. When n=5-8, the longest branch of CnN3- is polyacetylenelike. When n=5 or 7, the longest branch is connected to the central sp2 carbon in a nonlinear manner. The CnN3- (n=1-8) with an even number of carbon atoms are more stable than those with odd numbers, matching the peak pattern observed in laser-induced mass spectra of CnN3-. The trend of such odd/even alternation is explained based on concepts of bonding characteristics, electron affinities, and incremental binding energies.     

Structural evolution and stability of hydrogenated Li(n) (n = 1-30) clusters: a density functional study

The structure and electronic and optical properties of hydrogenated lithium clusters Li(n)H(m) (n = 1-30, m ≤ n) have been investigated by density functional theory (DFT). The structural optimizations are performed with the Becke 3 Lee-Yang-Parr (B3LYP) exchange-correlation functional with 6-311G++(d, p) basis set. The reliability of the method employed has been established by excellent agreement with computational and experimental data, wherever available. The turn over from two- to three-dimensional geometry in Li(n)H(m) clusters is found to occur at size n = 4 and m = 3. Interestingly, a rock-salt-like face-centered cubic structure is seen in Li(13)H(14). The sequential addition of hydrogen to small-sized Li clusters predicted regions of regular lattice in saturated hydrogenated clusters. This led us to focus on large-sized saturated clusters rather than to increase the number of hydrogen atoms monotonically. The lattice constants of Li(9)H(9), Li(18)H(18), Li(20)H(20), and Li(30)H(30) calculated at their optimized geometry are found to gradually approach the corresponding bulk values of 4.083. The sequential addition of hydrogen stabilizes the cluster, irrespective of the cluster size. A significant increase in stability is seen in the case of completely hydrogenated clusters, i.e., when the number of hydrogen atoms equals Li atoms. The enhanced stability has been interpreted in terms of various electronic and optical properties like adiabatic and vertical ionization potential, HOMO-LUMO gap, and polarizability.     

Doping-enhanced hyperpolarizabilities of silicon clusters: a global ab initio and density functional theory study of Si10 (Li, Na, K)n (n=1, 2) clusters

A global theoretical study of the (hyper)polarizabilities of alkali doped Si(10) is presented and discussed. First, a detailed picture about the low lying isomers of Si(10)Li, Si(10)Na, Si(10)K, Si(10)Li(2), Si(10)Na(2), and Si(10)K(2) has been obtained in a global manner. Then, the microscopic first (hyper)polarizabilities of the most stable configurations have been determined by means of ab initio methods of high predictive capability such as those based on the M?ller-Plesset perturbation and coupled cluster theory, paying extra attention to the (hyper)polarizabilities of the open shell mono-doped systems Si(10)Li, Si(10)Na, Si(10)K, and the influence of spin contamination. These results were used to assess the performance of methods of low computational cost based on density functional theory (DFT) in the reliable computation of these properties in order to proceed with an in-depth study of their evolution as a function of the alkali metal, the cluster composition, and the cluster structure. The most interesting outcomes of the performed (hyper)polarizability study indicate that while alkali doping leaves the per atom polarizability practically unaffected, influences dramatically the hyperpolarizabilities of Si(10). The lowest energy structures of the mono-doped clusters are characterized by significantly enhanced hyperpolarizabilities as compared to the analogue neutral or charged bare silicon clusters Si(10) and Si(11), while, certain patterns governed by the type and the number of the doping agents are followed. The observed hyperpolarizability increase is found to be in close connection with specific cluster to alkali metal charge transfer excited states and to the cluster structures. Moreover, an interesting correlation between the anisotropy of the electron density, and the hyperpolarizabilities of these systems has been observed. Finally, it is important to note that the presented method assessment points out that among the various DFT functionals used in this work, (B3LYP, B3PW91, BhandHLYP, PBE0, CAM-B3LYP, LC-BLYP, LC-BPW91) only B3PW91 and PBE0 out of the seven provided a consistent quantitative performance for both polarizabilities and hyperpolarizabilities with respect to the ab initio methods utilized here. On the other hand, the long range corrected functionals LC-(U)BLYP and LC-(U)BPW91 (μ = 0.47) failed to supply quantitatively accurate hyperpolarizability results in all the studied clusters while the CAM-(U)B3LYP functional performs satisfactory only in the case of the Na and K doped systems.     

Geometrical and electronic properties of neutral and charged cesium clusters Cs(n) (n=2-10): a theoretical study

We have investigated the structure and electronic properties of cesium clusters following all electron ab initio theoretical methods based on configuration interaction, second-order Moller-Plesset (MP2) perturbation theory, and density-functional theory. Becke's three-parameter nonlocal hybrid exchange-correlation functional (B3LYP) is found to perform best on the present systems with a split valence 3-21G basis function. We have calculated the optimized geometries of neutral and singly charged cesium clusters having up to ten atoms, their binding energy per atom, ionization potentials (IPs), and adiabatic electron affinity (EA). Geometry optimizations for all the clusters are carried out without imposing any symmetry restriction. The neutral clusters having up to six atoms prefer planar structure and three-dimensional structure is preferred only when the number of atoms in a cluster is more than six. There is a good agreement between the present theoretical and reported experimental IP values for the neutral clusters with cluster size n相似文献   

Electron scattering on OH-(H2O)n clusters (n = 0-4)

The cross sections for electron scattering on OH-(H2O)n for n = 0-4 were measured from threshold to approximately 50 eV. All detachment cross sections were found to follow the classical prediction given earlier [Phys. Rev. Lett. 74, 892 (1995)] with a threshold energy for electron-impact detachment that increased upon sequential hydration, yielding values in the range from 4.5 eV +/- 0.2 eV for OH- to 12.10 eV +/- 0.5 eV for OH-(H2O)4. For n > or = 1, we found that approximately 80% of the total reaction events lead to electron detachment plus total dissociation of the clusters into the constituent molecules of OH and H2O. Finally, we observed resonances in the cross sections for OH-(H2O)3 and for OH-(H2O)4. The resonances were located at approximately 15 eV and were ascribed to the formation of dianions in excited states.     

Collision-induced dissociation studies of protonated ether--(H2O)n (n = 1-3) clusters

We have studied the protonated ether-(H2O)n (n = 1-3) complexes containing tetrahydrofuran, dimethyl, diethyl, dibutyl, and butylmethyl ethers using a flowing afterglow triple-quadrupole mass spectrometer. Collision-induced dissociation, CID, of all clusters with n = 1, 2 shows sequential water loss. The n = 3 cluster of dimethyl ether shows sequential water loss, while all other ether clusters display selective product formation. The CID spectra are interpreted based on known energetics, and theoretical studies of the dimethyl and diethyl ether systems.     

In search of CS2(H2O)(n=1-4) clusters

Gaussian-3 and MP2/aug-cc-pVnZ methods have been used to calculate geometries and thermochemistry of CS(2)(H2O)n, where n=1-4. An extensive molecular dynamics search followed by optimization using these two methods located two dimers, six trimers, six tetramers, and two pentamers. The MP2/aug-cc-pVDZ structure matched best with the experimental result for the CS(2)(H2O) dimer, showing that diffuse functions are necessary to model the interactions found in this complex. For larger CS(2)(H2O)n clusters, the MP2/aug-cc-pVDZ minima are significantly different from the MP2(full)6-31G* structures, revealing that the G3 model chemistry is not suitable for investigation of sulfur containing van der Waals complexes. Based on the MP2/aug-cc-pVTZ free energies, the concentration of saturated water in the atmosphere and the average amount of CS(2) in the atmosphere, the concentrations of these clusters are predicted to be on the order of 10(5) CS(2)(H2O) clusters.cm(-3) and 10(2) CS(2)(H2O)(2) clusters.cm(-3) at 298.15 K. The MP2/aug-cc-pVDZ scaled harmonic and anharmonic frequencies of the most abundant dimer cluster at 298 K are presented, along with the MP2/aug-cc-pVDZ scaled harmonic frequencies for the CS(2)(H(2)O)(n) structures predicted to be present in a low-temperature molecular beam experiment.     

A density functional study on beryllium-doped carbon dianion clusters CnBe2- (n = 4-14)

Making use of the software of molecular graphics, we designed numerous models of C(n)()Be(2-) (n = 4-14). We carried out geometry optimization and calculation on vibration frequency by means of the B3LYP density functional method. After comparison of structure stability, we found that the ground-state isomers of C(n)()Be(2-) (n = 4-14) are linear with the beryllium atom located inside the C(n)() chain. When a side carbon chain is with an even number of carbon atoms, it is polyacetylene-like, whereas when a side chain is with an odd number of carbon atoms, it is cumulene-like. The C(n)Be(2-) (n = 4-14) clusters with an even number of carbon atoms are more stable than that with an odd number of carbon atoms, matching the peak pattern observed in accelerator mass spectrometry (AMS) and Coulomb Explosion Imaging (CEI) investigations of C(n)()Be(2-) (n = 4-14). The trend of such odd/even alternation is explained based on concepts of bonding characteristics, electronic configuration, electron detachment, and incremental binding energy.     

7-羟基喹啉-(H2O)n(n=1,2)复合物氢键相互作用的密度泛函理论研究

许多生理过程都通过分子间相互作用来实现。氢键则是最基本的化学作用力之一。具有碱性和酸性双官能团的芳香族化合物能与水作用形成氢键网络,对于实现生物体系的物质转移（质子转移、离子转移）起着十分重要的作用。在非水溶剂中,通过氢键发生质子转移反应动力学实验特征也己进行了广泛的研究。本文用密度泛函B3LYP方法在6-311G^＊基组水平上对7-羟基喹啉-水复合物相互作用进行了研究,从成键特征及氢键复合物的稳定关系方面进行了理论计算。     

Growth and structural properties of Mg(N) (N = 10-56) clusters: density functional theory study

Using the minima hopping global geometry optimization method on density functional potential energy surface, we have studied the structural and electronic properties of magnesium clusters for a size range of Mg(N) where N = 10-56. Our exhaustive search reveals that most of our global minima are nonsymmetric in the size range above N = 20. We elucidate the evolutionary trend of the entire series and present more details about the peculiar growth of the clusters. For N > 20, it is possible to divide the cluster into two regions: the core region and the surface region. It turns out that the growth follows a peculiar cyclic pattern where the core and surface grow alternatively. The surface energy, as a function of number of atoms shows a clear signature as the number of atoms in the core increases by one. We have also carried out stability analysis and the stable sizes(magic numbers) agree very well with the experimental magic numbers reported by Diederich [J. Chem. Phys. 2011, 134, 124302]. We point out the similarities and differences between our results and sodium clusters.     

Theoretical study on (Al2O3)n (n = 1-10 and 30) fullerenes and H2 adsorption properties

The structures and stabilities of (Al2O3)n (n = 1-10 and 30) clusters were studied by means of first principles calculations. The calculated results reveal that the global minima of small (Al2O3)n (n = 1-5) clusters are cage structures with high symmetries, in which Al and O atoms are three- and two-coordinated, respectively, and are linked to neighbors via single bonds. Beyond (Al2O3)5, we calculated both cage and cage-dimer structures for (Al2O3)n (n = 6-10), and the results show that, at this size range, cage-dimer structures are more stable than cage structures. Furthermore, an onion-like motif for (Al2O3)10 was studied, and it is interesting to find that, at this size, the onion structure is more favorable than cage and cage-dimer structures. For large clusters, a shell structure of Al60O90 is suggested. Electronic properties and calculations on hydrogen adsorption of these aluminum oxide structures are reported, and we discuss their possible use as hydrogen storage materials.     

Structural and electronic properties of TaSi(n) (n=1-13) clusters: a relativistic density functional investigation

The TaSi(n) (n=1-13) clusters with doublet, quartet, and sextet spin configurations have been systematically investigated by a relativistic density functional theory with the generalized gradient approximation available in Amsterdam density functional program. The total bonding energies, equilibrium geometries, Mulliken populations as well as Hirshfeld charges of TaSi(n) (n=1-13) clusters are calculated and presented. The emphasis on the stabilities and electronic properties is discussed. The most stable structures of the small TaSi(n) (n=1-6) clusters and the evolutional rule of low-lying geometries of the larger TaSi(n) (n=7-13) clusters are obtained. Theoretical results indicate that the most stable structure of TaSi(n) (n=1-6) clusters keeps the similar framework as the most stable structure of Si(n+1) clusters except for TaSi(3) cluster. The Ta atom in the lowest-energy TaSi(n) (n=1-13) isomers occupies a gradual sinking site, and the site moves from convex, to flatness, and to concave with the number of Si atom varying from 1 to 13. When n=12, the Ta atom in TaSi(12) cluster completely falls into the center of the Si frame, and a cagelike TaSi(12) geometry is formed. Meanwhile, the net Mulliken and Hirsheld populations of the Ta atom in the TaSi(n) (n=1-13) clusters vary from positive to negative, manifesting that the charges in TaSi(n) (n>/=12) clusters transfer from Si atoms to Ta atom. Additionally, the contribution of Si-Si and Si-Ta interactions to the stability of TaSi(n) clusters is briefly discussed. Furthermore, the investigations on atomic averaged binding energies and fragmentation energies show that the TaSi(n) (n=2,3,5,7,10,11,12) clusters have enhanced stabilities. Compared with pure silicon clusters, a universal narrowing of highest occupied molecular orbital-lowest unoccupied molecular orbital gap in TaSi(n) clusters is found.