Pairing of weakly correlated electrons in the platinum-based centrosymmetric superconductor SrPt3P

Abstract: We report a study of the normal- and superconducting-state electronic properties of the centrosymmetric compound SrPt3P via 31P nuclear-magnetic-resonance (NMR) and magnetometry investigations. Essential features such as a sharp drop of the Knight shift at T < Tc and an exponential decrease of the NMR spin-lattice relaxation ratio 1/(T1T) below Tc are consistent with an s-wave electron pairing in SrPt3P, although a direct confirmation in the form of a Hebel-Slichter-type peak is lacking. Normal-state NMR data at T < 50 K indicate conventional features of the conduction electrons, typical of simple metals such as lithium or silver. Our data are finally compared with available NMR results for the noncentrosymmetric superconductors LaPt$_3$Si and CePt$_3$Si, which adopt similar crystal structures.

\section{Introduction}\enlargethispage{8pt}
The continuing search for new superconductors has recently resulted in the identification of a new family of
superconducting phosphide compounds with the chemical composition \textit{A}Pt$_3$P (\textit{A} = Sr, Ca, or La).\cite{Takayama12}
These materials, whose superconductivity has been claimed to be driven by conventional electron-phonon
interactions, adopt a distorted anti-perovskite structure, resembling the structure
of several noncentrosymmetric superconductors, such as LaPt$_3$Si or
CePt$_3$Si.\cite{Bauer2005}
Contrary to the latter, however,
SrPt$_3$P exhibits an inversion center due to the staggered arrangement of Pt octahedra.
Non\-centrosymmetric superconductors, characterized by antisymmetric spin-orbit couplings, are currently
the subject of intense research.\cite{Yogi04,Yogi06,Mukuda09,Bauer12}
Therefore, a comparative study of their centrosymmetric counterparts is of particular interest.
Indeed, the Sr-based compound of the new family, with a critical temperature $T_c = 8.4$\,K, has
already been investigated in some detail.
The first study included measurements of the specific heat  $C(T)$ and Hall resistivity
$\rho_\mathrm{H}(T)$.\cite{Takayama12}
The specific-heat data were interpreted as indicating a strong electron-phonon coupling with a relatively large
ratio $2 \Delta_0 / k_{\mathrm{B}}T_c \sim 5 $, and a fully-gapped excitation spectrum below $T_c$ with
a zero-temperature gap value $\Delta_0 = 1.85$ meV. The nonlinear magnetic field dependence of
Hall resistivity was attributed to the presence of multiple Fermi-surface pockets.
Less details are known on LaPt$_3$P, whose critical temperature $T_c \simeq 1.5$\,K\cite{Takayama12}
is distinctly lower than that of the Sr compound, the difference being most likely due to the
different valencies of the cations.

There exist several theoretical interpretations of the SrPt$_3$P data, which differ in their conclusions.
In Ref.~\onlinecite{Chen12},
the proximity to a dynamical charge-density wave instability is claimed to favor the occurrence of
superconductivity in SrPt$_3$P.
In another approach, based on a Mig\-dal-Elia\-sh\-berg-ty\-pe analysis, it is argued that a conventional
phonon-mediated superconductivity is observed.\cite{Subedi13}
Density-functional-theory calculations\cite{Kang13}
suggest that the onset of superconductivity is due to a strong coupling between the $pd\pi$-hybridized bands with
low-energy phonon modes confined in the $ab$ plane of the crystal lattice. The 2D character of these modes seems
essential to preserve the antipolar arrangement of the distorted Pt octahedra in SrPt$_3$P, thus enhancing both the electron-phonon coupling
constant $\lambda_{\mathrm{ep}}$ and, consequently, the critical temperature $T_c$.

The only experimental investigation at a microscopic level cited in the literature focused on muon-spin rotation
($\mu$SR) measurements.\cite{Khasanov14} The established temperature dependence of the penetration depth
in the form of $\lambda^{-2}(T)$ again indicated a single gap with a value of  $\Delta_0 = 1.55$\,meV.  However,
an upward curvature of the upper critical field $H_{c2}(T)$ just below $T_c$ was interpreted as reflecting two-band
superconducivity with equal gaps but differing values of the coherence lengths.

In this work we report on the microscopic physical properties of SrPt$_3$P,
investigated by means of nuclear magnetic resonance.\cite{Walstedt08,Curro09}
The measurements described below are based on $^{31}$P NMR, and
probe both static (line widths and -shifts), as well as dynamic (spin-lattice
relaxation) properties of the material. Overall the data indicate that
SrPt$_3$P is a rather simple metal and that the pairing configuration
is most likely of spherically-symmetric ($s$-wave) character.
A comparison with the non\-centro\-symmetric superconductors LaPt$_3$Si and CePt$_3$Si
is intended to clarify the peculiarities of the latter. The lack of inversion symmetry of the
crystal lattice may partly account for the much lower $T_c$ of LaPt$_3$Si, while the
rare-earth cations seem responsible for the unusual electronic properties of
CePt$_3$Si.\cite{Yogi04,Mukuda09}

\vspace{-5pt}
\section{Experimental details\label{sec:methods}}
Polycrystalline samples of SrPt$_3$P were synthesized under conditions of high
pressure and high temperature. High-purity (99.99\%) coarse powders of Sr, Pt,
and P were mixed in the stoichiometric 1:3:1 ratio, then thoroughly ground, and
finally enclosed in a boron-nitride container. To avoid a possible degradation
of the material due to air exposure, all the preparatory steps were done in a
glovebox under argon atmosphere. For the heat treatment, the BN crucible
was placed into a pyrophyllite cube containing a graphite heater. An external
pressure of 2 GPa was applied at room temperature and kept constant during
the two-hour ramping up of the temperature to 1050$^{\circ}$C. The final
temperature was kept constant for 40 hours and then lowered back to room
temperature within approximately one hour. After releasing the pressure, the
polycrystalline samples were extracted and kept in closed capsules filled with
Ar gas during the corresponding experiments.
Subsequently-recorded x-ray powder diffraction patterns agreed with those
reported in Ref.~\onlinecite{Takayama12}. Measurements of the low-temperature
magnetization revealed the onset of a diamagnetic response at $T_c = 8.5$\,K,
again in good agreement with the originally reported value.\cite{Takayama12}
The relevant magnetic susceptibility data $\chi(T)$ are shown in Fig.~\ref{fig:magnetization}
where, at low temperatures, we observe a significant diamagnetic response,
typical for the onset of superconductivity.\footnote{Since the magnetic response of a
superconductor is influenced by many factors, such as the sample granularity, its compactness,
flux pinning, the demagnetization factor, the remnant field, etc., the proximity
to the ideal diamagnetic response is rather coincidental in our case.}
\begin{figure}[t]
\includegraphics[width=0.85\columnwidth]{fig1}
\caption{\label{fig:magnetization}Temperature dependence of the zero field-cooled (ZFC) and field-cooled (FC) dc magnetization measured at  $\mu_0H = 0.2$ mT. The sample shows a sizeable
diamagnetic response with $T_c = 8.5$\,K.
Inset: structure of SrPt$_3$P showing the Pt$_6$ octahedra
(adapted from Ref.~\onlinecite{Kang13}). The asymmetric positions of the apical Pt(II) ions
give rise to alternating electric polarization vectors (blue arrows), hence preserving the
centrosymmetric nature of  SrPt$_3$P.
}
\end{figure}

Although all the constituent elements of SrPt$_3$P feature NMR-active nuclei, the  ${}^{31}$P nucleus is clearly
the best choice with respect to isotopic abundance (100\%) and gyromagnetic ratio ($\gamma/2\pi = 17.2356$ MHz/T).
In addition, the $I = \frac{1}{2}$ nuclear spin of  ${}^{31}$P helps avoiding
complications due to quadrupole effects.
In order to probe an adequate region of the superconducting part of the $H$-$T$ phase diagram above 1.6 K
(i.e., accessible with a $^4$He flow cryostat),
the external field for the NMR studies was set to  $\mu_0 H_{\mathrm{ext}} = 2$ T, corresponding
to a critical temperature of approximately 5\,K (see, e.g., Fig.~1 in Ref.~\onlinecite{Khasanov14}).
Analogous measurements probing the normal state in the same temperature regime were made in a
field of $\mu_0 H_{\mathrm{ext}} = 7$ T, which exceeds the zero-temperature
upper critical field $H_{c2}(0) \simeq 5.5$\,T,\cite{Khasanov14}
and, for monitoring the superconducting transition, in a 0.6-T field
across a restricted temperature range around the relevant $T_c$.

