Nanothermodynamics of iron clusters: Small clusters, icosahedral and fcc-cuboctahedral structures.

The study of the thermodynamics and structures of iron clusters has been carried on, focusing on small clusters and initial icosahedral and fcc-cuboctahedral structures. Two combined tools are used. First, energy intervals are explored by the Monte Carlo algorithm, called σ-mapping, detailed in the work of Soudan et al. [J. Chem. Phys. 135, 144109 (2011), Paper I]. In its flat histogram version, it provides the classical density of states, gp(Ep), in terms of the potential energy of the system. Second, the iron system is described by a potential which is called "corrected EAM" (cEAM), explained in the work of Basire et al. [J. Chem. Phys. 141, 104304 (2014), Paper II]. Small clusters from 3 to 12 atoms in their ground state have been compared first with published Density Functional Theory (DFT) calculations, giving a complete agreement of geometries. The series of 13, 55, 147, and 309 atom icosahedrons is shown to be the most stable form for the cEAM potential. However, the 147 atom cluster has a special behaviour, since decreasing the energy from the liquid zone leads to the irreversible trapping of the cluster in a reproducible amorphous state, 7.38 eV higher in energy than the icosahedron. This behaviour is not observed at the higher size of 309 atoms. The heat capacity of the 55, 147, and 309 atom clusters revealed a pronounced peak in the solid zone, related to a solid-solid transition, prior to the melting peak. The corresponding series of 13, 55, and 147 atom cuboctahedrons has been compared, underscoring the unstability towards the icosahedral structure. This unstability occurs clearly in several steps for the 147 atom cluster, with a sudden transformation at a transition state. This illustrates the concerted icosahedron-cuboctahedron transformation of Buckminster Fuller-Mackay, which is calculated for the cEAM potential. Two other clusters of initial fcc structures with 24 and 38 atoms have been studied, as well as a 302 atom cluster. Each one relaxes towards a more stable structure without regularity. The 38 atom cluster exhibits a nearly glassy relaxation, through a cascade of six metastable states of long life. This behaviour, as that of the 147 atom cluster towards the amorphous state, shows that difficulties to reach ergodicity in the lower half of the solid zone are related to particular features of the potential energy landscape, and not necessarily to a too large size of the system. Comparisons of the cEAM iron system with published results about Lennard-Jones systems and DFT calculations are made. The results of the previous clusters have been combined with that of Paper II to plot the cohesive energy Ec and the melting temperature Tm in terms of the cluster atom number Nat. The Nat-1/3 linear dependence of the melting temperature (Pawlow law) is observed again for Nat > 150. In contrast, for Nat < 150, the curve diverges strongly from the Pawlow law, giving it an overall V-shape, with a linear increase of Tm when Nat goes from 55 to 13 atoms. Surprisingly, the 38 atom cluster is anomalously below the overall curve.

Nanothermodynamics of large iron clusters by means of a flat histogram Monte Carlo method.

The thermodynamics of iron clusters of various sizes, from 76 to 2452 atoms, typical of the catalyst particles used for carbon nanotubes growth, has been explored by a flat histogram Monte Carlo (MC) algorithm (called the σ-mapping), developed by Soudan et al. [J. Chem. Phys. 135, 144109 (2011), Paper I]. This method provides the classical density of states, gp(Ep) in the configurational space, in terms of the potential energy of the system, with good and well controlled convergence properties, particularly in the melting phase transition zone which is of interest in this work. To describe the system, an iron potential has been implemented, called "corrected EAM" (cEAM), which approximates the MEAM potential of Lee et al. [Phys. Rev. B 64, 184102 (2001)] with an accuracy better than 3 meV/at, and a five times larger computational speed. The main simplification concerns the angular dependence of the potential, with a small impact on accuracy, while the screening coefficients S(ij) are exactly computed with a fast algorithm. With this potential, ergodic explorations of the clusters can be performed efficiently in a reasonable computing time, at least in the upper half of the solid zone and above. Problems of ergodicity exist in the lower half of the solid zone but routes to overcome them are discussed. The solid-liquid (melting) phase transition temperature T(m) is plotted in terms of the cluster atom number N(at). The standard N(at)(-1/3) linear dependence (Pawlow law) is observed for N(at) >300, allowing an extrapolation up to the bulk metal at 1940 ±50 K. For N(at) <150, a strong divergence is observed compared to the Pawlow law. The melting transition, which begins at the surface, is stated by a Lindemann-Berry index and an atomic density analysis. Several new features are obtained for the thermodynamics of cEAM clusters, compared to the Rydberg pair potential clusters studied in Paper I.

A new Monte Carlo method for getting the density of states of atomic cluster systems.

A novel Monte Carlo flat histogram algorithm is proposed to get the classical density of states in terms of the potential energy, g(E(p)), for systems with continuous variables such as atomic clusters. It aims at avoiding the long iterative process of the Wang-Landau method and controlling carefully the convergence, but keeping the ability to overcome energy barriers. Our algorithm is based on a preliminary mapping in a series of points (called a σ-mapping), obtained by a two-parameter local probing of g(E(p)), and it converges in only two subsequent reweighting iterations on large intervals. The method is illustrated on the model system of a 432 atom cluster bound by a Rydberg type potential. Convergence properties are first examined in detail, particularly in the phase transition zone. We get g(E(p)) varying by a factor 10(3700) over the energy range [0.01 < E(p) < 6000 eV], covered by only eight overlapping intervals. Canonical quantities are derived, such as the internal energy U(T) and the heat capacity C(V)(T). This reveals the solid to liquid phase transition, lying in our conditions at the triple point. This phase transition is further studied by computing a Lindemann-Berry index, the atomic cluster density n(r), and the pressure, demonstrating the progressive surface melting at this triple point. Some limited results are also given for 1224 and 4044 atom clusters.

Molecular dynamics study of the hydration of lanthanum(III) and europium(III) including many-body effects.

Lanthanides complexes are widely used as contrast agents in magnetic resonance imaging (MRI) and are involved in many fields such as organic synthesis, catalysis, and nuclear waste management. The complexation of the ion by the solvent or an organic ligand and the resulting properties (for example the relaxivity in MRI) are mainly governed by the structure and dynamics of the coordination shells. All of the MD approaches already carried out for the lanthanide(III) hydration failed due to the lack of accurate representation of many-body effects. We present the first molecular dynamics simulation including these effects that accounts for the experimental results from a structural and dynamic (water exchange rate) point of view.

Can we understand the different coordinations and structures of closed-shell metal cation-water clusters?

We present a model potential for studying M(q+)(H(2)O)(n=1,9) clusters where M stands for either Na(+), Cs(+), Ca(2+), Ba(2+), or La(3+). The potential energy surfaces (PES) are explored by the Monte Carlo growth method. The results for the most significant equilibrium structures of the PES as well as for energetics are favorably compared to the best ab initio calculations found in the literature and to experimental results. Most of these complexes have a different coordination number in cluster compared to experimental results in solution or solid phase. An interpretation of the coordination number in clusters is given. In order to well describe the transition between the first hydration sphere and the second one we show that an autocoherent treatment of the electric field is necessary to correctly deal with polarization effects. We also explore the influence of the cation properties (charge, size, and polarizability) on both structures and coordination number in clusters, as well as the meaning of the second hydration sphere. Such an approach shows that the leading term in the interaction energy for a molecule in the second hydration sphere is an electrostatic attraction to the cation and not a hydrogen bond with the water molecules in the first hydration sphere.