Density-functional calculations for small iron clusters: Fen, Fen+, and Fen- for n <= 5

Cryo spectroscopy of N2 on cationic iron clusters

Infrared photodissociation (IR-PD) spectra of iron cluster dinitrogen adsorbate complexes [Fe_n(N₂)_m]⁺ for n = 8-20 reveal slightly redshifted IR active bands in the region of 2200-2340 cm^-1. These bands mostly relate to stretching vibrations of end-on coordinated N₂ chromophores, a μ_1,end end-on binding motif. Density Functional Theory (DFT) modeling and detailed analysis of n = 13 complexes are consistent with an icosahedral Fe₁₃ ⁺ core structure. The first adsorbate shell closure at (n,m) = (13,12)-as recognized by the accompanying paper on the kinetics of N₂ uptake by cationic iron clusters-comes with extensive IR-PD band broadening resulting from enhanced couplings among adjacent N₂ adsorbates. DFT modeling predicts spin quenching by N₂ adsorption as evidenced by the shift of the computed spin minima among possible spin states (spin valleys). The IR-PD spectrum of (17,1) surprisingly reveals an absence of any structure but efficient non-resonant fragmentation, which might indicate some weakly bound (roaming) N₂ adsorbate. The multiple and broad bands of (17,m) for all other cases than (17,1) and (17,7) indicate a high degree of variation in N₂ binding motifs and couplings. In contrast, the (17,7) spectrum of six sharp bands suggests pairwise equivalent N₂ adsorbates. The IR-PD spectra of (18,m) reveal additional features in the 2120-2200 cm^-1 region, which we associate with a μ_1,side side-on motif. Some additional features in the (18,m) spectra at high N₂ loads indicate a μ_1,tilt tilted end-on adsorption motif.

Kinetics of stepwise nitrogen adsorption by size-selected iron cluster cations: Evidence for size-dependent nitrogen phobia

We present a study of stepwise cryogenic N₂ adsorption on size-selected Fe_n ⁺ (n = 8-20) clusters within a hexapole collision cell held at T = 21-28 K. The stoichiometries of the observed adsorption limits and the kinetic fits of stepwise N₂ uptake reveal cluster size-dependent variations that characterize four structural regions. Exploratory density functional theory studies support tentative structural assignment in terms of icosahedral, hexagonal antiprismatic, and closely packed structural motifs. There are three particularly noteworthy cases, Fe₁₃ ⁺ with a peculiar metastable adsorption limit, Fe₁₇ ⁺ with unprecedented nitrogen phobia (inefficient N₂ adsorption), and Fe₁₈ ⁺ with an isomeric mixture that undergoes relaxation upon considerable N₂ uptake.

Evolution of Ferromagnetic and Antiferromagnetic States in Iron Nitride Clusters FenN and FenN2 (n = 1-10)

First-principles density functional theory calculations on neutral and singly negatively and positively charged iron clusters Fe_n and iron nitride clusters Fe_nN and Fe_nN₂ (n = 1-10) in the range of 1 ≤ n ≤ 10 revealed that there is a strong competition between ferromagnetic and antiferromagnetic states especially in the Fe_nN₂^0,±1 cluster series. This phenomenon was related to superexchange via a bridging N atom between two iron atoms in the Fe_nN₂^0,±1 cluster series and to a double superexchange effect via a Fe atom shared by two N atoms in the Fe_nN₂^0,±1 series. A thorough examination of the structure-energy-spin state relationships in these clusters is conducted, leading to new insights and confirmation of available experimental results on structural parameters and dissociation energetics. The bond energies of both nitrogen atoms in the Fe_nN₂ series are approximately the same. They weakly depend on the charge of the host cluster and fluctuate around 5.5 eV when moving along the series. The energy of N₂ desorption is relatively small; it varies by about 1.0 eV and depends on the charge of the cluster. The experimental finding that N₂ dissociates on the Fe_n⁺ clusters beginning with n = 4 was supported by the results of our computations. Our computed values of the Fe_n⁺-N bonding energies agree with the experimental data within the experimental uncertainty bars. It was found that the attachment of one or two N atoms does not seriously affect the polarizability, electron affinity, or ionization energy of the host iron clusters independent of the charge.

Dissociation of dinitrogen on iron clusters: a detailed study of the Fe16 + N2 case

The coalescence of two Fe8N as well as the structure of the Fe16N2 cluster were studied using density functional theory with the generalized gradient approximation and a basis set of triple-zeta quality. It was found that the coalescence may proceed without an energy barrier and that the geometrical structures of the resulting clusters depend strongly on the mutual orientations of the initial moieties. The dissociation of N2 is energetically favorable on Fe16, and the nitrogen atoms share the same Fe atom in the lowest energy state of the Fe16N2 species. The attachment of two nitrogen atoms leads to a decrease in the total spin magnetic moment of the ground-state Fe16 host by 6 μB due to the peculiarities of chemical bonding in the magnetic clusters. In order to gain insight into the dependence of properties on charge and to estimate the bonding energies of both N atoms, we performed optimizations of Fe16N and the singly charged ions of both Fe16N2 and Fe16N. It was found that the electronic properties of the Fe16N2 cluster, such as electron affinity and ionization energy, do not appreciably depend on the attachment of nitrogen atoms but that the average binding energy per atom changes significantly. The lowering in total energy due to the attachment of two N atoms was found to be nearly independent of charge. The IR and Raman spectra were simulated for Fe16N2 and its ions, and it was found that the positions of the most intense peaks in the IR spectra strongly depend on charge and therefore present fingerprints of the charged states. The chemical bonding in the ground-state Fe16N20,±1 species was described in terms of the localized molecular orbitals.

