Binding modes of carboxylic acids on cobalt nanoparticles

Owing to their high saturation magnetisation, cobalt nanoparticles hold significant potential for the hyperthermia treatment of tumours. Covalent binding of carboxylic acids to the nanoparticles can induce biocompatibility, whilst also preventing the formation of surface oxides which reduce the magnetic properties of cobalt. Understanding the origin of the acid–metal interaction is key, yet probably the most experimentally challenging step, for the rational design of such entities. In this density functional theory study, we use static calculations to establish that a 57-atom Co cluster is the smallest model able to reproduce the adsorption behaviour of carboxylic acids, and ab initio metadynamics to obtain the structure and the free energy landscape for its interaction with valeric acid. Our simulations show that a bridging bidentate binding mode has a stronger affinity compared to monodentate binding, with energetically high transition barriers between the two. A chelate interaction mode of two carboxyl oxygen atoms can be formed as an intermediate. These results clarify the organic–inorganic interactions in the cobalt–acid sys-tem, providing a basis for the rational design of biocompatible metallic nanoparticles.

S e e h t t p://o r c a .cf. a c. u k/ p olici e s. h t ml fo r u s a g e p olici e s.Co py ri g h t a n d m o r al ri g h t s fo r p u blic a tio n s m a d e a v ail a bl e in ORCA a r e r e t ai n e d by t h e c o py ri g h t h ol d e r s .

Introduction
With increasing numbers of cancer cases suffering from drug resistance and severe side effects as corollaries of commonly used chemotherapy, the development of materials for hyperthermia as a treatment, where a moderate degree of tissue warming results in malignant cell destruction with minimal collateral damage and no drug consumption, is becoming a priority. 1The high saturation magnetisation and large anisotropy field of pure cobalt nanoparticles give them an intrinsic advantage within the strict requirements of hyperthermia over commonly tested iron oxides and ferrites. 2,3However, biocompatibility is the key for their integration into biomedicine, and while titanium implants owe their biocompatibility to thin oxide layers formed on the surface in the presence of air and water, 4,5 the same phenomenon is a substantial drawback when magnetic cobalt nanoparticles transform into non-magnetic oxide compounds. 6][9] Understanding the surface chemistry and the interaction between the surfactant and the nanoparticle therefore plays a major role in tuning the properties and the performance of the final nanocomposite. 10,11While many chemical reactions at nanoparticle surfaces are very well understood experimentally, with strategies for the control of the surfactant architecture developed accordingly, atomic-level details about bond formation and charge transfer remain difficult to explore.Atomic understanding of the underlaying interaction mechanisms requires for the comprehension of structure, morphology, and size, but also physiochemical properties and reactivity of both metal of interest and capping agent.3][14] With an hcpfcc phase transition occurring at temperatures around 450 1C, 15,16 t h eh c pp h a s ei st h eo n eo f interest from a biomedical standpoint.8][19][20][21][22] Similarly, the shape and size of the nanoparticles can strongly affect the material uptake and biodistribution as well as their reactivity.Since the main interest in hyperthermia is to maximise magnetisation and mobility while minimising oxidation, small spherical nanoparticles are preferred, as after ligand functionalisation the overall diameter should be enough to circulate throughout the blood vessel network for a sustained period of time with minimal surface-to-volume ratio to avoid reacting with the surroundings.][25] A variety of ligands have so far been tested to enhance biocompatibility, improve growth control, stabilise the particle dispersions, prevent agglomeration, provide functionality, and limit the surface reactions.By now it is known that metallic nanoparticles experience strong van der Waals and magnetic dipole interactions, hence only the surfactants which covalently bond to the metal can considerably promote their stability.Carboxylic acids can bind strongly to the surface metal atoms through carboxyl group, and they are therefore an excellent choice as capping agents.Among them, oleic acid (OA), which has a C18 tail with a cis-double-bond in the middle, is the most common, since it can form the kink necessary for effective stabilisation.OA forms a dense protective monolayer, thereby producing nanocomposites in a highly uniform size range. 7,26,279][30] If one metal atom binds with one oxygen atom, the mode is called monodentate.If both oxygen atoms are involved, they can form complexes with either one or two different metal atoms, referred to as chelate and bidentate, respectively.The final possible interaction is of ionic nature.Which motif is dominating for a specific ligand-metal pair depends on the properties of the metal surface and the structure of the adsorbed carboxylic acid.Experimentally, the actual binding interaction may be distinguished by infrared (IR) spectroscopy due to the characteristically shifted carbonyl band.
The bonding nature of OA can have a critical effect on the magnetisation 31 and a detailed knowledge of the type of interaction between OA and cobalt nanoparticles hence allows fine-tuning of the size, morphology, and magnetic properties.It is therefore a crucial step in the characterisation of this nanocomposite for biomedical applications.However, experimental results are, so far, inconclusive on the binding mode between OA and cobalt nanoparticles.Older investigations gave a uniform conclusion on oleic 12,27,32 and elaidic 33 acid chemisorption onto cobalt surfaces, yielding nanoparticles in a narrow distribution, but they did not provide the exact binding mode.To the best of our knowledge, there exist only two such works that suggested a type of interaction between the cobalt atoms and carboxylate head of oleic acid.According to Wu et al. 34 oleic acid is chemisorbed as a carboxylate onto the Co nanoparticles, and the two oxygen atoms are coordinated symmetrically to the Co atoms.They interpreted the shift of the carbonyl's IR band to be due to a bridging bidentate adsorption.Similarly, Ansari et al. 35 presented an IR spectrum that demonstrated the chemisorption of the acid as a carboxylate onto the nanoparticles.However, they attributed the shift of the carbonyl band to a monodentate interaction.
In this work, we have first performed static density functional theory (DFT) calculations to establish a 57-atom cluster as the minimum size that reproduces the adsorption properties of larger Co nanoparticles.Next, we have performed metadynamics simulations of valeric acid (VA), a convenient model for OA, on the Co cluster.We have efficiently explored the free energy landscape of the system, evaluating in particular the most relevant energy minima and the barriers appearing in the corresponding transition pathways.

