On the brink of self-hydration: the water heptadecamer

For pure, neutral, isolated molecular clusters, (H2O)17 marks the transition from structures with all water molecules on the cluster surface to water self-hydration, i.e., cluster structures around one central water molecule. Getting this right with water model potentials turns out to be challenging. Even the best water potentials currently available, which reproduce collective properties very well, still deliver contradicting results for (H2O)17, when different low-energy isomers from global structure optimizations are examined. Interestingly, ab initio quantum chemistry also struggles with the only seemingly simple question if (H2O)17 is all-surface or water-centered. Hence, although the long history of water potential development may be entering its final phase, it is not quite finished yet.

Exploring self-organization of molecular tether molecules on a gold surface by global structure optimization

We employ nondeterministic global cluster structure optimization, based on the evolutionary algorithms paradigm, to model the self-assembly of complex molecules on a surface. As a real-life application example directly related to many recent experiments, we use this approach for the assembly of triazatriangulene "platform" molecules on the Au(111) surface. Without additional restrictions like spatial discretizations, coarse-graining or precalculated adsorption poses, and despite the proof-of-principle character of this study, we achieve satisfactory qualitative agreement with several experimental observations and can provide answers to questions that experiments on these species had left open so far. © 2019 Wiley Periodicals, Inc.

Are zinc clusters really amorphous? A detailed protocol for locating global minimum structures of clusters

We report the results of a conjoint experimental/theoretical effort to assess the structures of free-standing zinc clusters with up to 73 atoms. Experiment provides photoemission spectra for ZnN- cluster anions, to be used as fingerprints in structural assessment, as well as mass spectra for both anion and cation clusters. Theory provides both a detailed description of a novel protocol to locate global minimum structures of clusters in an efficient and reliable way, and its specific application to neutral and charged zinc clusters. Our methodology is based on the well-known hybrid EP-DFT (empirical potential-density functional theory) approach, in which the approximate potential energy surface generated by an empirical Gupta potential is first sampled with unbiased basin hopping simulations, and then a selection of the isomers so identified is re-optimized at a first-principles DFT level. The novelty introduced in our paper is a simple but efficient new recipe to obtain the best possible EP parameters for a given cluster system, with which the first step of the EP-DFT method is to be performed. Our method is able to reproduce experimental measurements at an excellent level for most cluster sizes, implying its ability to locate the true global minimum structures; meanwhile, if exactly the same method is applied based on the existing Gupta potential (fitted to bulk properties), it leads to wrong predicted structures with energies between 1 and 2 eV above the correct ones. Opposite to what was claimed in the past, our work unequivocally demonstrates that Zn clusters are not amorphous, and they rather adopt high symmetry structures for most sizes. We show that Zn clusters have a number of exotic, unprecedented structural and electronic properties which are not expected for clusters of a metallic element, and describe them in detail.

Cluster structures influenced by interaction with a surface

Systematic investigation of how cluster structures change upon interaction with a surface, using global structure optimization by evolutionary algorithms.

Global structure search for molecules on surfaces: Efficient sampling with curvilinear coordinates

Efficient structure search is a major challenge in computational materials science. We present a modification of the basin hopping global geometry optimization approach that uses a curvilinear coordinate system to describe global trial moves. This approach has recently been shown to be efficient in structure determination of clusters [C. Panosetti et al., Nano Lett. 15, 8044-8048 (2015)] and is here extended for its application to covalent, complex molecules and large adsorbates on surfaces. The employed automatically constructed delocalized internal coordinates are similar to molecular vibrations, which enhances the generation of chemically meaningful trial structures. By introducing flexible constraints and local translation and rotation of independent geometrical subunits, we enable the use of this method for molecules adsorbed on surfaces and interfaces. For two test systems, trans-β-ionylideneacetic acid adsorbed on a Au(111) surface and methane adsorbed on a Ag(111) surface, we obtain superior performance of the method compared to standard optimization moves based on Cartesian coordinates.

Theoretical vibrational spectra of OH(-)(H2O)2: the effect of quantum distribution and vibrational coupling

We performed ab initio path integral molecular dynamics simulations for the hydroxide-water cluster, OH(-)(H2O)2, at 50 K, 100 K, and 150 K to investigate its flexible structure. From our simulations, we found that nuclear quantum effects enhance hydroxide hydrogen atom inversion and the conformational change between isomers occurs by simultaneous rotation of the free hydrogen atom. We propose the importance of including the transition state conformer with C2 symmetry, for the description of this system at temperatures realized in predissociation experiments. Temperature dependence of relative populations of each conformer along with multidimensional vibrational calculations were used to simulate the vibrational spectra and compare with the experimental spectra of Johnson and coworkers. We assign the doublet peaks seen in the experiment at 2500 to 3000 cm(-1), as the mixture of the ionic hydrogen bonded OH stretching overtone, ionic hydrogen bonded OH bending overtone, and the combination band of the ionic hydrogen bonded OH stretch and bend, which are modulated by the van der Waals OO vibrations. We concluded that for OH(-)(H2O)2, the vibrational couplings between the ionic hydrogen bonded motion and floppy modes contribute to the broadening of peaks observed in the 2500 to 3000 cm(-1) region.