Support effects on conical intersections of Jahn-Teller fluxional metal clusters on the sub-nanoscale.

The concept of fluxionality has been invoked to explain the enhanced catalytic properties of atomically precise metal clusters of subnanometer size. Cu3 isolated in the gas phase is a classical case of a fluxional metal cluster where a conical intersection leads to a Jahn-Teller (JT) distortion resulting in a potential energy landscape with close-lying multiminima and, ultimately, fluxional behavior. In spite of the role of conical intersections in the (photo)stability and (photo)catalytic properties of surface-supported atomic metal clusters, they have been largely unexplored. In this work, by applying a high-level multi-reference ab initio method aided with dispersion corrections, we analyze support effects on the conical intersection of Cu3 considering benzene as a model support molecule of carbon-based surfaces. We verify that the region around the conical intersection and the associated Jahn-Teller (JT) distortion is very slightly perturbed by the support when the Cu3 cluster approaches it in a parallel orientation: Two electronic states remain degenerate for a structure with C3 symmetry consistent with the D3h symmetry of unsupported Cu3 at the conical intersection. It extends over a one-dimensional seam that characterizes a physisorption minimum of the Cu3-benzene complex. The fluxionality of the Cu3 cluster, reflected in large fluctuations of relaxed Cu-Cu distances as a function of the active JT mode, is kept unperturbed upon complexation with benzene as well. In stark contrast, for the energetically favored perpendicular orientation of the Cu3 plane to the benzene ring plane, the conical intersection (CI) is located 12 100 cm-1 (∼1.5 eV) above the chemisorption minimum, with the fluxionality being kept at the CI's nearby and lost at the chemisorption well. The first excited state at the perpendicular orientation has a deep well (>4000 cm-1), being energetically closer to the CI. The transition dipole moment between ground and excited states has a significant magnitude, suggesting that the excited state can be observed through direct photo-excitation from the ground state. Besides demonstrating that the identity of an isolated Jahn-Teller metal cluster can be preserved against support effects at a physisorption state and lifted out at a chemisorption state, our results indicate that a correlation exists between conical intersection topography and fluxionality in the metal cluster's Cu-Cu motifs.

Stability and Reversible Oxidation of Sub-Nanometric Cu5 Metal Clusters: Integrated Experimental Study and Theoretical Modeling.

Sub-nanometer metal clusters have special physical and chemical properties, significantly different from those of nanoparticles. However, there is a major concern about their thermal stability and susceptibility to oxidation. In situ X-ray Absorption spectroscopy and Near Ambient Pressure X-ray Photoelectron spectroscopy results reveal that supported Cu5 clusters are resistant to irreversible oxidation at least up to 773 K, even in the presence of 0.15 mbar of oxygen. These experimental findings can be formally described by a theoretical model which combines dispersion-corrected DFT and first principles thermochemistry revealing that most of the adsorbed O2 molecules are transformed into superoxo and peroxo species by an interplay of collective charge transfer within the network of Cu atoms and large amplitude "breathing" motions. A chemical phase diagram for Cu oxidation states of the Cu5 -oxygen system is presented, clearly different from the already known bulk and nano-structured chemistry of Cu.

High-level ab initio evidence of bipyramidal Cu5 clusters as fluxional Jahn-Teller molecules.

The front cover artwork is provided by María Pilar de Lara-Castells, Head of the AbinitFot Group at IFF-CSIC (Madrid), Coordinator of the National Project "COSYES", and Chair of the COST Action CA21101 "COSY", and Alexander O. Mitrushchenkov from the Université Paris-Est. The image shows the connection between the Jahn-Teller effect featured by bypiramidal Cu5 clusters and the property of fluxionality. Cover design by Katarzyna Krupka. Read the full text of the Research Article at 10.1002/cphc.202300317.

Superfluid helium droplet-mediated surface-deposition of neutral and charged silver atomic species.