In order to achieve an acceptable signal-to-noise (S/N) ratio,
the sample was powdered to a grain size of the order of 10\,$\mu$m and kept
in the form of a loose powder, hence reducing the electrical contact between the grains.
NMR measurements at the above-mentioned external magnetic fields were made
at different temperatures between 1.6 and 300\,K by employing a conventional
phase-coherent spectrometer.
The $^{31}$P spectra were obtained by fast Fourier transformation (FFT) of the spin-echo signals
generated by $\pi/2 - \pi$ rf pulses with a 200-$\mu$s delay between pulses (the measured
transverse relaxation time is $T_2$ of ca.\ 500 $\mu$s). Since the employed ``hard'' rf pulses
($t_{\pi/2} \sim 2.2$ $\mu$s) excite non-selectively the entire resonance line, no difference
between fixed- and swept-frequency lineshapes was observed. The nuclear spin-lattice
relaxation times $T_1$ were measured at the peak position of each NMR line by means of
the inversion-recovery method, using spin-echo detection at variable delays.
The magnetic field was calibrated via ${}^{27}$Al NMR of elemental aluminum,
whose gyromagnetic ratio and Knight shift are known to high precision. This data was
subsequently used to calculate the ${}^{31}$P NMR shifts.

\vspace{-5pt}
\section{Experimental results and analysis}\label{sec:results}
\vspace{-5pt}
\subsection{NMR line shapes and shifts}
In Fig.~\ref{fig:lines} we show the evolution of the  ${}^{31}$P lines with varying
temperature at $\mu_0H = 2.0$ T, both above and below $T_c$. Similar results
with respect to the line widths and -shifts were obtained for $\mu_0H = 7.0$\,T
(see Fig.~\ref{fig:shifts_2_vs_7T}).
\begin{figure}[t]
\includegraphics[width=0.75\columnwidth]{fig2}
\caption{\label{fig:lines}Representative ${}^{31}$P NMR line shapes in SrPt$_3$P
at $\mu_0H = 2.0$ T and temperatures in the range 1.75 to 290 K. The enhanced
line widths and asymmetry below 5\,K reflect the onset of the superconducting phase.}
\end{figure}

We focus first on the normal phase.
As can be seen in Fig.~\ref{fig:shifts}, the line width, of the order of 7\,kHz at
room temperature, increases gradually by approximately 50\% down to ca.\ 12\,K
and subsequently it stays constant between 12\,K and $T_c$.
As for the line shifts in the normal phase, a moderate increase upon lowering the
temperature is observed.
This increase is mostly linear down to approximately 50\,K, below which the line
shift remains constant down to $T_c$ (see inset of Fig.~\ref{fig:shifts_2_vs_7T}).
The constancy of the shift above the superconducting transition is a feature often observed in simple metals
with an only weakly-varying electronic susceptibility $\chi_{\mathrm{P}}(T)$. The relative shift, of the
order of 0.12\%, is close to that observed in simple metals such as Ag and Na and is the same for
the 2- and 7-T external fields, respectively (see Fig.~\ref{fig:shifts_2_vs_7T}). The line widths measured in a
field of 7\,T, instead, are approximately twice as large as those recorded at $\mu_0H = 2$\,T,
which, in turn, are three times larger than those at 0.6\,T. Thus the linewidth
scales almost linearly with the applied field.

The most interesting results are obtained below $T_c$, where the inhomogeneous
field distribution in the mixed state of a type-II superconductor significantly broadens
the NMR line. Fig.~\ref{fig:shifts} shows the steep increase of the line width, from 10
to ca.\ 18 kHz, between $T_c$ and 3.5\,K (and from 3.2 to 8\,kHz at
0.6\,T - data not shown).
Upon a further reduction of $T$, the width remains constant down to the lowest
reached temperatures. At the same time the $^{31}$P Knight shift exhibits a
sudden decrease below $T_c$, a strong indication of spin-singlet superconductivity,
since the pairing of electrons with opposite spins occurring in type $s$- (or $d$-wave)
superconductors implies a significant decrease of the local spin susceptibility, as detected by
the probe nuclei. Of the two possibilities, the experimental data rule out the occurrence
of a $d$-wave SC pairing, since in that case, the drop of the Knight shift would be
less abrupt and the line broadening would reach its asymptotic $T = 0$\,K value
only gradually.\cite{Walstedt08}
From these static NMR results, an $s$-type pairing seems to be the most plausible
configuration adopted by electrons in the superconducting state of SrPt$_3$P.
\begin{figure}[thb]
\includegraphics[width=0.9\columnwidth]{fig3}
\caption{\label{fig:shifts}${}^{31}$P NMR shift (left scale) and linewidth (right scale)
vs.\ temperature, as measured at an applied field of 2\,T.
The drop in shift and the increase in width below $T_c$ are followed by a constant behavior
at intermediate temperatures and by a final decrease at higher $T$ (see text for details).
The increase of data scattering below $T_c$ reflects the decreased S/N ratio
in the superconducting phase.}
\end{figure}

\vspace{-5pt}
\subsection{NMR relaxation rates and nature of superconductivity}
The conclusion about the $s$-wave pairing configuration is also supported by the
$^{31}$P spin-lattice $1/T_1$ relaxation-rate data displayed in Fig.~\ref{fig:relaxations}.
The $1/(T_1T)$ vs.\ $T$ plot exhibits what is considered to be a typical
signature of conventional BCS superconductivity, an exponential decrease with decreasing
$T$ well within the superconducting phase, but no clear indication of a Hebel-Slichter-type
coherence peak, which is masked by a distinct increase of $1/(T_1T)$  in a narrow regime
just above $T_c$.

\begin{figure}[b]
\includegraphics[width=0.8\columnwidth]{fig4}
\caption{\label{fig:T1_curves_2T}${}^{31}$P NMR spin-lattice relaxation data of
SrPt$_3$P at selected temperatures measured at 2\,T. All fit curves exhibit a simple
exponential ($\beta = 1$) recovery of the magnetization (see text).}
\end{figure}
The spin-lattice relaxation times $T_1$ were evaluated from magnetization-recovery
curves such as those shown in Fig.~\ref{fig:T1_curves_2T}. The nuclear magnetization data
$M_z(t)$ were fitted by using the  single-exponential equation valid for spin-1/2 nuclei:
\begin{equation}
\label{eq:T1-IR}
M_z(t) = M_z^0 \left[1-f\exp(-t/T_1)^{\beta}\right].
\end{equation}
Here $M_z^0$ represents the saturation value of magnetization at thermal equilibrium,
$f$ is the inversion factor (exactly 2 for a complete inversion), and $\beta$ is a
stretching exponent which accounts for possible distributions
of relaxation rates.\cite{Shiroka2011}
Except at the lowest temperatures, the $\beta$ values were always found to be
close to 1, with no difference in fit quality whether $\beta$ was fixed or set free.
Considering the purely magnetic relaxation of $I=1/2$ nuclear probes, this implies
the same local environment for all the $^{31}$P nuclei, i.e., disorder of either
structural or magnetic nature is insignificant in our case.

\begin{figure}[t]
\includegraphics[width=0.9\columnwidth]{fig5}
\caption{\label{fig:relaxations}Temperature dependence of $(T_1T)^{-1}$ measured
in a magnetic field of 2\,T. A similar cusp-like feature at $T_c$ is seen also at 0.6\,T.
The decrease below the maximum occurring at $T_c$
is compatible with
a fully-gapped superconductor (see text). Inset:  $(T_1T)^{-1}$ data
over the full temperature range show a broad maximum at $T^*=45$\,K
(see also Fig.~\ref{fig:suscept}).}
\end{figure}

Up to $T^* = 45$\,K, the  normal-state spin-lattice relaxation rate varies approximately linearly with
temperature (see inset in Fig.~\ref{fig:relaxations}), as is usually observed
in conventional simple metals.
The progressive deviation from the linear-in-$T$ behavior above $T^*$ reflects a
decrease in the dynamical spin susceptibility, as discussed in Sec.~\ref{ssec:korringa}
and \ref{ssec:susceptibility}.

Focusing back on the superconducting state, just below $T_c$ we do not
observe a distinct Hebel-Slichter (coherence) peak, but only a cusp-like feature peaking at $T_c$.
Analogous data taken in a 0.6-T field showed essentially the same type of feature in the vicinity
of $T_c$.
We recall that
the Hebel-Slichter anomaly is due to pair coherence and to a quasiparticle density-of-states anomaly
at the gap edge and is more pronounced at lower fields, where the pair-breaking effect of the field
is less important.\cite{Masuda69} %
Both these effects can be weakened by different causes, including interaction
anisotropies and the broadening of the quasiparticle states due to strong-coupling effects.
Indeed, while \textit{weakly}-coupled conventional $s$-wave superconductors,
typically exhibit a Hebel-Slichter peak, this is not the case in the \textit{strong}-coupling limit,
where the enhanced electron-phonon coupling implies reduced quasiparticle lifetimes.\cite{Ohsugi92}