Optimising the size of the cluster
The use of atomically flat surfaces can be appropriate when exploring the arrangement and properties of small molecules since the effects of irregular, curved, stepped, porous, and otherwise deformed nanoparticle surfaces can be largely excluded.However, changes in the reactivity become notable for large surfactants, i.e., when the nanoparticle diameter is less than 20 times the length of the surfactants in their extended conformation.For example, the length of a C18 alkyl ion is about 3 nm, and the nanoparticle diameter then needs to be at least 60 nm for curvature effects to become marginal.However, the maximum size of a nanoparticle allowed for biomedical applications should not exceed 20 nm, including the protective coating, which significantly magnifies the effects of surface curvature, edges, and vertices.To determine the optimal cluster size in terms of reactivity and computational cost, the adsorption characteristics of one oxygen atom, formic acid, and acetic acid have been investigated on hcp-shaped cobalt clusters, as the most stable morphology of larger clusters, going from 6 to 323 atoms, and compared with those on the extended slabs.
These calculations were carried out using the DFT Vienna Ab Initio Simulation Package (VASP) code, 36 with a spin-polarized formalism and the exchange-correlation functional developed by Perdew, Burke, and Ernzerhof (PBE). 37The core electrons up to and including the 3p levels of cobalt were kept frozen and their interaction with the valence electrons was described by the projector augmented wave (PAW) method. 38The DFT-D3 method with Becke-Johnson damping was used to include the long-range dispersion interactions. 39The kinetic energy cutoff of the wave functions was set to 400 eV.The Monkhorst-Pack k-point grid for the cluster calculations included only G-point.For the (0001) and (10% 11) surfaces, 3 Â 3 and 3 Â 2 supercells with 4 layers of atoms (top two relaxed, bottom two constrained in the bulk-relaxed positions) were used, with k-point meshes of 5 Â 5 Â 1 and 5 Â 4 Â 1, respectively.Structural optimisations were carried out without any constraints, with convergence criteria of 1.0 Â 10 À6 eV per atom and 1.0 Â 10 À2 eV Å À1 for the electronic energy and the atomic forces, respectively.The vacuum space in the unit cell was set to be 12 Å in the z-direction to avoid interactions between periodic images of the surface slabs, while clusters were positioned in the centre of a box with sizes adjusted such that there is at least 12 Å of vacuum between neighbouring clusters in each direction.The adsorption energy, E ads , was calculated as follows: where E surf+mol , E surf , and E mol are the DFT energies of the system with the adsorbed molecule, the clean surface or cluster, and the isolated molecule, respectively.

