「C-S-H系高圧超伝導は本物なのか論争」を自分でもチェックしてみる。

【イントロ】
　最近、人類最初の室温での超伝導状態の実現として、R. Dias等によるC-S-H系高圧超伝導論文[1]がNatureに掲載され注目されています。超伝導物質の特定、他グループによる追試、結晶構造の特定などまだまだ明らかにする問題が残っています。しかし、注目されている原因の１つとして、「そもそも最初の論文で報告されたデータは本物なのか？」という批判にさらされている点があります。
　この批判は、アメリカの物理学者、J. E. Hirschによって主になされています。Hirschによれば、電気抵抗の転移が磁場で幅広くならない、磁化率の振る舞いが過去に報告されたEu高圧超伝導の結果に似ている、検証の生データをDiasらに要求しても提供されないなど様々な批判を行っています[2]。一方でDiasらも超伝導の標準理論であるBCS理論を批判するようなトロールであるHirschに提供するデータはない、高圧下の実験をHirschは理解していないなど様々な反論を行っています[3]。この論争は、Hirschの批判論文がDiasらとの私信を勝手に批判論文に含めたことで掲載取り下げになったり、Diasらの共同研究者が関わったEu高圧超伝導論文[4]が撤回されたり、激しい戦いとなっています。
　本記事では、このHirschの最新批判論文[5]について取り上げたいと思います。

【最新の批判論文の内容】
　C-S-H系超伝導の報告以来、Hirschはそのデータの解釈について批判を続けています。その妥当性は別として、批判に答える形としてDias等は実験の詳細や、磁化率の生データを載せた反論論文[6]を2021年末にArxivに投稿しています。最新の批判論文では、Hirschはこの生データを解析し、そのデータが操作(Manipulate)されたものである可能性を指摘しています。

　その内容をみてみます。まず、Nature論文に報告され磁化率と、反論論文に掲載された生データによる磁化率の結果が以下の図です。磁化率の振る舞いは良く似ており確かに生データが掲載されていることが見て取れます。
　ただし、この再現についてヤバいところが１つあります。それは反論論文に掲載された生データが、表形式ではなく図として貼り付けられているため、Hirschらは１つ１つの数値を読み取りデータに直している点です。興味深いのは、この160GPaのデータだけ図で貼り付けられていて、他の圧力のデータはちゃんと表形式として掲載されコピペが可能になっている点です。明らかに意味不明なDiasらの対応の理由はよくわかっていません。
図、磁化率の生データ(左)とNatureに掲載されたデータ(右)


図、これは表形式

図、コレは画像、なんでやねん


　さて、この生データについて、Hirschらは隣の温度との差分をとってみることで、データに不自然な飛びがあることを指摘しています。更にこのトビが、0.16555という値の整数倍になっていることから、Hirschらはこれが生データからよくわからないバックグラウンドの引き方をした結果であり、実際のデータを「Manipulate」したものである可能性を指摘しています。つまり、反論論文で主張されている直線バックグラウンド除去以外の操作を行っているのではないか？、ということのようです。
図、謎のデータ操作が行われていると主張されている160GPaの磁化率(左上)と、
Hirschらがそこから再現した本来の生データ(左下)

　さて、このHirschらの指摘が妥当なのか、試しにDiasらが報告している他の圧力の生データで、同じ様に磁化率の温度差分を調べてみました。もし他のデータでもみえれば、装置の特性の可能性もあるかと考えるからです。
　下の図(上)はDiasらの反論論文[6]に掲載されたコピペ可能な圧力、189GPaの磁化率の生データです。この生データから直線のバックグラウンドを差し引いたのが磁化率となると論文中に記載されています。図(中)は図(上)から直線BGを引いた結果で、たしかに160GPaに似た超伝導転移ぽくなっており、Tc=198Kくらいに読み取れます。このデータについて、温度差分をとったのが図(下)となっています。図からは、Hirschの指摘するようなトビは確認できません。つまり、データのトビは装置の特性ではなさそうです。
　今回の結果は、すべての圧力で飛びが存在するわけではない可能性を意味しています。160GPaのデータはHirschらが頑張って図から読み取って手入力した値なのでミスが含まれているのかもしれません。自分でも追試してみればよいのですが、手入力が面倒臭すぎるので断念です。まさかこれがDiasらの狙い・・・？
　いずれにせよ、なぜか読み取りづらくした160GPaのデータだけ謎の操作(Manipulate)がされている可能性は否定できませんでした。
図(上)、189GPaの磁化率の生データ
図(中)、そこから直線のバックグラウンドを差し引いたデータ
図(下)、さらに横の温度の値との差分をとったデータ

