Within the family of cubic YbCu$_4$X compounds ($X$ = Ni, Au and Zn), we have investigated the YbCu$_{5-x}$Zn$_x$ ($1\geq x \geq 0.7$) alloys by means of structural, magnetic, thermal and transport measurements. In the $Ï1-$ YbCu$_{5-x}$Zn$_x$ (cubic AuBe$_5$ type, $0.7 \leq x \leq 1.5$) structural phase, Yb ion is in its Yb$^{3+}$ magnetic configuration. However, by increasing Zn content the unit cell grows faster than a reference computed as a Cu by Zn atoms substitution, which indicates a shift of Yb ions towards the larger Yb$^{2+}$ configuration. The magnetic behavior confirms such tendency with a clear decrease of the saturation magnetization and effective moment between $x=0.7$ and $x=1$. The specific heat at low temperature shows a logarithmic dependence characteristic for a non-fermi-liquid behavior. The characteristic energies of all studied parameters, including magneto resistivity, show notably low values as an indication that these alloys are close to a quantum critical point, which is approached from the non-magnetic side as the Zn content decreases.
- 論文ID: 2511.00277
- タイトル: I. Magnetic, thermal and transport properties of YbCu5−xZnx alloys
- 著者: I. Čurlik, F. Akbar, S. Gabani, M. Giovannini, J.G. Sereni
- 分類: cond-mat.str-el(強相関電子系)
- 発表日: 2025年11月4日
- 論文リンク: https://arxiv.org/abs/2511.00277
本研究は立方晶YbCu4X化合物族(X = Ni, Au, Zn)に属するYbCu5−xZnx(1 ≥ x ≥ 0.7)合金の構造的、磁気的、熱的および輸送特性を系統的に測定した。τ1-YbCu5−xZnx構造相(立方AuBe5型、0.7 ≤ x ≤ 1.5)において、Yb イオンはYb3+磁気配置に位置する。しかし、Zn含量の増加に伴い、格子定数の増加速度がCuをZn原子で置換する理論的予測を上回り、Ybイオンがより大きなYb2+配置への転移を示唆している。磁気的挙動はこの傾向を確認し、飽和磁化と有効磁気モーメントはx = 0.7からx = 1の間で著しく減少する。低温比熱は非フェルミ液体挙動の特徴である対数依存性を示す。すべての研究パラメータ(磁気抵抗を含む)の特性エネルギーは著しく低い値を示し、これらの合金が量子臨界点に接近し、非磁性側からこの臨界点に接近していることを示唆している。
本研究はYbCu5−xZnx合金系における電子濃度変化が量子臨界挙動に及ぼす影響を探索することを目的とし、特に以下の点を対象とする:
- 化学ポテンシャル変化と電子濃度変化が量子臨界点(QCP)に及ぼす影響の区別
- 価数不安定性と価数変動がこの系における役割の理解
- 非フェルミ液体挙動の微視的機構の探索
- 量子臨界物理:量子臨界点近傍の物理は凝聚系物理の中心課題であり、非従来型超伝導、重フェルミ子挙動などの重要な現象と密接に関連している
- 幾何学的フラストレーション:この材料の立方構造は共有辺四面体ネットワークとして見なすことができ、Yb磁気イオンが頂点に位置し、三角格子の3D類似体を構成し、磁気フラストレーション効果を研究するための理想的なプラットフォームを提供する
- 価数制御:Yb基化合物の価数不安定性は磁気と電子相関を制御するための独特の手段を提供する
- YbCu5−xAux系:深く研究されているが、Au置換はCuの化学ポテンシャルを主に変化させ、電子構造をあまり変えないため、異なる物理機構の寄与を分離することが困難である
- 系統的研究の欠如:YbCu5−xZnx系は合成されているが、包括的な物理特性の表現が欠けている
- 理論的理解の不足:非フェルミ液体挙動の微視的起源と量子臨界点の性質に関しては、さらなる実験データが必要である
YbCu5−xZnx系を研究することにより、電子濃度を独立に制御でき(CuとZnは異なる価電子数を持つ)、YbCu5−xAux系と対比させることで、量子臨界挙動を駆動する微視的機構をより深く理解することができる。
- 系統的実験表現:YbCu5−xZnx(0.7 ≤ x ≤ 1)合金の構造、磁気、比熱および電気輸送の初めての包括的測定
- 価数転移の実験的証拠:格子定数分析を通じて、格子体積がVegard則の予測を超過し、磁気測定と組み合わせてYbイオンのYb3+からYb2+への価数転移を明らかにした
- 非フェルミ液体挙動の確認:低温比熱がC4f/T ∝ -log(T)の特性依存性を示し、特性エネルギーTq = 7.5 Kで非フェルミ液体基底状態を証実した
- 結晶場エネルギー準位構造の決定:比熱と磁化率の結合分析を通じて、結晶場分裂を正確に決定した:第一励起状態Δ1 = 30 ± 2 K、第二励起状態Δ2 = 46 ± 3 K
- 量子臨界点の位置決定:複数の独立した物理量(磁化率、比熱、磁気抵抗)の特性エネルギーがすべて7~8.7 K範囲内にあり、系が量子臨界点に接近し、非磁性側から接近していることを示す
- 濃度依存挙動の明示:Zn濃度が磁気モーメント、特性温度および非フェルミ液体挙動の強度に及ぼす影響を系統的に研究し、完全な相図を確立した
本研究は複数の実験技術を通じてYbCu5−xZnx合金の物理特性を表現する:
- 入力:異なるZn濃度(x = 0.7, 0.8, 0.9, 1.0)の多結晶試料
- 出力:結晶構造パラメータ、磁化、磁化率、比熱、電気抵抗率の温度および磁場依存性
- 目標:基底状態特性、結晶場エネルギー準位、特性エネルギースケール、量子臨界点の位置決定
- 結晶構造:立方AuBe5型構造、空間群F4ˉ3m
- 格子体積分析:理論的Vegard則との比較、参照式:
VVegard(x)=V0+αx
ここでCuとZnの金属半径の差異を考慮
- 価数分析:格子超過膨張ΔV = Vexp - VVegardはYb価数変化を示す
磁化M(B,T):
- 磁場範囲:0-9 T
- 温度範囲:2-50 K
- データ拟合:Brillouin関数B1/2(y)を使用、ここでy = geffμBJB/kBT
- 飽和磁化抽出:
M(B,T=2K)=Msattanh(B/B0)
磁化率χ(T):
- 温度範囲:2-300 K
- 外場:B = 1 T
- 高温領域拟合(T > 200 K):Curie-Weiss則
χ(T)1=CCT+θP
- 低温領域改良モデル:
χ(T,x)1=CCT+θGSe−Δ/(T+θP)
結晶場分裂効果を含む
高温比熱CP(T):
- 温度範囲:2-300 K
- Debye温度:TD = 265 K
- 4f電子寄与抽出:
C4f/T=CP/T−CD/T
- エントロピー計算:
S4f(T)=∫0TT′C4f(T′)dT′
低温非フェルミ液体挙動:
複合拟合関数(式1)を採用:
TC4f=a+blog(TTq)+∑i=12cie−(T−Δi)2/(2TKi2)
各項の物理的意味:
- a = 0.15 J/molK2:電子帯寄与
- blog(Tq/T):非フェルミ液体項、Tq = 7.5 Kは量子変動エネルギー
- 指数項:結晶場励起状態のKondo展開Lorentzian関数
電気抵抗率ρ(T,B):
- 温度範囲:2-300 K
- 磁場:0-9 T
- 正規化表示:R/R300K
磁気抵抗スケーリング分析:
Schlottmann スケーリング理論を採用:
r(B,T)=R(0)ΔR(B)=−Aζe−B0/ζ
ここでスケーリング変数:
ζ=T+T∗B
異なる温度曲線の collapse を通じて特性温度T∗を決定
- 改良された磁化率分析モデル:従来のCurie-Weiss則は結晶場効果を正確に記述できず、本論文は結晶場分裂とKondo相互作用を含む指数修正項を導入し、基底状態パラメータをより正確に抽出する
- 比熱データの多スケール拟合:Kondo展開のGaussian関数を使用して励起状態寄与を記述する革新的方法を採用し、標準Schottky関数ではなく、Kondo相互作用がエネルギー準位に及ぼす展開効果をより良く反映する
- 複数物理量の交差検証:磁化率、比熱、磁気抵抗の3つの独立測定を通じて抽出された特性エネルギーを相互検証(7-8.7 K)し、結論の信頼性を強化する
- 価数分析の定量的方法:「ad hoc」Vegard則参照を確立し、格子超過膨張を定量的に評価し、価数転移の直接的な構造的証拠を提供する
- 合金組成:YbCu5−xZnx、x = 0.7, 0.8, 0.9, 1.0
- 構造相:τ1相、立方AuBe5型構造
- 試料形態:多結晶試料
- 製備方法:参考文献9を参照
- 磁気測定:SQUID磁力計またはVSM
- 比熱測定:
- T < 3 K:PPMS(Quantum Design)He3挿入ユニット
- T > 2 K:He4ポンプシステム
- データマッチング区間:2-3 K
- 電気輸送:標準四探針法
- 温度範囲:0.4 K~300 K(異なる測定技術により異なる)
- 磁場範囲:0~9 T
- 磁化率測定場:B = 1 T
- 背景除去:
- 比熱:Debye音声子寄与を除去 CD(TD = 265 K)
- 磁化率:高温線形外挿によりθPを決定
- 拟合方法:
- 誤差分析:拟合パラメータが不確定性を与える(例:Δ1 = 30 ± 2 K)
格子体積進化(図1):
- YbCu4Zn(x=1):V ≈ 87.5 Å3
- YbCu4.3Zn0.7(x=0.7):V ≈ 85 Å3
- 重要な発見:実験値は系統的にVegard則予測(破線)より高く、偏差はxの増加に伴い増大
- 物理的解釈:格子超過膨張はYb3+(1.00 Å)→ Yb2+(1.16 Å)価数転移を示す
飽和磁化:
- x = 0.7:Msat = 1.23 μB/f.u.(拟合関数:1.23×tanh(B/10.2))
- x = 1.0:Msat = 0.75 μB/f.u.(拟合関数:0.75×tanh(B/9.3))
- 濃度依存性:MsatはZn含量増加に伴い著しく低下(低下幅約40%)
基底状態決定:
- 理論値:Γ6基底状態Msat(Γ6)= 1.33 μB、Γ7基底状態Msat(Γ7)= 1.72 μB
- 実験値1.23 μBはΓ6により近く、基底状態がΓ6二重項であることを確認
高温挙動(T > 200 K):
- 常磁性Curie温度:θP = 26 ± 2 K(すべての濃度で一致)
- Kondo温度推定:TK = θP/√2 ≈ 15 K
- 注記:1/χ(T)は中温領域で直線から逸脱し、結晶場励起状態の退布居を反映
低温基底状態パラメータ(図3挿図):
- 有効磁気モーメント:
- x = 0.7:μeff = 4.57 μB
- x = 1.0:μeff = 3.95 μB
- 理論値(Yb3+、J=7/2):μeff = 4.54 μB
- 基底状態常磁性温度:
- x = 0.7:θGS = 4.7 K
- x = 1.0:θGS = 7 K
- 結晶場分裂:Δ = 38-42 K
自己無撞着性検証:
μeff(基底状態)/Msat(2K)= 2.13/1.23 ≈ √3、S=1/2二重項理論予測と一致
高温比熱(図4a):
- Debye温度:TD = 265 K(CD/T最大値位置0.28TDにより決定)
- LuCu5−xAgxと一致し、音声子スペクトルの合理性を検証
- 4fエントロピー(挿図):S4f → Rln2(基底状態二重項)+ Rln2(第一励起状態)+ Rln4(Γ8四重項)
低温非フェルミ液体挙動(図4b):
- YbCu4Zn主曲線:C4f/Tは0.4-10 Kで-log(T)依存性を示す
- 拟合パラメータ(式1):
- 電子項:a = 0.15 J/molK2
- NFL係数:b = 0.26 J/molK2
- 量子変動エネルギー:Tq = 7.5 K
- 第一励起状態:Δ1 = 30 ± 2 K、展開TK1 = 15 ± 1 K
- 第二励起状態:Δ2 = 46 ± 3 K、展開TK2 = 40 ± 2 K
濃度依存性(図4b挿図):
- TQ = 9 Kはすべての濃度で一定に保たれる
- CNFL/T強度はZn含量低下に伴い増強し、QCPへの接近を示す
零磁場電気抵抗率(図5b):
- 半対数スケールで単一不純物散乱の特性を示す
- 低温で飽和傾向、コヒーレンスピークなし(フェルミ液体と異なる)
- 濃度依存性:定性的挙動は類似、定量値に差異
磁場依存性(図5a、YbCu4Zn):
- 磁場は散乱を抑制し、電気抵抗を低下させる
- 低温曲率は磁場に伴い系統的に変化
負磁気抵抗(図6a):
- r(B) = ΔR(B)/R(0) < 0
- T = 2 Kで、9 T場下でr ≈ -18%
- 高温効果は減弱
Schlottmannスケーリング(図6b):
- スケーリング変数:ζ = B/(T + T∗)
- 特性温度:T∗ = 7 K
- スケーリング関数:r(B,T) = -29ζ exp(-0.3/ζ)
- 意義:2, 5, 10, 20 Kデータが単一曲線に collapse し、単一エネルギースケール挙動を証実
- 複数エネルギースケールの一致性:
- 磁化率:θGS = 4.7-7 K
- 比熱:Tq = 7.5-9 K
- 磁気抵抗:T∗ = 7 K
- エントロピー半高:8.7 K
- 結論:特性エネルギーは7-9 Kに集中し、典型的な重フェルミ子系(~100 K)より遠く低い
- 非フェルミ液体基底状態:
- C/T ∝ -log(T)は広い温度領域(0.4-10 K)で成立
- フェルミ液体のC/T = constと区別
- 価数転移の証拠チェーン:
- 構造:格子超Vegard膨張
- 磁気:μeffとMsat低下
- 熱力学:エントロピー減少
- 傾向:Zn含量↑ → Yb2+成分↑ → 磁気↓
- 量子臨界点の位置決定:
- 系はQCPの非磁性側に位置
- Zn含量低下はQCPに接近
- 幾何学的フラストレーションが長距離秩序を抑制
- YbCu5−xAux5-8:
- 最も深く研究された体系
- 40 mKまで磁気秩序なし
- QCPに接近するが機構が異なる(化学ポテンシャル制御 vs 電子濃度制御)
- YbCu4Ni3,19:
- T = 0.4 K磁気秩序
- 極めて高い電子比熱係数
- 重フェルミ子挙動
- Parks理論11:価数不安定性の初期フレームワーク
- Wohlleben経験則12:価数変動の判定基準
- EuPd2(Si2−xGex)13:価数変動→長距離秩序の実例
- Desgranges-Rasul モデル17:簡併Kondo模型における結晶場
- YbCu5−xAgx18:結晶場とKondoエネルギースケールの競争
- 本論文の革新:Kondo展開のGaussian vs 標準Schottky関数
- Schlottmannスケーリング21:磁気抵抗の普遍的スケーリング挙動
- 電核結合20:QCPでのYbCu4Znの超微細相互作用
- 本論文の貢献:広い温度領域での対数比熱の実験的確認
- 系統性:YbCu5−xZnxの物理特性の初めての包括的表現
- 補完性:YbCu5−xAuxとの対比により、化学ポテンシャルと電子濃度効果を分離
- 精度:複数技術の交差検証により、パラメータ抽出がより信頼性が高い
- 完全性:構造から輸送までの全方位測定
- 基底状態特性の確認:
- YbCu5−xZnx基底状態はΓ6二重項
- 低温は非フェルミ液体挙動を示す(C/T ∝ -log T)
- 特性エネルギースケールは異常に低い(7-9 K)
- 価数進化の法則:
- Zn含量増加がYb3+ → Yb2+転移を駆動
- 磁気性を減弱させる:μeffは4.57から3.95 μBに低下
- Msatは1.23から0.75 μBに低下
- 量子臨界点の特性:
- 系はQCPの非磁性側に位置
- Zn含量低下がQCPに接近
- 幾何学的フラストレーションが長距離磁気秩序を阻止
- 結晶場構造:
- 第一励起状態:30 K(二重項)
- 第二励起状態:46 K(四重項)
- Kondo相互作用は著しいエネルギー準位展開を引き起こす(15-40 K)
- 試料限定:
- 0.7 ≤ x ≤ 1範囲のみ研究
- より低いZn濃度(磁気側)データの欠如
- 多結晶試料、異方性研究不可
- 温度ウィンドウ:
- 最低温度0.4 K、真の量子臨界領域(T → 0)に達していない
- 低温磁気秩序の観測不可
- 希釈冷却機温度(< 50 mK)が必要
- 理論分析の深さ:
- 結晶場拟合は唯象モデル(Gaussian展開)を採用
- 第一原理計算の支持がない
- 非フェルミ液体起源が完全に明確化されていない
- 欠落した測定:
- 中性子散乱なし(磁気構造と変動確定)
- 光電子分光なし(価数直接測定)
- 圧力依存なし(QCP位置制御)
- 濃度範囲の拡張:
- x < 0.7:磁気秩序領域探索
- x > 1:深い非磁性領域
- 完全な相図作成
- 極低温測定:
- T < 50 mK:隠れた相転移探索
- 交流磁化率:動力学情報
- μSR:局所磁気探測
- 単結晶研究:
- 異方性磁気
- Fermi面測定(量子振動)
- 中性子散乱(磁気励起スペクトル)
- 理論協力:
- 動力学平均場(DMFT)計算
- 量子Monte Carlo シミュレーション
- 繰り込み群分析
- 外場制御:
- 高圧:4f-導帯ハイブリッド化調整
- 強磁場(> 30 T):磁場誘起相転移
- 化学ドーピング:無秩序効果導入
- 複数技術の組み合わせ:構造、磁気、熱学、輸送の四位一体、相互検証
- 広い温度範囲:0.4 Kから300 Kまで、すべての特性エネルギースケールをカバー
- 高いデータ品質:低温比熱がHe3とHe4システム間で滑らかに遷移、磁化曲線は滑らか
- 改良された磁化率モデル:指数修正項e−Δ/(T+θP)は結晶場効果を正確に捉える
- 複合比熱拟合:NFL項+Kondo展開励起状態、物理像が明確
- Schlottmannスケーリング検証:異なる温度データの collapse により単一エネルギースケール仮説を検証
- 価数転移の証拠チェーン:構造(格子膨張)→磁気(磁気モーメント低下)→熱力学(エントロピー減少)、論理が自己無撞着
- エネルギースケール一致性:4つの独立測定が7-9 Kを与え、結論の信頼性を強化
- 濃度依存の系統性:x = 0.7, 0.8, 0.9, 1.0の4点が明確な傾向を描く
- 図表情報が豊富(6つの主図+複数の挿図)
- 拟合パラメータと誤差が完全に与えられている
- 文献との比較が十分(21の参考文献)
- 唯象モデル:Kondo展開にGaussian関数を使用することは経験的であり、微視的導出が欠ける
- NFL機構:C/T ∝ -log Tのみ確認されるが、起源は議論されていない(Kondo破裂?価数変動?)
- QCP性質:スピン密度波型、価数転移型、その他のいずれかが明確でない
- 濃度間隔が大きい:Δx = 0.1は微細構造を見落とす可能性
- x < 0.7のデータ欠如:QCP正確位置の決定不可
- 単一磁場方向:多結晶試料は異方性研究を制限
- Vegard則参照:「ad hoc」であることを認め、Lu基参照化合物が欠ける
- 電気抵抗率分析が浅い:図5bは定性比較のみ、特性温度抽出なし
- 磁気フラストレーション効果:序論で言及のみ、データに反映されていない
- 中核動機:化学ポテンシャル vs 電子濃度の区別であるが、系統的対比図が欠ける
- 定量比較:両体系の特性エネルギー、NFL係数などのパラメータ並列がない
- 中程度以上:強相関電子系分野において、新しい材料体系の基準データを提供
- 方法論的価値:改良された分析モデルは他のYb基化合物に適用可能
- 示唆性:価数制御+幾何学的フラストレーションの組み合わせはQCP研究に新しい視点を提供
- 材料科学:YbCu5−xZnxの相図と物性データベース
- 応用前景は限定的:基礎研究志向であり、直接応用なし
- 技術示範:複数技術の組み合わせの範例
- 試料製備:参考文献9に依存し、本論文では詳述されていない
- 測定条件:パラメータが完全に与えられている(温度、磁場、拟合関数)
- データ可用性:論文ではデータ共有に言及していない、著者への連絡が必要な可能性
- 量子臨界物理研究:
- QCPに接近するが有序状態に入らないモデル系として
- 磁気フラストレーションが量子相転移に及ぼす影響を研究
- 価数物理:
- 4f電子の局所-巡回転移を研究
- 価数変動理論をテスト
- 重フェルミ子材料:
- 特性エネルギーは低いが基底状態がNFL挙動を示す
- 典型的な重フェルミ子系との対比
- 方法論的参照:
- 複数技術の組み合わせ表現の標準プロセス
- 結晶場+Kondo相互作用の分析方法
- 不適用シーン:
- 単結晶異方性を必要とする研究
- 極低温(< 0.4 K)の量子臨界挙動
- 微視的磁気構造情報(中性子散乱必要)
- Sarrao et al., PRB 59, 6855 (1999):YbCu4X族の基礎的研究
- Banda et al., arXiv:2308.15294 (2023):YbCu5−xAuxの電核量子臨界性
- Čurlik et al., PRB 90 (2014):YbCu5−xAuxの磁気励起高密度状態
- Akbar et al., J. Alloys Compd. 976, 173195 (2024):YbCu5−xZnxの相平衡と結晶構造
- Schlottmann, Z. Phys. B 51, 233 (1983):磁気抵抗スケーリング理論
総合評価:これは実験技術が堅牢で、データ品質が高く、物理像が明確な強相関電子系研究論文である。YbCu5−xZnx合金の系統的表現を通じて、価数転移、非フェルミ液体挙動および量子臨界特性を成功裏に明らかにし、Yb基化合物の複雑な物理を理解するための重要な実験基礎を提供している。主な不足は理論分析の深さと特定の測定の完全性にあるが、瑕疵は全体を損なわず、この分野の重要な貢献である。後続の単結晶研究と極低温測定により理解を深めることを推奨する。