Since the lack of a coherence peak cannot rule out
the possibility of an $s$-type electron pairing,\cite{Parker08} further analysis is required.
The temperature dependence of $T_1$ at $T < T_c$ can provide information about the symmetry
of the superconducting gap. In case of an anisotropic gap with nodes, a power-law dependence, $1/T_1 \sim T^{n}$,
is expected since very low-energy electronic excitations around gap nodes can still contribute to relaxation.
On the other hand, for standard fully-gapped superconductors with no or a weak gap anisotropy,
the relaxation rate depends
exponentially on temperature, $1/T_1 \sim \exp(-\Delta_0/T)$, with $\Delta_0$ the
superconducting gap at $T = 0$.
In either case, non-negligible $1/(T_1T)$ residual values may be observed at the lowest
temperatures.\footnote{Since samples of different quality show similar
residual values, the observed offset in $1/(T_1T)$ is likely of intrinsic origin.}
Most likely these are due to the presence of vortex-core or of thermally-excited vortex-motion
relaxation mechanisms,\cite{MacLaughlin76,Corti96} which can be accounted for by an
additional constant term.
While both models provide reasonable fits to the data, the power-law expression
shows a somehow poorer fit quality and, most importantly, the resulting exponent varies significantly
(from 3.6 for the 2-T dataset to 2.7 for the data collected at 0.6\,T).
On the other hand, as shown in  Fig.~\ref{fig:relaxations}, an exponential curve with
$\Delta_0 = 1.47(4)$\,meV fits the SrPt$_3$P relaxation data satisfactorily well. Also the data
taken at 0.6\,T can be consistently fitted with a similar gap value, $\Delta_0 = 1.51(4)$\,meV,
thus providing a good indication
that the superconducting gap of SrPt$_3$P is nodeless.
The value of the gap parameter,
although slightly modified by the presence of the constant term,
is in fair agreement with the 1.55\,meV value reported in
Ref.~\onlinecite{Khasanov14}, but distinctly lower than $\Delta_0 = 1.85$\,meV reported in
Ref.~\onlinecite{Takayama12}.
In any case, since the gap value from our analysis is clearly higher than 1.27 meV
(the gap expected from the BCS formula $2 \Delta_0 = 3.52 k_{\mathrm{B}}T_c$),
it is compatible with a strong-coupling scenario.

Besides the vortex-core relaxation mechanism mentioned above,
the almost constant value
of $1/(T_1T)$ well below $T_c$
might also reflect a certain degree of superconducting-gap anisotropy, but we
consider this as a less likely possibility.

\vspace{-5pt}
\subsection{Korringa relation and the degree of electronic correlation\label{ssec:korringa}}
Some insight into the degree of electronic correlations in the normal state of a
material can be gained by considering the so-called Korringa relation\cite{Korringa50}
\begin{equation}
\label{eq:Korringa_r}
T_1TK_s^2 = S_0, \quad \text{with} \quad
S_0 = \frac{\gamma_e^2}{\gamma_n^2} \frac{\hbar}{4 \pi k_{\mathrm{B}}},
\end{equation}
Here $\gamma_e$ and $\gamma_n$ are the electronic and the nuclear gyromagnetic ratios, respectively,
while $K_s$ is the line shift due to the polarized conduction electrons. Eq.~(\ref{eq:Korringa_r})
reflects the fact that both the readily accessible experimental parameters,
the Knight shift $K_s$ and the spin-lattice relaxation time $T_1$, depend on the same
electron-nucleus hyperfine interaction (assumed to be mostly of Fermi-contact character).
However, while for a simple (Fermi-gas) metal the parameter $S_0$ is a constant (depending
only on $\gamma_n$ of the probe-nucleus), for real materials the generalized
relation\cite{Narath68} $T_1TK_s^2 = \alpha S_0$, with $\alpha$ a measure of
the strength of the (Fermi-liquid) quasiparticle interactions,\cite{Pines54} is more appropriate.
Since $1/(T_1T)$ probes the dynamical susceptibility averaged over the Brillouin zone, either
FM or AFM spin correlations will enhance it, but only FM correlations can significantly enhance
also the shift.
Thus, $\alpha > 1$ ($< 1$) indicates the presence of ferro- (antiferro-)magnetic electronic correlations.

From the almost constant value of $S_0/(T_1TK_s^2)$ vs.\ $T$
(see inset in Fig.~\ref{fig:relaxations}) and
considering the negligible orbital shift contribution (see below), we find
$1/\alpha \simeq 0.55$ and, hence, $\alpha = 1.8$ in our case.
The resulting value, slightly larger than unity, seems to suggest
the presence of ferromagnetic correlations in the normal state. However, the ratio
$1/\alpha \simeq 0.55$ is similar to 0.58, 053, and 0.50 observed in typical
metals such as lithium, cesium, and silver, respectively.\cite{Walstedt08} In these nearly-free-electron
metals the Moriya theory of exchange enhancement, whereby the exchange fluctuations
also enhance the rate of the $T_1$ process,\cite{Moriya63,Walstedt08} explains the
experimental data satisfactorily. Consequently, we conclude that also in SrPt$_3$P the conduction
electrons experience rather weak correlation effects, in good agreement with
the results reported in Ref.~\onlinecite{Takayama12}.

\begin{figure}[t]
\includegraphics[width=0.9\columnwidth]{fig6}
\caption{\label{fig:korringa}Korringa relation [reported in
$S_0$ units, see Eq.~(\ref{eq:Korringa_r})] shows only a weak temperature dependence,
with no distinct features at $T^*=45$\,K (see also Fig.~\ref{fig:suscept}).}
\end{figure}

\vspace{-5pt}
\subsection{Magnetic susceptibility and NMR\label{ssec:susceptibility}}
The magnetization $M(T)$ was measured in an applied magnetic field of 2 T in specimens from the
same batch as that used for NMR experiments. The resulting values of the susceptibility,
$\chi = M/H$, for temperatures above 10\,K are shown in the main panel of Fig.~\ref{fig:suscept}.
The measured net susceptibility seems unusually small for a metal, because the contribution
from the conduction electrons is almost exactly canceled by the large core diamagnetism mainly
from Pt ($1.18\times10^{-4}$\,emu/mol per formula unit, according to standard tables). This
compensation effect makes it possible to measure even subtle variations of the electronic
contribution with $T$.

In order to relate $\chi(T)$ with the NMR data, the relative ${}^{31}$P
NMR shifts in the two chosen magnetic fields (2 and 7\,T, respectively)
are shown in Fig.~\ref{fig:shifts_2_vs_7T}.
\begin{figure}[t]
\includegraphics[width=0.9\columnwidth]{fig7}
\caption{\label{fig:suscept}Temperature dependence of the SrPt$_3$P molar magnetic
susceptibility $\chi_{\mathrm{mol}}$ in an applied field of 2\,T. Data were corrected by the
addition of $1.18\times10^{-4}$\,emu/mol to account for the core diamagnetism
(see text). Inset: ${}^{31}$P Clogston-Jaccarino ($K$-$\chi$) plot at
$\mu_0 H = 2$\,T. The line is a fit to the $T > 50$\,K data.}
\end{figure}
As expected for simple metals, the shifts are independent of the applied field and hence
coincide perfectly. Below $T_c$ the high-field (7-T) data reflect the normal phase: both
the line shift and width remain practically constant upon decreasing temperature.
\begin{figure}[t]
\includegraphics[width=0.8\columnwidth]{fig8}
\caption{\label{fig:shifts_2_vs_7T}Comparison of ${}^{31}$P NMR line shifts vs.\
temperature at 2 (circles) and 7\,T (squares) on logarithmic (main panel) and on
linear scale (inset). The nature of the shift changes from constant to linearly
decreasing at 45 K (arrows). The dashed lines highlight both regimes and serve
as guides to the eye.}
\end{figure}

From the $\chi(T)$ data (shown in Fig.~\ref{fig:suscept}) and the NMR shifts (see inset
of Fig.~\ref{fig:shifts_2_vs_7T}), one recognizes immediately very similar features.
In either case, a practically linear increase upon lowering the temperature,
from 300\,K down to ca.\ $T^* = 45$\,K, is followed by a $T$-independent part.
This clear break in the temperature dependence of both $K_s(T)$ and $\chi(T)$
coincides with a maximum of the Hall constant $R_{\mathrm{H}}$ reported in Ref.~\onlinecite{Takayama12}.
The reason for this might be that the trend to $R_{\mathrm{H}} = -1/ne$ starts to dominate
the $R_{\mathrm{H}} > 0$ contribution.

