Interlaced crystals having a perfect Bravais lattice and complex chemical order revealed by real-space crystallography

Abstract

The search for optimal thermoelectric materials aims for structures in which the crystalline order is disrupted to lower the thermal conductivity without degradation of the electron conductivity. Here we report the synthesis and characterisation of ternary nanoparticles (two cations and one anion) that exhibit a new form of crystalline order: an uninterrupted, perfect, global Bravais lattice, in which the two cations exhibit a wide array of distinct ordering patterns within the cation sublattice, forming interlaced domains and phases. Partitioning into domains and phases is not unique; the corresponding boundaries have no structural defects or strain and entail no energy cost. We call this form of crystalline order 'interlaced crystals' and present the example of hexagonal CuInS₂ . Interlacing is possible in multi-cation tetrahedrally bonded compound with an average of two electrons per bond. Interlacing has minimal effect on electronic properties, but should strongly reduce phonon transport, making interlaced crystals attractive for thermoelectric applications.

ARTICLE

Received 18 Feb 2014 | Accepted 30 Sep 2014 | Published 14 Nov 2014

Xiao Shen1,*, Emil A. Hernndez-Pagan2,3,*, Wu Zhou4,*, Yevgeniy S. Puzyrev1, Juan-Carlos Idrobo5, Janet E. Macdonald2,3, Stephen J. Pennycook6 & Sokrates T. Pantelides1,3,4,7

The search for optimal thermoelectric materials aims for structures in which the crystalline order is disrupted to lower the thermal conductivity without degradation of the electron conductivity. Here we report the synthesis and characterisation of ternary nanoparticles (two cations and one anion) that exhibit a new form of crystalline order: an uninterrupted, perfect, global Bravais lattice, in which the two cations exhibit a wide array of distinct ordering patterns within the cation sublattice, forming interlaced domains and phases. Partitioning into domains and phases is not unique; the corresponding boundaries have no structural defects or strain and entail no energy cost. We call this form of crystalline order interlaced crystals and present the example of hexagonal CuInS2. Interlacing is possible in multi-cation tetrahedrally bonded compound with an average of two electrons per bond. Interlacing has minimal effect on electronic properties, but should strongly reduce phonon transport, making interlaced crystals attractive for thermoelectric applications.

1 Department of Physics and Astronomy, Vanderbilt University, Nashville, Tennessee 37235, USA. 2 Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37235, USA. 3 Vanderbilt Institute for Nanoscale Science and Engineering, Vanderbilt University, Nashville, Tennessee 37235, USA. 4 Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA. 5 Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA. 6 Department of Materials Science and Engineering, The University of Tennessee, Knoxville, Tennessee 37996, USA. 7 Department of Electrical Engineering and Computer Science, Vanderbilt University, Nashville, Tennessee 37235, USA. * These authors contributed equally to this work. Correspondence and requests for materials should be addressed to X.S. (email: mailto:[email protected]

Web End [email protected] ) or to W.Z. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 5:5431 | DOI: 10.1038/ncomms6431 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

& 2014 Macmillan Publishers Limited. All rights reserved.

DOI: 10.1038/ncomms6431

Interlaced crystals having a perfect Bravais lattice and complex chemical order revealed by real-space crystallography

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6431

We discovered the novel crystalline structures during an investigation of CuInS2 semiconductor nanoparticles that are actively studied for solar energy harvesting

and other applications15. They often have the chalcopyrite structure of bulk CuInS2, that is, a zinc-blende structure in which

Cu and In share the cation sublattice1,2. Some methods of synthesis, however, yield nanoparticles that have a hexagonal structure, the wurtzite analogue of chalcopyrite1,35. These hexagonal CuInS2 nanoparticles are believed to be single crystals with an ordered sulfur sublattice, but it has not been established whether the Cu and In atoms are randomly distributed or ordered (it should be noted that the words lattice and structure are used differently by the physics and crystallography community; we use the physics terminology where lattice is a mathematical abstraction).

