Properties of Small Bimetallic Ni−Cu Clusters

Exploring structures of small anionic nickel-ethanol clusters with infrared spectroscopy

Small anionic nickel clusters with ethanol are investigated with a combination of mass-selective infrared photodissociation spectroscopy in a molecular beam and density functional theory simulations at the BLYP/6-311g(d,p) and TPSSh/def2-TZVPP level. In this context, the O-H stretching vibration of the ethanol is analyzed to obtain information about the structural motif, the geometry of the metal core, and the spin state of the clusters. For the [Ni2(EtOH)]- and [Ni3(EtOH)]- clusters, we assign quartet states of motifs with a hydrogen bond from the ethanol to the linear nickel core. The aggregation of a further ethanol molecule, yielding the [Ni3(EtOH)2]- cluster, results in the formation of a cooperative hydrogen bond network between the nickel core and the two ethanol molecules.

A Perspective on Modelling Metallic Magnetic Nanoparticles in Biomedicine: From Monometals to Nanoalloys and Ligand-Protected Particles

The focus of this review is on the physical and magnetic properties that are related to the efficiency of monometallic magnetic nanoparticles used in biomedical applications, such as magnetic resonance imaging (MRI) or magnetic nanoparticle hyperthermia, and how to model these by theoretical methods, where the discussion is based on the example of cobalt nanoparticles. Different simulation systems (cluster, extended slab, and nanoparticle models) are critically appraised for their efficacy in the determination of reactivity, magnetic behaviour, and ligand-induced modifications of relevant properties. Simulations of the effects of nanoscale alloying with other metallic phases are also briefly reviewed.

Dressing of Cu Atom over Nickel Cluster Stimulating the Poisoning-Free CO Oxidation: An Ab Initio Study

In this work using first-principles calculations based on spin-polarized density functional theory (DFT), the role of the Cu atom in degrading the poisoning of carbon monoxide (CO) over Ni_nCu clusters is unveiled. The search has been initiated with the examination of structural, magnetic, and electronic properties of Ni_n+1 and Ni_nCu clusters (1 ≤ n ≤ 12). X-ray absorption near-edge structure (XANES) spectra of Ni K-edge are computed to extract the information on the oxidation states and coordination environment of metal sites of the clusters. This study is operated with the two forms of dispersion corrections, i.e., D2 and D3, with standard DFT (with LDA and GGA functionals) for the consideration of van der Waals interactions during CO adsorption. The PBE and PBE-D3 approaches are found to be capable of yielding the experimentally observed preferential site for CO adsorption. The effect of spin-polarization on the reactivity of transition metals (TMs) toward CO adsorption is crucially assessed by the electronic reactivity descriptors such as d-band center, d-band width, and fractional filling of d-band using a spin-polarized d-band center model. The effective charge transfer from Cu to Ni atoms makes Ni atoms more efficient of charge and is attributed to the degrading adsorption of CO over Ni_nCu clusters. The Ni₁₂Cu cluster stands out with good CO oxidation activity for the Langmuir-Hinshelwood (L-H) reaction pathway.

DFT insights into structural effects of Ni-Cu/CeO2 catalysts for CO selective reaction towards water-gas shift

The water-gas shift (WGS) reaction is a key step in hydrogen production, particularly to meet the high-purity H2 requirement of PEM fuel cells. The catalysts currently employed in large-scale WGS plants require a two-step process to overcome thermodynamic and kinetic limitations. Ni-Cu/CeO2 solids are promising catalysts for the one-step process required for small-scale applications, as the addition of Cu hinders undesired methanation reactions occurring on Ni/CeO2. In this work, we performed calculations on Ni4-xCux/CeO2(111) systems to evaluate the influence of cluster conformation on the selectivity towards water-gas shift. The structure and miscibility of CeO2-supported Ni4-xCux clusters were investigated and compared with those of gas-phase clusters to understand the effect of metal-support interactions. The adsorption of CO onto apical Ni and Cu atoms of Ni4-xCux/CeO2(111) systems was studied, and changes in the C-O bond strength were confirmed at the electronic level by investigating shifts in the 3σ and 1π orbitals. The selectivity towards WGS was evaluated using Brønsted-Evans-Polanyi relations for the C-O activation energy. Overall, a strengthening of the C-O bond and an increase in CO dissociation energy were verified on Cu-containing clusters, explaining the improvement in selectivity of Ni4-xCux/CeO2(111) systems.

