Molecular Orientation of Intercalants Stabilized in the Interlayer Space of Layered Ceramics: 1-D Electron Density Simulation
Article information
Abstract
In this review, an attempt is made to calculate one-dimensional (1-D) electron density profiles from experimentally determined (00l) XRD intensities and possible structural models as well in an effort to understand the collective intracrystalline structures of intercalant molecules of two-dimensional (2-D) nanohybrids with heterostructures. 2-D ceramics, including layered metal oxides and clays, have received much attention due to their potential applicability as catalysts, electrodes, stabilizing agents, and drug delivery systems. 2-D nanohybrids based on such layered ceramics with various heterostructures have been realized through intercalation reactions. In general, the physico-chemical properties of such 2-D nanohybrids are strongly correlated with their heterostructures, but it is not easy to solve the crystal structures due to their low crystallinity and high anisotropic nature. However, the powder X-ray diffraction (XRD) analysis method is thought to be the most powerful means of understanding the interlayer structures of intercalant molecules. If a proper number of well-developed (00l) XRD peaks are available for such 2-D nanohybrids, the 1-D electron density along the crystallographic c-axis can be calculated via a Fourier transform analysis to obtain structural information about the orientations and arrangements of guest species in the interlayer space.
1. Introduction
Layered ceramic materials such as cationic clays, anionic clays such as layered double hydroxides (LDHs), layered metal oxides, layered metal chalcogenides, layered zirconium phosphates and carbonaceous layered compounds can be transformed into novel heterostructured materials through an intercalation reaction.1–7) Intercalation is defined as a chemical process in which guest species are reversibly inserted into the interlayer gallery space of layered hosts to form new 2-D heterostructured materials without any structural reconstruction.8–13) As widely studied, guest species including atoms, molecules, clusters and their charged species, regardless of whether they are inorganic or organic or bio species, have been successfully intercalated into layered ceramic materials to form various inorganic/ceramic, organic/ceramic and bio/ceramic 2-D nanohybrids with heterostructures. More recently, such hybrids have attracted considerable attention due to their multiplex potential for application in various fields such as energy and the environment, and in medical healthcare industries.14–21)
In general, the physico-chemical properties of 2-D nanohybrids are strongly related to their heterostructures; therefore, X-ray diffraction (XRD) analyses are required to solve their crystal structures in order to predict their structure-property relationships. However, when they are poor in terms of crystallinity due to their porous and/or amorphous nature resulting from lattice strains upon intercalation, such as structural defects and stacking faults, their structures can be analyzed not only by conventional XRD and small-angle X-ray scattering analyses but also by X-ray absorption spectroscopy (XAS), nuclear magnetic resonance spectroscopy (NMR), molecular dynamics simulations, and electron microscopic tools such as SEM and TEM.6,22–26)
To characterize such intercalative nanohybrids, powder XRD is the most common means of solving their structures along with one-dimensional (1-D) electron density map calculations, as well-ordered intercalative nanohybrids generally show a series of well-ordered (00l) XRD peaks corresponding to the basal spacing (d), of which the values increase upon intercalation.27) If a large number of (00l) XRD peaks were available for an intercalated nanohybrid, its 1-D electron density map along the crystallographic c-axis could be determined via the Fourier transform of the calculated structure factor, which reflects the intracrystal-line structure, such as the molecular orientations and arrangements, of the guest species in the interlayer space.
In this review, X-ray characterizations of intercalative nanohybrids are highlighted with some examples of layered hosts, such as metal oxides, cationic clays and LDHs as reported thus far, on the basis of the powder X-ray diffraction method with the 1-D electron density mapping technique. In addition, an attempt is also made to discuss experimental results with theoretical principles in detail in order to assist with an understanding of the structure-property relationships of such intercalative nanohybrids.
2. Structures of Layered Ceramic Materials
To predict the physico-chemical properties of 2-D nanohybrids with a heterostructure, their crystal structures should be well defined immediately after the synthesis process.
As well-known ceramic raw materials, smectites such as montmorillonite (MMT), saponite and hectorite are typical layered compounds belonging to cationic clays of the 2: 1 (the ratio of tetrahedral sheets to octahedral sheets) type. According to the crystal structure of smectite, infinite two-dimensional sheets are sequentially stacked along the c-axis to form a multilayer structure, where each sheet consists of a AlO6- or MgO6-octahedral (Oh) layer covalently sandwiched between two SiO4-tetrahedral (Td) layers (Fig. 1(a)).6,23,27) If Si4+ in the tetrahedral sheets and/or Al3+ or Mg2+ in the octahedral sheets are partially substituted by lower valent Al3+ or Fe3+, and/or Mg2+ or Fe2+ or Li+, respectively, a negative layer charge can be generated. In order to compensate for the generated layer charge, interlayer cations such as hydrated Na+ and Ca2+ ions must be stabilized between the two-dimensional clay lattice. Such interlayer cations are, however, readily exchangeable with other cationic species, regardless of whether they are inorganic, organic, or bioactive molecules.

