Density Functional Study of Neutral and Charged Silver Clusters Agn with n = 2–22. Evolution of Properties and Structure

Electronic properties and collision cross sections of AgOkHm± (k, m = 1-4) aerosol ionic clusters

Experimental evidence shows that hydroxylated metal ions are often produced during cluster synthesis by atmospheric pressure spark ablation. In this work, we predict the ground state equilibrium structures of AgOkHm± clusters (k and m = 1-4), which are readily produced when spark ablating Ag, using the coupled cluster with singles and doubles (CCSD) method. The stabilization energy of these clusters is calculated with respect to the dissociation channel having the lowest energy, by accounting perturbative triples corrections to the CCSD method. The interatomic interactions in each of the systems have been investigated using the frontier molecular orbital (FMO), natural bond orbital (NBO) and quantum theory of atoms in molecules (QTAIM) methods. Many of the ground states of these ionic clusters are found to be stable, corroborating experimental observations. We find that clusters having singlet spin states are more stable in terms of dissociation than the clusters that have doublet or triplet spin states. Our calculations also indicate a strong affinity of the ionic and neutral Ag atom towards water and hydroxyl radicals or ions. Many 3-center, 4-electron (3c/4e) hyperbonds giving rise to more than one resonance structure are identified primarily for the anionic clusters. The QTAIM analysis shows that the O-H and O-Ag bonds in the clusters of both polarities are respectively covalent and ionic. The FMO analysis indicates that the anionic clusters are more reactive than the cationic ones. Using the cluster structures predicted by the CCSD method, we calculate the collision cross sections of the AgOkHm± family, with k and m ranging from 1 to 4, by the trajectory method. In turn, we predict the electrical mobilities of these clusters when suspended in helium at atmospheric pressure and compare them with experimental measurements.

Reactivity of cationic silver clusters with O2: a probe of interplay between clusters' geometric and electronic structures

We explored the size-dependent reactivity of Agn+ (n = 2-22) with O2 under mild conditions and found that only a few sizes of Agn+, with even values of n = 4, 6, 12, 16, 18, and 22, are reactive. Possible structures of Agn+ (n = 2-22) were determined using a genetic algorithm with incomplete local optimizations at the DFT level, and the calculated bonding strengths of O2 on these structures are consistent with experimental observations. Analyses revealed a close relationship between the reactivity of Agn+ with O2 and its HOMO-LUMO gap: cationic silver clusters with a small HOMO-LUMO gap are reactive, which can be rationalized by the covalent character of chemical bonds between Agn+ and O2 involving their frontier orbitals. The peculiar size-dependent HOMO-LUMO gaps and reactivity with O2 correlate with the subtle interplay between the electronic configurations and geometric structures of these silver cluster cations.

Effect of Ligand Attachment at Ag11 for CO Oxidation: A Computational Investigation

Heterogeneous CO oxidation is a demanding reaction at room temperature due to the high activation energy required to break the O=O bond. While several metal clusters are reported to oxidize CO successfully, they fall short of their selectivity for the reaction and recyclability. In this regard, there is a need for economic catalysts with high catalytic activity, low activation barrier, and reusability. In this study, we have investigated the catalytic activity of the neutral pristine and ligated Ag11 cluster toward CO oxidation. We investigated the attachment effect of three organic donor ligands: trimethylphosphine, triethylphosphine, and N-ethyl pyrrolidone to the Ag11 cluster. Our results show that including donor ligands on the Ag11 cluster surface can significantly reduce the barrier heights for CO oxidation. The minimum barrier heights with the system coordinated with triethylphosphine showed the lowest activation barrier of 1.06 kcal/mol compared to the high activation barrier of 14.77 kcal/mol recorded for the pristine cluster. Exploration of the reaction mechanism and charge analysis showed that the electron donor ligands activate O2 via charge donation, thereby reducing the barrier heights of CO oxidation.

Structural and optical properties of the Agn-tyrosine complexes (n = 3-12): a density functional theory study

We study the optical properties of Agn (n = 3-12) neutral clusters and their coordination with a tyrosine (Tyr) molecule. A global search strategy coupled with density functional theory (DFT) computations explored the potential energy surface. Adsorption energy calculations predicted that Tyr coordination stabilizes the metal clusters, favouring the Agn-Tyr complexes with an even number of silver atoms. For the Agn low-lying isomers, the general shape and the major transitions of the calculated time dependent-DFT (TD-DFT) absorption spectra align with those of previous reports measured in an argon environment. We use the analysis of non-covalent interactions to identify the specific interactions between each silver cluster and functional groups of Tyr. The TD-DFT absorption spectra for the Agn-Tyr complexes showed that Tyr significantly modifies the optical properties of the coordinated silver clusters and affects the smaller systems to a greater extent. The optical absorption results of the bare Agn clusters and the Agn-Tyr complexes are compared and discussed in detail.

Understanding of the Effect of the Adsorption of Atom and Cluster Silver on Chitosan: An In Silico Analysis

