1110■ 研究紹介DOSUE-RR: ミリ波を用いたダークマター探索実験
京都大学

安達 俊介
adachi.shunsuke.5d@kyoto-u.ac.jp

隅田 土詞
toshi.sumida@cern.ch

鈴木 惇也
suzuki.junya.4r@kyoto-u.ac.jp
2022 年 (令和 4 年) 11 月 1 日

はじめに

1

DOSUE-RR は「どすえ-ダブルアール」と読む。 京風ダークフォトンの“隠れた” 電場

キラキラネームだが, “Dark-photon dark-matter Ob-

金属板

serving System for Un-Explored Radio-Range” の略である。DOSUE-RR 実験は, 冷たいダークマター (ColdDark Matter, CDM) を説明するようなダークフォトンを電波 (特にミリ波) 計測によって検出することを目指している。Simons Observatory や GroundBIRD といった宇宙マイクロ波背景放射 (CMB) 望遠鏡の開発研究で培った技術をもとに 2020 年から京都大学で実験を進めている。本記事では， DOSUE-RR の最初の探索結果 [1] と今後の展望について解説する。

ダークフォトンとその探索方法

22.1

ダークフォトンとは?

ダークフォトンは標準模型のラグランジアンに新たなU(1) 対称性を導入することで, 光子と弱く相互作用する未知粒子として現れる [2]。ダークフォトンは質量を持ち，WIMP (Weakly Interacting Massive Particle) やアクシオンとは異なる CDM 候補粒子になる。DOSUE-RR は，ダークマターの正体がダークフォトンだとして地上で検出しようとする実験である。

2.2

ダークフォトン CDM の転換光とその特徴

ダークフォトンは光子との結合定数 χ を持ち，χ に比例する “隠れた” 電場を持っている。そのダークフォトンの持つ電場 EDP が金属表面のような電磁気的境界面を横切る際に実光子 (“転換光”) を放出する [3]。転換転換光θ ≲ 0.06∘
図 1: ダークフォトンの隠れた電場が金属と空気との境界面で転換光に変換される。

光を放出する理由は、金属表面での電磁気的境界条件(表面に平行方向の電場が 0, EDP,|| + E実光子,|| = 0) と,ダークフォトンと転換光の間のエネルギー保存を満たす必要があるからである。 これら 2 つの条件から、転換光は 2 つの特徴を持つ。

1 つ目の特徴は，放出される転換光の方向が金属板に対してほぼ垂直となることである。電磁気的境界条件からダークフォトンの運動量の平行成分が保存され，かつエネルギー保存則からダークフォトンの質量エネルギーが転換光の運動量の垂直成分に変換される。ダークフォトンが CDM であれば，その運動エネルギーは質量エネルギーより小さいため1 転換光の運動量の平行成分より垂直成分が大きくなる。よって，図 1 に示すように転換光はほぼ垂直に放出される [5]。

2 つ目の特徴は，転換光の周波数がダークフォトンの質量に一対一対応することである。エネルギー保存則から導かれる転換光のエネルギーは，ダークフォトンの1 地球と CDM の相対速度はおよそ 220 km/s である [4]。

2111Δν ∼ 10−6ν0∼ 20 kHz

2.3

ダークフォトン転換光の 2 つの特徴によって，金属板から垂直に放出される転換光を捉えて分光し，得られ信号強度 [ 10-19 W ]信号のシミュレーション@ mDP=80 μeV/c2た周波数スペクトルの中から信号のピークを探すこと

PDP ∝ χ 2+

ν0 =++

周波数 [kHz + 20 GHz]
mDPc 2h+で，ダークフォトンを探索することができる。その質量mDP は検出される周波数で決まり，ダークフォトンと+光の結合定数 χ の 2 乗は得られた信号強度の大きさに比例する。図 3 に 2 つのパラメータ (mDP , χ) に対する+制限を示す。広い範囲に渡って宇宙観測から制限がつけられている。10–300 GHz にわたるミリ波の周波数帯域= 20 GHzでは宇宙観測からの制限が弱く, これまでに探索実験

図 2: ダークフォトンの転換光を測定した場合の周波がほとんど無い。つまり, 未開拓領域が広大にある。こ数スペクトルのシミュレーション。質量が 80µeV/c2 のの周波数領域は CMB の観測周波数帯と合致している。
DOSUE-RR は CMB 望遠鏡で利用されている技術を応