【まとめ】
　HirschらのC-S-H系高圧超伝導に対する批判論文でなされている検証を、自分でもやってみました。Hirschらが疑問視している、なぜか図で貼り付けられた圧力以外の生データの1つでは、Hirschらが見つけたデータの謎のトビは確認できませんでした。これだけでは何も言えませんが、C-S-H系高圧超伝導の著者たちがどのような再反論を示すのか、今後の動向が愉しみです。

【参考文献】
[1] Snider, E., Dasenbrock-Gammon, N., McBride, R. et al. Room-temperature superconductivity in a carbonaceous sulfur hydride. Nature 586, 373–377 (2020).
[2] Jorge E Hirsch, On the ac magnetic susceptibility of a room temperature superconductor: anatomy of a probable scientific fraud, arXiv:2110.12854
[3] Breakthrough or bust? Claim of room-temperature superconductivity draws fire, doi: 10.1126/science.acx9428
[4] M. Debessai, et al., Retraction: Pressure-Induced Superconducting State of Europium Metal at Low Temperatures [Phys. Rev. Lett. 102, 197002 (2009)], Phys. Rev. Lett. 127, 269902 – Published 23 December 2021
[5] Dirk van der Marel, Jorge E Hirsch, Comment on Nature 586, 373 (2020) by E. Snider et al., arXiv:2201.07686
[6] Ranga P. Dias, Ashkan Salamat, Standard Superconductivity in Carbonaceous Sulfur Hydride, arXiv:2111.15017





コメント

コメントを投稿

このブログの人気の投稿

交代磁性があるのかないのかどっちなんだい！～RuO2の場合～

イメージ
 【イントロ】 博士「はぁ～交代じゃ、交代じゃ」 学生「どうしたんですか？キャバクラですか？」 博士「違う違う違う違う今日は違う！交代は交代でも、交代磁性のことを調べておったんじゃ」 学生「交代磁性？」 博士「なんじゃ知らんのか？交代磁性とは、強磁性と反強磁性に次ぐ第三の磁性として最近発見された、新しい磁気状態のことなんじゃ。」 学生「知ってますよ、 Šmejkal（流暢な発音）が提案 してるやつでしょ？強磁性体のようにバンド構造にスピン分裂が生じるけれど、反強磁性体のようにNetの磁性はゼロとなっている、副格子間の反強磁性状態と解釈できる磁性体のことですよね。」 博士「そうじゃ、Šmejkal（流暢な発音）達の提案が大きなきっかけとなった研究分野じゃな。なんじゃ、知っとるじゃないか」 学生「当然です。中学入試程度です。交代磁性って、 MnTe や CrSb 、それに RuO2 で確立された概念じゃないんですか？」 よく見ると(b)がない（ 引用元 ） 博士「 どこ中出身 じゃ。。。うむ、それらの物質でARPESが行われ、交代磁性の特徴の1つであるバンド構造の巨大なスピン分裂が観測されて Nature や PRB に報告されるなど研究が加熱しておる。」 学生「ならいいんじゃないですか？流行物質どんとこい。ビッグウェーブに乗って論文を量産するタイミングですよ。h-index稼ぎましょう」 博士「しかし、そのうちの１つ、RuO2が実は交代磁性体ではなく、非磁性体なのではないかという報告がなされ議論を呼んでおるのじゃ。」 学生「どういうことですか？？？RuO2は 異常ホール効果 が見えたり、 REX（共鳴弾性X線散乱） や 中性子散乱実験 で反強磁性秩序がすでに観測されたんじゃないんですか？？？陰謀ですか？？？？」 博士「それでは今回は、独ヴュルツブルク大学のPhilipp Keßler博士らによる論文、 Philipp Keßler, et al. , Absence of magnetic order in RuO2: insights from μSR spectroscopy and neutron diffraction , arXiv:2405.10820 を読んでいくかのお」 学生「Absence論文、盛り上がるやつですね」 【どんな実験したの？】
続きを読む

狂気！室温超伝導体は脳内にあった！！！