Crystal structures of various layered materials; (a) montmorillonite, (b) KTiNbO5, (c) Bi2Sr2CaCu2O8, (d) LDH and (e) MoS2.
On the other hand, anionic clays, known as LDHs with anion-exchange properties, are often denominated as hydrotalcite- like materials due to similarities in terms of the crystal structure and layer charge of the hydrotalcite mineral (Mg6Al2(OH)16CO3·4H2O) in nature. The chemical formula of LDH is described as [M2+1−xM3+x(OH)2]x+(Am−)x/m·nH2O, where M2+ is a divalent metal, M3+ is a trivalent metal, and Am− is an interlayer anion.6,7,22,27) To understand the crystal structure of LDH in detail, it is necessary to explain the basic structure of brucite, Mg(OH)2, in which divalent Mg2+ ions are coordinated with six hydroxo ligands to form octahedra, which are then bound to one another by sharing edges to form an infinite sheet. In brucite, these sheets are stacked on top of each other to build up the layered structure through hydrogen bonding. In the LDH structure, however, divalent cations such as Mg2+, Zn2+, Ni2+, Ca2+, and Fe2+ are partially substituted by higher valent cations such as Al3+, Cr3+and Fe3+, resulting in a positive layer charge in the lattices. To compensate for the positive layer charge, solvated anionic species are simultaneously located in the interlayer space of the LDH lattice, which can be exchanged with other anionic molecules, as shown in Fig. 1(d).28)
Among the various types of layered metal oxides, the crystal structures of layered titanates11,14,29) and Bi-based cuprate high-Tc superconductors30,31) are briefly introduced in this section. A series of alkali metal titanates with a lepidocrocite (g-FeOOH)-related structure with the general formula of AxTi2−yMyOz (A = Na, K, Cs, Rb and M = Mg, Co, Ni, Fe, Ti, Nb, etc.) has a layer structure built up by sharing the edges and corners of TiO6 /MO6 octahedra, as shown in Fig. 1(b). The alkali metal ions located between the layers are exchangeable with other cations such as protons, other metal ions or organic cations. Their protonic forms are known to be expandable and, in consequence, easily exfoliated by intercalating bulky cationic molecules into metal titanate nanosheets.29,32–34) These exfoliated layers can further be used as building blocks for preparing heterostructured nanohybrids based on the lattice engineering strategy.
Bismuth cuprate high-Tc superconductors are also an example of a layered metal oxide, with the general formula of Bi2Sr2Can−1CunOy; they are often denoted as Bi2201, Bi2212 and Bi2223, where n equals 1, 2 or 3 in this case, respectively.30,31) For example, a single sheet of Bi2212 is composed of BiO, SrO, CuO2, Ca, CuO2, SrO and BiO layers, and the sheets are sequentially stacked on the top of each other to form a multilayered structure via weak van der Waals interaction, as demonstrated in Fig. 1(c). It is, therefore, not surprising that such layered Bi-based cuprates are also applicable as intercalative host materials.
Transition-metal dichalcogenide with the formula of MX2 is an another example of a layered material with a covalent bonding characteristic, where M is a transition-metal atom (e.g., Mo, W) and X is a chalcogen (S, Se, or Te).5,6) In the MX2 structure, the transition-metal atom is covalently bonded with six chalcogen ligands to form a prismatic unit, and the prismatic units formed in this manner are bounded to one another by the sharing of edges to form an infinite sheet (X-M-X). Specifically, these layers consist of a hexagonal plane of transition-metal atoms sandwiched between two hexagonal chalcogen planes (Fig. 1(e) that are stacked on top of each other to form a lamellar structure via weak Van der Waals interaction.
3. Computation and Examples of the 1-D Electron Density of 2-D Nanohybrids with a Heterostructure
XRD is a basic and powerful tool for characterizing the crystal structures of various solids. Therefore, the crystal structure of polycrystalline samples with good crystallinity can be solved via Rietveld refinement analysis (the entire XRD pattern fitting procedure to refine the structure).35–38) However, 2-D nanohybrids with a heterostructure, such as intercalation compounds, are not good enough in terms of the X-ray crystallinity for a Rietveld refinement analysis, as their XRD peaks are often broad and asymmetric due to the lattice strain upon intercalation, such as structural defects and stacking faults.39–43) Therefore, the collected XRD data are often not of a good quality; i.e., the (hkl) peaks of 2-D nanohybrids are poor as a whole, but a few sharp and strong (00l) peaks are frequently observed due to the effect of the preferred orientation of the 2-D particles. This is likely why only a few papers on Rietveld refinement have been reported thus far for layered materials intercalated with small organic and inorganic molecules, such as dimethylsulfoxide (DMSO), Cl−, and CO32−.44–47)
However, if a proper number of well-developed (00l) peaks were available for 2-D nanohybrids with a heterostructure, the structural arrangement of guest molecules in the interlayer space could be determined via Fourier transformation of the structure factors calculated from a series of (00l) peak intensities, resulting in a 1-D electron density map along the c-axis.12,42,48–53) If there were various possible interlayer structural models for 2-D nanohybrids, one could determine a best fit structural model as physically and chemically meaningful for 2-D nanohybrids by comparing the calculated electron density profile with that evaluated from the experimental XRD pattern.
3.1. Evaluation of 1-D Electron Density Map
The 1-D electron density profile can be calculated with Eq. (1) from the positons and the intensities of a large number of (00l) XRD peaks for intercalative nanohybrids.54–58)
This equation for centrosymmetric structures can be transformed to the following one (Eq. (2)),
where l, F(00l), c and z are the order of the (00l) diffraction peak, the structure factor, the unit cell parameter of the c-axis, and the atomic coordinate along the z-axis, respectively.
The structure factors corresponding to the expected structural model can also be calculated from the atomic scattering factors and positions in the structural model. The structural factor for centrosymmetric materials is calculated using the following equation (Eq. (3)),
where fi is the atomic scattering factor for each atom in the unit cell. For a practical calculation, the atomic scattering factors for x-ray as functions of the Bragg angle are represented by the following equation (Eq. (4)):
Here, the scattering factors of various neutral atoms and ions can be expressed as functions of nine coefficients for each element (c0, a1 – a4, b1 –b4) and sin θ/l, which are also found in the literature such as International Tables for Crystallography.59–60) The temperature factor, B, depends on the nature of the material. The typical range of the temperature factor (B) is ~0.5 to ~1 Å2 for ceramic ionic crystals and intermetallic compounds, and from ~1 to ~3 Å2 for other inorganic and many coordination compounds, while this range extends from ~3 to ~10 Å2 or higher for organic and organometallic compounds and for solvents or other intercalated non-bonded molecules or atoms.61) The atomic positions in the host layer can be obtained from the crystal structure of pristine layered materials before intercalation, and the initial interlayer structural model of guest species is established based on the Patterson diagram, the energy-minimized chemical structure, and the results from other analyses, including NMR, XAFS and TEM analyses.31,51,52)
The experimentally observed structure factor, Fo(00l), can be derived from the intensities of the (00l) XRD peaks, which are obtained from the integrated peak areas. Here, Fo(00l) should be corrected with the Lorentz-Polarization (LP) factor, expressed as Eq. (5), as the absorption factor can generally be negligible when the intensities are calculated in Debye-Scherrer method.58,61)
where
The sign (phase) of the observed structure factor can be obtained from the calculated result from the predicted structural model, as they should be identical. With Eq. (2), the electron density profiles can be plotted from z/c = −0.5 to 0.5 for 2-D nanohybrids.
From these computational processes, the observed 1-D electron density (ro) and the calculated one (rc) can be derived from the experimental XRD intensity (the integral peak area) and the predicted structural model, respectively. It then becomes possible to refine the reasonable orientation and arrangement of guest molecules in the interlayer by comparing these profiles. The reliability factor (R) can be calculated with the following equation (Eq. (6)): 57,62)
3.2. Examples of 1-D Electron Density Mapping for Inorganic/Inorganic 2D Nanohybrids
There are several reports on 1-D electron density mapping analyses of inorganic-inorganic 2-D nanohybrids such as the HgBr2-Bi2212 nanohybrid,48) the AgI-Bi2212 nanohybrid,31,63) and the Ni(OH)2-clay nanohybrid62) which can be used to confirm the structural model of intercalated species established by NMR, TEM and XAFS analyses. For example, Choy et al. reported a inorganic/inorganic 2-D heterostructured nanohybrid, the AgI-intercalated Bi2212 nanohybrid (Ag1.17I1.54-Bi2212), which demonstrated a superionic conductivity as well as a superconducting property.31,63) From the EXAFS and Raman analyses, the structural model of the Ag-I sub-lattice was proposed considering the available interlayer space and the volume required for a two-dimensional array of AgI4 tetrahedra. The expanded basal spacing for Ag-I-intercalated Bi2212 compared to that of iodine-intercalated Bi2212 demonstrated that an iodine double layer was present in the interlayer space and that the Ag-I sub-lattice should consist of a monolayer of AgI4 tetrahedra. Fig. 2 demonstrates the two most probable lattice types for AgI4 tetrahedra together with a schematic of the structural models upon intercalation, where one face (model A) or two edges (model B) are parallel to the AB plane of Bi2212. From the (Ag-I) bond distance determined by EXAFS analysis, the distance of I-I projected along the c-axis was estimated to be 3.75 Å and 3.24 Å for models A and B, respectively. As shown in Fig. 2(a), the height of the Ag-I layer in model A was consistent with the gallery height (~ 3.7 Å) evaluated by subtracting the basal spacing of iodine-intercalated Bi2212 (~ 3.6 Å) from that of AgI-intercalated Bi2212 (~7.3 Å), indicating the suitability of model A. To confirm the reliability of this structure, 1-D electron density profiles along the c-axis were calculated in the following two ways. First, the c-axial electron density profiles [r(z/c)] of AgI-Bi2212 can be obtained from the summation of structure factors as determined from the (00l) XRD peak intensities, where the magnitude of the structure factor, |F00l|, is the square root of the intensity corrected with the Lorentz polarization factor, as previous indicated in Eq. (5). The phase problem can be solved by assuming the conservation of intra-crystalline structure of host blocks during the intercalation process and by adopting the sign of calculated structure factor. Alternatively, |F00l| and r(z/c) can be calculated theoretically for the given intra-crystalline structural models of the Ag-I sub-lattice between Bi2212 layers. The 1-D electron density profiles simulated on the basis of structural models A and B are compared in Fig. 2 (a) and (b), respectively, with that obtained from the experimental XRD data. As expected from the two different layer stacking sequences of •••/BiO/I/Ag/Ag/I/BiO/••• (model A) and •••/BiO/I/Ag/I/BiO/••• (model B), the simulated patterns were found to be significantly distinguished in the region of the guest layer. The calculated (z/c) profile with model A was in good agreement with the experimental case, which clarified the reliability of model A, where the reliability factor was determined to be less than 0.05. In the r(z/c) profiles based on XRD and on structural model A, the peaks at z/c = ± 0.08 are attributed to the presence of Ag and I. Because both ions were located at highly adjacent positions on the c-axis, their individual contributions could not be separated but be observed as a single component, which provided an explicit explanation of the double-layered HRTEM pattern of the interlayer Ag-I lattice. As shown in Fig. 2(c), the dark spots in the magnified HRTEM image demonstrate good coincidence with peaks in the electron density profile evaluated from model A. This Ag1.17I1.54-Bi2212 nanohybrid exhibited fast ionic conductivity as high as 10−2 – 10−1 Ω−1cm−1 in the range of 80 - 270°C and high-temperature superconductivity as well with a TC value of 67K, which is thought to be slightly lower than that of the pristine Bi2212 (80 K) due to the charge transfer between the host and the guest.

