Home | About | Journals | Submit | Contact Us | Français |

**|**IUCr Open Access Articles**|**PMC3098556

Formats

Article sections

Authors

Related links

Acta Crystallogr B. 2011 June 1; 67(Pt 3): 205–217.

Published online 2011 May 14. doi: 10.1107/S0108768111013814

PMCID: PMC3098556

Liang Li,^{a} Alexander Wölfel,^{a} Andreas Schönleber,^{a} Swastik Mondal,^{a} Antoine M. M. Schreurs,^{b} Loes M. J. Kroon-Batenburg,^{b} and Sander van Smaalen^{a,}^{*}

Correspondence e-mail: smash/at/uni-bayreuth.de

Received 2011 January 25; Accepted 2011 April 12.

Copyright © Liang Li et al. 2011

This is an open-access article distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.

A combination of structure refinements, analysis of the superspace MEM density and interpretation of difference-Fourier maps has been used to characterize the incommensurate modulation of rubidium tetrachlorozincate, Rb_{2}ZnCl_{4}, at a temperature of *T* = 196 K, close to the lock-in transition at *T*
_{lock-in} = 192 K. The modulation is found to consist of a combination of displacement modulation functions, modulated atomic displacement parameters (ADPs) and modulated third-order anharmonic ADPs. Up to fifth-order Fourier coefficients could be refined against diffraction data containing up to fifth-order satellite reflections. The center-of-charge of the atomic basins of the MEM density and the displacive modulation functions of the structure model provide equivalent descriptions of the displacive modulation. Modulations of the ADPs and anharmonic ADPs are visible in the MEM density, but extracting quantitative information about these modulations appears to be difficult. In the structure refinements the modulation parameters of the ADPs form a dependent set, and *ad hoc* restrictions had to be introduced in the refinements. It is suggested that modulated harmonic ADPs and modulated third-order anharmonic ADPs form an intrinsic part, however small, of incommensurately modulated structures in general. Refinements of alternate models with and without parameters for modulated ADPs lead to significant differences between the parameters of the displacement modulation in these two types of models, thus showing the modulation of ADPs to be important for a correct description of the displacive modulation. The resulting functions do not provide evidence for an interpretation of the modulation by a soliton model.

The construction of a model-independent electron-density map from phased structure factors is an important application of the maximum entropy method (MEM) in crystallography (Takata, 2008 ; van Smaalen & Netzel, 2009 ). Analysis of the electron density after such a reconstruction provides the locations of the atoms in the unit cell. As such, the MEM has been used to determine the locations of the metal atoms in endohedral fullerenes (Nishibori *et al.*, 2006 ), to obtain information about disorder (multiple positions) in crystal structures (Dinnebier *et al.*, 1999 ; Wang *et al.*, 2001 ; Samy *et al.*, 2010 ) and about anharmonic atomic displacements (Kumazawa *et al.*, 1995 ; Bagautdinov *et al.*, 1998 ).

The MEM has been generalized towards the determination of the generalized electron density in -dimensional superspace (*d* = 1, 2, 3,…) of aperiodic crystals (van Smaalen *et al.*, 2003 ). Again, the MEM provides information about the locations of the atoms, which then result in a description of the modulation functions of incommensurately modulated crystals or incommensurate composite crystals (Palatinus & van Smaalen, 2004 ; van Smaalen & Li, 2009 ). Alternatively, the MEM in superspace has been used to determine the occupation domains of the atoms in quasicrystals (Yamamoto *et al.*, 1996 ). Here we will use the MEM to obtain information about the modulation functions of incommensurately modulated RbZnCl.

Many isostructural compounds of the -KSO structure type undergo phase transitions on cooling. Several compounds exhibit at least two phase transitions, first forming an incommensurately modulated structure which then becomes commensurate at lower temperatures (lock-in transition; Cummins, 1990 ).

Rubidium tetrachlorozincate, RbZnCl, is one of these compounds (Fig. 1 ). RbZnCl undergoes a phase transition from a periodic to an incommensurately modulated phase at = 303 K. The incommensurate modulation wavevector is = (). The lock-in transition towards a threefold superstructure () takes place at = 192 K (Sawada *et al.*, 1977 ).

The modulation of RbZnCl increasingly deviates from a sinusoidal shape on approaching the lock-in transition, as shown by the growth of the intensities of higher-order satellite reflections in the X-ray diffraction of this compound on cooling toward (Aramburu *et al.*, 1997 ). The results of structure refinements of a model of displacive modulation functions with up to fifth-order Fourier coefficients have been reported by Aramburu *et al.* (2006 ). The latter authors interpreted this structure model as providing evidence for a soliton shape of the incommensurate modulation wave.