ダークフォトンを仮定した時， 20 GHz にピークをも15ち, 幅 20 kHz 程度の信号となる。

IPNS Seminar, 2022/06/17用して，ダークフォトンの検出を目指す。

この転換光の特徴を利用した探索方法として, 元々パ運動量を無視するとその質量で決まり，hν0 = mDP c2ラボラアンテナを用いることで, 垂直に放出される転となる。例えば, mDP = 80 µeV/c2 であれば, ν0 =換光を一点に集めて転換光を検出する方法が考案され20 GHz となる。もちろんダークフォトンの運動量はゼた [3]。我々の実験では金属平板を置き，そこから出てロではないので, 転換光の周波数はその運動量分だけくる転換光をホーンアンテナで集めて測定する (図 4)。

ν0 からずれ, また運動量分布に起因するわずかな幅を金属表面のどのくらいの面積からくる転換光が集められもった周波数分布となる。ダークフォトンの運動量分布るかは，アンテナの実効的な開口面積 Aeﬀ [cm2 ] で決まは一般的に Maxwell-Boltzmann 分布 (熱平衡状態気体)

る。得られる信号の強度 PDP [W] は，Aeﬀ だけでなく，で仮定され，∆ν/ν0 ∼ (vDP /c)2 ∼ 10−6 になる2 。つま

結合定数 χ, ダークマターのエネルギー密度 ρDM (=り, 転換光はかなり細いピーク信号になる。図 2 に期待

0.39 ± 0.03 GeV/c2 [8]) から推定できる [3]:される周波数スペクトルを図示した。
2v

PDP =1.1 × 10−19 W()( χ )2ρCDM
Aeﬀ×. (1)10−1017.3 cm2 0.39 GeV/cm3

DP はダークフォトン CDM の速さ。

パラボラアンテナを用いた実験

光子との結合定数 χ

転換光の周波数 ν0 [GHz]Axion 共鳴探索実験による制限

よる制限

CMB にダークフォトン CDM 探索の現状

金属板
太陽寿命による制

広大な限未開拓領域

信号光(ミリ波)
ホーンアンテナ
低ノイズアンプ

熱放射

受信機雑音

(アンプなどから)スペクトラムダークフォトンの質量 mDP [μeV/c2]

図 3: ダークフォトンの質量 mDP と光子との結合定数 χ に対する制限。

その質量は転換光の周波数に一対一対応するので，x 軸下部に質量 mDPを，x 軸上部に周波数 ν0 を記載した。10 GHz 以上では宇宙観測よりも強い制限を与える実験結果はほぼ皆無であった。唯一ある細い縦線の制限は,KUMODeS と呼ばれる CMB 観測技術の気象応用で開発されたミリ波受信機を用いた先行研究である (といっても周波数幅は ∆ν = 10 MHz しかない)。[6, 7] を参照されたい。

アナライザ

図 4: DOSUE-RR の実験概念図。金属板からの転換光をアンテナで集光し， アンプで増幅したのちにスペクトラムアナライザで分光測定する。雑音には 2 種類あり, アンテナに入ってくる熱放射雑音と, アンプなどの機器で生じる受信機雑音がある。

3112

後述する我々のセットアップでは，Aeﬀ = 17.3 cm2 であり，目指す χ の精度を 10−1910−10金属板と設定すると，PDP ∼W となる。これは非常に小さな信号強度である。

窓

例えば, ダークフォトン転換光と同じ帯域幅（20 kHz）

真空層

で CMB の強度を計算すると, 10−18 W である。つま30 K 層

り, CMB よりも微弱な信号を従来の観測装置よりも桁違いに高い周波数分解能で検出する必要がある。

3K層
1.5 m
DOSUE-K
クライオスタット

初号機: DOSUE-K

3K-band 受信機の開発と探索

3.1

ホーンアンテナ受信機の開発

前述の図 4 のように, 金属平板から垂直に放出される信号光をホーンアンテナで集光し電気信号に変換する。
微弱な信号を読むために電気信号をアンプで増幅し，その後にスペクトラムアナライザでフーリエ変換し，周波数スペクトルを測定する。測定の際に問題となる雑音は2 種類あり，ホーンアンテナからスペクトラムアナライザまでの過程で混入する雑音（受信機雑音）と，ホーンアンテナに入ってくる信号光以外の光 (熱放射) による雑音がある（図 4）。受信機雑音は電子回路内での熱雑音が主原因であるため，低雑音なアンプを極低温 (3 K) で動作することで抑制する。一方，熱放射は金属板やホーンアンテナをすべて極低温クライオスタット内に入れてしまえば， 外界からの熱放射を遮断できるが，装置のゲイン較正を高頻度でおこなうことが困難になる。そこ

