Discrete-continuous hybrid models are a popular means for describing elastic membrane-mediated particle interactions in and on lipid bilayers. Here, the continuous part is usually given by an approximation of the lipid membrane by an infinitely thin and sufficiently smooth hypersurface, whose elastic energy is determined by a Canham-Helfrich type functional. The discrete component results from modeling non-membrane particles as rigid discrete entities, which, depending on their configuration, induce local constraints on the membrane along the membrane-particle interfaces.
In this context, the interaction potential describes the optimal elastic energy of such hybrid systems with a fixed particle configuration. Correspondingly, the energy minimization principle yields that stationary particle configurations are given by the local minima of the interaction potential.
The main goal of this work is the proof of differentiability of the interaction potential for a selected class of models. This is accomplished using a variational approach that is already established in the literature in order to develop and apply robust numerical optimization methods for computing stationary particle configurations. Correspondingly, an additional focus is the derivation of a numerically accessible representation of the gradient, including its discretization and relevant numerical analysis.
The proof of differentiability is brought forward by an application of the implicit function theorem. The basis for this is so-called boundary preserving domain transformations, which are induced by suitable families of vector fields and which locally admit the reformulation of the minimization problem that is implicitly defined by the interaction potential with respect to a fixed particle configuration. This subsequently enables the representation of the gradient as a volume integral using matrix analysis methods.
The discretization of the partial differential equations for describing optimal membrane shapes is done via finite element methods. For particle methods with so-called curve restrictions a fictitious domain stabilized Nitsche method is developed, and for models with point value restrictions a conforming Galerking discretization is made possible by local QR transformations of the nodal finite element basis. For both cases suitable a priori error estimates are proven, and in addition also error estimates for the volume representation of the gradient are shown within that context.
These developed methods open up the domain of efficient simulation of macro structures by isotropic and anisotropic particles, which is illustrated with the aid of various example applications and by means of perturbed gradient methods.
Diskretkontinuierliche hybride Modelle sind ein populäres Mittel zur Beschreibung elastisch Membran-mediierter Interaktionen von Partikeln in und auf Doppellipidschichten. Der kontinuierliche Anteil ist hierbei üblicherweise durch eine Näherung der Lipidmembran als unendlich dünne und hinreichend glatte Hyperfläche gegeben, deren elastische Energie durch ein Canham-Helfrich-artiges Funktional bestimmt ist. Die diskrete Komponente ergibt sich aus der Modellierung von Nicht-Membran-Partikeln als diskrete und starre Objekte, die in Abhängigkeit ihrer Konfiguration lokale Restriktionen der Membran entlang der Membran-Partikel-Schnittstellen induzieren. In diesem Zusammenhang beschreibt das Interaktionspotential die optimale elastische Energie eines solchen hybriden Systems in Abhängigkeit einer festen Partikelkonfiguration. Dementsprechend folgt aus dem Energieminimierungsprinzip, dass stationäre Partikelkonfigurationen durch die lokalen Minima des Interaktionspotentials gegeben sind. Das Hauptziel dieser Arbeit ist der Beweis der Differenzierbarkeit des Interaktionspotentials für eine ausgewählte Klasse von Modellen im Sinne eines in der Literatur bereits eingeführten variationellen Ansatzes zum Zwecke der Entwicklung und Anwendung robuster numerischer Optimierungsverfahren zur Berechnung stationärer Partikelkonfigurationen. Entsprechend steht darüber hinaus ebenso die Herleitung einer numerisch zugänglichen Darstellung des Gradienten im Mittelpunkt, inklusive deren Diskretisierung und der relevanten numerischen Analysis. Der Differenzierbarkeitsbeweis wird erbracht durch eine Anwendung des Satzes über implizite Funktionen. Die Grundlage hierfür sind sogenannte randwerterhaltende Gebietstransformationen, die von geeigneten Familien von Vektorfeldern induziert werden und die lokal die Umformulierung der im Interaktionspotential implizit definierten Minimierungsprobleme bezüglich einer festen Referenzkonfiguration erlauben. Dies ermöglicht im Anschluss die Darstellung des Gradienten als Volumenintegral mittels Methoden der Matrix-Analysis. Die Diskretisierung der Differentialgleichungen zur Beschreibung optimaler Membranformen erfolgt über Finite-Elemente-Methoden. Für Partikelmodelle mit sogenannten Kurvenrestriktionen wird ein stabilisiertes Nitscheverfahren mit fiktiven Gebietsanteilen entwickelt, und für Modelle mit Punktwertrestriktionen wird durch lokale QR-Transformationen der nodalen Finite-Element-Basis eine konforme Galerkin-Diskretisierung ermöglicht. In beiden Fällen werden geeignete a priori Fehlerabschätzungen bewiesen, und ebenso werden in diesem Rahmen auch Fehlerabschätzungen für die Volumendarstellung des Gradienten gezeigt. Durch die entwickelten Verfahren wird die effiziente Simulation von Makrostrukturen durch isotrope und anisotrope Partikel erschlossen, was anhand verschiedener Beispielanwendungen und mit Hilfe gestörter Gradientenverfahren illustriert wird.