Choice of carboxylic acid
Earlier studies of the interactions of primary carboxylic and mercapto-carboxylic acids 40,41 with zeolites and ZnO surfaces have shown that adsorption modes and dielectric properties of the hybrid systems are rather insensitive to the chain length, since the interaction state is localised on the specific atom (either carboxyl oxygen or mercapto-sulphur).Structural and electronic properties of OA, as the experimental compound of interest, and VA, as a segment of OA cut at the fourth carbon atom of the alkyl chain where the distance from the carboxyl group is larger than three bonds thus diminishing the inductive effect, have been compared by DFT calculations using the computational details specified in Section 2.1.We anticipate that the carbon chain size does not affect significantly the properties of interest and the computational cost of the metadynamics can therefore be substantially reduced by using the shorter VA.

Metadynamics simulations
Ab initio metadynamics 42,43 of a single molecule of VA on a 57-atom cobalt cluster was performed using CP2K 44,45 with the PBE exchange-correlation functional and the Goedecker, Teter, and Hutter (GTH) Gaussian basis set. 46The whole system was placed in a cubic box of 25 Å a side and equilibrated in the NVT ensemble during 3 ps of Born-Oppenheimer molecular dynamics.

AN o s e
´-Hoover thermostat 47 was set to 500 K, in accordance with the temperature of experiments (160-250 1C), 27,34,35 and a time step of 1 fs was used for the integration of the equations of motion.Fluctuations in temperature, potential energy, and Co-O bonds during the equilibration are shown in Fig. 1.
The most recent experimental study conducted on the OA-cobalt nanoparticle system by Ansari et al. 35 suggested a monodentate-type interaction, which was hence taken as the starting structure.During the equilibration, we observed that the oxygen atom which was not bound to the nanoparticle surface attached to a neighbouring cobalt atom within less than 0.25 ps (Fig. 1, top), yielding a bidentate interaction with both oxygen atoms atop two neighbouring cobalt surface atoms.
After equilibration, metadynamics was performed.The two collective variables for the projection of the metadynamics free energy landscape were chosen to be the coordination numbers (CNs) of the two oxygen atoms of VA.Both Co-O distances during the equilibration period oscillated around 2.0 Å, as shown in Fig. 1, top.Accordingly, this value was used in the definition of CNs, together with exponents n = 8 and m = 14. 42aussian hills with height of 0.01 Hartree and width of 0.05 were deposited every 50 fs.The metadynamics simulation was carried out for 15 ps, after which one of the oxygen atoms bound to a new Co atom that was not included in the binding states observed thus far, hence reproducing previously formed interaction modes on a different site of the nanoparticle.This was considered to be a terminating event, as it indicates that the molecule has exhausted all the binding possibilities with the initially assigned pair of cobalt atoms.VMD was used for postprocessing. 48