Proposed structures of (a) model A and (b) model B for the interlayer AgI lattice between Bi2212 layers, together with their corresponding 1-D electron density profiles, where the calculated values (solid line) are compared with the experimentally observed values (dotted line), and (c) HRTEM image of the Ag1.17I1.54-Bi2212 nanohybrid together with the 1-D electron density profile calculated based on the crystal structural model A. Reproduced and modified with permission from reference 63. Copyright 2000 American Chemical Society.
In the case of mercury bromide intercalated Bi2212, (HgBr2)0.5Bi2Sr2CaCu2Oy, an unusual heterostructure consisting of a superconducting layer and an insulating one was suggested.48) The interlayer structure of HgBr2 inbetween superconducting layers (Bi2212) was determined by the EXAFS fitting and XRD results. In order to confirm the reliability of this structure, 1-D electron density profiles along the c-axis obtained in two different ways were compared; one was calculated from the (00l) XRD peak intensities for the well-oriented HgBr2-intercalated Bi2212 single crystal, and the other was calculated from the structure factor based on the present structural model, including the atomic positions and their occupation characteristics. According to Fig. 3, there was excellent accordance between the calculated and experimental 1-D electron density profiles, indicating the reliability of the present structural model in which linear HgBr2 molecules were present between the (Bi-O) double layers, as shown in Fig. 3. From this result, it was concluded that the HgBr2 linear molecule was tilted by 65° with respect to the c-axis.

