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 ============================================================================== * FULVEN * * * 12 12 0 0 0 0 0 0 0 0 0 0 * -0.7600 0.0000 -0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 * 0.1252 -1.1805 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 * 1.4065 -0.7382 -0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 * 1.4065 0.7382 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 * 0.1252 1.1805 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 * -2.1053 0.0000 -0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 * -0.2245 2.2081 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0 * -0.2246 -2.2081 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0 * -2.6730 0.9273 -0.0001 H 0 0 0 0 0 0 0 0 0 0 0 0 * -2.6730 -0.9273 0.0001 H 0 0 0 0 0 0 0 0 0 0 0 0 * 2.3034 -1.3516 -0.0001 H 0 0 0 0 0 0 0 0 0 0 0 0 * 2.3034 1.3516 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0 * 1 2 1 0 0 0 0 * 1 5 1 0 0 0 0 * 1 6 2 0 0 0 0 * 2 3 2 0 0 0 0 * 2 8 1 0 0 0 0 * 3 4 1 0 0 0 0 * 3 11 1 0 0 0 0 * 4 5 2 0 0 0 0 * 4 12 1 0 0 0 0 * 5 7 1 0 0 0 0 * 6 9 1 0 0 0 0 * 6 10 1 0 0 0 0 Keine weitere Zeile gefunden. Keine weitere Zeile gefunden. ============================================================================== Ende der Eingabezeilen. ------------------------------------------------------------------------------ Summenformel fuer das Gesamtsystem H 6 C 6 ------------------------------------------------------------------------------ ---- Selektion der Pi-Atome 6 Pi-Atome gefunden 6 Pi-Bindungen gefunden 0 Gesamt-Ladung eingegeben 6 Pi-Elektronen berechnet Beta-Variation bis zum Atom BET= 6 ------------------------------------------------------------------------------ Projektion des gefundenen pi-Systems ACHTUNG! Die Achsen koennen vertauscht sein! Projektion auf die X - Y Ebene Skalierungsfaktor 1.00, 1 Angstrom = 2.50 cm 5 . Y 4 . . . . 6 . . . . . 1 . . . . . . . X . . . . . . . 3 . 2 ------------------------------------------------------------------------------ FULVEN ------------------------------------------------------------------------------ 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= 6 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 = 3 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.76000 0.00000 0.00000 C pi1 1. 11.42000 0.58000 0.00000 2 6 0.12520 -1.18050 0.00000 C pi1 1. 11.42000 0.58000 0.00000 3 6 1.40650 -0.73820 0.00000 C pi1 1. 11.42000 0.58000 0.00000 4 6 1.40650 0.73820 0.00000 C pi1 1. 11.42000 0.58000 0.00000 5 6 0.12520 1.18050 0.00000 C pi1 1. 11.42000 0.58000 0.00000 6 6 -2.10530 0.00000 0.00000 C pi1 1. 11.42000 0.58000 0.00000 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 1 - 6 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 - 5 C pi1 - C pi1 1.40000 2.31800 10.84000 1.00000 Abstandsmatrix R aus xyz-Koordinaten 1 2 3 4 5 6 1 0.0000 2 1.4755 0.0000 3 2.2888 1.3555 0.0000 4 2.2888 2.3072 1.4764 0.0000 5 1.4755 2.3610 2.3072 1.3555 0.0000 6 1.3453 2.5236 3.5885 3.5885 2.5236 0.0000 Alpha-Korrekturen und Beta-Werte vor den SCF-Iterationen (in eV) 1 2 3 4 5 6 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 2.3180 0.0000 6 2.3180 0.0000 0.0000 0.0000 0.0000 0.0000 Atomare Wechselwirkungsintegrale Gamma 1 2 3 4 5 6 1 10.8400 2 5.1354 10.8400 3 3.9807 5.3651 10.8400 4 3.9807 3.9606 5.1338 10.8400 5 5.1354 3.9028 3.9606 5.3651 10.8400 6 5.3855 3.7380 2.9284 2.9284 3.7380 10.8400 Beginn SCF-Rechnung mit Elektronen-Wechselwirkung ========================== 1.Iteration, maximale Dichteaenderung= 0.071613 2.Iteration, maximale Dichteaenderung= 0.004755 3.Iteration, maximale Dichteaenderung= 0.002170 4.Iteration, maximale Dichteaenderung= 0.001055 5.Iteration, maximale Dichteaenderung= 0.000538 6.Iteration, maximale Dichteaenderung= 0.000280 7.Iteration, maximale Dichteaenderung= 0.000147 8.Iteration, maximale Dichteaenderung= 0.000078 9.Iteration, maximale Dichteaenderung= 0.000042 10.Iteration, maximale Dichteaenderung= 0.000022 11.Iteration, maximale Dichteaenderung= 0.000012 12.Iteration, maximale Dichteaenderung= 0.000018 13.Iteration, maximale Dichteaenderung= 0.000015 14.Iteration, maximale Dichteaenderung= 0.000015 15.Iteration, maximale Dichteaenderung= 0.000010 Werte nach der SCF-Konvergenz mit 16 Iterationen ========================== Alpha-Korrekturen und Beta-Werte nach den SCF-Iterationen (in eV) 1 2 3 4 5 6 1 0.0000 2 1.9872 0.0000 3 0.0000 2.5489 0.0000 4 0.0000 0.0000 2.0358 0.0000 5 1.9872 0.0000 0.0000 2.5489 0.0000 6 2.5613 0.0000 0.0000 0.0000 0.0000 0.0000 Atomare Wechselwirkungsintegrale Gamma 1 2 3 4 5 6 1 10.8400 2 5.1354 10.8400 3 3.9807 5.3651 10.8400 4 3.9807 3.9606 5.1338 10.8400 5 5.1354 3.9028 3.9606 5.3651 10.8400 6 5.3855 3.7380 2.9284 2.9284 3.7380 10.8400 HARTREE-FOCK-Matrix 1 2 3 4 5 6 1 -6.0568 2 -2.7438 -5.9192 3 0.0661 -4.9524 -5.8206 4 0.0661 0.0886 -2.9420 -5.8206 5 -2.7438 0.5937 0.0886 -4.9524 -5.9192 6 -5.0057 -0.2257 0.3853 0.3853 -0.2257 -6.1866 SCF-Eigenwerte (eV) und zugehoerige SCF-Eigenvektoren HOMO = 3 LUMO = 4 1 2 3 4 5 6 -13.5407 -10.7774 -10.0543 -2.8684 0.6628 0.8549 1 0.5170 -0.4789 0.0000 -0.3263 0.0000 -0.6300 2 0.4067 0.1314 -0.5786 -0.3695 0.4065 0.4252 3 0.3797 0.4445 -0.4065 0.3419 -0.5786 -0.2034 4 0.3797 0.4445 0.4065 0.3419 0.5786 -0.2034 5 0.4067 0.1314 0.5786 -0.3695 -0.4065 0.4252 6 0.3371 -0.5839 0.0000 0.6219 0.0000 0.3983 Summe ( MO-Energien * Anzahl Elektronen ) -68.7448 (eV) Dichte-Matrix (Elektronendichten und Bindungsordnungen) 1 2 3 4 5 6 1 0.9932 2 0.2947 1.0349 3 -0.0332 0.8960 1.0139 4 -0.0332 -0.0448 0.3530 1.0139 5 0.2947 -0.3042 -0.0448 0.8960 1.0349 6 0.9078 0.1208 -0.2632 -0.2632 0.1208 0.9091 ------------------------------------------------------------------------------ Molekulare Elektronen-Wechselwirkungs-Integrale G(k,l,j,i) MEW N= 6 M=HOMO= 3 LUMO= 4 UU= 3 VV= 3 J G( 3, 3, 4, 4)= 4.9632 eV K G( 3, 4, 3, 4)= 0.6875 eV FULVEN ========================================================================= WFPPP Plot log(Oscill) 10.0 . . . . . . . . . 1.00 . . * * . * * . * * . * * * * . * * * * . * * * * . * * * * . * * * * 0.10 * * * * . * * * * . * * * * . * * * * . * * * * * . * * * * * . * * * * * . * * * * * . * * * * * . * * * * * 0.01 * * * * * . * * * * * . * * * * * . * * * * * . * * * * * . * * * * * . * * * * * . * * * * * . 5 4 3 2 1 . ---- cm-1 --50000----------40000-----------30000----------20000-----------10000 Singulett CI mit 9 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 3.5380 28544 350 0.038 0.000 -0.353 0.000 0.353 2 3 - 4 98.6% 3.5977 0.000 -0.424 0.000 0.424 4 3 - 6 1.3% 7.2599 0.000 0.649 0.000 0.649 2 5.0659 40871 244 0.603 1.166 0.000 0.000 1.166 3 3 - 5 8.5% 7.3889 0.852 0.000 0.000 0.852 5 2 - 4 89.4% 5.3025 1.501 0.000 0.000 1.501 7 2 - 6 1.9% 8.0495 0.028 0.000 0.000 0.028 3 6.7853 54743 182 0.264 -0.666 0.000 0.000 0.666 3 3 - 5 25.9% 7.3889 0.852 0.000 0.000 0.852 5 2 - 4 3.2% 5.3025 1.501 0.000 0.000 1.501 7 2 - 6 1.8% 8.0495 0.028 0.000 0.000 0.028 8 1 - 4 68.8% 7.2094 0.020 0.000 0.000 0.020 4 7.2808 58741 170 0.264 0.000 0.643 0.000 0.643 2 3 - 4 1.3% 3.5977 0.000 -0.424 0.000 0.424 4 3 - 6 97.0% 7.2599 0.000 0.649 0.000 0.649 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Grundzustand Singulettzustand Nr. 0 mu Ladung mu - nu mu - nu mu - nu 1 0.007 1 - 2 0.295 1 - 5 0.295 1 - 6 0.908 2 -0.035 2 - 3 0.896 3 -0.014 3 - 4 0.353 4 -0.014 4 - 5 0.896 5 -0.035 6 0.091 Dipolmoment Komponenten X= 1.174 Y= 0.000 Z= 0.000 Gesamt= 1.174 Projektion auf die X - Y Ebene Skalierungsfaktor 1.00, 1 Angstrom = 2.50 cm 5 . Y 4 . . . . 6 . . . . . 1 . . . . . . . X . . . . . . . 3 . 2 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Angeregter Singulettzustand Nr. 1 mu Ladung mu - nu mu - nu mu - nu 1 -0.149 1 - 2 0.421 1 - 5 0.421 1 - 6 0.646 2 0.198 2 - 3 0.560 3 0.051 3 - 4 0.617 4 0.051 4 - 5 0.560 5 0.198 6 -0.348 Dipolmoment Komponenten X= -4.991 Y= 0.000 Z= 0.000 Gesamt= 4.991 Projektion auf die X - Y Ebene Skalierungsfaktor 1.00, 1 Angstrom = 2.50 cm 5 . Y 4 . . . . 6 . . . . . 1 . . . . . . . X . . . . . . . 3 . 2 Aenderungen Grundzustand -> Angeregter Zustand Nr. 1 der pi-Ladungen und Bindungsordnungen mu mu - nu mu - nu mu - nu 1 -0.156 1 - 2 0.126 1 - 5 0.126 1 - 6 -0.262 2 0.233 2 - 3 -0.336 3 0.065 3 - 4 0.264 4 0.065 4 - 5 -0.336 5 0.233 6 -0.439 Summe der Betraege der E-Dichte-Aenderungen CT= 1.190 Summe der Betraege der Bindungs-Aenderungen CB= 1.450 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Angeregter Singulettzustand Nr. 2 mu Ladung mu - nu mu - nu mu - nu 1 0.145 1 - 2 0.448 1 - 5 0.448 1 - 6 0.553 2 -0.166 2 - 3 0.649 3 0.044 3 - 4 0.266 4 0.044 4 - 5 0.649 5 -0.166 6 0.099 Dipolmoment Komponenten X= 1.138 Y= 0.000 Z= 0.000 Gesamt= 1.138 Projektion auf die X - Y Ebene Skalierungsfaktor 1.00, 1 Angstrom = 2.50 cm 5 . Y 4 . . . . 6 . . . . . 1 . . . . . . . X . . . . . . . 3 . 2 Aenderungen Grundzustand -> Angeregter Zustand Nr. 2 der pi-Ladungen und Bindungsordnungen mu mu - nu mu - nu mu - nu 1 0.138 1 - 2 0.154 1 - 5 0.154 1 - 6 -0.354 2 -0.131 2 - 3 -0.247 3 0.058 3 - 4 -0.087 4 0.058 4 - 5 -0.247 5 -0.131 6 0.008 Summe der Betraege der E-Dichte-Aenderungen CT= 0.525 Summe der Betraege der Bindungs-Aenderungen CB= 1.242 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Angeregter Singulettzustand Nr. 3 mu Ladung mu - nu mu - nu mu - nu 1 0.223 1 - 2 0.208 1 - 5 0.208 1 - 6 0.585 2 -0.016 2 - 3 0.593 3 -0.095 3 - 4 0.347 4 -0.095 4 - 5 0.593 5 -0.016 6 -0.001 Dipolmoment Komponenten X= 2.116 Y= 0.000 Z= 0.000 Gesamt= 2.116 Projektion auf die X - Y Ebene Skalierungsfaktor 1.00, 1 Angstrom = 2.50 cm 5 . Y 4 . . . . 6 . . . . . 1 . . . . . . . X . . . . . . . 3 . 2 Aenderungen Grundzustand -> Angeregter Zustand Nr. 3 der pi-Ladungen und Bindungsordnungen mu mu - nu mu - nu mu - nu 1 0.216 1 - 2 -0.087 1 - 5 -0.087 1 - 6 -0.322 2 0.019 2 - 3 -0.303 3 -0.081 3 - 4 -0.006 4 -0.081 4 - 5 -0.303 5 0.019 6 -0.091 Summe der Betraege der E-Dichte-Aenderungen CT= 0.509 Summe der Betraege der Bindungs-Aenderungen CB= 1.109 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Angeregter Singulettzustand Nr. 4 mu Ladung mu - nu mu - nu mu - nu 1 -0.326 1 - 2 0.024 1 - 5 0.024 1 - 6 0.711 2 0.080 2 - 3 0.551 3 0.086 3 - 4 0.566 4 0.086 4 - 5 0.551 5 0.080 6 -0.006 Dipolmoment Komponenten X= -2.502 Y= 0.000 Z= 0.000 Gesamt= 2.502 Projektion auf die X - Y Ebene Skalierungsfaktor 1.00, 1 Angstrom = 2.50 cm 5 . Y 4 . . . . 6 . . . . . 1 . . . . . . . X . . . . . . . 3 . 2 Aenderungen Grundzustand -> Angeregter Zustand Nr. 4 der pi-Ladungen und Bindungsordnungen mu mu - nu mu - nu mu - nu 1 -0.333 1 - 2 -0.271 1 - 5 -0.271 1 - 6 -0.197 2 0.115 2 - 3 -0.345 3 0.100 3 - 4 0.213 4 0.100 4 - 5 -0.345 5 0.115 6 -0.097 Summe der Betraege der E-Dichte-Aenderungen CT= 0.858 Summe der Betraege der Bindungs-Aenderungen CB= 1.641 ------------------------------------------------------------------------------ FULVEN ENDE WFPPP Version 2011.11.29