W W FFFFFF PPPPP PPPPP PPPPP W W F P P P P P P W W F P P P P P P W W W FFF PPPPP PPPPP PPPPP W W W F P P P W W W F P P P WW WW F P P P Version 2011.11.29 ------------------------------------------------------------------------------ HMO und Pi-SCF-CI-Berechnungen fuer Singulett-Zustaende. Entwickelt aus dem PPP-CSVCH-Programm von E.Pfeil, A.Mehlhorn und J.Fabian, TU Dresden. ------------------------------------------------------------------------------ Referenzzitat: Programm WFPPP, PPP-Pi-SCF-Singulett-CI, Version 2011.11.29 Juergen Fabian, Technische Universitaet Dresden und Hans-Ulrich Wagner, Ludwig-Maximilians-Universitaet Muenchen. ------------------------------------------------------------------------------ FORTRAN-Versionen: DOS/Lahey, SGI, IBM, Linux Maximaldimensionierung: Anzahl Atomzentren ND2 = 200 Besetzte CI-MOs U = 15 Unbesetzte CI-MOs V = 15 Configurationen ND3 = 226 Kuroda Vector NKU = 15 ------------------------------------------------------------------------------ * = Ausgabe der Eingabezeilen = Eingabe-Datei ============================================================================== * MOLECULE CHARGE=1 * GENERATED BY GABEDIT 2.4.1 * TIME = MON NOV 28 15:46:11 2011 * 16 15 0 0 0 0999 V2000 * -0.1735 -3.6386 -1.4486 C 0 0 0 0 0 0 0 0 0 0 0 0 * 0.3465 -2.4406 -0.9116 C 0 0 0 0 0 0 0 0 0 0 0 0 * -0.4027 -1.2591 -0.7201 C 0 0 0 0 0 0 0 0 0 0 0 0 * 1.3995 -2.4253 -0.6290 H 0 0 0 0 0 0 0 0 0 0 0 0 * 0.0924 -0.0838 -0.2136 C 0 0 0 0 0 0 0 0 0 0 0 0 * 0.5453 -4.7886 -1.6572 C 0 0 0 0 0 0 0 0 0 0 0 0 * 0.0079 -6.0039 -2.0872 N 0 0 0 0 0 0 0 0 0 0 0 0 * -0.6702 1.0492 0.0771 N 0 0 0 0 0 0 0 0 0 0 0 0 * 1.1382 -0.0073 0.0747 H 0 0 0 0 0 0 0 0 0 0 0 0 * 1.6028 -4.8306 -1.4064 H 0 0 0 0 0 0 0 0 0 0 0 0 * 0.6320 -6.5716 -2.6469 H 0 0 0 0 0 0 0 0 0 0 0 0 * -0.9060 -5.9312 -2.5185 H 0 0 0 0 0 0 0 0 0 0 0 0 * -0.1857 1.9306 -0.0411 H 0 0 0 0 0 0 0 0 0 0 0 0 * -1.5791 1.0696 -0.3695 H 0 0 0 0 0 0 0 0 0 0 0 0 * -1.2328 -3.6500 -1.7164 H 0 0 0 0 0 0 0 0 0 0 0 0 * -1.4600 -1.2865 -0.9944 H 0 0 0 0 0 0 0 0 0 0 0 0 * 1 2 2 0 0 0 0 * 1 6 1 0 0 0 0 * 1 15 1 0 0 0 0 * 2 3 1 0 0 0 0 * 2 4 1 0 0 0 0 * 3 5 2 0 0 0 0 * 3 16 1 0 0 0 0 * 5 8 1 0 0 0 0 * 5 9 1 0 0 0 0 * 6 7 1 0 0 0 0 * 6 10 1 0 0 0 0 * 7 11 1 0 0 0 0 * 7 12 1 0 0 0 0 * 8 13 1 0 0 0 0 * 8 14 1 0 0 0 0 * M END ============================================================================== Ende der Eingabezeilen. ------------------------------------------------------------------------------ Summenformel fuer das Gesamtsystem H 9 C 5 N 2 ------------------------------------------------------------------------------ ---- Selektion der Pi-Atome 7 Pi-Atome gefunden 6 Pi-Bindungen gefunden 1 Gesamt-Ladung eingegeben 8 Pi-Elektronen berechnet Beta-Variation bis zum Atom BET= 7 ------------------------------------------------------------------------------ Projektion des gefundenen pi-Systems ACHTUNG! Die Achsen koennen vertauscht sein! Projektion auf die Y - Z Ebene Skalierungsfaktor 1.00, 1 Angstrom = 2.50 cm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Y 7 4 . . . 3 . . 2 . . . 1 . 5 . . 6 . ------------------------------------------------------------------------------ MOLECULE CHARGE=1 ------------------------------------------------------------------------------ Start PPP-Programm WFPPP Version 2011.11.29 PA = 4 Quantenchemische Berechnung nach dem pi-SCF-Verfahren von Pople mit Elektronenrepulsion nach Mataga mit Beta-Variation Betaij= Betaij * exp(BETA*(1.397-(1.517-0.18 * Pij))) BETA = 2.30000 Beta-Variation HET = 0 nur fuer C=C-Bindungen GAMMA = 1.00000 Parameter Gamma Matrix Gammaij = 14.3986/( (28.7972/2*Gtyp)*NIL + GAMMA*Rij ) Weitere Steuerparameter BMX= 7 maximale Atomnummer bei Alpha-Beta-Variation NIL= 0.0 Nishimoto L in K=(0.33*L+0.48), NIL=0 -> K=1 GAM(K) = 14.3986/( (28.7972/GIJ)*K + GAMMA*R ) U = 4 besetzte MO'S in der Konfigurations-Wechselwirkung V = 3 unbesetzte MO'S in der Konfigurations-Wechselwirkung UE = 4 zu druckende Uebergaenge IS = 4 zu druckende Ladungsverteilungen im angeregten Singulettzustand IT = 99 Iterationen maximal BA = 0 Ausgabe fuer BABA-Konfigurations-Analyse KU = 0 KURODA-Analyse MIM= 0 MIM-Analyse AMA= 0 Aufteilungsmatrix (z.Zt. nicht installiert) KDR= 0 STOP nach Hueckel-Matrix-Ausdruck) A = 0 Druck Zwischenergebnisse und Rechenkontrollen D = 0 Druck Elektronenwechselwirkungsintegrale Gamma E = 0 Druck Reaktivitaets-Indices F = 0 Druck Dipolmomente G = 0 Druck Hartree-Fock-, Dichte- und FCS-Matrix H = 0 Druck Konfigurations-WW-Matrix und Ladungsverteilungen 0.000010 Iterationsschranke 11.420000 SCF-Alfa in eV fuer Standard C 2.318000 SCF-Beta in eV fuer Standard C=C Atomparameter Nr. KLZ X Y Z At Typ El. IP(eV) EA(eV) HMO-alfa 1 6 -0.17350 -3.63860 -1.44860 C pi1 1. 11.42000 0.58000 0.00000 2 6 0.34650 -2.44060 -0.91160 C pi1 1. 11.42000 0.58000 0.00000 3 6 -0.40270 -1.25910 -0.72010 C pi1 1. 11.42000 0.58000 0.00000 4 6 0.09240 -0.08380 -0.21360 C pi1 1. 11.42000 0.58000 0.00000 5 6 0.54530 -4.78860 -1.65720 C pi1 1. 11.42000 0.58000 0.00000 6 7 0.00790 -6.00390 -2.08720 N pi2 2. 21.22000 8.24000 1.37000 7 7 -0.67020 1.04920 0.07710 N pi2 2. 21.22000 8.24000 1.37000 Gefundene Bindungen mu-nu und Bindungsparameter mu - nu --- Bindung --- DistTyp BetaSCF Gamma HMO-beta 1 - 2 C pi1 - C pi1 1.40000 2.31800 10.84000 1.00000 1 - 5 C pi1 - C pi1 1.40000 2.31800 10.84000 1.00000 2 - 3 C pi1 - C pi1 1.40000 2.31800 10.84000 1.00000 3 - 4 C pi1 - C pi1 1.40000 2.31800 10.84000 1.00000 4 - 7 C pi1 - N pi2 1.40000 2.31800 11.91000 0.89000 5 - 6 C pi1 - N pi2 1.40000 2.31800 11.91000 0.89000 Abstandsmatrix R aus xyz-Koordinaten 1 2 3 4 5 6 7 1 0.0000 2 1.4121 0.0000 3 2.4991 1.4121 0.0000 4 3.7726 2.4711 1.3722 0.0000 5 1.3721 2.4715 3.7728 4.9421 0.0000 6 2.4567 3.7675 4.9549 6.2101 1.3967 0.0000 7 4.9548 3.7669 2.4567 1.3963 6.2101 7.4088 0.0000 Alpha-Korrekturen und Beta-Werte vor den SCF-Iterationen (in eV) 1 2 3 4 5 6 7 1 0.0000 2 2.3180 0.0000 3 0.0000 2.3180 0.0000 4 0.0000 0.0000 2.3180 0.0000 5 2.3180 0.0000 0.0000 0.0000 0.0000 6 0.0000 0.0000 0.0000 0.0000 2.0630 0.0000 7 0.0000 0.0000 0.0000 2.0630 0.0000 0.0000 0.0000 Atomare Wechselwirkungsintegrale Gamma 1 2 3 4 5 6 7 1 10.8400 2 5.2543 10.8400 3 3.7620 5.2543 10.8400 4 2.8228 3.7897 5.3318 10.8400 5 5.3320 3.7893 2.8226 2.2963 10.8400 6 3.9280 2.8934 2.3360 1.9408 5.5260 12.9800 7 2.3360 2.8937 3.9280 5.5267 1.9408 1.6904 12.9800 Beginn SCF-Rechnung mit Elektronen-Wechselwirkung ========================== 1.Iteration, maximale Dichteaenderung= 0.010639 2.Iteration, maximale Dichteaenderung= 0.001308 3.Iteration, maximale Dichteaenderung= 0.000520 4.Iteration, maximale Dichteaenderung= 0.000203 5.Iteration, maximale Dichteaenderung= 0.000076 6.Iteration, maximale Dichteaenderung= 0.000029 7.Iteration, maximale Dichteaenderung= 0.000010 Werte nach der SCF-Konvergenz mit 8 Iterationen ========================== Alpha-Korrekturen und Beta-Werte nach den SCF-Iterationen (in eV) 1 2 3 4 5 6 7 1 0.0000 2 2.2820 0.0000 3 0.0000 2.2821 0.0000 4 0.0000 0.0000 2.2986 0.0000 5 2.2987 0.0000 0.0000 0.0000 0.0000 6 0.0000 0.0000 0.0000 0.0000 2.0630 0.0000 7 0.0000 0.0000 0.0000 2.0630 0.0000 0.0000 0.0000 Atomare Wechselwirkungsintegrale Gamma 1 2 3 4 5 6 7 1 10.8400 2 5.2543 10.8400 3 3.7620 5.2543 10.8400 4 2.8228 3.7897 5.3318 10.8400 5 5.3320 3.7893 2.8226 2.2963 10.8400 6 3.9280 2.8934 2.3360 1.9408 5.5260 12.9800 7 2.3360 2.8937 3.9280 5.5267 1.9408 1.6904 12.9800 HARTREE-FOCK-Matrix 1 2 3 4 5 6 7 1 -9.2907 2 -3.9342 -8.8692 3 0.2505 -3.9343 -9.2907 4 0.2899 -0.3043 -4.0219 -9.6791 5 -4.0221 -0.3042 0.2899 0.1385 -9.6786 6 0.5620 0.3714 -0.1939 -0.1524 -3.9440 -12.3140 7 -0.1939 0.3714 0.5620 -3.9443 -0.1524 0.1405 -12.3143 SCF-Eigenwerte (eV) und zugehoerige SCF-Eigenvektoren HOMO = 4 LUMO = 5 1 2 3 4 5 6 7 -16.6021 -15.9202 -14.1361 -11.8109 -6.5873 -4.0821 -2.2979 1 0.3965 -0.2640 -0.3037 -0.4965 0.0819 -0.4287 -0.4937 2 0.4059 -0.0002 -0.5053 0.0000 -0.5273 -0.0001 0.5495 3 0.3967 0.2637 -0.3038 0.4965 0.0819 0.4288 -0.4937 4 0.4063 0.4704 0.1673 0.1789 0.4655 -0.4966 0.3005 5 0.4061 -0.4706 0.1675 -0.1789 0.4654 0.4966 0.3006 6 0.3086 -0.4572 0.5021 0.4706 -0.3710 -0.2636 -0.1224 7 0.3089 0.4572 0.5020 -0.4706 -0.3710 0.2636 -0.1224 Summe ( MO-Energien * Anzahl Elektronen ) -116.9386 (eV) Dichte-Matrix (Elektronendichten und Bindungsordnungen) 1 2 3 4 5 6 7 1 1.1315 2 0.6289 0.8401 3 -0.1332 0.6289 1.1314 4 -0.2054 0.1606 0.6464 0.8928 5 0.6464 0.1606 -0.2054 -0.1206 0.8927 6 -0.2862 -0.2567 0.1660 0.1571 0.6808 1.5558 7 0.1660 -0.2567 -0.2862 0.6808 0.1571 -0.1662 1.5558 ------------------------------------------------------------------------------ Molekulare Elektronen-Wechselwirkungs-Integrale G(k,l,j,i) MEW N= 7 M=HOMO= 4 LUMO= 5 UU= 4 VV= 3 J G( 4, 4, 5, 5)= 4.3694 eV K G( 4, 5, 4, 5)= 1.1177 eV MOLECULE CHARGE=1 ========================================================================= WFPPP Plot log(Oscill) 10.0 . . . . . . . * . * . * 1.00 * . * . * . * . * . * . * . * . * . * * 0.10 * * . * * . * * * . * * * . * * * * . * * * * . * * * * . * * * * * . * * * * * . * * * * * 0.01 * * * * * . * * * * * . * * * * * . * * * * * . * * * * * * . * * * * * * . * * * * * * . * * * * * * . 7 5 4 3 2 1 . ---- cm-1 --50000----------40000-----------30000----------20000-----------10000 Singulett CI mit 12 angeregten Configurationen Liste fuer UE= 4 angeregte Singulett-Zustaende Configurationen sind nur dann angegeben, wenn der Anteil groesser als 1.0% ist Zustandsenergie relativ Grundzustand Oszillatorst. Uebergangsmomente # E(eV) (cm-1) (nm) f X Y Z gesamt Conf. i -> j Ant.% 1 2.9730 23986 416 1.692 -0.233 2.427 0.745 2.550 2 4 - 5 96.9% 3.0896 -0.234 2.432 0.746 2.555 6 3 - 6 2.1% 6.9060 -0.031 0.329 0.101 0.345 2 5.0096 40417 247 0.059 -0.356 -0.007 -0.088 0.366 3 4 - 6 58.9% 5.1510 -0.137 0.007 -0.065 0.152 5 3 - 5 40.4% 5.1974 0.396 0.021 0.056 0.400 3 5.2377 42257 236 0.019 0.202 0.019 0.000 0.203 3 4 - 6 39.4% 5.1510 -0.137 0.007 -0.065 0.152 5 3 - 5 57.1% 5.1974 0.396 0.021 0.056 0.400 11 1 - 5 2.3% 6.9606 0.146 0.013 0.002 0.147 4 6.3811 51483 194 0.003 -0.007 0.073 0.023 0.077 4 4 - 7 90.9% 6.4180 -0.010 0.107 0.034 0.113 8 2 - 5 8.0% 6.6041 0.015 -0.162 -0.050 0.170 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Grundzustand Singulettzustand Nr. 0 mu Ladung mu - nu mu - nu mu - nu 1 -0.131 1 - 2 0.629 1 - 5 0.646 2 0.160 2 - 3 0.629 3 -0.131 3 - 4 0.646 4 0.107 4 - 7 0.681 5 0.107 5 - 6 0.681 6 0.444 7 0.444 Dipolmoment fuer Ionen bezogen auf den Schwerpunkt der pi-Kernladungen. Dipolmoment Komponenten X= -0.035 Y= 0.059 Z= -0.204 Gesamt= 0.215 Projektion auf die Y - Z Ebene Skalierungsfaktor 1.00, 1 Angstrom = 2.50 cm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Y 7 4 . . . 3 . . 2 . . . 1 . 5 . . 6 . ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Angeregter Singulettzustand Nr. 1 mu Ladung mu - nu mu - nu mu - nu 1 0.118 1 - 2 0.571 1 - 5 0.576 2 -0.124 2 - 3 0.571 3 0.118 3 - 4 0.576 4 -0.054 4 - 7 0.605 5 -0.054 5 - 6 0.605 6 0.498 7 0.498 Dipolmoment fuer Ionen bezogen auf den Schwerpunkt der pi-Kernladungen. Dipolmoment Komponenten X= 1.793 Y= 0.105 Z= 0.222 Gesamt= 1.810 Projektion auf die Y - Z Ebene Skalierungsfaktor 1.00, 1 Angstrom = 2.50 cm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Y 7 4 . . . 3 . . 2 . . . 1 . 5 . . 6 . Aenderungen Grundzustand -> Angeregter Zustand Nr. 1 der pi-Ladungen und Bindungsordnungen mu mu - nu mu - nu mu - nu 1 0.250 1 - 2 -0.058 1 - 5 -0.070 2 -0.284 2 - 3 -0.058 3 0.250 3 - 4 -0.070 4 -0.161 4 - 7 -0.076 5 -0.161 5 - 6 -0.076 6 0.053 7 0.053 Summe der Betraege der E-Dichte-Aenderungen CT= 1.213 Summe der Betraege der Bindungs-Aenderungen CB= 0.406 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Angeregter Singulettzustand Nr. 2 mu Ladung mu - nu mu - nu mu - nu 1 -0.052 1 - 2 0.529 1 - 5 0.504 2 0.141 2 - 3 0.529 3 -0.052 3 - 4 0.504 4 -0.087 4 - 7 0.563 5 -0.087 5 - 6 0.563 6 0.568 7 0.568 Dipolmoment fuer Ionen bezogen auf den Schwerpunkt der pi-Kernladungen. Dipolmoment Komponenten X= 1.205 Y= 0.120 Z= -0.005 Gesamt= 1.211 Projektion auf die Y - Z Ebene Skalierungsfaktor 1.00, 1 Angstrom = 2.50 cm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Y 7 4 . . . 3 . . 2 . . . 1 . 5 . . 6 . Aenderungen Grundzustand -> Angeregter Zustand Nr. 2 der pi-Ladungen und Bindungsordnungen mu mu - nu mu - nu mu - nu 1 0.079 1 - 2 -0.100 1 - 5 -0.143 2 -0.018 2 - 3 -0.100 3 0.079 3 - 4 -0.143 4 -0.194 4 - 7 -0.118 5 -0.194 5 - 6 -0.118 6 0.124 7 0.124 Summe der Betraege der E-Dichte-Aenderungen CT= 0.814 Summe der Betraege der Bindungs-Aenderungen CB= 0.721 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Angeregter Singulettzustand Nr. 3 mu Ladung mu - nu mu - nu mu - nu 1 -0.031 1 - 2 0.493 1 - 5 0.561 2 0.231 2 - 3 0.493 3 -0.031 3 - 4 0.561 4 -0.127 4 - 7 0.520 5 -0.127 5 - 6 0.520 6 0.543 7 0.543 Dipolmoment fuer Ionen bezogen auf den Schwerpunkt der pi-Kernladungen. Dipolmoment Komponenten X= 1.156 Y= 0.112 Z= -0.003 Gesamt= 1.161 Projektion auf die Y - Z Ebene Skalierungsfaktor 1.00, 1 Angstrom = 2.50 cm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Y 7 4 . . . 3 . . 2 . . . 1 . 5 . . 6 . Aenderungen Grundzustand -> Angeregter Zustand Nr. 3 der pi-Ladungen und Bindungsordnungen mu mu - nu mu - nu mu - nu 1 0.100 1 - 2 -0.136 1 - 5 -0.085 2 0.071 2 - 3 -0.136 3 0.100 3 - 4 -0.085 4 -0.234 4 - 7 -0.161 5 -0.234 5 - 6 -0.161 6 0.099 7 0.099 Summe der Betraege der E-Dichte-Aenderungen CT= 0.937 Summe der Betraege der Bindungs-Aenderungen CB= 0.764 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Angeregter Singulettzustand Nr. 4 mu Ladung mu - nu mu - nu mu - nu 1 -0.124 1 - 2 0.377 1 - 5 0.421 2 -0.138 2 - 3 0.377 3 -0.124 3 - 4 0.421 4 0.055 4 - 7 0.691 5 0.055 5 - 6 0.691 6 0.638 7 0.638 Dipolmoment fuer Ionen bezogen auf den Schwerpunkt der pi-Kernladungen. Dipolmoment Komponenten X= 1.257 Y= 0.131 Z= -0.030 Gesamt= 1.264 Projektion auf die Y - Z Ebene Skalierungsfaktor 1.00, 1 Angstrom = 2.50 cm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Y 7 4 . . . 3 . . 2 . . . 1 . 5 . . 6 . Aenderungen Grundzustand -> Angeregter Zustand Nr. 4 der pi-Ladungen und Bindungsordnungen mu mu - nu mu - nu mu - nu 1 0.007 1 - 2 -0.252 1 - 5 -0.225 2 -0.298 2 - 3 -0.252 3 0.007 3 - 4 -0.225 4 -0.052 4 - 7 0.010 5 -0.052 5 - 6 0.010 6 0.194 7 0.194 Summe der Betraege der E-Dichte-Aenderungen CT= 0.804 Summe der Betraege der Bindungs-Aenderungen CB= 0.975 ------------------------------------------------------------------------------ MOLECULE CHARGE=1 ENDE WFPPP Version 2011.11.29