Crystal Growth of the Quasi-2D Quarternary Compound AgCrP₂S₆ by Chemical Vapor Transport

Abstract

We report optimized crystal growth conditions for the quarternary compound AgCrP${}_2$S${}_6$ by chemical vapor transport. Compositional and structural characterization of the obtained crystals were carried out by means of energy-dispersive X-ray spectroscopy and powder X-ray diffraction. AgCrP${}_2$S${}_6$ is structurally closely related to the $M_2$P${}_2$S${}_6$ family, which contains several compounds that are under investigation as 2D magnets. As-grown crystals exhibit a plate-like, layered morphology as well as a hexagonal habitus. AgCrP${}_2$S${}_6$ crystallizes in monoclinic symmetry in the space group $P2/a$ (No. 13). The successful growth of large high-quality single crystals paves the way for further investigations of low dimensional magnetism and its anisotropies in the future and may further allow for the manufacturing of few-layer (or even monolayer) samples by exfoliation.

1. Introduction

Among the magnetic quasi-two-dimensional materials that have recently moved in the focus of (quasi-)two-dimensional (2D) materials research [1,2,3], the $M_2$P${}_2$S${}_6$ class of layered materials offers a plenitude of isostructural compounds with different magnetic properties depending on M [4,5]. Thus, $M_2$P${}_2$S${}_6$ allows to investigate fundamental aspects of low dimensional magnetism and several members may be promising for future applications, e.g., complementing non-magnetic (quasi-)2D materials in heterostructures or in spintronic devices [6,7]. Furthermore, future applications in the field of catalysis are conceivable due to the structural similarity to the non-magnetic 2D materials such as graphene or the transition metal dichalcogenide compounds for which such applications are already discussed [8,9,10].

Regarding the crystal structure, the $M_2$P${}_2$S${}_6$ family consists of van der Waals layered compounds which share a honeycomb network of $M^{2+}$ and, most prominently, a dominantly covalent ${\left[{\mathrm{P}}_2{\mathrm{S}}_6\right]}^{4-}$ anion located in the voids of the honeycomb [4,5]. In the bulk, such layers are stacked on top of each other only interacting via weak van der Waals forces. Consequently, these compounds can be easily exfoliated potentially down to a single layer [11,12].

Several isovalent substitution series of $M^{2+}$ by another $M^{{}^{\prime }2+}$ (e.g., (Mn${}_{1-x}$Fe${}_x$)${}_2$P${}_2$S${}_6$ [13], (Mn${}_{1-x}$Ni${}_x$)${}_2$P${}_2$S${}_6$ [14], (Fe${}_{1-x}$Ni${}_x$)${}_2$P${}_2$S${}_6$ [15,16] and (Zn${}_{1-x}$Ni${}_x$)${}_2$P${}_2$S${}_6$ [4]) are reported to exhibit solid solution behavior and, thus, imply a random distribution of the substituents on the honeycomb network, as illustrated in Figure 1a. Beyond isovalent substitution, Colombet et al. [17,18,19,20] demonstrated that a substitution of $M_2^{2+}$ by $M^{1+}{M}^{{}^{\prime }3+}$ also yields several stable compounds. In contrast to the isovalent substitution series however, $M^{1+}$ and $M^{{}^{\prime }3+}$ do not randomly occupy the M positions in the lattice but order either in an alternating or in a zig-zag stripe-like arrangement on the honeycomb, as illustrated in Figure 1b,c, respectively. The former arrangement is attributed to a minimization of repulsive Coloumb interactions (i.e., charge ordering). The latter is observed for compounds for which $M^{1+}$ and $M^{{}^{\prime }3+}$ have notably different sizes and, thus, is dominantly driven by a minimization of lattice distortion and steric effects [4,17].

With $M^{{}^{\prime }3+}$ being a magnetic ion (e.g., V${}^{3+}$ or Cr${}^{3+}$) and $M^{1+}$ being non-magnetic (e.g., Cu${}^{1+}$ or Ag${}^{1+}$), the magnetic sublattices formed in $M^{1+}{M}^{{}^{\prime }3+}$P${}_2$S${}_6$ extend the magnetic structures of the usually magnetically hexagonal $M_2$P${}_2$S${}_6$ compounds by an alternating/triangular and a zig-zag stripe-like magnetic arrangement [4,17,19]. The stripe-like magnetic structure is especially notable, as each stripe of magnetic ions is well isolated from the adjacent magnetic stripes by a stripe of non-magnetic ions. Although the corresponding compound still has a (quasi-)2D layered crystal structure, the magnetic structure can be expected to exhibit 1D magnetic characteristics. Indeed, several indications for such low dimensional magnetism are reported for $M^{1+}{M}^{{}^{\prime }3+}$P${}_2$S${}_6$ with $M=\mathrm{Ag}$ and $M^{{}^{\prime }}=\mathrm{Cr}$ [19,21], making it an interesting compound for further studies.