Experimental and theoretical work has delivered evidence of the helium nanodroplet-mediated synthesis and soft-landing of metal nanoparticles, nanowires, clusters, and single atoms on solid supports. Recent experimental advances have allowed the formation of charged metal clusters into multiply charged helium nanodroplets. The impact of the charge of immersed metal species in helium nanodroplet-mediated surface deposition is proved by considering silver atoms and cations at zero-temperature graphene as the support. By combining high-level ab initio intermolecular interaction theory with a full quantum description of the superfluid helium nanodroplet motion, evidence is presented that the fundamental mechanism of soft-deposition is preserved in spite of the much stronger interaction of charged species with surfaces, with high-density fluctuations in the helium droplet playing an essential role in braking them. Corroboration is also presented that the soft-landing becomes favored as the helium nanodroplet size increases.

Carbon vacancy-assisted stabilization of individual Cu5 clusters on graphene. Insights from ab initio molecular dynamics.

Recent advances in synthesis and characterization methods have enabled the controllable fabrication of atomically precise metal clusters (AMCs) of subnanometer size that possess unique physical and chemical properties, yet to be explored. Such AMCs have potential applications in a wide range of fields, from luminescence and sensing to photocatalysis and bioimaging, making them highly desirable for further research. Therefore, there is a need to develop innovative methods to stabilize AMCs upon surface deposition, as their special properties are lost due to sintering into larger nanoparticles. To this end, dispersion-corrected density functional theory (DFT-D3) and ab initio molecular dynamics (AIMD) simulations have been employed. Benchmarking against high-level post-Hartree-Fock approaches revealed that the DFT-D3 scheme describes very well the lowest-energy states of clusters of five and ten atoms, Cu5 and Cu10. AIMD simulations performed at 400 K illustrate how intrinsic defects of graphene sheets, carbon vacancies, are capable of confining individual Cu5 clusters, thus allowing for their stabilization. Furthermore, AIMD simulations provide evidence on the dimerization of Cu5 clusters on defect-free graphene, in agreement with the ab initio predictions of (Cu5)n aggregation in the gas phase. The findings of this study demonstrate the potential of using graphene-based substrates as an effective platform for the stabilization of monodisperse atomically precise Cu5 clusters.

Meta-stability through intermolecular interactions protecting the identity of atomic metal clusters: ab initio evidences in (Cu5-Cu5)n (n < 3) cases.

Recent developments in new synthesis techniques have allowed the production of precise monodisperse metal clusters composed of a few atoms. These atomic metal clusters (AMCs) often feature a molecule-like electronic structure, which makes their physical and chemical properties particularly interesting in nanotechnology. Regarding potential applications, there is a major concern about the sintering of AMCs in nanoparticles due to the loss of their special properties. In this work, multireference ab initio theory is applied to demonstrate the formation of coupled AMC-AMC clusters in which the AMC partners maintain their 'identity' to a large extent in terms of their initial structures and atomic Mulliken charges, and their further oligomerization.

Aggregation and support effects in the oxidation of fluxional atomic metal clusters. The paradigmatic Cu5 case.

The recent development of new synthesis techniques has allowed the production of monodisperse metal clusters composed of a few atoms. Follow-up experimental spectroscopic characterization has indicated the stability of these atomic metal clusters (AMCs). Despite the common assumption that the occurrence of an irreversible oxidation becomes more likely as the cluster size decreases, its quenching and reversible nature has been experimentally identified in the particular case of Cu5 clusters, making them paradigmatic. This work aims to address the influence of aggregation and the effects of a chemically inert carbon-based support on the oxidation of AMCs, considering the case of Cu5 as a model system. For this purpose, we present an extended first-principles study of the oxidation of Cu5-Cu5 and circumpyrene-supported Cu5, comparing it with that of unsupported Cu5, and combine dispersion-corrected density-functionals, first principles thermochemistry, and ab initio molecular dynamics (AIMD) simulations within an adiabatic approach. Our results indicate that a molecular chemisorption/desorption model is sensible upon consideration of aggregation and support effects in such a way that the predicted (p-T)-phase diagrams do not differ significantly from those obtained for unsupported Cu5. We also provide insights into the decoupling of the Cu5-Cu5 dimer into Cu5 sub-units through activated fluxional rotational motion, upon heating, as well as the adsorption of multiple O2 molecules at high oxygen gas pressures. Furthermore, numerical evidence shows the likelihood of a support-mediated mechanism leading to the dissociation of chemisorbed peroxo O22- species, delivering states with very similar energies to those characterized by molecular chemisorption. A Boltzmann-weighted average of the free energies of formation is computed as well, coming up with a diagram of the dominant copper oxidation states as a function of temperature and oxygen gas pressure.