Such a comparison of micro- (NMR) vs.\ macroscopic (magnetometry) results provides
further information about the normal-state electronic properties of SrPt$_3$P. In our case, given
the relative insensitivity to the applied field (see Fig.~\ref{fig:shifts_2_vs_7T}),
it makes sense to compare the NMR shifts with the susceptibility data, as shown in
the Clogston-Jaccarino\cite{Clogston64} (or $K$-$\chi$) plot in the inset of
Fig.~\ref{fig:suscept}. Here both data sets refer to the same (2-T) field, with
the temperature as an implicit parameter.
From the slope $\mathrm{d}K/\mathrm{d}\chi$ of the $K$-$\chi$ plot for $T > 50$\,K,
and by assuming a standard electronic $g$-factor of 2.0, the relation
$K = g A_{\mathrm{hf}} \chi + K_{\mathrm{orb}}$ gives the hyperfine coupling
constant, $A_{\mathrm{hf}} = 4.4$\,T/$\mu_{\mathrm{B}}$,
while the zero intercept with the $K$-axis gives the orbital shift, $K_{\mathrm{orb}} = -0.003$\%.
The latter is perfectly compatible with known orbital-shift values for
a light nucleus such as ${}^{31}$P.
Clearly, the orbital shift is negligible in comparison with the dominant magnetic shift.
We recall that hyperfine coupling constants are known to be strongly material-dependent, so it
is not unusual that our value is an order of magnitude higher than, e.g.,
$A_{\mathrm{hf}} \sim 0.2$\,T/$\mu_{\mathrm{B}}$ found in another
phosphorus-based compound, Pb$_2$(VO)(PO$_4$)$_2$.\cite{Nath09}

\vspace{-5pt}
\section{Discussion}\label{sec:discussion}
Below we discuss in more detail some aspects of the results presented in the previous sections.

\emph{Line shapes.}
For all the applied magnetic fields the normal-state ${}^{31}$P lines show a slight
asymmetry, usually attributed to the crystallite surfaces or to minor impurities.
However, the persistence of the line asymmetry even in samples of different quality,
most likely indicates a (uniaxial) anisotropy of the hyperfine coupling tensor
(whose study is beyond the scope of this work), which can be
observed also in powdered samples.\cite{Kawamoto95}
The uniaxial nature of the HF tensor is to be expected, if the symmetry
of the distorted Pt octahedra probed by the phosphorus nuclei is considered
(see structural unit in the inset of Fig.~\ref{fig:magnetization}).

The line shapes, in either the normal or the superconducting phase, indicate
also the structural and chemical homogeneity
of the sample.
An essentially single line with a relatively narrow width of only
$\sim 7$\,kHz (see Fig.~\ref{fig:lines}) is typical of defect-free crystal structures,
especially when considering the random polycrystalline nature of the sample.
Similarly, the line asymmetry observed below $T_c$, again reflects the high
quality of the sample, since in a powdered superconductor with a significant
degree of disorder the lines would appear symmetric.

\emph{Residual shift.}
A key prediction of the BCS theory for superconductors with singlet $s$-wave pairing
is the vanishing of the Knight shift for $T \ll T_c$, reflecting the vanishing
susceptibility due to the pairing of quasiparticles with antiparallel spins. Yet, even for elemental superconductors
(as, e.g., Sn or Hg), a residual shift is observed even close to $T = 0$.\cite{MacLaughlin76}
The presence of a significant orbital Knight shift and/or of spin-orbit scattering (known as Ferrell-Anderson effect)
have been identified
as responsible for its occurrence.
The SrPt$_3$P case, however, is a bit more complex, since here $^{31}$P is a low-$Z$
nucleus
which interacts with superconducting carriers arising mostly from in-plane
$pd\pi$-hybridized states between Pt and P ions.\cite{Kang13}
Since orbital shifts are negligible for the light probe-nuclei,
the presence of a low-temperature nonzero Knight shift in SrPt$_3$P
is, therefore, most likely due to spin-orbit scattering effects.

\emph{Core-polarization effects and $T^*$ anomaly.}
Another contribution to the nonvanishing shift at $T = 0$ can arise from
core-polarization effects,\cite{Freeman65,MacLaughlin76} whereby the exchange
interaction of localized core electrons with field-polarized conduction band electrons
(partially of $d$ character in our case), implies an additional nonzero contact interaction
in the $^{31}$P nuclei.
The observation of large negative  $^{195}$Pt-NMR Knight shifts (data not shown)
seems to confirm the presence of core-polarization effects in SrPt$_3$P.\cite{Clogston64}
While this contribution to the Knight shift would disappear entirely if purely $d$-band electrons
were coupled into Cooper pairs, a residual shift would still arise considering the hybridized
(i.e., not exclusively $d$) character of the conduction band quasiparticles.

The consequences of core-polarization effects are more important, however, in the
normal state. The electronic contribution to the magnetic susceptibility of a $d$-type
band is generally temperature dependent, thus implying a temperature-dependent
Knight shift, as actually observed (for an ideal metal
one would have expected a constant shift). While core-polarization effects might explain
the similar $T$-dependence of $K$ and$\chi_s$, the abrupt change in slope observed at $T^*$,
as well as the value of $T^*$ itself, remain to be investigated.

\emph{Weak electron correlation.}
Our NMR data of SrPt$_3$P, indicate a fairly simple metallic behavior in the normal state,
followed by standard BCS features in the superconducting state. These results are
in remarkable agreement with the transport and magnetometry measurements reported in
Ref.~\onlinecite{Takayama12} in two important aspects.
Firstly, both NMR (via the Korringa product) and the macroscopic measurements (via the
Sommerfeld-Wilson ratio\footnote{The Sommerfeld-Wilson ratio, the dimensionless ratio
of the zero-temperature magnetic susceptibility and the linearly varying term of the
temperature dependence of specific heat, can provide insight into the degree of
electronic correlation.})
demonstrate the \textit{absence of significant electron correlations}.
Secondly, there is a distinct feature at $T^* = 45$\,K in the $^{31}$P Knight
shift and in the $1/(T_1T)$ relaxation behavior, respectively, both of which
change from a constant-in-$T$ to a linear-in-$T$ decrease as the temperature
changes across $T^*$.
This temperature coincides with the above-mentioned maximum in
$R_{\mathrm{H}}(T)$, the latter reflecting anisotropies of the electron mean
free path across different parts of the Fermi surface.
The latter fact is consistent with conclusions drawn in Ref.~\onlinecite{Khasanov14}.
The coincidence of $T^*$ values, derived from results of very different types of
experiments, is unclear at present.

As for the Korringa relation, we recall that the absolute value of $\alpha$ might also
be influenced by extraneous factors, such as the presence of strong disorder,\cite{Shastry94}
or effects related to the hyperfine form factor (i.e., $q$-space filtering). Nevertheless, as
mentioned above, these can be excluded in our case and the observation of a constant
Korringa ratio vs.\ temperature represents a valid proof of Fermi-liquid behavior in SrPt$_3$P,
where $1/\alpha \simeq 0.55$ is constant over a fairly wide temperature range.

\emph{Centro- vs.\ noncentro-symmetric superconductors.}
We conclude with a comparison of SrPt$_3$P with two noncentrosymmetric superconductors,
the non-magnetic LaPt$_3$Si ($T_c = 0.6$\,K) and the magnetic CePt$_3$Si ($T_c = 0.75$\,K),
both of which were studied via NMR, as well.\cite{Yogi04,Yogi06,Mukuda09}
The latter not only share the same structure, but show also similar antisymmetric spin-orbit couplings.
Hence, a comparison of the three is relevant to understand whether it is the lack
of inversion symmetry, or the presence of rare-earth ions to account for the unusual
properties of noncentrosymmetric superconductors.

While SrPt$_3$P and LaPt$_3$Si both follow a perfect Korringa law down to
$T_c$, in CePt$_3$Si the $1/(T_1T)$ ratio is enhanced upon cooling due to
the development of 4$f$-derived magnetic fluctuations.
Most importantly, unlike in SrPt$_3$P and LaPt$_3$Si, where the superconducting
transition is of conventional type and not dominated by correlation effects, in
CePt$_3$Si, the superconductivity emerges from a unique heavy-electron state
which coexists with antiferromagnetic order.

From these facts, we conclude that the key differences between the considered
centro- and non-centrosymmetric superconductors, are most likely related to
the presence of Ce$^{3+}$ ions and the ensuing heavy-fermion phenomena,
while the parity mixing alone (common to both LaPt$_3$Si and CePt$_3$Si)
cannot explain the peculiar behavior of the latter.
Despite the many similarities between SrPt$_3$P and LaPt$_3$Si, their belonging
to different classes of crystal symmetry may possibly explain their distinctly
different critical temperatures.
Incidentally, from the available data in literature, also the centrosymmetric LaPt$_3$P
compound seems most likely to be a conventional BCS type-I superconductor.

\vspace{-5pt}
\section{Conclusions\label{sec:conclusions}}
We presented the results of ${}^{31}$P NMR measurements in the recently discovered
SrPt$_3$P superconducting compound. In the normal state, the nuclear spin-relaxation
rate obeys a Korringa-type relation $(T_1T)^{-1}_n = 0.522$ (sK)$^{-1}$ and
the line shift is independent of temperature. In the superconducting state,
the lack of a clear coherence peak in $T_1^{-1}$ just below $T_c$,
most likely due to a strong electron-phonon coupling, is followed by an exponential,
thermally-activated decrease upon further cooling. A similar behavior is shown also
by the line shift.
The overall results of our NMR investigation make us conclude that SrPt$_3$P, which
behaves as a standard metal in the normal state, is characterized by singlet pairing
of the electrons, and a gap function with conventional $s$-wave-type symmetry.