で，我々のセットアップではポリエチレン製の窓をクライオスタットに搭載し， 窓の外にアルミ製の金属板を置くことによって, 金属板部分をミリ波光源に変更しておこなうゲイン較正を可能にした。常温に置かれた金属板自身からの熱放射を気にされる読者もいるかもしれな
いが, 金属は一般に放射率 (要するに完全な黒体に近いかどうか) が低いため, 金属板自身からの熱放射は少なくて済む。

クライオスタットは円筒形であり，円筒上側に信号を入れるための窓，下側に冷凍機を設置している。内部構
造としては外側から真空層，30 K 層，3 K 層という 3 層構造になっており，3 K 層の内壁は電波吸収体で覆われている。つまり, 窓のある極低温電波暗室のようになっている。アンテナ，低雑音アンプは最内層に配置され，3 K まで冷やされるようになっている (図 5)。このクライオスタットは，縦に長い円筒形（約 1.5 m）に設計した。ホーンアンテナが取り込むことができる信号光はアンテナ開口の面積で制限されるので， アンテナが見込む金属板の立体角を大きくしても転換光の信号強度は増えない。この事を利用して，外来の熱放射をさらに抑制する。窓の位置をアンテナから十分離すことで，アンテ信号光

GM 冷凍機図 5: 受信機全体の写真 (左図) と断面図 (右図)。クライオスタット上部の窓の上に設置された金属板の底面でダークフォトンが光に転換され，その転換光をクライオスタット下部にあるホーンアンテナで捉える。金属板からホーンアンテナまでの光路の大部分が 3 K 層で覆われており, 熱放射を抑制する構造となっている。

ナが望む窓の立体角を小さくし, 転換光の強度を維持しつつ熱放射を抑制した。

クライオスタットの組み立て手順としては，冷凍機に真空層, 30 K, 3 K 層のそれぞれの底面を取り付けた上で， 3 K 層，30 K 層，真空層の円筒を内側から順に組み付けていった (図 6)。3 K 層，真空層は上下の 2 段で分解できるようにしていたが， 30 K 層のみ 1 m ほ

どの高さの円筒を一体物で作ってしまったために, 手で持ち上げるのが難しいという困難を生んでしまったが,力を合わせて持ち上げてなんとか組み立てることができた。(現在は，クレーンが使える場所に移動して困難は解消されている。)

3.2

初探索と結果ダークフォトンの初探索は 2021 年 11 月 29 日から 12月 10 日の 12 日間にわたっておこなった。当時修士 2 年の小高駿平氏 (現在は民間企業に就職) が主導した。ス

ペクトラムアナライザは 1 回の測定でわずか 2 MHz の範囲しか測定できないので，測定周波数を 18 GHz から26.5 GHz まで少しづつずらして測定をおこなった。各周波数帯域に対する実効測定時間は 24 秒であった。おおよそ 40 分かけて 100 MHz の範囲を測定し， その前後にゲイン較正をおこない， その較正を元にスペクトラムアナライザでの測定値をアンテナに入ってきた信号の強度に変換した。12 日間かけて，18–26.5 GHz の範4113

30 K 層3K層30 K 層真空フランジ
図 6: DOSUE-K クライオスタットを組み立てている様子。30 K 層は長さ 1 m 程もあり, 組み立てに人手が必要である。

囲に渡る周波数スペクトルを得た。

取得したスペクトルデータに対してピークの高さ PDPをフリーパラメータとしたフィッティングを周波数毎におこなった。解析をおこなう際のスペクトルデータのビン幅の決定や誤差の推定などに, CMB のデータ解析手法を適用した (例えば, ヌルサンプル法と呼ばれるノイズレベルや系統誤差を評価する手法)。図 7 に PDP のフィッティング値とそれらに対する null-signal-hypothesis (信号が無い仮定) における p-value を示す。各フィッティングの p-value が低い周波数も見受けられたが， フィッティングの試行回数が非常に多い (4 × 106 回) ため， 全フィッティングに対する p-value (global p-value, いわゆ

る look-else where eﬀect を考慮した p-value) は 68 %であった。つまり，信号がないという仮定と矛盾がない
ため，今回の探索では信号の兆候はなかったと結論づけた。

