[スポンサーリンク]

一般的な話題

計算化学者は見下されているのか? Part 1

[スポンサーリンク]

近年、反応機構をサポートするために計算化学を用いた論文が多数発行されています。しかし、ひとりで実験も計算も両方できる研究者は少ないため、その多くが共同研究であると思われます。

こういった形の共同研究で論文を出す際に問題となりやすいのが、オーサーシップです。

最近の傾向では、実際に実験した人が 筆頭著者になることが一般的です。

では、

 

「実際に計算した人」は 筆頭著者にはなれないのでしょうか??

…計算化学者は、便利屋さんでも有機合成化学者の奴隷でもありません!計算した人も評価されるべきだと考えます。しかし、計算化学というのは、実験化学よりも下に見られることが多いです。(同様のことがバイオインフォマティストと生物学者の間にも言えると思います。…トランスクリプトーム解析とか簡単じゃないんだぜ?

 

こういった問題の背景にあるのは、

多くの人が計算化学の大変さを理解していない

という事実に尽きると思います。

実際に複雑な計算に取り組んだ人ならば、その大変さは実験と同程度と感じていると思います。(ワンステップで終わるような低分子触媒反応は除く)

しかし、未経験者にはその苦労はわからないでしょう!実際、わたしも若くて魅力的だった頃は、「計算ってワンクリックで終わるんでしょ!全自動なんでしょ?」と思ってました(笑)。

 

今回は、知られざる計算化学の大変さの一部分を紹介したいと思います(注:オーサーシップの順番は、研究に対する貢献度であり、大変であるということとは直接因果関係がないです。)

書いていたら、止まらなくなり非常に長くなってしまったので、3回に分けてお届けします。

Part 1 計算と実験には共通点が多い!

Part 2 コントリビューションの決め方(予定)

Paer 3  計算または実験。なぜ二者択一?(予定)

計算と実験には共通点が多い!

計算化学というと、全自動化・洗練されていて、実験のように地道な検討、個人のテクニックなど必要ないと思われがち。しかし、実際には違います。個々人のセンスが要求されますし、毎日ひたすら検討に次ぐ検討!です。しかし、こういった苦労は論文には載りません。一発で正解が出るわけではなく、試行錯誤しながら求めていく点において、苦労があり、実験化学と通じるものがあります。

 

有機化学者にとって最も身近なのは、反応機構解析の論文だと思いますので、まずは反応経路探索を例としてお話しましょう。

 

  1. 初期検討

反応経路探索では、遷移状態の構造を求め作業が研究全体の中で大きな比重を占めています。反応の種類にもよりますが、そんなに簡単には求まりません。逆に、遷移状態さえ求められてしまえば、あとはルーティーンワークとも言えます。

まずは軽い計算レベルで初期検討します。 遷移状態構造を最適化する際の初期構造を探索します。いきなり重たい基底関数での DFT 計算や MP2 なんて使いません。

反応点周りの結合長や結合角度などを少しずつ変えて、どのようにエネルギーが変化するかなどを調べます。この時点で、この後の計算がうまくいくかどうかほぼ分かります。この段階では、正確性よりも手数の多さ、計算量の数・スピードが問われるため、パラメーターも若干落とします。RMSDP の threshold 下げるとかね

このような初期検討は、有機合成でも同様に行われると思います。

 

  1. 計算レベルを上げる

続いて、高い計算レベルで計算します。ここからは、一つ一つの計算に時間がかかります。原子数や計算レベルにもよりますが、一つの計算に半日以上〜数日かかりますので、同時に数十個計算を投げます。この点は、反応開発などに近いと思います。半日以上かかる反応を同時に何個も仕込むような感覚です。でも計算レベルを変えた瞬間に遷移状態が求められなくなったり、構造が変わってしまったりすることも時々おこるので、その場合また初期構造を求め直さなくてはなりません。実験でもスケールアップして反応がおかしくなってしまうことってありますよね。

 

  1. 確認作業

遷移状態構造が最適化できたら、それが望みの遷移状態構造であるかを確認しなければなりません。まずは振動計算で虚振動があることを確認した後、IRC計算を行います(確かめ算のようなものだと思ってください) 。有機合成でも合成・精製できたっ!と思ったらMSや各種スペクトルで構造を確認しますよね。

その後、一点計算を行ったり、DFT計算の場合様々な汎関数を検討したりします。ちなみにアルゴリズムによっては、 IRC 計算には 1 週間以上かかる場合もあります。

 

  1. さらなる検討

また、紙の上で構造式を書く場合、2次元になってしまいますが、実際の3次元の構造は非常に複雑です。微妙なコンフォメーションの違い、置換基の角度、向きなどちょっとした構造の違いにより、全く異なる計算結果が出てしまいます。

例えば、Kendal HoukDiels-Alder 反応に対し、200 通りもの経路を計算することもあると講演で述べてました。[1] どんだけ Diels Alder 好きなんだよって感じですよね。

 

考えられる構造全てを計算する、これが計算化学の強みでもあり、苦労する点でもあります。

 

また、よくある間違いとしては、初心者が予想中間体と遷移状態をそれぞれ構造最適化して求めて、IRC計算をせずにエネルギーを求めることです。
また、構造最適化 ひとつとっても指定すべきキーワードはたくさんありますし、エネルギー計算でも指定するキーワードにより値が変わってきます。どのようなキーワードを指定して計算をしているかは論文には記載する必要がないため、正直言って、計算をきちんと学んだことのない人の計算結果は疑わしいです。

以上述べた、反応遷移状態の求め方は、すぐできる量子化学計算ビギナーズマニュアルに書いてあります。

 

大きな誤解

計算投げたら、あとは計算機がやってくれるから暇」というのは、間違っています。計算結果(過程)は、数時間おきに確認し、変な方向に進んでいないか確かめる必要があります。これは、有機合成で TLC 取るような感覚ですね。

また、初期検討を行う場合などは、一時間数十個のジョブを投げます。さらに、その後数十個のlog ファイルを解析しなくてはなりません。

実験で同時に何個も反応を仕込むのと同じように、同時に何個も計算を投げるため、だんだんと頭が混乱してきます。要領の良い人は、常時数十個の計算を投げ、それらを解析し、まとめていきます。一方で、要領の悪い人は、2つジョブを投げて2つ解析してといった感じで、解析中はジョブ回ってないィィィイ!!みたいな感じになつため進みが遅いです。

これも実験と同じですよね。

また、解析が終わって時間に余裕があるときは、自分でプログラム書いたりしています。市販のソフトウェア(Gaussian)使って「結果出ましたぁ!」なんて言っているのは、計算化学者ではないです。自分の使っているソフトウェアの内容をあまり理解せずに、ブラックボックス的に使用し、出てきた結果を鵜呑みにするなんてのは、あまりよくないですよね。

計算したことない人は、「ワンクリックで終わるんでしょ!」と言いますが、「ワンクリックで解析できているようにその解析プログラム書いたのわたしですよぉぉぉ!!」っていつも思っています。

実験だって、「試薬混ぜて一定温度で放置するだけだから、暇だろ!」なんて言われたら気分よくないですよね?

 

参考文献

  1.  Kersey Black, Peng Liu, Lai Xu, Charles Doubleday, Kendall N. Houk, Proc. Natl. Acad. Sci. 2012, 109, 12860–12865. doi: 10.1073/pnas.1209316109
ゼロ

ゼロ

投稿者の記事一覧

女の子。研究所勤務。趣味は読書とハイキング ♪
ハンドルネームは村上龍の「愛と幻想のファシズム」の登場人物にちなんでま〜す。5 分後の世界、ヒュウガ・ウイルスも好き!

関連記事

  1. DNAが絡まないためのループ
  2. 親水性ひも状分子を疎水性空間に取り込むナノカプセル
  3. J-STAGE新デザイン評価版公開 ― フィードバックを送ろう
  4. 有機分子触媒ーChemical Times特集より
  5. ケムステイブニングミキサー2018へ参加しよう!
  6. 医薬品への新しい合成ルートの開拓 〜協働的な触媒作用を活用〜
  7. その構造、使って大丈夫ですか? 〜創薬におけるアブナいヤツら〜
  8. イミンを求核剤として反応させる触媒反応

コメント、感想はこちらへ

注目情報

ピックアップ記事

  1. 投票!2019年ノーベル化学賞は誰の手に!?
  2. リチャード・シュロック Richard R. Schrock
  3. 液体ガラスのフシギ
  4. フォルスター・デッカー アミン合成 Forster-Decker Amine Synthesis
  5. ODOOSをリニューアル!
  6. 新たな特殊ペプチド合成を切り拓く「コドンボックスの人工分割」
  7. 【日産化学】新卒採用情報(2022卒)
  8. 2009年8月人気化学書籍ランキング
  9. トリス(2,4-ペンタンジオナト)鉄(III):Tris(2,4-pentanedionato)iron(III)
  10. 「イスラム国(ISIS)」と同名の製薬会社→社名変更へ

関連商品

ケムステYoutube

ケムステSlack

注目情報

注目情報

最新記事

転職を成功させる「人たらし」から学ぶ3つのポイント

転職活動を始めた場合、まずは自身が希望する職種、勤務地、年収などの条件を元にインターネットで求人を検…

mRNAワクチン(メッセンジャーRNAワクチン)

病原体のタンパクをコードしたmRNAをベースとしたワクチン。従来のワクチンは、弱毒化・不活化した病原…

第139回―「超高速レーザを用いる光化学機構の解明」Greg Scholes教授

第139回の海外化学者インタビューはグレッグ・ショールズ教授です。トロント大学化学科(訳注:現在はプ…

分子の対称性が高いってどういうこと ?【化学者だって数学するっつーの!: 対称操作】

群論を学んでいない人でも「ある分子の対称性が高い」と直感的に言うことはできるかと思います。しかし分子…

非古典的カルボカチオンを手懐ける

キラルなブレンステッド酸触媒による非古典的カルボカチオンのエナンチオ選択的反応が開発された。低分子触…

CEMS Topical Meeting Online 機能性材料の励起状態化学

1月28日に毎年行われている理研の無料シンポジウムが開催されるようです。事前参加登録が必要なので興味…

カルボン酸に気をつけろ! グルクロン酸抱合の驚異

 カルボン酸は、カルボキシ基 (–COOH) を有する有機化合物の一群です。カルボン…

第138回―「不斉反応の速度論研究からホモキラリティの起源に挑む」Donna Blackmond教授

第138回の海外化学者インタビューはドナ・ブラックモンド教授です。2009年12月現在、インペリアル…

Chem-Station Twitter

PAGE TOP