Finally, a comparison with the analogous non-centrosymmetric superconductors
LaPt$_3$Si and CePt$_3$Si, strongly suggests that the peculiarities of the latter
are closely related to the presence of rare-earth cations. The lack of inversion
symmetry may cause a much lower $T_c$, however.

\vspace{-5pt}
\section*{Acknowledgments}
\vspace{-5pt}

The authors thank R.\ Khasanov (Paul Scherrer Institut) for sharing his
data prior to publication and for useful discussions.
This work was financially supported in part by the Schweizerische Nationalfonds
zur F\"{o}rderung der Wissenschaftlichen Forschung (SNF) and the NCCR research
pool MaNEP of SNF.


Abstract: The crystal and magnetic structure of the magnetocaloric compound FeMnP0.5Si0.5 have been studied by means of neutron and X-ray powder diffraction. Single phase samples of nominal composition FeMnP0.5Si0.5 have been prepared by the drop synthesis method. The compound crystallizes in the Fe2P-type structure (P-62m) with the magnetic moments aligned along the aaxis. It is found that the Fe atoms are mainly situated in the tetrahedral 3g site while the Mn atoms prefer the pyramidal 3f position. The material is ferromagnetic (TC = 382 K) and at 296 K the total magnetic moment is 4.4 \mu B/f.u. It is shown that the magnetic moment in the 3f site is larger (2.5 \mu B) than in the 3g site (1.9 \mu B).

\section{Introduction}
\label{}
Magnetocaloric compounds have gained an increased interest since the middle 1990's due to environmental and energy benefits from magnetic refrigeration and the discovery of the Giant Magnetocaloric Effect (GMCE) \cite{gme_1}. In GMCE compounds rare earth metals are common alloying elements and due to their rareness in nature and high costs, they will be a problem in a future large-scale production. Hence, there is a need to find more common and cheaper compounds not based on rare earth metals that possess the GMCE.

A large number of studies on the compound Fe$_2$P have been performed since the 1960's and its crystallographic and magnetic properties have been well investigated \cite{rundqvist_59, wappling2, lundgren_fe2p}.  The relatively large saturation magnetization, first order nature of the transition and readily tunable transition temperature with various substitutions make the Fe$_2$P system a candidate compound for magnetocaloric applications. Numerous compounds based on Fe$_2$P have been fabricated during the years and compounds of the FeMnP$_{1-x}$M$_x$-type (M = Si, Ge and/or As) have shown improved  magnetocaloric properties.

In this investigation, samples of nominal composition FeMnP$_{0.5}$Si$_{0.5}$ have been synthesized and the crystal and magnetic structure and the magnitude of the magnetic moments at different temperatures have been determined.  There is a controversy about the exact properties of the FeMnP$_{1-x}$Si$_x$-system.  A previous study by Cam Thanh et al. \cite{cam_tanh_1} on samples of nominally the same composition, FeMnP$_{0.5}$Si$_{0.5}$, prepared by a ball milling technique, had a transition temperature of 332 K with higher preserved magnetocaloric effect.   There is a significant difference in the structural and magnetic properties presented by Cam Thanh et al. \cite{cam_tanh_1} compared to our results.  This dichotomy may be explained by the fact that their samples contained about 18\% of a second phase, identified as Fe$_{2}$MnSi.

Our magnetization measurements on the FeMnP$_{0.5}$Si$_{0.5}$  sample show a ferromagnetic transition temperature of 382 K and the magnetic entropy change $-\Delta S_M$ has from magnetization experiments been estimated to about 8 J/kgK in a magnetic field change of 1.8 T \cite{mossbauer_uu}. The magnetocaloric effect is slightly lower compared to similar compounds \cite{bruck_fe2p_comp, thanh_mnfepsige} but due to its content of only common, non-toxic, elements FeMnP$_{0.5}$Si$_{0.5}$ is a promising magnetocaloric compound within the  FeMnP$_{1-x}$M$_x$ alloy system.

\section{Experiments}
\label{}
\subsection{Sample preparation}
\label{}
FeMnP$_{0.5}$Si$_{0.5}$ samples were prepared by the drop synthesis method \cite{hf_ugn} using a high frequency induction furnace at 1623-1673 K in an Ar atmosphere of 40 kPa. Stoichiometric amounts of iron (Leico Industries, purity 99.995\%. Surface oxides were reduced in H$_2$-gas.),  manganese (Institute of Physics, Polish Academy of Sciences, purity 99.999\%),  phosphorus (Cerac, purity 99.999\%) and silicon (Highways International, purity 99.999\%) were used as raw materials. All samples were crushed, pressed into pellets and sealed in evacuated fused silica tubes. Subsequently, the samples were sintered at 1373 K for 1 h, annealed at 1073 K for 65 h and finally quenched in cold water.

\subsection{X-ray powder diffraction}
Phase analysis and crystal structure characterizations were performed using X-ray powder diffraction (XRD) with a Bruker D8 diffractometer equipped with a V\aa ntec position sensitive detector (PSD, 4 degree opening) using Cu K$\alpha_1$ radiation, $\lambda$ = 1.540598 \AA. The measurements were made using a 2$\theta$-range of 20-90\textdegree\ at  296 K and 403 K and a  2$\theta$-range of 35-60\textdegree\ in the temperature range 373 K to 393 K.

\subsection{Neutron powder diffraction}
\label{}
Neutron powder diffraction data were collected on the instrument MEREDIT at the Nuclear Physics Institute in Rez, Czech Republic. The neutron beam was monochromatised by a copper mosaic monochromator (reflection 220) giving a wavelength of $\lambda$ = 1.46 \AA. Samples were studied in a 2$\theta$-range of 4-148\textdegree\ at 296 K and 450 K.

\subsection{Refinements of the crystal and magnetic structure}
Structure refinements were performed on the neutron powder profiles by the Rietveld method \cite{rietveld} using the software \textsc{fullprof} \cite{fullprof} and unit cell parameters from XRD data were refined using the software \textsc{unitcell} \cite{unitcell}. The neutron wavelength was refined by using unit cell parameters determined from XRD data as starting point. The refined wavelength was used to refine the unit cell parameters at 450 K from neutron powder diffraction data. The peak shape was described by a psuedo-Voight profile function and the background was determined by a linear interpolation between chosen points. The following parameters in the 450 K data set were varied: peak shape, unit cell parameters, scale factor, half width parameters, zero point, background, atom occupancies, isotropic temperature parameter and atomic coordinates. The same parameters of the 296 K data set were varied  as well as the parameters for the magnetic moments. The magnetic form factors of Fe and Mn were set as shown in Refs. \cite{fuji-79, takei-63}, respectively. The occupancies of the P/Si sites were kept at a 50/50 ratio, thus not allowing the actual P and Si content to be determined from these refinements.

\section{Results}
\label{results}

\subsection{Phase analysis and crystal structure}
\label{structure}
 The XRD investigation confirms that FeMnP$_{0.5}$Si$_{0.5}$ crystallizes in the hexagonal Fe$_2$P-type structure, space group P\=62m and unit cell parameters a=6.2090(3) \AA, c=3.2880(3) \AA. The XRD pattern for FeMnP$_{0.5}$Si$_{0.5}$ at 296 K is shown in Figure \ref{xray1} which reveals a pattern of a single phase sample.  XRD-patterns in the range 363 to 403 K are shown in Figure \ref{xray_thermo} where it can be seen that FeMnP$_{0.5}$Si$_{0.5}$  undergoes a structural transition (within the space group) between $\sim$373 K and 393 K. The a-axis has decreased $\sim$2\% while the c-axis has increased $\sim$5\% compared to 296 K, see Table \ref{table:fe2p}. The factor c/a and the volume have increased $\sim$7\% and 1\% respectively. The structural transition occurs in the same region as the Curie temperature why it is likely that the transition originate from magnetostriction effects. Structure refinements of powder neutron diffraction data shows that the Fe and Mn atoms are preferably situated in the 3g and 3f sites respectively which are based on the interatomic distances (see Table \ref{table:inter-atom-dist}) and the refined occupancies of the Fe and Mn atoms (see Table \ref{table:occupancy}).