However, until now only details on the synthesis of AgCrP${}_2$S${}_6$ via solid state synthesis are reported (although Mutka et al. [21] mention CVT grown crystals, they do not report any further details or conditions regarding the crystal growth) [19]. Although small crystals in the $\mathsf{\mu }$m scale could be obtained by solid state synthesis, which allowed for a structural solution based on single crystal X-ray diffraction, significantly larger crystals are needed for detailed investigations of the physical properties including anisotropies.

As a crystal growth method of choice, for macroscopic AgCrP${}_2$S${}_6$ single crystals, the chemical vapor transport (CVT) technique is suitable due to the contained volatile elements such as S and P. Phosphorus and sulfur are both volatile and readily evaporate at elevated temperatures. The generation of volatile intermediate transition metal species for the vapor transport using so-called transport agents is well established [22]. CVT is the crystal growth technique of choice for virtually all ternary $M_2$P${}_2$S${}_6$ compounds [4,5]. For example, Taylor et al. [23] and Nitsche et al. [24] report the successful crystal growth of $M_2$P${}_2$S${}_6$ with $M=\mathrm{Mn},\mathrm{Fe},\mathrm{Ni},\mathrm{Cd} \mathrm{and} \mathrm{Sn}$ by CVT using either chlorine or iodine as transport agent and we present the crystal growth of mixed transition metal phosphorus sulfides of the substitution series (Fe${}_{1-x}$Ni${}_x$)${}_2$P${}_2$S${}_6$ [16] and (Mn${}_{1-x}$Ni${}_x$)${}_2$P${}_2$S${}_6$ [14] by the same technique with iodine as agent. To determine a suitable temperature gradient for the CVT growth of the quarternary compound AgCrP${}_2$S${}_6$, several growth experiments with different temperature profiles were conducted. The temperature profile, which is reported hereafter, resulted in the best crystal size and quality as well as regarding impurity contributions and opens up access to macroscopic AgCrP${}_2$S${}_6$ single crystals. In addition to the crystal growth of AgCrP${}_2$S${}_6$, we also present a comprehensive compositional and structural characterization of the obtained crystals.

2. Materials and Methods

The elemental educts for the crystal growth of AgCrP${}_2$S${}_6$, as listed in

Table 1. Elemental educts used for the CVT growth of AgCrP${}_2$S${}_6$.

Chemical	Specification	Purity
Silver	Powder, APS 4–7 micron	99.9%
Chromium	Powder, −100 + 325 mesh	99.99%
Red phosphorus	Lumps	99.999%
Sulfur	Pieces	99.999%
Iodine	Resublimed crystals	99.9985%

, were obtained from Alfa Aesar and kept in an argon filled glove box for storage and handling.

The crystals obtained from the CVT crystal growth experiments were thoroughly characterized by scanning electron microscopy (SEM) regarding their morphology and topography using a secondary electron (SE) detector and regarding chemical homogeneity via the chemical contrast obtained from a back scattered electron (BSE) detector. For this, a ZEISS EVO MA 10 scanning electron microscope was used. The chemical composition of the crystals was investigated by energy dispersive X-ray spectroscopy (EDX), which was measured in the same SEM device with an accelerating voltage of 30 kV for the electron beam and using an energy dispersive X-ray analyzer.

The crystal structure of the obtained crystals was investigated by powder X-ray diffraction (pXRD), which was measured on a STOE STADI laboratory diffractometer in transmission geometry with Cu-K${}_{\alpha 1}$ radiation from a curved Ge(111) single crystal monochromator and detected by a MYTHEN 1K 12.5${}^{\circ }$-linear position sensitive detector manufactured by DECTRIS. The pXRD patterns were initially analyzed by pattern matching using the HighScore Plus program suite [25]. After the crystallographic phase was identified, a structural refinement of the crystal structure model was performed based on our experimental patterns using the Rietveld method in Jana2006 [26].

3. Crystal Growth via Chemical Vapor Transport