Infrared spectroscopy
Infrared spectra of free VA and VA adsorbed on a cobalt cluster have been calculated within the perturbation theory at the PBE level.The PHONOPY code 49 was used to derive the associated intensities based on the Born effective charges tensor which yields information on the change in the dipole moment under atomic displacements.The intensity of a given mode can be evaluated as a square of all transition moments of the mode and expressed in terms of the Born effective charge matrix and eigenvectors of the mass-weighted Hessian.It is well-known that vibrational frequencies obtained by quantum-chemical calculations are typically larger than those evaluated experimentally, and thus, empirical scaling factors are often used for obtaining better match.However, since the current work focuses on the changes in the ligand spectra induced by binding to the nanoparticle rather than on the absolute values of the frequencies, the scaling was not carried out.

Adsorption behaviour with the nanoparticle size
The Co model clusters studied for the size-dependent adsorption consisted of 6, 13, 19, 30, 57, 76, 153, and 323 atoms arranged in an hcp-like atomic packing to mimic the behaviour of larger nanoparticles that are expected to perform well in cancer treatments.All clusters, except for the smallest 6 and 13-atom ones, exhibit facets of (0001) and (10% 11) orientations to minimise the overall surface energy.Size-dependent adsorption of O, HCOOH, and CH 3 COOH was therefore examined on both surfaces.
3.1.1Oxygen adsorption.Fig. 2, top, shows the optimised structures of clusters with a single adsorbed O atom and the corresponding adsorption energies, E ads , which are plotted in Fig. 3, top.For comparison, the structures of the periodic (0001) and (10% 11) surfaces are also shown (N = N).The oxygen atom, regardless of the cluster size or the surface stacking, adsorbs in the hollow position which allows binding to three different cobalt atoms of the (0001) surface.On the (10% 11) surface, clusters with more than 19 atoms can accommodate oxygen in the centre of four cobalt atoms, while the smaller surface area of the 19-atom cluster limits the interaction to three cobalt atoms.However, despite the similarities in the adsorption sites, there are still noticeable differences in the energetical aspects of the single oxygen atom adsorption on different clusters.On the (0001) surface, the interaction becomes stronger with shrinking cluster size but begins to weaken from 30-atom cluster downwards, with an exception for the adsorption on 13-atom cluster, which is the most favoured (E ads = À4.33 eV).A similar behaviour can be observed for the adsorption on the (10% 11) surface of the clusters, with lessening of E ads already triggered on the 76-atom cluster.The E ads values of the largest clusters considered, N = 153 and 323, differ only slightly from those of the periodic slab models.The initial trend of a strengthening interaction with decreasing cluster size prevails for larger particles only and reflects the decrease in the average CN of cobalt atoms as a consequence of the higher surface-tovolume ratio.Small clusters experience the influence of quantum effects, with a resulting non-scalable regime which often makes it hard to find a detailed quantitative dependence of the adsorption   energy on the cluster size.However, it is evident that even with a single atom adsorbate the two surfaces behave differently, which might have a more pronounced influence when it comes to the adsorption of larger molecular entities.In order to capture the main features of the cobalt cluster-oxygen atom systems, small 19-or 30-atom clusters may be used efficiently, although this may not be the case for carboxylic acids.
3.1.2HCOOH and CH 3 COOH adsorption.Fig. 2 illustrates the optimised structures and Fig. 3 the energy trends as a function of the cluster size for the adsorption of HCOOH and CH 3 COOH (middle and bottom, respectively).Compared to the single oxygen atom, where electronic effects make for the biggest change in the adsorption capabilities of differently sized clusters, the adsorption of acid molecules depends on the geometry of the clusters as well as on their electronic properties.On the (0001) surface, we observed dissociation only on the 6-atom cluster, where the preferred binding is with an oxygen atop one cobalt atom and the dissociated hydrogen atom in the bridging position over the two surface cobalt atom closest to the oxygen binding site.When the carboxyl oxygen was situated in the hollow as the most favourable adsorption site for a single oxygen atom adsorption, the interaction was weaker, probably due to the inconvenience of accommodating the -OH hydrogen atom.Adsorption of a carboxylic acid on the extended slab of the (10% 11) surface, on the other hand, shows a clear tendency toward dissociation.One of the oxygen atoms of the carboxyl group as well as the dissociated hydrogen tend to be located in the hollow sites on the largest clusters, while the remaining oxygen binds in the bridge position between two edge Co atoms on the 57-and 76-atom clusters.However, due to the limited (10% 11) surface areas of smaller clusters with 19 and 30 atoms, this dissociative behaviour is inhibited.It is important to note that dissociation was not manually forced and all initial structures of the acids prior to the adsorption process were undissociated.
For large clusters, the energy of HCOOH and CH 3 COOH adsorption on both the (0001) and (10 % 11) surface slowly progresses towards the asymptotic value for the infinite surface slabs.When adsorbed on the (0001) surface on clusters with less than 30 atoms, the adsorption energy reduces quickly with decreasing cluster size, yielding values well below the extended surface value.The opposite trends seen in large and small cobalt clusters meet at an intermediate size range of 15-50 atoms for the (0001) surface, with the 19-atom cluster exhibiting the least negative adsorption energies of À0.52 eV and À0.47 eV for HCOOH and CH 3 COOH, respectively.On the (10% 11) surface, the 57-atom cluster is the smallest of the tested sizes that showed the dissociative behaviour of an infinite slab model and, consequently, the adsorption energies of both carboxylic acids grow considerablywiththereductioninsizeofthe cluster.This trend is only interrupted with the 6-and 13-atom clusters, where the two surfaces cannot be distinguished.
As illustrated previously in the case of the oxygen atom, the strengthening of the adsorption as the cluster size shrinks is governed by the upward shift in the d-band centre throughout the whole range of clusters.However, for HCOOH and CH 3 COOH, this occurs only for large clusters, whereas the resemblance disappears once the clusters become smaller than 57 atoms, when the increase in the number of corner and edge sites with lower CNs leads towards stronger adsorptions.It can therefore be concluded that clusters with B50 atoms are precisely at the border where the domination of the electronic effects of large clusters vanishes and the ascendency of geometric effects of small clusters onsets.

