Diese Arbeit befasst sich mit zeitaufgelösten Prozessen in molekularen Systemen. Dabei wurde sowohl die Wellenpaketdynamik nach Photoanregung betrachtet als auch spektrale Eigenschaften mittels Absorptions- und zweidimensionaler Spektroskopie untersucht. Zunächst widmet sich die Arbeit der Wellenpaket- und Populationsdynamik in zwei diabatischen, gekoppelten Zuständen. Nach impulsiver Anregung aus dem zu Beginn besetzten Zustand treten in der Populationsdynamik zwei deutlich verschiedene Oszillationen auf. Der langsamer variierende Populationstransfer besitzt die Periodendauer der Vibrationsbewegung und ist auf einen Wechsel der Zustände beim Durchlaufen des Wellenpakets durch die Kreuzungsregion der diabatischen Potentiale zurückzuführen. Die ultraschnelle Komponente mit einer Periodendauer von etwa 4 fs lässt sich als eine Art Rabi-Oszillation beschreiben, die durch die (zeitunabhängige) Kopplung hervorgerufen wird. Sie wurde mit Hilfe von analytischen Berechnungen ausführlich charakterisiert. Damit dieser Prozess auftreten kann müssen mehrere Bedingungen erfüllt werden: Das Wellenpaket muss über die Dauer der Oszillationen annähernd örtlich lokalisiert bleiben; dies ist an den Umkehrpunkten der Wellenpaketsbewegung der Fall. Die Amplitude der Oszillationen in den Populationen ist proportional zum Verhältnis der Kopplung zum Energieabstand der Zustände. Deshalb muss an den stationären Stellen die Kopplung groß im Vergleich zum Energieabstand sein. Die Amplitude der Oszillationen hängt außerdem von dem Populationsverhältnis und den Phasen der Komponenten des Wellenpakets in den beiden Zuständen ab. Die ultraschnellen Oszillationen bleiben auch in mehrdimensionalen Systemen mit unterschiedlichen Vibrationsfrequenzen je Freiheitsgrad erhalten. Das gleiche Modell wurde benutzt, um Ladungstransferprozesse mittels linearer und 2D-Spektroskopie zu untersuchen. Eine Kopplung an die Umgebung wurde, aufbauend auf einer Quanten-Master-Gleichung in Markov-Näherung, wellenfunktionsbasiert mittels eines Quantum-Jump-Algorithmus mit expliziter Dephasierung beschrieben. Dabei findet mit vorher definierten Wahrscheinlichkeiten zu jedem Zeitschritt einer von drei stochastischen Prozessen statt. Neben kohärenter Propagation können Sprünge in einen anderen Eigenzustand des Systems und Dephasierungen auftreten. Zwei Dissipationsparameter spielen dabei eine Rolle. Dies ist zum einen die Stärke der System-Bad-Kopplung, welche die Gesamtrate der Energierelaxation beschreibt. Weiterhin beeinflusst die Dephasierungskonstante den Verlust kohärenter Phasen ohne Energieänderung. Fallenzustände wurden identifiziert, die durch sehr geringe Sprungraten in niedrigere Zustände charakterisiert sind. Die Langlebigkeit kann durch die Form der Eigenfunktionen erklärt werden, die eine deutlich andere Wahrscheinlichkeitsverteilung als die der Nicht-Fallenzustände besitzen. Dadurch werden die in die Sprungraten eingehenden Matrixelemente klein. Das Absorptionsspektrum zeigt Peaks an der Stelle der Fallenzustände, da nur die Eigenfunktionen der Fallenzustände große Franck-Condon-Faktoren mit der Anfangswellenfunktion besitzen. Verschiedene Kombinationen der Dissipationsparameter führen zu Änderungen der relativen Peakintensitäten und der Peakbreiten. Die 2D-Spektren des Ladungstransfersystems werden störungstheoretisch über die Polarisation dritter Ordnung berechnet. Sie zeigen viele eng nebeneinander liegende Peaks in einer schachbrettmusterförmigen Anordnung, die sich auf Übergänge unter Mitwirkung der Fallenzustände zurückführen lassen. Höhere System-Bad-Kopplungen führen aufgrund der effizienten Energiedissipation zu einer Verschiebung zu kleineren Energien. Peaks, die mit schneller zerfallenden Fallenzuständen korrespondieren, bleichen schneller aus. Höhere Dephasierungskonstanten resultieren in verbreiterten Peaks. Um den Einfluss der Dissipation genauer zu charakterisieren, wurden gefilterte 2D-Spektren betrachtet. Dazu wurden Ausschnitte der Polarisation dritter Ordnung zu verschiedenen Zeiten fouriertransformiert. Längere Zeiten führen zu einer effektiveren Energierelaxation entlang der entsprechenden Zeitvariablen. Die Entvölkerung der höher liegenden Zustände lässt sich somit zeit- und energieaufgelöst betrachten. Weiterhin wurde gezeigt, dass sich der Zerfall eines einzelnen Peaks mit dem Populationsabfall des damit korrespondierenden Eigenzustandes in Einklang bringen lässt, obwohl die Zuordnung der Peaks im 2D-Spektrum zu Übergängen zwischen definierten Eigenzuständen nicht eindeutig ist. Mit dem benutzten eindimensionalen Modell können auch Ladungstransferprozesse in organischen gemischtvalenten Verbindungen beschrieben werden. Es wurde die Frage untersucht, welche Prozesse nach einem optisch induzierten Energietransfer in solchen Systemen ablaufen. Experimentelle Daten (aufgenommen im Arbeitskreis von Prof. Lambert) deuten auf eine schnelle interne Konversion (IC) gefolgt von Thermalisierung hin. Um dies theoretisch zu überprüfen, wurden Absorptionsspektren bei verschiedenen Temperaturen berechnet und mit den gemessenen transienten Spektren verglichen. Es findet sich, abhängig von der Stärke der elektronischen Kopplung, eine sehr gute bis gute Übereinstimmung, was die Annahme eines schnellen ICs stützt. Im letzten Teil der Arbeit wurden vibronische 2D-Spektren von molekularen Aggregaten betrachtet. Dazu wurde die zeitabhängige Schrödingergleichung für ein Monomer-, Dimer- und Trimersystem mit der Multi-Configuration Time-Dependent Hartree-Methode gelöst und die Polarisation nicht-störungstheoretisch berechnet. Der Hamiltonoperator des Trimers umfasst hierbei sieben gekoppelte elektronische Zustände und drei bzw. sechs Vibrationsfreiheitsgrade. Der betrachtete Photonenecho-Beitrag der Polarisation wurde mittels phasencodierter Laserpulse extrahiert. Die resultierenden Spektren sind geometrieabhängig, ein Winkel zwischen den Übergangsdipolmomenten der Monomere von 0° (180°) resultiert in einem H-Aggregat (J-Aggregat). Die Lage und Intensität der Peaks im rein elektronischen Trimer wurde analytisch erläutert. Die Spektren unter Einbeziehung der Vibration zeigen eine ausgeprägte vibronische Struktur. Es wurde gezeigt, wie die Spektren für höhere Aggregationsgrade durch die höhere Dichte an vibronischen Zuständen komplexer werden. Im J-Aggregat ist mit zunehmender Aggregation eine stärkere Rotverschiebung zu sehen. Das Spektrum des H-Aggregats zeigt eine im Vergleich zum J-Aggregat kompliziertere Struktur. Die Verwendung zweier Vibrationsfreiheitsgrade je Monomer führt zu Spektren mit überlappenden Peaks und einer zusätzlichen vibronischen Progression. Der Vergleich von Spektren verschiedener Mischungen von Monomer, Dimer und Trimer, entsprechend einem von Temperatur und Konzentration abhängigen Aggregationsgrad, zeigt den Einfluss dieser experimentellen Faktoren. Schließlich wurden mögliche Ansätze aufgezeigt, anhand der Spektren auf den Aggregationsgrad zu schließen.