Die schnelle, computergestützte Berechnung thermodynamischer Eigenschaften ist die Basis zur in-silico-Untersuchung supramolekularer Verbindungen in Lösung. So können Freie Bindungsenthalpien durch statistisches Sampling der intermolekularen Freien Energiefläche bestimmt werden. Diese Technik ist im Rahmen der Molekulardynamik und Monte-Carlo-Methodik bereits etabliert, jedoch stets mit einem hohen numerischen Aufwand verbunden, der auch auf die explizite Behandlung von Solvensmolekülen zurückgeht. Durch Kombination von Simulations- und Integralgleichungsmethoden kann dieser Aufwand maßgeblich reduziert werden. Im Rahmen der vorliegenden Arbeit werden hybride Simulations-/Integralgleichungsmethoden entwickelt und ihre Eignung zur Berechnung struktureller und thermodynamischer Eigenschaften von Kronenetherkomplexen in Lösung gezeigt. Das Lösungsmittel wird dabei jeweils durch radiale Verteilungsfunktionen auf der Basis der RISM-Integralgleichungsmethode modelliert. Im einzelnen werden Verfahren zur lokalen und globalen Komplexstruktur-Optimierung und zur thermodynamischen Simulation vorgestellt. Die Kombination von Monte-Carlo-Simulation mit RISM kann zur schnellen Berechnung Freier Energieflächen in Lösung dienen. Auf der Basis dieses Verfahrens wird die Komplexbildung von Kronenethern mit Alkalimetallen und molekularen Ammoniumionen in den Lösungsmitteln Wasser, Methanol und Acetonitril untersucht. Dabei finden auch lokale Aspekte der Energiehyperflächen in diesen Systemen Erwähnung. Beim Einbringen chiraler Gruppen in das beschriebene System Kronenether-Ammoniumion kann oft eine lösungsmittelabhängige chirale Differenzierung beobachtet werden. Zur Berechnung dieser differenziellen Größe wird RISM mit einem genetischen Algorithmus und einem lokalen Optimierer verknüpft. Auf diese Weise gelingt die Berechnung der chiralen Selektion stereoisomerer Komplexpaare im Lösungsmittel Methanol.