In this work, the structural, electronic, and optical stability properties of the chitosan monomer (M-Ch) and atomic silver complex are reported, as well as a unitary cell of a silver cluster in the gas phase and acetic acid. The generalized gradient approximation HSEh1PBE/def2-TZVPP50 results established the structures' anionic charge (Q = -1|e|) and the doublet state (M = 2). The high cohesive energy indicates structural stability, and the quantum-mechanical descriptors show a high polarity and low chemical reactivity. Also, the quantum-mechanical descriptors present a low work function that shows the structures are suitable for applications in light-emitting diodes. Finally, the electronic behavior observed by the |HOMO-LUMO| gap energy changes depending on the atomic silver incorporated into the complex.

Adsorption of air pollutants onto silver and gold atomic clusters: DFT and PNO-LCCSD-F12 calculations

We provide a comprehensive investigation of intermolecular interactions between atmospheric gaseous pollutants, including CH₄, CO, CO₂, NO, NO₂, SO₂, as well as H₂O and Ag_n (n = 1-22) or Au_n (n = 1-20) atomic clusters. The optimized geometries of all the systems investigated in our study were determined using density functional theory (DFT) with M06-2X functional and SDD basis set. The PNO-LCCSD-F12/SDD method was used for more accurate single-point energy calculations. Compared to their isolated states, the structures of the Ag_n and Au_n clusters undergo severe deformations upon adsorption of the gaseous species, which become more significant as the size of the clusters decreases. Considering that, in addition to adsorption energy, we have determined the interaction and deformation energy of all the systems. All our calculations consistently show that among the gaseous species examined, SO₂ and NO₂ exhibit a higher preference for adsorption on both types of clusters, with a slightly higher preference for the Ag clusters compared to the Au clusters, with the SO₂/Ag₁₆ system exhibiting the lowest adsorption energy. The type of intermolecular interactions was investigated through wave function analyses, including natural bond orbital (NBO) and quantum theory of atoms in molecules (QTAIM), showing that NO₂ and SO₂ are chemisorbed on the Ag_n and Au_n atomic clusters, whereas the other gas molecules exhibit a much weaker interaction with them. The reported data can be used as input parameters for molecular dynamics simulations to study the selectivity of atomic clusters towards specific gases under ambient conditions, as well as to design materials that take advantage of the studied intermolecular interactions.

The global minimum of Ag30: a prolate spheroidal structure predicted using a genetic algorithm with incomplete local optimizations at the DFT level

Genetic algorithms have been widely used to explore global minimum points of atomic clusters, and their incorporation with ab initio calculations (including density functional theory methods) as local optimization approaches increases their ability to accurately locate the global minimum points on complicated potential energy surfaces. However, the local optimizations using ab initio calculations significantly increase the computational cost relative to those based on empirical or semi-empirical calculations. Herein, we develop a genetic algorithm program with an incomplete local optimization strategy at the DFT level. Using several representative clusters as test examples, this program showed high efficiency in locating their global minimum points. The low-lying isomers of Ag₃₀ were explored using this program, and the determined global minimum is a prolate spheroidal structure. The elongated spheroidal shape causes degeneracy lifting of the free electron shells, and endows Ag₃₀ with a large HOMO-LUMO gap. The sharp increase of silver clusters' reactivity around the sizes with 30 valence electrons observed in our previous experiments could be correlated with this theoretical figure.

Properties of Naked Silver Clusters with Up to 100 Atoms as Found with Embedded-Atom and Density-Functional Calculations

The structural and energetic properties of small silver clusters Agn with n = 2-100 atoms are reported. For n = 2-100 the embedded atom model for the calculation of the total energy of a given structure in combination with the basin-hopping search strategy for an unbiased structure optimization has been used to identify the energies and structures of the three energetically lowest-lying isomers. These optimized structures for n = 2-11 were subsequently studied further through density-functional-theory calculations. These calculations provide additional information on the electronic properties of the clusters that is lacking in the embedded-atom calculations. Thereby, also quantities related to the catalytic performance of the clusters are studied. The calculated properties in comparison to other available theoretical and experimental data show a good agreement. Previously unidentified magic (i.e., particularly stable) clusters have been found for n>80. In order to obtain a more detailed understanding of the structural properties of the clusters, various descriptors are used. Thereby, the silver clusters are compared to other noble metals and show some similarities to both copper and nickel systems, and also growth patterns have been identified. All vibrational frequencies of all the clusters have been calculated for the first time, and here we focus on the highest and lowest frequencies. Structural effects on the calculated frequencies were considered.

Electronic Structure, Stability, and Electrical Mobility of Cationic Silver Oxide Atomic Clusters

Silver oxide cluster cations (Ag_nO_m⁺) can readily be produced by a number of methods including atmospheric-pressure spark ablation of pure silver electrodes when trace amounts of oxygen are present in the carrier gas. Here we determine the equilibrium geometries of Ag_nO_m⁺ clusters (n = 1-4; m = 1-5) using accurate coupled cluster with singles and doubles (CCSD) method, while the stabilization energies are calculated with additional perturbative triples correction (CCSD(T)). Although a number of stable states have been identified, our results show that the Ag_nO_m⁺ clusters with m = 1 are more stable than those with m ≥ 2 due to the absence of the terminally attached O₂ molecule, corroborating recent observations by mass spectrometry. Using the computed structures, we calculate the electrical mobilities of the Ag_nO_m⁺ clusters and label the values on a respective experimentally determined spectrum in an attempt to better interpret the occurrence of the peaks and troughs in the measurements.

