The concept of spin symmetries has gained renewed interest as a valuable tool for classifying unconventional magnetic phases, including altermagnets and recently identified p-wave magnets. In this work, we show that in compounds with weak spin-orbit coupling, the dominant spin and charge photoresponse is determined by spin group rather than the conventional magnetic group symmetry. As a concrete realization we consider the nonlinear shift photocurrent in Mn$_5$Si$_3$, a material that features the two possible classes of unconventional p-wave magnetism in the form of two competing spin structures, a coplanar and non-coplanar one. While both are predicted to generate shift currents based on magnetic symmetry considerations, only the non-coplanar configuration survives the spin symmetry requirements. This is numerically confirmed by our $\textit{ab-initio}$ calculations, providing a protocol to experimentally identify the spin configuration of this promising material in photogalvanic or transport measurements.
論文ID : 2406.19842タイトル : Optical signatures of spin symmetries in unconventional magnets著者 : Javier Sivianes, Flaviano José dos Santos, Julen Ibañez-Azpiroz分類 : cond-mat.mes-hall(凝縮系物理-メゾスコピック・ホール効果)発表日時 : 2025年3月10日(v3版)論文リンク : https://arxiv.org/abs/2406.19842 本論文は、スピン対称性が非従来型磁性体の光学応答に果たす役割を研究している。弱いスピン軌道相互作用を持つ材料について、著者らは従来の磁気群対称性ではなく、スピン群対称性が主導的なスピンおよび電荷光応答を決定することを実証した。Mn₅Si₃を例として、この材料は低温で2つの可能なp波磁気構成(共面および非共面)を持つ。磁気対称性解析は両構成がシフト光電流を生成できることを予測するが、非共面構成のみがスピン対称性要件を満たす。第一原理計算がこの予測を確認し、光電圧または輸送測定実験を通じて材料のスピン構成を決定するための実行可能なスキームを提供する。
従来の磁性理論は強磁性と反強磁性が支配的なパラダイムであったが、近年出現した交替磁性 (altermagnetism)およびp波磁性 (p-wave magnetism)などの非従来型磁性相がこのフレームワークに異議を唱えている。これらの新型磁性相を効果的に分類・表現する方法が急務である。
新しい物理メカニズム :非従来型磁性体は独特の輸送特性(異常ホール効果、トポロジカルホール効果など)を示すが、その微視的メカニズムはまだ不明確である実験的鑑別の困難性 :複数の可能な磁気構成を持つ材料(Mn₅Si₃など)について、従来の中性子散乱などの方法では明確な区別が困難であるデバイス応用の可能性 :これらの材料はスピントロニクスおよび光電圧デバイスに重要な応用前景を持つ磁気群対称性解析 :従来の方法は磁気点群に基づくが、弱いスピン軌道相互作用系ではスピンと格子自由度がほぼ解耦し、磁気群解析は過度に緩い対称性制約を与える可能性がある実験的表現手段 :Mn₅Si₃の低温(T<60K)非共線磁気構造には議論があり、既存の実験技術では正確な構成を決定することが困難である著者らは、スピン群対称性 (spin group symmetry)解析を非線形光学応答、特にシフト光電流(shift photocurrent)に適用することを提案し、磁気構成を探査する新しい手段として利用する。この方法は、スピン対称性が光応答テンソルに追加の制約を与えることを利用し、磁気群対称性では区別できない磁気構造を区別することができる。
理論フレームワークの確立 :スピン群対称性を光電導率テンソル解析に初めて体系的に適用し、弱いスピン軌道相互作用系では有効点群がスピン対称性の空間部分によって決定されることを証明したモデル検証 :二重交換モデルの数値計算を通じて、スピン対称性で禁止された光電導成分が許可された成分より3~4桁小さいことを確認した材料予測 :Mn₅Si₃の2つの候補磁気構造に対する第一原理計算を実施し、以下を発見した:非共面構造のシフト光電流は~10⁻⁴ A/V²に達する(Weyl半金属と同等) 共面構造の応答は2桁小さい 測定可能な角度依存光電流特性(~20 nA)を提供 実験スキーム :偏光角依存の直流光電流測定を通じてMn₅Si₃の真の磁気構成を鑑別する実験プロトコルを提案した普遍性への拡張 :スピン光電流への分析を拡張し、スピン対称性制約が同様に適用されることを証明した入力 :磁性材料の結晶構造と磁気構成出力 :2次非線形光電導率テンソルσᵃᵇᶜ(ω)およびその対称性制約制約 :弱いスピン軌道相互作用近似(スピンと格子自由度の解耦)
スピン群操作は Rₛ||Rᵣ|τ と表現される。ここで:
Rₛ:スピン空間に作用する操作 Rᵣ|τ:従来の空間群操作(回転+平行移動) 重要な洞察:電荷応答は純粋なスピン操作Rₛの下で不変であるため、光電導テンソルは空間操作Rᵣから構成される有効点群によってのみ制約される 。
長さゲージでのシフト光電導率:
σ a b c ( ω ) = − i π e 3 2 ℏ 2 ∫ d k ( 2 π ) 3 ∑ m n f n m ( I m n a b c + I m n a c b ) [ δ ( ω n m − ω ) + δ ( ω m n − ω ) ] \sigma^{abc}(\omega) = -\frac{i\pi e^3}{2\hbar^2}\int\frac{dk}{(2\pi)^3}\sum_{mn}f_{nm}(I^{abc}_{mn} + I^{acb}_{mn})[\delta(\omega_{nm}-\omega) + \delta(\omega_{mn}-\omega)] σ ab c ( ω ) = − 2 ℏ 2 iπ e 3 ∫ ( 2 π ) 3 d k ∑ mn f nm ( I mn ab c + I mn a c b ) [ δ ( ω nm − ω ) + δ ( ω mn − ω )]
ここで遷移行列要素:
I m n a b c = r m n b r n m c ; a I^{abc}_{mn} = r^b_{mn}r^{c;a}_{nm} I mn ab c = r mn b r nm c ; a
は双極子項 r m n b r^b_{mn} r mn b と一般化導関数 r n m c ; a r^{c;a}_{nm} r nm c ; a を含む。
Wannier関数構築における対称性破れの問題を回避するため、著者らは速度ゲージ を採用した:
σ a b c ( ω ) = π ∣ e ∣ 3 2 ℏ 2 V ∑ k , m , n f m , n ω m , n 2 Im [ v m , n b ∑ l ≠ m , n ( v n , l a v l , m c ω n , l − v n , l c v l , m a ω l m ) + ( b ↔ c ) ] δ ( ω m , n − ω ) \sigma^{abc}(\omega) = \frac{\pi|e|^3}{2\hbar^2V}\sum_{k,m,n}\frac{f_{m,n}}{\omega^2_{m,n}}\text{Im}\left[v^b_{m,n}\sum_{l\neq m,n}\left(\frac{v^a_{n,l}v^c_{l,m}}{\omega_{n,l}} - \frac{v^c_{n,l}v^a_{l,m}}{\omega_{lm}}\right) + (b\leftrightarrow c)\right]\delta(\omega_{m,n}-\omega) σ ab c ( ω ) = 2 ℏ 2 V π ∣ e ∣ 3 ∑ k , m , n ω m , n 2 f m , n Im [ v m , n b ∑ l = m , n ( ω n , l v n , l a v l , m c − ω l m v n , l c v l , m a ) + ( b ↔ c ) ] δ ( ω m , n − ω )
スピン対称性 Rₛ||Rᵣ に対して、電場と電流は以下のように変換される:
E → g R r E , j → g R r j \mathbf{E} \xrightarrow{g} R_r\mathbf{E}, \quad \mathbf{j} \xrightarrow{g} R_r\mathbf{j} E g R r E , j g R r j
ノイマン原理(物理的性質は結晶対称性を持たなければならない)を利用して:
σ a ′ b ′ c ′ = R r a ′ a R r b ′ b R r c ′ c σ a b c = σ a b c \sigma^{a'b'c'} = R^{a'a}_r R^{b'b}_r R^{c'c}_r \sigma^{abc} = \sigma^{abc} σ a ′ b ′ c ′ = R r a ′ a R r b ′ b R r c ′ c σ ab c = σ ab c
これは従来の点群制約と同じ形式であるが、R_rは磁気点群ではなくスピン群から来ている 。
スピンシフト電流テンソル:
ξ shift s , a b c = i π q 3 2 ℏ 2 V Ω ∑ k , m , n f m n I m n s , a b c ( k ) δ ( Ω + ω m n ) \xi^{s,abc}_{\text{shift}} = \frac{i\pi q^3}{2\hbar^2V\Omega}\sum_{k,m,n}f_{mn}I^{s,abc}_{mn}(k)\delta(\Omega+\omega_{mn}) ξ shift s , ab c = 2 ℏ 2 V Ω iπ q 3 ∑ k , m , n f mn I mn s , ab c ( k ) δ ( Ω + ω mn )
変換規則は追加のスピン回転を含む:
ξ s ′ , a ′ b ′ c ′ = Det ( R s ) R s s ′ s R r a ′ a R r b ′ b R r c ′ c ξ s , a b c \xi^{s',a'b'c'} = \text{Det}(R_s)R^{s's}_s R^{a'a}_r R^{b'b}_r R^{c'c}_r \xi^{s,abc} ξ s ′ , a ′ b ′ c ′ = Det ( R s ) R s s ′ s R r a ′ a R r b ′ b R r c ′ c ξ s , ab c
ハミルトニアン:
H = ∑ ⟨ i , j ⟩ , α t 1 ( c i α † c j α + h . c . ) + ∑ ⟨ ⟨ i , j ⟩ ⟩ , α t 2 ( c i α † c j α + h . c . ) + J ∑ i , α , α ′ m i ⋅ c i α † S α α ′ c i α ′ H = \sum_{\langle i,j\rangle,\alpha}t_1(c^\dagger_{i\alpha}c_{j\alpha} + h.c.) + \sum_{\langle\langle i,j\rangle\rangle,\alpha}t_2(c^\dagger_{i\alpha}c_{j\alpha} + h.c.) + J\sum_{i,\alpha,\alpha'}\mathbf{m}_i\cdot c^\dagger_{i\alpha}\mathbf{S}_{\alpha\alpha'}c_{i\alpha'} H = ∑ ⟨ i , j ⟩ , α t 1 ( c i α † c j α + h . c . ) + ∑ ⟨⟨ i , j ⟩⟩ , α t 2 ( c i α † c j α + h . c . ) + J ∑ i , α , α ′ m i ⋅ c i α † S α α ′ c i α ′
パラメータ設定:
格子定数 a=2Å 遷移積分 t₁=1 eV, t₂=0.25 eV 交換定数 J=0.4 eV/μ²ᵦ 方法 :VASPソフトウェアパッケージ、PAW擬ポテンシャル、制約局所磁気モーメント法パラメータ :エネルギーカットオフ440 eV、k網格9×5×11、収束基準10⁻⁷ eV磁気モーメント :Mn原子最大磁気モーメント~2.6 μᵦエネルギー差 :共面と非共面構造のエネルギー差5 meV/atom2次元正方格子 (図1a):
磁気空間群:P1.1' (T τₘ対称性のみを含む) スピン対称性:T||E|τₘ およびE||Mₓ|0 有効点群:m (鏡面対称) 対称性許可成分:σʸʸʸ, σʸˣˣ, σˣˣʸ 磁気空間群:PS₁ = {E, Tτₓᵧ} スピン対称性有効点群:m 許可成分:z指標が奇数の成分は禁止される 磁気空間群:PC2221 = {E, C₂ʸ, C₂ˣτz, C₂ʸτz} ⊗ {E, Tτₓᵧ} スピン対称性:C₂ˣT||E|0 (磁気モーメントはyz平面内に共面) 有効点群:mmm(空間反転対称性を回復) すべての非線形光応答は非相対論的極限で禁止される 光子エネルギー:ℏω = 0.15 eV 光強度:480 nW スポットサイズ:w = 5 μm 偏光:E(ω) = E₀(ω)(0, cos(θ), sin(θ)) 許可成分σʸˣˣ :ピーク値~10⁻⁴ A/V² (|m|=0.75μᵦ時)禁止成分σˣˣˣ :許可成分より3桁小さい 数値結果はスピン対称性予測と完全に一致 重要な発見:
シフト光電流は|m|に単調に増加 |m| < 1 μᵦ:応答は弱い |m| ≈ 1.5 μᵦ:~10⁻⁴ A/V²に達する(Weyl半金属TaAsと同等) |m|→0極限では光電流は消失する (従来の強誘電体と異なる)ピーク位置の進化:
ω₁ ≈ 1.4 eV主導 (|m| < 0.33μᵦ) ω₂ ≈ 1.9 eV主導 (|m| > 0.38μᵦ) スピン対称性で予測された許可成分(ξˣ'ˣˣʸ, ξʸ'ˣˣʸ, ξᶻ'ˣˣʸ)と禁止成分の差は4桁に達する 。
非共面構造 :
金属性、フェルミ面は全方向スピン分裂を示す Sˣ, Sʸ, Sᶻ成分はすべて有限 スピン分裂~0.5 eV(非相対論的起源) 共面構造 :
金属性、フェルミ面スピン分裂 Sʸ, Sᶻはほぼゼロ(Γ点付近を除く) k空間スピン偏極は実空間磁化平面に垂直 理論的説明:共面構造のC₂ˣT||E|0 対称性は以下を要求する:
S n y , z ( k ) = ⟨ ψ k n ∣ C 2 x S y , z C 2 x − 1 ∣ ψ k n ⟩ = − S n y , z ( k ) S^{y,z}_n(k) = \langle\psi_{kn}|C_{2x}S^{y,z}C^{-1}_{2x}|\psi_{kn}\rangle = -S^{y,z}_n(k) S n y , z ( k ) = ⟨ ψ kn ∣ C 2 x S y , z C 2 x − 1 ∣ ψ kn ⟩ = − S n y , z ( k )
したがって平行成分はゼロでなければならない。
非共面構造 :
ピーク値:10⁻⁴ A/V² (ℏω0.5 eV) 許可成分(σˣʸʸ, σˣᶻᶻ):10⁻⁴ A/V²のオーダー 禁止成分(σˣʸᶻ):~10⁻⁶ A/V² スピン対称性予測と完全に一致 共面構造 :
すべての成分:~10⁻⁶ A/V² 非共面構造より2桁小さい スピン対称性予測の完全な抑制とほぼ一致 比較:
非共面Mn₅Si₃ > Weyl半金属TaAs/TaIrTe₄ 非共面Mn₅Si₃ >> 半導体GaAs/BC₂N (1桁高い) 非共面構造 :
振幅:~20 nA(実験で測定可能) 角度依存性:J x = G x y ′ y ′ ( 1 − R ) ( 1 − e − α d ) w I J_x = G_{xy'y'}(1-R)(1-e^{-\alpha d})wI J x = G x y ′ y ′ ( 1 − R ) ( 1 − e − α d ) w I
ここで σ x y ′ y ′ = cos 2 θ σ x y y + sin 2 θ σ x z z + sin 2 θ σ x y z \sigma_{xy'y'} = \cos^2\theta\sigma_{xyy} + \sin^2\theta\sigma_{xzz} + \sin2\theta\sigma_{xyz} σ x y ′ y ′ = cos 2 θ σ x yy + sin 2 θ σ x zz + sin 2 θ σ x yz パターン:主にcos²θ + sin²θ形式 共面構造 :
振幅:~0.2 nA(2桁小さい) 角度依存性:sin(2θ)形式(σˣʸᶻから来ている) 実験的実行可能性 :
光子エネルギー0.15 eVは実験的に到達可能な範囲 20 nA電流は測定ノイズレベルをはるかに超える 2つの構造の振幅と角度分布の差は顕著 |m|を0.13から2.0 μᵦまで系統的にスキャン:
光電流と磁気モーメントの単調関係を検証 非相対論的メカニズムの支配的役割を確認 顕著な応答に必要な磁気モーメント閾値(~1 μᵦ)を定量化 速度ゲージ計算は合計を切り詰める必要がある:
300本のバンドを使用して収束を確保 k網格19×15×21(Wannier補間より粗いが対称性を保証) 非共面Mn₅Si₃の誘電テンソル:
許可される非対角成分εˣʸ:20ε₀ (ℏω1 eV) 禁止成分εᶻˣ, εʸᶻ:<1ε₀ スピン対称性制約の有効性を再度検証 非共面:18個の独立成分中、12個がスピン対称性で許可 共面:18個の成分すべてがスピン対称性で禁止 数値結果は理論予測と成分ごとに一致 交替磁性 (Altermagnets) Šmejkal et al., PRX 2022 補償共線磁性 + {T, R}結合対称性 偶数回転対称のスピン分裂フェルミ面 代表材料:RuO₂, CrSb p波磁性 (p-wave magnets) Hellenes et al., 2024 補償非共線磁性 + P破れ + T τ対称性 2つのクラス:I類(共面)とII類(非共面) 本論文が初めて系統的に光学特性を研究 基礎理論 :Litvin & Opechowski (1974) - スピン群概念現代的発展 :Liu et al. (PRX 2022) - 弱SOC材料におけるスピン群対称性系統的分類 :Xiao et al. (PRX 2024) - 5万以上のスピン群タイプを識別、近2000種の化合物をカバーシフト光電流理論 von Baltz & Kraut (PRB 1981) - 理論的基礎 Sipe & Shkrebtii (PRB 2000) - 2次光学応答 Tan et al. (npj Comput. Mater. 2016) - トポロジカル特性との関連 高応答材料 Weyl半金属TaAs:~10⁻⁴ A/V² Osterhoudt et al., Nat. Mater. 2019 強誘電体BaTiO₃:~10⁻⁵ A/V² Fei et al., PRB 2020 WS₂ナノチューブ:増強された本質的光電圧効果 Zhang et al., Nature 2019 磁気構造の議論 :
Brown et al. (1992):非共面構造 Biniskos et al. (PRB 2022):共面構造 輸送特性 :
大きな異常ホール効果 Sürgers et al., AIP Adv. 2016 トポロジカルホール効果 Sürgers et al., Nat. Commun. 2014 異常ネルンスト効果 Badura et al., 2024 初めて スピン群対称性を光応答解析に系統的に適用初めて p波磁性体の大きなシフト光電流を予測初めて 光学測定に基づいて磁気構成を区別する実験スキームを提案スピン光電流に拡張し、完全な理論フレームワークを確立 理論的貢献 :弱SOC材料では、スピン群対称性は磁気群対称性より光応答に対してより厳格な制約を提供 電荷応答はスピン対称性の空間部分から構成される有効点群によって決定される スピン応答はスピンと空間操作の両方を同時に考慮する必要がある 材料予測 :Mn₅Si₃の非共面構造は巨大なシフト光電流(~10⁻⁴ A/V²)を生成 共面構造の応答はスピン対称性によってほぼ完全に抑制される 2つの構造の光電流差は実験的鑑別手段として利用可能 物理的洞察 :非共線磁性体の光応答は磁気モーメント大きさと強く相関 p波磁性体はWeyl半金属に匹敵する光電圧効果を生成可能 非相対論的メカニズムが弱SOC系の光応答を支配 理論的仮定 :弱SOC近似に依存(スピン-格子解耦) 強いSOC材料(5d遷移金属化合物など)には適用不可 励子効果と多体相互作用を考慮していない 計算面 :速度ゲージは多くのバンド(300本)の収束を必要とする k網格は相対的に粗い(19×15×21)で精度に影響する可能性 反射率R(ω)を計算していないため、実際の電流が過大評価される可能性 実験的検証 :Mn₅Si₃低温光電流の実験データはまだない サンプル品質と表面状態が測定に影響する可能性 低温(T<60K)と中赤外光源が必要 材料的限界 :1種類のp波磁性体のみを研究 他の軽元素磁性化合物を系統的に探索していない エネルギー差(5 meV/atom)はDFT誤差範囲に近い 理論的拡張 :線形応答(異常ホール、ネルンスト効果)への推広 集団励起(プラズモン、マグノン)が光学特性に与える影響の研究 スピンゼーベック効果とスピンネルンスト効果の考慮 材料探索 :5万以上のスピン群化合物から高応答材料を探索 他のp波磁性体(Fe₃Sn₂, Mn₃Geなど)の研究 特定のスピン対称性を持つ人工構造の設計 実験的提案 :低温偏光角依存光電流測定 2次高調波生成実験による対称性検証 超高速動力学の時間分解分光研究 応用前景 :スピン対称性に基づく光電圧デバイス設計 純粋なスピン光電流の生成と制御 不揮発性光学メモリ 概念的革新 :スピン群対称性を光応答に適用する先駆的な取り組みで、新しい分析パラダイムを確立 非従来型磁性体における対称性-応答関係の新しい法則を明らかに 弱SOC磁性材料に従来の磁気群を超える理論的ツールを提供 方法の厳密性 :単純なモデルから複雑な材料への段階的な論証 速度ゲージ計算によるWannier関数の対称性破れ問題の回避 完全なテンソル解析(18個の独立成分を1つずつ検証) 結果の説得力 :理論予測と数値計算の高い一致(3~4桁差) 明確な実験判定基準を提供(20 nA vs 0.2 nA) 合理的なパラメータ設定で実験的に到達可能な範囲 科学的意義 :Mn₅Si₃磁気構造争議の鑑別問題を解決 p波磁性体の巨大な光電圧効果を予測 2000以上の化合物の光学特性研究に理論的基礎を提供 執筆品質 :論理が明確で、動機から結論まで層次が分明 図表が精美で、核心情報を効果的に伝達 補足資料が詳細(完全な導出、すべてのテンソル成分) 実験的検証の欠如 :核心予測(20 nA光電流)はまだ実験的に検証されていない Mn₅Si₃の真の磁気構造はまだ議論がある 実験グループとの協力による検証作業を推奨 普遍性の限界 :詳細に研究した材料は1種類のみ 強いSOC系には適用不可(5d/4f材料を制限) 温度効果と相転移挙動について議論していない 計算精度 :k網格が粗い可能性があり、鋭い特徴に影響する可能性 GW補正または励子効果を考慮していない エネルギー差(5 meV/atom)は方法誤差に近い 物理的メカニズム :磁気モーメント-光応答関係の微視的起源の説明が不足 具体的な遷移チャネルの寄与を分析していない 実験的に観測可能な量(スペクトル重み)との直接比較が欠ける 実用性の考慮 :低温(60K)と中赤外条件が実験難度を増加 サンプル製造と表面効果について議論していない 既存技術(MOKE)との比較が欠ける 短期的影響 :
Mn₅Si₃実験研究に明確な指導を提供 スピン群対称性の輸送と光学における役割への関心を喚起 p波磁性体の実験的探索を推進 長期的影響 :
非従来型磁性体の光学特性の新しい分類フレームワークを確立 スピン対称性に基づく材料設計方法を触発する可能性 5万以上のスピン群化合物の特性予測にツールを提供 潜在的応用 :
磁性光電圧デバイス(効率はWeyl半金属に匹敵) 全光スピン操作技術 不揮発性光学ストレージ 再現性 :
理論導出は完全で再現が容易 計算詳細が充分(VASPパラメータ、k網格など) コードはオープンソースツール(WannierBerri)に基づく 補足資料がすべての原子座標を提供 理想的な応用シーン :
軽元素(3d)非共線磁性体の光応答予測 複数の候補磁気構造の実験的鑑別 弱SOC材料の対称性解析 p波磁性体と交替磁性体の分類研究 不適用シーン :
強いSOC系(Ir、Pt基化合物など) 高温応用(方法は秩序磁気構造に依存) 量子化応答が必要なトポロジカルシステム(本論文は古典応答を研究) 超高速動力学(本論文は定常状態理論) 推広の可能性 :
他の応答関数(ホール電導、熱電係数)への拡張が可能 反強磁性スピントロニクスに適用可能 機械学習による高スループット探索と組み合わせ可能 Šmejkal et al., PRX 12, 031042 (2022) - 交替磁性理論の基礎Liu et al., PRX 12, 021016 (2022) - 弱SOC材料におけるスピン群対称性Hellenes et al., arXiv:2309.01607 (2024) - p波磁性体の分類Sipe & Shkrebtii, PRB 61, 5337 (2000) - 2次光学応答理論Ibañez-Azpiroz et al., PRB 97, 245143 (2018) - Wannier補間によるシフト電流計算Brown et al., J. Phys.: Condens. Matter 4, 10025 (1992) - 非共面構造Biniskos et al., PRB 105, 104404 (2022) - 共面構造予測Sürgers et al., Nat. Commun. 5, 3400 (2014) - トポロジカルホール効果総合評価 :これは概念、方法、応用レベルで重要な革新を持つ高品質な理論物理論文である。実験的検証が欠けているが、理論フレームワークは厳密で、予測は明確であり、非従来型磁性体研究に新しい方向を切り開いている。高影響力ジャーナル(PRX、Nature Physicsなど)への掲載を推奨し、核心予測の検証のための実験協力を早期に開始することを提案する。