それでも万が一, 「もうちょっと統計精度を上げたら発見があった！」なんてことがあったら、悔やんでも悔やみきれない。そこで, さらなる努力もおこなった。具

体的には，p-value が小さい 27 箇所の周波数領域に対して, 10 倍の測定時間をかけて追測定をおこなった。そして，あらためて信号がないことを確認した。図 7 には追測定後の結果も赤色（白黒印刷では濃い色）で重ね書きした。

以上の測定結果を式 (1) を用いて χ への制限に焼き直したものを， 図 8 に示す。この質量領域で, 世界最高感度の探索結果である。また, CMB による χ への制限より 1 桁以上小さい制限を与えることができた。

図 7: 各周波数でのフィッティングから得られたピークの高さ (PDP , 上図) と local p−value (plocal , 下図)。

plocal < 10−6 であった 27 点に対しては，追測定を 10 倍時間をかけておこなった。追測定の結果を赤色 (白黒印刷では濃い色) で示す。

今後の戦略

44.1

質量領域の拡大:DOSUE-Q/V/W/D/Y

DOSUE-RR の探索はこれで終わらない。さらなる未開の周波数領域を開拓することを計画している。各周
波数帯域でミリ波の規格が Q-バンド, V-バンド, W-バ

ンド, D-バンド, Y-バンドと名付けられており，それに合わせて DOSUE-Q, DOSUE-V, DOSUE-W, DOSUE-

D, DOSUE-Y と受信機を改造し， 探索を行っていく予定である3 。110 GHz を境に受信機に利用する技術が異なるのでそれぞれ説明していきたい。

DOSUE-Q/V/W
K-バンドよりも一つ上の Q-バンド (33–50 GHz) での測定は，基本的には K-バンドと同様の手法で行える。ただし, スペクトラムアナライザの対応周波数が 26.5 GHz

以下に限られているため4 , 図 9 に示すようなミキサを

3 東北大の K 氏に DOSUE-RR/Q, DOSUE-RR/V などという言い方にすべきではないかと言われたが，RR=radio-range は非常に広い範囲の周波数帯域を表す言い方であるので，ここを各バンドに置き換えて DOSUE-Q = “Dark Photon Observing System for
Un-Explored Q-band” というように言っていると理解してほしい。
4 高い周波数に対応したスペクトラムアナライザも存在するが非常に高価になる。

5114

よる制限

太陽

DOSUE-V
DOSUE-W
DOSUE-D
DOSUE-Y

寿命

DOSUE-Q

DOSUE-K

CMB にDOSUE-J

光子との結合定数 χパラボラアンテナを用いた実験

転換光の周波数 ν0 [GHz]による制限

今回の探索結果

Axion 共鳴探索実験
による制限

ダークフォトンの質量 mDP [μeV/c2]
図 8: ダークフォトンの質量 mDP と光子との結合定数 χ に対する制限。DOSUE-RR の初探索である 18–26.5 GHzでの測定 (DOSUE-K) で CMB より強い制限を与えた。今後は，帯域ごとに受信機を改造して, より広い帯域に渡って探索していく。

実現を目指している。

ダークフォトン
金属板

さらなる高みへ... DOSUE-D/Y
入力帯域33-50 GHz

ホーン
アンテナ

低ノイズ
低温アンプ

DOSUE-RR では D-バンド (110–170 GHz) と Y-バ

スペクトラムアナライザ

ンド (170–260 GHz) での探索も計画している。しかし同軸ケーブル

ミキサ
常温アンプ (ダウンコンバータ)

導波管

クライオスタット

ながら， 110 GHz 以上の周波数になると，市販の低ノ差周波帯域
3-20 GHz

基準信号(LO)
30 GHz
信号発生器

イズアンプの雑音レベルが劇的に (典型的には一桁程度)悪くなり，価格も非常に高価になる。つまり，アンプで増幅してミキサでダウンコンバートする手法は望ましくない。

この解決策として電波天文学で導入されたのが, SISミキサと呼ばれる超伝導デバイスである。SIS ミキサは

図 9: Q-バンドで測定のダイアグラム。

