Global geometry optimization of atomic clusters using a modified genetic algorithm in space‐fixed coordinates

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.

Structure elucidation and construction of isomerisation pathways in small to moderate-sized (6-27) MgO nanoclusters: an adaptive mutation simulated annealing based analysis with quantum chemical calculations

Determination of global minimum structures and elucidation of reaction paths or minimum energy paths between low-lying minima are of great chemical importance. To that end, we have used our own Adaptive Mutation Simulated Annealing method to determine the global minimum and the minimum energy paths for various isomerisation reactions for small to moderate-sized (MgO)n (n = 6-27) clusters, using the Born-Mayer potential with suitable parameter values. The minimum energy structures obtained by us match well with previously reported data and are used as guess structures for further optimisation at the DFT level (using the B3LYP functional and DGDZVP basis set). Our optimised structures are found to match very well with the further DFT optimised structures, where the comparison is done by determining the root mean square deviation values as well as the radial distribution function profiles. A scheme is proposed to determine the minimum energy paths for isomerisation reactions for some cluster sizes where the transition state/s obtained by us, at very low computational cost, match well with those obtained from further optimisation using DFT calculations. We have shown the efficacy of our method in determining the reaction pathways, even for cases that involve multi-step reactions.

A charge optimized many-body potential for iron/iron-fluoride systems

A classical interatomic potential for iron/iron-fluoride systems is developed in the framework of the charge optimized many-body (COMB) potential. This interatomic potential takes into consideration the effects of charge transfer and many-body interactions depending on the chemical environment. The potential is fitted to a training set composed of both experimental and ab initio results of the cohesive energies of several Fe and FeF2 crystal phases, the two fluorine molecules F2 and the F2-1 dissociation energy curve, the Fe and FeF2 lattice parameters of the ground state crystalline phase, and the elastic constants of the body centered cubic Fe structure. The potential is tested in an NVT ensemble for different initial structural configurations as the crystal ground state phases, F2 molecules, iron clusters, and iron nanospheres. In particular, we model the FeF2/Fe bilayer and multilayer interfaces, as well as a system of square FeF2 nanowires immersed in an iron solid. It has been shown that there exists a reordering of the atomic positions for F and Fe atoms at the interface zone; this rearrangement leads to an increase in the charge transfer among the atoms that make the interface and put forward a possible mechanism of the exchange bias origin based on asymmetric electric charge transfer in the different spin channels.

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.

Global Materials Structure Search with Chemically Motivated Coordinates

Identification of relevant reaction pathways in ever more complex composite materials and nanostructures poses a central challenge to computational materials discovery. Efficient global structure search, tailored to identify chemically relevant intermediates, could provide the necessary first-principles atomistic insight to enable a rational process design. In this work we modify a common feature of global geometry optimization schemes by employing automatically generated collective curvilinear coordinates. The similarity of these coordinates to molecular vibrations enhances the generation of chemically meaningful trial structures for covalently bound systems. In the application to hydrogenated Si clusters, we concomitantly observe a significantly increased efficiency in identifying low-energy structures and exploit it for an extensive sampling of potential products of silicon-cluster soft landing on Si(001) surfaces.

Use of density functional theory method to calculate structures of neutral carbon clusters Cn (3 ≤ n ≤ 24) and study their variability of structural forms

In this work, we present modifications to the well-known basin hopping (BH) optimization algorithm [D. J. Wales and J. P. Doye, J. Phys. Chem. A 101, 5111 (1997)] by incorporating in it the unique and specific nature of interactions among valence electrons and ions in carbon atoms through calculating the cluster's total energy by the density functional tight-binding (DFTB) theory, using it to find the lowest energy structures of carbon clusters and, from these optimized atomic and electronic structures, studying their varied forms of topological transitions, which include a linear chain, a monocyclic to a polycyclic ring, and a fullerene/cage-like geometry. In this modified BH (MBH) algorithm, we define a spatial volume within which the cluster's lowest energy structure is to be searched, and introduce in addition a cut-and-splice genetic operator to increase the searching performance of the energy minimum than the original BH technique. The present MBH/DFTB algorithm is, therefore, characteristically distinguishable from the original BH technique commonly applied to nonmetallic and metallic clusters, technically more thorough and natural in describing the intricate couplings between valence electrons and ions in a carbon cluster, and thus theoretically sound in putting these two charged components on an equal footing. The proposed modified minimization algorithm should be more appropriate, accurate, and precise in the description of a carbon cluster. We evaluate the present algorithm, its energy-minimum searching in particular, by its optimization robustness. Specifically, we first check the MBH/DFTB technique for two representative carbon clusters of larger size, i.e., C60 and C72 against the popular cut-and-splice approach [D. M. Deaven and K. M. Ho, Phys. Rev. Lett. 75, 288 (1995)] that normally is combined with the genetic algorithm method for finding the cluster's energy minimum, before employing it to investigate carbon clusters in the size range C3-C24 studying their topological transitions. An effort was also made to compare our MBH/DFTB and its re-optimized results carried out by full density functional theory (DFT) calculations with some early DFT-based studies.