A nuclear spin and spatial symmetry-adapted full quantum method for light particles inside carbon nanotubes: clusters of 3He, 4He, and para-H2.

We present a new nuclear spin and spatial symmetry-adapted full quantum method for light fermionic and bosonic particles under cylindrical carbon nanotube confinement. The goal is to address Fermi-Dirac and Bose-Einstein nuclear spin statistics on an equal footing and to deliver excited states with a similar accuracy to that of the ground state, implementing ab initio-derived potential models as well. The method is applied to clusters of up to four (three) 4He atoms and para-H2 molecules (3He atoms) inside a single-walled (1 nm diameter) carbon nanotube. Due to spin symmetry effects, the bound states energy landscape as a function of the angular momentum around the tube axis becomes much more complex and rich as the number of 3He atoms increase compared to the spinless 4He and para-H2 counterparts. Four bosonic 4He and para-H2 particles form pyramidal-like structures which are more compact as the particle mass and the strength of the inter-particle interaction increases. They feature stabilization of the collective rotational motion as bosonic quantum rings bearing persistent rotational motion and superfluid flow. Our results are brought together with two key experimental findings from the group of Jan-Peter Toennies: (1) the congestion of spectral profiles in doped 3He droplets as opposed to the case of 4He droplets (S. Gebenev, J. P. Toennies and A. F. Vilesov, Science, 1998, 279, 2083); (2) the onset of microscopic superfluidity in small doped clusters of para-H2 molecules (S. Grebenev, B. G. Sartakov, J. P. Toennies and A. F. Vilesov, Science, 2000, 289, 1532), but at the reduced dimensionality offered by the confinement inside carbon nanotubes.

Nonadiabatic Effects in the Molecular Oxidation of Subnanometric Cu5 Clusters.

The electronic structure of subnanometric clusters, far off the bulk regime, is still dominated by molecular characteristics. The spatial arrangement of the notoriously undercoordinated metal atoms is strongly coupled to the electronic properties of the system, which makes this class of materials particularly interesting for applications including luminescence, sensing, bioimaging, theranostics, energy conversion, catalysis, and photocatalysis. Opposing a common rule of thumb that assumes an increasing chemical reactivity with smaller cluster size, Cu5 clusters have proven to be exceptionally resistant to irreversible oxidation, i.e., the dissociative chemisorption of molecular oxygen. Besides providing reasons for this behavior in the case of heavy loading with molecular oxygen, we investigate the competition between physisorption and molecular chemisorption from the perspective of nonadiabatic effects. Landau-Zener theory is applied to the Cu5(O2)3 complex to estimate the probability for a switching between the electronic states correlating the neutral O2 + Cu5(O2)2 and the ionic O2- + (Cu5(O2)2)+ fragments in a diabatic representation. Our work demonstrates the involvement of strong nonadiabatic effects in the associated charge transfer process, which might be a common motive in reactions involving subnanometric metal structures.

A Path Integral Molecular Dynamics Simulation of a Harpoon-Type Redox Reaction in a Helium Nanodroplet.

We present path integral molecular dynamics (PIMD) calculations of an electron transfer from a heliophobic Cs2 dimer in its (3Σu) state, located on the surface of a He droplet, to a heliophilic, fully immersed C60 molecule. Supported by electron ionization mass spectroscopy measurements (Renzler et al., J. Chem. Phys.2016, 145, 181101), this spatially quenched reaction was characterized as a harpoon-type or long-range electron transfer in a previous high-level ab initio study (de Lara-Castells et al., J. Phys. Chem. Lett.2017, 8, 4284). To go beyond the static approach, classical and quantum PIMD simulations are performed at 2 K, slightly below the critical temperature for helium superfluidity (2.172 K). Calculations are executed in the NVT ensemble as well as the NVE ensemble to provide insights into real-time dynamics. A droplet size of 2090 atoms is assumed to study the impact of spatial hindrance on reactivity. By changing the number of beads in the PIMD simulations, the impact of quantization can be studied in greater detail and without an implicit assumption of superfluidity. We find that the reaction probability increases with higher levels of quantization. Our findings confirm earlier, static predictions of a rotational motion of the Cs2 dimer upon reacting with the fullerene, involving a substantial displacement of helium. However, it also raises the new question of whether the interacting species are driven out-of-equilibrium after impurity uptake, since reactivity is strongly quenched if a full thermal equilibration is assumed. More generally, our work points towards a novel mechanism for long-range electron transfer through an interplay between nuclear quantum delocalization within the confining medium and delocalized electronic dispersion forces acting on the two reactants.

