Far‐Infrared Spectra of Yttrium‐Doped Gold Clusters AunY (n=1–9)

The geometric, spectroscopic, and electronic properties of neutral yttrium‐doped gold clusters Au_nY (n=1–9) are studied by far‐infrared multiple photon dissociation (FIR‐MPD) spectroscopy and quantum chemical calculations. Comparison of the observed and calculated vibrational spectra allows the structures of the isomers present in the molecular beam to be determined. Most of the isomers for which the IR spectra agree best with experiment are calculated to be the energetically most stable ones. Attachment of xenon to the Au_nY cluster can cause changes in the IR spectra, which involve band shifts and band splittings. In some cases symmetry changes, as a result of the attachment of xenon atoms, were also observed. All the Au_nY clusters considered prefer a low spin state. In contrast to pure gold clusters, which exhibit exclusively planar lowest‐energy structures for small sizes, several of the studied species are three‐dimensional. This is particularly the case for Au₄Y and Au₉Y, while for some other sizes (n=5, 8) the 3D structures have an energy similar to that of their 2D counterparts. Several of the lowest‐energy structures are quasi‐2D, that is, slightly distorted from planar shapes. For all the studied species the Y atom prefers high coordination, which is different from other metal dopants in gold clusters.     

Structures and electronic properties of the small rubidium‐doped silicon RbSin (n = 1–12) clusters

The geometries, relative stabilities, and electronic properties of small rubidium‐doped silicon clusters RbSi_n (n = 1–12) have been systematically investigated using the density functional theory at the B3LYP/GENECP level. The optimized structures show that lowest‐energy isomers of RbSi_n are similar with the ground state isomers of pure Si_n clusters and prefer the three‐dimensional for n = 3–12. The relative stabilities of RbSi_n clusters have been analyzed on the averaged binding energy, fragmentation energy, second‐order energy difference, and highest occupied molecular orbital‐lowest unoccupied molecular orbital energy gap. The calculated results indicate that the doping of Rb atom enhances the chemical activity of Si_n frame and the magic number is RbSi₂. The Mulliken population analysis reveals that the charges in the corresponding RbSi_n clusters transfer from the Rb atom to Si atoms. The partial density of states and chemical hardness are also discussed. © 2014 Wiley Periodicals, Inc.     

Li2B12 and Li3B12: Prediction of the Smallest Tubular and Cage‐like Boron Structures

An intriguing structural transition from the quasi‐planar form of B₁₂ cluster upon the interaction with lithium atoms is reported. High‐level computations show that the lowest energy structures of LiB₁₂, Li₂B₁₂, and Li₃B₁₂ have quasi‐planar (C_s), tubular (D_6d), and cage‐like (C_s) geometries, respectively. The energetic cost of distorting the B₁₂ quasi‐planar fragment is overcompensated by an enhanced electrostatic interaction between the Li cations and the tubular or cage‐like B₁₂ fragments, which is the main reason of such drastic structural changes, resulting in the smallest tubular (Li₂B₁₂) and cage‐like (Li₃B₁₂) boron structures reported to date.     

Stability,Infrared Spectra and Electronic Structures of (ZrO2)n(n=3-6)Clusters:DFT Study

The stability, infrared spectra and electronic structures of (ZrO₂)_n (n=3–6) clusters have been investigated by using density‐functional theory (DFT) at B3LYP/6‐31G* level. The lowest‐energy structures have been recognized by considering a number of structural isomers for each cluster size. It is found that the lowest‐energy (ZrO₂)₅ cluster is the most stable among the (ZrO₂)_n (n=3–6) clusters. The vibration spectra of Zr? O stretching motion from terminal oxygen atom locate between 900 and 1000 cm^?1, and the vibrational band of Zr? O? Zr? O four member ring is obtained at 600–700 cm^?1, which are in good agreement with the experimental results. Mulliken populations and NBO charges of (ZrO₂)_n clusters indicate that the charge transfers occur between 4d orbital of Zr atoms and 2p orbital of O atoms. HOMO‐LUMO gaps illustrate that chemical stabilities of the lowest‐energy (ZrO₂)_n (n=3–6) clusters display an even‐odd alternating pattern with increasing cluster size.     

Density functional study of structural and electronic properties of AlnN (1 ≤ n ≤ 12) clusters