Quantum simulations of the hydrogen molecule on ammonia clusters

Mixed ammonia-hydrogen molecule clusters [H2-(NH3)n] have been studied with the aim of exploring the quantitative importance of the H2 quantum motion in defining their structure and energetics. Minimum energy structures have been obtained employing genetic algorithm-based optimization methods in conjunction with accurate pair potentials for NH3-NH3 and H2-NH3. These include both a full 5D potential and a spherically averaged reduced surface mimicking the presence of a para-H2. All the putative global minima for n ≥ 7 are characterized by H2 being adsorbed onto a rhomboidal ammonia tetramer motif formed by two double donor and two double acceptor ammonia molecules. In a few cases, the choice of specific rhombus seems to be directed by the vicinity of an ammonia ad-molecule. Diffusion Monte Carlo simulations on a subset of the species obtained highlighted important quantum effects in defining the H2 surface distribution, often resulting in populating rhomboidal sites different from the global minimum one, and showing a compelling correlation between local geometrical features and the relative stability of surface H2. Clathrate-like species have also been studied and suggested to be metastable over a broad range of conditions if formed.

The structure and properties of small Pd clusters

The zero-temperature minimal energy structure of small free-standing Pd clusters (14≤N≤21, where N is the number of atoms in the cluster), their characteristics and their magnetic configurations are investigated. Results obtained using five different phenomenological many-body potentials (implemented in combination with a genetic algorithm search) are refined by means of various density functional theory (DFT) techniques. The agreement and differences between the results obtained with our procedure, using these five potentials, are displayed in detail. While phenomenological potentials yield values that approach the minimal energies of larger clusters, as compared with DFT results, they fail to predict the right symmetry group for some of the clusters with N>14. We find that the minimal energy configurations are not necessarily associated with high symmetry of the atomic arrangement. Actually, several cases of previously overlooked low symmetry structures turn out to have lower energies than more symmetric ones.

Global optimization of ionic Mg(n)F(2n) (n=1-30) clusters

The global optimization basin-hopping (BH) method has been used to locate the global minima (GM) of Mg(n)F(2n) (n=1-30) clusters using a Born-Mayer-type potential. Some of the GM were particularly difficult to find, requiring more than 1.5 x 10(4) BH steps. We have found that both the binding energy per MgF2 unit and the effective volume of the GM isomers increase almost linearly with n, and that cluster symmetry decreases with cluster size. The data derived from the BH runs reveal a growing density of local minima just above the GM as n increases. Despite this, the attraction basin around each GM is relatively large, since after all their atomic coordinates are randomly displaced by values as high as 2.0 bohrs, the perturbed structures, upon reoptimization, relax back to the GM in more than 50% of the cases (except for n=10 and 11). The relative stabilities derived from energy second differences suggest that n=8,10,13,15, and 20 are probably the magic numbers for these systems. Mass spectrum experiments would be very useful to clarify this issue.

Global minimum structure searches via particle swarm optimization

Novel implementation of the evolutionary approach known as particle swarm optimization (PSO) capable of finding the global minimum of the potential energy surface of atomic assemblies is reported. This is the first time the PSO technique has been used to perform global optimization of minimum structure search for chemical systems. Significant improvements have been introduced to the original PSO algorithm to increase its efficiency and reliability and adapt it to chemical systems. The developed software has successfully found the lowest-energy structures of the LJ(26) Lennard-Jones cluster, anionic silicon hydride Si(2)H(5) (-), and triply hydrated hydroxide ion OH(-) (H(2)O)(3). It requires relatively small population sizes and demonstrates fast convergence. Efficiency of PSO has been compared with simulated annealing, and the gradient embedded genetic algorithm.

Alternative search strategy for minimal energy nanocluster structures: The case of rhodium, palladium, and silver

An alternative strategy to find the minimal energy structure of nanoclusters is presented and implemented. We use it to determine the structure of metallic clusters. It consists in an unbiased search, with a global minimum algorithm: conformational space annealing. First, we find the minima of a many-body phenomenological potential to create a data bank of putative minima. This procedure assures us the generation of a set of cluster configurations of large diversity. Next, the clusters in this data bank are relaxed by ab initio techniques to obtain their energies and geometrical structures. The scheme is successfully applied to magic number 13 atom clusters of rhodium, palladium, and silver. We obtained minimal energy cluster structures not previously reported, which are different from the phenomenological minima. Moreover, they are not always highly symmetric, thus casting some doubt on the customary biased search scheme, which consists in relaxing with density functional theory global minima chosen among high symmetry structures obtained by means of phenomenological potentials.

