Ar n HF van der Waals clusters revisited. I. New low-energy isomeric structures for n=6-13

New low-lying isomeric structures of Ar(n)HF clusters are reported for n=6-13. They were determined using simulated annealing and evolutionary programming, for pairwise additive intermolecular potential energy surfaces. New global minima were found for the clusters with n=7, 10, 11. The new lowest-energy structure of Ar(7)HF and several new local minima for n=6, 7 clusters have the HF bound on a threefold surface site, consistent with the recent spectroscopic data for Ar(n)HF clusters in helium nanodroplets. A new type of low-energy local minima were determined for n=9-13 clusters.     

Evolution of the properties of Al(n)N(n) clusters with size

A global optimization of stoichiometric (AlN)(n) clusters (n = 1-25, 30, 35, ..., 95, 100) has been performed using the basin-hopping (BH) method and describing the interactions with simple and yet realistic interatomic potentials. The results for the smaller isomers agree with those of previous electronic structure calculations, thus validating the present scheme. The lowest-energy isomers found can be classified in three different categories according to their structural motifs: (i) small clusters (n = 2-5), with planar ring structures and 2-fold coordination, (ii) medium clusters (n = 6-40), where a competition between stacked rings and globular-like empty cages exists, and (iii) large clusters (n > 40), large enough to mix different elements of the previous stage. All the atoms in small and medium-sized clusters are in the surface, while large clusters start to display interior atoms. Large clusters display a competition between tetrahedral and octahedral-like features: the former lead to a lower energy interior in the cluster, while the latter allow for surface terminations with a lower energy. All of the properties studied present different regimes according to the above classification. It is of particular interest that the local properties of the interior atoms do converge to the bulk limit. The isomers with n = 6 and 12 are specially stable with respect to the gain or loss of AlN molecules.     

A DFT-Based Study of the Low-Energy Electronic Structures and Properties of Small Gold Clusters

Gold clusters Au_n of size n = 2–12 atoms were studied by the density-functional theory with an ab-initio pseudopotential and a generalized gradient approximation. Geometry optimizations starting from a number of initial candidate geometries were performed for each cluster size, so as to determine a number of possible low-energy isomers for each size. Along with the lowest-energy structures, metastable structures were obtained for many cluster sizes. Interestingly, a metastable planar zigzag arrangement of Au atoms was obtained for every cluster size n ≥ 5. The stable electronic structure, binding energy, relative stability and HOMO–LUMO gap for the lowest-energy isomer were calculated for each cluster size. Variation of the electronic properties with size is investigated in this paper and compared with experimental results and other calculations.     

Density-functional study of structural and electronic properties of Si(n)C(n) (n=1-10) clusters

Density-functional theory with generalized gradient approximation for the exchange-correlation potential has been used to calculate the structural and electronic structure of Si(n)C(n) (n=1-10) clusters. The geometries are found to undergo a structural change from two dimensional to three dimensional when the cluster size n equals 4. Cagelike structures are favored as the cluster size increases. A distinct segregation between the silicon and carbon atoms is observed for these clusters. It is found that the C atoms favor to form five-membered rings as the cluster size n increases. However, the growth motif for Si atoms is not observed. The Si(n)C(n) clusters at n=2, 6, and 9 are found to possess relatively higher stability. On the basis of the lowest-energy geometries obtained, the size dependence of cluster properties such as binding energy, HOMO-LUMO gap, Mulliken charge, vibrational spectrum, and ionization potential has been computed and analyzed. The bonding characteristics of the clusters are discussed.     

Ab initio study of structural stability of Mo-S clusters and size specific stoichiometries of magic clusters

Using first-principles calculations with ultrasoft pseudopotential formalism and the generalized gradient approximation for the exchange-correlation functional, we study the stability of MonSm (n =1-6 and m ranging from n to 3n) clusters and obtain the optimal stoichiometry for each n corresponding to the magic cluster. It is found that in this size range, the lowest-energy structures favor a core of metal atoms, which is covered by sulfur. In particular, we observe that for Mo6S14 isolated clusters, a 3D structure is significantly lower in energy as compared to platelet structures found recently on Au (111) surface. The composition ratio between S and Mo in the magic clusters is less than 2 for n=3 and greater than 2 for n<3. The structural stability of the magic clusters arises from the optimization of the Mo-Mo and Mo- S bonding as well as the symmetry of the cluster. Addition of a terminal sulfur in a magic cluster generally lowers its binding energy. The presence of partially occupied d-orbitals in Mo atoms contributes to Mo-Mo bonding and for higher S concentration it leads to sulfur-sulfur bond formation. The variation in energy due to a change in the sulfur composition suggests that sulfurization of the magic clusters is generally more favorable than desulfurization.     