Here we present a re-analysis of the incommensurate structure of RbZnCl close to the lock-in transition, employing a more extensive data set of Bragg reflections than has been used by Aramburu *et al.* (2006 ). The purpose of this work is twofold. The first aim is to investigate the nature of modulations by means of the maximum entropy method (MEM). As we will show, modulations of atomic displacement parameters (ADPs) and modulations of anharmonic ADPs form an intrinsic and important part of the modulation. Secondly we show that the modulation functions do not provide evidence for a soliton character of the modulation in this compound.

Single crystals of RbZnCl have been grown from aqueous solution (Sawada *et al.*, 1977 ). RbCl (2.73 g, Aldrich, 99.99%) and ZnCl (1.54 g, Aldrich, 99.999%) were dissolved in 4.5 g of ultra pure water (from a Simplicity UV system by Millipore) at *T* = 323 K. The solution was slowly cooled to *T* = 313 K, and crystals were obtained by slow evaporation at this temperature.

A suitable single crystal was glued to a thin glass fibre mounted on a copper pin. X-ray diffraction experiments were performed at beamline F1 of Hasylab, DESY, Hamburg, employing the radiation of a wavelength of 0.5000 Å and a MAR-CCD area detector. The temperature of the sample was maintained at *T* = 196 K, employing a nitrogen-flow cryostat. A large crystal-to-detector distance of 225 mm was chosen, in order to be able to resolve closely spaced reflections.

With the aid of the four-circle kappa diffractometer at beamline F1, diffraction data were collected by and scans with a scan step of 0.3° per image. Several values were chosen for the off-set of the detector and for the orientation of the crystal, thus allowing the measurement of a nearly complete data set up to a high resolution of = 0.86 Å^{−1}. With the purpose of increasing the effective dynamic range of the experiment, runs with a zero detector off-set were repeated with exposure times of 2 and 8 s, and runs at higher scattering angles were repeated with 8 and 64 s exposure. The long exposure times resulted in overexposed strong (main) reflections, while they allowed weak reflections (mostly higher-order satellite reflections) to be measured.

Integrated intensities of Bragg reflections were extracted from the measured images by the software *EVAL*15 (Schreurs *et al.*, 2010 ). Absorption correction was performed with *SADABS* (Sheldrick, 2008 ). A fraction of the area of the CCD detector was not properly cooled during parts of the experiment. This is a technical problem that occurred for experiments of long durations (Paulmann, 2009 ). As a result several pixels of the detector always gave a large intensity, which could negatively affect data quality. Therefore, the coordinates of these pixels have been determined by inspection of the images, and they were excluded from the integration. Experimental data and crystallographic information are summarized in Table 1 .^{1} The observed volume of the unit cell is significantly smaller than reported by Aramburu *et al.* (2006 ), who gave = 844.04 Å with *a* = 7.241 (3), *b* = 12.648 (5) and *c* = 9.216 (3) Å. Since lattice parameters from point-detector measurements are much more accurate than from area detectors, we have employed the lattice parameters from Aramburu *et al.* (2006 ) in the present refinements.

The resulting data set of intensities of Bragg reflections — including satellite reflections up to fifth order — was used for structure determination, structure refinements and maximum entropy calculations.

Aramburu *et al.* (2006 ) have kindly supplied the diffraction data from their publication. These data will be denoted as the Aramburu data. Various models have also been tested by calculation of the values of indices on these data.

A peculiar property of the Aramburu data is that a selection of the reflections were measured, which included all main reflections and only the strongest satellite reflections as expected on the basis of a soliton model. Satellite reflections up to order five, except fourth order, have been measured in this way by Aramburu *et al.* (2006 ). The result is a data set that consists of many fewer reflections than available in the present data. On the other hand, CCD detectors have a limited dynamic range so that the lower bound on measurable intensities is relatively high, resulting in the number of high-order ‘observed’ satellite reflections being comparable in the two data sets (Table 1 ).

Structure models of different complexity have been refined against the diffraction data. They involve the basic structure coordinates and the harmonic atomic displacement parameters (ADPs) for each of the six crystallographically independent atoms (Fig. 1 ). Depending on the complexity of the model, they may include Fourier coefficients for displacement modulation ( and for the sine and cosine Fourier coefficients of the order along the direction ); anharmonic ADPs of third () and fourth () order; Fourier coefficients for the modulation of the ADPs (, for the sine and cosine Fourier coefficients of order ) as well as and (Table 2 ; van Smaalen, 2007 ).

