1 / 20

Coupled Clusters (СС) и Quadratic CI (QCI) Методи

Coupled Clusters (СС) и Quadratic CI (QCI) Методи. Малонов алдехид. Аминов радикал. Дублет. Въглероден оксид. Тестови молекули. Синглети. Кога смятате с Coupled Clusters?. Винаги ... ако имате персонален суперкомпютър!.

luyu
Download Presentation

Coupled Clusters (СС) и Quadratic CI (QCI) Методи

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. Coupled Clusters (СС) и Quadratic CI (QCI) Методи

  2. Малонов алдехид Аминов радикал Дублет Въглероден оксид Тестови молекули Синглети

  3. Кога смятате с Coupled Clusters? Винаги ... ако имате персонален суперкомпютър! Получавате изключително точна енергия и вълнова функция с високо качество При неясна изчислителна стратегия и липса на експериментални данни се прави “юстиране” по СС-изчисления QCI е с почти същото качество

  4. Формулите ... CC/SD(TQ)/

  5. #p ccsd(t)/d95 units=au test Singles и doubles се включват, а приносът от triples просто се добавя към енергията  #p bd(t,e4t)/d95 units=au test Пресмятане с Brueckner’s doubles, приносът от triples просто се добавя към енергията … ... пресмятат се и МР4triples Какво въвеждате ...

  6. #P uccsd(t,e4t)/gen test tran=iabc Неограниченият вариант ...  #p uqcisd/sto-3g guess=alter force test ... и на Quadratic CI Какво въвеждате ... …………………………………… Z-матрица a842=116.18239988 18 20

  7. ... и какво получавате  ccsd(t)/d95 SCF Done: E(RHF) = -112.394512472 A.U. after 16cycles Convg = 0.3799D-08 -V/T = 2.0055 ................. Range of M.O.s used for correlation: 3 18 ................. Compute canonical integrals, LenV= 2051653. ................. E2= -0.2307371915D+00 EUMP2= -0.11262524966369D+03 ................. CCSD(T) ======= ................. T4(CCSD)= -0.99821730D-01 T5(CCSD)= 0.59006510D-01 CCSD(T)= -0.11271993664D+03 Job cpu time: 0 days 0 hours 0 minutes 7.0 seconds

  8. ... и какво получавате  bd(t,e4t)/d95 SCF Done: E(RHF) = -112.394512472 A.U. after 16cycles Convg = 0.3799D-08 -V/T = 2.0055 ................. Range of M.O.s used for correlation: 3 18 ................. E2= -0.2307371915D+00 EUMP2= -0.11262524966369D+03 ................. Brueckner Doubles with Triples (BD(T)) ====================================== ................. T4(BD)= -0.49017607D-01 BD(T)= -0.11270565566D+03 Job cpu time: 0 days 0 hours 0 minutes 18.0 seconds

  9. ... и какво получавате  uccsd(t,e4t)/gen General basis read from cards: (5D, 7F) ................... SCF Done: E(UHF) = -55.5771823051 A.U. after 12 cycles Convg = 0.4772D-08 -V/T = 2.0007 S**2 = 0.7579 Annihilation of the first spin contaminant: S**2 before annihilation 0.7579, after 0.7500 .................... E2= -0.1432673103D+00 EUMP2= -0.55720449615407D+02 .................... CCSD(T) ======= ................. T4(CCSD)= -0.27104732D-02 T5(CCSD)= 0.45953846D-04 CCSD(T)= -0.55742073339D+02 Job cpu time: 0 days 0 hours 0 minutes 16.0 seconds

  10. ... и какво получавате  uqcisd/sto-3g Pairs of Alpha orbitals switched: 18 20 No Beta orbitals switched. ................... SCF Done: E(UHF) = -262.011304022 A.U. after 15 cycles Convg = 0.6283D-08 -V/T = 2.0058 S**2 = 1.2113 ................... E2= -0.2567064288D+00 EUMP2= -0.26226801045131D+03 ................... Quadratic Configuration Interaction (QCISD) =========================================== ................... DE(Corr)= -0.32066872E(CORR)= -262.33197274 Delta=-2.00D-09 Job cpu time: 0 days 0 hours 0 minutes 39.0 seconds

  11. ccsd(t)/d95  E(RHF) = -112.39451247 a.u. E(MP2) = -112.62524966 a.u. E(CCSD)= -112.67912142 a.u. E(CCSD(T))= -112.71993664 a.u. DE=0.231 a.u.=144.79 kcal/mol  DE=0.054 a.u.= 33.81kcal/mol  DE=0.041 a.u.= 25.61kcal/mol bd(t,e4t)/d95  E(RHF) = -112.39451247 a.u. E(MP2) = -112.62524966 a.u. E(BD)= -112.65663805 a.u. E(BD(T))= -112.70565566 a.u. DE=0.231 a.u.=144.79 kcal/mol  DE=0.031 a.u.= 19.70kcal/mol  DE=0.049 a.u.= 30.76 kcal/mol Количествена оценка

  12. uccsd(t,e4t)/gen  E(UHF) = -55.57718231 a.u. E(MP2) = -55.72044962 a.u. E(CCSD)= -55.73940891 a.u. E(CCSD(T))= -55.74207334 a.u. DE=0.143 a.u.= 89.90kcal/mol  DE=0.019 a.u.= 11.90 kcal/mol  DE=0.003 a.u.= 1.67 kcal/mol Количествена оценка

  13. uqcisd/sto-3g  E(UHF) = -262.01130402 a.u. E(MP2) = -262.26801045 a.u. E(QCI)= -262.33197274 a.u. DE=0.257 a.u.=161.09kcal/mol  DE=0.064 a.u.= 40.14 kcal/mol uccsd(t)/sto-3g  E(RHF) = -262.01130402 a.u. E(MP2) = -262.26801045 a.u. E(CCSD)= -262.32751820 a.u. E(CCSD(T))= -262.34577626 a.u. DE=0.257 a.u.=161.09kcal/mol  DE=0.060 a.u.= 37.34 kcal/mol  DE=0.018 a.u.= 11.46 kcal/mol Количествена оценка

  14. Енергия на връзката в Н—F при базис 6-31++G(3df,3pd) Енергия на ковалентни връзки

  15. Геометрия на озонпри базис 6-31G* Геометрия Геометрията се счита за ‘точна’ при отклонения във връзките < 0.010.02 Å и във валентните и торзионните ъгли < 120 спрямо експериментална структура

  16. Зависимост от размера на задачата Ресурсите ... N – брой базисни ф-ции O – брой заети МО V – брой виртуални МО

  17. Ресурсите ... Зависимост от метода и базиса

  18. CC2 – комбинация между CC и MP2 Много по-бърз – позволява пресмятане на системи със стотици атоми (фулерени,порфирини)!!!

  19. Как става? Resolution of the Identity (RI) Ускорява изчислението на интегралите многократно, без да нараства грешката

  20. Ефективност Абсолютна и средна грешка при използване на едноелектронен базис и RI: възбудено състояние с CC2/aug-cc-pVTZ

More Related