Low‐lying equilibrium geometric structures of Al_nN (n = 1–12) clusters obtained by an all‐electron linear combination of atomic orbital approach, within spin‐polarized density functional theory, are reported. The binding energy, dissociation energy, and stability of these clusters are studied within the local spin density approximation (LSDA) and the three‐parameter hybrid generalized gradient approximation (GGA) due to Becke–Lee–Yang–Parr (B3LYP). Ionization potentials, electron affinities, hardness, and static dipole polarizabilities are calculated for the ground‐state structures within the GGA. It is observed that symmetric structures with the nitrogen atom occupying the internal position are lowest‐energy geometries. Generalized gradient approximation extends bond lengths as compared with the LSDA lengths. The odd–even oscillations in the dissociation energy, the second differences in energy, the highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO–LUMO) gaps, the ionization potential, the electron affinity, and the hardness are more pronounced within the GGA. The stability analysis based on the energies clearly shows the Al₇N cluster to be endowed with special stability. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

The Evolution of Bonding and Thermodynamic Properties of Boron‐Doped Small Carbon Clusters: An Ab Initio Study

Theoretical studies on BC_n (n=1–6) clusters are carried out using density functional theory, Møller–Plesset second‐order perturbation theory (MP2), coupled‐cluster calculations including up to triple excitations (CCSD(T)), and higher‐level approaches. All possible isomers depending on the positions of the boron atom are generated and the lowest‐energy isomers are determined for doublet and quartet electronic states. The three potential evolution paths of the clusters are determined as a function of their size. The energetic and electronic consequences for the increased size of structures differ significantly, which leads to representatives of the ground electronic state from different structural groups. The ab initio calculated thermal functions allow enhancements to the available atomization energies and improve the agreement between the calculated and experimental heat content.     

Structure and stability of coinage metal fluoride and chloride clusters (MnFn and MnCln,M = Cu,Ag, or Au; n = 1–6)

Calculations in the framework of the density functional theory are performed to study the lowest‐energy isomers of coinage metal fluoride and chloride clusters (M_nF_n, M_nCl_n, M = Cu, Ag, or Au, n = 1–6). For all calculated species starting from the trimers the most stable structures are found to be cyclic arrangements. However, planar rings are favored in the case of metal fluorides whereas metal chlorides prefer nonplanar cycles. Calculated bond lengths and infrared frequencies are compared with the available experimental data. The nature of the bonding, involving both covalent and ionic contributions, is characterized. The stability and the fragmentation are also investigated. Trimers are found to be particularly stable when considering the Gibbs free energies. © 2012 Wiley Periodicals, Inc.     

硼碳团簇BnC2 (n＝1～6)的理论研究

应用密度泛函理论在B3LYP/6-311＋G(d)水平上研究了硼碳团簇B_nC₂ (n＝1～6)的几何结构、生长机制和相对稳定性. 计算结果表明, 对于n＝2～6的簇, 平面多环状构型为最稳定的结构, 其中C原子分布于环的顶点、有尽可能多的三配位硼原子和尽可能多的B—C键. 碳原子作为杂原子倾向掺杂于团簇的顶点位置, 它的掺杂不改变硼团簇的主体结构. 与平面多环状结构相比, 随着簇尺寸的增大, 三维结构和线性链结构更不稳定. 在低能线性结构中, C原子位于链两侧的第二个位置. 计算的碎片分裂能、递增键能以及HOMO-LUMO能隙表明, B₄C₂为幻数簇.     

Optoelectronic properties of benzotrithiophene isomers: A density functional theory study

Benzotrithiophene (BTT) isomers were investigated using density functional theory (DFT) and time‐dependent DFT (TD‐DFT) with the aim to explore their structures, linear optical properties, vertical and adiabatic ionization potentials (IP_v and IP_a), electron affinities (EA_v and EA_a), and reorganization energies (λ). The computed bond lengths and bond angles at the B3LYP/6–311+G (d, p) level of theory are in good agreement with experimental crystal structures of the known BTTs. These molecules are planar with zero dihedral angle, making them an ideal backbone for high charge mobility. The UV–visible spectra of BTT isomers are in the range 280–360 nm. All BTT isomers have low hole/electron reorganization energies, which is the main characteristic of good hole/electron transporting materials, and these isomers in turn have potential applications in the field of organic materials.     

The effects of atom pyramidalization and square distribution on the stability of F4F6‐(BN)n polyhedrons

The structures and stability of F₄F₆‐(BN)_n polyhedrons (n = 20–30) with the alternation of B and N atoms were studied with DFT method. The calculation results reveal that the atoms at square–square fusions with large pyramidalization angles are remarkably extruded out of the surfaces of (BN)_n polyhedrons. The energetically favorable isomers do not contain square–square bonds and the energies of those isomers containing square–square bonds increase with the number of square–square bonds linearly, demonstrating that the energetically favorable structures of F₄F₆‐(BN)_n polyhedrons satisfy the isolated square rule and square adjacency penalty rule. The atom pyramidalization determines the stability of the isomers. The binding energy is fitted to the numbers of vertices formed from different faces and a model is proposed to predict the relative stability of these polyhedral molecules. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Double icosahedron‐based motif of Nin (n = 20−30)