Structure refinements were performed with the computer program *JANA*2006 (Petricek *et al.*, 2006 ). The model published by Aramburu *et al.* (2006 ) involves displacement modulation parameters of orders 1, 2, 3 and 5. Refinement of these parameters against the Aramburu data reproduced the published model within one standard uncertainty () of all parameters.

Model A was created to resemble the published structure model as much as possible. It includes all Fourier coefficients up to fifth order for the displacement modulation, because the availability of fourth-order satellite reflections in the present data allows the refinement of the fourth-order Fourier coefficients of the displacement modulations. Refinements were initiated with the values of the published structure model as starting parameters. Values of the refined parameters are similar to those of the published structure model, with only 12 out of 140 parameters having differences larger than and with a maximum difference of for of atom Rb1 (*cf.* Table 3 with Table VI in Aramburu *et al.*, 2006 ).

Model B is an extension of model A, where the first- and second-order Fourier coefficients of the modulation of the harmonic ADPs have been incorporated. Refinements with model A as starting values for the parameters gave a smooth convergence and led to a considerable improvement of the fit to all orders of reflections (Table 4 ).

Model B was used to create the phased observed diffraction data from the measured intensities for the MEM calculations (see §2.3). Analysis of the MEM-derived electron-density map suggested that the next important feature was the modulation of the third-order anharmonic ADPs, while their average structure values remained zero. Model C includes, in addition to the parameters of model B, the Fourier coefficients up to for the modulation of the third-order anharmonic ADPs, . This refinement suffered from large correlations between parameters. Therefore, a reduced model, model C_{r}, was defined, in which those Fourier coefficients were set to zero that had values less than in the refinement of model C. This reduced the number of coefficients from 244 to 132 (Table 2 ), while models C and C_{r} fitted the data almost equally well (Table 4 ).

Difference-Fourier maps based on the observed structure factors and those calculated for a model indicate the improvement of the fit to the data for increasing complexity of the model (Fig. 2 and Table 4 ). The difference-Fourier map of model B compared with that of model C_{r} confirms the importance of modulated third-order anharmonic ADPs, as it has been derived based on MEM density maps. The difference-Fourier map of model C_{r} displays structure around the Rb2 atom which, to a first approximation, is independent of the phase of the modulation. It has the signature of unmodulated fourth-order anharmonic ADPs, as they are missing in model C_{r}. The inclusion of fourth-order anharmonic ADPs for all atoms led to highly nonphysical values of these parameters, that is, large negative values of the joint probability distribution function for the resulting model. Model D_{r} was then constructed to include fourth-order anharmonic ADPs for atoms Rb1, Rb2 and Cl3 only. The improvement, compared with model C_{r}, of the fit to the data, in particular to the main reflections, is apparent (Table 4 ). Refinements of the extinction coefficient led to a negative value for this parameter, so it was fixed to zero.

(, )-sections of difference-Fourier maps centered on the Rb2 atom (, and ) from different models, showing two periods along . Solid lines represent positive values, dashed lines are negative values, and long dashed lines represent the zero contour. The **...**

The remaining discrepancies between calculated and observed structure factors can be attributed in part to the incompleteness of the model. As indicated above, the introduction of more parameters leads to nonphysical values and high correlations between them, while these additional parameters would have been required for a full characterization of the modulation. A second reason for the rather high final values of the higher-order satellite reflections lies in the less than optimal accuracy of the present data due to limited counting statistics. This interpretation becomes apparent when the values are considered for model D_{r} on the stronger reflections of the present data [reflections with ; column in Table 4 ]. In particular, the partial values of the higher-order satellite reflections are considerably lower than on the full data set (compare columns D and D_{r} in Table 4 ).

The fit of the models A, B, C, C_{r} and D_{r} to the Aramburu data has been tested by refinement of the basic structure parameters of each model against these data, while the modulation parameters and anharmonic ADPs were kept fixed to the values determined from the present data. The fit to the main reflections and first-order satellite reflections is reasonable, but it becomes worse on the introduction of modulation parameters for the (an)harmonic ADPs (models B–D_{r}; Table 5 ). On the other hand, the latter models lead to an improvement of the fit to the third- and fifth-order satellites of the Aramburu data, but with values that are considerably higher than those on the present data. These discrepancies can be attributed to different qualities of the sample and especially different temperatures, which will affect the shapes of the modulation functions and the contributions of modulated and anharmonic ADPs to it.

