Die Hydratation hat großen Einfluss auf biomolekulare Vorgänge in wässriger Phase wie zum Beispiel die Ligandenbindung oder die Struktur und Funktion von Proteinen. Die molekulare Simulation bietet hier ein vielseitiges Werkzeug zur Berechnung makroskopischer struktureller und thermodynamischer Eigenschaften von Flüssigkeiten und Solute/Solvent Systemen, beispielsweise radiale Paarverteilungsfunktionen oder freie Solvatationsenthalpien. Eine klassische Molekulardynamik Simulation auf atomarem Level ist jedoch sehr rechenaufwendig. Die Effizienz läßt sich durch den Einsatz empirischer Wechselwirkungspotenziale mit sphärischen Abschneidebedingungen und kleinen Simulationszellen erheblich steigern. Diese effizienten Simulationsbedingungen rufen aber systematische Fehler in den berechneten Observablen hervor. Die begrenzte Anzahl von Molekülen in der Simulationszelle erzeugt einen finte size Effekt, während die Näherung der Wechselwirkungspotenziale einen Abschneidefehler einbringt. Finite size Effekte können nur durch eine Vergrößerung der Simulationszelle vermieden werden. Abschneidefehler hingegen lassen sich durch den Einsatz geeigneter Simulationsalgorithmen minimieren oder korrigieren. Das Ziel dieser Arbeit ist die Steigerung der Effizienz molekularer Simulationen von Flüssigkeiten durch Anwendung möglichst keiner sphärischer Abschneideradien und Minimierung oder Korrektur der Abschneidefehler. Die erste verfolgte Strategie ist auf polare Flüssigkeiten wie Wasser beschränkt. Unter Berücksichtigung der dielektrischen Eigenschaften des Mediums wird ein effektives Coulomb Potenzial formuliert und für einige Wassermodelle parametrisiert. Es wird gezeigt, dass zahlreiche mit kleinen Abschneideradien berechnete Observablen dennoch Werte annehmen, die nah bei entsprechenden duch Ewald Summation erhaltener Referenzwerte liegen. Die zweite eingesetzte Technik basiert auf einer Kombination von molekularer Simulation und eindimensionaler RISM Integralgleichungstheorie. Ein Vergleich der unabhängig voneinander mit Simulation und Integralgleichung bestimmten radialen Paarverteilungsfunktionen von Wasser offenbart, dass beide Methoden exakt die gleichen Abschneideeffekte aufweisen. Über die RISM Integralgleichungstheorie können daher die Abschneidefehler von Simulationsergebnissen korrigiert werden. Mit Hilfe der thermodynamischen Integration wird die absolute freie Hydrationsenthalpie von Argon in Wasser und der elektrostatische Anteil des chemischen Potenzials von Wasser bei verschiedenen Abschneideradien bestimmt. Es wird gezeigt, dass die RISM Korrektur die Abschneidefehler zuverlässig verhersagen kann. Die neue RISM Korrekturmethode wird erfolgreich angewendet bei der effizienten Simulation von radialen Paarverteilungsfunktionen, potenzieller Energie und schließlich freien Hydratationsenthalpien.