1 Introduction
High-temperature metal solutions as flux are indispensable methods to grow single crystals with diverse exotic physical properties [
1,
2]. Among a variety of metal flux, Al flux has some unique features. For example, it can dissolve a large number of elements, especially boron. Thus, Al flux is an important method to grow metal borides single crystals, such as transition metal and rare earth borides [
3,
4]. Furthermore, Al can be dissolved readily not only in non-oxidizing acids but also in bases. This property is very important to synthesis single crystals that are easy to decompose in acids but not in bases. In addition, Al flux can also be used as a reaction flux to grow complex (ternary or quaternary) aluminides, especially alumosilicides and alumogermanides, because the phase diagrams of the Al-Si and Al-Ge show low eutectic points without forming any binary phases. A large amount of intermetallic alumosilicides and alumogermanides have been prepared from liquid Al, such as RENiAl
4Ge
2 (RE = Sm, Tb, Y) [
5], Sm
2Ni(Ni, Si)Al
4Si
6 [
6], and RE
8Ru
12Al
49Si
9 (Al
xSi
12-x) (RE = Sm, Pr) etc [
7].
In this work, we discovered two new quaternary compounds Gd
2CoAl
4T
2 (T = Si, Ge) using Al flux. They are isostructural to Tb
2NiAl
4Ge
2 and Er
2CoAl
4Ge
2 [
8,
9]. Detailed study on physical properties shows that these materials are metals with long-range antiferromagnetic (AFM) transition at low temperature. These AFM orderings can be mainly ascribed to the magnetic interaction of local moments of Gd
3+ ions.
2 Experiment
Gd2CoAl4T2 (T = Si, Ge) single crystals were grown by the Al flux method. Gd chunks (99.99%), Co lumps (99.99%), Al grains (99.99%) and Si or Ge granular (99.9999%) were mixed together in a molar ratio of and placed in an alumina crucible. Then, the crucible was sealed in the quartz ampoule under partial argon atmosphere. The sealed quartz ampoule was heated up to 1323 K for 12 h and then kept at that temperature for 4 h. Finally the ampoule was cooled down slowly to 1023 K with 5 K/h. At this temperature, the ampoule was decanted with a centrifuge to extract single crystals from the Al flux. X-ray diffraction (XRD) patterns were collected using a Bruker D8 X-ray diffractometer with Cu radiation ( nm) at room temperature. The elemental analysis was performed using energy-dispersive X-ray spectroscopy (EDX) analyses. Electrical transport measurements were carried out in a Quantum Design physical property measurement system (PPMS-14). The electrical resistivity was measured using a four-probe method on rectangular-shaped single crystals. The current flows in the ab plane of the crystal. Magnetization measurements were performed in a Quantum Design magnetic property measurement system (MPMS3).
3 Results and discussion
Figure 1(a) shows the powder XRD pattern of ground Gd
2CoAl
4Ge
2 single crystals. It can be fitted well by using the structure of Er
2CoAl
4Ge
2 with the space group I4/mmm (No. 139). The fitted
a- and
c-axial lattice parameters are 4.1392(3) Å and 19.543(2) Å for Gd
2CoAl
4Ge
2 (also shown in Table 1), which are larger than those in Er
2CoAl
4Ge
2 (
Å and
Å) [
9]. It can be ascribed to the larger ionic radius of Gd than Er. On the other hand, when exploring the Gd-Co-Al-Si system we found that the powder XRD pattern of ground single crystals can also be fitted using above structure. The obtained lattice parameters from fitting give
Å and
Å for Gd
2CoAl
4Si
2, which are smaller than those in Gd
2CoAl
4Ge
2, because of smaller ionic radius of Si when compared to Ge. Correspondingly, the volume of unit cell
V for Gd
2CoAl
4Si
2 is also smaller than that for Gd
2CoAl
4Ge
2 (Table 1). The composition analysis further confirms the formation of Gd
2CoAl
4Si
2. The atomic ratio of
determined from the EDX measurement is
when setting Gd as 2, close to the chemical formula of Gd
2CoAl
4Si
2. Figures 1(c) and (d) show the XRD patterns of Gd
2CoAl
4T
2 single crystals. There are only the (00
l) diffraction peaks appearing, indicating that the crystal surface is parallel to the
ab-plane and the diffraction peaks of Gd
2CoAl
4Ge
2 shift to lower angles compared to the Gd
2CoAl
4Si
2 single crystal, consistent with the larger
c axial lattice parameter of the former. The thick plate-like shapes of crystals [insets of Figs. 1(c) and (d)] are consistent with the crystal orientation and the layered structure of Gd
2CoAl
4T
2 mentioned below.
As shown in Fig. 2(a), the structure of Gd
2CoAl
4T
2 is built up by stacking CoAl
4T
2 slabs and bilayer of Gd atoms alternatively along the
c axis. Similar to the NiAl
4Ge
2 unit in Tb
2NiAl
4Ge
2 [
8], the CoAl
4T
2 slab has an antiflourite (anti-CaF
2) structure type, where T and Co occupy the site of Ca and F is replaced by Al. Thus, each Al atom is coordinated by two T and two Ni atoms, forming a tetrahedron. Two layers of these AlCo
2T
2 tetrahedra connect each other via edge sharing and form a two-dimensional
slab [Fig. 2(b)]. On the other hand, the local structure of CoAl
4T
2 slab can also be described as following [
8]. Each Co atom is surrounded by eight Al atoms, consisting of a slightly distorted CoAl
8 cube. These distorted CoAl
8 cubes connect each other by edge sharing and form an infinite layer with one half of the cubes filled with a Co center. The T atoms then are located above and below the Al
8 cubes which are not centered by Co atoms [Fig. 2(c)]. The fitted atomic positions for Gd, Co, Al and T are shown in Table 1. The values are close to those in Tb
2NiAl
4Ge
2 and Er
2CoAl
4Ge
2 [
9,
8]. As shown in Table 1, the
,
and
(in the
ab plane) for Gd
2CoAl
4Si
2 are slightly smaller than those for Gd
2CoAl
4Ge
2which can be ascribed to the smaller lattice parameters of the former. But the
and
(out of
ab plane) for the former are abnormally larger than those for the latter. It suggests that the distorted CoAl
8 cubes in Gd
2CoAl
4Si
2 expand along the
c axis and one of possible explanation is that the stronger bonding strength of Si-Al weakens those of Co-Al and Al-Al because of larger electronegativity of Si than Ge.
Figure 3(a) shows the temperature dependence of magnetic susceptibility
of the Gd
2CoAl
4Ge
2 single crystal at
kOe along the
ab plane and the
c axis with zero-field-cooling (ZFC) and field-cooling (FC) modes. There are sharp drops below 20 K in both ZFC and FC curves, and it can be ascribed to the onset of long-range AFM ordering. The Weiss temperature
TN determined from the peak of d
/d
T(
T) curves with ZFC mode is 20.1 and 20.0 K for
and
[inset of Fig. 3(a)] [
10]. These values are close to that in Tb
2NiAl
4Ge
2 [
8]. The decrease of
below
TN is more obvious than
, strongly suggesting that the spin is ordered antiferromagnetically in the
ab-plane (easy plane of magnetization). When
, the
decrease with increasing temperature for both field directions and they can be fitted very well using the Curie–Weiss (CW) law
, where
C is Curie constant and
θ is the CW temperature. The linear fits for the curves of
vs.
T between 50 and 300 K with ZFC mode (solid lines in Fig. 3(c)) give
emu·g
-1·Oe
-1·K and
K for
, and
emu·g
-1·Oe
-1·K and
K for
. The negative values of
θ indicate the AFM interaction in this material, consistent with the AFM ordering below
TN. The values of
C correspond to the effective moment
and 8.567(1)
μB/Gd for
and
. These values are close to the theoretical value of free Gd
3+ ion (7.94
μB/Gd). It suggests that Co ions may not have local moments in Gd
2CoAl
4Ge
2 as in the case in Tb
2NiAl
4Ge
2 and other rare earth-cobalt-aluminum compounds [
8,
11]. For Gd
2CoAl
4Si
2, there is also a drop in the
curve at low temperature [Fig. 3(b)], corresponding to the AFM transition. The determined
TN from the peak of d(
χT)/d
T(
T) curves is 26.4 and 26.8 K for
and
[inset of Fig. 3(b)], slightly higher than those in Gd
2CoAl
4Ge
2. It could be due to the smaller
in the former. The easy plane of magnetization in Gd
2CoAl
4Si
2 is also the
ab-plane because of the faster drop of
than
. From the linear fits of
curves [solid lines in Fig. 3(c)], the obtained
C and
θ are 0.031690(7) emu·g
-1·Oe
-1·K and
K for
, when
emu·g
-1·Oe
-1·K and
K for
. Again, the negative values of
θ indicate that the dominant magnetic interaction at high temperature is AFM, which evolves into the long-range AFM ordering below
TN. On the other hand, the corresponding
is 8.255(9) and 8.069(1)
μB/Gd for
and
, respectively, comparable with those values in Gd
2CoAl
4Ge
2 and theoretical value of Gd
3+ ions. In addition, the nearly perfect overlapping between the ZFC and FC
curves indicate the absence of magnetic glassy state in these compounds. Figure 3(d) shows the
loops at 2 K for both field directions. For Gd
2CoAl
4Ge
2, the
loops exhibit the linear field dependence without the trend of saturation and hysteresis, irrespective of the field directions. It is consistent with the AFM orders in these materials. For Gd
2CoAl
4Si
2, the
loop also shows similar linear field dependence when
. Interestingly, the
curve exhibits a sharp jump at
H ~ 29 kOe when
, implying that there is a spin-flop transition, i.e., the spin directions are almost perpendicular to the field direction and nearly antiparallel to each other [
12-
14].
Figure 4 shows the temperature dependence of the in-plane resistivity
of Gd
2CoAl
4T
2 single crystals. These compounds exhibit metallic behavior in the whole temperatures range. The most notable feature of
is a sharp drop at low temperature. The determined temperature
from the peak of d
ρab/d
T curve [Inset (a) of Fig. 4] is 20.0 and 25.3 K for Gd
2CoAl
4Ge
2 and Gd
2CoAl
4Si
2, respectively, consistent with the transition temperature observed in the magnetization and heat capacity data (shown as below). The drop in the
curve originates from the rapid decrease in the spin disorder scattering below AFM transition. when taking into account the scattering of conduction electrons by phonons and antiferromagnetic magnons, the temperature dependence of
for
can be expressed by [
15,
16]
where is residual resistivity, A and B is the constants of the electron-phonon and electron-magnon scattering term at low temperature, respectively, kB is Boltzmann constant, Δ is the magnon gap energy. Both ρab curves can be fitted very well [Inset (b) of Fig. 4] and the obtained , A, B, and Δ is 35.3(2) μΩ·cm, 0.0107(2) nΩ·cm·K-5, 0.148(2) μΩ·cm·K-2, and 20.2(3) K for Gd2CoAl4Ge2, and 22.9(1) μΩ·cm, 0.00211(5) nΩ·cm·K-5, 0.0720(7) μΩ·cm·K-2, and 25.3(3) K for and Gd2CoAl4Si2, respectively. The fitted values of Δ are close to their AFM ordering temperature.
Figure 5 shows the temperature dependence of heat capacity for Gd2CoAl4Ge2 single crystal. There is a λ-type anomaly with the peak position located at 19.95 K. It is the evidence of bulk long-range AFM ordering and the transition temperature is also consistent with that temperature obtained from d/dT curve [inset of Fig. 3(a)]. Moreover, at kOe the curve is perfectly overlapped with that at zero field, indicating that that the field of 50 kOe has a negligible influence on the AFM transition of Gd2CoAl4Ge2. Assuming that the total heat capacity consists of electronic (Ce), lattice () and magnetic (Cm) components, Cm can be estimated by the subtraction of Ce and . Because the absence of non-magnetic counterpart of this material, we fitted the heat capacity data in the high-temperature range ( K) well away from TN when only considering Ce and and then extrapolated to 2 K. Using the Debye function, it has
where γ is the Sommerfeld coefficient, R is the gas constant ( J·mol-1·K-1), is Debye temperature, C is a fitting parameter. The fitted γ and is 355(2) J·mol-1·Gd K-2 and 321.9(6) K, respectively. The fitted curve is shown in Fig. 5(a) (red solid line) and after subtracting the electronic and lattice contributions from , we obtain the curve [Fig. 5(b)]. Then, the magnetic entropy can be calculated using the integral . As shown in Fig. 5(b), the increases gradually when increasing temperature and then tends to saturate above TN. The obtained Sm up to 40 K is about 10.8 J·mol-1·Gd K-1, which is about 62% of theoretical value [Rln(2J + 1) = 17.3 J·mol-1·Gd K-1 for Gd3+ ions with ]. This discrepancy may originate from an overestimation of the electronic and lattice contributions which may cause a reduction in . Similar to Gd2CoAl4Ge2, there is also a peak at 26.97 K for Gd2CoAl4Si2 at zero field [Fig. 4(c)], consistent with the results of magnetization measurement. In contrast to Gd2CoAl4Ge2, at high field (H=50 kOe), the peak of shifts to lower temperature slightly ( K), implying the suppression of AFM order by field. This could be related to the metamagnetic transition in Gd2CoAl4Si2. After subtracting the electronic and lattice contributions from , the obtained is shown in Fig. 4(d) and it exhibits similar behavior as in Gd2CoAl4Ge2. The estimated Sm up to 40 K is 8.9 J·mol-1·Gd K-1, which is about 51% of theoretical value for Gd3+ ions. Again, this discrepancy may also be due to the overestimation of the Ce and .
4 Conclusion
In summary, we report the physical properties of two new quaternary compounds Gd2CoAl4T2 (T = Si, Ge) single crystals grown by using Al flux. Gd2CoAl4T2 belong to the tetragonal space group I4/mmm and the fitted a- and c-axial lattice parameters are 4.1392(3) Å and 19.543(2) Å for Gd2CoAl4Ge2 and 4.1099(1) Å and 19.4146(8) Å for Gd2CoAl4Si2. The crystal structure consists of layers of CoAl4T2 slab separated by bilayers Gd atoms. Both compounds are AFM metals and the TN is about 20 and 27 K for Gd2CoAl4Ge2 and Gd2CoAl4Si2, respectively. Heat capacity measurements confirm the bulk feature of these AFM transitions. The long-range antiferromagnetism in these materials may solely originate from the magnetic interaction between Gd3+ ions and the Co ions should have negligible local moments. Moreover, the Gd2CoAl4Ge2 exhibits metamagnetic transition when the field is along the ab plane and larger than 29 kOe.
Higher Education Press and Springer-Verlag GmbH Germany, part of Springer Nature 2019