From theory, ordered cations are favoured over random distributions of cations according to the GrimmSommerfeld (GS) or valence-octet rule6, which has been conrmed by density functional theory (DFT) calculations for bulk cubic I-III-VI27,8 and related quaternary materials9. In cubic CuInS2 thin lms, X-ray diffraction (XRD) and transmission electron microscopy (TEM)10 have established the coexistence of two ordered phases: chalcopyrite, which is the ground state, and a CuAu-like phase, which, according to theory8, is only slightly higher in energy. For nanoparticles, however, reciprocal-space X-ray crystallographic techniques yield patterns with broad, weak peaks. Thus, any small peaks and subtle changes in the pattern that would be diagnostic of the presence of one or several ordered phases are not discernible. Instead, only the average lattice (hexagonal in this case) is determined and the question of cation ordering in the wurtzite phase remains unanswered. Although the octet rule, which is backed by DFT calculations, enshrines the notion that ordered cations are energetically preferred over a random distribution, it does not specify if there is a preferred ordering or set of orderings or any potential distribution. The aim of the present paper is to present the experimental evidence for the manifestation of interlaced crystals and a theoretical rational for their existence.

Here we use real-space crystallography in the form of Z-contrast imaging on an aberration-corrected scanning TEM (STEM), to directly distinguish the In and Cu atoms in the nanoparticles. We also employ DFT to determine an entire family of ordered hexagonal CuInS2 phases (polytypes) that are essentially degenerate in energy and are favoured over phases with disordered cations in accord with the GS rule, as is the case for cubic CuInSe2 (ref. 7). The two approaches independently yield the same result: the Cu and In cations are not randomly distributed in the cation sublattice of a wurtzite crystal structure. Instead the abstract cation sublattice retains its perfect structure while the two cation species exhibit a range of ordered structures. The key experimental nding is that the Z-contrast images show a continuum of interlaced phases and domains that can be demarcated in a variety of ways as illustrated in Fig. 1, with an uninterrupted, perfect Bravais lattice. At the same time, the calculations found that phase and domain boundaries are defect-and strain-free and cost no energy, and conrm that the global Bravais lattice remains intact. As a further evidence, strain analysis of the Z-contrast images reveals that there is no measurable strain. We will further demonstrate that the partitioning into domains and phases is not unique and that domain boundaries are twins by either rotational or reection symmetry (see Supplementary Information), while phase boundaries have no twinning symmetry operation. It is the combination of a perfect, uninterrupted global Bravais lattice and the nonuniqueness of partitioning, with strain-free, defect-free domain and phase boundaries that characterises what we call

interlaced crystals. Such crystallographic features are likely to exist in many similar compounds.

ResultsNanoparticle characterisation. Nanoparticles of CuInS2 were synthesized by a solution-phase colloidal synthesis at 245 C 11. Characterisation by TEM, XRD and UVvis absorption spectrum analysis were performed. TEM shows that the nanoparticles have a pancake shape, with a diameter of 182 nm and thickness of 82 nm (Supplementary Fig. 1). Powder XRD (Fig. 2a) veried that the structure is hexagonal, as it would be in a binary compound with the wurtzite structure. Scherrer broadening calculations of the strong (100)-like reection at 26.8 yields a crystallite size of 20 nm (Supplementary Fig. 2 and Supplementary Note 1), which is consistent with the TEM analysis. The XRD provides supporting evidence that there is an uninterrupted global Bravais lattice in the interlaced crystals, as the 20-nm size calculated from Scherrer broadening on wurtzite reections is far larger than the 35-nm lateral size of cation-ordered domains observed in the STEM image (Fig. 1a).

While the wurtzite reections are clearly observed in the XRD pattern of nanoparticles, the reections due to cation ordering are too low in intensity to be discerned. This is due to the combined effect of several factors: (1) the different ordered phases do not necessarily generate superlattice reections at the same 2y (Fig. 2), whereby the presence of multiple cation-ordered phases reduces the signal from each phase proportionally (in addition to the WZ1WZ4 structures shown in Fig. 3, other cation-ordered phases with similar band gaps are possible and would also weaken the superlattice reections); (2) the natural intensities of these reections are intrinsically low, as can be seen from Fig. 2cf;(3) Scherrer broadening is greater for the cation-ordered reections than for the overarching wurtzite reections (Supplementary Note 1).

Indirect information about cation ordering is provided by the UVvis absorption spectrum, which shows that the nanoparticles have a band gap of 1.55 eV (Supplementary Fig. 3), close to that of the chalcopyrite CuInS2 (1.54 eV)12. Indeed, DFT calculations nd that a disordered cation sublattice results in a much smaller band gap than that of chalcopyrite CuInS2 (Supplementary

Table 1).

