A Digital Twin of the FPGA Digital Signal Processing Chain for MKIDs Readout: Root-Cause Analysis and Mitigation of Spurs

宇宙が、目に見える光ではなく、ミリ波という隠れた電波の言語でその秘密をささやいている様子を想像してみてください。この電波は、星がいかにして誕生し、銀河がいかに進化するかという物語を運んできます。これらのささやきに耳を傾けるために、科学者たちはMKID（マイクロ波・キネティック・インダクタンス検出器）と呼ばれる超高感度な検出器を使用します。これらの検出器を、数千個の小さな、超冷却された音叉だと考えてみてください。光のフォトンのひとつが当たると、その音叉のピッチがごくわずかに変化します。これら数千の音叉のピッチを同時に聴き取ることで、天文学者は宇宙の地図を描くことができるのです。

しかし、数千もの音叉の音を同時に聴き取ることは、歓声に包まれたスタジアムの中で、たった一挺のバイオリンの音を聞き取ろうとするようなものです。これを可能にするために、エンジニアたちは「周波数多重化」と呼ばれる巧妙なトリックを使います。それぞれの音叉に固有の無線トーンを割り当て、それらをすべて一本のワイヤーで送り出すのです。課題は、これらのトーンを生成し、変化を聴き取るための電子機器が非常に複雑であり、電光石火のスピードで行われるデジタル信号処理を伴うことです。もしデジタルシステムにわずかなでも不具合があれば、それは本物の宇宙のイベントのように見える偽の「ゴースト信号」を作り出してしまうことがあり、実際にはただのノイズなのです。これが、研究者たちが解決しようとした問題です。彼らは、なぜ自分たちの聴取機器が幽霊を聞いているのか、そしてどうすればそれらを沈黙させられるのかを知る必要がありました。

機械の中の幽霊

アンデス高地の望遠上に設置されているCONCERTO観測装置に取り組んでいたチームは、奇妙なことに気づきました。400個のMKID検出器から送られてくるデータを見たとき、ラジオのノイズのような、鋭くて厄介なスパイク（突起）が見られたのです。これらのスパイクは特定の周波数（763 Hzと1526 Hz）に現れ、宇宙からの微かな信号を見分ける能力を妨害していました。問題は、これらの「スパー（Spur：不要な信号）」がどこから来ているのか誰も分からなかったことです。それらは望遠鏡からでも、冷やされた検出器からでも、ケーブルからでもありませんでした。それらはコンピュータのコードの中に隠れていたのです。

これらのデジタル上の幽霊を追い詰めるために、研究者たちは「デジタルツイン」を構築しました。もし、自動車のエンジンの完璧な仮想コピーをコンピュータの中に作り上げ、最後のボルトや火花に至るまで再現し、現実の世界よりも100万倍速く動かせるとしたら、と想像してみてください。それが、彼らがMKID信号を処理しているFPGA（特殊なコンピュータチップ）に対して行ったことです。実際のハードウェア上でシミュレーションを実行するのに何日も待つ代わりに、彼らはチップの挙動をビット単位、サイクル単位で正確に模倣するPythonベースのモデルを作成しました。この仮想的な双子のおかげで、彼らはわずか1時間で数千回の実験を行うことができました。これは、実際の機械では数日かかる作業です。

探偵の仕事

デジタルツインを用いて、チームはシステムを「ループバック」モードで実行しました。これは、生成した信号をデジタル処理チェーン全体に通し、それを再び入力へと戻すことで、マシンが自分自身の声に対してどのような挙動を示すかを確認するという意味です。これを実行したところ、仮想的な双子は、実際の望遠鏡のデータで見られた763 Hzと1526 Hzの幽霊のようなスパイクを完璧に再現しました。これは大きな手がかりとなりました。つまり、問題は物理的な世界ではなく、完全にデジタルロジックの内部にあるということが証明されたのです。

調査の結果、2つの主要な原因が協力してノイズを作り出していることが判明しました。

1. 不一致なリズム： システムには主に2つの部分があります。「励起（Excitation）」チェーン（トーンを作成するもの）と、「解析（Analysis）」チェーン（それらを聴き取るもの）です。励起チェーンは、16ビットコンピュータの標準である65,536ステップごとに繰り返されるデジタルカウンタを使用してトーンを作成します。しかし、信号がシフトされ処理される際、「バンドシフター（Band-shifter）」が40という数字に基づいた新しいリズムを導入します。65,536のリズムと40のリズムを混ぜ合わせると、数学的に混乱が生じます。結合された信号は、最終的にずっと長く、奇妙に複雑な繰り返しパターン（具体的には、5 × 2の19乗サンプルというサイクル）を持つことになりました。この不一致により、信号はもはや「クリーン」ではなくなり、隠れた複雑な揺らぎを持つようになったのです。
2. 漏れ出すフィルター： 問題の第2の部分は、不要な周波数をフィルタリングするためのふるいとして機能する「デジタル・ダウン・コンバータ（DDC）」でした。このふるいは、65,536サンプルにわたって信号を平均化するように設計されていました。しかし、信号のリズムが上述の奇妙に長いパターンになってしまったため、「幽霊」の周波数がふるいの穴と一致しなくなりました。ブロックされる代わりに、これらの不要な周波数がフィルターをすり抜け、メイン信号へと折り返され、厄介なスパイクを作り出したのです。

解決策：デジタルエンジンのチューニング

リズムの不一致を理解したことで、解決策は明らかになりました。チームは、励起チェーンと解析チェーンが同じビートで踊るようにする必要がありました。

彼らはシンプルかつ巧妙な変更を提案しました。

• カウンタの変更： デジタルカウンタが標準の65,536でラップアラウンド（一周）させるのではなく、65,520でラップアラウンドするように強制しました。なぜこの特定の数字なのか？ それは、65,520が40で完全に割り切れるからです。これにより、トーンジェネレータのリズムとバンドシフターのリズムがついに同期することになりました。信号は今や、クリーンで予測可能な繰り返しパターンを持つようになりました。
• ふるいのサイズ調整： また、DDC内の平均化フィルタを、65,536サンプルではなく65,520サンプルを見るように調整しました。これにより、フィルターの「穴」が幽霊信号の周波数と完全に一致し、それらが完全にブロックされるようになりました。

結果

チームはこの新しい設計を実際のFPGAチップでテストしました。デジタル・ループバックを設定し（一時的に望遠鏡をバイパスして）、古いファームウェアと新しいファームウェアを比較しました。

結果は劇的でした。古いファームウェアでは、ノイズスペクトルに763 Hzと1526 Hzの明確でギザギザした2つのスパイクが現れていました。しかし、新しいファームウェアでは、それらのスパイクは消失しました。ノイズフロアは完全に平坦で静かになり、-240 dBc/Hzに達しました。「幽霊」は退治されたのです。

この勝利のために支払った代償は、コンピュータチップのメモリ使用量がわずかに増加したことだけでした。チップのロジックリソースは約3.2%増加しましたが、チップにはまだ十分な余裕があったため、これは宇宙をクリアに観測するための非常に小さな代償でした。

デジタルツインを構築することで、研究者たちは問題の場所を推測しただけでなく、それを証明し、修正し、検証しました。これにより、CONCERTO観測装置は、自らのデジタルの鼓動による静電気に邪魔されることなく、宇宙の声を聴き続けることができるようになったのです。