イメージ
 【イントロ】 博士「これはすごい・・・これはすごいぞ！！！」 学生「どうしたんですか？自分の愚かさのことですか？」 博士「そうではない！ワシは出会ったんじゃ！科学史上最高の論文に！」 学生「どんな論文ですか一体」 博士「なんと 室温超伝導体 が見つかったという話じゃ！」 学生「え！室温超伝導体といえば、室温で電気抵抗がゼロになる物質で、世界を支えるエレクトロニクスに不可避なエネルギーの散逸、発熱の問題を一気に解決し、 世界のエネルギー問題を救うキーテクノロジーになる物質 じゃないですか！」 博士「その通りじゃ！」 学生「気になってきました、早く教えて下さい、その論文！」 博士「うむ…驚くなよ、なんと、 生物の脳の中では超伝導状態が実現しておる のじゃ！！！」 学生「・・・」 博士「ふふふ、あまりの衝撃に言葉もないようじゃな。」 学生「あの」 博士「では、今回は、 P. Mikheenko, Possible Superconductivity in the Brain , Journal of Superconductivity and Novel Magnetism, 2018, No 5, p. 1121-1134 （ Arxiv版 ） を読んでいくぞい！」 学生「帰っていいですか？」 【論文誌がすごい！】 博士「まずは、この論文のなにが凄いかといえば、査読を通って論文誌に掲載されておるのじゃ」 学生「ほんとですよ、どうなってんですか」 博士「そして、掲載された Journal of Superconductivity and Novel Magnetismの2022年度のImpact Factor は1.8なんじゃ。これがどういうことかわかるか？」 学生「そ、それは日本の物理学会が発行するJournal of the Physical Society of Japan、いわゆる JPSJのImpact Factor 、1.7を超えています・・・！」 博士「そうじゃな、つまり論文誌としての影響力としてはJPSJよりも大きな雑誌に掲載された論文ということじゃ」 学生「JPSJもっとがんばってくださいよ」 博士「しかもこの論文は8回も引用されておる」 学生「 みなしPRL じゃないですか・・・」 博士「これだけで如何にこの論文が本物であるか納得してもらえ
続きを読む

まだ最小二乗法で消耗してるの？ ～Passing Bablok法について～

イメージ
【イントロ】 与えられたデータの傾向を見る、つまり回帰分析の一つとして、直線でフィッティングする方法が挙げられます。 ココで問題となるのが、「どうやって直線の傾き(と切片)を求めるか」という点です。 代表的な方法が最小二乗法です。目的の直線を $y=ax$ と仮定します。データ点$y$と求めたい直線との差の二乗を誤差として、すべての点に対する誤差の和を最小化する傾き$a$を求める方法です。 具体的には誤差の和を $E=\sum_{i=1}^{n}(y_i－ax_i)^2$ とすると、$a$について最小値が欲しいので、 $\frac{\partial E}{\partial a}=0$ となる$a$を求めれば良いです。整理すると、 $a=\frac{\sum_{i=1}^{n}(x_i－x_{average})(y_i－y_{average})}{\sum_{i=1}^{n}(x-x_{average})^2}$ となります。 　最小二乗法は直観的にもちょうどよい直線を見つけてくれそうで、最強の方法に見えます。しかし、 いくつかの課題 があります。まず一つが、外れ値に弱いという点です。えっちですね。つまり、ノイズ等の影響でデータ点の集合から大きくハズレた点が存在すると、その点に引きずられて傾き$a$が大きく変化してしまうという点です。 　また、別の問題として、最小二乗法は暗にデータが正規分布に従った誤差を持つことを仮定している点が挙げられます。今回最小化すべき誤差をデータ点と直線の差の二乗和として与えましたが、これはデータ点の周りに分散一定の誤差が正規分布している場合に導かれる最尤推定解と対応しています。詳細は PRML を参照。 さて、それではこうした最小二乗法の呪縛から解き放たれた、 回帰モデルは存在するでしょうか？そうしたモデルとしては 主成分分析 等があります 本記事では、そうした手法の中でもノンパラメトリック回帰モデルの一つである”Passing-Bablok法”について、 原論文 に従った解説とPythonでの実装に取り組みました。 【方法】 Passing-Bablok法はその名の通り、PassingとBablok(PB)により提案された、一次の比例関係を比較するノンパラメトリック回
続きを読む