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 ============================================================================== * HEXATRIENE * * * 14 13 0 0 0 0 0 0 0 0 0 0 * -1.2076 2.8328 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 * -0.2814 3.4038 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0 * -2.1342 3.3991 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0 * -1.2027 1.4889 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 * -2.1534 0.9506 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0 * 0.0000 0.6766 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 * 0.9553 1.2079 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0 * -0.0001 -0.6767 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 * -0.9554 -1.2079 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0 * 1.2027 -1.4889 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 * 2.1533 -0.9505 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0 * 1.2077 -2.8327 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 * 0.2816 -3.4038 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0 * 2.1343 -3.3990 0.0000 H 0 0 0 0 0 0 0 0 0 0 0 0 * 1 2 1 0 0 0 0 * 1 3 1 0 0 0 0 * 1 4 2 0 0 0 0 * 4 5 1 0 0 0 0 * 4 6 1 0 0 0 0 * 6 7 1 0 0 0 0 * 6 8 2 0 0 0 0 * 8 9 1 0 0 0 0 * 8 10 1 0 0 0 0 * 10 11 1 0 0 0 0 * 10 12 2 0 0 0 0 * 12 13 1 0 0 0 0 * 12 14 1 0 0 0 0 Keine weitere Zeile gefunden. Keine weitere Zeile gefunden. ============================================================================== Ende der Eingabezeilen. ------------------------------------------------------------------------------ Summenformel fuer das Gesamtsystem H 8 C 6 ------------------------------------------------------------------------------ ---- Selektion der Pi-Atome 6 Pi-Atome gefunden 5 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 . . . . . . . . 2 . . . . Y 3 . . . . . . . . . . . . X . . . . . 4 . . . . . 5 . . . . . . . 6 ------------------------------------------------------------------------------ HEXATRIENE ------------------------------------------------------------------------------ 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 -1.20760 2.83280 0.00000 C pi1 1. 11.42000 0.58000 0.00000 2 6 -1.20270 1.48890 0.00000 C pi1 1. 11.42000 0.58000 0.00000 3 6 0.00000 0.67660 0.00000 C pi1 1. 11.42000 0.58000 0.00000 4 6 -0.00010 -0.67670 0.00000 C pi1 1. 11.42000 0.58000 0.00000 5 6 1.20270 -1.48890 0.00000 C pi1 1. 11.42000 0.58000 0.00000 6 6 1.20770 -2.83270 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 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.3439 0.0000 3 2.4713 1.4513 0.0000 4 3.7114 2.4771 1.3533 0.0000 5 4.9484 3.8280 2.4771 1.4513 0.0000 6 6.1589 4.9484 3.7113 2.4713 1.3438 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 0.0000 0.0000 0.0000 0.0000 2.3180 0.0000 Atomare Wechselwirkungsintegrale Gamma 1 2 3 4 5 6 1 10.8400 2 5.3883 10.8400 3 3.7895 5.1801 10.8400 4 2.8570 3.7837 5.3694 10.8400 5 2.2940 2.7925 3.7838 5.1800 10.8400 6 1.9231 2.2940 2.8571 3.7896 5.3885 10.8400 Beginn SCF-Rechnung mit Elektronen-Wechselwirkung ========================== 1.Iteration, maximale Dichteaenderung= 0.041349 2.Iteration, maximale Dichteaenderung= 0.009391 3.Iteration, maximale Dichteaenderung= 0.004643 4.Iteration, maximale Dichteaenderung= 0.002333 5.Iteration, maximale Dichteaenderung= 0.001183 6.Iteration, maximale Dichteaenderung= 0.000603 7.Iteration, maximale Dichteaenderung= 0.000308 8.Iteration, maximale Dichteaenderung= 0.000158 9.Iteration, maximale Dichteaenderung= 0.000081 10.Iteration, maximale Dichteaenderung= 0.000042 11.Iteration, maximale Dichteaenderung= 0.000021 12.Iteration, maximale Dichteaenderung= 0.000011 13.Iteration, maximale Dichteaenderung= 0.000013 Werte nach der SCF-Konvergenz mit 14 Iterationen ========================== Alpha-Korrekturen und Beta-Werte nach den SCF-Iterationen (in eV) 1 2 3 4 5 6 1 0.0000 2 2.6072 0.0000 3 0.0000 1.9983 0.0000 4 0.0000 0.0000 2.5617 0.0000 5 0.0000 0.0000 0.0000 1.9983 0.0000 6 0.0000 0.0000 0.0000 0.0000 2.6072 0.0000 Atomare Wechselwirkungsintegrale Gamma 1 2 3 4 5 6 1 10.8400 2 5.3883 10.8400 3 3.7895 5.1801 10.8400 4 2.8570 3.7837 5.3694 10.8400 5 2.2940 2.7925 3.7838 5.1800 10.8400 6 1.9231 2.2940 2.8571 3.7896 5.3885 10.8400 HARTREE-FOCK-Matrix 1 2 3 4 5 6 1 -6.0000 2 -5.1686 -6.0000 3 0.0000 -2.7967 -6.0000 4 0.4048 0.0000 -4.9998 -6.0000 5 0.0000 0.0477 0.0000 -2.7966 -6.0000 6 -0.1212 0.0000 0.4048 0.0000 -5.1687 -6.0000 SCF-Eigenwerte (eV) und zugehoerige SCF-Eigenvektoren HOMO = 3 LUMO = 4 1 2 3 4 5 6 -12.9932 -11.6326 -9.7126 -2.2874 -0.3674 0.9932 1 0.2739 -0.4674 -0.4544 -0.4544 -0.4674 0.2739 2 0.4042 -0.5085 -0.2794 0.2794 0.5085 -0.4042 3 0.5115 -0.1516 0.4641 0.4641 -0.1516 0.5115 4 0.5115 0.1515 0.4642 -0.4642 -0.1515 -0.5115 5 0.4042 0.5085 -0.2794 -0.2794 0.5085 0.4042 6 0.2739 0.4674 -0.4544 0.4544 -0.4674 -0.2739 Summe ( MO-Energien * Anzahl Elektronen ) -68.6768 (eV) Dichte-Matrix (Elektronendichten und Bindungsordnungen) 1 2 3 4 5 6 1 1.0000 2 0.9507 1.0000 3 0.0000 0.3082 1.0000 4 -0.2833 0.0000 0.9081 1.0000 5 0.0000 -0.0342 0.0000 0.3082 1.0000 6 0.1261 0.0000 -0.2833 0.0000 0.9507 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.0404 eV K G( 3, 4, 3, 4)= 1.2388 eV HEXATRIENE ========================================================================= WFPPP Plot log(Oscill) 10.0 . . . . . . . . * . * 1.00 * . * . * . * . * . * . * . * . * . * 0.10 * . * . * . * . * . * . * . * . * . * 0.01 * . * . * . * . * . * . * . * . 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 4.8324 38987 256 1.333 -0.458 1.715 0.000 1.775 2 3 - 4 99.5% 4.8624 -0.440 1.738 0.000 1.793 2 5.9433 47950 208 0.000 0.000 0.000 0.000 0.000 3 3 - 5 47.5% 6.4732 0.000 0.000 0.000 0.000 5 2 - 4 47.5% 6.4732 0.000 0.000 0.000 0.000 7 2 - 6 2.5% 9.2498 0.000 0.000 0.000 0.000 9 1 - 5 2.5% 9.2498 0.000 0.000 0.000 0.000 3 6.7715 54632 183 0.000 0.000 0.000 0.000 0.000 3 3 - 5 49.0% 6.4732 0.000 0.000 0.000 0.000 5 2 - 4 49.0% 6.4732 0.000 0.000 0.000 0.000 7 2 - 6 1.0% 9.2498 0.000 0.000 0.000 0.000 9 1 - 5 1.0% 9.2498 0.000 0.000 0.000 0.000 4 6.8525 55285 180 0.000 0.000 0.000 0.000 0.000 4 3 - 6 50.0% 7.5246 0.041 -0.067 0.000 0.079 8 1 - 4 50.0% 7.5246 0.041 -0.067 0.000 0.079 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Grundzustand Singulettzustand Nr. 0 mu Ladung mu - nu mu - nu mu - nu 1 0.000 1 - 2 0.951 2 0.000 2 - 3 0.308 3 0.000 3 - 4 0.908 4 0.000 4 - 5 0.308 5 0.000 5 - 6 0.951 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 . . . . . . . . 2 . . . . Y 3 . . . . . . . . . . . . X . . . . . 4 . . . . . 5 . . . . . . . 6 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Angeregter Singulettzustand Nr. 1 mu Ladung mu - nu mu - nu mu - nu 1 0.000 1 - 2 0.696 2 0.000 2 - 3 0.564 3 0.000 3 - 4 0.479 4 0.000 4 - 5 0.564 5 0.000 5 - 6 0.696 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 . . . . . . . . 2 . . . . Y 3 . . . . . . . . . . . . X . . . . . 4 . . . . . 5 . . . . . . . 6 Aenderungen Grundzustand -> Angeregter Zustand Nr. 1 der pi-Ladungen und Bindungsordnungen mu mu - nu mu - nu mu - nu 1 0.000 1 - 2 -0.255 2 0.000 2 - 3 0.256 3 0.000 3 - 4 -0.429 4 0.000 4 - 5 0.256 5 0.000 5 - 6 -0.255 6 0.000 Summe der Betraege der E-Dichte-Aenderungen CT= 0.000 Summe der Betraege der Bindungs-Aenderungen CB= 1.451 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Angeregter Singulettzustand Nr. 2 mu Ladung mu - nu mu - nu mu - nu 1 0.000 1 - 2 0.643 2 0.000 2 - 3 0.334 3 0.000 3 - 4 0.610 4 0.000 4 - 5 0.334 5 0.000 5 - 6 0.643 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 . . . . . . . . 2 . . . . Y 3 . . . . . . . . . . . . X . . . . . 4 . . . . . 5 . . . . . . . 6 Aenderungen Grundzustand -> Angeregter Zustand Nr. 2 der pi-Ladungen und Bindungsordnungen mu mu - nu mu - nu mu - nu 1 0.000 1 - 2 -0.307 2 0.000 2 - 3 0.026 3 0.000 3 - 4 -0.298 4 0.000 4 - 5 0.026 5 0.000 5 - 6 -0.307 6 0.000 Summe der Betraege der E-Dichte-Aenderungen CT= 0.000 Summe der Betraege der Bindungs-Aenderungen CB= 0.965 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Angeregter Singulettzustand Nr. 3 mu Ladung mu - nu mu - nu mu - nu 1 0.000 1 - 2 0.623 2 0.000 2 - 3 0.348 3 0.000 3 - 4 0.647 4 0.000 4 - 5 0.347 5 0.000 5 - 6 0.623 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 . . . . . . . . 2 . . . . Y 3 . . . . . . . . . . . . X . . . . . 4 . . . . . 5 . . . . . . . 6 Aenderungen Grundzustand -> Angeregter Zustand Nr. 3 der pi-Ladungen und Bindungsordnungen mu mu - nu mu - nu mu - nu 1 0.000 1 - 2 -0.327 2 0.000 2 - 3 0.039 3 0.000 3 - 4 -0.261 4 0.000 4 - 5 0.039 5 0.000 5 - 6 -0.327 6 0.000 Summe der Betraege der E-Dichte-Aenderungen CT= 0.000 Summe der Betraege der Bindungs-Aenderungen CB= 0.994 ------------------------------------------------------------------------------ ------------------------------------------------------------------------------ pi-Ladungen und Bindungsordnungen Angeregter Singulettzustand Nr. 4 mu Ladung mu - nu mu - nu mu - nu 1 0.000 1 - 2 0.713 2 0.000 2 - 3 0.231 3 0.000 3 - 4 0.431 4 0.000 4 - 5 0.231 5 0.000 5 - 6 0.713 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 . . . . . . . . 2 . . . . Y 3 . . . . . . . . . . . . X . . . . . 4 . . . . . 5 . . . . . . . 6 Aenderungen Grundzustand -> Angeregter Zustand Nr. 4 der pi-Ladungen und Bindungsordnungen mu mu - nu mu - nu mu - nu 1 0.000 1 - 2 -0.238 2 0.000 2 - 3 -0.077 3 0.000 3 - 4 -0.477 4 0.000 4 - 5 -0.077 5 0.000 5 - 6 -0.238 6 0.000 Summe der Betraege der E-Dichte-Aenderungen CT= 0.000 Summe der Betraege der Bindungs-Aenderungen CB= 1.106 ------------------------------------------------------------------------------ HEXATRIENE ENDE WFPPP Version 2011.11.29