Quantum-Chemical Design of Molecular Structures of Tetra-, Penta- and Hexanuclear Metal Clusters Containing Aluminum and 3d-Element Atoms

Various data on the structural and thermodynamic characteristics of polynuclear metal clusters containing atoms of aluminum and various d-elements with the general formula Al_nM_m where (n + m) is 4, 5, or 6, and which can be precursors for the formation of nanoparticles of elemental metals or intermetallic compounds, have been systematized and discussed. It has been noted that each of these metal clusters in principle is able to exist in very diverse structural isomers, differing significantly among themselves in terms of the total energy and spin multiplicity of the ground state, the number of which is determined by both the specific values of n and m, and the nature of d-elements in their compositions. The presence of very complex dynamics with respect to the changes of the individual thermodynamic characteristics of the metal clusters under consideration as well as the thermodynamic parameters of the reactions of their formation, depending on the nature of the d-element, were also ascertained. In the main, the given review is devoted to the authors’ works published over the last 10 years. Bibliography – 96 references.

Uniaxial magnetic anisotropy energy of bimetallic Co-Ni clusters from a first-principles perspective

Along with the growing precision in the control of matter at increasingly smaller size scales, a field of research, based onto magnetic materials of technical interest, such as bimetallic clusters, has been developed in very recent years. Thereby, here, we report on a complete study of bimetallic clusters composed of cobalt and nickel with up to 7 atoms using ab initio methods in the GGA approach. We applied an unrestricted search method based on the tensor of inertia eigenvalues to find the most stable configurations of the clusters, obtaining a diverse set of structures with different geometric properties. We explored the effect of composition on the structural properties, the chemical stability, the magnetization and the magnetic anisotropy energy (MAE) of the so-obtained systems. Our results indicate that the behavior of the clusters is mainly governed by the Co-Co interaction and to a lesser extent by the Co-Ni and Ni-Ni interactions. Furthermore, for a given cluster size the magnetic moment increases by 2 μB/Co-substitution plus 1 μB/Ni-substitution coming from the cobalt and nickel core d-states, while in some cases unpaired hybrid s-electrons can also give rise to itinerant magnetism. These features have been analyzed with the help of a Jellium model and have important consequences for the magnetism and the magnetic anisotropy of the clusters. The magnetic behavior and MAE present complex and intriguing landscapes, which suggests the possibility of finely controlling the magnetic states, by tuning the cluster composition, aiming at technical implementation in fields such as molecular magnetism or quantum computation. In particular, cases such as Co6Ni, Co4Ni3, CoNi6 Co2Ni5 and Co3Ni3 present high relative stability and enhanced magnetic moments (around 10 μB), what makes them promising candidates for applications such as subnanometer magnetic information storage.

Insights into geometries, stabilities, electronic structures, reactivity descriptors, and magnetic properties of bimetallic Nim Cun-m (m = 1, 2; n = 3-13) clusters: Comparison with pure copper clusters

A long-range corrected density functional theory (LC-DFT) was applied to study the geometric structures, relative stabilities, electronic structures, reactivity descriptors and magnetic properties of the bimetallic NiCu_n-1 and Ni₂ Cu_n-2 (n = 3-13) clusters, obtained by doping one or two Ni atoms to the lowest energy structures of Cu_n , followed by geometry optimizations. The optimized geometries revealed that the lowest energy structures of the NiCu_n-1 and Ni₂ Cu_n-2 clusters favor the Ni atom(s) situated at the most highly coordinated position of the host copper clusters. The averaged binding energy, the fragmentation energies and the second-order energy differences signified that the Ni doped clusters can continue to gain an energy during the growth process. The electronic structures revealed that the highest occupied molecular orbital and the lowest unoccupied molecular orbital energies of the LC-DFT are reliable and can be used to predict the vertical ionization potential and the vertical electron affinity of the systems. The reactivity descriptors such as the chemical potential, chemical hardness and electrophilic power, and the reactivity principle such as the minimum polarizability principle are operative for characterizing and rationalizing the electronic structures of these clusters. Moreover, doping of Ni atoms into the copper clusters carry most of the total spin magnetic moment. © 2018 Wiley Periodicals, Inc.