The composition based on the refined occupancies extracted from the neutron powder diffraction data  (however assuming fixed ratio P/Si=1) indicates that the acquired composition of the metallic atoms in the sample is (Fe$_{1.018(3)}$Mn$_{0.982(3)}$P$_{0.5}$Si$_{0.5}$. Also, the synthetic process with stoichiometric amounts of the raw materials showed minor (less than 0.5\%) losses.

\begin{figure}[t!]
\centering
\includegraphics[width = 75mm]{figure1.eps}
\caption{X-ray powder diffraction pattern of FeMnP$_{0.5}$Si$_{0.5}$ at 296 K. The tick marks indicate the Bragg positions of FeMnP$_{0.5}$Si$_{0.5}$. $\lambda$ = 1.540598 \AA.}
\label{xray1}
\end{figure}

\begin{table}[b!]
\centering
\caption{Unit cell parameters of  FeMnP$_{0.5}$Si$_{0.5}$ at 296 K and 450 K refined from XRD data, $\lambda$ = 1.540598 \AA.}
\begin{tabular}{ l c c c c}
\toprule
T (K) & a (\AA)  & c (\AA) & V (\AA$^3$) & c/a \\
\midrule
296 & 6.2090(3)  & 3.2880(2) & 109.78(2)  & 0.5296(1) \\
450 & 6.0830(8)  & 3.4507(9) & 110.58(4)  & 0.5672(1) \\
\bottomrule
\label{table:fe2p}
\end{tabular}
\end{table}

\begin{table}[b!]
\centering
\caption{Interatomic distances in  FeMnP$_{0.5}$Si$_{0.5}$ at 296 K.}
\begin{tabular}{l l l l l l}
\toprule
Atoms & & Distance (\AA) & Atoms & & Distance (\AA) \\
\midrule
Fe(1) -- & 2 P/Si(2)  & 2.292(2) & Mn(2) -- & 1 P/Si(2)  & 2.502(6) \\
&  2 P/Si(1)&  2.342(2) &  & 4 P/Si(1)  & 2.506(2) \\
& 2 Mn(2)  & 2.675(5) &  & 2 Fe(1)  & 2.675(5) \\
& 4 Mn(2)  & 2.742(1) &  & 4 Fe(1)  & 2.742(4) \\
& 2 Fe(1)  & 2.766(3) &  & 4 Mn(2)  & 3.275(6) \\
\\
P/Si(1) --  & 3 Fe(1)  & 2.342(2) & P/Si(2)  -- & 6 Fe(1)  & 2.292(2) \\
& 6 Mn(2)  & 2.506(5) &  & 3 Mn(2)  & 2.502(3) \\
\bottomrule
\label{table:inter-atom-dist}
\end{tabular}
\end{table}

\begin{figure}[b!]
  \centering
\includegraphics[width = 60mm]{figure2.eps}
\caption{XRD-patterns of FeMnP$_{0.5}$Si$_{0.5}$  showing the structural transition occurring at approx. 385 K. $\lambda$ = 1.540598 \AA.}
\label{xray_thermo}
\end{figure}

\begin{figure}[t]
\centering
\includegraphics[width = 70mm]{figure3.eps}
\caption{Structure refinements from neutron powder diffraction data of  FeMnP$_{0.5}$Si$_{0.5}$ at 296 K and 450 K. The differences in intensity and position of the Bragg peaks are due to the ferro- to paramagnetic transition at $\sim$390 K and the structural transition at $\sim$385 K. The peaks with the highest magnetic intensity are marked with an arrow.}
\label{neutron_diff}
\end{figure}

\subsection{Magnetic structure}
Refinements of the neutron powder diffraction data at 296 K and 450 K are shown in Figure \ref{neutron_diff}. The magnetic contribution from the ferromagnetic phase at 296 K is distinguished by the higher peak-intensities at lower 2-theta angles. The magnetic moments are aligned in the a-direction and are 1.9(1) $\mu_B$ and 2.5(1) $\mu_B$ in the M(1) and M(2) site respectively which gives a total magnetic moment of 4.4(2) $\mu_B$. The magnetic symmetry group was found to be Cm2m.

\begin{table*}[t!]
  \centering
  \caption{Placement and occupancy of the Fe and Mn atoms in FeMnP$_{0.5}$Si$_{0.5}$ at 296 K and 450 K. Derived from refinements of neutron powder diffraction data. (0, 0, 0) was chosen as an origin. }
\begin{tabular}{lllllllllllll}
\toprule
 &  &  & 296 K & & & & & 450 K & & &\\
Atom & Site & & x & y & z &  Occ. &  & x & y & z &  Occ.\\
\midrule
Fe(1) & 3g & & 0.2572(4) & 0 & 1/2  & 0.2353(4) & & 0.2548(3) & 0 & 1/2   & 0.2353(4)\\
Mn(1) & 3g & & 0.2572(4) & 0 & 1/2  & 0.0147(4) & & 0.2548(3) & 0 & 1/2  & 0.0147(4)\\
Fe(2) & 3f & & 0.597(1) & 0 & 0  & 0.0193(2) & & 0.591(1) & 0 & 0  & 0.0193(2) \\
Mn(2) & 3f & & 0.597(1) & 0 & 0  &  0.2307(2) & & 0.591(1) & 0 & 0  & 0.2307(2) \\
P/Si(1) & 2d & & 1/3 & 2/3 & 1/2  & 0.1667 & & 1/3 & 2/3 & 1/2 & 0.1667\\
P/Si(2) & 1a & & 0 & 0 & 0  & 0.0833 & & 0 & 0 & 0 &  0.0833 \\
\midrule
 & & & \multicolumn{4}{l}{R$_p$ = 2.37\%, R$_{wp}$ = 3.05\%, $\chi^2$ = 3.03} & & \multicolumn{4}{l}{R$_p$ = 2.59\%, R$_{wp}$ = 3.36\%, $\chi^2$ = 3.44} & \\
 & & & \multicolumn{4}{l}{R$_{Bragg}$ = 5.26\%, R$_{mag}$ = 6.08\% }  & & \multicolumn{4}{l}{R$_{Bragg}$ = 4.53\%} & \\
\bottomrule
\label{table:occupancy}
\end{tabular}
\end{table*}

\section{Discussion}
\label{}

\begin{figure}[b!]
\centering
\includegraphics[width = 80mm]{figure4.eps}
\caption{The magnetic structure of FeMnP$_{0.5}$Si$_{0.5}$, a), and Fe$_2$P, b). The magnetic moments in Fe$_2$P are aligned in the c-direction while the moments in FeMnP$_{0.5}$Si$_{0.5}$ are aligned in the a-direction. The length of the arrows correspond to the magnitude of the magnetic moments.}
\label{magn_str}
\end{figure}

The magnitude of the magnetic moments of the Fe atoms and the total magnetic moment per formula unit is confirmed by recent  M\" ossbauer and magnetization studies of  FeMnP$_{0.5}$Si$_{0.5}$ \cite{mossbauer_uu}. The total magnetic moment has been reported to be approx. 4.4 $\mu_B$/f.u which is in good agreement with our result, where the total magnetic moment is 4.4(2) $\mu_B$/f.u. The size of the total magnetic moment is higher than the corresponding value for Fe$_{2}$P (2.9(1) $\mu_B$) \cite{fujii-87, koumina, scheerlinck, tobola}. Earlier studies have also shown that there are significant differences in magnitude of the moments between the two Fe sites in  Fe$_{2}$P. The Fe atom in the pyramidal 3f site possesses a smaller magnetic moment than the Fe atom in the tetrahedral 3g site, which is illustrated in Figure \ref{magn_str}b.

The refinements of the neutron powder diffraction intensities of FeMnP$_{0.5}$Si$_{0.5}$ indicate that the Mn atoms prefer to be situated in the pyramidal 3f site in the Fe$_{2}$P-structure. The substitution of Fe with Mn in the 3f site increases the magnetic moment in the site with $\sim$0.8 $\mu_B$  compared to Fe$_{2}$P. The magnetic moment of the Fe atom in the 3g site of FeMnP$_{0.5}$Si$_{0.5}$ is also shown to increase $\sim$0.6 $\mu_B$.

The magnetic structures of FeMnP$_{0.5}$Si$_{0.5}$ and Fe$_{2}$P are shown in Figure \ref{magn_str}. The magnetic moments are aligned in the a-direction while the moments of Fe$_{2}$P are aligned in the c-direction \cite{lundgren_fe2p, fujii_c_axis}. A similar alignment of the magnetic moments as in FeMnP$_{0.5}$Si$_{0.5}$ has also been reported to occur in FeMnP$_{0.5}$As$_{0.5}$ \cite{bacmann-94}. It is of interest to note that in FeMnP$_{0.7}$As$_{0.3}$ the moments deviate from the c-axis by 50\textdegree.

Magnetoelastic transitions have been found for the hexagonal system FeMnP$_{1-y}$As$_{y}$ (0.15 $\le$ y $\le$ 0.66) isostructural with FeMnP$_{0.5}$Si$_{0.5}$. A marked increase both in magnetic hyperfine field and in magnetic moment take place on the Fe tetrahedral site (Fe$_t$) as compared to the tetrahedral site in Fe$_{2}$P \cite{bacmann-94, malaman-1996}. Table \ref{table:mag-state} displays the interatomic distances in paramagnetic and ferromagnetic Fe$_{2}$P, FeMnP$_{0.7}$As$_{0.3}$ and the presently studied compound FeMnP$_{0.5}$Si$_{0.5}$.

\begin{table}[t!]
  \centering
\caption{Magnetic state, Fe saturation magnetic moment and interatomic average distances for tetrahedral Fe$_t$ in FeMX with Fe$_2$P structure, M$_p$ = pyramidal Fe or Mn and X = P, Si and As.}
\begin{tabular}{ l c c c c c c c}
\toprule
 & Magn.  &  Fe$_t$ sat.  &  &  Fe$_t$-X &  Fe$_t$- Fe$_t$ &  Fe$_t$- M$_p$  &  \\