Structure of the Na(x)Cl(x+1) (-) (x=1-4) clusters via ab initio genetic algorithm and photoelectron spectroscopy

The application of the ab initio genetic algorithm with an embedded gradient has been carried out for the elucidation of global minimum structures of a series of anionic sodium chloride clusters, Na(x)Cl(x+1) (-) (x=1-4), produced in the gas phase using electrospray ionization and studied by photoelectron spectroscopy. These are all superhalogen species with extremely high electron binding energies. The vertical electron detachment energies for Na(x)Cl(x+1) (-) were measured to be 5.6, 6.46, 6.3, and 7.0 eV, for x=1-4, respectively. Our ab initio gradient embedded genetic algorithm program detected the linear global minima for NaCl(2) (-) and Na(2)Cl(3) (-) and three-dimensional structures for the larger species. Na(3)Cl(4) (-) was found to have C(3v) symmetry, which can be viewed as a Na(4)Cl(4) cube missing a corner Na(+) cation, whereas Na(4)Cl(5) (-) was found to have C(4v) symmetry, close to a 3x3 planar structure. Excellent agreement between the theoretically calculated and the experimental spectra was observed, confirming the obtained structures and demonstrating the power of the developed genetic algorithm technique.

Stereographic projection path-integral simulations of (HF)n clusters

We perform several quantum canonical ensemble simulations of (HF)(n) clusters. The HF stretches are rigid, and the stereographic projection path-integral method is employed for the simulation in the resulting curved configuration space. We make use of the reweighted random series techniques to accelerate the convergence of the path-integral simulation with respect to the number of path coefficients. We develop and test estimators for the total energy and heat capacity based on a finite difference approach for non-Euclidean spaces. The quantum effects at temperatures below 400 K are substantial for all sizes. We observe interesting thermodynamic behaviors in the quantum simulations of the octamer and the heptamer.

A genetic algorithm based determination of the ground and excited (1Lb) state structure and the orientation of the transition dipole moment of benzimidazole

The structure of benzimidazole has been determined in the electronic ground and excited states using rotationally resolved electronic spectroscopy. The rovibronic spectra of four isotopomers and subsequently the structure of benzimidazole have been automatically assigned and fitted using a genetic algorithm based fitting strategy. The lifetimes of the deuterated isotopomers have been shown to depend on the position of deuteration. The angle of the transition dipole moment with the inertial a-axis could be determined to be -30 degrees. Structures and transition dipole moment orientation have been calculated at various levels of theory and were compared to the experimental results.

Parallel tempering simulations of the 13-center Lennard-Jones dipole-dipole cluster (muD=0-->0.5 a.u.)

We investigate the thermodynamic behavior of the thirteen center uniform Lennard-Jones dipole-dipole cluster [(LJDD)(13)] for a wide range of dipole moment strengths. We find a relatively wide range of potential parameters where solid-solid coexistence manifests itself. Using structural characterization methods we determine the shape of the few isomers that contribute to the solid-solid coexistence region. The thermal distributions of the size of the net dipole moment are broad even at the coldest temperatures of the simulation where the (LJDD)(13) cluster is solid.

Parameter space minimization methods: Applications to Lennard-Jones–dipole-dipole clusters

The morphology of the uniform Lennard-Jones-dipole-dipole cluster with 13 centers (LJDD)13 is investigated over a relatively wide range of values of the dipole moment. We introduce and compare several necessary modifications of the basin-hopping algorithm for global optimization to improve its efficiency. We develop a general algorithm for T=0 Brownian dynamics in curved spaces, and a graph theoretical approach necessary for the elimination of dissociated states. We find that the (LJDD)13 cluster has icosahedral symmetry for small to moderate values of the dipole moment. As the dipole moment increases, however, its morphology shifts to an hexagonal antiprism, and eventually to a ring.

On the encapsulation of nickel clusters by molecular nitrogen

The structures and energetic effects of molecular nitrogen adsorbates on nickel clusters are investigated using an extended Huckel model coupled with two models of the adsorbate-nickel interaction. The potential parameters for the adsorbates are chosen to mimic experimental information about the binding strength of nitrogen on both cluster and bulk surface phases of nickel. The first model potential is a simple Lennard-Jones interaction that leads to binding sites in holes defined by sets of near-neighbor nickel atoms. The second model potential has a simple three-body form that forces the model nitrogen adsorbates to bind directly to single nickel atoms. Significant rearrangement of the core nickel structures are found in both models. A disconnectivity graph analysis of the potential energy surfaces implies that the rearrangements arise from low transition state barriers and the small differences between available isomers in the nickel core.