1-D electron density profiles along the c-axis of (HgBr2)0.5Bi2Sr2CaCu2Oy with its structural model and the atomic z/c coordinates used. The dotted and solid lines represent 1-D electron density profile calculated from the intensities of the (00l) XRD peaks of HgBr2-intercalated Bi2212 nanohybrid and that calculated from the present structural model, respectively. Reproduced and modified with permission from reference 48. Copyright 1997 American Chemical Society.
Uehara et al. reported another example of 1-D electron density mapping for inorganic-inorganic 2-D nanohybrids.62) They were successful in intercalating Ni(OH)2 sheets into clay layers and suggested a model with a monolayer orientation of Ni(OH)2 simply by comparing the calculated and observed 1-D electron density profiles with a reasonable R value of 0.10. Argüelles et al. also suggested the positions of interlayer water molecules in vermiculite through 1-D electron density profile fitting.64)
3.3. Examples of 1-D Electron Density Mapping for Organic/Inorganic 2-D Nanohybrids
There have been several reports of 1-D electron density mapping for organic-inorganic 2-D nanohybrids, for example layered ceramic materials intercalated with organic molecules, in an effort to determine the orientation and arrangement of organic molecules stabilized in the interlayer space.65–68) For example, Marangoni et al. synthesized Zn2Al-LDH nanohybrids intercalated with a series of blue dye molecules such as Evans blue (EB), Chicago sky blue (CB) and Niagara blue (NB) via a coprecipitation method, where interlayer structures of dye molecules were proposed by 1-D electron density mapping. Due to the relatively large number of (00l) diffraction peaks (up to 7 (00l) peaks), they were able to suggest interlayer structures of those dye molecules projected on the c-axis via Fourier transform analysis. The 1-D electron density profiles for CB-Zn2Al, EB-Zn2Al and NB-Zn2Al LDH nanohybrids are displayed in Fig. 4, where the integral intensities of XRD peaks were determined by the peak profile fitting method with the Fullprof program.69) The 1-D electron density profiles for CB-Zn2Al and EB-Zn2Al LDH nanohybrids were similar, while that NB-Zn2Al LDH nanohybrid was different from the others. For all the dye-LDH nanohybrids, the two strong peaks at z = 0 Å and 20.4 Å/20.3 Å/23.4 Å were attributed to the LDH layers, which contained metal cations, and the other lower electron density arose from the interlayer dye molecules. In all of the nanohybrids, relatively strong peaks at a distance of 4.0 Å from the center of LDH layer were attributed to the sulfonate groups, as there were hydrogen bonding interactions between sulfonate groups and OH ones in LDH layers (-S-O···H-O-Metal). Hence, they proposed that the dye molecules such as CB and EB were intercalated into LDH layers with a tilting angle of 55–60° with respect to the basal plane of LDH by considering the structural dimensions of dye molecules with a minimized level of steric energy. The methoxy groups in CB and EB provided medium electron density in the middle range of the interlayer of LDH. However, 1-D electron density of NB-Zn2Al LDH nanohybrid was minimized at the center, which is different from those of other dye-LDH nanohybrids, indicating that NB molecules are most likely oriented with an arrangement of parallel bilayer inbetween the LDH layers due to the positions of sulfonate groups in NB; that is, all sulfonate groups in NB were not present at the opposite ends of the molecule, as in CB and EB, but on the same side of the molecule. It is, therefore, highly expected to have such a bilayer arrangement resulting in a surface area of ~ 17 Å2 per unit charge, which is sufficiently small comparable to the equivalent area of the host lattice (24.8 Å2e−1 for Zn2Al LDH).
As an another example of 1-D electron density mapping for an organic/inorganic 2-D nanohybrid, Fujita et al. proposed an appropriate molecular orientation of rhodamine 6G (R6G) in the interlayer space of Li-fluortaeniolite,66) a smectite-type clay, by calculating 1-D electron density along the z-direction from ten well-ordered (00l) XRD peaks as well as from the possible model structures. Based on the expansion of basal spacing and the molecular dimension of R6G, they proposed two possible orientation models of R6G in the interlayer for Fourier transformation; one is a vertical model with a perpendicular monolayer orientation where the longest axis of xanthene ring in R6G is perpendicular to the basal plane of clay, and the other is a parallel model where the longest axis of xanthene ring is parallel to the basal plane, that is, R6G molecules are oriented with a bilayer stacking arrangement. The 1-D electron density profile calculated from the vertical model and the difference Fourier synthesis result are presented in Fig. 5 together with the proposed vertical model, including the positional parameters for R6G between the clay layers. As can be seen clearly in Fig. 5, the values from the difference Fourier synthesis along the c-axis are very small, indicating that the 1-D electron density profile calculated from this vertical structural model is in excellent agreement with that calculated from the intensities of the (00l) XRD peaks. The R-factor for the vertical model was determined to be ~0.01, while that for the parallel model was ~0.20. Based on this 1-D electron density mapping analysis, they suggested the structural model of the R6G intercalated clay, where R6G was intercalated with a monolayer perpendicular orientation to the ab plane of the clay layer, and the positively charged N atoms of R6G were bonded to the negatively charged clay layers by strong electrostatic interactions.