A genetic algorithm (GA) coupled with a tight‐binding (TB) interatomic potential was used to search for the low‐energy structures of the medium‐sized Ni_n (n = 20?30) clusters. The low‐energy candidate structures from the GA/TB search were further optimized by using the density functional theory calculations with the Perdew, Burke, and Ernzerhof exchange‐correlation energy functional. The obtained lowest‐energy structures of the medium‐sized Ni_n (n = 20?30) clusters are shown to exhibit double icosahedron‐based motif. The properties of the nickel clusters including binding energies, second differences in energy, and especially magnetic properties have also been studied. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Excited state properties,fluorescence energies,and lifetimes of a poly(fluorene‐phenylene), based on TD‐DFT investigation

The structural and electronic properties of fluorene‐phenylene copolymer (FP)_n, n = 1–4 were studied by means of quantum chemical calculations based on density functional theory (DFT) and time dependent density functional theory (TD‐DFT) using B3LYP functional. Geometry optimizations of these oligomers were performed for the ground state and the lowest singlet excited state. It was found that (FP)_n is nonplanar in its ground state while the electronic excitations lead to planarity in its S₁ state. Absorption and fluorescence energies were calculated using TD‐B3LYP/SVP and TD‐B3LYP/SVP+ methods. Vertical excitation energies and fluorescence energies were obtained by extrapolating these values to infinite chain length, resulting in extrapolated values for vertical excitation energy of 2.89 and 2.87 eV, respectively. The S₁ ← S₀ electronic excitation is characterized as a highest occupied molecular orbital to lowest unoccupied molecular orbital transition and is distinguishing in terms of oscillator strength. Fluorescence energies of (FP)_n calculated from TD‐B3LYP/SVP and TD‐B3LYP/SVP+ methods are 2.27 and 2.26 eV, respectively. Radiative lifetimes are predicted to be 0.55 and 0.51 ns for TD‐B3LYP/SVP and TD‐B3LYP/SVP+ calculations, respectively. These fundamental information are valuable data in designing and making of promising materials for LED materials. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Electronic structure and stability of BP clusters: theoretical calculations for (BP)n (n = 2–4)

The geometry, electronic configurations, harmonic vibrational frequencies, and stability of the structural isomers of boron phosphide clusters have been investigated using density functional theory (DFT). CCSD(T) calculations show that the lowest‐energy structures are cyclic (IIt, IVs) with D_nh symmetry for dimers and trimers. The caged structure for B₄P₄ lie higher in energy than the monocyclic structure with D_2d symmetry (VIs). The B–P bond dominates the structures for many isomers, so that one preferred dissociation channel is loss of the BP monomer. The hybridization and chemical bonding in the different structures are also discussed. Comparisons with boron nitride clusters, the ground state structures of B_nP_n (n = 2, 3) clusters are analogous to those of their corresponding B_nN_n (n = 2, 3) counterparts. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Optimization of bimetallic Cu–Au and Ag–Au clusters by using a modified adaptive immune optimization algorithm

A modified adaptive immune optimization algorithm (AIOA) is designed for optimization of Cu–Au and Ag–Au bimetallic clusters with Gupta potential. Compared with homoatom clusters, there are homotopic isomers in bimetallic cluster, so atom exchange operation is presented in the modified AIOA. The efficiency of the algorithm is tested by optimization of Cu_nAu_38‐n (0 ≤ n ≤ 38). Results show that all the structures with the putative global minimal energies are successfully located. In the optimization of Ag_nAu_55‐n (0 ≤ n ≤ 55) bimetallic clusters, all the structures with the reported minimal energies are obtained, and 36 structures with even lower potential energies are found. On the other hand, with the optimized structures of Cu_nAu_55‐n, it is shown that all 55‐atom Cu–Au bimetallic clusters are Mackay icosahedra except for Au₅₅, which is a face‐centered cubic (fcc)‐like structure; Cu₅₅, Cu₁₂Au₄₃, and Cu₁Au₅₄ have two‐shell Mackay icosahedral geometries with I_h point group symmetry. © 2009 Wiley Periodicals, Inc. J Comput Chem 2009     

A conformational analysis of the six isomers of oxybispyridine