「トンネル効果」, 「粒子性」(単一光子による Photon-用いた周波数のダウンコンバージョンが必要となる。具

いう量子的性質を統合した機能を提供するデバイスであ

体的には, 入力帯域に対応したアンプで信号を増幅した後にミキサに入力し, 信号発生器で設定した基準信号との差周波をスペクトラムアナライザで読み出す。このような手法を用いることで，より高い周波数 (V-バンド:50–75 GHz, W-バンド: 75–110 GHz) においても，適切な周波数の基準信号を使用して, 測定可能な低い周波数に変換することができる。

高周波になるほど同軸ケーブルでの減衰が大きいため，ミキサの手前まで導波管を用いて信号を取り回す必要があるが， 低温層の中での断熱を適切に行いながらこれを行う。現在，断熱導波管，フレキシブル導波管などを組み合わせて, 低温層内の綿密な設計をおこない,温度による収縮なども考慮した減衰の少ない信号測定のAssisted Tunneling 電流), 「波動性」(周波数変換) とり, 低ノイズ（雑音温度換算で 50–100 K）に周波数をダウンコンバートできる [9]。ダウンコンバートできれば，低周波数用の低ノイズアンプを用いることができ，高感度に測定できる。現在は，Y-バンド向けのプロトタ
イプ受信機の準備をおこなっている (図 10)。今年度中に，SIS ミキサを初めて用いた探索を Y-バンドでおこなう計画である。

6115で 4 GHz 以下の信号にダウンコンバートした後，開発す金属板

るスペクトロメータでフーリエ変換・データ蓄積の処理

信号光
窓 φ8 cmを行う。その後イーサネット経由でコンピュータにデータを送り，解析処理を行う。

信号発生器(~200 GHz)

m40 c

クライオスタット
内部へ

このシステムの心臓部であるフーリエ変換／データ蓄積処理部は Xilinx 社製の RFSoC で構築する。RFSoC

は高速 ADC・DAC，FPGA，CPU などが一体になったチップで，5G 通信などを背景に生まれた新しい技術である。4 GSPS を超えるサンプリング速度の ADC で広い測定帯域幅を実現し，FPGA で並列データ処理を行う。

DOSUE-Y 用
クライオスタット

この開発で最もチャレンジングなのは， RFSoC の

図 10: DOSUE-Y 受信機のプロトタイプ。上の窓から転換光を導入し，中にあるホーンアンテナで集光して計測する。SIS ミキサには基準となる高周波信号を入れる必要があり，それを右側の信号発生器で生成し，クライオスタット内に導入している。

4.2

FPGA 部分のファームウェア開発である。4 GSPS での超高レートサンプリングによりデータ量が膨大になるため， dSpec では FPGA でフーリエ変換・蓄積処理を行いデータ量を削減した後，PC へとデータ転送する。

FPGA は CPU のような高速動作はできないため，帯域 4 GHz のデータをさばく工夫が必要となる。この問題を解決するため，高速フーリエ変換の回路を 16 個並列に用意し，最後に並列回路の結果を組み合わせること統計精度の向上戦略によって, 広帯域のフーリエ変換を実現した。組み合わRFSoC エレクトロニクス「dSpec」の開発

ミリ波周波数領域でのダークフォトンやアクシオン

ダークマターの探索では，図 2 に示したように周波数スペクトルを測定して特徴的なピークを見つける探索手法が有効である。ただし, 一般的な市販のスペクトラムア

ナライザでは, デッドタイムが 99% を占め， 一気に測定できる周波数帯域が狭いという致命的な欠点がある。
この欠点により, 測定の効率が著しく低い。さらに, 市
販品は必要以上に機能が多く (つまり, イラナイ機能がせ部分と後段のデータ蓄積部分では高位合成 (HLS) を駆使し，迅速なファームウェア開発を成し遂げることができる。

本稿の執筆時点では, 修士 1 年の竹内広樹氏の努力によって, dSpec のファームウェアが整った。性能の評価をおこなった後に, DOSUE-RR 実験に導入する予定である。今後の発展計画として, 帯域をさらに倍化することも検討している。

RFSoC

多く)， 価格も高い。
スペクトラム
アナライザ

IQ
ミキサ

フーリエ変換／データ蓄積

RFSoC

置き換え
10 - 40GHz~イーサネット< 4GHzPCdSpec信号発生器

図 11: デッドタイムなし・広帯域測定を可能にするデータ取得システム dSpec の概略図。ダークフォトン由来の信号をミキサで ダウンコンバートした後， RFSoCでフーリエ変換・データ蓄積の処理を行う。