All procedures for the preparation were performed in a glove box under argon atmosphere. The elemental educts silver, chromium, red phosphorus and sulfur were weighed out in a molar ratio of Ag:Cr:P:S = 1:1:2:6 and homogenized in an agate mortar. 0.5 g of reaction mixture were loaded in a quartz ampule (6 mm inner diameter, 2 mm wall thickness) together with approx. 50 mg of the transport agent iodine. Immediately prior to use, the ampule was cleaned by washing with distilled water, rinsing with isopropanol and, subsequently, baking out at 800 ${}^{\circ }$C for at least 12 h in an electric tube furnace. This is done to avoid contamination of the reaction volume with (adsorbed) water. The filled ampule was then transferred to a vacuum pump and evacuated to a residual pressure of ${10}^{-8}$ bar. To suppress the unintended sublimation of the transport agent during evacuation, the end of the ampule containing the material was cooled with a small Dewar flask filled with liquid nitrogen. After reaching the desired internal pressure, the valve to the vacuum pump was closed, the cooling was stopped and the ampule was sealed under static pressure at a length of approximately 12 cm.

The ampule was carefully placed in a two-zone tube furnace in such a way that the reaction mixture was only at one side of the ampule which is referred to as the charge region. As illustrated in Figure 2a, the furnace was initially heated homogeneously to 750 ${}^{\circ }$C at 100 ${}^{\circ }$C/h. The charge region was kept at this temperature for 274.5 h while the other side of the ampule, which is the sink region (see Figure 2b), was initially heated up to 800 ${}^{\circ }$C at 100 ${}^{\circ }$C/h, dwelled at this temperature for 24 h and then cooled back to 750 ${}^{\circ }$C at 1 ${}^{\circ }$C/h. An inverse transport gradient is formed, i.e., transport from sink to charge, to clean the sink region of particles which stuck to the walls of the quartz ampule during the previous preparation steps. This ensures improved nucleation conditions in the following step. Then the sink region was cooled to 690 ${}^{\circ }$C at 0.5 ${}^{\circ }$C/h to gradually form the thermal transport gradient resulting in a controlled nucleation. With a final gradient of 750 ${}^{\circ }$C (charge) to 690 ${}^{\circ }$C (sink), the ampule was dwelled for 100 h. After this period of time, the charge region was cooled to the sink temperature in 1 h before both regions were furnace cooled (i.e., the heating elements were turned off) to room temperature.

Shiny plate-like crystals of AgCrP${}_2$S${}_6$ in the size of approximately $2\mathrm{mm}\times 2\mathrm{mm}\times 100\mathsf{\mu }\mathrm{m}$ were obtained. As example, as-grown single crystals are shown in Figure 2c,d. These crystals exhibit a layered morphology and are easily exfoliated, which is typical for bulk crystals of (quasi-)2D materials.

As shown in Figure 3a, some as-grown crystals exhibit small spherical particles (likely solidified droplets) of a secondary phase (bright) and pieces of a second secondary phase (dark). As the secondary phases are found only on the surface, exfoliating the crystals is sufficient to remove the secondary phases and results in crystals with clean surfaces, as illustrated in Figure 3b.

4. Crystal Morphology and Compositional Analysis

The topographical SE image of a AgCrP${}_2$S${}_6$ crystal in Figure 4a exhibits a flat crystal surface and sharp edges. The terrace close to the upper edge of the crystal is a typical feature of layered systems. Furthermore, on the upper edge of the crystal, some steps can be seen, which form 120${}^{\circ }$ angles, indicative of a hexagonal crystal habitus. The SE image with BSE detector in in Figure 4b shows an overall homogeneous contrast over the surface of the crystal demonstrating that it is chemically homogeneous. At some small areas, a change in contrast is observed. In comparison with the SEM image in SE mode, these spots can be clearly attributed to impurity particles on top of the crystal and not to any region of intergrowth with a secondary phase.

By EDX measurements on multiple spots on several crystals, the mean elemental composition of the crystals was obtained as Ag${}_{1.03\left(2\right)}$Cr${}_{1.06\left(2\right)}$P${}_{2.03\left(2\right)}$S${}_{5.88\left(1\right)}$. This composition is in ideal agreement with the expected composition of AgCrP${}_2$S${}_6$ and the small standard deviations (given in parentheses) indicate a homogeneous elemental distribution and composition.

5. Structural Analysis

The pXRD pattern obtained from exfoliated AgCrP${}_2$S${}_6$ crystals, as shown in Figure 5a, was indexed in the space group $P2/a$ (No. 13), in agreement with the literature [19]. No additional reflection were observed demonstrating the intrinsic phase purity of our crystals.