Parametrization of the PM7 Semiempirical Quantum Mechanical Method for Silver Nanoclusters

Semiempirical quantum mechanical methods (SEQMs) are widely used in computational chemistry because of their low computational cost, but their accuracy depends on the quality of the parameters. The neglect of diatomic differential overlap method PM7 is among the few SEQMs that contain parameters for Ag, but the experimental reference data was insufficient to obtain reliable parameters in the original parametrization. In this work, we reparametrize the PM7 parameters for Ag to accurately reproduce the ground-state potential energy surfaces of Ag clusters. Since little experimental data is available, we use reference data obtained from the ab initio method CCSD(T). The resulting parameters significantly reduce the errors in binding energies, energies required to displace clusters along their normal modes, and relative energies of isomers compared to the default PM7 Ag parameters.

Interactions of Nitric Oxide Molecules with Pure and Oxidized Silver ClustersAgn{plus minus}/AgnO{plus minus} (n=11-13). A Computational Study

In this work we studied, within DFT, the interaction of NO with pure and oxidized Agn, both anionic and cationic, composed from11 to 13 Ag atoms. In that size interval, shell closing effects are not expected, and structural and electronic odd-even effects will determine the strength of interaction. We obtained that species Agn{plus minus} and AgnO{plus minus} with odd number of electrons (n=12) adsorb NO with higher energy than their neighbours. This result agrees with the facts observed in recent mass spectroscopy measurements, which were performed at finite temperature. The adsorption energy is about twice for oxidized clusters compared to pure ones, and higher for anions than for cations. The adsorption of another NO molecule on AgnNO{plus minus} forms Agn(NO)2{plus minus}, with the dimer (NO)2 in cis configuration, and binding the two N atoms with two neigbour Ag atoms. The n=12 show the higher adsorption energy again. In absence of reaction barriers, Agn(NO)2{plus minus} dissociate spontaneously into AgnO{plus minus} and N2O, except the n= 12 anion. The máximum high barrier along the dissociation path of Ag13(NO)2- is about 0.7 eV. Further analysis of PDOS for Ag11-13 (NO)x{plus minus} (x=0,1,2) molecules shows that bonding between NO and Agn mainly occurs in the range between -3.0 eV and 3.0 eV. The overlap between 4 d of Ag and 2 p of N and O is larger for Ag12(NO)2{plus minus} than for neighbour sizes. For n=12, the d bands are close to the (NO)2 2π orbital, leading to extra back-donation charge from the 4 d of Ag to the closer 2π orbital of (NO)2.

Coordination of Ethylamine on Small Silver Clusters: Structural and Topological (ELF, QTAIM) Analyses

Amine ligands are expected to drive the organization of metallic centers as well as the chemical reactivity of silver clusters early growing during the very first steps of the synthesis of silver nanoparticles via an organometallic route. Density functional theory (DFT) computational studies have been performed to characterize the structure, the atomic charge distribution, and the planar two-dimensional (2D)/three-dimensional (3D) relative stability of small-size silver clusters (Ag_n, 2 ≤ n ≤ 7), with or without an ethylamine (EA) ligand coordinated to the Ag clusters. The transition from 2D to 3D structures is shifted from n = 7 to 6 in the presence of one EA coordinating ligand, and it is explained from the analysis of the Ag-N and Ag-Ag bond energies. For fully EA saturated silver clusters (Ag_n-EA_n), the effect on the 2D/3D transition is even more pronounced with a shift between n = 4 and 5. Subsequent electron localization function (ELF) and quantum theory of atoms in molecules (QTAIM) topological analyses allow for the fine characterization of the dative Ag-N and metallic Ag-Ag bonds, both in nature and in strength. Electron transfer from ethylamine to the coordinated silver atoms induces an increase of the polarization of the metallic core.

SERS of Dopamine: Computational and experimental studies

Computational and experimental studies have been carried out on Dopamine. The calculated Raman spectra of Dopamine with and without Silver clusters (Ag_n (n = 1-4)) are compared with each other and it is shown that the intensity of the Raman activity increases with increasing number of silver atoms. The SERS effect shown by this system is further supported by calculating the Global electrophilicity index ω, the static mean polarizability α₀, and the anisotropy of the polarizabilities Δα. Stabilities of the complexes are analysed using the charge transfer, stabilization energies, and interaction energies. The reactive parameters for these complexes were further supported by looking at the molecular electrostatic potential (MESP) surfaces. SERS substrates were fabricated by sintering silver nanoparticle paste onto a fused silica substrate, using a femtosecond laser. Detection of Dopamine up to 1 μM is reported using the SERS substrates.

Low-energy structures and electronic properties of small titanium nitride nanoclusters

Some small nanoclusters of plasmonically superior titanium nitride are generated usingab initiomolecular dynamics simulation under the regime of density functional theory. The global minima structures of TiN lack symmetry in the local environment as compared to the bulk counterpart. Electronic properties, namely Bader charge, electron localization function, and density of states, help us to have a deeper understanding of these nanoclusters. In all the calculations, bulk TiN has been taken as a reference to compare the properties. It has been observed that with an increasing number of atoms in the nanoclusters, they tend to show properties similar to bulk TiN, indicating that they might as well be used as plasmonic materials.