Chain length dependence on the adsorption behaviour
To assess the suitability of VA as a model for OA when adsorbed on cobalt nanoparticles, both molecules were adsorbed on the 57-atom Co cluster and their structural, energetic, and electronic properties have been compared.

3.2.1
Structures and adsorption energies.Fig. 5 depicts the structures and adsorption energies of VA and OA adsorbed on the 57-atom cluster.Qualitatively, the positioning of the shorter chain of VA resembles the shape of the first few carbon atoms of OA.The adsorption energies on the ð0001Þ surface of the cluster differ only by 0.03 eV, while undissociative and dissociative modes on the (10% 11) surface give adsorption energies of OA lower by 0.10 and 0.17 eV than VA, respectively.3.2.2Electronic properties.For comparison of the electronic properties, the densities of states (DOS) of VA and OA adsorbed on the (0001) and (10% 11) surfaces of the 57-atom cluster are presented in Fig. 6, together with their projection on the carbon chains and the carboxyl oxygen atoms.In all three cases of adsorption, VA and OA show similar electronic structures, with bonding hybridisation mainly occurring at BÀ5.0 eV for undissociative adsorption and at BÀ3.5 eV for dissociative adsorption.
To estimate the share of carbon-chain charge transfer in the overall bonding process, the Bader charge difference 51 was calculated as where q ads and q vac are the number of electrons of all the atoms in the carbon-chain of the adsorbate on the cluster and in vacuum, respectively.Dq c represents the charge gained or lost in the process of adsorption.If it is positive, the electrons are lost from the carbon-chain, and vice versa, if it is negative, the carbon-chain drains electrons from the atoms of the cluster.In thecaseofadsorptiononthe(0001)surface,thecarbonchainofVA, namely -(CH 2 ) 3 CH 3 , experiences a charge transfer of Dq c = +0.062eV, which increases to +0.086 eV for the -(CH 2 ) 7 CQC(CH 2 ) 7 CH 3 carbonchain of OA.For undissociative adsorption on the (10% 11) surface, carbon-chain charge transfers are equal to +0.060 and +0.057 eV, while for dissociative adsorption they are equal to À0.011 and À0.008 eV for VA and OA, respectively.No appreciable electron transfer occurs during the adsorption, with carbon chains of both acid molecules losing or, in the case of dissociation, gaining a negligible amount of charge during the process.Importantly, the differences between the two acids are minimal.Combined with the findings on the energetics of adsorption, the comparison of the electronic properties suggests that the carboxyl group plays the dominant role and interaction with the surface is not significantly affected by the size of the carbon-chain.Accordingly, VA has been used in the continuation of the study.