Therefore, independent refinements were performed against the Aramburu data, now varying all parameters, and resulting in models A, B, C, C and D, which differ from the corresponding models A, B, C, C_{r} and D_{r} in the values of the parameters. The fit to the Aramburu data is dramatically improved in this way (see supplementary material), resulting in values comparable to values on the present data. Exceptions are the main reflections, which are much better fitted for the present data, indicating the higher accuracy of these data compared with the Aramburu data.

Despite convergence of the refinements against the Aramburu data and the resulting low values, the primed models suffer from high correlations between parameters and large standard uncertainties. For example, none of the modulation parameters for ADPs in model B exceed , which prevents a meaningful analysis of the modulation on the basis of model B, as has already been noted by Aramburu *et al.* (2006 ).

The standard uncertainties of modulation parameters and anharmonic ADPs are a multiple of the standard uncertainties of these parameters in the corresponding unprimed models (refinements against the present data). Therefore, we refrain from a further consideration of the primed models.

Phased observed structure factors corrected for anomalous scattering and scaled to the scattering power of one unit cell were obtained from the observed data and model B according to published procedures (Bagautdinov *et al.*, 1998 ). These data were used for the calculation of a maximum-entropy-optimized generalized electron density in -dimensional superspace [MEM density or ] with the computer program *BAYMEM* (van Smaalen *et al.*, 2003 ). A uniform prior, the Cambridge algorithm and the weights of type F2 have been used (Li *et al.*, 2010 ). The MEM calculation converged in 69 iterations (see Table 1 for more information on the MEM calculation).

The -dimensional electron-density map has been analyzed with the computer program *EDMA* (van Smaalen *et al.*, 2003 ). Physical space sections of have been obtained for 100 equally spaced values within one period along the fourth axis, *i.e.* for . Atoms in the crystal correspond to local maxima in the physical space sections of the generalized electron density.

The position of each local maximum as a function of then provides an estimate for the modulated position of an atom. Alternatively, the center-of-charge has been determined for the atomic basins surrounding each local maximum. The dependence on of the positions of the center-of-charge provides an alternative measure for the atomic positions. Modulation functions have been extracted from by taking the difference between the modulated atomic position and the basic structure position as obtained from model B (Fig. 3 ).

Modulation functions of the crystallographically independent atoms of RbZnCl. Displacements along , and are given in Å. Open circles represent the center of the charge and filled circles are the maxima of the MEM electron density. Lines represent **...**

Two-dimensional sections of have been visualized by the plotting option of the computer program *JANA*2006 (Petricek *et al.*, 2006 ). The () section centered on the Rb2 atom clearly shows the modulated position of this atom (Fig. 4 ).

The incommensurate modulation of RbZnCl at a temperature of *T* = 196 K, close to the lock-in transition at K, has been determined to be comprised of atomic displacement modulation functions that contain contributions of Fourier coefficients up to fifth order. This finding is in agreement with previous studies on this compound (Aramburu *et al.*, 2006 ).

Modulations of the ADPs are found to be an intrinsic part of the incommensurate crystal structure. That is, the harmonic ADPs are modulated with up to second-order Fourier coefficients and the third-order anharmonic ADPs are modulated with Fourier coefficients up to fifth order, while the basic structure components of the third-order anharmonic ADPs are zero.

The finding of modulated ADPs is in agreement with studies on several other compounds, such as incommensurately modulated NaCO (Dusek *et al.*, 2003 ), and the composite crystals [LaS][NbS] and [SrO][CrO] (Jobst & van Smaalen, 2002 ; Castillo-Martínez *et al.*, 2008 ). Modulated ADPs have also been found necessary in cases where only first-order satellite reflections were available in the diffraction data, like LaC, NaSiO and PbNiVO (Dusek *et al.*, 2000 ; Krüger *et al.*, 2006 ; Roussel *et al.*, 2009 ).

The necessity of modulation functions for third-order anharmonic ADPs has been revealed in our previous studies with the MEM on (NH)BeF (Palatinus *et al.*, 2004 ) and CrPO (Li *et al.*, 2010 ). Many incommensurate crystal structures have been published, where values are higher than they should have been for the perfect structure model. It can thus be speculated that the fit to diffraction data might be improved for many compounds by the inclusion of modulated ADPs and modulated third-order anharmonic ADPs. On the other hand, correlations between modulation parameters, as shown here for RbZnCl, might prevent their determination by structure refinements. This problem especially exists for high-order Fourier coefficients of modulation functions. Meaningful values are almost always limited to coefficients of orders equal to and less than the maximum order of observed satellite reflections.