Compound & state (T$_C$)   & mom. ($\mu_B$) & V (\AA$^3$) & (\AA) & (\AA) & (\AA)  & Ref. \\
\midrule
Fe$_{2}$P & PM (295 K)  & - & 103.1  & 2.255 & 2.610 & 2.682 & \cite{hf_ugn} \\
 & FM (77 K) & 1.03 & 102.9 & 2.253 & 2.597 & 2.682 & \cite{scheerlinck}\\
FeMnP$_{0.7}$As$_{0.3}$ & PM (250 K)  & - & 110.7  & 2.311 & 2.638 & 2.755 & \cite{bacmann-94} \\
 & FM (100 K) & 1.25 & 110.2 & 2.312 & 2.747 & 2.743 & \cite{bacmann-94}\\
FeMnP$_{0.5}$Si$_{0.5}$ & PM (450 K) & - & 109.8 & 2.312 & 2.685 & 2.726 & This work\\
 & FM (296 K) & 1.65 & 110.6 & 2.317 & 2.766 & 2.708 & This work\\
\bottomrule
\label{table:mag-state}
\end{tabular}
\end{table}

As can be seen from the table the cell volume, the average near Fe$_t$--X distances and Fe$_t$--M$_p$ do not change significantly between the paramagnetic (PM) and the ferromagnetic (FM) state. However the crystal a- and c-axis decreases and increases, respectively, when passing the first order ferromagnetic transition from lower temperature. A large interatomic difference is however observed for the Fe$_t$--Fe$_t$ distances for the  high moment cases of FeMnP$_{0.7}$As$_{0.3}$ and FeMnP$_{0.5}$Si$_{0.5}$. The cell volume expansion between the three different compounds is obvious. The overall crystal expansion and the increase in Fe$_t$--Fe$_t$  distances yields a stronger Fe$_t$  electron localization and results in a larger Fe$_t$  magnetic moment.

\section{Conclusions}
\label{conclusions}

The magnetocaloric compound FeMnP$_{0.5}$Si$_{0.5}$ has been synthesized and studied regarding the crystallographic and magnetic structure. X-ray and neutron powder diffraction experiments show that the sample is single phase and reveals the magnetic structure of FeMnP$_{0.5}$Si$_{0.5}$. An isotructural phase transition has been observed at about the same temperature as the magnetic phase transition (382 K). The Fe atoms are mainly situated in the tetrahedral 3g site while the Mn atoms prefer the pyramidal 3f position. The magnetic moments derived from neutron powder diffraction are shown to be coordinated along the a-axis with a total moment of 4.4 $\mu_B$. This high value of the magnetic moment goes along with our M\"ossbauer results \cite{mossbauer_uu} and is in accord with a strong increase in the Fe$_t$-Fe$_t$ distances. The high magnetic moment and a readily tunable transition temperature make (slightly) off-stoichiometric FeMnP$_{0.5}$Si$_{0.5}$ a promising alloy system for magnetocaloric applications.

\section*{Acknowledgments}
This work was financed by the Swedish Research Council and the Swedish Energy Agency, which is gratefully acknowledged.

\bibliographystyle{model1a-num-names}


Magnetic, electronic, and structural investigation of the strongly correlated Y$_{1-x}$Sm$_x$Co$_5$ system

\section*{Abstract}
{\bf
SmCo$_5$ and YCo$_5$ are isostructural compounds both showing large magnetocrystalline anisotropy, where the former originates mainly from the crystal-electric field and magnetic interactions. We investigate the contribution of both interactions by partially substituting Y by Sm in as-cast polycrystalline Y$_{1-x}$Sm$_{x}$Co$_5$ (with $x$=0, 0.1, 0.2, 0.3, and 0.4) and measuring their structural, magnetic, and electrical properties through X-ray diffraction, magnetization, and electrical transport measurements. Our results suggest an interplay between microstructure strain in as-cast samples and the electronic and magnetic interactions.
}

\section{Introduction}
\label{sec:intro}
SmCo$_5$ alloys are well-known high-performance rare-earth permanent magnets which have been intensively studied in the last five decades due to their high magnetization saturation, high Curie temperature, and high coercivity \cite{strnat}. SmCo$_5$ has a hexagonal CaCu$_5$-type crystal structure (space group P6/mmm) and its excellent magnetic properties result from the large magnetocrystalline anisotropy originated mainly from the Sm localized 4$f$-electrons interaction with the hexagonal crystal field (CF) \cite{larson,anisotropy}. On the other hand, the isostructural compound YCo$_5$ also has a considerably high magnetocrystalline anisotropy but it originates from the spin-orbit coupling of the large orbital angular momentum of the Co 3$d$-electrons \cite{anisotropy,steinbeck}.

Y$_{1-x}$Sm$_{x}$Co$_5$ is a promising candidate system to understand the role of the rare-earth localized $f$-electrons and the contribution of the transition metal itinerant $d$-electrons to the magnetocrystalline anisotropy and magnetic properties. Moreover, the different atomic ratios of Sm and Y allow to investigate the effect of changing the atomic spacing in the packed hexagonal crystal structure on the physical properties \cite{steinbeck}. We have successfully synthesized polycrystalline Y$_{1-x}$Sm$_{x}$Co$_5$ ($x$=0, 0.1, 0.2, 0.3, and 0.4) solid solutions and investigated, for the first time, the influence of the partial substitution of Y by Sm in the structural, magnetic, and electrical properties.

\section{Experimental details}
\label{sec:exp}
Polycrystalline Y$_{1-x}$Sm$_{x}$Co$_5$ samples with $x$=0, 0.1, 0.2, 0.3, and 0.4 were synthesized by arc-melting the stoichiometric amounts of high-purity parent materials under an Argon atmosphere. The X-ray diffraction (XRD) of the samples were performed with the Cu-K$_\alpha$ radiation using a Rigaku Ultima III diffractometer with a monochromator.

The magnetization vs magnetic field (M-H) hysteresis of each sample was measured using a vibrating-sample magnetometer (VSM) at room temperature with a maximum magnetic field of 15 kOe. For the electrical transport measurements, the samples were polished to acquire a regular shape and their resistances were measured using the standard four-probe method at temperatures from 10 K to 290 K. The electrical contacts on the samples were made using silver conducting paint and 50 $\mu$m gold wires.

\section{Results and discussion}
\label{sec:results}

\subsection{Structural properties}
\label{sec:structural}

Figure \ref{fig:rietvelds} shows the XRD pattern of the as-cast Y$_{1-x}$Sm$_{x}$Co$_5$ ($x$=0, 0.1, 0.2, 0.3, and 0.4) samples and their respective Rietveld refinements. We can successfully recognized the diffraction peaks associated with the CaCu$_5$-type hexagonal crystal structure, besides a small Y$_2$O$_3$ phase (see Fig. \ref{fig:rietvelds}). Fig. \ref{fig:lattice} shows the lattice parameters and the relative volume change as function of Sm substitution $x$, obtained by the Rietveld refinement of our XRD patterns.

\begin{figure}[h]
	\centering
	\includegraphics[width=\textwidth]{Figures/rietvelds.png}
	\caption{XRD patterns (circles) and Rietveld refinements (solid lines) of the polycrystalline Y$_{1-x}$Sm$_{x}$Co$_5$ samples (with x=0, 0.1, 0.2, 0.3 and 0.4). The Miller indexs correspond to the main peaks of the YCo$_5$ phase and the black dots indicate the peaks from the Y$_2$O$_3$ phase.}
	\label{fig:rietvelds}
\end{figure}

\begin{figure}[h]
	\centering
	\includegraphics[width=0.32\textwidth]{Figures/a.png}
	\includegraphics[width=0.32\textwidth]{Figures/c2.png}
	\includegraphics[width=0.32\textwidth]{Figures/DeltaV.png}
	\caption{\label{fig:lattice} Lattice parameters a and c, and unit cell volume of the Y$_{1-x}$Sm$_{x}$Co$_5$ samples for different Sm doping $x$ obtained from the Rietveld refinement method. Since this phase has the hexagonal P6/mmm space group symmetry, the lattice parameter a and b are the same. Dashed lines are guide to eyes.}
\end{figure}

We observe that the lattice parameter $c$ (Fig. \ref{fig:lattice}b) monotonically increases with Sm substitution $x$, except for $x$<0.1. On the other hand, the lattice parameter $a$ (Fig. \ref{fig:lattice}a) does not show the linear variation expected from the Vergard's law, remaining almost constant up to $x\approx0.2$, with a subsequent sizable rapid increase around $x$=0.3. This change in the lattice between $x$=0.2 and $x$=0.3 becomes noticeable also in the variation of volume as a function of doping (see Fig. \ref{fig:lattice}c), indicating a structural deformation.

\subsection{Magnetic properties}
\label{sec:magnetic}