The $C2/m$ space group, which is typically observed for compounds of the $M_2$P${}_2$S${}_6$ family [23,27], including $M^{2+}{M}^{{}^{\prime }2+}$P${}_2$S${}_6$ compounds of isovalent substitution series (e.g., MnFeP${}_2$S${}_6$ [13], MnNiP${}_2$S${}_6$ [14] and FeNiP${}_2$S${}_6$ [16]), can be ruled out. Assuming a monoclinic unit cell, several observed reflections correspond to Laue indices that are systematically absent for C centering, as they violate the reflection condition $hkl$: $h+k=2n$. Examples are the reflections at $2\theta ={22.94}^{\circ }$ corresponding to 120 and at $2\theta ={31.92}^{\circ }$ corresponding to $\overline{2}11$, as shown in Figure 5b.

This implies that Ag and Cr indeed arrange as zig-zag stripes in AgCrP${}_2$S${}_6$ and do not just randomly occupy the corners of the structural honeycomb network, as it is the case for isovalent substitutions. While the former scenario breaks the mirror symmetry of the $C2/m$ space group of the Fe${}_2$P${}_2$S${}_6$ aristotype [28], which results in a $P2/a$ space group, the latter scenario would not. Furthermore, a $C2/c$ space group, as reported, e.g., for CuCrP${}_2$S${}_6$ [17] with a triangular arrangement of the two transition element cations, can be ruled out based on the same considerations.

Starting from the crystal structure model proposed by Colombet et al. [19], a refined crystal structure model is obtained using the Rietveld method which is sufficient to describe our experimental pattern with good agreement, as shown in Figure 5a. The obtained lattice parameter and reliability factors are summarized in

Table 2. Top: Summary of experimental parameters of the pXRD experiment on AgCrP${}_2$S${}_6$, extracted lattice parameters and reliability factors of the structural model obtained by the Rietveld method. Bottom: Refined crystal structure model of AgCrP${}_2$S${}_6$ and isotropic displacement parameters with standard deviations given in parentheses. All sites were treated as fully occupied.

Experiment
Temperature (K)				293(2)
Radiation				Cu-K${}_{{\alpha }_1}$
Wavelength (Å)				1.54059
${\theta }_{\min }$ (${}^{\circ }$)				10.00
${\theta }_{\mathrm{step}}$ (${}^{\circ }$)				0.03
${\theta }_{\max }$ (${}^{\circ }$)				120.13
Crystal Data
Crystal System				monoclinic
Space Group				$P2/a$
a (Å)				5.8832(1)
b (Å)				10.6214(2)
c (Å)				6.7450(3)
$\beta$ (${}^{\circ }$)				106.043(2)
Refinement
Goodness-Of-Fit				2.13
$R_{\mathrm{p}}$ (%)				2.08
$w{R}_{\mathrm{p}}$ (%)				3.11
$R_{\mathrm{F}}$ (%)				5.45
Label	Type	Wyck	x	y	z	$U_{\mathrm{iso}}$
						($\times {10}^{-3}{\mathsf{\AA }}^2$)
Ag1	Ag	$2e$	0.75	0.4364(2)	0	34(1)
Cr1	Cr	$2e$	0.25	0.9229(4)	0	20(2)
P1	P	$4g$	0.2979(6)	0.2466(4)	0.1659(6)	2(1)
S1	S	$4g$	0.9792(7)	0.2309(5)	0.2336(8)	18(2)
S2	S	$4g$	0.9880(5)	0.9233(4)	0.2165(8)	7(1)
S3	S	$4g$	0.4777(7)	0.3947(4)	0.2802(7)	34(1)

(top) and the refined structural model is given in the same table on the bottom and is illustrated in Figure 6. The strongest disagreement between model and experiment is observed for the high intensity 001 reflection at $2\theta ={13.67}^{\circ }$. As this reflection corresponds to the stacking of layers, it is most prominently affected by any kind of disorder or defects influencing the stacking. Due to the weak structural interaction between layers, which are only based on weak van der Waals forces, the $M_2$P${}_2$S${}_6$ compounds are prone to stacking faults and twinning between layers. In the presence of such defects, the shape of the corresponding 001 reflection is altered, which may be a reason for the observed deviation between the experiment and the model without defects.

Additionally, the experimental pattern exhibits significantly altered reflection intensities compared to an initial model, which are attributed to a strongly preferred orientation of the crystallites in the investigated sample. Due to the layered structure with only weak van der Waals interactions between layers, the powder particles obtained from grinding AgCrP${}_2$S${}_6$ crystals are plate-like and tend to lie flat on the sample holder. Thus, reflections with a dominant l component (e.g., 001) exhibit higher intensities than expected for spherical crystallites in transmission geometry. To adjust for this effect in the model, the method proposed by March [29] and extended by Dollase [30] was used. However, the preferred orientation in AgCrP${}_2$S${}_6$ is strongly pronounced, such that it might be beyond the limit of what the semi-empirical March-Dollase model is capable of describing accurately. This may furthermore contribute to the deviation between model and experiment around the 001 reflection.

