In dieser Dissertation diskutiere ich theoretische Methoden zur Simulation von Grenzflächen zwischen Metallen und Molekülen auf atomarer Maßstabsebene, die für die Speicherung und Erzeugung "sauberer" Energie von Bedeutung sind, und wende sie an. Wir verwenden die Dichtefunktionaltheorie für das elektronische Subsystem und verschiedene Methoden wie die (quasi-)harmonische Näherung und die Pfadintegral-Molekulardynamik, um die Quanteneigenschaften des nuklearen Subsystems zu berücksichtigen. Wir berechnen den Isotopeneffekt auf die Arbeitsfunktion von Cyclohexan, das an der Rh(111)-Oberfläche adsorbiert wird, ein Effekt, der sich aus der Elektron-Phonon-Kopplung nur dann ergibt, wenn die nuklearen Freiheitsgrade quantenmechanisch behandelt werden. Deuteriertes Cyclohexan C6D12 hat einen größeren Adsorptionsabstand als gewöhnliches Cyclohexan. Pfadintegral-Molekulardynamiksimulationen zeigen auch eine temperaturabhängige Renormierung der elektronischen Zustandsdichte in diesem System. Schließlich befassen wir uns mit Oberflächenreaktionen auf einer geladenen metallischen Oberfläche. Wir stellen unsere Implementierung der Nudged-Elastic-Band-Methode (NEB) im i-PI-Paket vor und diskutieren ihre Leistungsfähigkeit. Anschließend setzen wir die Methode ein, um die Energiebarriere der Wasserspaltungsreaktion auf einer Pd(111)-Oberfläche zu berechnen, die einem elektrischen Feld unterschiedlicher Intensität ausgesetzt ist. Wir zeigen, dass die niedrigste Dissoziationsbarriere auftritt, wenn das Feld eine Stärke erreicht, die eine geometrische Frustration des auf der Oberfläche adsorbierten Wassermoleküls hervorruft, und dass die Nullpunktenergiebeiträge zur Barriere dieser Reaktion über den weiten Bereich der auf das System angelegten elektrischen Feldstärken nahezu konstant bleiben. Wir erklären dies durch eine gegenseitige Aufhebung der Rot- und Blauverschiebungen einzelner Schwingungsmoden zwischen Reaktant und Übergangszustand. / In this thesis, I discuss and apply theoretical methods for simulating interfaces between metals and molecules of relevance to "clean" energy storage and production on an atomistic scale. We use density-functional theory for the electronic subsystem and various methods such as (quasi-)harmonic approximation and path integral molecular dynamics to account for quantum properties of the nuclear subsystem, determining which methods are sufficient to grasp the essential phenomena while remaining computationally affordable. We calculate isotope effect on the work function of cyclohexane adsorbed on Rh(111) surface, an effect that emerges from electron-phonon coupling only when the nuclear degrees of freedom are treated quantum-mechanically. Deuterated cyclohexane C6D12 has larger adsorption distance than ordinary cyclohexane. Path integral molecular dynamics simulations also show a temperature-dependent renormalization of the electronic density of states in this system, induced by both thermal and quantum fluctuations of nuclei. Finally, we address surface reactions on a charged metallic surface. We present our implementation of the nudged elastic band (NEB) method in i-PI package and discuss its performance. We then employ the method to calculate the energy barrier of water splitting reaction on a Pd(111) surface subjected to electric fields of different strengths. We show that the lowest dissociation barrier takes place when the field reaches a strength that induces a geometric frustration of the water molecule adsorbed on the surface, and that the zero-point energy contributions to the barrier of this reaction remain nearly constant across the wide range of electric field strengths applied to the system. We explain this by a mutual cancellation of the red and blue shifts of individual vibrational modes between reactant and transition states.