Computational Appraisal of Silver Nanocluster Evolution on Epitaxial Graphene: Implications for CO Sensing

Early stages of silver nucleation on a two-dimensional (2D) substrate, here, monolayer epitaxial graphene (MEG) on SiC, play a critical role in the formation of application-specific Ag nanostructures. Therefore, it is of both fundamental and practical importance to investigate the growth steps when Ag adatoms start to form a new phase. In this work, we exploit density functional theory to study the kinetics of early-stage nuclei Ag n (n = 1-9) assembly of Ag nanoparticles on MEG. We find that the Ag1 monomer tends to occupy hollow site positions of MEG and interacts with the surface mainly through weak dispersion forces. The pseudoepitaxial growth regime is revealed to dominate the formation of the planar silver clusters. The adsorption and nucleation energies of Ag n clusters exhibit evident odd-even oscillations with cluster size, pointing out the preferable adsorption and nucleation of odd-numbered clusters on MEG. The character of the interaction between a chemisorbed Ag3 cluster and MEG makes it possible to consider this trimer as the most stable nucleus for the subsequent growth of Ag nanoparticles. We reveal the general correlation between Ag/MEG interaction and Ag-Ag interaction: with increasing cluster size, the interaction between Ag adatoms increases, while the Ag/MEG interaction decreases. The general trend is also supported by the results of charge population analysis, according to which the average charge per Ag adatom in a Ag n cluster demonstrates a drastic decrement with cluster size increase. 2D-3D structural transition in Ag n clusters was investigated. We anticipate that the present investigation is beneficial by providing a better understanding of the early-stage nucleation of Ag nanoparticles on MEG at the atomic scale. Specific interaction between odd-numbered Ag clusters preadsorbed onto the MEG surface and carbon monoxide (CO) as well as clusters' stability at 300 K is discussed in terms of sensing applications.

Reply to 'Comment on "Structural characterization, reactivity, and vibrational properties of silver clusters: A new global minimum for Ag16"' by P. V. Nhat, N. T. Si, L. V. Duong and M. T. Nguyen, Phys. Chem. Chem. Phys., 2021, , DOI: D1CP00646K

Recently, P. V. Nhat et al., have discussed and commented on our article (DOI: 10.1039/D0CP04018E) for the case of the most stable structure of Ag15. They have found a new most stable structure (labeled as 15-1) in comparison to the putative global minimum reported by us, which is a four layered 1-4-6-4 stacking structure with a C2v point group (15-2). In this reply, we have performed a larger structure search which allowed us to confirm the results of Nhat et al. The results show the existence of multiple isoenergetic isomers with similar structure motifs for the Ag15 system, increasing the problem complexity to locate the global minimum. The results in regard to the structure and electronic properties of the new lowest energy structure are discussed.

Comment on 'Structural characterization, reactivity and vibrational properties of silver clusters: a new global minimum for Ag16' by P. L. Rodríguez-Kessler, A. R. Rodríguez-Domínguez, D. MacLeod Carey and A. Muñoz-Castro, Phys. Chem. Chem. Phys., 2020, 22, 27255, DOI: D0CP04018E

A recent paper by Rodríguez-Kessler et al., Phys. Chem. Chem. Phys., 2020, 22, 27255-27262, reported not only results of quantum chemical computations (using the PW91 density functional) on Ag₁₆ clusters as emphasized in the article's title, but also on the Ag₁₅ size. These authors confirmed previous results obtained by McKee and Samokhvalov (J. Phys. Chem. A, 2017, 121, 5018-5028 using the M06 density functional) that the most stable isomer of Ag₁₅ is a C_2v structure. We wish to point out that two low symmetry isomers of Ag₁₅ that have a similar energy content are even lower in energy than their reported C_2v global minimum. The relative energies between low-lying Ag₁₅ isomers were again found to be method-dependent, and within the expected accuracy of DFT and CCSD(T) methods they could be considered as energetically degenerate, and likely coexist in a molecular beam. The new lower-energy Ag₁₅ isomers appear to fit more consistently within the structural evolution of small silver clusters.

Theoretical Study of the Binding of the Thiol-Containing Cysteine Amino Acid to the Silver Surface Using a Cluster Model

Computational approaches within the framework of density functional theory (DFT) were employed to elucidate the binding mechanism of the cysteine amino acid on silver nanoparticles using several small silver clusters Ag_n with n = 2-10 as surface models. The long-range corrected LC-BLYP functional and correlation consistent basis sets cc-pVTZ-PP and cc-pVTZ were used to determine the structural features, energetics, and spectroscopic and electronic properties of the resulting complexes. In vacuum and highly acidic conditions, cysteine molecules prefer to adsorb on silver clusters via their amine group. In aqueous solution, the thiolate head turns out to be the most energetically favorable binding site. The cysteine affinity of silver clusters is greatly altered in different conditions, i.e., acidic solution < vacuum < aqueous solution, and is strongly dependent on the cluster size. As compared to free clusters, the frontier orbital energy gap of the ones capped by cysteine is significantly improved, which corresponds to stronger stability, especially in aqueous solution. The analysis of frontier orbitals also reveals that both forward and backward electron donations exhibit comparable contributions to the enhancement of stabilizing interactions. As for an application, a chemical enhancement mechanism of the surface-enhanced Raman scattering (SERS) procedure of cysteine by silver clusters was also analyzed.