Activation of Methane by the Iron Dimer Cation. A Theoretical Study

A detailed investigation of the reaction mechanisms underlying the observed reactivity of the iron dimer cation with respect to methane has been performed by density functional hybrid (B3LYP) and nonhybrid (BPW91) calculations. Minima and transition states have been fully optimized and characterized along the potential energy surfaces leading to three different exit channels for both the ground and the first excited states of the dimer. A comparison with our previous work covering the reactivity of the Fe(+) monomer was made to underline similarities and differences of the reaction mechanisms. Results show that geometric arrangements corresponding to bridged positions of the ligands with respect to iron atoms are always favored and stabilize intermediates, transition states and products, facilitating their formation. Binding energies of reaction products have been computed and compared with experimental measurements, and ELF analysis of the bond has been performed to rationalize trends as a function of the structure.

DFT study of the monomers and dimers of 2-pyrrolidone: equilibrium structures, vibrational, orbital, topological, and NBO analysis of hydrogen-bonded interactions

A computational study of the monomers and hydrogen-bonded dimers of 2-pyrrolidone was executed at different DFT levels and basis sets. The above dimeric complexes were treated theoretically to elucidate the nature of the intermolecular hydrogen bonds, geometry, thermodynamic parameters, interaction energies, and charge transfer. The processes of dimer formation from monomers and concerted reactions of double proton transfer were considered. The evolution of geometry, vibrational frequencies, charge distribution, and AIM properties in going from monomers to dimers was systematically followed. The solvent effects upon dimer formation were investigated in terms of the self-consistent reaction field (SCRF Onsager model). For the monomers and three dimers, vibrational frequencies were calculated and the changes in frequencies of the vibrations most sensitive to complexation were discussed. The orbital interactions were shown to lengthen the X-H (X = N, O) bond and lower its vibrational frequency (a red shift). To better understand the nature of the corresponding intermolecular interactions, we performed natural bond orbital (NBO) analysis. Topological analysis of electron density at bond critical points (BCP) was executed for complex molecules using the Bader's atoms in molecules (AIM) theory. The interaction energies were calculated, and the basis set superposition errors (BSSE) were estimated systematically. Satisfactory correlations between the structural parameters, interaction energies, and electron density characteristics at BCP were found.

Homonuclear transition-metal trimers

Density-functional theory has been used to determine the ground-state geometries and electronic states for homonuclear transition-metal trimers constrained to equilateral triangle geometries. This represents the first application of consistent theoretical methods to all of the ten 3d block transition-metal trimers, from scandium to zinc. A search of the potential surfaces yields the following electronic ground states and bond lengths: Sc3(2A1',2.83 A), Ti3(7E',2.32 A), V3(2E",2.06 A), Cr3(17E',2.92 A), Mn3(16A2',2.73 A), Fe3(11E",2.24 A), Co3(6E",2.18 A), Ni3(3A2",2.23 A), Cu3(2E',2.37 A), and Zn3(1A1',2.93 A). Vibrational frequencies, several low-lying electronic states, and trends in bond lengths and atomization energies are discussed. The predicted dissociation energies DeltaE(M3-->M2+M) are 49.4 kcal mol(-1)(Sc3), 64.3 kcal mol(-1)(Ti3), 60.7 kcal mol(-1)(V3), 11.5 kcal mol(-1)(Cr3), 32.4 kcal mol(-1)(Mn3), 61.5 kcal mol(-1)(Fe3), 78.0 kcal mol(-1)(Co3), 86.1 kcal mol(-1)(Ni3), 26.8 kcal mol(-1)(Cu3), and 4.5 kcal mol(-1)(Zn3).

Evolution of small copper clusters and dissociative chemisorption of hydrogen

The structural evolution of small copper clusters of up to 15 atoms and the dissociative chemisorption of H2 on the minimum energy clusters are studied systematically using density functional theory. The preferred copper sites for chemisorption are identified and the transition state structures and activation barriers for clusters four to nine atoms are determined and found to be inconsistent with the empirical Bronsted-Evans-Polanyi relationship. The physicochemical properties of the clusters are computed and compared with the bulk and surface values. The results indicate that a phase transition must occur in the going from cluster to bulk.

Reactions of cationic iron clusters with ammonia, models of nitrogen hydrogenation and dehydrogenation

The gas phase reactivities of small cationic iron clusters, Fen+ (n = 1-20), towards ammonia were investigated using Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry. Sequential addition of ammonia molecules to the clusters was observed to be the dominating process for n > 4. In the case of n = 4 we observed addition of ammonia accompanied by dehydrogenation. This reaction was modelled using hybrid density functional theory. Clusters with n < 4 do not react with ammonia. Clusters Fen+ (n = 1-20) react with neither N2 nor H2 at around 10(-8) mbar. When dinitrogen was seeded into the expanding helium, mixed clusters of the type FenNm+ were observed. These ions react with H2, either by addition, or by substitution of N2. The clusters with m = 1 were isolated in separate experiments and reacted with H2, which showed that mixed clusters with n = 5-13 add up to 5 molecules of dihydrogen in successive slow reactions.