Low-Lying Electronic States of CuAu

Coinage metal diatomic molecules are building blocks for nanostructured materials, electronic devices, and catalytically or photochemically active systems that are currently receiving lively interest in both fundamental and applied research. The theoretical study presented here elucidates the electronic structure in the ground and several low-lying excited states of the diatomic molecule CuAu that result from the combination of the atoms in their ground states nd(10)(n + 1)s(1 2)S and lowest excited d-hole states nd(9)(n + 1)s(2 2)D (n = 3 for Cu, n = 5 for Au). Full and smooth potential energy curves, obtained at the multireference configuration interaction (MRCI) level of theory, are presented for the complete set of the thus resulting 44 Λ-S terms and 86 Ω terms. Our approach is based on a scalar relativistic description using the Douglas-Kroll-Hess (DKH) Hamiltonian, with subsequent perturbative inclusion of spin-orbit (SO) coupling via the spin-orbit terms of the Breit-Pauli (BP) Hamiltonian. The Ω terms span an energy interval of about 7 eV at the ground state's equilibrium distance. Spectroscopic constants, calculated for all terms, are shown to accurately reproduce the observation for those nine terms that are experimentally known.

A DFT study of the chemical and optical properties of 7-atom Ag–X [X = Li, Na] nanoalloys for potential applications in opto-electronics and catalysis

In this paper, Ag atoms are substituted by X (Li, Na) atoms to form Ag_mX_(7−m) clusters to explore their electronic, chemical and optical properties in the framework of density functional theory (DFT).

Potential energy curves for the ground and low-lying excited states of CuAg

The ground and low-lying excited states of heteronuclear diatomic CuAg are examined by multi-reference configuration interaction (MRCI) method. Relativistic effects were treated and probed in two steps. Scalar terms were considered using the spin-free DKH Hamiltonian as a priori and spin-orbit coupling was calculated perturbatively via the spin-orbit terms of the Breit-Pauli Hamiltonian based on MRCI wavefunctions. Potential energy curves of the spin-free states and their corresponding Ω components correlating with the separated atom limits (2)S(Cu) + (2)S(Ag) and (2)D(Cu) + (2)S(Ag) are obtained. The results are in fine agreement with the experimental measurements and tentative conclusions for the ion-pair B0(+) state are confirmed by our theoretical calculations. Illustrative results are presented to reveal the relative importance and magnitude of the scalar and spin-orbit effects on the spectroscopic properties of this molecule. Time dependent density functional theory calculations, using the LDA, BLYP, B3LYP, and SAOP functionals have been carried out for CuAg and the accuracy of TD-DFT has been compared with ab initio results.

An optimized interatomic potential for Cu-Ni alloys with the embedded-atom method

We have developed a semi-empirical and many-body type model potential using a modified charge density profile for Cu-Ni alloys based on the embedded-atom method (EAM) formalism with an improved optimization technique. The potential is determined by fitting to experimental and first-principles data for Cu, Ni and Cu-Ni binary compounds, such as lattice constants, cohesive energies, bulk modulus, elastic constants, diatomic bond lengths and bond energies. The generated potentials were tested by computing a variety of properties of pure elements and the alloy of Cu, Ni: the melting points, alloy mixing enthalpy, lattice specific heat, equilibrium lattice structures, vacancy formation and interstitial formation energies, and various diffusion barriers on the (100) and (111) surfaces of Cu and Ni.

Effects of trimethylaluminium and tetrakis(ethylmethylamino) hafnium in the early stages of the atomic-layer-deposition of aluminum oxide and hafnium oxide on hydroxylated GaN nanoclusters