(a) The orientation of rhodamine 6G (R6G) in clay layers and its 1-D electron density profile along the z-direction, and (b) positional parameters of the R6G intercalated clay nanohybrid. From [Fujita et al. (1997, Fig. 8 and Table 4)] (reference 66). Reproduced and modified with kind permission of The Clay Minerals Society, publisher of Clays and Clay Minerals. Copyright 1997 The Clay Minerals Society.
A similar study of 1-D electron density mapping was also carried out on the PEO-V2O5 nanohybrid by Liu et al.55) They successfully intercalated poly(ethylene oxide) (PEO) into a layered V2O5 xerogel simply by mixing an aqueous PEO solution with an aqueous V2O5 solution, after which they slowly evaporated the water molecules. According to the 1-D electron density mapping results, it was confirmed that the intercalated PEO chains were not in the form of a helical coil structure but with a zigzag conformation; that is, the PEO chains were fully extended and arranged horizontally with a bilayer conformation stacked on top of each other with respect to the (002) plane. The interaction between PEO and V2O5 in this organic-inorganic nanohybrid appeared to be the type of van der Waals and hydrogen bonding. As a similar example of organic/inorganic nanohybrid reported by Adams et al., tetrahydrofuran (THF) was intercalated into MMTs containing various cations such as NH4+, Na+, Sr2+, Co2+, Ni2+, and Cu2+.67,68) According to the 1-D electron density calculation, it was found that THF molecules were oriented perpendicularly to the basal plane of clay layers. Kanamaru et al. also proposed a parallel bi-layer orientation of 6-amino hexanoic acid in vermiculite through an 1-D electron density mapping analysis.56) Bauer et al. also suggested the possible interlayer orientation of perylene dye, N,N′-di(phenyl-3,5- disulfonic acid)perylene-3,4:9,10-tetracarboxyldiimide (PBITS) stabilized inbetween LDH layers based on 1-D electron density mapping analysis.70) Lerner et al. was also successful in preparing the bis(oxalato)borate anion (BOB) intercalated graphite in the presence of a solution containing BOB anion molecules in anhydrous hydrofluoric acid, and in mapping its 1-D electron density, and eventually in proposing the orientation of BOB anions in the basal plane of graphite.71) For porphyrin-LDH 2-D nanohybrids,72) porphyrin molecules such as Pd(II)-5,10,15,20-tetrakis(4- carboxyphenyl) porphyrin (PdTPPC) and Zn-5,10,15,20- tetrakis(4-sulfonateophenyl)porphyrin (ZnTPPS) were intercalated into LDH to form 2-D organic-inorganic nanohybrids. The interlayer orientation of porphyrin molecules was also confirmed to be perpendicular to the LDH layers on the basis of 1-D electron density mapping.
3.4. Examples of 1-D Electron Density Mapping for Bio-organic/Inorganic 2-D Nanohybrids
It has been reported that a large number of bio- and bioactive molecules, such as RNA, DNA, anti-oxidant molecules, and drug molecules, can be stabilized in layered materials such as clays and LDHs to form bio-inorganic 2-D nanohybrids. In this section, several examples of 1-D electron density mapping results of drug-LDH nanohybrids are proposed to confirm the interlayer structure of drug molecules between LDH layers.42,43,73) Recently, Kang et al. reported a ferulic acid (FA) intercalated LDH nanohybrid with sustained release property. They proposed two possible interlayer structural models of FA on the basis of XRD and FT-IR analyses and the molecular dimension of FA. In order to confirm the more reliable interlayer structure of FA in the LDH layers, they computed the 1-D electron density from the integral intensities of the (00l) XRD peaks and the proposed structural models shown in Fig. 6; the model-1 was based on the perpendicular packing of FA in a zig-zag manner with alternating phenyl–acrylate π–π interaction (Fig. 6(a)), whereas the model-2 was made under an assumption of tilted packing (tilting angle ~ 61°) of FA moiety to maximize the π–π interaction between phenyl groups (Fig. 6(b)). The 1-D electron density profile computed from the observed (00l) XRD peaks was compared with those computed from the two structural models. As shown in Fig. 6, the 1-D electron density profile computed from the observed (00l) XRD peaks was better fitted with that of the model-1 than that of the model-2, where the reliability factor for the former was 0.10, while that for the latter was 0.19. This result revealed that FA molecules were intercalated between the LDH layers with a perpendicular zig-zag packing orientation rather than with a tilted orientation of π–π packing between phenyl rings. This FA-LDH nanohybrid showed a controlled release behavior in a simulated body fluid medium.