Impact of metals on (star)dust chemistry: a laboratory astrophysics approach

Laboratory experiments are essential in exploring the mechanisms involved in stardust formation. One key question is how a metal is incorporated into dust for an environment rich in elements involved in stardust formation (C, H, O, Si). To address experimentally this question we have used a radiofrequency cold plasma reactor in which cyclic organosilicon dust formation is observed. Metallic (silver) atoms were injected in the plasma during the dust nucleation phase to study their incorporation in the dust. The experiments show formation of silver nanoparticles (~15 nm) under conditions in which organosilicon dust of size 200 nm or less is grown. The presence of AgSiO bonds, revealed by infrared spectroscopy, suggests the presence of junctions between the metallic nanoparticles and the organosilicon dust. Even after annealing we could not conclude on the formation of silver silicates, emphasizing that most of silver is included in the metallic nanoparticles. The molecular analysis performed by laser mass spectrometry exhibits a complex chemistry leading to a variety of molecules including large hydrocarbons and organometallic species. In order to gain insights into the involved chemical molecular pathways, the reactivity of silver atoms/ions with acetylene was studied in a laser vaporization source. Key organometallic species, Ag n C2H m (n=1-3; m=0-2), were identified and their structures and energetic data computed using density functional theory. This allows us to propose that molecular Ag-C seeds promote the formation of Ag clusters but also catalyze hydrocarbon growth. Throughout the article, we show how the developed methodology can be used to characterize the incorporation of metal atoms both in the molecular and dust phases. The presence of silver species in the plasma was motivated by objectives finding their application in other research fields than astrochemistry. Still, the reported methodology is a demonstration laying down the ground for future studies on metals of astrophysical interest such as iron.

Photocatalytic Surface Restructuring in Individual Silver Nanoparticles

Light absorption and scattering by metal nanoparticles can drive catalytic reactions at their surface via the generation of hot charge carriers, elevated temperatures, and focused electromagnetic fields. These photoinduced processes can substantially alter the shape, surface structure, and oxidation state of surface atoms of the nanoparticles and therefore significantly modify their catalytic properties. Information on such local structural and chemical change in plasmonic nanoparticles is however blurred in ensemble experiments, due to the typical large heterogeneity in sample size and shape distributions. Here, we use single-particle dark-field and Raman scattering spectroscopy to elucidate the reshaping and surface restructuring of individual silver nanodisks under plasmon excitation and during photocatalytic CO₂ hydrogenation. We show that silver nanoparticles reshape significantly in inert N₂ atmosphere, due to photothermal effects. Furthermore, by collecting the inelastic scattering during laser irradiation in a reducing gas environment, we observe intermittent light emission from silver clusters transiently formed at the nanoparticle surface. These clusters are likely to modify the photocatalytic activity of silver nanodisks and to enable detection of reaction products by enhancing their Raman signal. Our results highlight the dynamic nature of the catalytic surface of plasmonic silver nanoparticles and demonstrate the power of single-particle spectroscopic techniques to unveil their structure–activity relationship both in situ and in real time.

Thermal Stability of Ag13- Clusters Studied by Ab Initio Molecular Dynamics Simulations

Identification of the geometric structures of silver clusters is of great importance in future nanotechnologies due to their superior properties. Nevertheless, some ground-state structures are still in academic debate, partly because the experiments and theoretical calculations are not performed at the same temperatures. For example, silver clusters usually have compact configurations. However, a combined experimental and theoretical study proposed that the most stable structure of Ag₁₃^- had a two-coordinated atom. By using the CALYPSO approach for the global minima search followed by first-principles calculations, we discovered that a more compact trilayer Ag₁₃^- cluster was the ground state, in accordance with another three works published recently. In addition, its O₂ adsorption structure is also energetically favored. By tracing characteristic bond changes in ab initio molecular dynamics (MD) simulations, we confirmed that, compared with other isomers, this trilayer structure and its O₂ adsorption structure also had the highest thermal stability. This work emphasized the thermal stability concept in theoretical calculations, which may be a necessary supplement to explain the experimental observations on cluster science.

Dispersion-Corrected Density Functional Theory Study of the Noncovalent Complexes Formed with Imidazo[1,2-a]pyrazines Adsorbed onto Silver Clusters