Spectroscopy of a rotating hydrogen molecule in carbon nanotubes.

A first-principles study of the spectroscopy of a single hydrogen molecule rotating inside and outside of carbon nanotubes is presented. Density functional theory (DFT)-based symmetry-adapted perturbation theory (SAPT) is applied to analyze the influence of the rotation in the dispersionless and dispersion energy contributions to the adsorbate-nanotube interaction. A potential model for the H2-nanotube interaction is proposed and applied to derive the molecular energy levels of the rotating hydrogen molecule. The SAPT-based analysis shows that a subtle balance between the dispersionless and dispersion contributions is key in determining the angular dependence of the H2-nanotube interaction, which is strongly influenced by the diameter of the carbon nanotubes. As a consequence, the structure of molecular energy levels is very different in wide and narrow nanotubes with the diameter above and below 1 nanometer, respectively. Strong anisotropy effects lead to a rather constrained rotation of molecular hydrogen inside narrow nanotubes.

Ab initio design of light absorption through silver atomic cluster decoration of TiO2.

A first-principles study of the stability and optical response of subnanometer silver clusters Agn (n ≤ 5) on a TiO2(110) surface is presented. First, the adequacy of the vdW-corrected DFT-D3 approach is assessed using the domain-based pair natural orbital correlation DLPNO-CCSD(T) calculations along with the Symmetry-Adapted Perturbation Theory [SAPT(DFT)] applied to a cluster model. Next, using the DFT-D3 treatment with a periodic slab model, we analyze the interaction energies of the atomic silver clusters with the TiO2(110) surface. Finally, the hybrid HSE06 functional and a reduced density matrix treatment are applied to obtain the projected electronic density of states and photo-absorption spectra of the TiO2(110) surface, with and without adsorbed silver clusters. Our results show the stability of the supported clusters, the enhanced light absorbance intensity of the material upon their deposition, and the appearance of intense secondary broad peaks in the near-infrared and the visible regions of the spectrum, with positions depending on the size and shape of the supported clusters. The secondary peaks arise from the photo-induced transfer of electrons from intra-band valence 5s orbitals of the noble-metal cluster to 3d Ti band states of the supporting material.

Spatial quenching of a molecular charge-transfer process in a quantum fluid: the Csx-C60 reaction in superfluid helium nanodroplets.

A recent experimental study [Renzler et al., J. Chem. Phys., 2016, 145, 181101] on superfluid helium nanodroplets reported different reactivities for Cs atoms and Cs2 dimers with C60 fullerenes inside helium droplets. Alkali metal atoms and clusters are heliophobic, therefore typically residing on the droplet surface, while fullerenes are fully immersed into the droplet. In this theoretical study, which combines standard methods of computational chemistry with orbital-free helium density functional theory, we show that the experimental findings can be interpreted in the light of a quenched electron-transfer reaction between the fullerene and the alkali dopant, which is additionally hindered by a reaction barrier stemming from the necessary extrusion of helium upon approach of the two reactants.

Quantum confinement of molecular deuterium clusters in carbon nanotubes: ab initio evidence for hexagonal close packing.

An ab initio study of quantum confinement of deuterium clusters in carbon nanotubes is presented. First, density functional theory (DFT)-based symmetry-adapted perturbation theory is used to derive parameters for a pairwise potential model describing the adsorbate-nanotube interaction. Next, we analyze the quantum nuclear motion of N D2 molecules (N < 4) confined in carbon nanotubes using a highly accurate adsorbate-wave-function-based approach, and compare it with the motion of molecular hydrogen. We further apply an embedding approach and study zero-point energy effects on larger hexagonal and heptagonal structures of 7-8 D2 molecules. Our results show a preference for crystalline hexagonal close packing hcp of D2 molecules inside carbon nanotubes even at the cost of a reduced volumetric density within the cylindrical confinement.