Ar(n)HF van der Waals clusters revisited: II. Energetics and HF vibrational frequency shifts from diffusion Monte Carlo calculations on additive and nonadditive potential-energy surfaces for n=1-12

The ground-state energies and HF vibrational frequency shifts of Ar(n)HF clusters have been calculated on the nonadditive potential-energy surfaces (PESs) for n=2-7 and on the pairwise-additive PESs for the clusters with n=1-12, using the diffusion Monte Carlo (DMC) method. For n>3, the calculations have been performed for the lowest-energy isomer and several higher-lying isomers which are the closest in energy. They provide information about the isomer dependence of the HF redshift, and enable direct comparison with the experimental data recently obtained in helium nanodroplets. The agreement between theory and experiment is excellent, in particular, for the nonadditive DMC redshifts. The relative, incremental redshifts are reproduced accurately even at the lower level of theory, i.e., the DMC and quantum five-dimensional (rigid Ar(n)) calculations on the pairwise-additive PESs. The nonadditive interactions make a significant contribution to the frequency shift, on the order of 10%-12%, and have to be included in the PESs in order for the theory to yield accurate magnitude of the HF redshift. The energy gaps between the DMC ground states of the cluster isomers are very different from the energy separation of their respective minima on the PES, due to the considerable variations in the intermolecular zero-point energy of different Ar(n)HF isomers.     

Structures and stability of lithium monosilicide clusters SiLin (n = 4–16): What is the maximum number,magic number,and core number for lithium coordination to silicon?

In the coordination, hypervalent and cluster chemistry, three important characteristic properties are the maximum coordination number, magic number, and core coordination number. Yet, few studies have considered these three numbers at the same time for an ML(n) cluster with n larger than 8. In this article, we systematically studied the three properties of SiLi(n) (n = 4-16) clusters at the B3LYP/6-31G(d), B3LYP/6-311++G(2d), and CCSD(T)/6-311++G(3df)//B3LYP/6-311++G(2d) (for energy only) levels. Various isomeric forms with different symmetries were calculated. For each SiLi(n) (n = 4-9), silicon cohesive energy (cE) from SiLi(n) --> Si + Li(n) reaction, vertical ionization potential (vIP), and vertical electron affinity (vEA) were obtained for the lowest-energy isomer. We found that the maximum Li-coordination number of Si is 9, which is the largest number among the known MLi(n) clusters. All cE, vIP, and vEA values predicted that 6 is the magic Li-coordination number of Si. For small SiLi(n) (n < or = 6) clusters, Li atoms favor direct coordination to Si, whereas for larger SiLi(n) (n > or = 7) clusters, there is a core cluster that is surrounded by excessive Li atoms. The core Li-coordination number is 6 for SiLi(n) (n = 7,8), 7 for SiLi(n) (n = 9,10), 8 for SiLi(n) (n = 11-15) and 9 for SiLi(n) (n > or = 16). Through the calculations, we verified the relationship between the structure and stability of SiLi(n) with the maximum coordination number, magic number, and core coordination number.     

Structures and dissociation channels of protonated mixed clusters around a small magic number: infrared spectroscopy of ((CH3)3N)n-H(+)-H2O (n = 1-3)

The magic number behavior of ((CH(3))(3)N)(n)-H(+)-H(2)O clusters at n = 3 is investigated by applying infrared spectroscopy to the clusters of n = 1-3. Structures of these clusters are determined in conjunction with density functional theory calculations. Dissociation channels upon infrared excitation are also measured, and their correlation with the cluster structures is examined. It is demonstrated that the magic number cluster has a closed-shell structure, in which the water moiety is surrounded by three (CH(3))(3)N molecules. The ion core (protonated site) of the clusters is found to be (CH(3))(3)NH(+) for n = 1-3, but coexistence of an isomer of the H(3)O(+) ion core cannot be ruled out for n = 3. Large rearrangement of the cluster structures of n = 2 and 3 before dissociation, which has been suggested in the mass spectrometric studies, is confirmed on the basis of the structure determination by infrared spectroscopy.     