It is suggested here that modulations of ADPs are an intrinsic part of modulations in incommensurate crystals in general. Their presence can be rationalized by the fact that any displacive modulation defines a modulation of the environments of the atoms. Different environments require different ADPs, which can be achieved by a modulation of the ADPs. It is important to include at least up to second-order Fourier coefficients (Perez-Mato *et al.*, 1991 ).

The displacement modulation of an atom creates a tightening of its environment in the direction in which this atom is displaced. Since the modulation always involves displacements out of the average position into both the positive and negative directions, this general feature of modulations explains the presence of modulated third-order anharmonic ADPs, while their average values are zero. The correlation between displacement modulation and modulated is apparent from the values found for model D_{r} of RbZnCl. Both the displacement modulation and the modulated have their most important nonzero components along the axis for all five independent atoms on the mirror plane (Table 6 and supplementary material). The Cl3 atom, not on the mirror plane, also has contributions to its modulation for the other directions, which again affects both the displacement modulation and the modulated third-order anharmonic ADPs.

Although not perfectly matched, positive displacements along (Fig. 5 ) of the Rb2 atom match negative values of (Fig. 6 ), which implies a decreased probability (Figs. 7 and 8 ) for the presence of an atom in the direction of a tighter environment.

Modulation functions of the crystallographically independent atoms of RbZnCl. Displacements along , and are given in Å. Open circles are the center of the charge and filled circles are the maxima of the MEM electron density. Dashed lines represent **...**

Value in model D_{r} of the component of third-order anharmonic ADPs of the Rb2 atom as a function of . Minimum and maximum values are located at and 0.6.

The joint probability distribution function at the site of Rb2 for selected values for model D_{r}, but obtained with and only. The contour interval is 1 e Å with a maximum density of 17.7 e Å.

The joint probability distribution function at the site of Rb2 for selected values for model D_{r}. The contour interval is 1 e Å. Over the selected map region the minimum density is −0.66 e Å and **...**

The trace of the center-of-charge for each atom in the MEM density indicates smooth modulations, whereas the trace of the local maximum of the density varies around the position defined by the center-of-charge for each atom (Fig. 3 ). We take this variation as evidence for the presence of modulated third-order anharmonic ADPs. Similarly, several of the components of the modulation functions in model A exhibit variations (ripples) according to higher-order Fourier coefficients, while the same functions appear smooth in model D_{r} (Fig. 5 ). (The largest effect is visible for the components [Cl], [Cl] and [Rb].) These smoother functions seem more plausible and they match the trace of the center-of-charge of the MEM density very well. Both refinements (model D_{r}
*versus* model A) and the MEM thus provide evidence for modulated third-order anharmonic ADPs. Lastly, the reduction of values on the increasing complexity of the structure model provides strong evidence for modulated harmonic ADPs and modulated third-order anharmonic ADPs (Table 4 ).

As mentioned above, structure refinements without (model A) and with modulation functions for ADPs (models B–D_{r}) result in significantly different functions for the displacive modulation. The inclusion of modulation functions for ADPs thus appears to be necessary to achieve an accurate description of the displacive modulation, with concomitant implications for the interpretation of the modulation (§3.2). Alternatively, the center-of-charge of each atom in the MEM density also provides a good description of the displacement modulation functions.

Comparison of the two approaches, MEM analysis and structure refinements, shows the different limitations of the two methods. The MEM density gives evidence for the modulations of the ADPs as well as the presence of anharmonic ADPs. However, the finite size of the pixels in the MEM density (here 0.1 Å) limits the accuracy of the atomic positions to ~ 0.01 Å (van Smaalen *et al.*, 2003 ), while atoms on special positions may sometimes lead to more accurate values of the positions. An error of up to 0.01 Å is not small, if modulations are considered with amplitudes significantly below 0.1 Å. On the other hand, structure refinements readily lead to large dependencies between parameters, such that Fourier components of orders cannot be determined. Furthermore, a full *ab initio* determination appeared impossible for the third- and fourth-order anharmonic ADPs, and we had to resort to a method of selecting relevant parameters (compare models C, C_{r} and D_{r} and the discussion in §2.2).