The refined crystal structure model for AgCrP${}_2$S${}_6$ shows that the Ag–S bonds are notably longer than the Cr–S bonds, as expected based on the difference between the size of the transition element cations (e.g., ionic radii for octahedral coordination: $r\left({\mathrm{Ag}}^{1+}\right)=1.15$ Å and $r\left({\mathrm{Cr}}^{3+}\right)=0.62$ Å [31]). These different bond lengths result in a distortion of the structure compared to the aristotype Fe${}_2$P${}_2$S${}_6$, which can be clearly observed, e.g., in Figure 6a. In detail, the CrS${}_6$ coordination environment remains antiprismatic (i.e., close to octahedral with a slight trigonal elongation along $c^{*}$) with the faces above and below the shared plane of the transition elements being parallel to each other. However, the AgS${}_6$ coordination environment as well as for the P${}_2$S${}_6$ environment are distorted in such a way that the faces above and below the transition element plane are not parallel to each other. In the view along the $c^{*}$ direction in Figure 6c, this distortion manifests in Ag and P${}_2$ being shifted off-center in their respective sulfur coordination environments away from the closest Cr positions. Meanwhile, Cr is located exactly in the center of the CrS${}_6$ unit. The observation that the CrS${}_6$ unit is closer to an ideal octahedral coordination environment than the AgS${}_6$ unit can be understood considering the local charge density (i.e., ionic size and charge). Cr${}^{3+}$ is small and highly charged and, thus, interacts with the surrounding S atoms stronger than the comparable large and less charged Ag${}^{1+}$. Another notable structural aspect is the strong distortion of the [P${}_2$S${}_6$]${}^{4-}$ units that demonstrates how flexible this covalent complex anion is. This complex anion is a common and characteristic building unit in the $M_2$P${}_2$S${}_6$ family and its flexibility may indicate that several more compounds of the general formula $M^{1+}{M}^{{}^{\prime }3+}$P${}_2$S${}_6$ are stable but have not been synthesized yet.

6. Summary and Conclusions

We report optimized crystal growth conditions for the quarternary compound AgCrP${}_2$S${}_6$ via Chemical Vapor Transport (CVT). A temperature profile adapted from the CVT growth of ternary $M_2$P${}_2$S${}_6$ compounds is sufficient to yield crystals of the target AgCrP${}_2$S${}_6$ phase in the mm-size. On some crystals, traces of a superficial impurity phase is found which could be readily removed by exfoliation.

The as-grown crystals exhibit a plate-like, layered morphology as well as a hexagonal habitus and have the expected composition of AgCrP${}_2$S${}_6$ based on EDX spectroscopy. The pXRD pattern is indexed in the space group $P2/a$ in agreement with the literature [19]. The $P2/a$ space group, on which the zig-zag type arrangement of M and $M^{{}^{\prime }}$ is based on, can be well distinguished from, e.g., the $C2/m$ and $C2/c$ space groups due to reflections that are systematically absent for C centering. Starting from the model of Colombet et al. [19], a refined structural model is obtained using the Rietveld method. This model contains a notable distortion of the AgS${}_6$ and P${}_2$S${}_6$ coordination environments, while the CrS${}_6$ units remain antiprismatic with a slight trigonal distortion.

The zig-zag stripe-like arrangement in AgCrP${}_2$S${}_6$ and the alternating arrangement of M and $M^{{}^{\prime }}$, which is reported, e.g., for CuCrP${}_2$S${}_6$, are promising to yield interesting magnetic and electronic structures. While only few such quarternary phosphorus sulfide compounds have been synthesized until now, many more combinations of a 1+-ion and a 3+-ion can be expected to form analogous compounds. Furthermore, the fundamental idea of replacing $M^{X+}$ by $M_{0.5}^{(X-1)+}{M}_{0.5}^{{}^{\prime }(X+1)+}$ may be adoptable to the closely related structures such as, $M_2^{3+}$(Si,Ge)${}_2$Te${}_6$ compounds.

The single crystals of AgCrP${}_2$S${}_6$ that were obtained using the presented growth conditions allow for studies of the low dimensional magnetic interactions including the magnetic anisotropy of this compound in the future, which may lead to a better fundamental understanding of low dimensional magnetism. Furthermore, the van der Waals layered structure makes exfoliation easily possible and, thus, our successful growth of single crystals paves the way for further manufacturing of few-layer or even monolayer samples of AgCrP${}_2$S${}_6$.