We calculate the interactions of two atomic layer deposition (ALD) reactants, trimethylaluminium (TMA) and tetrakis(ethylmethylamino) hafnium (TEMAH) with the hydroxylated Ga-face of GaN clusters when aluminum oxide and hafnium oxide, respectively, are being deposited. The GaN clusters are suitable as testbeds for the actual Ga-face on practical GaN nanocrystals of importance not only in electronics but for several other applications in nanotechnology. We find that TMA spontaneously interacts with hydroxylated GaN; however it does not follow the atomic layer deposition reaction path unless there is an excess in potential energy introduced in the clusters at the beginning of the optimization, for instance, using larger bond lengths of various bonds in the initial structures. TEMAH also does not interact with hydroxylated GaN, unless there is an excess in potential energy. The formation of a Ga-N(CH3)(CH2CH3) bond during the ALD of HfO2 using TEMAH as the reactant without breaking the Hf-N bond could be the key part of the mechanism behind the formation of an interface layer at the HfO2/GaN interface.

Density functional study of structural and electronic properties of small binary Be(n)Cu(m) (n + m = 2~7) clusters

The geometrical structures, electronic properties and relative stabilities of small bimetallic Be n Cu m (n + m = 2-7) clusters have been systematically investigated by using a density functional method at the B3PW91 level. In the most stable structures of Be n Cu m , the Be atoms tend to gather together and construct similar configurations to those of pure Be n clusters. Meanwhile, there is a tendency for Cu atoms to segregate toward the Be n cluster surface. The successive binding energies, cohesive energies, second difference of energies, the highest occupied-lowest unoccupied molecular orbital energy gaps and chemical hardness of Be n Cu m are also investigated. All of them demonstrate that the clusters with even number of copper atoms present relatively higher stabilities. The natural population analyses on the Be n Cu m clusters reveal that, the charge transfers from Be to Cu when the average coordination numbers (Nc) of Be atom is less than 3, whereas the charge-transferring direction reverses when Nc(Be) increases.

Molecular electrostatic potentials of DNA base-base pairing and mispairing

An understanding of why adenine (A) pairs with thymine (T) and cytosine (C) with guanine (G) in DNA is very useful in the design of sensors and other related devices. We report the use of dissociation energies, geometries and molecular electrostatic potentials (MEPs) to justify the canonical (AT and CG) Watson-Crick pairs. We also analyze all mismatches in both configurations-cis and trans-with respect to their glycoside bonds. As expected, we found that the most stable pair configuration corresponds to CG, providing an energy criterion for that preferred configuration. The reason why A gets together with T is much more difficult to explain as the energy of this pair is smaller than the energy of some other mismatched pairs. We tested MEPs to see if they could shed light on this problem. Interestingly, MEPs yield a unique pattern (shape) for the two canonical cases but different shapes for the mismatches. A tunnel of positive potential surrounded by a negative one is found interconnecting the three H-bonds of CG and the two of AT. This MEP tunnel, assisted partially by energetics and geometrical criteria, unambiguously determine a distinctive feature of the affinity between A and T as well as that between G and C.

Light activation of the isomerization and deprotonation of the protonated Schiff base retinal

We perform an ab initio analysis of the photoisomerization of the protonated Schiff base of retinal (PSB-retinal) from 11-cis to 11-trans rotating the C10-C11=C12-C13 dihedral angle from 0° (cis) to -180° (trans). We find that the retinal molecule shows the lowest rotational barrier (0.22 eV) when its charge state is zero as compared to the barrier for the protonated molecule which is ∼0.89 eV. We conclude that rotation most likely takes place in the excited state of the deprotonated retinal. The addition of a proton creates a much larger barrier implying a switching behavior of retinal that might be useful for several applications in molecular electronics. All conformations of the retinal compound absorb in the green region with small shifts following the dihedral angle rotation; however, the Schiff base of retinal (SB-retinal) at trans-conformation absorbs in the violet region. The rotation of the dihedral angle around the C11=C12 π-bond affects the absorption energy of the retinal and the binding energy of the SB-retinal with the proton at the N-Schiff; the binding energy is slightly lower at the trans-SB-retinal than at other conformations of the retinal.