Post-Hartree-Fock studies of the He/Mg(0001) interaction: Anti-corrugation, screening, and pairwise additivity.

The adsorption of noble gases on metallic surfaces represents a paradigmatic case of van-der-Waals (vdW) interaction due to the role of screening effects on the corrugation of the interaction potential [J. L. F. Da Silva et al., Phys. Rev. Lett. 90, 066104 (2003)]. The extremely small adsorption energy of He atoms on the Mg(0001) surface (below 3 meV) and the delocalized nature and mobility of the surface electrons make the He/Mg(0001) system particularly challenging, even for state-of-the-art vdW-corrected density functional-based (vdW-DFT) approaches [M. P. de Lara-Castells et al., J. Chem. Phys. 143, 194701 (2015)]. In this work, we meet this challenge by applying two different procedures. First, the dispersion-corrected second-order Möller-Plesset perturbation theory (MP2C) approach is adopted, using bare metal clusters of increasing size. Second, the method of increments [H. Stoll, J. Chem. Phys. 97, 8449 (1992)] is applied at coupled cluster singles and doubles and perturbative triples level, using embedded cluster models of the metal surface. Both approaches provide clear evidences of the anti-corrugation of the interaction potential: the He atom prefers on-top sites, instead of the expected hollow sites. This is interpreted as a signature of the screening of the He atom by the metal for the on-top configuration. The strong screening in the metal is clearly reflected in the relative contribution of successively deeper surface layers to the main dispersion contribution. Aimed to assist future dynamical simulations, a pairwise potential model for the He/surface interaction as a sum of effective He-Mg pair potentials is also presented, as an improvement of the approximation using isolated He-Mg pairs.

Nuclear Bound States of Molecular Hydrogen Physisorbed on Graphene: An Effective Two-Dimensional Model.

The interaction potential of molecular hydrogen physisorbed on a graphene sheet is evaluated using the ab initio-based periodic dlDF+Das scheme and its accuracy is assessed by comparing the nuclear bound-state energies supported by the H2(D2/HD)/graphite potentials with the experimental energies. The periodic dlDF+Das treatment uses DFT-based symmetry-adapted perturbation theory on surface cluster models to extract the dispersion contribution to the interaction whereas periodic dispersionless density functional (dlDF) calculations are performed to determine the dispersion-free counterpart. It is shown that the H2/graphene interaction is effectively two-dimensional (2D), with the distance from the molecule center-of-mass to the surface plane and the angle between the diatomic axis and the surface normal as the relevant degrees of freedom. The global potential minimum is found at the orthogonal orientation of the molecule with respect to the surface plane, with an equilibrium distance of 3.17 Å and a binding energy of -51.9 meV. The comparison of the binding energies shows an important improvement of our approach over the vdW-corrected DFT schemes when we are dealing with the very weak H2/surface interaction. Next, the 2D nuclear bound-state energies are calculated numerically. As a cross-validation of the interaction potential, the bound states are also determined for molecular hydrogen on the graphite surface (represented as an assembly of graphene sheets). With the largest absolute deviation being 1.7 meV, the theoretical and experimental energy levels compare very favorably.

Argon Interaction with Gold Surfaces: Ab Initio-Assisted Determination of Pair Ar-Au Potentials for Molecular Dynamics Simulations.