Imidazo[1,2-a]pyrazines are cyclic amidine-type compounds composed of α-amino acid residues. A full structural identification of these molecules constitutes an analytical challenge, especially when imidazo[1,2-a]pyrazines are obtained from physical processes (e.g., sublimation and pyrolysis of amino acids). A valuable source of molecular information can be obtained from absorption spectroscopies and related techniques encompassing the use of metallic substrates. The aim of this study is to provide new knowledge and insights into the noncovalent intermolecular interactions between imidazo[1,2-a]pyrazines and two Ag _n (n = 4 and 20) clusters using density functional theory (DFT) methods. Semiempirical DFT dispersion (DFT-D) corrections were addressed using Grimme's dispersion (GD2) and Austin-Petersson-Frisch (APF) functionals in conjunction with the 6-31+G(d,p) + LANL2DZ mixed basis set. These DFT-D methods describe strong interactions; besides, in all cases, the APF dispersion (APF-D) energies of interaction appear to be consistently overestimated. In comparison with B3LYP calculations, the mean values for the difference in the energies of interaction calculated are 2.25 (GD2) and 6.24 (APF-D) kcal mol^-1 for Ag₄-molecules, and 2.30 (GD2) and 8.53 (APF-D) kcal mol^-1 for Ag₂₀-molecules. The effect of applying GD2 and APF-D corrections to the noncovalent complexes is nuanced in the intermolecular distances calculated, mainly in the Ag···N(amidine) bonding, which appears to play the most important role for the adsorptive process. Selective enhancement and considerable red shifts for Raman vibrations suggest strong interactions, whereas a charge redistribution involving the metallic substrate and the absorbate leads to a significant rearrangement of frontier molecular orbitals mainly in the Ag₂₀-molecule complexes. Finally, time-dependent DFT calculations were carried out to access the orbital contributions to each of the transitions observed in the absorption spectrum. The corresponding UV-vis spectra involve transitions in the visible region at around 400 and 550 nm for the Ag₄-molecule and the Ag₂₀-molecule complexes, respectively.

Dissociative chemisorption of O2 on Agn and Agn−1Ir (n = 3–26) clusters: a first-principle study

Lower dissociation barriers and higher reaction rates of O₂ on doped Ag_n−1Ir clusters, and a gradually weakened dopant effect.

The adsorption and growth of Agn (n = 1–4) clusters on cubic, monoclinic, and tetragonal ZrO2 surfaces: a first-principles study

The adsorption and growth of silver clusters on different zirconia surfaces.

Study of odd-even effects in physisorption and chemisorption of Ar, N2, O2 and NO on open shell Ag11-13+ clusters by means of self-consistent van der Waals density functional calculations

We have studied the adsorption and coadsorption properties of one or more X = Ar, N₂, O₂, and NO adsorbates on cationic silver clusters Ag_11-13⁺, whose sizes are in the open shell region of metal clusters, aiming to understand the observed odd-even effects in the abundance spectra of Ag_11-13⁺·mX complexes. All calculations were performed self-consistently using a non-local van der Waals correlation functional, covering the different nature of the interactions between the silver substrate and the several adsorbates, which range from dispersion (London) forces for Ar, non covalent π-π interactions for N₂, charge-transfer interactions for O₂ and NO, and the covalent Ag-Ag bond in the nude silver cluster. Despite the wide interval of adsorption energies, spanning two orders of magnitude, we have been able to explain the following experimental facts. For X = Ar, N₂, and O₂ reactions with Ag_11-13⁺, it was observed in the mass spectra an abundance peak at n = 12 [M. Schmidt, et al., ChemPhysChem, 2015, 16, 855]. In addition it was observed the competitive adsorption of two or more N₂ molecules, and the cooperative effect of adsorbing N₂ together with O₂ molecules. For X = NO, an abundance peak at n = 12 has been also observed [J. Ma, et al., Phys. Chem. Chem. Phys., 2016, 18, 12819]. We find that the main factors determining these properties are the different core motifs of the cluster geometry (pentagonal bipiramid for Ag₁₁⁺ and Ag₁₃⁺, but triangular prism for Ag₁₂⁺) and, on the other hand, the odd number of valence electrons for Ag₁₂⁺, leading to a smaller HOMO-LUMO gap than those of its neighbours. Further details about the preferred adsorption sites, dipole moments, and dipole polarizabilities are also discussed.

Molecular Docking as a Promising Predictive Model for Silver Nanoparticle-Mediated Inhibition of Cytochrome P450 Enzymes

Cytochrome P450 (CYP) enzymes are responsible for oxidative metabolisms of a large number of xenobiotics. In this study, we investigated interactions of silver nanoparticles (AgNPs) and silver ions (Ag⁺) with six CYP isoforms, namely, CYP1A2, CYP2C9, CYP2C19, CYP2D6, CYP2E1, and CYP3A4, within CYP-specific inhibitor-binding pockets by molecular docking and quantum mechanical (QM) calculations. The docking results revealed that the Ag₃ cluster, not Ag⁺, interacted with key amino acids of CYP2C9, CYP2C19, and CYP2D6 within a distance of about 3 Å. Moreover, the QM analysis confirmed that the amino acid residues of these CYP enzymes strongly interacted with the Ag₃ cluster, providing more insight into the mechanism of the potential inhibition of CYP enzyme activities. Interestingly, these results are consistent with previous in vitro data indicating that AgNPs inhibited activities of CYP2C and CYP2D in rat liver microsomes. It is suggested that the Ag₃ cluster is a minimal unit of AgNPs for in silico modeling. In summary, we demonstrated that molecular docking, together with QM analysis, is a promising tool to predict AgNP-mediated CYP inhibition. These methods are useful for deeper understanding of reaction mechanisms and could be used for other nanomaterials.

Electrostatics-Assisted Building-Up Procedure for Capturing Energy Minima of Metal Clusters: Test Case of Agn Clusters