A Density Functional Theory Investigation of the Bimetallic Clusters Nb2Rh and NbRh2 and the Complexes They Form with CO

The interaction of CO with the bimetallic clusters Nb2Rh and NbRh2 has been theoretically investigated using density functional theory. The lowest energy structure of Nb2Rh is found to be a doublet Cs scalene triangle and the global minimum of Nb2Rh–CO is a dissociative structure with C1 symmetry. The lowest energy minimum of NbRh2 is found to be a doublet C2v isosceles triangle and the global minimum of NbRh2–CO is a dissociative structure with Cs symmetry. In comparison with our previous work on Rh3 + CO (J. Comp. Chem., 2008, 29, 1497), these results show that substitution of a single Rh atom with Nb is sufficient to dissociate CO.

Transverse electronic transport in double-stranded DNA nucleotides

We calculate the transverse current through double-stranded DNA nucleotides using ab initio techniques in order to establish a protocol to recognize the type and sequence of double-stranded DNA nucleotides. The distinct current-voltage features between nucleotides are used as signatures for their characterization and sequencing. Extended bulk gold electrodes as well as extensions of the DNA backbones are tested as contacts for the electron transport, yielding currents 2 orders of magnitude larger for the former. The addition of Na or H positive counterions improves the signal levels, thus leading to a better discrimination, especially when sodium cations are added.

Geometry and stability of BenCm (n=1-10; m=1, 2, ..., to 11-n) clusters

Density functional theory B3PW91/6-31+G* calculations on BenCm (n=1-10; m=1, 2, ..., to 11-n) clusters have been carried out to examine the effect of cluster size, relative composition, binding energy per atom, HOMO-LUMO gap, vertical ionization potential, and electron affinity on their relative stabilities. The most stable planar cyclic conformations of these clusters always show at least a set of two carbon atoms between two beryllium atoms, while structures where beryllium atoms cluster together, or allow the intercalation of one carbon atom between two of them, generally seem to be the least stable ones. Clusters containing 1, 2, and 3 beryllium atoms (Be2C8, Be3C6, Be2C6, BeC6, Be2C4, BeC4, Be2C2, and BeC2) are identified as clusters of "magic numbers" in terms of their high binding energy per atom, high HOMO-LUMO gap, vertical ionization potential, and second difference in energy per beryllium atom.

Growth format, electronic architecture, magnetic, and optical properties of aromaticcyclo-Cu3Au3homotops

Bimetallic Cu(3)Au(3) clusters have been investigated using electronic structure calculation techniques (DFT) to understand their electronic, magnetic, and optical properties as well as the geometrical structures. The most stable homotop is the planar cyclo-[Cu(3)(micro-Au)(3)] form consisting of a triangular positively charged Cu(3) structural core with negatively charged Au atoms occupying exposed positions. This structure is characterized by the maximum number of heterobonds and peripheral positions of Au atoms. Possible growth formats of the cyclo-[Cu(3)(micro-Au)(3)] homotops have been explored following both the edge-capping and the stepwise metal atom substitution mechanism. The bonding pattern along with the density of states (DOS) plots of the cyclo-[Cu(3)(micro-Au)(3)] homotop are thoroughly analyzed and compared with those of the pure cyclo-[Cu(3)(micro-Cu)(3)] and cyclo-[Au(3)(micro-Au)(3)] clusters. Particular attention was paid on the stability of these bimetallic clusters in relation with the ring-shaped electron density distribution (aromaticity). It was found that all 3-membered metal rings exhibit significant aromatic character, which was verified by a number of established criteria of aromaticity, such as structural, energetic, magnetic (NICS profiles), and out-of-plane ring deformability criteria. The NICS (1) values correlate well with the out-of-plane ring deformation energy. Finally, a comprehensive analysis of the optical spectra of the CuAu, Cu(2), and Au(2) diatomics and the cyclo-[Cu(3)(micro-Au)(3)], cyclo-[Cu(3)(micro-Cu)(3)], and cyclo-[Au(3)(micro-Au)(3)] clusters placed the electronic assignments of the optical transitions on a firm footing.