Global potentials for the interaction between the Ar atom and gold surfaces are investigated and Ar-Au pair potentials suitable for molecular dynamics simulations are derived. Using a periodic plane-wave representation of the electronic wave function, the nonlocal van-der-Waals vdW-DF2 and vdW-OptB86 approaches have been proved to describe better the interaction. These global interaction potentials have been decomposed to produce pair potentials. Then, the pair potentials have been compared with those derived by combining the dispersionless density functional dlDF for the repulsive part with an effective pairwise dispersion interaction. These repulsive potentials have been obtained from the decomposition of the repulsive interaction between the Ar atom and the Au2 and Au4 clusters and the dispersion coefficients have been evaluated by means of ab initio calculations on the Ar+Au2 complex using symmetry adapted perturbation theory. The pair potentials agree very well with those evaluated through periodic vdW-DF2 calculations. For benchmarking purposes, CCSD(T) calculations have also been performed for the ArAu and Ar+Au2 systems using large basis sets and extrapolations to the complete basis set limit. This work highlights that ab initio calculations using very small surface clusters can be used either as an independent cross-check to compare the performance of state-of-the-art vdW-corrected periodic DFT approaches or, directly, to calculate the pair potentials necessary in further molecular dynamics calculations.

Combining density functional and incremental post-Hartree-Fock approaches for van der Waals dominated adsorbate-surface interactions: Ag2/graphene.

A combined density functional (DFT) and incremental post-Hartree-Fock (post-HF) approach, proven earlier to calculate He-surface potential energy surfaces [de Lara-Castells et al., J. Chem. Phys. 141, 151102 (2014)], is applied to describe the van der Waals dominated Ag2/graphene interaction. It extends the dispersionless density functional theory developed by Pernal et al. [Phys. Rev. Lett. 103, 263201 (2009)] by including periodic boundary conditions while the dispersion is parametrized via the method of increments [H. Stoll, J. Chem. Phys. 97, 8449 (1992)]. Starting with the elementary cluster unit of the target surface (benzene), continuing through the realistic cluster model (coronene), and ending with the periodic model of the extended system, modern ab initio methodologies for intermolecular interactions as well as state-of-the-art van der Waals-corrected density functional-based approaches are put together both to assess the accuracy of the composite scheme and to better characterize the Ag2/graphene interaction. The present work illustrates how the combination of DFT and post-HF perspectives may be efficient to design simple and reliable ab initio-based schemes in extended systems for surface science applications.

Transferability and accuracy by combining dispersionless density functional and incremental post-Hartree-Fock theories: Noble gases adsorption on coronene/graphene/graphite surfaces.

The accuracy and transferability of the electronic structure approach combining dispersionless density functional theory (DFT) [K. Pernal et al., Phys. Rev. Lett. 103, 263201 (2009)] with the method of increments [H. Stoll, J. Chem. Phys. 97, 8449 (1992)], are validated for the interaction between the noble-gas Ne, Ar, Kr, and Xe atoms and coronene/graphene/graphite surfaces. This approach uses the method of increments for surface cluster models to extract intermonomer dispersion-like (2- and 3-body) correlation terms at coupled cluster singles and doubles and perturbative triples level, while periodic dispersionless density functionals calculations are performed to estimate the sum of Hartree-Fock and intramonomer correlation contributions. Dispersion energy contributions are also obtained using DFT-based symmetry-adapted perturbation theory [SAPT(DFT)]. An analysis of the structure of the X/surface (X = Ne, Ar, Kr, and Xe) interaction energies shows the excellent transferability properties of the leading intermonomer correlation contributions across the sequence of noble-gas atoms, which are also discussed using the Drude oscillator model. We further compare these results with van der Waals-(vdW)-corrected DFT-based approaches. As a test of accuracy, the energies of the low-lying nuclear bound states supported by the laterally averaged X/graphite potentials (X = (3)He, (4)He, Ne, Ar, Kr, and Xe) are calculated and compared with the best estimations from experimental measurements and an atom-bond potential model using the ab initio-assisted fine-tuning of semiempirical parameters. The bound-state energies determined differ by less than 6-7 meV (6%) from the atom-bond potential model. The crucial importance of including incremental 3-body dispersion-type terms is clearly demonstrated, showing that the SAPT(DFT) approach effectively account for these terms. With the deviations from the best experimental-based estimations smaller than 2.3 meV (1.9%), the accuracy of the combined DFT and post-HF incremental scheme is established for all the noble-gas atoms. With relative deviations smaller than 4% and 11%, good agreement is also achieved by applying the vdW-corrected DFT treatments PBE-D3 and vdW-DF2 for noble-gas atoms heavier than neon.