Table Of ContentMagnetization and transport properties of single crystalline RPd P (R=Y, La–Nd,
2 2
Sm–Ho, Yb)
Gil Drachuck, Anna E. B¨ohmer, Sergey L. Bud’ko, and Paul C. Canfield
Department of Physics and Astronomy, Iowa State University, Ames, IA 50011,
USA and Ames Laboratory, Iowa State University, Ames, IA 50011, USA
(Dated: January 13, 2016)
SinglecrystalsofRPd P (R=Y,La–Nd,Sm–Ho,Yb)weregrownoutofahightemperaturesolu-
2 2
tionrichinPdandPandcharacterizedbyroom-temperaturepowderX-raydiffraction,anisotropic
temperature- and field-dependent magnetization and temperature-dependent in-plane resistivity
measurements. In this series, YPd P and LaPd P and YbPd P (with Yb+2) are non local-
2 2 2 2 2 2
moment bearing, whereas CePd P and PrPd P order at low temperature with a ferromagnetic
2 2 2 2
component along the crystallographic c-axis. The rest of the series manifest low temperature anti-
ferromagnetic ordering. EuPd P has Eu+2 ions and both EuPd P and GdPd P have isotropic
2 2 2 2 2 2
paramagnetic susceptibilities consistent with L=0 and J =S = 7 and exhibit multiple magnetic
2
transitions. For R=Eu-Dy, there are multiple, T > 1.8 K transitions in zero applied magnetic
field and for R=Nd, Eu, Gd, Tb, and Dy there are clear metamagnetic transitions at T=2.0 K
6
for H < 55 kOe. Strong anisotropies arising mostly from crystal electric field (CEF) effects were
1
observed for most magnetic rare earths with L (cid:54)= 0. The experimentally estimated CEF param-
0 eters B2 were calculated from the anisotropic paramagnetic θ and θ values and compared to
2 0 ab c
theoretical trends across the rare earth series. The ordering temperatures as well as the polycrys-
n talline averaged paramagnetic Curie–Weiss temperature, θ , were extracted from magnetization
ave
a and resistivity measurements, and compared to the de-Gennes factor.
J
2
1 I. INTRODUCTION surementsweredoneonpolycrystallinesamples,thusthe
anisotropicpropertieswereaveraged-outoverallcrystal-
] lographic directions.
l The RT X (R = Y, La-Lu; T = transition metal;
e 2 2 In the present work, a systematic study of the
- X = Si, Ge, P, As) family of intermetallic com-
r anisotropic magnetic properties and electrical resistiv-
pounds had been extensively studied over the past 50
t
s years [1]. Nearly all RT X compounds crystallize into ity of RPd2P2 single crystals is presented for R = Y,
. 2 2 La-Nd, Sm-Ho, Yb. The experimental techniques used
t the ThCr Si (space group I4/mmm), where the rare
a 2 2 in the crystal growth and characterization are described
m earth(R)ionoccupiesthe2(a)sitewhichhasatetragonal
in Section II. The experimental results are summarized
point symmetry [2]. Moreover, the transition metal ions
- and presented in Section III, starting with x-ray diffrac-
d in this structure, except for Mn [1] (and perhaps Fe in
tion followed by physical properties of R = Y and La
n LuFe Ge [3–5]),bearnomagneticmoments,meaningall
2 2 members, combined, and then separately for all other
o the magnetic properties are a consequence of the R local
c members. Discussions of trends along the series, such
moment. Therareearthionsinteractviathelongrange,
[ as ordering temperature, anisotropic Curie-Weiss (CW)
indirect, Ruderman-Kittel-Kasuya-Yosida (RKKY) type
temperatures and crystal electric field (CEF) effects will
1 interactions, mediated by the conduction electrons [1].
v Therefore,aninterplaybetweenFermisurfacenesting,or bepresentedinSectionIV,followedbyabriefconclusion
5 in Section V.
maximainthegeneralizedmagneticsusceptibility(χ)[6],
2
and local moment anisotropy is expected to lead to a
0
multitude of potential magnetic transitions and ground
3
0 states. Namely, incommensurate or commensurate mag- II. EXPERIMENTAL
. netic propagation vectors, multiple transitions from one
1
0 to the other and metamagnetism are expected. Growth of intermetallic compounds containing signifi-
6 In the specific case of the RPd P , limited work has cantamountsofvolatileelementsisoftenchallengingdue
2 2
1 been done, mainly due to the high cost of palladium and to the apparent conflict between accessible liquidus sur-
v: the difficulties associated with the volatility of phospho- facesandallowablevaporpressure. Overthepastseveral
Xi rusathightemperatures. TheRPd2P2serieswassynthe- years we have been developing mixed metal-chalcogen
sizedinthe1983[7],andexceptoneearlyphoto-emission and metal-pnictogen fluxes that alleviate this problem
r study on EuPd P [8], which revealed that Eu is in a by greatly reducing the partial vapor pressure over the
a 2 2
divalent state, the system has been mainly overlooked. melt. Solutiongrowthusingsulfur[13],nitrogen[14]and
Recently, work has been done on the CePd P com- phosphorous [15] have been possible by careful identifi-
2 2
pound, revealing its ferromagnetic (FM) Kondo-lattice cation and testing of binary melts for use.
nature[9–12]. ThemagneticpropertiesofGdPd P have Given that RPd P crystals have equal amounts of
2 2 2 2
been reported as well [12], but only to serve as reference palladium and phosphorous, the Pd P binary eutec-
67 33
to CePd P . Moreover, all the above mentioned mea- tic (with T (cid:39) 780 ◦C) was identified as a promising
2 2 eu
2
( a ) E u P d P
E u P d P 2 2
4 0 0 0 2 2 )
2
1
(1
)
s
t
3 0 0 0
n
u
o
C
(
y2 0 0 0
it ) )
Intenis1 0 0 0 (002) (101) (110)(103 (004 (200) (202)(114)(211)(105) (213) (204) (220)(116)(222)(215)(301)(312)(224) (008) ( b )
0
2 0 3 0 4 0 5 0 6 0 7 0 8 0
q (
2 d e g )
FIG. 1: (a) Powder X-ray diffraction pattern of EuPd P with (hkl) values for all peak positions. The arrow indicates a peak
2 2
from the Pd-P binary phases (b) Single crystal of EuPd P displayed on a 1 mm grid.
2 2
melt. We first tested the Pd P binary by placing sto- the residual resistivity ratio (ρ(300K)/ρ(2.0K)) mono-
67 33
ichiometric amounts of elemental P and Pd (powder) in tonicallydecreasesacrosstheR+3 membersoftheseries,
one side of a 2ml fritted crucible set [16]. The crucibles again suggesting that as heavier rare earth are used the
weresealedinamorphoussilicatube[17]under0.2atmo- single crystal growth becomes less stable.
spherespartialpressureofArgonandthenheatedover24
DC magnetization measurements were performed in a
hoursto1100◦Cfollowedbydecantingtheliquid. Given
Quantum Design Magnetic Property Measurement Sys-
that there was (i) no apparent phosphorus migration or
tem(MPMS),superconductingquantuminterferencede-
significant vapor pressure at high temperature and (ii)
vice (SQUID) magnetometer (T = 1.8 - 300 K, H =
max
noapparentcrucibleorampuleattack,thePd P melt
67 33 55kOe). Allsamplesweremanuallyalignedwithin5de-
was used for these growths.
greesofaccuracy,tomeasurethemagnetizationalongthe
desiredaxes. Thesamplesweretightlysqueezedbetween
Single crystals of the RPd P (R=Y, La–Nd, Sm–Ho,
2 2 twoplasticstrawsforH (cid:107)ab orientation. Inthisconfigu-
Yb)seriesweregrownoutofselfflux, byadding5−10%
rationthereisnoriskofsamplerotationduetotorquefor
rare earth into the Pd P eutectic [17, 18]. The ini-
67 33 samples with magnetic anisotropy, and no addendum to
tial elements where placed into the bottom 2 ml alu-
the magnetic signal. For the H (cid:107)c orientation, the sam-
mina crucible of a fritted crucible set [16], and sealed
ples where mounted between two strips of Teflon tape
in amorphous silica ampules under a partial argon at-
suspended over the edges of two internal straws inserted
mosphere. The ampules were heated to 300 ◦C in 3
into an external straw. Given that the signal of the mo-
hoursanddwelledtherefor6hours,inordertoallowthe
ment bearing samples was much larger than that from
phosphorous and palladium to react, therefore reducing
the addendum for this configuration, only the data from
the risk for explosions. Subsequently, the ampules were
LaPd P ,YPd P andYbPd P measurementwerecor-
heated over 10-12 hours to 1180 ◦C where they dwelled 2 2 2 2 2 2
rected for addendum contribution.
for 3 additional hours, then cooled, over 90-120 hours to
930 ◦C. At that point, the excess molten flux was de- The temperature-dependent magnetization (M(T)) of
themoment-bearingmembers,wasmeasuredwithanex-
cantedand, givenPd-content, recycled. Thegrowncrys-
ternal magnetic field of H =1 kOe. Due to the tetrago-
talshadplate-likemorphologywiththec-axisperpendic-
nalsymmetryoftheRPd P ,thepolycrystallineaverage
ular to the plate surface. An optical image of EuPd P 2 2
2 2 method, taking χ = 1χ + 2χ , could be applied to
is given as an example in Fig. 1(b). As R progressed ave 3 c 3 ab
eliminate CEF effects [19]. The transition temperatures
pastHo,theRPd P compoundsbecamehardertogrow
2 2
forallantiferromagneticallyorderedcompoundswerein-
in single crystal form. Attempts to grow ErPd P were
2 2
ferred from d(χ T)/dT [20].
unsuccessful,therefore,TmPd P andLuPd P growths ave
2 2 2 2
were not attempted. The fact that YbPd P could be Additional measurements of DC magnetization for
2 2
grown is very likely associated with the fact that Yb is TbPd P and DyPd P up to 140 kOe were performed
2 2 2 2
divalent and its unit cell volume is between GdPd P using an extraction magnetometer of the ACMS option
2 2
and TbPd P . As will be shown in the next section, of a Quantum Design Physical Property Measurement
2 2
3
and/or dρ/dT. The error bars due to mass uncertainty
1 0 . 0
) ( a ) and different ranges of CW fit are about 2% for effective
(Å 9 . 9 Y L a C e N d G d D y c - a x is mθ o.mTenhteaunndce1rt0a%intfyorinpatrhaemsaagtnuertaitcedCWmotmemenptevraatluueresis,
axis 9 . 8 P r S m T b H o esptimated to be about 2% as well.
- Y b
c E u
9 . 7
)
(Å 4 . 2 III. RESULTS
is
x 4 . 1
a a - a x is A. Powder X-ray Diffraction
-
a
4 . 0
1 7 5 ( b ) The unit cell parameters for the RPd2P2 compounds
L a weredeterminedatroom-temperature,usinggroundsin-
) glecrystals,withaRigakuMiniflexpowderX-raydiffrac-
(Å1 7 0 C e E u tometer(CuKαradiation). TheX-raydiffraction(XRD)
e P r pattern ofEuPd2P2 is shown inFig. 1(a)as an example.
m N d Allmajorpeakscanbeidentifiedandareconsistentwith
lu c e ll v o lu m e the reported ThCr Si (I4/mmm,139) tetragonal struc-
o S m 2 2
V1 6 5 ture. Insomecases,smallpeaksassociatedwiththePd-P
ll G d binaryphaseswerealsodetected. Thelatticeparameters,
e Y b
C T b a-axis and c-axis (Fig. 2(a)), were refined for all series
Y D y
1 6 0 H o members and are summarized in Table I. The unit cell
volume is shown in Fig. 2(b). The values are in excel-
lent agreement with previously published data [7]. The
3 6 4 0 6 0 6 5 7 0
trivalent rare earth members of the RPd P series show
2 2
A t o m ic N u m b e r
a standard lanthanide contraction in volume. It is worth
noting that the effects of divalency (for R=Eu and Yb)
are very anisotropic with an ∼ 0.15˚A increase in the a-
FIG. 2: (a) Powder x-ray diffraction a-axis and c-axis unit lattice parameter and a comparable ∼0.15˚A decrease in
cellparameter. (b)Unitcellvolumecalculatedfromtheabove
thec-latticeparameter. Thisisanunusualeffect,sincein
values.
relatedcompounds,forwhichEuisdivalent,suchbehav-
iorwasnotobserved. InEuRu P [23]forexample, only
2 2
c-lattice parameter is increased, in EuCu Ge [24] only
System (PPMS). For these measurements the samples 2 2
a-lattice parameter is increased and in EuCo Ge [25]
were glued to a Kel-F disk to ensure H (cid:107) c direction 2 2
and EuNi Ge [24] both a- and c-lattice parameter are
of the applied field. The signals from the samples were 2 2
increasedrelativetothetrivalentlanthanidecontraction.
significantly larger than the diamagnetic signal from the
disk [21], so no correction for the disk’s signal was used.
Resistivity measurements were performed within the RPd2P2 a (˚A) c (˚A) Volume (˚A3)
temperature-field environment of the MPMS system us- Y 4.05 9.84 161.4
La 4.12 9.89 174.1
ing a Linear Research Inc. LR-700 4-wire AC resistance
Ce 4.16 9.90 171.1
bridge. The samples were shaped into bars with typical
Pr 4.14 9.88 169.5
dimensions of 1.5 × 0.8 × 0.3 mm3 mm. Epotek-H20E
Nd 4.12 9.88 168.1
silverepoxywasusedtocontactPtwires(0.05mmdiam-
Sm 4.10 9.87 165.8
eter)tothesamples. Typicalcontactresistanceswere1–2
Eu 4.16 9.74 170.8
Ω. Theplate-likemorphologyofthecrystalshasonlyal-
Gd 4.08 9.86 164.2
lowedmeasurementswithcurrentflowingintheab-plane.
Tb 4.06 9.86 162.5
Theresistivetransitiontemperaturevalueswereinferred
Dy 4.02 9.85 161.2
from anomalies in dρ/dT [22]. The transition tempera- Ho 4.03 9.84 159.8
tures from magnetization and resistivity measurements Yb 4.10 9.72 163.0
are summarized below in Table II.
The uncertainty in absolute value of resistivity due to
TABLE I: Lattice parameters and unit cell volume of the
the measurement of the sample’s dimension and sample
RPd P series. The uncertainty is ∼0.2% for lattice param-
irregularity is estimated to be ∼ 20%. The uncertainty 2 2
eter values.
in determination of the transition temperature was de-
termined by half width at half maximum for d(χT)/dT
4
B. YPd P and LaPd P their band structures and Fermi surface topologies. For
2 2 2 2
YPd P an upturn in susceptibility is present because
2 2
YPd2P2 andLaPd2P2 exhibitmagneticandelectronic of trace amounts impurities. (e.g Y1−xGdx Pd2P2 with
properties consistent with the empty 4f-shells of Y and x=0.000025 would give a comparable CW tail). The
La ions. The zero-field resistivity (ρ(T)) in Fig.3(a) magnetic signal from the addendum used for the H (cid:107)c
demonstratescharacteristicmetallicbehaviorwithanal- was subtracted for both R=Y and La.
most linear increase of the resistivity with temperature
forT >75K,withnoobservedanomaliesdownto1.8K.
The residual resistivity ratios, RRR≡ ρ(300K)/ρ(2K)), C. CePd2P2
of YPd P and LaPd P are 2.4 and 7.8 respectively.
2 2 2 2
Similar to CeAgSb [26], CePd P is a rare example
2 2 2
of a low-temperature, Ce-based ferromagnet. Although
4 0
(a )
this compound was previously reported, the work has
3 5 beendoneonpolycrystallinesamples[9–12],thereforethe
)
anisotropic characteristics of CePd P were not studied.
m 3 0 2 2
c InFig.4(a)M(T)/H datameasuredabovethesaturation
hm 2 5 field (H=10 kOe) is shown. A ferromagnetic transition
O is apparent at T < 30 K with large anisotropy between
m ( 2 0 theeasy(H(cid:107)c)andhard-axis(H(cid:107)ab)directionsofapplied
istivity 1 5 LY aP Pd d2P2P2 2 fimealdgn.eTtihcefiinelsdetosfhHow=s1MkO(Te.)/H data,measuredinalower
s 1 0 TheH/M(T)dataareplottedinFig.4(b). Byplotting
e I || (a b )
R 5 H = 0 the data in such manner, high temperature anisotropies
becomepronounced,andCWlawbehaviorismoreeasily
0 identified. For the H (cid:107)ab direction, a non-linear feature
0 5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0
is apparent in the 150 < T < 50 K region, whereas the
0
ol) - 1 0 (b ) (c ) data T > 150 K, H/Mab(T) exhibit close to linear be-
havior with a CW temperature θ =−344 K. For H (cid:107)c,
/m - 2 0 L a P d P ab
u 2 2 below 150 K, H/Mc(T) deviates from the expected lin-
-60 em -- 43 00 H = 1 0 k O e H || c YH =P 1d02 Pk O2 e eaaCrWformtemaspewrealtlu.rAebθocv=e 16540KK,wHh/icMhci(sTc)oinssilsinteenatr wwiitthh
(1 - 5 0 H || a b the axial ferromagnetic ordering. From the H/Mave(T)
/H - 6 0 H || c dmaotma,eantCoWf µtemp=era2t.4u4rµeofwθaevree=inf−er1r3edK,aslnigdhatnlyesffmecatlilveer
M - 7 0 H || a b eff B
than expected for a free Ce+3 ion (2.54µ ). The exper-
- 8 0 B
0 1 0 0 2 0 0 3 0 0 1 0 0 2 0 0 3 0 0 imental values are in good agreement with the previous
report on polycrystalline sample [12].
T e m p e r a t u r e ( K )
Figure 4(c) presents the anisotropic magnetization
isotherm data measured at T = 2.0 K. The easy
axis (H(cid:107)c) M(H) curve saturates near H = 2 kOe.
FIG. 3: (a) Zero-field, in-plane resistivity of LaPd P and
2 2
YPd2P2 (b) Anisotropic M(T)/H of LaPd2P2 measured at The saturated moment of CePd2P2 is Msat = 1.64µB,
H =10 kOe. (c) Anisotropic M(T)/H of YPd P measured higher than the reported value for polycrystalline sam-
2 2
at H =10 kOe. ples[12],yetlowerthanthepredictmoment(M(T) =
sat
g J) value of Ce+3 (2.14µ ) for free ions. The hard-
J B
Magnetization measurements performed at H = axis(H(cid:107)ab)M(H)curveontheotherhand, showsakink
10 kOe are shown in Fig. 3(b) and (c). Both com- in the linear slope at H = 18 kOe, possibly indicating
pounds present a net diamagnetic susceptibility, im- weak meta-magnetism, but does not saturate up to the
plying that the sum of Landau and core diamagnetic maximal measurement field of H =55 kOe.
contributions to the magnetic susceptibility, is greater Finally,thezero-field,in-planeresistivityasafunction
than the Pauli paramagnetic contribution. In compar- of temperature is shown in Fig. 4(d). The residual resis-
ison, a recent study has revealed that YCo Ge and tivity ratio of CePd P is RRR=12.4. The 4f electronic
2 2 2 2
LaCo Ge are Pauli paramagnetic [25]. The low γ, ∼ contributionofCePd P wasdeducedbysubtractingthe
2 2 2 2
6 mJ/K2 [9] as opposed to 10 mJ/K2 for YCo Ge and resistivity of LaPd P , which should consist primarily of
2 2 2 2
14.6 mJ/K2 LaCo Ge [25], and the diamagnetism of electron-phonon and electron-impurity scattering. The
2 2
YPd P and LaPd P imply that they have a relatively characteristicsoftheresistivitycurve,areconsistentwith
2 2 2 2
small density of states (DOS) at the Fermi surface. The previously reported Kondo-lattice behavior [9], with a
compounds show different anisotropies, probably due to weak minimum in the 4f (ρ = ρ −ρ ) component
4f Ce La
unit cell contraction, which in turn, causes changes in oftheresistivityaroundT =45K.Theferromagnetic
min
5
1 .0 2 .0
( a ) 5 ( c )
ol) 00 ..68 (/Hemu/mol)1234 H = 1 k O e (mnt /f.u.)B1 .5 -1m K)00 ..57 05 ---001 ...680 )/dT HH |||| ca b
(M/Hemu/m 00 ..24 C e P d 2P 2 M00 T1e0m p He r2 a|0|t uc r e (3K0) 4 0 netic Mome01 ..50 rmW d/dT (c00 ..02 05 0T e 1m 0p e2r0a tu3 r0e 4(K0 ) 00-- 00..20..24 d(M/Have TC =e 2P dK 2P 2
H = 1 0 k O e H || a b g
a
M
0 .0 0 .0
0 5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0 0 1 0 2 0 3 0 4 0 5 0
T e m p e r a tu r e ( K ) M a g n e tic F ie ld ( k O e )
6 0 0
( b ) H || c 7 0
( d )
H || a b
5 0 0 A v e ra g e 6 0
u) )
(H/Mmol/em234 000 000 mWistivity ( cm2345 0000 C err Pt4oftd P
s 2 2
1 0 0 C e P d 2P 2 e H = 0
H = 1 k O e R1 0 I || (a b )
0 0
0 5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0 0 5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0
T e m p e r a tu r e ( K ) T e m p e r a tu r e ( K )
FIG. 4: Measurements of CePd P . (a) Anisotropic M(T)/H measured at H = 10 kOe. Inset: Low temperature blow
2 2
up of anisotropic M(T)/H measured at H = 1 kOe (b) Anisotropic and polycrystalline averaged H/M(T) (c) Anisotropic
magnetizationisothermmeasured atT =2.0K. Inset: low-temperatureblowupof dρ/dTand d(M /H)/dT. (d)Zero-field,
ave
in-plane resistivity (blue) and the 4f electronic contribution (ρ =ρ −ρ ) (dashed magenta)(see text for details).
4f Ce La
ordering is evident with a sharp drop in the resistivity Figure 5(b) depicts the inverse magnetization. The
below the Curie temperature T , nearly to zero 4f resis- paramagnetic anisotropy is less extreme than in
C
tivity at T = 2.0 K, on account of loss of spin-disorder CePd P . Both H/M (T) and H/M (T) show linear
2 2 ab c
scattering, as the ferromagnetic order sets in. T was behavior with a CW temperature θ = −55 K and
C ab
estimated from the temperature dependence of the re- θ = 32 K. The polycrystalline averaged H/M (T) is
c ave
sistivity, by taking the derivative of ρ(T). The jump in linear with a CW temperature of θ = −8 K and an
ave
dρ/dT (insetofFig.4(c))yieldsT =28.4±0.4Kwhich effectivemomentofµ =3.52µ . Thevalueofinferred
C eff B
is consistent with the temperature at which a jump in moment is consistent with the predicted value for a free
d(M(T) /H)/dT occurs. Pr+3 ion (3.58µ ).
ave B
Figure 5(c) presents the anisotropic magnetization
isotherm data measured at T = 2.0 K, confirming the
D. PrPd P FM nature of PrPd P . The M(H) curve for the easy
2 2 2 2
axis (H(cid:107)c) rapidly saturates at H = 5 kOe. The
hard-axis (H(cid:107)ab) M(H) curve on the other hand, lin-
As shown in Fig. 5(a), PrPd P is the other ferromag-
2 2
early increases up to the maximal measurement field of
neticmemberintheRPd P series. M(T)/H data,mea-
2 2
H = 55 kOe. The saturated moment of PrPd P is
sured at H =10 kOe, above the saturation of PrPd P , 2 2
2 2 M =2.10µ , smallerthanthepredictedvalueofPr+3
are strongly anisotropic with M (T)/H (cid:29)M (T)/H in sat B
c ab
(3.2µ ) for free ions.
the ferromagnetic state. The inset shows M(T)/H data B
measured at H = 1 kOe. The easy-axis (H(cid:107)c) shows a Thezero-field,in-planeresistivity(Fig.5(d))isroughly
clear transition just below 10 K after which the magne- linear from 300 K down to around 30 K. The residual
tization saturates. resistivity ratio of PrPd P is RRR=7. The FM order-
2 2
6
2 .5
1 .2 ( a ) ( c ) P rP d 2P 2
H || c .)
1 .0 P rP d 2P 2 H || a b /f.u2 .0
B
(/Hemu/mol)000 ...468 H = 1 0 k O e (M/Hemu/mol)1234 H = 1 k O e (metic Moment 11 ..05 HH |||| ca b rmW -1d/dT (cm K)0000 ....0123 0 5 1 0 1 5 2 0 2 50--01.0..50 )d(M/H/dTave
M0 .2 00 5 1 0 1 5 2 0 2 5 gn0 .5 T e m p e ra tu re (K )
T e m p e ra tu re (K ) a T = 2 K
M
0 .0 0 .0
0 5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0 0 1 0 2 0 3 0 4 0 5 0
T e m p e r a tu r e ( K ) M a g n e tic F ie ld ( k O e )
2 5 0
( b ) P rP d P 4 0 ( d )
2 2
2 0 0
H = 1 k O e
)
) m3 0
u1 5 0 c
m W P rP d P
(mol/e1 0 0 mtivity (2 0 2 2 W cm) 667 ...050
H/M5 0 HH |||| ca b esis1 0 H = 0 rm( 55 ..05
A v e ra g e R I || (a b ) 0 T e5m p1e0ra tu1 r5e (2K 0)
0 0
0 5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0 0 5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0
T e m p e r a tu r e ( K ) T e m p e r a tu r e ( K )
FIG. 5: Measurements of PrPd P . (a) Anisotropic M(T)/H measured at H = 10 kOe. Inset: Low temperature blow up
2 2
of anisotropic M(T)/H measured at H = 1 kOe. (b) Anisotropic and polycrystalline averaged H/M(T). (c) Anisotropic
magnetization isotherm measured at T =2.0 K. Inset: low temperature blow upof dρ/dT and d(M /H)/dT. (d) Zero-field,
ave
in-plane resistivity. Inset: low-temperature blow up of ρ(T).
ing is evident with a sharp drop in the resistivity below rections are identical. For H (cid:107)c, NdPd P shows a
2 2
T ,associatedwithlossofspin-disorderscatteringwhich clear metamagnetic transition at H = 3 kOe. The
C
is clearly seen in the inset. The Curie temperature es- magnetization in H (cid:107)ab monotonically increases, crosses
timated from the jump in dρ/dT (inset of Fig 5(c)) is H (cid:107)c curve at 30 kOe and reverse the anisotropy at
T =8.4±0.3 K, in agreement with the peak tempera- H = 30 kOe. Neither for H (cid:107)ab nor for H (cid:107)c direction,
C
ture in d(M (T)/H)/dT. the M(H) curve reaches the predicted saturated moment
ave
predicted for Nd+3 (3.27µ ) up to 55 kOe.
B
The zero-field, in-plane resistivity is depicted in
E. NdPd P
2 2 Fig. 6(d). The residual resistivity ratio of NdPd P is
2 2
RRR=5.4. ρ(T) decreases linearly down around 25 K.
The M(T)/H data of NdPd P (Fig. 6(a)), measured Below25K,theresistivitystartstosaturate. At15K,a
2 2
at H = 1 kOe, are only weakly anisotropic with χ > change of slope is observed. In dρ/dT a minimum is evi-
c
χ . The M(T)/H curves follow a CW law for both dent,followedbyapeakatT =6.0±1Kasdemonstrated
ab
measurement directions. From the H/M(T) data shown intheinsetofFig5(c). Ataclosetemperature,achange
in Fig. 6(b), θ =−25 K, θ =−2 K and θ =−19 K ofslopeisobservedind(χ T)/dT. FromM(T)/H and
ab c ave ave
CW temperatures were inferred. The effective moment ρ(T) data there appears to be a weak signature which
µ =3.64µ is in good agreement with the theoretical could be associated with magnetic ordering near 6 K.
eff B
value for Nd+3 (3.62µ ). M(H) data are consistent with an AFM ground state.
B
The anisotropic magnetization isotherm data mea- Given that the rest of the RPd P compounds also or-
2 2
sured at T = 2.0 K are shown in Fig. 6(c). Up to der antiferromagnetically (see below), we propose that
H = 3 kOe, the M(H) curves of H (cid:107)c and H (cid:107)ab di- NdPd P also adopts AFM order below T ≈6 K.
2 2 N
7
2 .0
( a ) -1)0 .1 5 0 .0 4 T ( c )
u/mol) 00 ..23 HN =d P1 dk 2OP e2 (M/Hemu/mol)000 ...123 0 5 1 0 1 5 2 0 (ment /f.u.)B11 ..05 rmW d/dT (cm K00 ..01 50 T 0e m 5p e r1a0tu1re5 (2K0) 00 ..00 02 cd(T)/dave
(em T e m p e ra tu re (K ) Mom HH |||| ca b
/H 0 .1 H || c tic0 .5 N d P d P
M H || a b ne 2 2
g T = 2 K
a
0 .0 M0 .0
0 5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0 0 1 0 2 0 3 0 4 0 5 0
T e m p e r a tu r e ( K ) M a g n e tic F ie ld ( k O e )
2 0 0
7 0
( b ) ( d )
6 0
1 5 0 N d P d P
N d P d P ) 2 2
(mol/emu)1 0 0 H = 1 k 2O e2 mWity ( cm 345 000 IH |=| (0 a b ) m)1 5
H/M 5 0 HH |||| ca b esistiv 2 0 rmW ( c111 234
A v e ra g e R 1 0 0 1 0 2 0 3 0
T e m p e ra tu re (K )
0 0
0 5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0 0 5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0
T e m p e r a tu r e ( K ) T e m p e r a tu r e ( K )
FIG. 6: Measurements of NdPd P . (a) Anisotropic M(T)/H measured at H = 1 kOe. Inset: low-temperature blow up
2 2
of M(T)/H. (b) Anisotropic and polycrystalline averaged H/M(T). (c) Anisotropic magnetization isotherm measured at
T = 2.0 K. Inset: low-temperature blow up of dρ/dT and d(χ T)/dT. (d) Zero-field, in-plane resistivity. Inset: low-
ave
temperature blow up of ρ(T).
F. SmPd P (Fig. 6(c)) show a linear increase up to H = 55 kOe
2 2
with a slight anisotropy in favor of H (cid:107)c.
The zero-field, in-plane resistivity is plotted in
The M(T)/H data of SmPd P are rather dif-
2 2
Fig. 7(d). The residual resistivity ratio of SmPd P is
ferent from the previous members. As shown in 2 2
RRR=3.8. ρ(T) decreases linearly down to 25 K. Below
Fig. 7(a), the M(T)/H data are anisotropic, with
50K,theresistivitysaturates,thendecreasesslightlybe-
a change of anisotropy around 30 K (below which
low10K.AtT =3.1±0.3Kapeakindρ/dTisapparent,
M (T)/H > M (T)/H), and another change to
ab c
and coincides with the AFM transition in d(χ T)/dT.
M (T)/H >M (T)/H below7K.AtT =3.50±0.3K, ave
c ab N
an anomaly is observed for both M (T)/H, M (T)/H
ab c
and in d(χ T)/dT data (insets of Fig 7(a,c)), which
ave
likely indicates antiferromagnetic (AFM) ordering. G. EuPd P
2 2
From H/M(T) data (Fig. 7(b)), it is clear that the
paramagnetic susceptibility of SmPd P does not follow Figure 8 presents measurements done on EuPd P .
2 2 2 2
a CW law up to 300K, and appears to saturate roughly The temperature-dependent susceptibility (Fig. 8(a)), is
at room temperature. Similar behavior has been re- isotropic down to ∼ 20 K. Around 19 K, a pronounced
ported for other Sm bearing compounds [25, 27–29]. A peak indicating AFM ordering of the Eu moments is ob-
likelyexplanationforthiscouldbethermalpopulationof served. Given that M (T)/H is essentially temperature
c
the Sm3+’s Hund’s rule excited states. The anisotropic independent for 10 K below the peak, the ordered mo-
magnetization isotherm data measured at T = 2.0 K ments are likely aligned perpendicular to the crystallo-
8
0 .0 0 4 0 .0 4
( a ) 0 .0 6 0 .0 6
u/mol)00 ..00 00 23 SH m= 1P kd O2 P e2 (M/Hemu/mol)0000 ....0000 0000 2334 5050 0 2 4 6 8 1 0 (ment /f.u.)B00 ..00 23 rmW -1d/dT (cm K)000 ...000 024 0 T e2 m p4 e r6a tu8re 1 (0K )000 ...000 024 cd(T)/dTave
(em T e m p e ra tu re (K ) Mom TS =m 2 P Kd 2 P 2
/H0 .0 0 1 tic0 .0 1 ( c )
M e H || c
H || c n
g H || a b
H || a b a
0 .0 0 0 M0 .0 0
0 5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0 0 1 0 2 0 3 0 4 0 5 0
T e m p e r a tu r e ( K ) M a g n e tic F ie ld ( k O e )
1 2 0 0 ( b ) 5 0 ( d ) S m P d P
2 2
1 0 0 0
)4 0
u) 8 0 0 cm H = 0
l/em 6 0 0 mW (3 0 I || ( a b ) 1 1 .6
(mo 4 0 0 tivity2 0 cm)11 11 ..24
H/M 2 0 0 SH m= 1P kd O2 P e2 HH |||| ca b esis1 0 rmW (1 1 .0
R 0 5 1 0 1 5 2 0
A v e ra g e
T e m p e ra tu re (K )
0 0
0 5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0 0 5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0
T e m p e r a tu r e ( K ) T e m p e r a tu r e ( K )
FIG. 7: Measurements of SmPd P . (a) Anisotropic M(T)/H measured at H = 1 kOe. Inset: low-temperature blow up
2 2
of M(T)/H. (b) Anisotropic and polycrystalline averaged H/M(T). (c) Anisotropic magnetization isotherm measured at
T = 2.0 K. Inset: low-temperature blow up of dρ/dT and d(χ T)/dT. (d) Zero-field, in-plane resistivity. Inset: low-
ave
temperature blow up of ρ(T).
graphic c-axis over this temperature range. On the en- is linear down to ∼50 K. Below 25 K, at least four tran-
larged temperature scale shown in the inset of Fig. 8(a), sitions can be clearly seen. By looking closer at dρ/dT
multiple features can be identified clearly. andd(χ T)/dT(insetofFig.8(c)),acascadeoftransi-
ave
From the polycrystalline average H/M(T) shown in tionsisapparentatT1 =18.2±0.3K,T2 =12.4±0.5K,
Fig. 8(b), an effective moment µeff = 7.60µB and an T3 = 9.8±0.3 K and T4 = 5.8±0.3 K. The first tran-
average CW temperature θ = −30 K were evaluated. sition can be associated with an opening of a superzone
ave
The size of the effective moment is consistent with the gap. Such a complex magnetic ground state has been
anomalousunitcellvolumeofEuPd P ,showninTableI previously reported for CeSb [31], where six transition
2 2
and Fig. 2(b), both suggesting that Eu is in a divalent at H =0 were identified. More advanced probes will be
state. This is in agreement with M¨ossbauer [8] and pho- necessary to determine the exact nature and number of
toemission spectroscopy [30]. the observed transitions.
The anisotropic magnetization isotherm data, mea-
sured at T = 2.0 K is shown in Fig. 10(c). The M(H)
curveforH (cid:107)c directionhasalinearfielddependenceup H. GdPd2P2
to 55 kOe, whereas the M(H) curve of H (cid:107) ab reveals a
metamagnetic transition at around H = 35 kOe. Above Figure 9 shows measurements done on GdPd P . The
2 2
thetransition,H (cid:107)abandH (cid:107)cM(H)curvesmerge. At M(T)/H data (Fig. 9(a)), similarly to EuPd P , are
2 2
the maximum applied field of 55 kOe the magnetization isotopic down to 20 K. GdPd P exhibits two mag-
2 2
reaches 1.4µB, far below the theoretical value of 7µB for netic transitions around 10 K, one peak in Mab(T)/H at
a free Eu+2 ion. ∼ 11 K followed by a second peak in M (T)/H around
c
In Fig. 8(d), the zero-field resistivity of EuPd P is 7 K, as can be clearly seen in the inset of Fig. 9(a).
2 2
shown. The residual resistivity ratio is RRR = 10. ρ(T) The inverse susceptibility is depicted in Fig. 9(b).
9
0 .1 5
ol)0 .1 0 ( a ) u/mol)00 ..11 24 HE =u 1P 0d 2kP O 2 e (mnt /f.u.)B11 ..05 mW -1/dT (cm K)024 000 ...011 505 cd(T)/dTave ( c )
u/m (em0 .1 0 H = 1 k O e me rd 0 T e m 1p0e ra tu2r0e (K )3 0 E u P d P
(/Hem0 .0 5 H || c M/H0 .0 8 0 T e m p2e0ra tu re (K4 )0 etic Mo0 .5 T = 2 H 2K ||2 c
M n
H || a b ag0 .0 H || a b
M
0 .0 0
5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0 0 1 0 2 0 3 0 4 0 5 0 6 0
T e m p e r a tu r e ( K ) M a g n e tic F ie ld ( k O e )
5 0 7 0
E u P d P
( b ) 6 0 (( dd )) 2 2
4 0
E u P d 2 P 2 )5 0 H = 0
) m
u H = 1 0 k O e c I || ( a b )
m 3 0 W 4 0
(H/Mmol/e 2 0 HH |||| ca b msistivity (23 00 rmW( cm)112 050
1 0 A ve ra g e Re1 0 5
0 5 1 0 1 5 2 0 2 5 3 0
T e m p e ra tu re (K )
0 0
5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0 5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0
T e m p e r a tu r e ( K ) T e m p e r a tu r e ( K )
FIG. 8: Measurements of EuPd P . (a) Anisotropic M(T)/H measured at H = 10 kOe. Inset: low-temperature blow up
2 2
of M(T)/H measured at H = 1 kOe. (b) Anisotropic and polycrystalline averaged H/M(T). (c) Anisotropic magnetization
isothermmeasuredatT =2.0K.Inset: low-temperatureblowupofdρ/dTandd(χ T)/dT. (d)Zero-field,in-planeresistivity.
ave
Inset: low-temperature blow up of ρ(T).
H/M (T)andH/M (T),andthereforeH/M (T),are both transitions can be clearly identified in dρ/dT and
ab c ave
indistinguishableintheparamagneticstate. Theinferred d(χ T)/dT.
ave
effective moment, µ = 8.01µ , is in agreement with
eff B
the theoretical prediction for Gd+3 (7.94µ ). The aver-
B
age CW temperature is θ = −26 K, comparable with I. TbPd P
ave 2 2
EuPd P . Both µ and θ are in agreement with
2 2 eff ave
previously reported polycrystalline data [12]. The M(T)/H data of TbPd P (Fig. 10(a)) shows ex-
2 2
In Figure 9(c), the magnetization isotherm at T = tremeaxialanisotropywithχ (cid:29)χ atlowtemperature.
c ab
2.0 K is presented. The M(H) curves of the H (cid:107)ab M (T)/H is monotonically increasing with decreasing
c
and H (cid:107)c directions are isotropic up to H = 10 kOe. temperature. Below 15 K there are 3 clear transitions.
The H (cid:107)ab curve shows a metamagnetic transition at M (T)/H follows a CW law above T > 150 K. Below
ab
H = 14 kOe. At H = 30 kOe, a metamagnetic tran- 150KthereisabroadmaximuminM (T)/H,centered
ab
sition takes place for the H (cid:107)c direction. At the high- around110K,consistentathermaldepopulationofCEF
est applied field of 55 kOe the magnetic moment reaches levels leading to the strong anisotropy.
1.9µB, which is lower than the theoretical value of 7µB The H/M(T) data are plotted in Fig. 10(b). Both
for a free Gd+3 ion. H /M(T) and H/M (T) are linear above 200K with a
ab c
The zero-field, in-plane resistivity, is shown in CW temperatures θ = −107 K and θ = 43 K. An av-
ab c
Fig. 9(d). ρ(T) is metallic down to 20 K, with RRR erageCWtemperatureofθ =−21Kwithaneffective
ave
= 3. At T =10.1±0.3 K an increase in ρ(T) is evident, moment of µ = 9.71µ were inferred for polycrys-
1 eff B
which can be explained by an opening of a superzone talline averaged H/M (T), consistent with the theo-
ave
gap at the AFM transition, followed by a sharp decrease retical value for Tb+3 (9.72µ ).
B
in ρ(T) at T = 7.0±0.3 K. In the inset of Fig. 9(c), The anisotropic magnetization isotherm data were
2
10
0 .2 0
u/mol)00 ..11 05 ( a ) GH =d 1P dk 2OP e2 (M/Hemu/mol)00 ..12 80 0 1 0 2 0 (mment /f.u.)B112 ...050 rmW -1d/dT (cm K)-0000 ....3036 0 T e m5 p e r1a0tu re1 5(K )2 0000 ...123 cd(T)/dTave
m T e m p e ra tu re (K ) o G d P d P
(e M 2 2
M/H0 .0 5 H || c netic0 .5 ( c ) T = 2 H K || c
H || a b g
a H || a b
M
0 .0 0 0 .0
0 5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0 0 1 0 2 0 3 0 4 0 5 0
T e m p e r a tu r e ( K ) M a g n e tic F ie ld ( k O e )
5 0
4 0
( b ) ( d )
4 0 G d P d P G d P d P
2 2 ) 2 2
) H = 1 k O e cm3 0 H = 0
u 3 0 m
m h )1 5 .5
(mol/e 2 0 mity (O2 0 Ohm cm111 445 ...050
H/M 1 0 HHA v ||e|| rcaa b g e Resistiv1 0 I || ( a b ) rm (11 33 ..05 0 T5e m p e1 r0a tu r1e 5 (K ) 2 0
0 0
0 5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0 0 5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0
T e m p e r a tu r e ( K ) T e m p e r a tu r e ( K )
FIG. 9: Measurements of GdPd P . (a) Anisotropic M(T)/H measured at H = 1 kOe. Inset: low-temperature blow up
2 2
of M(T)/H. (b) Anisotropic and polycrystalline averaged H/M(T). (c) Anisotropic magnetization isotherm measured at
T = 2.0 K. Inset: low-temperature blow up of dρ/dT and d(χ T)/dT. (d) Zero-field, in-plane resistivity. Inset: low-
ave
temperature blow up of ρ(T).
measuredatT =2.0KandareshowninFig.10(c). The value of the resistivity is rather small, all transitions
M(H)curvefortheH (cid:107)ab directionhasafeaturelesslin- can be clearly identified in the inset. Three transitons,
ear field dependence. In the case of H (cid:107)c direction, sev- T =12.0±0.4K,T =7.5±0.4KandT =4.6±0.4K,
1 2 3
eralmetamagnetictransitionwereobservedintheM(H) canbeinferredfromdρ/dT (insetofFig10(c))andcorre-
curve, where the magnetic moment manifests character- lated with the corresponding anomalies in d(χ T)/dT.
ave
isticstep-likebehavior. Forincreasedfield,threewellde-
finedplateauswithamomentsizeof3µ atH ∼=20kOe,
B
5.5µ atH ∼=60kOeand9µ aboveH ∼=75kOewere
B B J. DyPd P
observed. Thetwolower-fieldtransitionshaveasubstan- 2 2
tial hysteresis, clearly showing they are first order. The
magnetic moment at the last plateau reaches the the- MeasurementsperformedonDyPd2P2aresummarized
oretical saturated moment for Tb+3 (9.00µB), suggest- in Fig. 11. The M(T)/H data of DyPd2P2 (Fig. 11(a))
ing that more meta-magnetic transitions are unlikely at are similar to that of TbPd2P2, but with only two ob-
higher applied magnetic fields. Similarly rich metamag- served magnetic transitions.
netic behavior had been observed in TbNi Ge [27, 32]. Theanisotropicinversemagneticsusceptibilityisplot-
2 2
ted in Fig. 11(b). The CW temperatures inferred in the
The zero-field, in-plane resistivity of TbPd2P2 is plot- paramagnetic state are θab = -37 K, θc = 16 K and θave
ted in Fig. 10(d). The residual resistivity ratio of = -11 K. The effective moment is 10.62 µB, consistent
TbPd2P2 is RRR=2.5. ρ(T) is metallic down to 30 K, with predicted value for Dy3+ (10.64 µB).
with no clear signature for loss of spin-disorder scatter- The anisotropic magnetization isotherms were mea-
ing below the highest transition temperature. A feature sured at T = 2.0 K and are depicted in Fig. 11(c). For
more consistent with an opening of a superzone gap is the H (cid:107)c M(H) curve, two meta-magnetic transitions
observed instead. Although the effect on the absolute were observed. There is a plateau with a moment size of