[論文レビュー] PALMA, an improved algorithm for DOSY signal processing

Afef Cherni, Émilie Chouzenoux+1 more|PubMed|Aug 25, 2016

NMR spectroscopy and applications参考文献 43被引用数 6

ひとこと要約

本論文では、近接作用素を用いた分割スキームと最大エントロピーと ℓ₁-ノルムのハイブリッド正則化を組み合わせることで、逆ラプラス変換問題を解く新しいアルゴリズム、PALMA を提案する。この手法により、モノ分散およびポリ分散系の両方において、高速でロバストかつ高精度な拡散分布の再構成が可能となり、従来の手法に比べてノイズ耐性と収束速度に優れる。

ABSTRACT

NMR is a tool of choice for the measurement of diffusion coefficients of species in solution. The DOSY experiment, a 2D implementation of this measurement, has been proven to be particularly useful for the study of complex mixtures, molecular interactions, polymers, etc. However, DOSY data analysis requires to resort to the inverse Laplace transform, in particular for polydisperse samples. This is a known difficult numerical task for which we present here a novel approach. A new algorithm based on a splitting scheme and on the use of proximity operators is introduced. Used in conjunction with a Maximum Entropy and hybrid regularisation, this algorithm converges rapidly and produces results robust against experimental noise. This method has been called PALMA. It is able to reproduce faithfully monodisperse as well as polydisperse systems, and numerous simulated and experimental examples are presented. It has been implemented on the server where users can have their datasets processed automatically.

研究の動機と目的

• DOSY NMR データの従来の逆ラプラス変換手法に内在する不安定性とノイズ感受性を解消すること。
• 最小限のユーザ介入でポリ分散系の拡散分布を正確に再構成できる、ロバストで効率的なアルゴリズムの開発。
• 装置に起因するアーティファクトや信号の重なりがある状況でも、DOSY 分析の分解能と信頼性を向上させること。
• 分布の形状に関する事前仮定を必要とせず、モノ分散およびポリ分散系の両方に適用可能な汎用的ソリューションを提供すること。
• 公開済みのウェブサーバーとオープンソースコードを通じて、自動的でユーザーフrndlyな DOSY データセット処理を可能とすること。

提案手法

• アルゴリズムは、DOSY データ解析における逆問題を解くために、プライマル・デュアル分割スキームを採用する。
• 非滑らか性を扱うために ℓ₁-ノルム正則化と滑らか性を保証する最大エントロピー正則化の両方に近接作用素を用いる。
• 最大エントロピー（滑らかさのため）と ℓ₁-ノルム（スパarsityのため）の2つの正則化を、1つのパラメータ λ でバランスさせる。
• 最適化問題は、ハイブリッド正則化を組み込んだ制約付き最小二乗問題として定式化され、安定性と収束性を保証する。
• アルゴリズムは、特にノイズが多いまたは複雑なスピン系に対しても、最小限のチューニングで自動処理が可能に実装されている。
• 本手法は、http://palma.labo.igbmc.fr で利用可能なウェブサービスとして提供されており、ソースコードも公開済み。

実験結果

リサーチクエスチョン

• RQ1 実験的ノイズや装置アーティファクトが存在する中で、DOSY NMR における逆ラプラス変換をよりロバストに解く方法は何か？
• RQ2 最大エントロピーと ℓ₁-ノルムをハイブリッドに組み合わせた正則化戦略は、拡散分布再構成の精度と分解能を向上させ得るか？
• RQ3 再構成された拡散分布における滑らかさとスパarsityの最適なバランスは何か？また、結果は正則化パラメータ λ に対してどれほど感度を示すか？
• RQ4 分布の形状に関する事前知識がなくても、モノ分散およびポリ分散系の両方で高品質な結果を得られるか？
• RQ5 収束速度、ノイズ耐性、スペクトル分解能という観点から、既存手法に比べて本アルゴリズムはどの程度優れているか？

主な発見

• PALMA は、正則化パrameter λ が 0.01 から 0.05 の範囲で最適化された場合に、広範な λ 値において高速かつ安定した収束を達成する。
• ノイズが存在する状況でも、モノ分散およびポリ分散の両方の拡散分布を忠実に再現する高分解能の再構成が可能である。
• 古典的 MaxEnt（λ = 1）に比べ、PALMA は収束が速く、特にノイズが多い条件下でも優れた結果を示す。
• 純粋な ℓ₁-ノルム（λ = 0）に比べ、ハイブリッド正則化アプローチは誤った特徴を減少させ、信頼性を向上させる。ただし、ℓ₁-ノルム単体でもスパース解は得られる。
• 勾配の非線形性や温度漂移といった装置の不具合に対しても、本手法はロバストである。これらはしばしば仮想のポリ分散性を引き起こす要因となる。
• http://palma.labo.igbmc.fr に設置されたウェブサーバー実装により、最小限のユーザ入力で DOSY データセットの完全自動処理が可能となり、広範な利用可能性と再現可能性が実現された。