1 / 30

Anomalous non-additive dispersion interactions in systems of three one -dimensional wires

押山科研費 /10Mar'14/15+5min. ( 東大工 6). Anomalous non-additive dispersion interactions in systems of three one -dimensional wires. Ryo Maezono. rmaezono@mac.com. School of Information Science, JAIST, Ishikawa, Japan. 本科研費期間中の原著業績. DIB 分子結晶の量子拡散モンテカルロ法電子状態計算 /

temima
Download Presentation

Anomalous non-additive dispersion interactions in systems of three one -dimensional wires

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. 押山科研費/10Mar'14/15+5min (東大工6) Anomalous non-additive dispersion interactions in systems of three one-dimensional wires Ryo Maezono rmaezono@mac.com School of Information Science, JAIST, Ishikawa, Japan.

  2. 本科研費期間中の原著業績 DIB分子結晶の量子拡散モンテカルロ法電子状態計算/ K. Hongo and T. Iitaka, M.W, A.A, and R. Maezono, submitted to J. Chem. Theory Comput. (2014). 金属ナノワイヤ間の分散力相互作用と非加算性/ A.J. Misquitta, R. Maezono, N.D.D., A.J. Stone, R.J.N, Phys. Rev. B 89, 045140 (2014) 半導体二層膜の中密度域におけるバイエキシトン気体 R. Maezono, P.L. Rios, T. Ogawa, and R.J. Needs, Phys. Rev. Lett. 110, 216407 (2013). DNA塩基間引力の量子モンテカルロ法電子状態計算, K. Hongo N.T. Cuong, and R. Maezono. J. Chem. Theory Comput. 9, 1081 (2013). GPGPUに適した量子モンテカルロの新しい配位更新法/ Y. Uejima and R. Maezono. J. Comput. Chem. 34, 83 (2013). 化合物半導体構造相転移の量子モンテカルロ計算/ C.N.M. Ouma, M.Z.M., N.W.M., G.O.A., and R. Maezono. Phys. Rev. B 86 104115 (2012). チタン酸化物の密度汎関数計算/ M. Abbasnejad, M. R. M. and R. Maezono, Europhys. Lett. 97, 56003 (2012). チタン酸化物の量子モンテカルロ計算/ M. Abbasnejad, E.S., M.R.M., M.A., and R. Maezono, Appl. Phys. Lett. 100, 261902 (2012). クロム二量体結合の量子モンテカルロ計算/ K. Hongo and R. Maezono, Int. J. Quant. Chem. 112, 1243 (2012).

  3. vdW in Nano Wire Equilateral Triangle Geometry + + + - - - + + + - - -

  4. vdW in Nano Wire An interesting QMC challenge. ・Difficulty of conventional DFT XC functionals for Dispersion force ・Metallic Wire (Different Power law) → local polarization not well be defined... (Screening Length) ・Non-additivity Additive modeling ... valid only for local polarizations

  5. Previous Studies Between Insulating wires Metal wires/RPA J.F. Dobsonet.al., PRL96, 073201 (2006) Metal-Semiconductor/SAPT Huckel J. Spencer and A. Alavi, (thesis) SAPT（TB + Perturbation） A. Misquittaet.al., Phys. Rev. B, 82 , 075312, (2010). Metal/QMC N.D. Drummond et.al., Phys. Rev. Lett. 99, 247401 (2009). Much earlier works for 1-dim. - Coulson and Davies, Trans. Faraday Soc. (1952). - Longuet-Higgins and Salem, Proc. R. Soc. A. (1961).

  6. Modeling DMC (diffusion Monte Carlo) method ・Intra-Wire ; Anti-symetrized products of PW orbitals 1D HEG → No Nodal problem. Equilateral Triangle Geometry ・Inter-Wire ; treated as Distinguishable particles + + + - ・Odd number of electrons ; Real WF. - - + + + - - - N=55, 111 Bi-wire interaction Non-additive contribution ; Energy/wire @ Bi (Tri)-wire system

  7. Rs = 1.0

  8. Rs = 3.0 Diffusion Monte Carlo by RM (2012) α=2.541

  9. Rs = 10.0 09

  10. Kink in the exponent Huckel Note; correlation effect makes kink broader [Misquetta/SAPT, 2010] 22

  11. 2nd order perturbation inter-molecular interactions integral representation 18

  12. Longuet-Higgins Representation of dispersion energy H. C. Longuet-Higgins, Discuss. Faraday Soc. 40, 7 1965. non-local polarizability multi-pole expansion 19

  13. Exponent of decay Induced dipole-dipole int. Rank of tensor --> lowest possible contribution --> 20

  14. Inv.Sq. contribution Insulator (localized) ; ; localization length --> inv-sq. contribution vanishes ; integral region valid for multi-pole expansion 21 gets larger --> integral gets disappear

  15. Kink in the exponent Note; correlation effect makes kink broader [Misquetta/SAPT, 2010] 22

  16. Rs = 3.0 Diffusion Monte Carlo by RM (2012) α=2.541

  17. Exponents 2体相互作用 3体非加算性寄与 Long-ranged interactions ; α< 3

  18. Exponents of Decay Summing up interactions… We expect it not diverging. ; Num. of wires per unit area ；Energy per wire 、 → Finite → Diverge ! 15

  19. Exponents of Decay Possible mechanisms making it finite. two-body superposition non-additivity cancelled by and higher. １）Divergence of (c.f., Our equilateral triangle gives repulsive ) ２）Exponent falls off more rapidly at larger distances. (at which we cannot perform accurate calculations) ３）Relativistic retardation reducing the interaction at larger distances.

  20. Exponents 2体相互作用 3体非加算性寄与 Long-ranged interactions ; α< 3

  21. 今回の仕事 先行研究はいずれも、分極率評価に簡素な近似を適用したもの 分極率評価に第一原理計算を用いて何か新しい事が見つかるか？ 甲、 の範囲で、DFT摂動計算の枠組み(SAPT)で評価。 /A.J. Misquitta 乙、DMCで評価（こちらは摂動表式ではなく、第一原理で評価） /R. Maezono 甲は、ギャップのある系にしか適用出来ない。 乙はメタルに適用出来るが、は切り分けられない 乙は、一次元メタルに適用したが、甲は、の水素鎖に適用。

  22. 先行研究/非加算性寄与 ；n体相互作用から来る3次の寄与 を3つのダイポール間の相互作用で扱う。 [13] P. M. Axilrod and E. Teller, J. Chem. Phys. 11, 299 (1943). 昭和18年 [14] Y. Muto, Proc. Phys.-Math. Soc. Japan 17, 629 (1943). /非常に短距離な寄与を帰結 → "Axilrod-Teller-Muto contribution" 我々の今回の結論； 逆9乗の短距離ATM型ではなく、逆2.5乗程度の長距離型非加算寄与

  23. 非加算性寄与 摂動論による分散力寄与 Longuet-Higginsの一般化表式 /分散力に関わる分子内励起の「起こりやすさ」 ・重み ・分母/どの配置間での相互作用が生じているか 非加算性寄与

  24. 多重極展開 から、 → テンソルを込めた分極率 の上添え字； に関する非局所積分の事情を反映 下添え字；テンソルと組む事情、 すなわち異方性の事情

  25. 幾つかの極限状況 先行の結論を導くことが出来る 摂動による3体3次の分散力寄与 １）系のギャップが大きい場合 非局所積分の事情を反映するの上添え字はつぶれて、 → 非局所性弱くなる → Stogrynの結果 に帰着。 [D.E. Stogryn, Mol. Phys. 22, 81 (1971)] ギャップの大きな系のみで有効な結果

  26. 幾つかの極限状況 ２）双極子部分のみ考慮し、テンソルに関する異方性を方位平均化 → Axilrod–Teller–Mutoの形に帰着

  27. SAPT算定 分極率をコーンシャム摂動理論で評価 DFT計算はNWChem/PBE汎関数を用い、 KSPT計算はCamCASP/LDAの線形応答核 ダイマリゼーション・パラメタを含んだの水素鎖に適用。 基底系Sadlej-pVTZ 上記の分散エネルギーは「Dispersion」というコード実装で評価。

  28. SAPT解析からの結論 のSAPT/DFT摂動計算評価 ギャップを減らしていくと、ベキは上から3に近づく。 (c.f., DMC算定では、系が低密度となると冪は下から3に近づく） ギャップが小さくなると非加算寄与の増長 2体相互作用はファクター2倍だが、非加算性寄与は4倍の増長。

  29. SAPT解析からの結論 2体, 3体と展開していくと、どれもがr逆ベキの遅い収束、 且つ、交番級数的 →展開が良くないことを示唆。 代わる提案として、自己無撞着分極モデル - L. Silberstein, Phil. Mag. 33, 92 (1917). - J. Applequist et.al., J. Am. Chem. Soc. 94, 2952 (1972). (尚、DMCは全てを自己無撞着に扱う事に相当) 3体非加算性が正の寄与となる理由； +-に対して-+が誘起されて、更に其れにより+-が誘起されると、 元の+-に対して同符号・斥力となるため。

  30. まとめ 従来型の逆9乗と異なる、逆2.5乗程度 の非加算相互作用 Equilateral Triangle Geometry SAPTからも支持される。 + + + - - - + + + - - -

More Related