Real-space crystallography. Figure 1a shows an unprocessed STEM Z-contrast image from the pancake-shaped CuInS2 nanocrystals, taken along the c axis of the wurtzite lattice. Figure 1b shows a simulated image based on the atomic positions extracted from the image in Fig. 1a, assuming there is only one type of cation. It is clear from this image that the abstract cation sublattice of the structure is a perfect hexagonal Bravais lattice without any structural interfaces, such as stacking faults, twins or grain boundaries. As we will now demonstrate, however, the cations form several different ordered structures within the cation sublattice, creating interlaced domains that can be demarcated by domain boundaries, but the partitioning is not unique. The perfect Bravais anion and cation sublattices, and interlaced domains are the key features of what we call interlaced crystals.

The colourized version of the STEM Z-contrast image (Fig. 1c,d) shows clearly that there are distinct Cu and In columns (blue and yellow, respectively), which means that the majority of Cu and In atoms do not occupy random sites in the cation lattice, as that would result in a roughly uniform intensity in all atomic columns throughout the nanoparticle. The Cu and In columns form a complex pattern that can be decomposed into a set of fundamental cation-ordered phases (WZ1WZ4) shown in Fig. 3 and domains (as further conrmed by atomic-resolution

2 NATURE COMMUNICATIONS | 5:5431 | DOI: 10.1038/ncomms6431 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6431 ARTICLE

2 nm

WZ4*

WZ4

WZ1 WZ4

WZ2

WZ4*

0.35

WZ2

WZ3

WZ2*

WZ2* 0.17

Figure 1 | Direct imaging of cation ordering in an interlaced CuInS2 nanoparticle. (a) Unprocessed Z-contrast STEM annular dark eld (ADF) image of a CuInS2 nanoparticle. The image is taken along the c axis of the wurtzite lattice, which is also the vertical axis of the pancake-shaped nanocrystal. (b) Simulated image based on the atomic positions extracted from a, with the same intensity at all atomic sites. (c,d) Different ways to partition the nanoparticle, demonstrating the interlacing nature of the cation ordering. The partitioning is done on the colourized low by-pass ltered ADF image to highlight the difference between Cu and In columns. The colour scale bar is shown in d. The dark level of the images is set to 0.17 and the maximum level is set to 0.35 for contrast enhancement, so that columns containing the heavier In cations display as yellow dots while columns containing the lighterCu cations display as blue dots. (See Supplementary Fig. 4c,d for more details.) The notion WZ2 and WZ4 refer to phases of CuInS2 with different ways of cation ordering. The notation WZ2/WZ2* (and WZ4/WZ4*) is used to label different cation-ordering domains of the same phase butwith different orientations. Over two hundred high-resolution images were taken from three batches of nanoparticle samples, and the interlaced crystal structure was consistently observed.

Experimental

Disordered

WZ1

WZ2

Intensity (a.u.)

WZ3

WZ4

10 20 30 50

40 60 70 80

2

Figure 2 | Experimental and simulated XRD patterns. (a) XRD patterns of CuInS2 nanoparticles, (b) simulated pattern for a hexagonal wurtzite structure with disordered cation sublattice and (cf) simulated patterns for the four different ordered structures (see Fig. 3 and Methods).

elemental mapping shown in Supplementary Fig. 5). The coexistence of various cation-ordered phases is also conrmed by secondary spots in the simulated electron diffraction pattern (Supplementary Fig. 6) and Fourier transform of the STEM image (Supplementary Fig. 4b). In addition to the chemical ordering, strain analysis of the Z-contrast images reveals that there is no measurable strain (Supplementary Fig. 7).

DFT calculations. The peculiar structure of the interlaced crystals described in this paper originates from the bonding requirements in tetrahedrally coordinated I-III-VI2 compounds such as CuInS2.

To have an average of four electrons associated with each atom, which is optimal for tetrahedral bonds, the group-VI atom (S) must have two group I (Cu) neighbours and two group III (In) neighbours according to the GS or valence-octet rule9. In contrast, in isovalent ternary compounds like AlGaAs, the GS rule is always satised regardless of the distribution of cations. In these alloys, a disordered cation sublattice is favoured by entropy, as long as the size mismatch between the cations is not large13.