この問題を解決するため，デッドタイムなしで広帯域測定ができる測定装置 “dSpec”(ディースペック,darkmatter Spectrometer) の開発プロジェクトを立ち上げた。図 11 に開発するエレクトロニクスの概略図を示す。ダークフォトン／アクシオン由来の信号をミキサ

LO用
信号発生器

テスト用
信号発生器

IQミキサ

図 12: dSpec のテスト中の様子。二つの信号発生器とIQ ミキサを組み合わせて模擬信号を生成し， RFSoCの搭載されたボードに入力している。

7116

アクシオン

磁場

S

信号光(ミリ波)
DOSUE-RR
受信機
図 13: DOSUE-RR 受信機を用いてアクシオン探索をする場合のセットアップの概念図。

アクシオン探索への拡張展望

4.3

仲間を常に募集している。周波数バンドの数だけ新たな実験結果を出すことができ, それを主導できるチャンス

金属板N

しく，全く人手が足りていないのが実情であり，新たなDOSUE-RR の測定セットアップに磁場を組み込むことでアクシオン探索への転用も可能である。アクシオンもダークフォトンと同様に軽いダークマターの候補である。磁場が存在すると，アクシオンと磁場が結合して金属境界面からわずかに転換光を生成する。アクシオンから生成される転換光もエネルギー保存によって同様のミリ波信号になる。

従来のアクシオン探索実験では，主に共振器を用いて信号を増幅して探索する方式が取られており，数 GHz以下での探索が精力的におこなわれてきた。共振周波数が高くなるほど, 共振器には高い工作精度を必要とがある。今後，どんな実験しようか考えている学生，若手研究者にはぜひ興味をもってもらいたい。そして, 筆者達まで気軽にご連絡いただきたい。
本 記 事 で 紹 介 し た 研 究 は,科 研 費 20K14486,20K20427, 21H01093, 21H05460, 22K18713 と村田財団, 住友財団, 京都大学白眉プロジェクトの助成のもとにおこなっています。

参考文献[1] S. Kotaka et al., “Search for dark photon cold
dark matter in the mass range 74–110 µeV/c2
with a cryogenic millimeter-wave receiver”,
arXiv:2205.03679 (2022).
[2] P. Arias et al., “WISPy cold dark matter”, Journal of Cosmology and Astroparticle Physics 06(2012) 013.
[3] D. Horns et al., “Searching for WISPy cold dark
matter with a dish antenna”, Journal of Cosmology and Astroparticle Physics 04 (2013) 016.
[4] F. Mayet et al., “A review of the discovery reachof directional Dark Matter detection”, Physics
Reports 627, (2016) 1になっても同じ実験手法で探索することが可能である。

[5] J. Jaeckel and S. Knirck, “Directional resolution of dish antenna experiments to search for
WISPy dark matter”, Journal of Cosmology andAstroparticle Physics 01 (2016) 005.

増す。この強みを活かした実験の準備を東北大学の岸本[6] N. Tomita et al., “Search for hidden-photon cold
dark matter using a K-band cryogenic receiver”,

し, 一般に性能が低下する。一方で，DOSUE-RR のように金属板表面からの信号を検出する場合は高い周波数
つまり, 高い周波数ほど従来手法と比較した際の強みが康宏氏と共に進めている。図 13 に DOSUE-RR 受信機を用いたアクシオン探索実験のセットアップの概念図を示す。

5まとめ

DOSUE-RR の初探索結果と今後の戦略について解説した。DOSUE-RR は，着実な成果を上げている小規模

Journal of Cosmology and Astroparticle Physics
09 (2020) 012.
[7] 田島治, “ミリ波でみる大気と宇宙 ビッグバンとゲリラ豪雨, そしてダークマター”, 応用物理 91 巻(2022) 1 号

数でおこなっており，自身のアイディアで様々に工夫す[8] R. Catena and P. Ullio, “A novel determinationof the local dark matter density”, Journal ofCosmology and Astroparticle Physics 08 (2010)004.など大規模，中規模の実験も経験しているが，テーブル[9] 中井直正, 坪井昌人, 福井康雄, 「シリーズ現代の実験である。まだまだ始まったばかりの実験でかつ少人ることができる。筆者達は ATLAS 実験や CMB 実験

トップ実験ならではの面白さを日々感じている。少人数なだけでなく，各々が DOSUE-RR 以外の大実験で忙

天文学 第 16 巻 電波天文の観測 II–電波天文学」
日本評論社 (2020)