Structures of bimetallic clusters

We report an optimization algorithm for studying bimetallic nanoclusters. The algorithm combines two state-of-the-art methods, the genetic algorithm and the basin hopping approach, widely employed in the literature for predicting structures of pure metallic and nonmetallic clusters. To critically test the present algorithm and its use in determining the lowest-energy structures of bimetallic nanoclusters, we apply it to study the bimetallic clusters Cu(n)Au(38-n) (0< or =n< or =38). It is predicted that the Au atoms, being larger in size than the Cu atoms, prefer to occupy surface sites showing thus the segregating behavior. As the atom fraction of Cu increases, the bimetallic cluster Cu(n)Au(38-n), as a whole, first takes on an amorphous structure and is followed by dramatic changes in structure with the Cu atoms revealing hexagonal, then assuming pentagonal, and finally shifting to octahedral symmetry in the Cu-rich range.     

Structural and energetic properties of Ni-Cu bimetallic clusters

The lowest-energy structures for all compositions of Ni n Cu m bimetallic clusters with N = n + m up to 20 atoms, N = 23, and N = 38 atoms have been determined using a genetic algorithm for unbiased structure optimization in combination with an embedded-atom method for the calculation of the total energy for a given structure. Comparing bimetallic clusters with homoatomic clusters of the same size, it is shown that the most stable structures for each cluster size are composed entirely of Ni atoms. Among the bimetallic clusters in the size range N = 2-20, the Ni N-1 Cu 1 clusters possess the highest stability. Further, it has been established that most of the bimetallic cluster structures have geometries similar to those of pure Ni clusters. The size N = 38 presents a special case, as the bimetallic clusters undergo a dramatic structural change with increasing atom fraction of Cu. Moreover, we have identified an icosahedron, a double, and a triple icosahedron with one, two, and three Ni atoms at the centers, respectively, as particularly stable structures. We show that in all global-minimum structures Ni atoms tend to occupy mainly high-coordination inner sites, and we confirm the segregation of Cu on the surface of Ni-Cu bimetallic clusters predicted in previous studies. Finally, it is observed that, in contrast to the bulk, the ground-state structures of the 15-, 16-, and 17-atom bimetallic clusters do not experience a smooth transition between the structures of the pure copper and the pure nickel clusters as a function of the relative number of the two types of atoms. For these sizes, the concentration effect on energy is more important than the geometric one.     

Clusters of glycolic acid with three to six water molecules

Semiempirical, ab initio, and density functional theory calculations are used to locate many low-energy minima on the potential energy surfaces of the CH2OHCOOH-(H2O)n complexes with n = 3,4,5,6. In the clusters with three, four, and five water molecules, the lowest-energy structure consists of a (H2O)n complex, not necessarily of lowest energy, hydrogen bonded to the carboxylic group of the glycolic acid. The lowest-energy structure for n = 6 is similar except that the water hexamer is hydrogen bonded to both the carboxylic and alpha-hydroxyl groups of the acid. In all the lowest-energy clusters, the intramolecular hydrogen bond remains intact in the glycolic acid.     

Electronic and infrared spectroscopy of [benzene-(methanol)(n)](+) (n = 1-6)

The microsolvation structure of the [benzene-(methanol)(n)](+) (n = 1-6) clusters was analyzed by electronic and infrared spectroscopy. For the n = 1 and 2 clusters, further spectroscopic investigation was carried out by Ar atom attachment, which has been know as a useful technique for discriminating isomers of the clusters. The coexistence of multiple isomers was confirmed for the n = 1 and 2 clusters, and remarkably, preferential production of the specific isomers occurred in the Ar attachment. The most stable isomer of the n = 1 cluster was suggested to be of the "on-ring" structure where the nonbonding electrons of the methanol moiety directly interact with the pi orbital of the benzene cation moiety. This is a sharp contrast to [benzene-(H(2)O)(1)](+), exhibiting the "side" structure, where the water moiety is bound to the C-H sites of the benzene cation moiety. The structure of the n = 2 cluster was discussed with the help of density functional theory calculations. Spectral signatures of the intracluster proton-transfer reaction were found for n > or = 5. The intracluster electron-transfer reaction leading to the (methanol)(m)()(+) fragment was also seen upon vibrational and electronic excitation of n > or = 4.     

Isomerization dynamics and thermodynamics of ionic argon clusters