For cubic CuInS2, the GS rule is satised in two known cubic structures: one is the chalcopyrite structure that is the ground state for CuInS2; the other is the CuAu-like ordered structure, which is slightly higher in energy8 and has been observed in thin lms10. Wei et al.7 identied a family of cubic polytypes for the similar CuInSe2 whose energies are only slightly higher than that of the chalcopyrite ground state. We found that for the wurtzite-like hexagonal CuInS2, the GS rule results in a family of structures, including the four structures that are depicted in detail in Fig. 3 and were independently predicted by DFT calculations. The rst two structures, which we call WZ1 and WZ2, are constructed by hexagonal analogies to the chalcopyrite and CuAu-like structures, respectively. The third structure, WZ3, was obtained by running a reverse Monte Carlo simulation14, which ensured that S atom has exactly 50% probability of nding Cu atom as its nearest neighbour. The common feature of structures WZ1, WZ2 and WZ3 is that their lattices are constructed by repeating a Cu1In1 unit in one direction and a CuxIny unit in another direction (with a conjugated InxCuy arrangement a

NATURE COMMUNICATIONS | 5:5431 | DOI: 10.1038/ncomms6431 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6431

half-cell below). Therefore, there exists a family of an innite number of orderings, characterised by non-negative integers x and y, that satises the GS/valence-octet rule. For the WZ1 structure, x y 2. For the WZ2 structure, x N and y 0. In

the case of the WZ3 structure, x 3 and y 1. The next simplest

member of the family is x 2 and y 1, which we label as WZ4.

(For the structure data, see Supplementary Note 2). DFT total-energy calculations show that these structures are very close in energy (enthalpy), with WZ1 being slightly lower by o4 meV per

CuInS2 unit. As in the case of the cubic CuInSe2 polytypes7, all four structures have signicantly lower energy than structures with disordered cation sublattices (Supplementary Table 2), where each unit of deviation from the GS rule (one Cu-S/In-S bond pair) costs 0.3B0.4 eV of energy.

It is possible to interpret some of the structures as derivatives of others. For example, the WZ3 structure can be viewed as a derivative of WZ2 with twinning domains while WZ4 can be viewed as a derivative of WZ2 with slips. However, they are better described as distinct, ordered structures, because the derivative descriptions imply structural defects and do not fully capture the translational symmetries.

One can extend the family of WZn structures. While the repetition of a Cu1In1 unit in one direction is always required, the repetition of a CuxIny unit in the second direction is not necessary. We found that any combination of Cu and In, Cux1

Iny1Cux2Iny2y, can result in a structure in which every S has two Cu and two In neighbours if there is a conjugated Inx1Cuy1 Inx2Cuy2y arrangement a half-cell below. This feature makes the arrangement of Cu and In columns very exible, allowing the coexistence of various phases of cation ordering and different orientations of cation-ordering domains, as observed experimentally in the nanoparticles. The calculated energies of domain and phase boundaries are essentially zero (Supplementary Table 3).

a2

a1 a1

a1

a2

a2

a2

a1

Figure 3 | Projection view of the four structures along the [001] direction with ordered cation sublattices. (a) WZ1, (b) WZ2, (c) WZ3 and (d) WZ4. Cu atoms are in blue, In atoms are in yellow and S atoms are in grey. The repeating Cu1In1 and CuxIny units along the direction of the lattice vectors in the wurtzite structure are shown by red arrows. The conjugated InxCuy units are shown as the dotted red arrows. The primitive unit cells are shown by green dashed lines.

WZ2 WZ4

WZ4

WZ2

Figure 4 | Nonunique partition of ordered domains in the interlaced crystal. (a) Z-contrast STEM image of a nanoparticle, with the area of interest marked. (b) Enlarged view of the area of interest. (c,d) Two different ways to partition the area into WZ2 and WZ4 structures, demonstrating the interlaced lattice. The solid white line shows the position of the phase boundary. The dashed white arrow marks the In5CuIn2CuIn2 arrangement of Cu and In atoms discussed in the text. Each green and red box shows a unit cell of the repeating WZ4 structure on the right.

4 NATURE COMMUNICATIONS | 5:5431 | DOI: 10.1038/ncomms6431 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6431 ARTICLE