Aramburu *et al.* (2006 ) have shown that a soliton model for the modulation leads to displacement modulation functions with Fourier components of first and fifth (and higher) order. They introduced a measure, , for the soliton density, which can be interpreted as the ratio between the width and the separation of the discommensurations (Babkevich & Cowley, 1999 ), and which describes the shape of the modulation functions with for a sinusoidal shape and for non-overlapping discommensurations. Within this approach, the phase of the fifth-order Fourier coefficient depends on the phase of the first-order Fourier coefficient by a simple relation and it is independent from the soliton density. The ratio between amplitudes of fifth- and first-order Fourier coefficients should be the same for all atoms, while its value is a measure for the soliton density. Aramburu *et al.* (2006 ) have found these relations to be approximately valid for their structure model for RbZnCl, and they proposed that the modulation of RbZnCl can be described by a soliton wave, with a soliton density of at the temperature of their experiment.

Here we have shown that significant differences exist for the displacement modulation functions in cases of a pure displacive modulation model (model A) and a model including modulated (an)harmonic ADPs (model D_{r}). Since model D_{r} is the more accurate model and the Aramburu model resembles model A, this finding questions the interpretation by Aramburu *et al.* (2006 ) concerning the soliton shape of displacive modulations in their model. The data from Aramburu *et al.* (2006 ) are re-plotted in Figs. 9 and 10 . The relation between the phases of the first- and fifth-order Fourier coefficients matches that of a soliton wave much better than the observed standard uncertainties would suggest. This excellent agreement might be an artefact resulting from the fact that the only reflections included in the data set were those which possessed high intensities in the soliton model. On the other hand, one out of eight data points has a phase relation that is significantly different from the soliton model, again suggesting that the Aramburu data do not necessarily provide evidence of the soliton model.

Both models A and D_{r} appear to be at variance with the soliton model, as follows most prominently from the ratio of amplitudes of fifth- and first-order Fourier coefficients (Fig. 10 ). The standard uncertainties of the phases of the fifth Fourier coefficients are much larger for some functions in models A and D_{r} than in the Aramburu model. While standard uncertainties of refined parameters are of comparable magnitude in the different models, this discrepancy can be ascribed to the much smaller amplitudes of some fifth-order coefficients in model D_{r} than in the Aramburu model, thus leading to a poorer estimate of the phases of these functions. Nevertheless, large discrepancies are found between the calculated and experimental phases of the fifth-order Fourier coefficients, again amounting to several standard uncertainties. This indicates that the present data fail to provide evidence for a soliton shape of the modulation functions.

Another feature of the modulation arguing against the soliton model is the presence of Fourier coefficients of orders two, three and four, with magnitudes comparable to the magnitudes of the fifth-order coefficients (Table 6 ). The interpretation favored by Aramburu *et al.* (2006 ) is that the third-order Fourier coefficient would represent a secondary mode, while they have not determined the fourth-order Fourier coefficients and the second-order Fourier coefficients are much smaller in the model by Aramburu *et al.* (2006 ) than presently obtained (Table 6 ). The interpretation of Aramburu *et al.* (2006 ) would thus imply that secondary modes are more important than the distortion (magnitude of fifth harmonics) toward the supposed soliton-shaped wave, a situation that is not necessarily likely.

An alternative interpretation of the observed structure model is that of a squaring of the modulation wave, then involving all harmonics of the modulation functions (Leist *et al.*, 2008 ). This interpretation is not at variance with the interpretation of the incommensurate structure by discommensurations. It only questions the structure of the discommensurations as a structure given by the solution of the sine-Gordon equation. Since the sine-Gordon equation is only a simple model for discommensurations, a more advanced theory might be able to describe the observed shapes of the modulation functions.

Finally, we have presently established that modulated harmonic ADPs and modulated third-order anharmonic ADPs are an important part of the modulation. However, these features have not been incorporated into the soliton model considered by Aramburu *et al.* (2006 ). Establishing the relation between modulated ADPs and a possible soliton property of the modulation wave will require further theoretical analysis that is beyond the scope of the present work.

The origin of the modulations in RbZnCl and in *A*
*BX*
-type compounds in general lies in the incompatibility between the observed orthorhombic packing of ZnCl and Rb ions and the nearly tetrahedral symmetry of the ZnCl complex ions. This incompatibility results in one short distance between the 11-coordinated A atom (Rb1 in the present models) and an *X* atom (Cl1 in the present models) in the same mirror plane. The bond strength of this short bond in the unmodulated high-temperature structure has been taken as a measure for the propensity of the compound to form modulated structures at low temperatures (Fabry & Perez-Mato, 1994 ).

Analysis of the interatomic distances of model D_{r} shows that they are in agreement with previous studies on similar compounds (Friese *et al.*, 2000 ), and that they support the interpretation given by Fabry & Perez-Mato (1994 ); see -plots of distances and bond angles in the supplementary material.