The dynamics and thermodynamics of small Ar(n) (+) clusters, n=3, 6, and 9, are investigated using molecular dynamics (MD) and exchange Monte Carlo (MC) simulations. A diatomic-in-molecule Hamiltonian provides an accurate model for the electronic ground state potential energy surface. The microcanonical caloric curves calculated from MD and MC methods are shown to agree with each other, provided that the rigorous conservation of angular momentum is accounted for in the phase space density of the MC simulations. The previously proposed projective partition of the kinetic energy is used to assist MD simulations in interpreting the cluster dynamics in terms of inertial, internal, and external modes. The thermal behavior is correlated with the nature of the charged core in the cluster by computing a dedicated charge localization order parameter. We also perform systematic quenches to establish a connection with the various isomers. We find that the Ar(3) (+) cluster is very stable in its linear ground state geometry up to about 300 K, and then isomerizes to a T-shaped isomer in which a quasineutral atom lies around a charged dimer. In Ar(6) (+) and Ar(9) (+), the covalent trimer core is solvated by neutral atoms, and the weakly bound solvent shell melts at much lower energies, occasionally leading to a tetramer or pentamer core with weakly charged extremities. At high energies the core itself becomes metastable and the cluster transforms into Ar(2) (+) solvated by a fluid of neutral argon atoms.     

Size resolved infrared spectroscopy of Na(CH3OH)n (n = 4-7) clusters in the OH stretching region: unravelling the interaction of methanol clusters with a sodium atom and the emergence of the solvated electron

Size resolved IR action spectra of neutral sodium doped methanol clusters have been measured using IR excitation modulated photoionisation mass spectroscopy. The Na(CH(3)OH)(n) clusters were generated in a supersonic He seeded expansion of methanol by subsequent Na doping in a pick-up cell. A combined analysis of IR action spectra, IP evolutions and harmonic predictions of IR spectra (using density functional theory) of the most stable structures revealed that for n = 4, 5 structures with an exterior Na atom showing high ionisation potentials (IPs) of ~4 eV dominate, while for n = 6, 7 clusters with lower IPs (~3.2 eV) featuring fully solvated Na atoms and solvated electrons emerge and dominate the IR action spectra. For n = 4 simulations of photoionisation spectra using an ab initio MD approach confirm the dominance of exterior structures and explain the previously reported appearance IP of 3.48 eV by small fractions of clusters with partly solvated Na atoms. Only for this cluster size a shift in the isomer composition with cluster temperature has been observed, which may be related to kinetic stabilisation of less Na solvated clusters at low temperatures. Features of slow fragmentation dynamics of cationic Na(+)(CH(3)OH)(6) clusters have been observed for the photoionisation near the adiabatic limit. This finding points to the relevance of previously proposed non-vertical photoionisation dynamics of this system.     

Ab initio study of neutral and charged SinNap(+) (n <or= 6, p <or= 2) clusters