DiscussionOne important consequence of the existence of a perfect global Bravais lattice is that the structure can be partitioned into domains and phases in a nonunique way as illustrated in Fig. 1c,d, where we show two different ways of partitioning part of the image into the fundamental ordered structures. In Fig. 1c, we divide the structure into six parts, two WZ2 domains, two WZ4 domains, a WZ1 phase and a WZ3 phase. In Fig. 1d, we divide the same particle into only four parts: two WZ2 and two WZ4 domains. The nonuniqueness of partitioning occurs because the different phases are not stacked along the z-direction (c axis), which would yield a superlattice of polytypes. In such stacking, the phase boundaries would not obey the GS rule. In these nanoparticles, however, different phases are primarily juxtaposed as columns in the xy plane, that is, the nanoparticles can be described as textured15 and the GS rule is obeyed without interruption.

The 60 boundary between two WZ2 domains is equivalent to a piece of WZ3 structure, and the boundary between two WZ4 domains is equivalent to a piece of WZ1 structure. A boundary between two domains of the same phase yielding a different stable phase is a feature that has not been observed in other material systems.

The interlaced nature of the crystal structure not only shows up at domain boundaries, but also at phase boundaries. The boundary between WZ2 and WZ4 phases is assigned differently in Fig. 1c,d. To illustrate the phenomenon more clearly, we zoom in to the area around the boundary (Fig 4a,b). An In5CuIn2

CuIn2 arrangement is marked with a dotted white arrow that crosses the two phases. The position of the phase boundary could be assigned at multiple locations with the extremes illustrated in Fig. 4c,d. Therefore, the middle section belongs to both the WZ2 and the WZ4 structures, making the two phases interlaced. We note that, in addition to the I-III-VI2 and related materials, there have been other cases where two cations share the same sublattice and cation ordering occurs in domains, for example, in perovskite-structure transition-metal oxides16,17. However, in these systems the local changes in stoichiometry increase the Coulombic repulsions between cations, perturb the anion sublattice and raise the free energy16. Such distortions and associated energy costs are typical in all cases of structural boundaries, such as grain boundaries, twin boundaries and stacking faults. Even coherent twin boundaries were recently reported to entail distortions and energy costs18. In contrast, in what we call interlaced crystals, the local bonding at both domain and phase boundaries is the same as in the bulk crystal, without any distortions, disruptions or energy cost.

Although the domain boundaries in CuInS2 can be identied as a form of twinning19,20, they exhibit an unbroken global Bravais lattice, which is not a feature of any known twin boundaries (Supplementary Fig. 8 and Supplementary Note 3). The phase boundaries in CuInS2, on the other hand, are denitely not twins as they lack a twinning symmetry operation. All boundaries are dened strictly through the ordering of the cations in the otherwise perfect Bravais lattice, which underlies the nonuniqueness of the partitioning. It is for these reasons that the order observed in CuInS2 nanoparticles is a prototype of a new form of crystalline order.

Although the size and shape of the nanoparticles have facilitated this discovery, (the small thickness allows for clear Z-contrast between the In and Cu columns, and the pancake-shaped particles lie at so that the c axis is conveniently observed with Cu and In columns easily distinguished), each nanoparticle contains B50,000 atoms, and thus surface energy contributions are likely low. Therefore, the observed crystalline behaviour is bulk-like and not intrinsic to nanoparticles. We expect that

interlaced crystals should exist in thin lms and bulk materials as well. If the size of ordering phase and domain in these materials is B5 nm as observed in the present study, XRD analysis is likely to be insufcient to identify the cation-ordered phases. The very weak reections (Fig. 2cf) due to cation ordering will have a line broadening of 1.8 (Supplementary Note 1), which makes them difcult to be observed. In addition to hexagonal CuInS2 and other I-III-VI2 compounds, interlacing is likely to occur in cubic versions of these materials, in II-IV-V2 compounds such as ZnGeAs2 or ZnGeP2 and in I2-II-IV-VI4 quaternaries. Such materials have been synthesized and have been found to exhibit ordering in the cation sublattice in accord with the GS rule and DFT calculations21,22. In addition, anti-phase boundaries in ZnGeAs2 were observed by TEM 25 years ago21, and could be described by the interlacing phenomena. In all cases modern Z-contrast imaging could verify the presence of interlaced crystal ordering.