Free energy landscape of valeric acid on a Co cluster
Now that we have identified a system that can accurately represent the OA-cobalt nanoparticle interaction at affordable computational cost, the stability and the likelihood of the possible binding modes can be explored.To provide insight in the mechanisms and barriers of the transitions while taking into account temperature effects, we have employed metadynamics, which enforces the process of sampling between different minima by adding a small repulsive Gaussian potential in a steady manner.
The most prominent differences between the well-known experimentally defined modes of interaction are the CNs of the two oxygen atoms, O a and O b , of carboxyl group with respect to surface metal atoms.By going from a bidentate mode, for which both CN (O a -Co) and CN (O b -Co) equal one (1,1), to a monodentate, the CN of one oxygen decreases to zero, while the other oxygen can bind to one, two, or three cobalt atoms, yielding combinations of (0,1), (0,2), and (0,3), respectively.Transition to a chelate interaction induces binding of both oxygen atoms to the same cobalt atom, which would at the first sight give exactly the same combination of CNs of bidentate binding, namely (1,1).However, the oxygen-oxygen repulsion results in an increased distance between the oxygen atoms in the chelate configuration, making the chelating Co-O bonds longer than in the bidentate, giving an approximate CN combination of (0.5,0.5).Thus, being able to distinguish between the most relevant structures, CN (O a -Co) and CN (O b -Co) make a good choice for collective variables.
Co-O distances and CNs of O a and O b have been presented for the whole metadynamics trajectory in Fig. 7, alongside the most representative structures.The starting structure was bridging bidentatewitheachoxygenbondedtooneCoatom,ascontinued from the ab initio molecular dynamics equilibration.After 2.5 ps of metadynamics run, O a formed a bridge in between neighbouring cobalt atoms, resulting in a chelating form which was stable for a  The type of interaction between the carboxylate head and the metal atoms can be confirmed by the wavenumber separation (D) between the asymmetric (n as (COO À )) and symmetric (n s (COO À )) stretches of COO À : 34,52,53 the largest D (B200-320 cm À1 ) corresponds to the monodentate interaction, medium D (B140-190 cm À1 ) to bridging bidentate interaction, and the smallest D (Bo110 cm À1 ) to chelate interaction.We have selected three geometries corresponding to the main binding modes resulting from the metadynamics trajectory and simulated their IR spectra, together with that of VA in vacuum (Fig. 8).The IR spectrum of the pure acid shows good agreement with the positions of the fundamental frequencies of carboxylic acids.After adsorption, the intense peak at 1749 cm À1 , typical of the CQO stretch, disappears and two new peaks as features of carboxylate are formed at B1200 and 1650 cm À1 instead.The obtained wavenumber separations D between the asymmetric and symmetric stretches are 320, 108, and 42 cm À1 , respectively, which, following the experimental trends discussed above, confirms our assignment of the type of interactions to monodentate, bridging bidentate,  The next question that arises is to what extent does the change from bidentate to monodentate interaction happen; i.e., what are the pathways and barriers taking place, considering finite-temperature effects at the appropriate thermodynamic conditions.To investigate the underlying mechanisms in pertinent detail, the reconstructed free energy surface projected onto the plane spanning the two collective variables has been plotted in Fig. 9. Several minima beyond the one that defines the starting structure are readily identified and denoted using Roman numbers.Minimum I represents the bridging bidentate and is by far the most stable structure with an energy of À139 kJ mol À1 .The minima denoted as IIa and IIb are atop monodentate structures with O a and O b attached to one surface cobalt atom, respectively.They are higher in energy with respect to bridging bidentate and sit at around À115 kJ mol À1 .Chelates are represented by IVa and IVb and are energetically less favourable at À90 kJ mol À1 .Bridging monodentate, III, is least likely to be formed with an energy of around À60 kJ mol À1 .A wide area in energy up to À70 kJ mol À1 around the global minimum can be assigned to exhaustive efforts of the system to stay in the bidentate mode with both oxygen atoms of carboxyl group bonded to a single cobalt atom, even when one or both bonds do not reach the ideal value of 2.0 Å.What clearly emerges from Fig. 9 is that the bridging bidentate is the most favourable mode of interaction between carboxylic acids and cobalt nanoparticles.barrier of 77 kJ mol À1 .However, the inclusion of the chelate interaction in the transition separates the process into two major steps (B and C), for which the two barriers of 58 and 41 kJ mol À1 are both lower than in the direct route D.
Other possible pathways include bridging monodentate as a starting point, with two possible final destinations -atop monodentate (A, B), or the global minimum -bridging bidentate (A, C).Starting with a bridging monodentate and going towards the atop analogue requires an energy of 24 kJ mol À1 to overcome the rate-determining step, which increases to 41 kJ mol À1 if bridging bidentate is the final structure.In general, all transitions between minima are moderately energetically costly, with the rate-determining steps depending on the pathway.