Two possible interlayer structural models of ferulic acid (FA) intercalated LDH nanohybrid and their corresponding 1-D electron density profiles: (a) Model-1: perpendicular zig-zag packing of FA molecules with an alternating phenol–acrylate π–π arrangement. (b) Model-2: tilted packing of FA molecules with maximized π–π interaction between phenyl groups. The solid (Dobs) and dashed (Dcal) lines represent the 1-D electron density profiles from the XRD intensity and structural model, respectively. Reproduced and modified with permission from reference 73. Copyright 2008 Elsevier Ltd.
For the flurbiprofen-LDH nanohybrid, flurbiprofen (FB), a non-steroidal anti-inflammatory drug, was intercalated into Zn2Al-LDH as a new transdermal drug delivery system, and its 1-D electron density was calculated in order to understand the interlayer arrangement of FB drug molecules between the LDH layers. The integral intensities of six well-developed (00l) peaks were used for Fourier transformation, and the positions of Zn, Al, and OH in the LDH layer were obtained from the crystal structure of Zn2Al-CO3 LDH. It was found that the rational structure of interlayer FB molecules was an interdigitated bilayer structure oriented with a tilting angle of ~21° by comparing the calculated and observed electron density profiles, those which were made on the basis of different interlayer orientation models with a tilting angle of the interdigitated bi-layer ranging from 0° to 31.6°. A strong peak at the z-parameter of 0 Å in the 1-D electron density map corresponds to Zn and Al in the LDH layer, and six additional peaks with lower electron densities could be attributed to the intercalated FB and water molecules in the interlayer. Another example of a bio/inorganic 2-D nanohybrid was the artesunate (AS), anti-malarial drug, intercalated layered zinc hydroxide (ZBS).43) The rational structure of the interlayer AS molecules in the AS-ZBS nanohybrid was also confirmed by 1-D electron density mapping analysis together with TEM analysis.
Recently, a time-resolved in situ energy dispersive X-ray diffraction (EDXRD) experiment was developed using a synchrotron X-ray device.26,74–76) This is a very useful tool for investigating the intercalation mechanism and its kinetics, since the intercalation reactions can be simultaneously traced as soon as the guest species is added to the host suspension. From this time-resolved in situ EDXRD, Taviot- Gueho et al. clarified the second-stage phenomena (the interstratified structure) during the intercalation of tartrate into Zn2Cr-Cl LDH.26) In addition, Williams et al. demonstrated the staging structure during phosphonic acid intercalation into LiAl2-Cl LDH via this time-resolved EDXRD analysis.75) This in-situ measurement would be very useful to elucidate the interlayer structures of 2-D nanohybrids together with 1-D electron density mapping analysis.
4. Summary
Various types of guest species, such as molecules with cationic charge or anionic one, metal ions, metal oxide clusters, polymers, and even 2D oxide nanosheets have been incorporated into the inter-gallery of layered ceramic materials such as clays, LDHs, graphite, layered metal oxides and layered transition metal chalcogenides to form 2-D nanohybrids with various heterostructures. It is always an important issue to understand the interlayer structure of guest species between host ceramic layers, not only to trace intercalation reactions but also to understand their physicochemical properties in terms of the structure-property relationships. In this regard, the 1-D electron density mapping analysis is very useful tool to elucidate the interlayer structure in 2-D nanohybrids together with x-ray spectroscopy (XANES and EXAFS), NMR and TEM analyses. We hope that some of the highlighted examples in this review will provide useful information to readers who are currently working on new 2-D nanohybrids based on ceramic lamellar materials such as inorganic/ceramic-, organic/ceramic- and bio/ceramic-2D nanohybrids with various heterostructures.

1-D electron density profiles along the z-direction for flurbiprofen stabilized between the Zn2Al-LDH layers and its corresponding structural model, where the solid and dashed lines indicate the observed and calculated electron density profiles, respectively. Reproduced and modified with permission from the reference 42. Copyright 2013 Elsevier Ltd.
Acknowledgments
This work was supported by a grant from the National Research Foundation of Korea (NRF) funded by the Korean government (MSIP) (2005-0049412 and NRF-2013R1A1A 2062239).