Interlaced crystals are expected to have contrasting electronic and phonon mobilities due to the nonunique assignment of phase boundaries and domain structure. Unlike grain boundaries in polycrystalline materials, the local bonding at domain and phase boundaries in an interlaced crystal is the same as in the bulk, with no mismatched or distorted bonds. Therefore, the domain and phase boundaries of interlaced crystals should have minimal effect on electron transport. In contrast, phonon transport should be reduced signicantly. Phonon scattering happens at not only phase boundaries because of their different phonon spectra, but also at domain boundaries, as it is impossible to conserve both the phonon wavevector and vibrational eigenvector due to the different orientation of the domains. The existence of such a phonon scattering mechanism at domain boundaries has been demonstrated at ferroelectric domain walls23. As a result, interlaced crystals should have high electron conductivities but low phonon conductivities, which make them ideal for thermoelectric applications.

Recently, a high gure of merit of 1.4 was reported for bulk chalcopyrite CuGaTe2 (ref. 24). Interlaced crystalline order is possible in such materials and is the likely source of the high thermoelectric performance. Future efforts to control crystallisation in the mixed cation copper chalcogenides and encourage a high density of interlaced crystal domains, should yield new thermoelectric materials with even higher gures of merit.

Methods

Nanocrystal synthesis. CuCl (99.99%, STREM), InCl3 (99.99%, Alfa Aesar), thiourea (ACS reagent Zreagen Sigma-Aldrich), oleylamine (technical grade 70%

Aldrich) were used as received. The CuInS2 nanoparticles were synthesized following the procedure reported in ref. 11, CuCl (0.5 mmol), InCl3 (0.5 mmol), thiourea (1.0 mmol) and oleylamine (10.0 ml) were loaded on a 25-ml 3-neck ask containing a stir bar. The ask was equipped with a condenser connected to a vacuum/purging adaptor, as well as a thermocouple and two septa. The mixture was placed under vacuum at 60 C for 30 min, then heated to 245 C under nitrogen and held at this temperature for 1 h after which the heating mantle was removed. After cooling to room temperature, the particles were precipitated with ethanol followed by centrifugation and removal of the supernatant. Three cycles were performed of suspension in hexanes:oleylamine (20:1, by volume) followed by precipitation by excess ethanol, centrifugation and removal of the supernatant. Finally, the nanocrystals were stored in hexanes:oleylamine.

Low-resolution TEM, XRD and UVvis absorption spectrum. Low-resolution TEM was carried out on a Phillips CM20 operating at 200 kV. The powder XRD pattern was collected over 43 h on an Inel powder diffractometer equipped with a curved position-sensitive detector (CPS 120) and a Cu Ka1 radiation source operated at 40 kV and 20 mA. The software Powder Cell was employed for the simulation of the diffraction patterns from DFT results of ordered structures. To t the positions of major peaks to that of experiment, the DFT lattice constants are adjusted (the a and c parameters of the underlying wurtzite lattice are scaled to 3.86 and 6.37 , respectively). UVvis absorption spectrum was collected on a Jasco V-670 spectrophotometer.

NATURE COMMUNICATIONS | 5:5431 | DOI: 10.1038/ncomms6431 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6431

Reverse Monte Carlo. The WZ3 structure and two structures with disordered cation sublattices were generated using a reverse Monte Carlo algorithm in a 128-atom supercell, which corresponds to a 4 4 2 supercell of the wurzite structure.

The cost function in reverse Monte Carlo simulation is 1

N

Psites

2, where nI is the number of group I atoms and N is coordination number (which equals four in this material). The starting state of the simulation is a random distribution of 50% of Cu and 50% of In atoms in the cation site. For generating the WZ3 structure, the minimisation of cost function stops when it equals zero, which corresponds to every S atom in the simulation cell having exactly 50% probability of nding Cu as nearest neighbour (two Cu atoms out of four nearest neighbours). For generating the two structures with disordered cation sublattices, the mini-misation stopped when the cost function is close to zero, therefore ensuring that each group-VI atom had close to 50% probability of nding a group I atom as its nearest neighbour, while retaining some randomness in the cation sublattice, which is characterised by an increased number of same-type neighbour pairs in the cation sublattice as compared with the ordered congurations.

DFT calculations. The band gaps are obtained from hybrid density functional theory calculations using the PBE0 version of the exchange-correlation functional25. The cohesive energies are obtained using DFT calculations with the PBE version of the exchange-correlation functional26. We employed PAW potentials27 and plane-wave basis as implemented in the VASP code28. We used 400 eV for the plane-wave cutoff. The geometrical optimizations are converged to 10 4 eV for the energy difference between two ionic steps. The energies of domain and phase boundaries are obtained using DFT calculations with the PBE functional. The calculations are done in 12 2 1 supercell of the wurzite