Conclusion
We have presented a theoretical investigation which clarifies the interaction modes of cobalt nanoparticles with carboxylic acids.On the basis of DFT calculations, we suggest that a 57-atom cluster is the smallest cluster that can be adopted to model the interaction of carboxylic acids with Co nanoparticles.In addition, our ab initio metadynamics simulations have provided important insights into the free energy landscape around the two experimentally observed binding modes, namely monodentate and bridging bidentate.In particular, we found that the free energies of the two modes differ by 24 kJ mol À1 in favour of the bridging bidentate, with the barrier allowing for the direct transition being moderately high, i.e. 77 kJ mol À1 .A chelating interaction was detected as a possible intermediate state, in a transition mechanism with an energy barrier of 58 kJ mol À1 for the rate-limiting step.We also observed the formation of a fourth possible mode, bridging monodentate, but its high free energy makes it even more unlikely to be present under experimental conditions.Thus, it is conceivable that the majority of the carboxylic acid molecules attached to cobalt nanoparticles would maintain a bridging bidentate interaction on the surface under synthesis conditions.Determination of the binding mode between the two modalities was a first step towards simulating the organic acid coatings on cobalt nanoparticles and their behaviour in an in vivo environment.The next set of issues to be addressed are the possible density of the adsorbed molecules and whether the introduction of the coating would be sufficient to prevent further chemical reactions, such as oxidation, to take place at the nanoparticles' surface.Work along these directions is in progress.
Pl e a s e n o t e: C h a n g e s m a d e a s a r e s ul t of p u blis hi n g p r o c e s s e s s u c h a s c o py-e di ti n g, fo r m a t ti n g a n d p a g e n u m b e r s m a y n o t b e r efl e c t e d in t his ve r sio n.Fo r t h e d efi nitiv e ve r sio n of t hi s p u blic a tio n, pl e a s e r ef e r t o t h e p u blis h e d s o u r c e.You a r e a d vis e d t o c o n s ul t t h e p u blis h e r's v e r sio n if yo u wi s h t o cit e t hi s p a p er. Thi s v e r sio n is b ei n g m a d e a v ail a bl e in a c c o r d a n c e wit h p u blis h e r p olici e s.