Density Functional Study of Structural and Electronic Properties of Small Bimetallic Silver−Nickel Clusters

Theoretical study on the structure and electronic properties of small AgmNip (m + p < or = 6) clusters has been carried out in the framework of density functional theory. Structural features, cohesive energies, vertical ionization potentials, and charge transfers are evaluated for each Ag/Ni ratio. In all the AgmNip clusters, the nickel atoms are brought together, yielding a maximum of Ni-Ni bonds, and the silver atoms are located around a Ni core with a maximum of Ag-Ni bonds. The ionization potential and the highest occupied molecular orbital shape are directly related to the two- or three-dimensional character of the cluster's geometry. A very low electronic charge transfer from Ni to Ag is found, and the magnetic moment is located on Ni atoms but with a low spin polarization on silver in the Ni-rich clusters.

Periodic table of 3d-metal dimers and their ions

The ground states of the mixed 3d-metal dimers TiV, TiCr, TiMn, TiFe, TiCo, TiNi, TiCu, TiZn, VCr, VMn, VFe, VCo, VNi, VCu, VZn, CrMn, CrFe, CrCo, CrNi, CrCu, CrZn, MnFe, MnCo, MnNi, MnCu, MnZn, FeCo, FeNi, FeCu, FeZn, CoNi, CoCu, CoZn, NiCu, NiZn, and CuZn along with their singly negatively and positively charged ions are assigned based on the results of computations using density functional theory with generalized gradient approximation for the exchange-correlation functional. Except for TiCo and CrMn, our assignment agrees with experiment. Computed spectroscopic constants (r(e),omega(e),D(o)) are in fair agreement with experiment. The ground-state spin multiplicities of all the ions are found to differ from the spin multiplicities of the corresponding neutral parents by +/-1. Except for TiV, MnFe, and MnCu, the number of unpaired electrons, N, in a neutral ground-state dimer is either N(1)+N(2) or mid R:N(1)-N(2)mid R:, where N(1) and N(2) are the numbers of unpaired 3d electrons in the 3d(n)4s(1) occupation of the constituent atoms. Combining the present and previous results obtained at the same level of theory for homonuclear 3d-metal and ScX (X=Ti-Zn) dimers allows one to construct "periodic" tables of all 3d-metal dimers along with their singly charged ions.

Effect of Ni and Pd on the Geometry, Electronic Properties, and Active Sites of Copper Clusters

The geometry, electronic properties, and active sites of copper clusters doped with Ni or Pd atoms, Cu(n)()(-)(1)M (n = 2-6; M = Ni, Pd) have been investigated using first-principles methods. Planar structures are energetically favorable in Cu(n)()(-)(1)Ni (n = 2-6). However, for Pd-doped clusters, three-dimensional structures are competitive in energy, and for n = 6, the most stable structure is not planar. Several properties of doped copper clusters present odd-even oscillations as the number of copper atoms grow. The different atomic ground-state configuration of Ni and Pd determines the bonding and electronic properties of doped copper clusters. The interaction between impurities and copper atoms can modify the chemical hardness and active sites of doped copper clusters markedly inducing directionality in the reactivity. This effect is relevant to the behavior of catalysts as well as in the growth of metallic films.

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).

A density-functional study of Al-doped Ti clusters: TinAl (n=1–13)

Equilibrium geometries, stabilities, and electronic properties of TinAl (n = 1-13) clusters have been studied by using density-functional theory with local spin density approximation and generalized gradient approximation. The ground-state structures of TinAl clusters have been obtained. The resulting geometries show that the aluminum atom remains on the surface of clusters for n < 9, but is slowly getting trapped beyond n = 9, meanwhile, the Al atom exhibits a valent transition from monovalent to trivalent. The geometric effects and electronic effects clearly demonstrate the Ti4Al cluster to be endowed with special stability. The studies on the bonds indicate the change from ionic to metalliclike.