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 ============================================================================== * BENZOL R=1.397 * MARVIN 01071018003D * * 12 12 0 0 0 0 999 V2000 * 0.0000 1.3970 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 * 1.2098 0.6985 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 * 1.2098 -0.6985 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 * 0.0000 -1.3970 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 * -1.2098 -0.6985 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 * -1.2098 0.6985 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 * 0.0000 2.3970 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0 * 2.0759 1.1985 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0 * 2.0759 -1.1985 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0 * 0.0000 -2.3970 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0 * -2.0759 -1.1985 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0 * -2.0759 1.1985 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0 * 1 7 8 0 0 0 0 * 1 6 8 0 0 0 0 * 1 2 2 0 0 0 0 * 2 8 8 0 0 0 0 * 2 3 8 0 0 0 0 * 3 9 8 0 0 0 0 * 3 4 2 0 0 0 0 * 4 10 8 0 0 0 0 * 4 5 8 0 0 0 0 * 5 11 8 0 0 0 0 * 5 6 2 0 0 0 0 * 6 12 8 0 0 0 0 * M END 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 1 . . Y 6 . 2 . . . . . . . . . . . . X . . . . . 5 . 3 . . . 4 ------------------------------------------------------------------------------ BENZOL R=1.397 ------------------------------------------------------------------------------ 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.00000 1.39700 0.00000 C pi1 1. 11.42000 0.58000 0.00000 2 6 1.20980 0.69850 0.00000 C pi1 1. 11.42000 0.58000 0.00000 3 6 1.20980 -0.69850 0.00000 C pi1 1. 11.42000 0.58000 0.00000 4 6 0.00000 -1.39700 0.00000 C pi1 1. 11.42000 0.58000 0.00000 5 6 -1.20980 -0.69850 0.00000 C pi1 1. 11.42000 0.58000 0.00000 6 6 -1.20980 0.69850 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 - 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 5 - 6 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.3970 0.0000 3 2.4197 1.3970 0.0000 4 2.7940 2.4197 1.3970 0.0000 5 2.4197 2.7939 2.4196 1.3970 0.0000 6 1.3970 2.4196 2.7939 2.4197 1.3970 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 0.0000 0.0000 0.0000 2.3180 0.0000 6 2.3180 0.0000 0.0000 0.0000 2.3180 0.0000 Atomare Wechselwirkungsintegrale Gamma 1 2 3 4 5 6 1 10.8400 2 5.2834 10.8400 3 3.8417 5.2833 10.8400 4 3.4929 3.8417 5.2834 10.8400 5 3.8417 3.4929 3.8418 5.2834 10.8400 6 5.2834 3.8418 3.4929 3.8417 5.2833 10.8400 Beginn SCF-Rechnung mit Elektronen-Wechselwirkung ========================== Werte nach der SCF-Konvergenz mit 0 Iterationen ========================== Alpha-Korrekturen und Beta-Werte nach 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 0.0000 0.0000 0.0000 2.3180 0.0000 6 2.3180 0.0000 0.0000 0.0000 2.3180 0.0000 Atomare Wechselwirkungsintegrale Gamma 1 2 3 4 5 6 1 10.8400 2 5.2834 10.8400 3 3.8417 5.2833 10.8400 4 3.4929 3.8417 5.2834 10.8400 5 3.8417 3.4929 3.8418 5.2834 10.8400 6 5.2834 3.8418 3.4929 3.8417 5.2833 10.8400 HARTREE-FOCK-Matrix 1 2 3 4 5 6 1 -6.0000 2 -4.0791 -6.0000 3 0.0000 -4.0791 -6.0000 4 0.5821 0.0000 -4.0791 -6.0000 5 0.0000 0.5822 0.0000 -4.0791 -6.0000 6 -4.0791 0.0000 0.5822 0.0000 -4.0791 -6.0000 SCF-Eigenwerte (eV) und zugehoerige SCF-Eigenvektoren HOMO = 3 LUMO = 4 1 2 3 4 5 6 -13.5761 -10.6613 -10.6613 -1.3387 -1.3387 1.5761 1 0.4082 -0.5771 0.0184 0.0039 -0.5773 0.4082 2 0.4082 -0.3044 -0.4906 0.4980 0.2920 -0.4082 3 0.4082 0.2726 -0.5089 -0.5019 0.2853 0.4082 4 0.4082 0.5771 -0.0184 0.0039 -0.5773 -0.4082 5 0.4082 0.3044 0.4906 0.4980 0.2920 0.4082 6 0.4082 -0.2726 0.5089 -0.5019 0.2853 -0.4082 Summe ( MO-Energien * Anzahl Elektronen ) -69.7973 (eV) Dichte-Matrix (Elektronendichten und Bindungsordnungen) 1 2 3 4 5 6 1 1.0000 2 0.6667 1.0000 3 0.0000 0.6667 1.0000 4 -0.3333 0.0000 0.6667 1.0000 5 0.0000 -0.3333 0.0000 0.6667 1.0000 6 0.6667 0.0000 -0.3333 0.0000 0.6667 1.0000 ------------------------------------------------------------------------------ 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)= 5.8640 eV K G( 3, 4, 3, 4)= 1.4759 eV BENZOL R=1.397 ========================================================================= WFPPP Plot log(Oscill) 10.0 . . . . . . . . . * 1.00 * . * . * . * . * . * . * . * . * . * 0.10 * . * . * . * . * . * . * . * . * . * 0.01 * . * . * . * . * . * . * . * . 4 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 4.7600 38403 260 0.000 0.000 0.000 0.000 0.000 3 3 - 5 50.0% 5.7912 -0.987 -0.038 0.000 0.988 5 2 - 4 50.0% 5.7912 -0.987 -0.038 0.000 0.988 2 5.9998 48406 206 0.000 0.000 0.000 0.000 0.000 2 3 - 4 50.0% 6.4104 0.038 -0.987 0.000 0.988 6 2 - 5 50.0% 6.4104 -0.038 0.987 0.000 0.988 3 6.8216 55036 181 1.165 -1.397 0.000 0.000 1.397 3 3 - 5 49.9% 5.7912 -0.987 -0.038 0.000 0.988 5 2 - 4 49.9% 5.7912 -0.987 -0.038 0.000 0.988 4 6.8216 55036 181 1.165 0.000 1.397 0.000 1.397 2 3 - 4 49.9% 6.4104 0.038 -0.987 0.000 0.988 6 2 - 5 49.9% 6.4104 -0.038 0.987 0.000 0.988 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Grundzustand Singulettzustand Nr. 0 mu Ladung mu - nu mu - nu mu - nu 1 0.000 1 - 2 0.667 1 - 6 0.667 2 0.000 2 - 3 0.667 3 0.000 3 - 4 0.667 4 0.000 4 - 5 0.667 5 0.000 5 - 6 0.667 6 0.000 Dipolmoment Komponenten X= 0.000 Y= 0.000 Z= 0.000 Gesamt= 0.000 Projektion auf die X - Y Ebene Skalierungsfaktor 1.00, 1 Angstrom = 2.50 cm 1 . . Y 6 . 2 . . . . . . . . . . . . X . . . . . 5 . 3 . . . 4 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Angeregter Singulettzustand Nr. 1 mu Ladung mu - nu mu - nu mu - nu 1 0.000 1 - 2 0.500 1 - 6 0.500 2 0.000 2 - 3 0.500 3 0.000 3 - 4 0.500 4 0.000 4 - 5 0.500 5 0.000 5 - 6 0.500 6 0.000 Dipolmoment Komponenten X= 0.000 Y= 0.000 Z= 0.000 Gesamt= 0.000 Projektion auf die X - Y Ebene Skalierungsfaktor 1.00, 1 Angstrom = 2.50 cm 1 . . Y 6 . 2 . . . . . . . . . . . . X . . . . . 5 . 3 . . . 4 Aenderungen Grundzustand -> Angeregter Zustand Nr. 1 der pi-Ladungen und Bindungsordnungen mu mu - nu mu - nu mu - nu 1 0.000 1 - 2 -0.167 1 - 6 -0.167 2 0.000 2 - 3 -0.167 3 0.000 3 - 4 -0.167 4 0.000 4 - 5 -0.167 5 0.000 5 - 6 -0.167 6 0.000 Summe der Betraege der E-Dichte-Aenderungen CT= 0.000 Summe der Betraege der Bindungs-Aenderungen CB= 1.000 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Angeregter Singulettzustand Nr. 2 mu Ladung mu - nu mu - nu mu - nu 1 0.000 1 - 2 0.500 1 - 6 0.500 2 0.000 2 - 3 0.500 3 0.000 3 - 4 0.500 4 0.000 4 - 5 0.500 5 0.000 5 - 6 0.500 6 0.000 Dipolmoment Komponenten X= 0.000 Y= 0.000 Z= 0.000 Gesamt= 0.000 Projektion auf die X - Y Ebene Skalierungsfaktor 1.00, 1 Angstrom = 2.50 cm 1 . . Y 6 . 2 . . . . . . . . . . . . X . . . . . 5 . 3 . . . 4 Aenderungen Grundzustand -> Angeregter Zustand Nr. 2 der pi-Ladungen und Bindungsordnungen mu mu - nu mu - nu mu - nu 1 0.000 1 - 2 -0.167 1 - 6 -0.167 2 0.000 2 - 3 -0.167 3 0.000 3 - 4 -0.167 4 0.000 4 - 5 -0.167 5 0.000 5 - 6 -0.167 6 0.000 Summe der Betraege der E-Dichte-Aenderungen CT= 0.000 Summe der Betraege der Bindungs-Aenderungen CB= 1.000 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Angeregter Singulettzustand Nr. 3 mu Ladung mu - nu mu - nu mu - nu 1 0.000 1 - 2 0.500 1 - 6 0.500 2 0.000 2 - 3 0.500 3 0.000 3 - 4 0.500 4 0.000 4 - 5 0.500 5 0.000 5 - 6 0.500 6 0.000 Dipolmoment Komponenten X= 0.000 Y= 0.000 Z= 0.000 Gesamt= 0.000 Projektion auf die X - Y Ebene Skalierungsfaktor 1.00, 1 Angstrom = 2.50 cm 1 . . Y 6 . 2 . . . . . . . . . . . . X . . . . . 5 . 3 . . . 4 Aenderungen Grundzustand -> Angeregter Zustand Nr. 3 der pi-Ladungen und Bindungsordnungen mu mu - nu mu - nu mu - nu 1 0.000 1 - 2 -0.167 1 - 6 -0.167 2 0.000 2 - 3 -0.167 3 0.000 3 - 4 -0.167 4 0.000 4 - 5 -0.167 5 0.000 5 - 6 -0.167 6 0.000 Summe der Betraege der E-Dichte-Aenderungen CT= 0.000 Summe der Betraege der Bindungs-Aenderungen CB= 1.000 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Angeregter Singulettzustand Nr. 4 mu Ladung mu - nu mu - nu mu - nu 1 0.000 1 - 2 0.500 1 - 6 0.500 2 0.000 2 - 3 0.500 3 0.000 3 - 4 0.500 4 0.000 4 - 5 0.500 5 0.000 5 - 6 0.500 6 0.000 Dipolmoment Komponenten X= 0.000 Y= 0.000 Z= 0.000 Gesamt= 0.000 Projektion auf die X - Y Ebene Skalierungsfaktor 1.00, 1 Angstrom = 2.50 cm 1 . . Y 6 . 2 . . . . . . . . . . . . X . . . . . 5 . 3 . . . 4 Aenderungen Grundzustand -> Angeregter Zustand Nr. 4 der pi-Ladungen und Bindungsordnungen mu mu - nu mu - nu mu - nu 1 0.000 1 - 2 -0.167 1 - 6 -0.167 2 0.000 2 - 3 -0.167 3 0.000 3 - 4 -0.167 4 0.000 4 - 5 -0.167 5 0.000 5 - 6 -0.167 6 0.000 Summe der Betraege der E-Dichte-Aenderungen CT= 0.000 Summe der Betraege der Bindungs-Aenderungen CB= 1.000 ------------------------------------------------------------------------------ BENZOL R=1.397 ENDE WFPPP Version 2011.11.29