The present model gives displacement modulations of Rb1 and Cl1 that are in-phase with each other (Fig. 5 ). Consequently, the very short Rb1–Cl1 distance hardly varies with phase of the modulation. Instead, the strain of this contact is resolved by the modulated third-order anharmonic ADPs.

A combination of structure refinements, analysis of the superspace MEM density and interpretation of difference-Fourier maps has been used to characterize the incommensurate modulation of RbZnCl at a temperature of *T* = 196 K, close to the lock-in transition at = 192 K. The basic characteristics of the modulation are a displacement modulation that contains contributions of Fourier coefficients up to fifth order.

A modulation of the ADPs is found to be an intrinsic part of the modulation. That is, the harmonic ADPs are modulated with up to second-order Fourier coefficients and the third-order anharmonic ADPs are modulated with Fourier coefficients up to fifth order, while the basic structure or average components of the third-order anharmonic ADPs are zero.

Model D_{r}, which includes modulated ADPs and modulated third-order anharmonic ADPs, provides different values for the parameters of the displacement modulation than model A, which lacks any modulation of ADPs. Modulations of ADPs are thus essential for the correct description of the displacement modulation functions.

The MEM density gives an excellent description of the displacement modulations of the atoms by means of the -dependencies of the traces of the centers-of-charge of the atoms. These traces coincide with the displacement modulation functions of the atoms in model D_{r} but not in model A, providing further support for the necessity of modulated ADPs in the structure model.

Modulations of the ADPs and anharmonic ADPs are visible in the MEM density as variations of the distributions of the density about their average value, as exemplified by the traces of the local maxima of the MEM density around the positions of the atoms. A quantitative interpretation of the MEM density is made difficult by the finite resolution of this map, which limits the accuracy of positions to ~ 0.01 Å.

Structure refinements may lead to a quantitative description of the modulation, but the introduction of the required model parameters readily leads to correlated parameters. Nevertheless, with the extensive data set available in the present study, we have been able to obtain significant values for higher-order Fourier coefficients of the displacive modulation functions and for modulated parameters of the harmonic ADPs and the third-order anharmonic ADPs.

The results suggest that modulated harmonic ADPs and modulated third-order anharmonic ADPs form an intrinsic part, however small, of incommensurately modulated structures.

For RbZnCl we could show that the modulation fails to provide clear evidence for a soliton wave as the principal shape of the modulation functions (Aramburu *et al.*, 2006 ). Instead, an extended theory will be necessary which includes the effects of modulated (an)harmonic ADPs to understand the modulations in *A*
*BX*
compounds.

Crystal structure: contains datablocks global, I, II. DOI: 10.1107/S0108768111013814/bp5035sup1.cif

Click here to view.^{(79K, cif)}

Structure factors: contains datablocks I. DOI: 10.1107/S0108768111013814/bp5035Isup2.hkl

Click here to view.^{(324K, hkl)}

Structure factors: contains datablocks II. DOI: 10.1107/S0108768111013814/bp5035IIsup3.hkl

Click here to view.^{(324K, hkl)}

Supplementary material file. DOI: 10.1107/S0108768111013814/bp5035sup4.pdf

Click here to view.^{(868K, pdf)}

The single crystal was grown by Alfred Suttner at the Laboratory of Crystallography in Bayreuth. We thank Vaclav Petricek, Michal Dusek and Lukas Palatinus (Praha) for help with *JANA*2006. We are grateful to I. Aramburu for providing a copy of their diffraction data of RbZnCl. Financial support was obtained from the German Science Foundation (DFG).

^{1}Supplementary data for this paper are available from the IUCr electronic archives (Reference: BP5035). Services for accessing these data are described at the back of the journal.