Fig. 1
Fig. 1 Co-O distances (top), potential energy (middle), and temperature (bottom, with horizontal line representing the temperature of the thermostat, T = 500 K) during the ab initio molecular dynamics equilibration of valeric acid on the 57-atom Co cluster.

Fig. 4
Fig. 4 shows the d-band energies calculated following the definition by Hammer et al., 50 for both the 3d spin-up and spindown electrons of cobalt atoms in differently sized clusters and the infinite surfaces.The spin-up (3 unpaired electrons per Co atom which are responsible for oxygen-binding) d-band energy sequence mirrors the trend of oxygen adsorption, with the 13-atom cluster having the most positive d-band energy and, accordingly, the strongest interaction.Small clusters experience the influence of quantum effects, with a resulting non-scalable regime which often makes it hard to find a detailed quantitative dependence of the adsorption

Fig. 3
Fig. 3 Calculated adsorption energies, E ads , of O (top), HCOOH (middle), and CH 3 COOH (bottom) adsorption on the (0001) and (10 % 11) surfaces of Co N clusters as a function of the cluster size, N. Results for surface slabs are shown as N = N. Lines are to guide an eye only.

Fig. 4
Fig. 4 d-Band energies for 3d spin-up and spin-down electrons of Co as a function of cluster size, N. Results for surface slabs are shown as N = N.

Fig. 5
Fig. 5 Structures of valeric (left panel) and oleic acid (right panel) undissociated on the (0001) surface (left) and both undissociated (middle) and dissociated (right) on the (10 % 11) surface of the 57-atom Co cluster.Adsorption energies are also shown below each structure.

Fig. 7
Fig. 7 Representative structures sampled during the metadynamics run (top).Changes in the coordination numbers CN (middle), and Co-O bond distances d Co-O (bottom) throughout the metadynamics trajectory.O a and O b CNs and distances, as well as windows across the simulation are coloured according to the colours associated in the top panel.

Fig. 8
Fig. 8 Simulated IR spectra of valeric acid in vacuum (top) and adsorbed on the 57-atom Co nanoparticle (bottom) as monodentate, bridging bidentate, and chelate, with accompanying wavenumber separations (D) between the asymmetric and symmetric stretches of COO À .Both axes are the same for all spectra.

Fig. 9
also illustrates that the transition of the system from the minima II to the global minimum I has two possible pathways, i.e. one with a chelating mode as an intermediate (B, C) and another direct pathway (D).The formation of the intermediate chelating state IV allows the molecule to access a lower energy pathway.In particular, the shift from atop monodentate directly to the bridging bidentate global minimum following route D, requires the system to overcome an energy

Fig. 9
Fig. 9 Free energy surface for the adsorption process of valeric acid on the 57-atom Co cluster; Roman numbers denote different minima, with accompanying letters a and b standing for the configurations differing only by the oxygen atom, O a or O b , to which the surface cobalt atom is bonded.All structures are visually represented in the middle panel.Upper case letters illustrate the lowest energy paths whose mechanisms are proposed together with the free energy barriers for individual processes in the bottom panel.All energy barriers are given in kJ mol À1 .
an intermediate state.No desorption of VA was observed on the time scale of the entire simulation.After 15 ps, instead of desorption of the molecule, O b bonded to the third different surface cobalt atom.Ultimately, four distinct interaction modes were captured: atop monodentate where one of carboxyl oxygen atoms binds one cobalt atom; bridging monodentate where one O atom has two bonds with a pair of neighbouring Co atoms; bridging bidentate where each carboxyl oxygen binds a single Co atom from a neighbouring pair; and, finally, chelate, where one O atom forms two bonds, one with each atom from Co pair, while the other O atom forms a single bond with one of the Co atoms from the same pair.The bridging bidentate mode was formed for more than half of the time, B54%, with the bridging monodentate, atop monodentate, and chelate mode having a similar share of B18, B15 and B13%, respectively.