Global geometry optimization of metal clusters is an important problem in nanophysics. The starting geometries of the clusters generated with empirical or other model potentials are generally optimized further by density functional theory (DFT)-based energy minimization. For this purpose, several algorithms such as simulated annealing, genetic algorithms, basin hopping, etc. are used. Our building-up procedure generates putative lower-energy structures of metal (M) clusters, M_n+1, M_n+2, etc., by anchoring one or more metal atoms in the vicinity of the minima of the molecular electrostatic potential (MESP) of M_n. Here, we report an application of this method to Ag_n clusters, for 5 ≤ n ≤ 20, followed up by DFT-based geometry optimization, generating several lower-energy structures than those reported in the literature. New low-energy isomers are obtained by applying the same procedure to the test case of mixed-metal clusters, Ni_nAg_m, for n + m = 4 and 5. In conclusion, our MESP-based building-up procedure offers a new general methodology for generating lower-energy geometries of metal clusters.

Comment on "Theoretical Investigations on Geometrical and Electronic Structures of Silver Clusters"

We would like to make a few comments on the results reported in a recent paper in this journal by Tsuneda, J. Comput. Chem, 2019, 40, 206. The structures of some pure silver clusters were not correctly assigned. © 2019 Wiley Periodicals, Inc.

Adsorption Kinetics of Nitrogen Molecules on Size-Selected Silver Cluster Cations

We present adsorption processes of dinitrogen on size-selected silver cluster cations, Ag n + (n = 1–10), studied by kinetics measurement using an ion trap. The cluster ions showed sequential adsorption of N2 molecules when the ion trap was cooled down to 105 K, excluding n = 8 and 9 that were exceptionally inactive at this temperature. Termolecular rate coefficients of each adsorption step are determined by analyzing time-dependent changes in the reactant and product ion signals. The first-step rate coefficients were found to increase exponentially from n = 1 to 7 due to increased internal degrees of freedom at larger sizes, which are favorable for accommodating the adsorption energy in a free cluster. In contrast, the adsorption rate turned to decrease for n > 7 due to weaker binding of dinitrogen as revealed by density-functional-theory (DFT) calculation. Adsorption sites on Ag n + are further discussed on the basis of the maximum number of adsorbing N2 molecules observed in the experiment.

First-principles modeling of GaN(0001)/water interface: Effect of surface charging

The accumulation properties of photogenerated carriers at the semiconductor surface determine the performance of photoelectrodes. However, to the best of our knowledge, there are no computational studies that methodically examine the effect of "surface charging" on photocatalytic activities. In this work, the effect of excess carriers at the semiconductor surface on the geometric and electronic structures of the semiconductor/electrolyte interface is studied systematically with the aid of first-principles calculations. We found that the number of water molecules that can be dissociated follows the "extended" electron counting rule; the dissociation limit is smaller than that predicted by the standard electron counting rule (0.375 ML) by the number of excess holes at the interface. When the geometric structure of the GaN/water interface obeys the extended electron counting rule, the Ga-originated surface states are removed from the bandgap due to the excess holes and adsorbates, and correspondingly, the Fermi level becomes free from pinning. Clearly, the excess charge has a great impact on the interface structure and most likely on the chemical reactions. This study serves as a basis for further studies on the semiconductor/electrolyte interface under working conditions.

Structural properties of sub-nanometer metallic clusters

At the nanoscale, the investigation of structural features becomes fundamental as we can establish relationships between cluster geometries and their physicochemical properties. The peculiarity lies in the variety of shapes often unusual and far from any geometrical and crystallographic intuition clusters can assume. In this respect, we should treat and consider nanoparticles as a new form of matter. Nanoparticle structures depend on their size, chemical composition, ordering, as well as external conditions e.g. synthesis method, pressure, temperature, support. On top of that, at finite temperatures nanoparticles can fluctuate among different structures, opening new and exciting horizons for the design of optimal nanoparticles for advanced applications. This article aims to overview geometrical features of transition metal clusters and of their various rearrangements.

Electronic processes in NO dimerization on Ag and Cu clusters: DFT and MRMP2 studies

Experimentally observed NO dimerization on Cu and Ag surfaces is surprising because binding energy of NO dimer is very small in gas phase. MRMP2, MP2 to MP4, CCSD(T), and DFT studies of NO dimerization on Ag₂ and Cu₂ clusters disclosed that the CCSD(T) method could be applied to this reaction on Ag₂ and Cu₂ unlike NO dimerization in gas phase which exhibits significantly large nondynamical electron correlation effect. Charge-transfer (CT) from Ag₂ and Cu₂ to NO moieties plays important role in NN bond formation between two NO molecules. This CT considerably decreases nondynamical correlation effect. Also, the DFT method could be applied to this NO dimerization, if appropriate DFT functional is used; all pure functionals examined here and most of the hybrid functionals underestimated the activation barrier (E_a ), while only ωB97X provided E_a similar to CCSD(T)-calculated value. NO dimerization on similar Cu₂ and Cu₅ needs moderately larger E_a than those on Ag₂ and Ag₅ , because frontier orbital participating in the CT exists at lower energy in Cu₂ and Cu₅ than in Ag₂ and Ag₅ . The E_a decreases in the order Ag₂ >> Ag₃₈ > Ag₇ ∼ Ag₅ and the reaction energy (ΔE) is positive (endothermic) in Ag₂ but significantly negative in Ag₃₈ , Ag₇ , and Ag₅ , indicating that various Ag clusters could be effective for NO dimerization except for Ag₂ . The decreasing order of E_a and increasing order of exothermicity are attributed to increasing order of the frontier orbital energy of Ag₂ < Ag₃₈ < Ag₇ ∼ Ag₅ . © 2018 Wiley Periodicals, Inc.