structure with 96 atoms, which contains a pair of the boundaries. The energies of the boundaries are calculated by subtracting the total energies of the respective bulk phases from the total energy of the boundary model.

High-resolution electron microscopy characterisation. The high-resolution STEM imaging and electron energy loss spectroscopy analysis were performed with an aberration-corrected Nion UltraSTEM-100 operated at 100 kV accelerating voltage. The convergence semi-angle for the incident probe was 31 mrad and the annular dark eld images were acquired with a collection semi-angle range of about 86200 mrad. EEL spectra were collected using a Gatan Enna spectrometer with an EELS collection semi-angle of 48 mrad. Energy dispersion of 1 eV per channel was used for the elemental mapping. The atomic position in the annular dark eld image was extracted using the standard Fine Maxima in ImageJ. Electron diffraction simulation was performed using SingleCrystal plug-in within Crystal-Maker (http://www.crystalmaker.com/

Web End =http://www.crystalmaker.com/). Strain eld analysis was performed using a GPA analysis plug-in (http://www.christophtkoch.com/strain/index.html

Web End =http://www.christophtkoch.com/strain/index.html).

References

1. Pan, D. et al. Synthesis of Cu In S ternary nanocrystals with tunable

structure and composition. J. Am. Chem. Soc. 130, 56205621 (2008).2. Castro, S. L., Bailey, S. G., Raffaelle, R. P., Banger, K. K. & Hepp, A. F. Synthesis and characterisation of colloidal CuInS2 nanoparticles from a molecular single-source precursor. J. Phys. Chem. B 108, 1242912435 (2004).

3. Qi, Y. et al. Synthesis and characterisation of nanostructured wurtzite CuInS2:

a new cation disordered polymorph of CuInS2. J. Phys. Chem. C 113, 39393944 (2009).4. Connor, S. T., Hsu, C.-M., Weil, B. D., Aloni, S. & Cui, Y. Phase transformation of biphasic Cu2S CuInS2 to monophasic CuInS2 nanorods. J. Am. Chem. Soc.

131, 49624966 (2009).5. Norako, M. E., Franzman, M. A. & Brutchey, R. L. Growth kinetics of monodisperse Cu In S nanocrystals using a dialkyl disulde sulfur source.

Chem. Mater. 21, 42994304 (2009).6. Miller, A., Mackinnon, A. & Weaire, D. in Solid State Physics vol. 36 (eds Seitz, F. & Turnbu, D.) 119175 (Academic Press, 1981).

7. Wei, S.-H., Ferreira, L. G. & Zunger, A. First-principles calculation of the order-disorder transition in chalcopyrite semiconductors. Phys. Rev. B 45, 25332536 (1992).

8. Wei, S.-H., Zhang, S. B. & Zunger, A. Band structure and stability of zinc-blende-based semiconductor polytypes. Phys. Rev. B 59, R2478R2481 (1999).

9. Chen, S., Walsh, A., Luo, Y., Gong, X. G. & Wei, S.-H. Wurtzite-derived polytypes of kesterite and stannite quaternary chalcogenide semiconductors. Phys. Rev. B 82, 195203 (2010).

10. Su, D. S. et al. CuAu-type ordering in epitaxial CuInS2 lms. Appl. Phys. Lett. 73, 785787 (1998).

11. Koo, B., Patel, R. N. & Korgel, B. A. Wurtzite chalcopyrite polytypism in

CuInS2 nanodisks. Chem. Mater. 21, 19621966 (2009).12. Yoshino, K., Ikari, T., Shirakata, S., Miyake, H. & Hiramatsu, K. Sharp band edge photoluminescence of high-purity CuInS2 single crystals. Appl. Phys. Lett.

78, 742744 (2001).

13. Wood, D. M., Wei, S.-H. & Zunger, A. Stability and electronic structure of ultrathin [001] (GaAs)m(AlAs)m superlattices. Phys. Rev. B 37, 13421363 (1988).

14. McGreevy, R. L. Reverse Monte Carlo modeling. J. Phys. Condens. Matter 13, R877R913 (2001).

15. Wenk, H.-R. & Van Houtte, P. Texture and anisotropy. Rep. Prog. Phys. 67, 1367 (2004).

16. Davis, P. K., Tong, J. & Negas, T. Effect of ordering-induced domain boundaries on low-loss Ba(Zn1/3Ta2/3)O3-BaZrO3 perovskite microwave

dielectrics. J. Am. Ceram. Soc. 80, 17271740 (1997).17. Sun, J., Liu, S., Newman, N. & Smith, D. J. Atomic resolution transmission electron microscopy of the microstructure of ordered Ba(Cd1/3Ta2/3)O3

perovskite ceramics. J. Am. Ceram. Soc. 89, 10471052 (2006).18. Wang, Y. M. et al. Defective twin boundaries in nanotwinned metals. Nat. Mater. 12, 697702 (2013).

19. Hahn, T. h. & Klapper, H. in International Tables for Crystallography vol. D (eds Authier, A.) 415 (Wiley, 2014).

20. Janovec, V. & Privratska, J. in International Tables for Crystallography vol. D (eds Authier, A.) 506508 (Wiley, 2014).

21. Soloman, G. S., Posthill, J. B. & Timmons, M. L. Antiphase domain boundary formation in singe-crystal chalcopyrite-structure ZnGeAs2 grown on GaAs.

Appl. Phys. Lett. 55, 15311533 (1989).22. Shimony, Y., Fledman, R., Dahan, I. & Kimmel, G. Anti-phase domain boundaries in ZnGeP2 (ZGP). Opt. Mater. 16, 119123 (2001).

23. Weilert, M. A., Msall, M. E., Anderson, A. C. & Wolfe, J. P. Phonon scattering from ferroelectric domain walls: phonon imaging in KDP. Phys. Rev. Lett. 71, 735738 (1993).

24. Plirdpring, T. et al. Chalcopyrite CuGaTe2: a high-efciency bulk thermoelectric material. Adv. Mater. 24, 36223626 (2012).

25. Adamo, C. & Barone, V. Toward reliable density functional methods without adjustable parameters: the PBE0 model. J. Chem. Phys. 110, 61586170 (1999).

26. Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 38653868 (1996).

27. Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 59, 17581775 (1999).

28. Kresse, G. & Furthmuller, J. Efcient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 1116911186 (1996).

Acknowledgements

This work is supported by a Wigner Fellowship through the Laboratory Directed Research and Development Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for DOE (W.Z.), by the National Science Foundation grant no. DMR-0938330 (W.Z., for work at Vanderbilt prior to 2013), by NSF EPS-1004083 (J.E.M., E.A.H.-P.), by NSF CHE-1253105 (J.E.M.), by the US-Israel Binational Foundation and the Bergmann Memorial Award (J.E.M.), by ORNLs Center for Nanophase Materials Sciences (CNMS), which is sponsored by the U.S. DOE Ofce of Basic Energy Sciences, Scientic User Facilities Division, (J.-C.I.), by the U.S. DOE Ofce of Science, Basic Energy Sciences, Materials Sciences and Engineering Division (S.J.P., X.S.), and by DOE grant DE-FG02-0946554 (S.T.P., Y.S.P.). The work used resources of the National Energy Research Scientic Computing Center, supported by the DOE Ofce of Science under Contract No. DE-AC02-05CH11231. We thank Lasse Norn, Ray Withers, Barry Carter, and Laurie Marks for useful discussions and Young-Min Kim and Ryo Ishikawa, and Christoph Koch for their help in strain analysis. We thank Jonathan G.C. Veinot and Tapas Purkait at the University of Alberta Chemistry Department for their assistance in obtaining the XRD pattern.

Author contributions

X.S. and S.T.P. coordinated the project and performed DFT calculations; Y.S.P. performed Monte Carlo simulations; E.A.H.-P. and J.E.M. synthesized and characterised the samples; W.Z., J.-C.I. and S.J.P. performed aberration-corrected scanning transmission electron microscopy; X.S.,W.Z., J.E.M. and S.T.P. wrote the manuscript.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Shen, X. et al. Interlaced crystals having a perfect Bravais lattice and complex chemical order revealed by real-space crystallography. Nat. Commun. 5:5431 doi: 10.1038/ncomms6431 (2014).

n

N

12

6 NATURE COMMUNICATIONS | 5:5431 | DOI: 10.1038/ncomms6431 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.