Figure \ref{fig:vsm} shows the M-H hysteresis loops of Y$_{1-x}$Sm$_{x}$Co$_5$ for each Sm substitution. The hysteresis loops suggest a ferromagnetic order for all compositions.

\begin{figure}[h]
	\centering
	\includegraphics[width=0.75\textwidth]{Figures/VSM_doping1.png}
	\caption{\label{fig:vsm} Bulk magnetization as a function of the magnetic field of the Y$_{1-x}$Sm$_{x}$Co$_5$ samples with different dopings of Sm measured at room temperature. The inset shows a zoom to better visualize the coercivity of each sample.}
\end{figure}

From the hysteresis loops we obtain the coercivity and the saturation magnetization for each sample. Fig. \ref{fig:magnetic} shows the evolution of these properties with Sm substitution.

\begin{figure}[h]
	\centering
	\includegraphics[width=0.49\textwidth]{Figures/saturation.png}
	\includegraphics[width=0.49\textwidth]{Figures/coercivity.png}
	\caption{\label{fig:magnetic} (a) Saturation magnetization and (b) coercivity as a function of Sm doping in the Y$_{1-x}$Sm$_{x}$Co$_5$ samples. In (a) the dashed line is a linear fit to the data and in (b) the dashed line is a guide to eyes.}
\end{figure}

We note that the saturation magnetization shows a linear decrease with the Sm substitution. On the other hand, we observe that the coercivity first slightly decreases from $x$=0 to $x$=0.1, and then starts to increase monotonically with the increasing of Sm substitution, in a very similar way to what we observed in the $c$ lattice parameter. This similar dependence of the coercivity and the $c$ lattice parameter suggests magneto-elastic coupling effects \cite{rosner,stillwell} and can also be related to an interplay between the change on the hexagonal crystal structure when substituting Y by Sm and the (Y,Sm)Co$_5$ large uniaxial magnetocrystalline anisotropy along the $c$-axis \cite{steinbeck}, revealing as the easy-axis magnetization.

\subsection{Electrical transport measurements}
\label{sec:electrical}

We compare the temperature dependence of the electrical resistance normalized to its room temperature value for all the five samples, as shown in Fig. \ref{fig:resistances}. We observe that for all samples the resistance shows a decrease with temperature very similar to the observed in other metals \cite{metals}.

\begin{figure}[h]
	\centering
	\includegraphics[width=0.75\textwidth]{Figures/Rnorm_doping2.png}
	\caption{\label{fig:resistances} Electrical resistances normalized at 290 K for polycrystalline Y$_{1-x}$Sm$_{x}$Co$_5$.}
\end{figure}

We can separate the contributions to the total resistivity as follows:

\begin{equation}
	\rho(T) = \rho_0 + \rho_{ee}(T) + \rho_{mag}(T) + \rho_{ph}(T)
\end{equation}

where $\rho_0$ is the residual resistivity, $\rho_{ee}$ is the contribution from the electron-electron interaction, $\rho_{mag}$ is the magnetic contribution from the electron-spin wave scattering, and $\rho_{ph}$ is the contribution coming from the electron-phonon scattering processes. The dominant scattering mechanism varies significantly among materials and with the temperature range. In the case of metals at low temperatures, the strongly interacting electron-electron scattering is usually well described by the Fermi liquid model which results in a T$^2$ dependence of the resistivity. On the other hand, for ferromagnets, the magnon modes quadratic dispersion relation $\omega(k)=\Delta+Dk^2 $ gives rise to a magnetic contribution to the resistivity with the following temperature dependence for $k_BT\ll \Delta$ \cite{magnon}:

\begin{equation}
	\rho_{mag}(T) = b\Delta T e^{-\Delta/T}(1+2\frac{T}{\Delta})
	\label{gap}
\end{equation}

where $\Delta$ is the spin-wave gap, $D$ is the spin-wave stiffness, and b is a constant such that $b\propto1/D^2$. Finally, the electron-phonon interactions usually has a $\rho\propto T^5$ contribution at lower temperatures, and a $\rho\propto T$ contribution at higher temperatures \cite{mannari,kowalczky}.

The attempt to include the $T^5$ contribution from the phonons in the fittings at lower temperatures (below 60 K) did not pay off, which indicates that phonon scattering is not significant. Indeed, it has already been suggested that the phonon contribution for RCo$_5$ systems in this temperature range is much smaller compared to the other contributions \cite{mannari} and, therefore, can be neglected. Furthermore, for the YCo$_5$ sample, a quadratic dependence of the resistivity up to 60 K was already observed previously \cite{kowalczky}. Therefore, we perform a fitting of the electrical resistivity below 50 K for each sample including only the Fermi liquid and magnetic contributions using the following equation:

\begin{equation}
	\rho(T)=\rho_0+AT^2+b\Delta T(1+\frac{2T}{\Delta})e^{-\Delta/T}
	\label{eq:quadratic}
\end{equation}

where $A$ is the quadratic parameter from the electron-electron scattering. Fig. \ref{quadratic-fit} shows that the fittings of eq. \ref{eq:quadratic} have a good agreement with the experimental data up to 50 K and Fig. \ref{rho-fit} shows the values of all parameters obtained from the fittings of the resistivity data as a function of the Sm substitution.

\begin{figure}[h]
	\centering
	\includegraphics[width=0.75\textwidth]{Figures/exponential-fittings.png}
	\caption{\label{quadratic-fit}Temperature dependence of the difference between the electrical resistivity and the residual electrical resistivity ($\rho_{0}$) for our polycrystalline samples Y$_{1-x}$Sm$_x$Co$_5$. Solid lines are fittings according to eq. 3. }
\end{figure}

\begin{figure}[h]
	\centering
	\includegraphics[width=0.45\textwidth]{Figures/RRR.png}
	\includegraphics[width=0.45\textwidth]{Figures/A_rho.png}
	\includegraphics[width=0.45\textwidth]{Figures/b.png}
	\includegraphics[width=0.45\textwidth]{Figures/gap.png}
	\caption{\label{rho-fit} Parameters obtained fitting the experimental with equation \ref{eq:quadratic} for each Y$_{1-x}$Sm$_{x}$Co$_5$ sample. (a) Residual-resistance ratio (RRR) obtained from $\rho_0$ using RRR=$\rho$(290 K)/$\rho_0$. (b) Quadratic parameter from the Fermi liquid contribution. (c) Magnetic contribution parameter where $b\propto1/D^2$. (d) Spin-wave gap from the magnetic contribution.}
\end{figure}

We can note in Fig. \ref{rho-fit}a that the residual-resistance ratio (RRR) initially decreases almost linearly until $x$=0.2 followed by an increase at $x$=0.3. Besides that, the $A$ and $b$ parameters (Figs. \ref{rho-fit}b and c) show an overall tendency of decrease except at $x$=0.3 where all of them present an abrupt increase. This abrupt change agrees with the in-plane lattice parameters (see Fig. \ref{fig:lattice}a) and unit cell volume (see Fig. \ref{fig:lattice}c) changes at $x$=0.3, which suggests strain effects on the electronic and magnetic scattering mechanisms \cite{stillwell}. The decrease of the $b$ parameter (Fig. \ref{rho-fit}c) means an increase of the spin-wave stiffness $D$ ($b\propto1/D^2$) and, therefore, a hardening of the magnon modes with increasing Sm substitution. Finally, in Fig. \ref{rho-fit}d we observe that the magnon gap has a maximum around $x$=0.1 and then starts to decrease monotonically, which might indicate a magnetic anisotropy loss.

\section{Conclusion}

Polycrystalline Y$_{1-x}$Sm$_{x}$Co$_5$ as-cast samples were successfully synthesized, where the physical properties were not dramatically affected by the presence of a small spurious yttrium/samarium oxide phase. Our results indicate an abrupt change in the lattice parameter $a$ and in the unit cell volume (Fig. \ref{fig:lattice}) at the same Sm substitution of $x$=0.3, which influences the electric and magnetic scattering process. In addition, similar dependence of lattice parameter $c$ and the coercivity with the Sm substitution suggests the important role that the $f$-electron orbitals plays on the magnetocrystalline anisotropy and crystal electrical field, both underlying mechanisms to understand the magnetic ground state in RCo$_5$ magnets. Finally, the electrical transport measurements indicated the presence of both electron-spin and electron-electron interactions and an apparent hardening of the magnon modes with Sm substitution.

Our results trigger further investigation at higher substitutions of Sm in order to understand the role that strong correlation between itinerant and localized electrons plays on setting the magnetic ground state in hard magnets.

\section*{Acknowledgements}

\paragraph{Funding information}

FP acknowledges FAPESP grant number 2019/24711-0. VM acknowledges the support of JP-FAPESP (2018/19420-3). JLJ acknowledges support of JP-FAPESP (2018/08845-3) and CNPq-PQ2 fellowship (310065/2021-6).

Strongly Correlated Electrons (cond-mat.str-el); Materials Science (cond-mat.mtrl-sci)