Open Access Article
Jingxiang
Cheng
ab,
Malsha
Udayakantha
ab,
Saul
Perez-Beltran
ab,
Luis
Carrillo
ab,
Wasif
Zaheer
ab,
Lucia
Zuin
c and
Sarbajit
Banerjee
*ab
aDepartment of Chemistry, Texas A&M University, College Station, TX 77843-3012, USA. E-mail: banerjee@chem.tamu.edu
bDepartment of Material Science and Engineering, Texas A&M University, College Station, TX 77843-3012, USA
cCanadian Light Source, University of Saskatchewan, Saskatoon, SK S7N 2V3, Canada
First published on 29th August 2024
Increasing interest in halide-ion batteries as a means of long-duration energy storage has led to intense recent interest in mechanisms of anion migration in periodic solids. Lanthanide oxyhalides crystallize in a richly diverse set of structures with well-separated halide-ion slabs and represent potential solid electrolytes of halide batteries. However, oxyhalides of later lanthanides are underexplored with regard to preparatory methods, structural preferences, defect energetics, and anion migration pathways. Here, we report a non-hydrolytic sol–gel condensation route to ligand-passivated single-crystalline ErOCl nanoplatelets, which crystallize in van der Waals layered SmSI and YOF structure types with reduced coordination of the lanthanide centers as compared to the PbFCl structure type preferred by early lanthanide oxyhalides. Pressure-less annealing in an argon ambient environment results in anisotropic sintering of the platelets with strongly preferred grain growth along the ab plane whilst retaining the thin platelet morphology. The preferred grain growth can be traced to strongly anisotropic chloride-ion migration along the ab plane of the bilayered van der Waals solid, which represents the primary mechanism of mass transport across grain boundaries. First-principles nudged elastic band simulations indicate that a low-energy intralayer chloride-ion migration pathway becomes accessible within the chloride-ion slab along the ab plane in the presence of chloride vacancies. The simulations demonstrate the promise of modulation anion conductivity through site-selective modification on the cation and anion lattices. These findings offer valuable insights into an important emerging class of halide electrolytes, illuminate fundamental mechanisms governing anion transport in 2D materials, and provide a blueprint for materials processing to achieve directional ion diffusion pathways.
Lanthanide oxyhalides exhibit a diverse array of crystal structures characterized by well-separated halide-ion slabs,16 which makes them promising candidates for solid electrolytes of halide-ion batteries. This promise is manifested in record halide ion diffusion coefficients recorded for lanthanide oxyiodides, oxybromides, and oxychlorides.17–21 The well-separated halide slabs of these compounds are amenable to vacancy engineering, which enables facile anion transport. In addition to their capability for halide-ion conduction, lanthanide oxyhalides demonstrate excellent thermal stability, which makes them particularly well suited for high-temperature applications.17,18,20,21 The large bandgap and low electronic conductivity of lanthanide oxyhalides are furthermore optimal for employment as solid electrolytes of halide batteries where it is imperative to separate ion and electron flows.3,22,23 Such solid electrolytes can be paired with insertion or conversion electrodes with transition metal redox centers.3,24
Aliovalent alloying on the lanthanide cation sublattice is used to engender point defect formation and facilitate anion-plane conduction.20 Recent work has demonstrated facile and versatile halide-ion mobility in the anion slabs of PbFCl-structured LaOI wherein topochemical halide-ion replacement is driven along a gradient of increasing hardness of anions from LaOI to LaOF whilst entirely preserving the tetragonal PbFCl crystal structure.25
In contrast to the tetragonal matlockite crystal structure (P4/nmm-D74h) preferred by early lanthanides (La–Ho), the later lanthanide (Er–Lu) oxyhalides tend to adopt a layered rhombohedral R
m–D53d structure.16,26,27 These structures have alternating hexagonally packed (LnO)+ bilayers separated by hexagonally packed X− layers (amenable even to cation insertion in their quasi van der Waals gap).28 Mechanisms of halide-ion diffusion in these open-framework structures remain almost entirely unexplored. In this work, we report the synthesis of ligand-passivated ErOCl nanoplatelets through a non-hydrolytic sol–gel condensation reaction, examine mechanisms of chloride-ion diffusion, and demonstrate a surprising anisotropic sintering mechanism mediated by constrained chloride-ion diffusion along two-dimensional anion planes.
ErOCl is a dimorphic compound and crystallizes both in tetragonal and rhombohedral polymorphs. The rhombohedral structure, layered along the c-axis, features two polymorphs with different stacking sequences: ABC stacking corresponds to the SmSI structure type, whereas ACB stacking corresponds to the YOF structure type. Synthetic outcomes are typically a mix of both stacking sequences. In general, the synthesis of oxyhalides of later lanthanides is relatively poorly explored. Hydrochlorination of lanthanide oxides followed by hydrolysis yielded poorly crystalline mixtures of rhombohedral (YOF)/tetragonal27 or SmSI/YOF structure types.28 Annealing these mixtures at temperatures >750 °C for extended periods of time (2–7 weeks) using halide fluxes (LiCl/KCl/ErCl3) to dissolve impurity phases has been reported to yield more crystalline materials. However, the products continued to be contaminated by Er2O3 and other unidentified phases.27,29 Attempts to synthesize ErOCl by precipitation of ethanol dispersions of ErCl3 with NH4OH and subsequent calcination at 500 °C yielded amorphous products even though this synthesis provided a reliable path to obtaining oxyhalides of the early lanthanides (La–Ho) crystallized in the PbFCl structure.30 ErOCl phases obtained with a similar approach where hydrated ErCl3 was sintered started to oxidize to Er2O3 at 550 °C,31 precluding further crystallization of SmSI/YOF-structured phases.
In this article, we devise a facile non-hydrolytic ligand-exchange and condensation route to prepare ErOCl nanocrystals crystallizing in the rhombohedral R
m space group without additional annealing. ErOCl nanoplatelets prepared by this method adopt a rhombohedral 2D layered R
m crystal structure with SmSI and YOF structure types. The prepared nanoplatelets are characterized by powder X-ray diffraction, transmission electron microscopy, Fourier transform infrared spectroscopy, and X-ray absorption near-edge structure (XANES) spectroscopy. High-temperature annealing of ErOCl results in anisotropic sintering into larger thin nanoplatelets, which we posit is mediated by 2D chloride-ion diffusion. Fundamental mechanisms of chloride-ion diffusion in these structures have been examined by first-principles density functional theory (DFT) calculations.
The nanoplatelets were sintered using an MTI GSL-1700 high-pressure tube furnace. As-prepared powders were placed in an alumina crucible, gently pressed to maintain contact, and then heated at 1050 °C for a duration of 12 h at a controlled heating rate of 5 min °C−1 under a consistent Ar flow of 70 mL min−1.
FTIR spectroscopy was performed using a Bruker Vertex-70 instrument with a Pike MIRacle single-reflection total attenuated reflectance (ATR) accessory.
Transmission electron microscopy (TEM) images were acquired using an FEI Tecnai G2 F20 ST instrument operated at an accelerating voltage of 200 kV. The ErOCl powder was dispersed in 2-propanol and then the solution was deposited onto 300-mesh carbon-coated Cu grids for high-resolution TEM (HRTEM).
Er N4,5 X-ray absorption near-edge fine structure (XANES) spectra were concurrently acquired at the Variable Line Spacing-Plane Grating Monochromator (VLS-PGM) beamline of the Canadian Light Source with an energy resolution E/ΔE >10
000.37 Powder samples were adhered to carbon tape and all spectra were recorded at room temperature, with a step size of 0.1 eV and a dwell time of 1 s in a sample chamber maintained below 1 × 10−8 Torr. In collecting XANES spectra, the total fluorescence yield (FLY)38 was recorded with an MCP detector. All spectra recorded were normalized with respect to the incident photon flux (I0), which in turn was recorded by monitoring the current emitted by a Ni mesh (90% transmission) located just upstream of the sample.
N K-edge XANES spectra were collected at the National Synchrotron Light Source II of Brookhaven National Laboratory beamline SST-1 operated by the National Institute of Standards and Technology. Measurements were performed in partial electron yield (PEY) mode with a nominal resolution of 0.1 eV. The PEY signal was normalized to the incident beam intensity of a clean gold grid to eliminate the effects of any incident beam fluctuations and optical absorption features.
Vacancy formation energies were determined by calculating the difference in total energy between the pristine structure and a modified structure with a vacancy. This vacancy was created by displacing the atom 3 Å from the topmost layer on the slab to avoid charge imbalance.45–47 The vacancy formation energy EF was calculated using the total energy per atom (n) from the pristine structure EP subtracted by the total energy per atom (n) of the vacant structure EV: EF = (EP − EF)/n.
A single-point energy calculation was performed to calculate the projected density of state (pDOS) with Local Orbital Suite Toward Electronic-Structure Reconstruction (LOBSTER).48,49 Bunge's description for the local basis functions was used for the pDOS calculation with 6s and 4f orbitals for erbium, 2s and 2p orbitals for oxygen, and 3s and 3p orbitals for chlorine.
The climbing-image nudged elastic-band (CI-NEB) method using a 288-atom supercell (14.88 Å × 12.88 Å × 27.15 Å) was used to calculate the migration barrier along the minimum energy pathway.50 The same simulation parameters from geometry optimizations were applied except that the Brillouin zone integration grid was set to the Γ-point (1 × 1 × 1); the ionic optimization loop was used in conjunction with a force optimizer with a convergence criterion of 0.03 eV Å−1.51
m–D53d polymorphs being accessible.26,27 In contrast to the tetragonal structure with 3D bonding mediated by La–Cl connectivity between the slabs (Fig. 1C), rhombohedral ErOCl is layered along the crystallographic c-axis (Fig. 1D). The separation between halide ions in adjacent layers is ca. 3.45 Å (which yields an interlayer spacing of ca. 2.56 Å). Er atoms reside at the centers of mono-capped trigonal antiprisms.27,52 The rhombohedral structure adopts two polymorphs differentiated by distinctive stacking sequences along the crystallographic c-axis as depicted in Fig. 1D. The same Cl–Er–O–O–Er–Cl connectivity is manifested with the only distinction being the sequence of stacks: ABC stacking corresponds to the SmSI structure type, whereas ACB stacking emulates the YOF structure.
First-principles DFT calculations of the total density of states (DOS) are plotted in Fig. 1E along with the atom-projected partial DOS. Fig. 1F plots the orbital-projected DOS (pDOS) for both stacking sequences. In both polymorphs, the conduction band is derived primarily from narrow erbium 4f contributions, whereas the valence band is derived primarily from chloride 3p and oxygen 2p-derived states. As is typical of a dimensionally confined van der Waals solid, the electronic structure is primarily governed by intralayer bonding in the bilayers with the bilayer stacking sequence having scarce little effect on the Fermi surface. The relative invariancy of the electronic structure points to the origins of the synthetic challenges in controlling the stacking sequence in homogeneous bulk synthesis given that the two polymorphs are closely spaced in energy.
As noted above, synthetic methods for the preparation of well-defined nanocrystals of rhombohedral ErOCl are rather sparse. Here, we adapt a synthetic methodology (Fig. 2A) based on non-hydrolytic sol–gel condensation to prepare well-defined nanoplatelets of ErOCl.25,32 Our proposed reaction mechanism is analogous to non-hydrolytic sol–gel processes for metal oxides which proceed with the elimination of alkyl halide or dialkyl ether.34,35 Evidence for ligand scrambling as an initial step is derived from previous work on tetravalent alkoxides where the chain length and branching of the alkoxide precursor were systematically modified.33 An initial erbium haloalkoxide is thought to be stabilized as per:
| ErCl3(soln.) + Er(OiPr)3(soln.) → ErCl(3−x)(OiPr)x(soln.) + ErClx(OiPr)(3−x)(soln.) | (1) |
| ErCl3−x(OiPr)x(soln.) + ErClx(OiPr)3−x(soln.) → (OiPr)(Cl)–Er–O–Er–(Cl)(OiPr)(soln.) + iPrCl(soln.) | (2) |
| ErCl3(soln.) + (OiPr)(Cl)–Er–O–Er–(Cl)(OiPr)(soln.) → 2ErOCl(s) + 3iPrCl(soln.) | (3) |
An alternative E1 elimination, ligand scrambling, and condensation route has also been suggested, which in this case would involve ErCl(OH)(OiPr) and ErCl2(OH) intermediates.34,53
Notably, in comparison with previous syntheses of PbFCl-structured early lanthanide oxyhalides,25,54 a high molar excess of the amine coordinating solvent leads to stabilization of ErOCl nanocrystals. Similarly, trials with trioctylphosphine oxide (TOPO) as the coordinating ligand led to similar preferential stabilization of oxides and mostly amorphous products (see Fig. S1†). The use of long-chain amines such as oleylamine as the coordinating solvent and passivating ligand at low molar excess (25 mmol to 5 mmol amine per 2 mmol ErCl3) yields rhombohedral ErOCl.
Fig. 2B contrasts the powder XRD patterns of ErOCl synthesized with oleylamine (OAm), hexadecylamine (HDA), and tetradecylamine (TDA) as the capping ligands. In stark contrast to previous preparatory methods involving hydrochlorination of oxide followed by pyrohydrolysis, which yielded amorphous products with relatively poor long-range ordering,27 our synthetic approach yields crystalline materials. The powder XRD patterns have been indexed to a trigonal R
m unit cell (PDF# 49-1800) corresponding to (a = b = 3.7 Å; c = 27.6 Å) for mixed SmSI/YOF structure types. Correspondence to the expected layered SmSI and YOF-type phases is seen in all the measured powder XRD patterns, with substantial splitting of the (003) reflection suggestive of significant layer disorder. No other crystalline products are recovered.
The ligand-passivated ErOCl nanocrystals have been annealed under an Ar atmosphere.55 The powder XRD patterns of the annealed samples are shown in Fig. 2C. The relative intensities are modified selectively for the (003) and (110) reflections with respect to the other reflections. The strong growth of these reflections is suggestive of strongly anisotropic sintering. Fig. S2† displays simulated powder XRD patterns for differently constructed supercells for both YOF and SmSI structure types. Based on these simulations, the only powder XRD patterns that are consistent with the strong amplification of the (003) and (110) reflections are those that correspond to selective sintering along the ab plane, which results in extended X-ray coherent domains along this orientation. Fig. 2C shows that upon sintering, the XRD patterns include contributions from several by-products of the sintering process, such as graphite (C, PDF# 25-0284), carbon nitride (C3N4, PDF# 50-1512), and carbon chloride (C6Cl6, PDF# 02-0443), which are derived from sintering of the passivating ligands. Nevertheless, the reflections attributed to the SmSI phase are indicative of the preferred orientation along the c-axis, as corroborated by the simulated powder XRD patterns in Fig. S2.†
In Fig. 2D–F, prior to sintering, nanoplatelets are observed as individual platelets or agglomerates. The platelets adopt a mixture of hexagonal and rectangular morphologies. Notably, compared to prior studies on amine-capped early lanthanide oxychlorides,32 the amine-capped ErOCl nanoparticles exhibited a more pronounced preferred orientation along the ab plane.
The TEM images obtained before and after annealing corroborate a similar trend, with the particle size growing significantly along the ab plane (Fig. 3A–F). Particles of ErOCl are fused along different directions,56 resulting in faceted irregular-shaped particles that nevertheless retain their thin platelet morphology. The lattice-resolved HRTEM images of OAm-capped ErOCl before sintering (Fig. 3G and H) reveal distinct (110) and (003) planes, consistent with the simulated XRD patterns in Fig. S2,† which are indicative of anisotropic sintering along the ab plane as schematically illustrated in Fig. 3I.
In pressure-less sintering, the sintering rates along different crystallographic directions are dependent on both structural and geometric characteristics of the platelets.57 In general, direction-dependent diffusion processes along particle surfaces modify the curvature of the particle surfaces, resulting ultimately in evolving grain boundary orientations in polycrystalline intermediates.58 The directionality of sintering is governed by the grain boundary mass transfer or migration rates, which depend on the relative inclination and misorientation at grain boundaries. Given the strongly anisotropic crystal structure, mass transfer underpinning anisotropic sintering is likely to be governed by migration of the most mobile species, which in ErOCl are the chloride ions.20 Chloride-ion diffusion in the ab plane thus yields extended thin platelets instead of 3D sintered solids. Notably, sintering hexagonal platelets manifest strongly preferred orientation of particles as desired for anisotropic ion diffusion. For instance, Youn and co-workers demonstrated strongly anisotropic proton migration in {0001} preferentially oriented hexagonal Ba5Er2Al2ZrO13 perovskite-structured solid electrolytes.59 Controlling sintering directionality and anion conduction requires a detailed evaluation of chloride-ion diffusion mechanisms in these layered structures, which we describe in subsequent sections using first-principles simulations.
The Fourier transform infrared (FTIR) spectrum obtained for oleylamine-capped ErOCl nanocrystals is plotted alongside the spectrum measured for the oleylamine ligand in Fig. 4A. The broad band centered at 3300 cm−1 corresponds to symmetric and asymmetric stretching vibrations of the NH2 group. Bands at 2926 and 2852 cm−1 are ascribed to asymmetric stretching modes of –CH3 and –CH2– moieties. FTIR bands at 900–1000 and 650–900 cm−1 correspond to N–H bending and wagging modes, respectively.28 Moreover, –C–N– bending (1078 cm−1) and –C
C– bending (1073 cm−1) modes can also be discerned.52,55 Similarly, the FTIR spectra of hexadecylamine- and tetradecylamine-capped ErOCl demonstrate characteristic modes of HDA and TDA (Fig. S3†).
The Er N4,5-edge XANES spectrum in Fig. 4B corresponds to excitation of core 4d3/2 (N4-edge) and 4d5/2 (N5-edge) electrons to the 4f levels of Er3+. Preceding the giant resonance, three distinct features are assigned to multiplet splitting of 4d → 4f excitations. The N4,5 edge absorption spectra of heavier lanthanide cations60,61 are dominated by a combination of discrete 4d104fn → 4d94fn+1 transitions.62–64
The N K-edge XAS spectra illustrated in Fig. 4C unveil multiple absorption features that are attributed to ligand-derived states.65,66 The broad resonance in the vicinity of 405 eV can be attributed to NH4+ and C–N σ* interactions, reminiscent of spectral features measured for allylamine and propylamine.65,67 The region spanning from 398 to 401 eV reflects amine–erbium, C–N π*, and N–H σ* interactions.68,69 The observation of N–Er features, which are indicative of dative interactions involving nitrogen lone-pairs, provides insight into the mode of surface passivation. Specifically, the hybridization of nitrogen 2p lone-pair-derived slightly bonding LUMO states with Er states gives rise to both low-energy π* features and higher energy σ* features.
In order to evaluate the mechanisms governing anisotropic sintering, we next turn our attention to examining defect energetics and pathways for chloride-ion migration in ErOCl. We have evaluated chloride-ion migration considering three distinctive materials, an ordered ErOCl lattice, ErOCl with 10 at% Cl vacancies, and ErOCl with 10 at% Cl, 5 at% Er, and 3 at% O vacancies. Considering defect formation energies, Cl vacancies (0.042 eV per atom) have a substantially lower formation energy as compared to Er (0.111 eV per atom) and oxygen (0.078 eV per atom) vacancies, respectively. A strong preference for chloride vacancies over oxygen vacancies was previously evidenced for PbFCl-structured Eu-doped LaOCl in neutron activation analysis upon aliovalent substitution of divalent europium ions on trivalent lanthanum sites.70
Fig. 5A and B show that chloride-ion migration in non-defective ErOCl (rhombohedral SmSI structure type) proceeds within the same chloride plane. The diffusing anion adopts a triangular pyramid configuration with the closest Er atoms at the transition state; the closest Er ion is only 2.26 Å apart, whereas the other two adjacent Er ions are each at a distance of 3.41 Å. No significant structural deformation is observed beyond the first coordination sphere of the diffusing chloride ion. However, the migration barrier along this path is quite substantial, estimated to be 2.15 eV.
In contrast, randomly distributed 10 at% Cl vacancies in ErOCl afford access to an alternative lower-energy interlayer chloride-ion migration pathway (Fig. 5C and D). The transition state configuration for this migration path has the diffusing chloride ion interacting with multiple Er atoms from both layers at distances between 2.96 and 4.24 Å. The activation energy barrier for anion migration is reduced to 0.96 eV. Finally, for an ErOCl structure with 10 at% Cl, 5 at% Er, and 3 at% O vacancies, an interlayer chloride-ion migration path becomes accessible with an energy barrier of 0.82 eV. The Er vacancy leads to undercoordinated chloride ions being able to diffuse along the same Cl layer (Fig. 5E). In this configuration, the mobile chloride ion interacts with three distinctive Er ions at a distance of 3.26, 2.64, and 3.42 Å. Fig. 5F plots the energy barriers for these three distinctive modes of chloride-ion diffusion; the distance between images plotted in the horizontal axis shows that the presence of vacancies results in longer albeit lower energy chloride-ion migration pathways. As such, point defects such as those introduced through aliovalent alloying on the cation sublattice are critical to ensuring high chloride-ion mobility.70 Videos S1A–S1C† provide animated views of the chloride-ion diffusion pathways. Chloride-ion diffusion along the ab-plane provides the primary mode of mass transfer during sintering and yields anisotropically sintered platelets.
Video S1† shows the top view of chloride-ion diffusion between layers A and B in non-defective ErOCl; Video S2† shows the side view of chloride-ion diffusion between layers A and B in ErOCl after 10 at% Cl vacancies; Video S3† shows the top view of diffusion on layer B in ErOCl with 10 at% Cl vacancies, 5 at% Er vacancies, and 3% at. O vacancies. In all three videos, Er atoms are colored green, Cl atoms are black, O atoms are red, and the diffusing Cl atom is distinguished as a dark orange sphere.
Pathways and mechanisms of chloride-ion diffusion have been examined using first-principles nudged elastic band simulations, which indicate that an interlayer low-energy pathway becomes accessible along the ab plane in the presence of point defects. The results show a promising synthetic route to prepare van der Waals materials that can mediate anion diffusion. The DFT simulations demonstrate the promise of enhancing halide-ion mobility through site-selective modification on the cation and anion sublattice. Future work will focus on systematic aliovalent substitution on the cation sublattice and elucidation of the temperature-dependent chloride-ion conductivity of these materials.
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4ce00585f |
| This journal is © The Royal Society of Chemistry 2024 |