- Aramburu, I., Friese, K., Pérez-Mato, J. M., Morgenroth, W., Aroyo, M., Breczewski, T. & Madariaga, G. (2006).
*Phys. Rev. B Condens. Matter*,**73**, 014112. - Aramburu, I., Madariaga, G., Grebille, D., Perez-Mato, J. M. & Breczewski, T. (1997).
*J. Phys. I*,**7**, 371–383. - Babkevich, A. Y. & Cowley, R. A. (1999).
*J. Phys. Condens. Matter*,**11**, 1639–1655. - Bagautdinov, B., Luedecke, J., Schneider, M. & van Smaalen, S. (1998).
*Acta Cryst.*B**54**, 626–634. - Castillo-Martínez, E., Schönleber, A., van Smaalen, S., Arévalo-López, A. M. & Alario-Franco, M. (2008).
*J. Solid State Chem.***181**, 1840–1847. - Cummins, H. Z. (1990).
*Phys. Rep.***185**, 211–409. - Dinnebier, R. E., Schneider, M., van Smaalen, S., Olbrich, F. & Behrens, U. (1999).
*Acta Cryst.*B**55**, 35–44. [PubMed] - Dusek, M., Chapuis, G., Meyer, M. & Petricek, V. (2003).
*Acta Cryst.*B**59**, 337–352. [PubMed] - Dusek, M., Chapuis, G., Schobinger-Papamantellos, P., Wilkinson, C., Petricek, V., Tung, L. D. & Buschow, K. H. J. (2000).
*Acta Cryst.*B**56**, 959–971. [PubMed] - Fabry, J. & Perez-Mato, J. M. (1994).
*Phase Transitions*,**49**, 193–229. - Friese, K., Madariaga, G. & Breczewski, T. (2000).
*Acta Cryst.*B**56**, 17–21. [PubMed] - Jobst, A. & van Smaalen, S. (2002).
*Acta Cryst.*B**58**, 179–190. [PubMed] - Krüger, H., Kahlenberg, V. & Friese, K. (2006).
*Acta Cryst.*B**62**, 440–446. [PubMed] - Kumazawa, S., Takata, M. & Sakata, M. (1995).
*Acta Cryst.*A**51**, 651–658. - Leist, J., Gibhardt, H., Hradil, K. & Eckold, G. (2008).
*J. Phys. Condens. Matter*,**20**, 415209. - Li, L., Schönleber, A. & van Smaalen, S. (2010).
*Acta Cryst.*B**66**, 130–140. [PubMed] - Nishibori, E., Terauchi, I., Sakata, M., Takata, M., Ito, Y., Sugai, T. & Shinohara, H. (2006).
*J. Phys. Chem. B*,**110**, 19215–19219. [PubMed] - Palatinus, L., Amami, M. & van Smaalen, S. (2004).
*Acta Cryst.*B**60**, 127–137. [PubMed] - Palatinus, L. & van Smaalen, S. (2004).
*Z. Kristallogr.***219**, 719–729. - Paulmann, C. (2009). Personal communication.
- Perez-Mato, J. M., Madariaga, G. & Elcoro, L. (1991).
*Solid State Commun.***78**, 33–37. - Petricek, V., Dusek, M. & Palatinus, L. (2006).
*JANA*2006. Institute of Physics, Praha, Czech Republic. - Roussel, P., Labidi, O., Huve, M., Drache, M., Wignacourt, J.-P. & Petricek, V. (2009).
*Acta Cryst.*B**65**, 416–425. [PubMed] - Samy, A., Dinnebier, R. E., van Smaalen, S. & Jansen, M. (2010).
*Acta Cryst.*B**66**, 184–195. [PubMed] - Sawada, S., Shiroishi, Y., Yamamoto, A., Takashige, M. & Matsuo, M. (1977).
*J. Phys. Soc. Jpn*,**43**, 2099–2100. - Schreurs, A. M. M., Xian, X. & Kroon-Batenburg, L. M. J. (2010).
*J. Appl. Cryst.***43**, 70–82. - Sheldrick, G. M. (2008).
*SADABS*, Version 2008/1. University of Göttingen, Germany. - Takata, M. (2008).
*Acta Cryst.*A**64**, 232–245. [PubMed] - van Smaalen, S. (2007).
*Incommensurate Crystallography.*Oxford University Press. - van Smaalen, S. & Li, L. (2009).
*Phys. Scr.***79**, 048305. - van Smaalen, S. & Netzel, J. (2009).
*Phys. Scr.***79**, 048304. - van Smaalen, S., Palatinus, L. & Schneider, M. (2003).
*Acta Cryst.*A**59**, 459–469. [PubMed] - Wang, C.-R., Kai, T., Tomiyama, T., Yoshida, T., Kobayashi, Y., Nishibori, E., Takata, M., Sakata, M. & Shinohara, H. (2001).
*Angew. Chem. Int. Ed.***40**, 397–399. [PubMed] - Yamamoto, A., Weber, S., Sato, A., Kato, K., Ohshima, K.-I., Tsai, A. P., Niikura, A., Hiraga, K., Inoue, A. & Masumoto, T. (1996).
*Philos. Mag. Lett.***73**, 247–254.

Articles from Acta Crystallographica Section B: Structural Science are provided here courtesy of **International Union of Crystallography**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's Canada Institute for Scientific and Technical Information in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |