Structure, thermodynamics, and energy content of aluminum-cyclopentadienyl clusters

We present quantum chemistry simulations of aluminum clusters surrounded by a surface layer of cyclopentadiene-type ligands to evaluate the potential of such complexes as novel fuels or energetic materials. Density functional theory simulations are used to examine the aluminum-ligand bonding and its variation as the size of the aluminum cluster increases. The organometallic bond at the surface layer arises mainly from ligand charge donation into the Al p orbitals balanced with repulsive polarization effects. Functionalization of the ligand and changes in Al cluster size are found to alter the relative balance of these effects, but the surface organometallic bond generally remains stronger than Al-Al bonds elsewhere in the cluster. In large clusters, such as the experimentally observed Al(50)Cp(12)*, this suggests that unimolecular thermal decomposition likely proceeds through loss of surface AlCp* units, exposing the strained interior aluminum core. The calculated heats of combustion per unit volume for these systems are high, approaching 60% that of pure aluminum. We discuss the possibility of using organometallic aluminum clusters as a means of achieving rapid combustion in propellants and fuels.     

Origin of the size-dependence of the polarizability per atom in heterogeneous clusters: The case of AlP clusters

An analysis of the atomic polarizabilities α in stoichiometric aluminum phosphide clusters, computed at the MP2 and density functional theory (DFT) levels, the latter using the B3LYP functional, and partitioned using the classic and iterative versions of the Hirshfeld method, is presented. Two sets of clusters are examined: the ground-state Al(n)P(n) clusters (n=2-9) and the prolate clusters (Al(2)P(2))(N) and (Al(3)P(3))(N) (N≤6). In the ground-state clusters, the mean polarizability per atom, i.e., α/2n, decreases with the cluster size but shows peaks at n=5 and at n=7. We demonstrate that these peaks can be explained by a large polarizability of the Al atoms and by a low polarizability of the P atoms in Al(5)P(5) and Al(7)P(7) due to the presence of homopolar bonds in these clusters. We show indeed that the polarizability of an atom within an Al(n)P(n) cluster depends on the cluster size and the heteropolarity of the bonds it forms within the cluster, i.e., on the charges of the atoms. The polarizabilities of the fragments Al(2)P(2) and Al(3)P(3) in the prolate clusters were found to depend mainly on their location within the cluster. Finally, we show that the iterative Hirshfeld method is more suitable than the classic Hirshfeld method for describing the atomic polarizabilities and the atomic charges in clusters with heteropolar bonds, although both versions of the Hirshfeld method lead to similar conclusions.     

Reactions of Al(n)I(x)- with methyl iodide: the enhanced stability of Al7I and the chemical significance of active centers

Al(n)I(x)- are reacted with methyl iodide, and the reaction mechanisms and products are discussed. The relevance of previous studies of the reactions between bare aluminum clusters and methyl iodide is addressed, and the chemical differences reported herein are explained. Particular attention is given to parallels with the known chemistry of alkyl halides on aluminum surfaces, where kinetically mediated etching reactions are prominent. The emergence of Al7I- as the dominant product in the present reactions is addressed via electronic structure calculations, which reveal that the cluster can be described in terms of an electron bound to a "jellium compound". Other significant products of the etching reaction include I-, I3-, and, importantly, the polyhalide-like Al13I2x- clusters. In the Al13I(x)- series, clusters with odd values for x are found to be reactive, and those with even x are far more stable. This observation is explained in terms of the presence or absence of active sites.     

Reactivity of molecular oxygen with aluminum clusters: Density functional and Ab Initio molecular dynamics simulation study

Dissociative adsorption of molecular oxygen (O₂) on aluminum (Al) clusters has attracted much interest in the field of surface science and catalysis, but theoretical predictions of the reactivity of this reaction in terms of barrier height is still challenging. In this regard, we systematically investigate the reactivity of O₂ with Al clusters using density functional theory (DFT) and atom‐centered density matrix propagation (ADMP) simulations. We also calculate potential energy surfaces (PESs) of the reaction between O₂ and Al clusters to estimate the barrier energy of this reaction. The M06‐2X functional gives the barrier energy in agreement with the one calculated by coupled cluster singles and doubles with perturbed triples (CCSD(T)) while the TPSSh functional significantly underestimates the barrier height. The ADMP simulation using the M06‐2X functional predicts the reactivity of O₂ with the Al cluster in agreement with the experimental findings, that is, singlet O₂ readily reacts with Al clusters but triplet O₂ is less reactive. We found that the ability of a DFT functional to describe the charge transfer appropriately is critical for calculating the barrier energy and the reactivity of the reaction of O₂ with Al clusters. The M06‐2X functional is relevant for investigating chemical reactions involving Al and O₂. © 2016 Wiley Periodicals, Inc.     

Hydration of alumina cluster anions in the gas phase

Hydration reactions of anionic aluminum oxide clusters were measured using a quadrupole ion trap secondary ion mass spectrometer, wherein the number of Lewis acid sites were determined. The extent of hydration varied irregularly as cluster size increased and indicated that the clusters possessed condensed structures where the majority of the Al atoms were fully coordinated, with a limited number of undercoordinated sites susceptible to hydrolysis. For maximally hydrated ions, the number of OH groups per Al decreased in an exponential fashion from 4.0 in Al(1) cluster to 1.4 in the Al(9) cluster, which was greater than that expected for a highly hydroxylated surface but less than that for solution phase alumina clusters.     

Al-H bond formation in hydrated aluminum oxide cluster anions

Quantum chemical calculations have been performed to investigate the interaction of a water molecule with gas phase aluminum oxide cluster anions. While oxygen-rich clusters (AlxOy-,xy) generate metal hydrides. These hydride species are, in many cases, 30-35 kcal/mol more stable than their hydroxide counterparts. Our observations on such competing reaction pathways may be useful to understand the catalytic role of alumina nanoparticles in many chemical reactions.     

Hydrogen dissociation on small aluminum clusters

Transition states and reaction paths for a hydrogen molecule dissociating on small aluminum clusters have been calculated using density functional theory. The two lowest spin states have been taken into account for all the Al(n) clusters considered, with n=2-6. The aluminum dimer, which shows a (3)Π(u) electronic ground state, has also been studied at the coupled cluster and configuration interaction level for comparison and to check the accuracy of single determinant calculations in this special case, where two degenerate configurations should be taken into account. The calculated reaction barriers give an explanation of the experimentally observed reactivity of hydrogen on Al clusters of different size [Cox et al., J. Chem. Phys. 84, 4651 (1986)] and reproduce the high observed reactivity of the Al(6) cluster. The electronic structure of the Al(n)-H(2) systems was also systematically investigated in order to determine the role played by interactions of specific molecular orbitals for different nuclear arrangements. Singlet Al(n) clusters (with n even) exhibit the lowest barriers to H(2) dissociation because their highest doubly occupied molecular orbitals allow for a more favorable interaction with the antibonding σ(u) molecular orbital of H(2).     

Formation and properties of halogenated aluminum clusters

The fast-flow tube reaction apparatus was employed to study the halogenation of aluminum clusters. For reactions with HX (X=Cl, Br, and I), acid-etching pathways are evident, and we present findings for several reactions, whereby Al(n)X(-) generation is energetically favorable. Tandem reaction experiments allowed us to establish that for Al(n)Cl(-), Al(n)I(-), and Al(n)I(2) (-), species with n=6, 7, and 15 are particularly resistant to attack by oxygen. Further, trends in reactivity suggest that, in general, iodine incorporation leaves the aluminum clusters' electronic properties largely unperturbed. Ab initio calculations were performed to better interpret reaction mechanisms and elucidate the characteristics of the products. Lowest energy structures for Al(13)X(-) were found to feature icosahedral Al(13) units with the halogen atom located at the on-top site. The charge density of the highest occupied molecular orbital in these clusters is heavily dependent on the identity of X. The dependence of reactivity on the clusters' charge state is also discussed. In addition, we address the enhanced stability of Al(13)I(-) and Al(13)I(2) (-), arguing that the superhalogen behavior of Al(13) in these clusters can provide unique opportunities for the synthesis of novel materials with saltlike structures.     

Catalytic reactions on neutral Rh oxide clusters more efficient than on neutral Rh clusters

Gas phase catalytic reactions involving the reduction of N(2)O and oxidation of CO were observed at the molecular level on isolated neutral rhodium clusters, Rh(n) (n = 10-28), using mass spectrometry. Sequential oxygen transfer reactions, Rh(n)O(m-1) + N(2)O → Rh(n)O(m) + N(2) (m = 1, 2, 3,…), were monitored and the rate constant for each reaction step was determined as a function of the cluster size. Oxygen extraction reactions by a CO molecule, Rh(n)O(m) + CO → Rh(n)O(m-1) + CO(2) (m = 1, 2, 3,…), were also observed when a small amount of CO was mixed with the reactant N(2)O gas. The rate constants of the oxygen extraction reactions by CO for m ≥ 4 were found to be two or three orders of magnitude higher than the rate constants for m ≤ 3, which indicates that the catalytic reaction proceeds more efficiently when the reaction cycles turn over around Rh(n)O(m) (m ≥ 4) than around bare Rh(n). Rhodium clusters operate as more efficient catalysts when they are oxidized than non- or less-oxidized rhodium clusters, which is consistent with theoretical and experimental studies on the catalytic CO oxidation reaction on a rhodium surface.     

Kinetics and mechanisms of CF3CHFOCH3, CF3CHFOC(O)H, and FC(O)OCH3 reactions with OH radicals

The kinetics and mechanism of oxidation of CF3CHFOCH3 was studied using an 11.5-dm3 environmental reaction chamber. OH radicals were produced by UV photolysis of an O3-H2O-He mixture at an initial pressure of 200 Torr in the chamber. The rate constant of the reaction of CF3CHFOCH3 with OH radicals (k1) was determined to be (1.77 +/- 0.69) x 10(-12) exp[(-720 +/- 110)/T] cm3 molecule(-1)(s-1) by means of a relative rate method at 253-328 K. The mechanism of the reaction was investigated by FT-IR spectroscopy at 298 K. CF3CHFOC(O)H, FC(O)OCH3, and COF2 were determined to be the major products. The branching ratio (k1a/k1b) for the reactions CF3CHFOCH3 + OH --> CF3CHFOCH2* + H2O (k1a) and CF3CHFOCH3 + OH --> CF3CF*OCH3 + H2O (k1b) was estimated to be 4.2:1 at 298 K from the yields of CF3CHFOC(O)H, FC(O)OCH3, and COF2. The rate constants of the reactions of CF3CHFOC(O)H (k2) and FC(O)OCH3 (k3) with OH radicals were determined to be (9.14 +/- 2.78) x 10(-13) exp[(-1190 +/- 90)/T] and (2.10 +/- 0.65) x 10(-13) exp[(-630 +/- 90)/T] cm3 molecule(-1)(s-1), respectively, by means of a relative rate method at 253-328 K. The rate constants at 298 K were as follows: k1 = (1.56 +/- 0.06) x 10-13, k2 = (1.67 +/- 0.05) x 10-14, and k3 = (2.53 +/- 0.07) x 10-14 cm3 molecule(-1)(s-1). The tropospheric lifetimes of CF3CHFOCH3, CF3CHFOC(O)H, and FC(O)OCH3 with respect to reaction with OH radicals were estimated to be 0.29, 3.2, and 1.8 years, respectively.     

Gas phase reactivity of thermal metal clusters

Reaction kinetics of metal cluster ions under well defined thermal conditions were studied using a flow tube reactor in combination with laser vaporization. Aluminum anions and cations were reacted with oxygen, and several species which are predicted jellium shell closings, were found to have special stability. Metal alloy cluster anions comprised of Al, V and Nb were also seen to react with oxygen. Alloy clusters with an even number of electrons reacted more slowly than odd electron species, and certain clusters appeared to be exceptionally unreactive. Copper cation clusters were observed to associate with carbon monoxide with reactivities that approach bulk behavior at surprisingly small cluster size. These reactions demonstrate how the rate of reaction changes with cluster size.     

Chemisorption of CO and mechanism of CO oxidation on supported platinum nanoclusters

Kinetic, isotopic, and infrared studies on well-defined dispersed Pt clusters are combined here with first-principle theoretical methods on model cluster surfaces to probe the mechanism and structural requirements for CO oxidation catalysis at conditions typical of its industrial practice. CO oxidation turnover rates and the dynamics and thermodynamics of adsorption-desorption processes on cluster surfaces saturated with chemisorbed CO were measured on 1-20 nm Pt clusters under conditions of strict kinetic control. Turnover rates are proportional to O(2) pressure and inversely proportional to CO pressure, consistent with kinetically relevant irreversible O(2) activation steps on vacant sites present within saturated CO monolayers. These conclusions are consistent with the lack of isotopic scrambling in C(16)O-(18)O(2)-(16)O(2) reactions, and with infrared bands for chemisorbed CO that did not change within a CO pressure range that strongly influenced CO oxidation turnover rates. Density functional theory estimates of rate and equilibrium constants show that the kinetically relevant O(2) activation steps involve direct O(2)* (or O(2)) reactions with CO* to form reactive O*-O-C*=O intermediates that decompose to form CO(2) and chemisorbed O*, instead of unassisted activation steps involving molecular adsorption and subsequent dissociation of O(2). These CO-assisted O(2) dissociation pathways avoid the higher barriers imposed by the spin-forbidden transitions required for unassisted O(2) dissociation on surfaces saturated with chemisorbed CO. Measured rate parameters for CO oxidation were independent of Pt cluster size; these parameters depend on the ratio of rate constants for O(2) reactions with CO* and CO adsorption equilibrium constants, which reflect the respective activation barriers and reaction enthalpies for these two steps. Infrared spectra during isotopic displacement and thermal desorption with (12)CO-(13)CO mixtures showed that the binding, dynamics, and thermodynamics of CO chemisorbed at saturation coverages do not depend on Pt cluster size in a range that strongly affects the coordination of Pt atoms exposed at cluster surfaces. These data and their theoretical and mechanistic interpretations indicate that the remarkable structure insensitivity observed for CO oxidation reactions reflects average CO binding properties that are essentially independent of cluster size. Theoretical estimates of rate and equilibrium constants for surface reactions and CO adsorption show that both parameters increase as the coordination of exposed Pt atoms decreases in Pt(201) cluster surfaces; such compensation dampens but does not eliminate coordination and cluster size effects on measured rate constants. The structural features and intrinsic non-uniformity of cluster surfaces weaken when CO forms saturated monolayers on such surfaces, apparently because surfaces and adsorbates restructure to balance CO surface binding and CO-CO interaction energies.     

Dipole polarizabilities of germanium clusters

Dipole polarizabilities of Ge_n clusters with 2–25 atoms are calculated using finite field (FF) method within density functional theory. The dipole moments and polarizabilities of clusters are sensitively dependent on the cluster geometries and electronic structures. The clusters with low symmetry and large HOMO–LUMO gap prefer to large dipole moments. The polarizabilities of the Ge_n clusters increase rapidly in the size range of 2–5 atoms and then fluctuate around the bulk value. The larger HOMO–LUMO gap may lead to smaller polarizability. As compared with the compact structure and diamond structure, the prolate cluster structure corresponds to a larger polarizability.     

Sequential reactions of silicon clusters with SiD4: constrained heterogeneous nucleation of deuterated silicon particles

The sequential clustering reactions of gas phase silicon cluster ions are ideally suited for studying heterogeneous nucleation of silicon particles. Stepwise reactions of silicon cluster ions with silane have been observed which lead to growth of larger deuterated silicon clusters. Extensive information is obtained about their exothermic gas phase reaction rates and product distributions as a function of cluster size and degree of hydrogenation. What is found is that each size of silicon cluster exhibits a unique pattern of growth. Furthermore, nearly all of the silicon clusters encounter chemical constraints to rapid nucleation. These constraints are a consequence of specific electronic and structural requirements in the chemical reactions. The nature of these requirements are deduced using experimental results in concert with ab initio electronic structure theoretical calculations of the energetics and structures of various species.     

Oxidation of CO by aluminum oxide cluster ions in the gas phase

Small aluminum oxide cluster cations and anions, produced by laser vaporization, were investigated regarding their reactivity toward CO and N2O employing guided-ion-beam mass spectrometry. Clusters with the same stoichiometry as bulk alumina, Al2O3, exhibited atomic oxygen transfer products when reacted with CO, suggesting the formation of CO2. Anionic clusters were less reactive than cations but showed higher selectivity towards the transfer of only a single oxygen atom. Cationic clusters, in contrast, exhibited additional products corresponding to the sequential transfer of two oxygen atoms and the loss of an aluminum atom. To determine if these stoichiometric clusters could be generated from oxygen-deficient species, clusters having a stoichiometry with one less oxygen atom than bulk alumina, Al2O2, were reacted with N2O. Cationic clusters were found to be selectively oxidized to Al2O3(+), while anionic clusters added both one and two oxygen atoms forming Al2O3(-) and Al2O4(-). The oxygen-rich Al2O4(-) cluster exhibited comparable reactivity to Al2O3(-) when reacted with CO.     

Small Ligated Organometallic $${{\hbox {Pd}}}_{n}$$ Clusters ($$n= 4{-}12$$n=4-12): A DFT Investigation

DFT calculations were carried out at the BP86/LANL2DZ level of theory to probe the structures and relative stabilities of small and medium size ligated palladium clusters. Results show that optimized geometries satisfactorily mimic experimental structural data. Observed cluster electron counts generally deviate from the expected Wade–Mingos electron counts because of the presence of several non-conical T-shaped \({\hbox {ML}}_{3}\) fragments in the clusters. These species are thermodynamically stable with substantial HOMO–LUMO gaps.     

Structural characterization of MAO and related aluminum complexes. 1. Solid-state (27)Al NMR with comparison to EFG tensors from ab initio molecular orbital calculations.

Experimental and ab initio molecular orbital techniques are developed for study of aluminum species with large quadrupole coupling constants to test structural models for methylaluminoxanes (MAO). The techniques are applied to nitrogen- and oxygen-containing complexes of aluminum and to solid MAO isolated from active commercial MAO preparations. (Aminato)- and (propanolato)aluminum clusters with 3-, 4-, and 6-coordinate aluminum sites are studied with three (27)Al NMR techniques optimized for large (27)Al quadrupole coupling constants: field-swept, frequency-stepped, and high-field MAS NMR. Four-membered (aminato)aluminum complexes with AlN(4) coordination yield slightly smaller C(q) values than similar AlN(2)C(2) sites: 12.2 vs 15.8 MHz. Planar 3-coordinate AlN(2)C sites have the largest C(q) values, 37 MHz. In all cases, molecular orbital calculations of the electric field gradient tensors yields C(q) and eta values that match with experiment, even for a large hexameric (aminato)aluminum cage. A D(3d) symmetry hexaaluminum oxane cluster, postulated as a model for MAO, yields a calculated C(q) of -23.7 MHz, eta = 0.7474, and predicts a spectrum that is too broad to match the field-swept NMR of methylaluminoxane, which shows at least three sites, all with C(q) values greater than 15 MHz but less than 21 MHz. Thus, the proposed hexaaluminum cluster, with its strained four-membered rings, is not a major component of MAO. However, calculations for dimers of the cage complex, either edge-bridged or face-bridged, show a much closer match to experiment. Also, MAO preparations differ, with a gel form of MAO having significantly larger (27)Al C(q) values than a nongel form, a conclusion reached on the basis of (27)Al NMR line widths in field-swept NMR spectra acquired from 13 to 24 T.     

Dramatic size effects and evidence of structural isomers in the reactions of rhodium clusters, Rhn+/-, with nitrous oxide

The reactions of gas phase rhodium clusters, Rhn+/- (n<30), with nitrous oxide, N2O, have been investigated under single collision conditions by Fourier transform ion cyclotron resonance mass spectrometry. The only significant reaction observed is the sequential generation of oxides. Absolute rate constants for the reactions of all clusters have been determined and, in the case of the cationic clusters especially, they exhibit large fluctuations as a function of cluster size with local minima observed for n=5, 19, 28. Striking similarities are observed with the variation in rate constants for these clusters in reactions with small hydrocarbons (C. Adlhart and E. Uggerud, J. Chem. Phys., 2005, 123, 214709). Corresponding size effects are also observed but are less marked in the reactions of the anionic clusters. The reactions of several clusters exhibit marked deviations from simple pseudo-first-order kinetics suggesting the presence of multiple isomeric forms: Rh11+, Rh12+ and Rh8- exhibit characteristic biexponential decays which are interpreted in terms of the existence of different structural forms of the cluster which have markedly different reactivity. By contrast, Rh6+, Rh7+ and Rh8+ show rates which apparently increase with time, probably due to collisional activation. Thermalisation of the clusters prior to reaction by exposure to pulses of argon results in changes to the kinetics of these anomalous systems which can be explained in terms of collision induced isomerisation.     

AlmN2 (m=1～8)团簇结构与稳定性的DFT研究

用密度泛函理论(DFT)的B3LYP方法，在6-31G*水平上对AlmN2(m=1～8) 团簇的几何构型、电子结构、振动频率和热力学性质进行了理论研究. 结果表明， AlmN2团簇的基态结构有两种基本构型， m≤2的结构是以N－N键为核心周围与Al原子相配位形成的; m >2的结构是由两个AlnN（n< m）分子碎片通过共用Al原子或Al－Al键相互结合形成的，这为较快找到AlmN2团簇基态结构提供了一条有效途径. 通过对基态结构能量二次差分的讨论，得到了m为偶数的AlmN2团簇比m为奇数的稳定.     

Nitric oxide decomposition on small rhodium clusters, Rh(n)+/-

The decomposition of nitric oxide on small charged rhodium clusters Rh(n)(+/-) (6 < n < 30) has been investigated by Fourier transform ion cyclotron resonance mass spectrometry. For both cationic and anionic naked clusters, the rates of reaction with NO increase smoothly with cluster size in the range studied without the dramatic size-dependent fluctuations often associated with the reactions of transition-metal clusters. The cationic clusters react significantly faster than the anions and both exhibit rate constants exceeding collision rates calculated by average dipole orientation theory. Both the approximate magnitude and the trends in reactivity are modeled well by the surface charge capture model recently proposed by Kummerl?we and Beyer. All clusters studied here exhibit pseudo-first-order kinetics with no sign of biexponential kinetics often interpreted as evidence for multiple isomeric structures. Experiments involving prolonged exposure to NO have revealed interesting size-dependent trends in the mechanism and efficiency of NO decomposition: For most small clusters (n < 17), once two NO molecules are coadsorbed on a cluster, N(2) is evolved, generating the corresponding dioxide cluster. By analogy with experiments on extended surfaces, this observation is interpreted in terms of the dissociative adsorption of NO in the early stages of reaction, generating N atoms that are mobile on the surface of the cluster. For clusters where n < 13, this chemistry, which occurs independently of the cluster charge, repeats until a size-dependent, limiting oxygen coverage is achieved. Following this, NO is observed to adsorb on the oxide cluster without further N(2) evolution. For n = 14-16 no single end-point is observed and reaction products are based on a small range of oxide structures. By contrast, no evidence for N(2) production is observed for clusters n = 13 and n > 16, for which simple sequential NO adsorption dominates the chemistry. Interestingly, there is no evidence for the production of N(2)O or NO(2) on any of the clusters studied. A simple general mechanism is proposed that accounts for all observations. The detailed decomposition mechanisms for each cluster exhibit size (and, by implication, structure) dependent features with Rh(13)(+/-) particularly anomalous by comparison with neighboring clusters.