13.07.2015 Views

Documento PDF - UniCA Eprints - Università degli studi di Cagliari.

Documento PDF - UniCA Eprints - Università degli studi di Cagliari.

Documento PDF - UniCA Eprints - Università degli studi di Cagliari.

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

1.3 theoretical modeling of hybrid interfaces 13acting atoms representing the material of interest by solvingthe Newton’s Equation of motions (F = ma). Forcesare derived from a suitable model potential of the atomicpositions that is calibrated in such a way to reproduce a setof physical properties of the material (see Appen<strong>di</strong>x A).The relatively low computational workload associated toMPMD, allows to obtain pre<strong>di</strong>ctive informations regar<strong>di</strong>ngthermo<strong>di</strong>namics and microcrystalline evolution over the 10ns timescale of systems as large as 10 nm.Furthermore, the molecular dynamics approach makespossible to easily take into account long range <strong>di</strong>spersiveinteractions, by using simple Lennard-Jones type potential.The accurate description of interatomic forces in hybridsis however challenging. A general model potential for thehybrid system is not available, but there are reliable potentialsfor the organic and inorganic phases separately. Organicpolymers can be described by means of the “triedand true” Amber force field [62], while in the case of metaloxide the modeling is slightly more complicated. In particular,the ZnO description must take into account its partiallyionic and partially covalent nature [63]. A simpleand succesfull solution is the use of pair interactions consistingof a short-range part (usually a Buckingham interaction[64]) and long-range Coulombic terms employingfixed charges. This method, however, does not take into accountthe charge re<strong>di</strong>stribution around a defect or at thesurface [63]. In the more advanced shell-model description[63, 65, 66], the electronic polarizability is includedad<strong>di</strong>ng an ad<strong>di</strong>tional charged site to each ion connected viaa spring [67]. The shell models however do not properly describethe covalent character of ZnO [63], problem that canbe aided by using higher-order terms in the many-bodyexpansion [68], or by neglecting also the ionic characterand using a bond-order potential [69], but these solutionshave also several drawbacks [63], inclu<strong>di</strong>ng larger computationalcosts. Finally, the reactive force field (ReaxFF)[63, 70, 71] is also a bond-order interaction model consistingof the two-body, three-body and four-body short-rangeinteraction terms. It allows the re<strong>di</strong>stribution of charges,can simulate the breaking and reforming of bonds and canreproduce the structures and mechanical properties of con-

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!