Regium bonds between Mn clusters (M = Cu, Ag, Au and n = 2-6) and nucleophiles NH3 and HCN

The most stable geometries of the coinage metal (or regium) atom (Cu, Ag, Au) clusters Mn for n up to 6 are all planar, and adopt the lowest possible spin multiplicity. Clusters with even numbers of M atoms are thus singlets, while those with odd n are open-shell doublets. Examination of the molecular electrostatic potential (MEP) of each cluster provides strong indications of the most likely site of attack by an approaching nucleophile, generally one of two positions. A nucleophile (NH3 or HCN) most favorably approaches one particular M atom of each cluster, rather than a bond midpoint or face. In the closed-shell clusters, the interaction energies are highly dependent upon the intensity of the MEP, but this correlation fades for the open-shell systems studied in this work. The strength of the interaction is also closely related to the basicity of the nucleophile. Regium bond energies can be more than 30 kcal mol-1 and tend to follow the Au > Cu > Ag order. These interaction energies are in large part derived from Coulombic attraction, with a smaller orbital interaction contribution.

Density-functional tight-binding approach for metal clusters, nanoparticles, surfaces and bulk: application to silver and gold

Density-functional based tight-binding (DFTB) is an efficient quantum mechanical method that can describe a variety of systems, going from organic and inorganic compounds to metallic and hybrid materials. The present topical review addresses the ability and performance of DFTB to investigate energetic, structural, spectroscopic and dynamical properties of gold and silver materials. After a brief overview of the theoretical basis of DFTB, its parametrization and its transferability, we report its past and recent applications to gold and silver systems, including small clusters, nanoparticles, bulk and surfaces, bare and interacting with various organic and inorganic compounds. The range of applications covered by those studies goes from plasmonics and molecular electronics, to energy conversion and surface chemistry. Finally, perspectives of DFTB in the field of gold and silver surfaces and NPs are outlined.

Elucidation of the molecular and electronic structures of some magic silver clusters Agn (n = 8, 18, 20)

Density functional theory (DFT) calculations were carried out to explore the geometric, spectroscopic, and electronic properties of three magic silver clusters Ag_n (n = 8, 18, and 20) in detail. The computed results show that the global minima of these clusters are compact, near-spherical structures, while other low-lying isomers exhibit oblate or prolate shapes. Vertical ionization energies for the low-lying isomers were also computed and assigned with respect to available experimental values. Although several isomers were predicted to have similar energies, their electronic and vibrational signatures were quite distinctive, meaning that they could be used as fingerprint signals to distinguish between isomers. In addition, the electronic structures of these systems were explored using the phenomenological shell model. Calculations for the coinage metal clusters M₂₀ (M = Cu, Ag, Au) indicated that the structures and properties of the Ag cluster are similar to those of the Cu cluster in that both Cu₂₀ and Ag₂₀ prefer a compact structure whereas Au₂₀ prefers to adopt a tetrahedral form. Graphical abstract Shell Orbitals of Ag₈ Cluster.

Probing the structural, electronic and magnetic properties of AgnSc (n = 1–16) clusters

The structural, electronic and magnetic properties of Ag_nSc (n = 1–16) clusters have been studied on the basis of density functional theory and the CALYPSO structure prediction method.

Probing the Structural, Electronic, and Magnetic Properties of Ag n V (n = 1-12) Clusters

The structural, electronic, and magnetic properties of Ag _n V (n = 1-12) clusters have been studied using density functional theory and CALYPSO structure searching method. Geometry optimizations manifest that a vanadium atom in low-energy Ag_nV clusters favors the most highly coordinated location. The substitution of one V atom for an Ag atom in Ag _n + 1 (n ≥ 5) cluster modifies the lowest energy structure of the host cluster. The infrared spectra, Raman spectra, and photoelectron spectra of Ag _n V (n = 1-12) clusters are simulated and can be used to determine the most stable structure in the future. The relative stability, dissociation channel, and chemical activity of the ground states are analyzed through atomic averaged binding energy, dissociation energy, and energy gap. It is found that V atom can improve the stability of the host cluster, Ag₂ excepted. The most possible dissociation channels are Ag _n V = Ag + Ag _n - 1V for n = 1 and 4-12 and Ag _n V = Ag₂ + Ag _n - 2V for n = 2 and 3. The energy gap of Ag _n V cluster with odd n is much smaller than that of Ag _n + 1 cluster. Analyses of magnetic property indicate that the total magnetic moment of Ag _n V cluster mostly comes from V atom and varies from 1 to 5 μ _B. The charge transfer between V and Ag atoms should be responsible for the change of magnetic moment.