Conformational energy contour maps of the six isomers of oxybispyridine have been constructed using the ab initio STO-3G molecular orbital method. The calculations (employing a partial rigid rotor) for all six isomers indicate that the minimum energy conformers are not planar and that energy barriers between 70–1000 kJ mol^?1 restrict interconversion to planar structures, thereby preventing conjugation between the p-electrons of the oxygen atom with the π system of the pyridine rings. It is postulated that of the three mechanisms used to explain conformer interconversion about the C? O bond, the disrotatory one-ring flip mechanism is the most appropriate, since the “Morino's” structures are all within 2.5 kJ mol^?1 of the minimum. Furthermore, room temperature accessibility of the “Morino's” structures suggests that the Smiles rearrangement would be possible for suitably substituted derivatives of these isomers.     

Search for the global minimum structures of AlB3H2n (n = 0 − 6) clusters

The global minimum structures of AlB₃H_2n (n = 0–6) clusters are determined using the stochastic search method at the B3LYP/6–31G level of theory. These initially specified geometries are recalculated using B3LYP and CCSD(T) methods using the 6–311++G^** basis set. The structural and electronic properties of the two lowest‐lying isomers are presented. The structural parameters obtained for aluminum borohydride are compared with the experimental and theoretical results. The H₂ fragmentation energies of the most stable isomers are investigated. Chemical bonding analyses for the global minimum of AlB₃H_2n (n = 0–6) clusters are performed using the adaptive natural density partitioning method. © 2014 Wiley Periodicals, Inc.     

Structural,spectroscopic aspects,and electronic properties of (TiO2)n clusters: A study based on the use of natural algorithms in association with quantum chemical methods

In this article, we propose a stochastic search‐based method, namely genetic algorithm (GA) and simulated annealing (SA) in conjunction with density functional theory (DFT) to evaluate global and local minimum structures of (TiO₂)_n clusters with n = 1–12. Once the structures are established, we evaluate the infrared spectroscopic modes, cluster formation energy, vertical excitation energy, vertical ionization potential, vertical electron affinity, highest occupied molecular orbital (HOMO)‐lowest unoccupied molecular orbital (LUMO) gaps, and so forth. We show that an initial determination of structure using stochastic techniques (GA/SA), also popularly known as natural algorithms as their working principle mimics certain natural processes, and following it up with density functional calculations lead to high‐quality structures for these systems. We have shown that the clusters tend to form three‐dimensional networks. We compare our results with the available experimental and theoretical results. The results obtained from SA/GA‐DFT technique agree well with available theoretical and experimental data of literature. © 2013 Wiley Periodicals, Inc.     

Structural and static electric response properties of highly symmetric lithiated silicon cages: Theoretical predictions

It is shown by density functional theory calculations that high symmetry silicon cages can be designed by coating with Li atoms. The resulting highly symmetric lithiated silicon cages (up to D_5d symmetry) are low‐lying true minima of the energy hypersurface with binding energies of the order of 4.6 eV per Si atom and moderate highest occupied molecular orbital–lowest unoccupied molecular orbital gaps. Moreover, relying on a systematic study of the electric response properties obtained by ab initio (Hartree–Fock, MP2, and configuration interaction singles (CIS)) and density functional (B3LYP, B2PLYP, and CAM‐B3LYP) methods, it is shown that lithium coating has a large impact on the magnitude of their second hyperpolarizabilities resulting to highly hyperpolarizable species. Such hyperpolarizable character is directly connected to the increase in the density of the low‐lying excited states triggered by the interaction between the Si cage and the surrounding Li atoms. © 2012 Wiley Periodicals, Inc.     

Effect of cage size on the third‐order optical properties of endohedral metallofullerenes Sc3N@C2n (2n = 68, 70, 78, and 80): A theoretical study

Geometrical structures of the investigated endohedral metallofullerenes Sc₃N@C_2n (2n = 68, 70, 78, and 80) were optimized at the B3LYP/6‐31G* level. The analyses of electronic structures display that the contribution of fullerene cage to the lowest unoccupied molecular orbital decreases as the cage size increases. Based on the optimized structures, the time‐dependent density functional theory combined with the sum‐over‐states method was used to investigate their nonlinear optical properties. Calculated third‐order polarizabilities γ and two‐photon absorption (TPA) cross‐section δ do not present the monotone variation with the size of fullerene cage, with largest γ of 0.48 × 10^?34 esu for Sc₃N@C₇₈ in static state, and largest δ of 12.374 GM for Sc₃N@C₇₀ in the wavelength of 902.5 nm. However, the obtained TPA resonant peaks shift red with the size of fullerene cage. By analyzing the electronic origin of the third‐order optical properties, it is found that the charge transfers from the fullerene cage to the encapsulated Sc₃N cluster make important contributions to the studied properties. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011