Ab initio calculations in the framework of the density functional theory, with B3LYP functional, are performed to study the lowest-energy isomers of silicon sodium clusters Si(n)Na(p)(+) (n 相似文献   

From ar clustering dynamics to Ar solvation for Na+-benzene

The gradual evolution from cluster rearrangement to solvation dynamics is discussed by considering the rearrangement of n (n = 1, ..., 19) Ar atoms around Na+-benzene clusters and using an atom-bond potential energy surface. The nature of the bonding is discussed on the basis of the decomposition of the interaction energy and of the formation of the possible conformers. The benzene molecule is found to remain strongly bound to Na+ independently of the number of solvating rare-gas atoms, although due to the anisotropy of the interaction potential, the Ar atoms solvate the Na+-benzene cluster preferentially on the side of the cation. Other specific features of the solvation process are discussed.     

Disparate effects of Cu and V on structures of exohedral transition metal-doped silicon clusters: a combined far-infrared spectroscopic and computational study

The growth mechanisms of small cationic silicon clusters containing up to 11 Si atoms, exohedrally doped by V and Cu atoms, are described. We find that as dopants, V and Cu follow two different paths: while V prefers substitution of a silicon atom in a highly coordinated position of the cationic bare silicon clusters, Cu favors adsorption to the neutral or cationic bare clusters in a lower coordination site. The different behavior of the two transition metals becomes evident in the structures of Si(n)M(+) (n = 4-11 for M = V, and n = 6-11 for M = Cu), which are investigated by density functional theory and, for several sizes, confirmed by comparison with their experimental vibrational spectra. The spectra are measured on the corresponding Si(n)M(+)·Ar complexes, which can be formed for the exohedrally doped silicon clusters. The comparison between experimental and calculated spectra indicates that the BP86 functional is suitable to predict far-infrared spectra of these clusters. In most cases, the calculated infrared spectrum of the lowest-lying isomer fits well with the experiment, even when various isomers and different electronic states are close in energy. However, in a few cases, namely Si(9)Cu(+), Si(11)Cu(+), and Si(10)V(+), the experimentally verified isomers are not the lowest in energy according to the density functional theory calculations, but their structures still follow the described growth mechanism. The different growth patterns of the two series of doped Si clusters reflect the role of the transition metal's 3d orbitals in the binding of the dopant atoms.     

Intracluster anionic oligomerization of acrylic ester molecules initiated by electron transfer from an alkali metal atom

Stabilities and intracluster reactions have been investigated by photoionization mass spectrometry for clusters composed of an alkali metal atom (M; Na and K) and acrylic ester molecules, CH(2)=CHCO(2)R, such as methyl acrylate (MA; R = CH(3)) and ethyl acrylate (EA; R = C(2)H(5)). The following two features are commonly observed in the photoionization mass spectra of M(CH(2)=CHCO(2)R)(n): (1) The ion with n = 3 is clearly observed as a magic number. (2) Fragmented cluster ions with the loss of ROH, [M(CH(2)=CHCO(2)R)(n) - ROH] are detected only for n = 3. These features are both explained by an intracluster oligomerization reaction initiated by electron transfer from the metal atoms. The magic number trimer is concluded to have the stable structure of cyclohexane derivatives as a result of oligomerization. The fragmentation reaction is explained by Dieckmann cyclization after anionic oligomerization to produce another isomer of the trimer. The intracluster electron transfer is also supported by theoretical calculation for Na(MA) based on density functional theory.     

Optimal covering of C(60) fullerene by rare gases

Putative global energy minima of clusters formed by the adsorption of rare gases on a C(60) fullerene molecule, C(60)X(N) (X=Ne, Ar, Kr, Xe; N ≤ 70), are found using basin-hopping global optimization in an empirical potential energy surface. The association energies per rare gas atom as a function of N present two noticeable minima for Ne and Ar and just one for Kr and Xe. The minimum with the smallest N is the deepest one and corresponds to an optimal packing monolayer structure; the other one gives a monolayer with maximum packing. For Kr and Xe, optimal and maximum packing structures coincide. By using an isotropic average form of the X-C(60) interaction, we have established the relevance of the C(60) surface corrugation on the cluster structures. Quantum effects are relevant for Ne clusters. The adsorption of these rare gases on C(60) follows patterns that differ significantly from the ones found recently for He by means of experimental and theoretical methods.     

IR+vacuum ultraviolet (118 nm) nonresonant ionization spectroscopy of methanol monomers and clusters: neutral cluster distribution and size-specific detection of the OH stretch vibrations

Small methanol clusters are formed by expanding a mixture of methanol vapor seeded in helium and are detected using vacuum UV (vuv) (118 nm) single-photon ionization/linear time-of-flight mass spectrometer (TOFMS). Protonated cluster ions, (CH3OH)(n-1)H+ (n=2-8), formed through intracluster ion-molecule reactions following ionization, essentially correlate to the neutral clusters, (CH3OH)n, in the present study using 118 nm light as the ionization source. Both experimental and Born-Haber calculational results clarify that not enough excess energy is released into protonated cluster ions to initiate further fragmentation in the time scale appropriate for linear TOFMS. Size-specific spectra for (CH3OH)n (n=4 to 8) clusters in the OH stretch fundamental region are recorded by IR+vuv (118 nm) nonresonant ion-dip spectroscopy through the detection chain of IR multiphoton predissociation and subsequent vuv single-photon ionization. The general structures and gross features of these cluster spectra are consistent with previous theoretical calculations. The lowest-energy peak contributed to each cluster spectrum is redshifted with increasing cluster size from n=4 to 8, and limits near approximately 3220 cm(-1) in the heptamer and octamer. Moreover, IR+vuv nonresonant ionization detected spectroscopy is employed to study the OH stretch first overtone of the methanol monomer. The rotational temperature of the clusters is estimated to be at least 50 K based on the simulation of the monomer rotational envelope under clustering conditions.