Modelling and Simulation in Materials Science and Engineering PAPER A 3D platform for the morphology modulation of materials: first principles calculations on the thermodynamic stability and surface structure of metal oxides: Co3O4, α-Fe2O3, and In2O3 To cite this article: M M Ferrer et al 2016 Modelling Simul. Mater. Sci. Eng. 24 025007   View the article online for updates and enhancements. Related content Effects of surface stability on the morphological transformation of metals and metal oxides as investigated by first- principles calculations Juan Andrés, Lourdes Gracia, Amanda Fernandes Gouveia et al. - Identifying and rationalizing the morphological, structural, and optical properties of -Ag2MoO4 microcrystals, and the formation process of Ag nanoparticles on their surfaces: combining experimental data and first-principles calculations Maria T Fabbro, Carla Saliby, Larissa R Rios et al. - Modelling the effect of particle shape on the phase stability of ZrO2 nanoparticles Amanda S Barnard, Rakesh R Yeredla and Huifang Xu - Recent citations Surface broken bonds: An efficient way to assess the surface behaviour of fluorite Zhiyong Gao et al - Inverse Wulff Construction for Surface Energies of Coexisting and Missing Surfaces of Crystal Particles Fuming Lai et al - Preparation of -CaSO 4 ·½H 2 O with tunable morphology from flue gas desulphurization gypsum using malic acid as modifier: A theoretical and experimental study Qingjun Guan et al - This content was downloaded from IP address 186.217.236.64 on 29/05/2019 at 18:47 https://doi.org/10.1088/0965-0393/24/2/025007 http://iopscience.iop.org/article/10.1088/0957-4484/26/40/405703 http://iopscience.iop.org/article/10.1088/0957-4484/26/40/405703 http://iopscience.iop.org/article/10.1088/0957-4484/26/40/405703 http://iopscience.iop.org/article/10.1088/0957-4484/26/40/405703 http://iopscience.iop.org/article/10.1088/1468-6996/16/6/065002 http://iopscience.iop.org/article/10.1088/1468-6996/16/6/065002 http://iopscience.iop.org/article/10.1088/1468-6996/16/6/065002 http://iopscience.iop.org/article/10.1088/1468-6996/16/6/065002 http://iopscience.iop.org/article/10.1088/1468-6996/16/6/065002 http://iopscience.iop.org/article/10.1088/1468-6996/16/6/065002 http://iopscience.iop.org/article/10.1088/0957-4484/17/12/038 http://iopscience.iop.org/article/10.1088/0957-4484/17/12/038 http://dx.doi.org/10.1016/j.mineng.2018.09.024 http://dx.doi.org/10.1016/j.mineng.2018.09.024 http://dx.doi.org/10.1016/j.jcrysgro.2018.12.006 http://dx.doi.org/10.1016/j.jcrysgro.2018.12.006 http://dx.doi.org/10.1016/j.jcrysgro.2018.12.006 http://dx.doi.org/10.1016/j.jcis.2018.06.068 http://dx.doi.org/10.1016/j.jcis.2018.06.068 http://dx.doi.org/10.1016/j.jcis.2018.06.068 http://dx.doi.org/10.1016/j.jcis.2018.06.068 http://dx.doi.org/10.1016/j.jcis.2018.06.068 https://oasc-eu1.247realmedia.com/5c/iopscience.iop.org/181771891/Middle/IOPP/IOPs-Mid-MSMSE-pdf/IOPs-Mid-MSMSE-pdf.jpg/1? 1 Modelling and Simulation in Materials Science and Engineering A 3D platform for the morphology modulation of materials: first principles calculations on the thermodynamic stability and surface structure of metal oxides: Co3O4, α-Fe2O3, and In2O3 M M Ferrer1, A F Gouveia1, L Gracia2, E Longo3 and J Andrés2 1 INCTMN-UFSCar, Universidade Federal de São Carlos, PO Box 676, 13565–905 São Carlos, SP, Brazil 2 Department of Analytical and Physical Chemistry, University Jaume I (UJI), Castelló 12071, Spain 3 INCTMN-UNESP, Universidade Estadual Paulista, PO Box 355, CEP 14801–907 Araraquara, SP, Brazil E-mail: andres@qfa.uji.es Received 16 July 2015, revised 16 November 2015 Accepted for publication 26 November 2015 Published 12 January 2016 Abstract Essentially, the exposed crystal planes of a given material, which primarily determine their morphology, tremendously affect its behavior. First principle calculations, based on the Wulff construction model and broken bonding density index, have been performed to calculate the equilibrium and their transformations for different metal oxides: Co3O4, α-Fe2O3, and In2O3. Present results point out that starting by surface thermodynamics is a helpful approach to predict and assess the morphology transformations of these materials. These complete set of morphologies may serve as a guide for researchers, when analyzing the images from electron microscopies, to gain further understanding of how to control crystal shape synthetically by tuning the surface chemistry and by controlling the relative values of surface energies. Keywords: morphology, Wulff construction, surface energy, Co3O4, α-Fe2O3, In2O3 (Some figures may appear in colour only in the online journal) M M Ferrer et al Printed in the UK 025007 MSMEEU © 2016 IOP Publishing Ltd 2016 24 Modelling Simul. Mater. Sci. Eng. MSMS 0965-0393 10.1088/0965-0393/24/2/025007 Paper 2 1 10 Modelling and Simulation in Materials Science and Engineering IOP 0965-0393/16/025007+10$33.00 © 2016 IOP Publishing Ltd Printed in the UK Modelling Simul. Mater. Sci. Eng. 24 (2016) 025007 (10pp) doi:10.1088/0965-0393/24/2/025007 mailto:: andres@qfa.uji.es http://crossmark.crossref.org/dialog/?doi=10.1088/0965-0393/24/2/025007&domain=pdf&date_stamp=2016-01-12 publisher-id doi http://dx.doi.org/10.1088/0965-0393/24/2/025007 2 1. Introduction A few years ago, McIlroy et al [1] stated: ‘Currently, there is a large divide between surface scientists and most catalysis researchers (who make and test new materials), caused by the inability of existing experimental and theoretical techniques to deal with the real-world nano- morphologies. Closing this divide would allow the entire surface science community to be much more effective.’ In recent years, there has been such speed and development in this area that this is no longer a valid assessment. The current combination of theory, simulation, and experiment is not only capable of developing new functional materials with predictable and desirable properties, but also enables us to understand the nature of micro and nano-structures with novel morphologies. Control of the morphology of micro and nanocrystals depends on factors such as the internal structure of the crystal and the external growth conditions. The surface facet ener- gies can be modified by the use of different synthetic methods, reaction times, the nature of surfactants etc. In this sense, the rational design of nanomaterials is one of the holy grails of nanoscience and nanotechnology. The precise regulation and controlled synthesis of specific morphologies is crucial to enhance their performance in practical applications of micro and nanomaterials [2]. Predicting the 3D morphologies of crystalline materials can be used to guide the analy- sis of field-emission scanning electron microscopy images of the materials, and could be a very useful tool to design the synthesis of new materials. However, there are many scientific and technological challenges to overcome; in particular, we must understand the relationship between properties and structure. Since various physical and chemical processes occur at surfaces, an understanding of the surface structure and the factors that affect surfaces is important. The shape of crystals grown under or near equilibrium conditions is governed mainly by surface energies, form- ing simple shapes with well-defined facets that can achieve minimum surface energy [3]. One method of predicting the morphology of a crystal is to calculate the energies of the different surfaces and these value can be associated to the fraction of surface atoms with unsaturated coordination shells, i.e. the number of under-coordinated surface atoms. The predominant surface will be the most stable surface; but measuring the geometries and energies of the surface facets is extremely challenging. In this context, surface energies are not readily acquired by experiment, and computer modeling and simulation are necessary to obtain these energies, being a powerful tool to explore the morphological mechanism at the atomistic/molecular level. Theory and simulation in this research area should be further strengthened, which focuses on understanding the morphology, can guide the efficient synthesis of outstanding materials. First principles investigations are well established in the study of the morphol- ogy and surface properties of materials, from simple to complex metal oxides [4–27]. There are several insightful reviews that report experimental results and theoretical calculations on the most important aspects that govern crystal shape modulation in semiconductors and in metal nanoparticles [12, 28–31]. Theoretical works, in which surface modification by explicit or implicit solvation/hydration/redox effect are taken into account, can be consid- ered the current state-of-the-art on theoretical methods and computing procedures. In this context, good illustrations have been published, as showed in these recent works: (1) the stabilization effect of surface impurities on the structure of ultrasmall ZrO2 nanoparticles [32], (2) on the basis of theoretical predictions, Yang et al have synthesized uniform anatase TiO2 single crystals with a high percentage of reactive facets [33], (3) theoretical investiga- tions into the structure and surface thermodynamics of the (1 0 0) facet of Co3O4 at various M M Ferrer et alModelling Simul. Mater. Sci. Eng. 24 (2016) 025007 3 redox conditions [34] as well as the morphological control of Co3O4 to enhace the oxida- tion of CO at low tempertaure [35], and (4) the mechanism of morphology transformation of tetragonal phase LaVO4 nanocrystals controlled by surface chemistry: experimental and theoretical insights, in which equilibrium geometries, surface energies, and surface tensions are calculated for selected low-index surfaces of tetragonal LaVO4 under the different sur- face passivated conditions [24]. Recently, we have determined the morphologies of several metal oxides such as SnO2, PbMoO4, and CaWO4 using ab initio quantum mechanical calculations [15–20]. First prin- ciples calculations were also used to study the the effects of surface stability on the morphological transformation of different metals and metal oxides (Ag, anatase TiO2, BaZrO3, and α-Ag2WO4) [36]. Morphological transformations are due to the geometric constraints imposed by the crystal structure and are associated with the relative values of the surface energy of each face. In this context, and we have analyzed theoretical morphology modu- lation for α-Ag2WO4 [37], Ag3PO4 [38] and Ag2MoO4 [39] materials in order to explain the surface energy changes observed in experimental conditions, and how the knowledge of surface-specific properties can be utilized to design a number of crystal morphologies that may offer improved performance in various applications. It is important to remark that along our calculations specific solvent and/or surfactants interactions or hydration/redox effects on each face are not considered. The main objective of this work is to show how the values of surface energies, obtained from first principle calculations in vacuo, can be used to obtain possible morphologies under thermodynamic equilibrium conditions for three metal oxides: Co3O4, α-Fe2O3, and In2O3 as well as to better understand the mechanism of morphological transformation at atomic level. The results of calculations have been compared with the experimental data. It is important to note the following aspects: (i) we have used the values of surface energy from literature, and we are confident that the proposed set of morphologies is not dependent on the calculation level; (ii) the values of surface energies are obtained in vacuo; (iii) from these values we are capable to predict the available morphologies for a crystalline material; (iv) by screening of the relative values of surface energies, a desired morphology can be obtained; (v) an analysis of how the idealized in vacuo surface results can be related to the in situ behavior, and this strategy allow us to find the target and control of crystal morphologies, as well as to rationalize the different channels starting from the ideal morphology. The whole paper is arranged as follows: (2) a description of theoretical method and com- putational procedure, (3) the application of the method and the presentation of results, and (4) our main conclusions and perspectives close this letter. 2. Theoretical method and computational procedure From the thermodynamic point of view, the free energy of the crystal facets (Esurf) determines the equilibrium shape by a classic Wulff construction, which minimizes the total surface free energy at a fixed volume [40]. The Wulff theorem provides a simple relationship between the Esurf of the (h k l) plane and the distance, rhkl, in the normal direction from the center of the crystallite. In addition, we have calculated the broken bonding density (Db) index proposed by Gao et al [41] Db, the number of bonds broken per area on creation of a new surface calculated from equation (1), where Nb is the number of broken bonds per unit cell area on a particular surface and A is the area of the surface unit cell. =D N A/b b (1) M M Ferrer et alModelling Simul. Mater. Sci. Eng. 24 (2016) 025007 4 3. Results and discussion 3.1. Co3O4 Co3O4 is a multifunctional material that has many applications and has received much theoret- ical and practical investigation [14, 42–45]. Co3O4 has a face centered cubic spinel structure, crystallizing in space group Fd-3m, and has two formula units (Z = 2) in the unit cell. Zasada et al [14] obtained the surface energies for the low index faces (1 0 0), (1 1 0), and (1 1 1) using calculations on a Co3O4 slab and compared the Wulff construction for this simulated material with the experimental morphology of a Co3O4 nanocrystal. Su et al [46] synthesized single crystalline Co3O4 nanocrystals with different exposed crystal planes and wide range of mor- phologies such as cubic, pseudo octahedral, nanosheets, and nanoplatelets, and so on. These authors found a correlation between the reducing charge − discharge overpotential with crys- tal planes of Co3O4 and established an order of catalytic activity: (1 1 1) > (1 1 0) > (1 1 2) > (1 0 0), and very recently, Liu et al [47] have confirmed the facet-dependent electrocatalytic performance of Co3O4 for rechargeable Li − O2 battery. The Wulff construction for Co3O4 structure has been obtained from the work of Su et al [46]. Table 1 lists the values of the surface energies obtained by Su et al [46] as well as the broken bonding density (Db) of the faces used in the Wulff construction. The Db values can be directly related to the order of surface stability because higher values represent the presence of a larger quantity of defects on the surface (broken bonds). For the Co3O4 simulated system, table 1 reveals that the order of stability of the faces is (1 0 0) > (1 1 0) > (1 1 2) > (1 1 1). The complete set of morphologies for Co3O4 structure has been depiected in figure  1, in which the transformations can be obtained by changing the ratio bettwen the values of Esurf for each surface. Transformations between the different morphologies are due to the geometric constraints imposed by the crystal structure and are associated with the relative surface energy value of each surface. This interpretation has the advantage that all faces grow from the initial Co3O4 crystal (ideal) as function of their surface energy value. It is important to note that these images will serve as a guide to better interpret the morphological changes according to the synthetic method used, and so on. In addition, the different degrees of saturation at the surface, and consequently the different Db values, in the simulated system play an important role in the predicted morphology of the crystal and these differences must be considered when compar- ing experimental and theoretical results. In a real system, the surfaces can adsorb species pre- sent during synthesis, and these adsorbed species will affect the surface energy and, thus, the final morphology, as shown in figure 1. Therefore, several research groups have been reported the synthesis of different morphologies of the Co3O4 indicated on figure 1, such as cubic (a) [46–49], octahedron (b) [47], and pseudo octahedron (c) [46], in which it is necessary to differ one of the (1 1 0)-surface group ( ( )E * surf 1 1 0 can be ( )Esurf 1 1 0 or ( )Esurf 0 1 1 or ( )Esurf 1 0 1 ). Table 1. Values of Esurf, number of broken bonds, area, and broken bonding density (Db) calculated for Co3O4. Co3O4 Surface Esurf (J m−2) Broken bonds Area (nm2) Db (nm−2) (1 0 0) 0.92 12 0.3258 36.83 (1 1 0) 1.31 20 0.4607 43.41 (1 1 1) 2.31 15 0.2821 53.17 (1 1 2) 1.46 36 0.7980 45.11 M M Ferrer et alModelling Simul. Mater. Sci. Eng. 24 (2016) 025007 5 3.2. α-Fe2O3 Iron nanomaterials are used in catalysts and electrodes, as medical imaging contrast agents, adsorbents, and magnetic recording devices. The morphology of individual particles is one of the determining factors in their performance, stability, and in the safety of both engineered nanoparticles as well as natural samples. For example, pseudocube-shaped hematite nanopar- ticles have a significantly larger permanent magnetization and coercivity than rice-shaped particles [50]. Hematite (α-Fe2O3) also has unique optical and electronic properties and is among the most important materials discussed in literature regarding energy-conversion processes [51–55]. Thus, there is considerable scientific and technological interest in the morphology of iron oxides, as it is demonstrated in various studies, reporting different morphologies synthe- sized via a template-free hydrothermal route [56], under the influence of −H PO2 4 ions [57], pre- pared by forced hydrolysis of ferric chloride [58], synthesized via an octadecylamine-assisted [59] and via atriphenyl-phosphine-assisted [60] hydrothermal method and with an oriented attachment mechanism controlled by different solvents [61]. Recently, a relationship between surface structure of hematite (0 0 1) and the macroscopic charging (zetapotential measure- ments in electrolyte solutions as a function of pH) has been reported [62], and the shape of Fe3O4 ferroferric oxide nanoparticles synthesized by a coprecipitation method as well as their electrochemical performances has been investigated [63]. Much effort is being directed towards the controlled synthesis of functional nanomaterials in which the shape and size can be engineered. α-Fe2O3 crystallizes in the trigonal space group R-3c, and is isostructural with corundum, Al2O3. It is one of the most thermodynamically stable and abundant phases among all of the iron oxides and oxyhydroxides. Each rhombohedral unit cell contains four Fe atoms, dis- tributed over two cation layers. The spins of the Fe atoms in the layers perpendicular to the [0 0 0 1] axis are aligned. The spin of the alternating layers is antiparallel, which leads to an antiferromagnetic material. Guo et al [64] performed surface energy calculations for α-Fe2O3. The authors stud- ied the most important surfaces of hematite: (0 0 1)-Fe, (0 1 2), (1 0 0) and (1 1 0), of which Figure 1. Crystallographic structure and morphologies of Co3O4 with crystal planes (1 0 0), (1 1 0), (1 1 1) and (1 1 2). Surface energy is in joule per⋅square meter. M M Ferrer et alModelling Simul. Mater. Sci. Eng. 24 (2016) 025007 6 (0 0 1) was the most extensively studied. The reported surface energies for these slabs and the morphologies generated by Wulff construction are shown in the center of figure 2. Central morphologies of figure 2 seem to be in agreement with scanning electron microscopy images of α-Fe2O3 particles formed in 5 mol% octadecylamine (ODA) for 48 h reported by Liu et al [59]. However, they reported a tetrakaidecahedral morphology increasing ODA to 12 mol% and enhancing the appearance of (1 1 2), (0 0 1) and (−1 1 2) facets over the other. The influ- ence in morphology of temperature as well as a mixture EtOH/H2O has been investigated [56] obtaining rhombohedra more or less truncated and cubic morphologies, which fit with some morphologies depicted in figure 2 involving different ratios of (0 1 2) and (1 0 0) surfaces. Table 2 shows that the value for Db index decreases in the order: (0 0 1)/Fe > (1 0 0) > (1 0 1) > (0 0 1)/O > (0 1 2) > (1 1 0), while the surface energy decreases in this order: (1 0 0) > (1 0 1) > (1 1 0) > (0 0 1)/Fe > (0 1 2) > (0 0 1)/O. Because H atoms were added to the (0 0 1)-O surface there is a large difference in the order of the (0 0 1)-O surface energy compared to the Db representation. For this reason, the (0 0 1)-O surface energy was not used in the Wulff construction to generate the hematite morphology shown in figure 2. From the analysis of the (0 0 1)-O surface, it is clear that H atoms stabilized the under-coordinated O atoms on the surface, and this decreased both the surface energy and the density of broken bonds. In this system, the Db values are not consistent with the simulated surface energy order. This can be related to the theoretical methodology used by Guo et al [64] to study the surface energy. That is, the surface model contained surface vacancies that are not fully explained. 3.3. In2O3 Indium oxide (In2O3) is an n-type highly degenerate wide band gap material that has many applications. Indium oxide is used in flat panel displays, photoelectric devices, windshield defrosting layers, heat-reflecting mirrors, high-transparency layers and gas sensor devices Figure 2. Crystallographic structure and morphologies of α-Fe2O3 with crystal planes (0 0 1), (0 1 2), (1 0 0), (1 0 1) and (1 1 0). Surface energy is in joule per⋅square meter. M M Ferrer et alModelling Simul. Mater. Sci. Eng. 24 (2016) 025007 7 [65]. The body-centered cubic polymorph of In2O3 crystallizes in the space group Ia-3 with eight formula units (Z = 8) in the unit cell. Walsh and Catlow [66] performed surface energy calculations for In2O3. We have used their surface values to generate the Wulff construction and the morphologies of In2O3. Figure 3 shows the energies, morphologies and the comparison with experimental results for In2O3 in our laboratory [67]. According to table 3, the order of stability using the Db method is (1 1 1) < (1 1 0) < (1 0 0), the same as the theoretical surface energies. Thus, In2O3 surfaces with higher energies possess the highest amount of unsaturation on the surface. 4. Conclusions It is mandatory to control the morphology of a micro/nanocrystal for for many scientific and technological applications. In this work, we used first principle calculations, based on the Wulff construction model and broken bonding density index, to investigate the effect of surface chemistry on the mechanism of morphology transformation of metal oxides: Co3O4, α-Fe2O3, Table 2. Values of Esurf, number of broken bonds, area, and broken bonding density (Db) calculated for α-Fe2O3. α-Fe2O3 Surface Esurf (J m−2) Broken bonds Area (nm2) Db (nm−2) (0 0 1)-Fe 1.145 12 0.2196 54.65 (0 1 2) 1.056 8 0.2733 29.27 (1 0 0) 1.369 28 0.6923 40.44 (1 0 1) 1.306 12 0.2421 49.57 (1 1 0) 1.230 8 0.3997 20.01 (0 0 1)-O 0.782 9 0.2196 40.99 Figure 3. Crystallographic structure and morphologies of In2O3 with crystal planes (1 0 0), (1 1 0) and (1 1 1). Surface energy is in joule per⋅square meter. M M Ferrer et alModelling Simul. Mater. Sci. Eng. 24 (2016) 025007 8 and In2O3. Using a methodology recently developed by our research group, a detailed con- figurational analysis of the different faces was performed. A correlation between the broken bonding density (Db) and the surface energies was found. The relaxed structures and surface energies were used to obtain a complete array of accessible morphologies. Manipulating the values of surface energy for each face plays a significant role, and it is exemplified in the aforementioned materials. This method provides an approach with both predictive and explanatory capabilities and the results presented here are an illustration of how first principle calculations can rationalize the mechanisms stabilizing the morphology of micro- and nano-particles at the atomic level. The calculated diagrams relate the crystal growth conditions with the observed morphologies in an attempt to rationalize the morphologies obtained under different conditions. The mor- phological transformation of Co3O4, α-Fe2O3, and In2O3 materials revealed by this model was consistent with the Wulff construction. Our results point out that starting by surface thermodynamics is a helpful approach to predict and assess the morphology transformation of materials and they are very useful in interpreting changes in crystal morphology. By this way, this method allow to understand how tuning the surface can control the shape of micro and nanocrystals synthetically. We expect that our work is beneficial to get a profound understanding of how to achieve shape-controlled materials synthesis than the normal trial and error synthetic procedure. Acknowledgments The authors are grateful to for Prometeo II/2014/022 and ACOMP/2014/270 projects (Generalitat Valenciana), Ministerio de Economía y Competitividad (Spain), CTQ2012- 36253-C03-02, FAPESP (2013/07296-2 and scholarship process 2013/26671–9 and 2012/14468-12014/04350-9), CAPES (process A104/2013 and 99999.002998/2014-09), CNPq INCTMN 573636/2008-7, and the Spanish-Brazilian program (PHB2009-0065-PC) for financially supporting this research. References [1] McIlroy A et al 2006 Basic Research Needs for Clean and Efficient Combustion of 21st Century Transportation Fuels (US Dept. of Energy, Basic Energy Sciences, 2007) (http://science.energy. gov./~/media/bes/pdf/reports/files/ctf_rpt.pdf) [2] Kuang Q, Wang X, Jiang Z, Xie Z and Zheng L 2014 Acc. Chem. Res. 47 308–18 [3] Chen K, Sun C, Song S and Xue D 2014 CrystEngComm. 16 5257–67 [4] Tompsett D A, Parker S C, Bruce P G and Islam M S 2013 Chem. Mater. 25 536–41 [5] Kim Y, Lee H and Kang S 2012 J. Mater. Chem. 22 12874–81 [6] Karim A, Fosse S and Persson K A 2013 Phys. Rev. B 87 075322 [7] Eunseok L and Kristin A P 2013 Nanotechnology 24 424007 Table 3. Values of Esurf, number of broken bonds, area, and broken bonding density (Db) calculated for In2O3. In2O3 Surface Esurf (J m−2) Broken bonds Area (nm2) Db (nm−2) (1 0 0) 1.759 48 1.024 46.90 (1 1 0) 1.070 24 0.7238 33.16 (1 1 1) 0.891 42 1.773 23.69 M M Ferrer et alModelling Simul. Mater. Sci. Eng. 24 (2016) 025007 http://science.energy.gov./~/media/bes/pdf/reports/files/ctf_rpt.pdf http://science.energy.gov./~/media/bes/pdf/reports/files/ctf_rpt.pdf http://dx.doi.org/10.1021/ar400092x http://dx.doi.org/10.1021/ar400092x http://dx.doi.org/10.1021/ar400092x http://dx.doi.org/10.1039/C4CE00339J http://dx.doi.org/10.1039/C4CE00339J http://dx.doi.org/10.1039/C4CE00339J http://dx.doi.org/10.1021/cm303295f http://dx.doi.org/10.1021/cm303295f http://dx.doi.org/10.1021/cm303295f http://dx.doi.org/10.1039/c2jm31145c http://dx.doi.org/10.1039/c2jm31145c http://dx.doi.org/10.1039/c2jm31145c http://dx.doi.org/10.1103/PhysRevB.87.075322 http://dx.doi.org/10.1103/PhysRevB.87.075322 http://dx.doi.org/10.1088/0957-4484/24/42/424007 http://dx.doi.org/10.1088/0957-4484/24/42/424007 9 [8] Hörmann N G and Groß A 2014 J. Solid State Electrochem. 18 1401–13 [9] Li Y et al 2014 Phys. Chem. Chem. Phys. 16 24604–9 [10] Barnard A S 2012 Acc. Chem. Res. 45 1688–97 [11] Barnard A S and Kirkland A I 2008 Chem. Mater. 20 5460–3 [12] Seyed-Razavi A, Snook I K and Barnard A S 2010 J. Mater. Chem. 20 416–21 [13] Gurlo A 2011 Nanoscale 3 154–65 [14] Zasada F et al 2011 J. Phys. Chem. C 115 6423–32 [15] Stroppa D G et al 2011 Chem. – Eur. J. 17 11515–9 [16] Beltrán A, Andrés J, Longo E and Leite E R 2003 Appl. Phys. Lett. 83 635–7 [17] Leite E R, Giraldi T R, Pontes F M, Longo E, Beltrán A and Andrés  J 2003 Appl. Phys. Lett. 83 1566–8 [18] Stroppa D G et al 2014 Phys. Chem. Chem. Phys. 16 1089–94 [19] Bomio M R D et al 2013 J. Phys. Chem. C 117 21382–95 [20] Longo V M et al 2011 J. Phys. Chem. C 115 20113–9 [21] Qi K et al 2014 J. Phys. Chem. C 118 23320–7 [22] Piskorz W et al 2012 J. Phys. Chem. C 116 19307–20 [23] Whiteside  A, Fisher  C  A  J, Parker  S  C and Saiful Islam  M 2014 Phys. Chem. Chem. Phys. 16 21788–94 [24] Li P et al 2012 Cryst. Growth Des. 12 5042–50 [25] Suleiman I et al 2015 Phys. Chem. Chem. Phys. 17 7038–45 [26] Kanaki E, Gohr S, Mueller C and Paulus B 2015 Surf. Sci. 632 158–63 [27] Barmparis G D, Lodziana Z, Lopez N and Remediakis I N 2015 Beilstein J. Nanotechnol. 6 361–8 [28] Lovette M A, Browning A R, Griffin D W, Sizemore  J P, Snyder R C and Doherty M F 2008 Ind. Eng. Chem. Res. 47 9812–33 [29] Chen J, Lim B, Lee E P and Xia Y 2009 Nano Today 4 81–95 [30] Jiang Z-Y, Kuang Q, Xie Z-X and Zheng L-S 2010 Adv. Funct. Mater. 20 3634–45 [31] Lee K, Kim M and Kim H 2010 J. Mater. Chem. 20 3791–8 [32] Grena R et al 2015 J. Phys. Chem. C 119 15618–26 [33] Yang H G et al 2008 Nature 453 638–41 [34] Zasada F, Piskorz W and Sojka Z 2015 J. Phys. Chem. C 119 19180–91 [35] Xie X, Li Y, Liu Z-Q, Haruta M and Shen W 2009 Nature 458 746–9 [36] Andrés J, Gracia L, Gouveia A F, Ferrer M M and Longo E 2015 Nanotechnology 26 405703–13 [37] Roca R A et al 2015 Catal. Sci. Technol. 5 4091–107 [38] Botelho G, Andres J, Gracia L, Matos L S and Longo E 2015 Chem Plus Chem. doi: 10.1002/ cplu.201500485 [39] Fabbro M T et al 2015 Sci. Technol. Adv. Mater. 16 065002 [40] Wulff G 1901 Z. Kristallogr. 34 449–530 [41] Gao Z-Y, Sun W and Hu Y-H 2014 Trans. Nonferr. Met. Soc. 24 2930–7 [42] Mulinari T A, La Porta F A, Andres J, Cilense M, Varela J A and Longo E 2013 CrystEngComm. 15 7443–9 [43] Sun H, Ang H M, Tade M O and Wang S 2013 J. Mater. Chem. A 1 14427–42 [44] Yuan C, Zhang L, Hou L, Pang G and Oh W-C 2014 RSC Adv. 4 14408–13 [45] Singh V, Kosa M, Majhi K and Major D T 2014 J. Chem. Theory Comput. 11 64–72 [46] Su D, Dou S and Wang G 2014 Sci. Rep. 4 5767 [47] Gao R, Zhu J, Xiao X, Hu Z, Liu J and Liu X 2015 J. Phys. Chem. C 119 4516–23 [48] Mu J, Zhang L, Zhao G and Wang Y 2014 Phys. Chem. Chem. Phys. 16 15709–16 [49] Xiao J, Kuang Q, Yang S, Xiao F, Wang S and Guo L 2013 Sci. Rep. 3 2300 [50] Cao H, Wang G, Zhang L, Liang Y, Zhang S and Zhang X 2006 ChemPhysChem. 7 1897–901 [51] Braun A et al 2012 J. Phys. Chem. C 116 16870–5 [52] Hisatomi T et al 2012 Faraday Discuss. 155 223–32 [53] Sivula K 2013 J. Phys. Chem. Lett. 4 1624–33 [54] McKone J R, Lewis N S and Gray H B 2013 Chem. Mater. 26 407–14 [55] Mayer M T, Lin Y, Yuan G and Wang D 2013 Acc. Chem. Res. 46 1558–66 [56] Li X, Lei Y, Li X, Song S, Wang C and Zhang H 2011 Solid State Sci. 13 2129–36 [57] Lv B, Xu Y, Wu D and Sun Y 2011 CrystEngComm. 13 7293–98 [58] Wang W, Howe J Y and Gu B 2008 J. Phys. Chem. C 112 9203–8 [59] Liu Z, Lv B, Wu D and Sun Y 2012 Particuology 10 456–61 [60] Liu Z, Lv B, Wu D, Sun Y and Xu Y 2013 Particuology 11 327–33 M M Ferrer et alModelling Simul. Mater. Sci. Eng. 24 (2016) 025007 http://dx.doi.org/10.1007/s10008-013-2189-x http://dx.doi.org/10.1007/s10008-013-2189-x http://dx.doi.org/10.1007/s10008-013-2189-x http://dx.doi.org/10.1039/C4CP03628J http://dx.doi.org/10.1039/C4CP03628J http://dx.doi.org/10.1039/C4CP03628J http://dx.doi.org/10.1021/ar3000184 http://dx.doi.org/10.1021/ar3000184 http://dx.doi.org/10.1021/ar3000184 http://dx.doi.org/10.1021/cm8001892 http://dx.doi.org/10.1021/cm8001892 http://dx.doi.org/10.1021/cm8001892 http://dx.doi.org/10.1039/B915383G http://dx.doi.org/10.1039/B915383G http://dx.doi.org/10.1039/B915383G http://dx.doi.org/10.1039/C0NR00560F http://dx.doi.org/10.1039/C0NR00560F http://dx.doi.org/10.1039/C0NR00560F http://dx.doi.org/10.1021/jp200581s http://dx.doi.org/10.1021/jp200581s http://dx.doi.org/10.1021/jp200581s http://dx.doi.org/10.1002/chem.201100972 http://dx.doi.org/10.1002/chem.201100972 http://dx.doi.org/10.1002/chem.201100972 http://dx.doi.org/10.1063/1.1594837 http://dx.doi.org/10.1063/1.1594837 http://dx.doi.org/10.1063/1.1594837 http://dx.doi.org/10.1063/1.1605241 http://dx.doi.org/10.1063/1.1605241 http://dx.doi.org/10.1063/1.1605241 http://dx.doi.org/10.1039/C3CP53427H http://dx.doi.org/10.1039/C3CP53427H http://dx.doi.org/10.1039/C3CP53427H http://dx.doi.org/10.1021/jp407416h http://dx.doi.org/10.1021/jp407416h http://dx.doi.org/10.1021/jp407416h http://dx.doi.org/10.1021/jp205764s http://dx.doi.org/10.1021/jp205764s http://dx.doi.org/10.1021/jp205764s http://dx.doi.org/10.1021/jp503891s http://dx.doi.org/10.1021/jp503891s http://dx.doi.org/10.1021/jp503891s http://dx.doi.org/10.1021/jp3050059 http://dx.doi.org/10.1021/jp3050059 http://dx.doi.org/10.1021/jp3050059 http://dx.doi.org/10.1039/C4CP02356K http://dx.doi.org/10.1039/C4CP02356K http://dx.doi.org/10.1039/C4CP02356K http://dx.doi.org/10.1021/cg3009927 http://dx.doi.org/10.1021/cg3009927 http://dx.doi.org/10.1021/cg3009927 http://dx.doi.org/10.1039/C4CP05340K http://dx.doi.org/10.1039/C4CP05340K http://dx.doi.org/10.1039/C4CP05340K http://dx.doi.org/10.1016/j.susc.2014.10.007 http://dx.doi.org/10.1016/j.susc.2014.10.007 http://dx.doi.org/10.1016/j.susc.2014.10.007 http://dx.doi.org/10.3762/bjnano.6.35 http://dx.doi.org/10.3762/bjnano.6.35 http://dx.doi.org/10.3762/bjnano.6.35 http://dx.doi.org/10.1021/ie800900f http://dx.doi.org/10.1021/ie800900f http://dx.doi.org/10.1021/ie800900f http://dx.doi.org/10.1016/j.nantod.2008.09.002 http://dx.doi.org/10.1016/j.nantod.2008.09.002 http://dx.doi.org/10.1016/j.nantod.2008.09.002 http://dx.doi.org/10.1002/adfm.201001243 http://dx.doi.org/10.1002/adfm.201001243 http://dx.doi.org/10.1002/adfm.201001243 http://dx.doi.org/10.1039/b921857b http://dx.doi.org/10.1039/b921857b http://dx.doi.org/10.1039/b921857b http://dx.doi.org/10.1021/acs.jpcc.5b00684 http://dx.doi.org/10.1021/acs.jpcc.5b00684 http://dx.doi.org/10.1021/acs.jpcc.5b00684 http://dx.doi.org/10.1038/nature06964 http://dx.doi.org/10.1038/nature06964 http://dx.doi.org/10.1038/nature06964 http://dx.doi.org/10.1021/acs.jpcc.5b05136 http://dx.doi.org/10.1021/acs.jpcc.5b05136 http://dx.doi.org/10.1021/acs.jpcc.5b05136 http://dx.doi.org/10.1038/nature07877 http://dx.doi.org/10.1038/nature07877 http://dx.doi.org/10.1038/nature07877 http://dx.doi.org/10.1088/0957-4484/26/40/405703 http://dx.doi.org/10.1088/0957-4484/26/40/405703 http://dx.doi.org/10.1088/0957-4484/26/40/405703 http://dx.doi.org/10.1039/C5CY00331H http://dx.doi.org/10.1039/C5CY00331H http://dx.doi.org/10.1039/C5CY00331H http://dx.doi.org/10.1002/cplu.201500485 http://dx.doi.org/10.1002/cplu.201500485 http://dx.doi.org/10.1088/1468-6996/16/6/065002 http://dx.doi.org/10.1088/1468-6996/16/6/065002 http://dx.doi.org/10.1016/S1003-6326(14)63428-2 http://dx.doi.org/10.1016/S1003-6326(14)63428-2 http://dx.doi.org/10.1016/S1003-6326(14)63428-2 http://dx.doi.org/10.1039/c3ce41215f http://dx.doi.org/10.1039/c3ce41215f http://dx.doi.org/10.1039/c3ce41215f http://dx.doi.org/10.1039/c3ta12960h http://dx.doi.org/10.1039/c3ta12960h http://dx.doi.org/10.1039/c3ta12960h http://dx.doi.org/10.1039/c4ra00762j http://dx.doi.org/10.1039/c4ra00762j http://dx.doi.org/10.1039/c4ra00762j http://dx.doi.org/10.1021/ct500770m http://dx.doi.org/10.1021/ct500770m http://dx.doi.org/10.1021/ct500770m http://dx.doi.org/10.1038/srep05767 http://dx.doi.org/10.1038/srep05767 http://dx.doi.org/10.1021/jp511363p http://dx.doi.org/10.1021/jp511363p http://dx.doi.org/10.1021/jp511363p http://dx.doi.org/10.1039/c4cp01326c http://dx.doi.org/10.1039/c4cp01326c http://dx.doi.org/10.1039/c4cp01326c http://dx.doi.org/10.1002/cphc.200600130 http://dx.doi.org/10.1002/cphc.200600130 http://dx.doi.org/10.1002/cphc.200600130 http://dx.doi.org/10.1021/jp304254k http://dx.doi.org/10.1021/jp304254k http://dx.doi.org/10.1021/jp304254k http://dx.doi.org/10.1039/C1FD00103E http://dx.doi.org/10.1039/C1FD00103E http://dx.doi.org/10.1039/C1FD00103E http://dx.doi.org/10.1021/jz4002983 http://dx.doi.org/10.1021/jz4002983 http://dx.doi.org/10.1021/jz4002983 http://dx.doi.org/10.1021/cm4021518 http://dx.doi.org/10.1021/cm4021518 http://dx.doi.org/10.1021/cm4021518 http://dx.doi.org/10.1021/ar300302z http://dx.doi.org/10.1021/ar300302z http://dx.doi.org/10.1021/ar300302z http://dx.doi.org/10.1016/j.solidstatesciences.2011.08.011 http://dx.doi.org/10.1016/j.solidstatesciences.2011.08.011 http://dx.doi.org/10.1016/j.solidstatesciences.2011.08.011 http://dx.doi.org/10.1039/c1ce05431g http://dx.doi.org/10.1039/c1ce05431g http://dx.doi.org/10.1039/c1ce05431g http://dx.doi.org/10.1021/jp800683j http://dx.doi.org/10.1021/jp800683j http://dx.doi.org/10.1021/jp800683j http://dx.doi.org/10.1016/j.partic.2011.11.012 http://dx.doi.org/10.1016/j.partic.2011.11.012 http://dx.doi.org/10.1016/j.partic.2011.11.012 http://dx.doi.org/10.1016/j.partic.2012.06.016 http://dx.doi.org/10.1016/j.partic.2012.06.016 http://dx.doi.org/10.1016/j.partic.2012.06.016 10 [61] Chakrabarty S and Chatterjee K 2013 J. Cryst. Growth 381 107–13 [62] Lützenkirchen J et al 2015 Faraday Discuss. 180 55–79 [63] Liang C et al 2015 New J. Chem. 39 2651–6 [64] Guo H and Barnard A S 2011 J. Mater. Chem. 21 11566–77 [65] Mehta B R and Singh V N 2005 Pramana – J. Phys. 65 949–58 [66] Walsh A and Catlow C R A 2010 J. Mater. Chem. 20 10438–44 [67] Motta F V, Lima R C, Marques A P A, Leite E R, Varela J A and Longo E 2010 Mater. Res. Bull. 45 1703–6 M M Ferrer et alModelling Simul. Mater. Sci. Eng. 24 (2016) 025007 http://dx.doi.org/10.1016/j.jcrysgro.2013.07.013 http://dx.doi.org/10.1016/j.jcrysgro.2013.07.013 http://dx.doi.org/10.1016/j.jcrysgro.2013.07.013 http://dx.doi.org/10.1039/C4FD00260A http://dx.doi.org/10.1039/C4FD00260A http://dx.doi.org/10.1039/C4FD00260A http://dx.doi.org/10.1039/C4NJ02032D http://dx.doi.org/10.1039/C4NJ02032D http://dx.doi.org/10.1039/C4NJ02032D http://dx.doi.org/10.1039/c1jm10381d http://dx.doi.org/10.1039/c1jm10381d http://dx.doi.org/10.1039/c1jm10381d http://dx.doi.org/10.1007/BF02704096 http://dx.doi.org/10.1007/BF02704096 http://dx.doi.org/10.1007/BF02704096 http://dx.doi.org/10.1039/c0jm01816c http://dx.doi.org/10.1039/c0jm01816c http://dx.doi.org/10.1039/c0jm01816c http://dx.doi.org/10.1016/j.materresbull.2010.06.056 http://dx.doi.org/10.1016/j.materresbull.2010.06.056 http://dx.doi.org/10.1016/j.materresbull.2010.06.056