The subject of this thesis is the rigorous passage from discrete systems to continuum models via variational methods. The first part of this work studies a discrete model describing a one-dimensional chain of atoms with finite range interactions of Lennard-Jones type. We derive an expansion of the ground state energy using \(\Gamma\)-convergence. In particular, we show that a variant of the Cauchy-Born rule holds true for the model under consideration. We exploit this observation to derive boundary layer energies due to asymmetries of the lattice at the boundary or at cracks of the specimen. Hereby we extend several results obtained previously for models involving only nearest and next-to-nearest neighbour interactions by Braides and Cicalese and Scardia, Schlömerkemper and Zanini. The second part of this thesis is devoted to the analysis of a quasi-continuum (QC) method. To this end, we consider the discrete model studied in the first part of this thesis as the fully atomistic model problem and construct an approximation based on a QC method. We show that in an elastic setting the expansion by \(\Gamma\)-convergence of the fully atomistic energy and its QC approximation coincide. In the case of fracture, we show that this is not true in general. In the case of only nearest and next-to-nearest neighbour interactions, we give sufficient conditions on the QC approximation such that, also in case of fracture, the minimal energies of the fully atomistic energy and its approximation coincide in the limit. / Der Gegenstand dieser Doktorarbeit ist der mathematisch strenge Übergang von diskreten Systemen zu kontinuierlichen Modellen. Im ersten Teil der Arbeit betrachten wir ein diskretes Modell für eine Kette von Atomen welche durch Lennard-Jones-artige Potentiale mit endlicher Reichweite miteinander wechselwirken. Wir leiten eine Entwicklung der Energie des Grundzustands mittels \(\Gamma\)-Konvergenz her. Dabei zeigen wir zunächst die Gültigkeit einer Variante der sogenannten Cauchy-Born Regel für unser Modell. Diese Beobachtung benutzen wir um Oberflächenenergien herzuleiten. Diese entstehen zum einen am Rand der Kette, aber auch an den Oberflächen die durch Brüche im kontinuierlichem Material gebildet werden. Dabei werden einige Resultate verallgemeinert, welche von Braides und Cicalese bzw. Scardia, Schlömerkemper und Zanini für Modelle bewiesen wurden in denen nur nächste und übernächste Nachbarn wechselwirken. Im zweiten Teil der Arbeit widmen wir uns der Untersuchung einer Quasi-continuum (QC) Methode. Wir betrachten die Energie aus dem ersten Teil als atomistisches Modellproblem und konstruieren eine Annäherung dieser Energie durch eine QC Methode. Im Falle von elastischen Verhalten zeigen wir, dass die Entwicklung mittels \(\Gamma\)-Konvergenz der atomistischen Energie und ihrer QC Approximation übereinstimmen. Falls allerdings Brüche auftreten gilt dies im Allgemeinen nicht mehr. Für Modelle in denen nur nächste und übernächste Nachbarn miteinander wechselwirken leiten wir hinreichende Bedingungen her die, auch im Falle von Brüchen, garantieren, dass das Minimum der atomistischen Energie und das Minimum der QC Approximation den gleichen Grenzwert haben.

The work in this thesis contains three main topics. These are the passage from discrete to continuous models by means of $\Gamma$-convergence, random as well as periodic homogenization and fracture enabled by non-convex Lennard-Jones type interaction potentials. Each of them is discussed in the following. We consider a discrete model given by a one-dimensional chain of particles with randomly distributed interaction potentials. Our interest lies in the continuum limit, which yields the effective behaviour of the system. This limit is achieved as the number of atoms tends to infinity, which corresponds to a vanishing distance between the particles. The starting point of our analysis is an energy functional in a discrete system; its continuum limit is obtained by variational $\Gamma$-convergence. The $\Gamma$-convergence methods are combined with a homogenization process in the framework of ergodic theory, which allows to focus on heterogeneous systems. On the one hand, composite materials or materials with impurities are modelled by a stochastic or periodic distribution of particles or interaction potentials. On the other hand, systems of one species of particles can be considered as random in cases when the orientation of particles matters. Nanomaterials, like chains of atoms, molecules or polymers, are an application of the heterogeneous chains in experimental sciences. A special interest is in fracture in such heterogeneous systems. We consider interaction potentials of Lennard-Jones type. The non-standard growth conditions and the convex-concave structure of the Lennard-Jones type interactions yield mathematical difficulties, but allow for fracture. The interaction potentials are long-range in the sense that their modulus decays slower than exponential. Further, we allow for interactions beyond nearest neighbours, which is also referred to as long-range. The main mathematical issue is to bring together the Lennard-Jones type interactions with ergodic theorems in the limiting process as the number of particles tends to infinity. The blow up at zero of the potentials prevents from using standard extensions of the Akcoglu-Krengel subadditive ergodic theorem. We overcome this difficulty by an approximation of the interaction potentials which shows suitable Lipschitz and Hölder regularity. Beyond that, allowing for continuous probability distributions instead of only finitely many different potentials leads to a further challenge. The limiting integral functional of the energy by means of $\Gamma$-convergence involves a homogenized energy density and allows for fracture, but without a fracture contribution in the energy. In order to refine this result, we rescale our model and consider its $\Gamma$-limit, which is of Griffith's type consisting of an elastic part and a jump contribution. In a further approach we study fracture at the level of the discrete energies. With an appropriate definition of fracture in the discrete setting, we define a fracture threshold separating the region of elasticity from that of fracture and consider the pointwise convergence of this threshold. This limit turns out to coincide with the one obtained in the variational $\Gamma$-convergence approach. / Diese Arbeit vereinigt im Wesentlichen drei Themen: Den Übergang von diskreten zu kontinuierlichen Modellen mittels $\Gamma$-Konvergenz, stochastische sowie periodische Homogenisierung, sowie Bruchmechanik, die durch nicht-konvexe Wechselwirkungspotentiale vom Lennard-Jones-Typ ermöglicht wird. Jedes dieser drei Themen wird im Folgenden diskutiert. Wir betrachten ein diskretes Modell, bestehend aus einer eindimensionale Kette von Teilchen mit zufällig verteilten Wechselwirkungspotentialen. Wir sind am Kontinuumsgrenzwert interessiert, welcher das effektive Verhalten des Systems widerspiegelt. In diesem Grenzwert läuft die Anzahl der Atome gegen unendlich, was einem verschwindenden Abstand zwischen den Teilchen entspricht. Ausgehend von einer Energie eines diskreten Systems erhalten wir den Kontinuumsgrenzwert durch die variationelle Methode der $\Gamma$-Konvergenz, welche den Übergang zum kontinuierlichen System liefert. Die $\Gamma$-Konvergenzmethoden werden im Rahmen der Ergodentheorie mit einem Homogenisierungsprozess kombiniert, wodurch die Betrachtung heterogener Systeme möglich wird. Einerseits werden Verbundwerkstoffe oder Materialien mit Verunreinigungen durch eine stochastische oder periodische Verteilung der Teilchen oder der Wechselwirkungspotentiale modelliert. Andererseits können Systeme einer Teilchenart als zufällig angesehen werden, wenn die Orientierung der Teilchen von Bedeutung ist. Nanomaterialien wie Ketten von Atomen, Molekülen oder Polymeren bieten eine Anwendung des Modells der heterogenen Ketten in den experimentellen Wissenschaften. Von besonderem Interesse ist das Auftreten von Brüchen in diesen heterogenen Systemen. Wir betrachten Wechselwirkungspotentiale vom Lennard-Jones Typ. Die nicht-standardisierten Wachstumsbedingungen und die konvex-konkave Struktur der Lennard-Jones Potentiale werfen mathematische Schwierigkeiten auf, ermöglichen jedoch das Auftreten von Brüchen. Die Wechselwirkungen gelten als langreichweitig in dem Sinne, dass ihr Betrag langsamer als exponentiell abfällt. Darüber hinaus betrachten wir Wechselwirkungen jenseits der nächsten Nachbarn, was ebenfalls als langreichweitig bezeichnet wird. Eine der größten mathematischen Schwierigkeiten besteht darin, die Wechselwirkungen vom Lennard-Jones Typ mit den Ergodensätzen zusammenzuführen. Die Singularität der Potentiale bei Null erlaubt keine Verwendung der Standardtechniken zur Erweiterung des subadditiven Ergodensatzes von Akcoglu-Krengel. Die Lösung dieses Problems ist eine Approximation der Wechselwirkungspotentiale, welche eine geeignete Lipschitz- und Hölder-Regularität besitzt. Darüber hinaus stellt die Verwendung von kontinuierlichen Wahrscheinlichkeitsverteilungen, anstelle von nur endlich vielen verschiedenen Potentialen, eine weitere Herausforderung dar. Das Integralfunktional im Grenzwert besteht aus einer homogenisierten Energiedichte und ermöglicht Brüche, jedoch ohne einen Beitrag dieser Brüche zur Energie. Um dieses Ergebnis zu verfeinern, skalieren wir unser Modell neu und betrachten dessen $\Gamma$-Grenzwert, der in Form einer Energie vom Griffith-Typ gegeben ist und aus einem elastischen Teil und einem Sprungbeitrag besteht. In einem weiteren Ansatz untersuchen wir Brüche auf Ebene der diskreten Energien. Mit einer geeigneten Definition des Bruchpunktes im diskreten System definieren wir eine Bruchschwelle, die den Elastizitätsbereich von dem Gebiet mit Brüchen trennt. Von diesem Schwellwert berechnen wir anschließend den punktweisen Grenzwert. Es stellt sich heraus, dass dieser Grenzwert mit dem durch die variationelle $\Gamma$-Konvergenz errechneten übereinstimmt.

Diese Dissertation beschäftigt sich mit EDP-Konvergenz. Dabei handelt es sich um einen Konvergenzbegriff auf dem Gebiet der verallgemeinerten Gradientensysteme und metrischen Gradientensysteme, der geeignet ist für Gradientenflüsse, die von einem kleinen Parameter abhängen. EDP-Konvergenz liefert einen Algorithmus, der es erlaubt in der Energie und dem Dissipationspotenzial zum Grenzwert überzugehen. Es ist die fundamentale Frage evolutionärer Γ-Konvergenz, wie das Limes-Dissipationspotenzial berechnet werden kann. Das Ziel dieser Arbeit ist es aufzuzeigen, dass EDP-Konvergenz das mikro- und das makroskopische Dissipationspotenzial in einer sinnvollen und eindeutigen Art und Weise in Beziehung setzt. Anhand von drei Beispielen wird der Konvergenzbegriff untersucht: die Diffusionsgleichung auf einem dünnen, dreischichtigen Gebiet, die Poröse-Medien-Gleichung mit einer dünnen Membran und ein Modell mit oszillierender Energie. Es wird die Definition von relaxierter EDP-Konvergenz und EDP-Konvergenz mit Kippung motiviert. EDP-Konvergenz basiert auf dem Prinzip, dass es ein Gleichgewicht zwischen Energie und Dissipation gibt – das Energie-Dissipations-Prinzip (EDP). Mittels Γ-Konvergenz wird sowohl in der Energie, als auch dem totalen Dissipationsfunktional zum Grenzwert übergegangen. Durch die zusätzliche Entkopplung von Zustand und Triebkraft wird die Dissipationslandschaft erkundet und die kinetische Beziehung des Limessystems ermittelt. Das Modell mit oszillierender Energie zeigt die Bedeutung der kinetischen Beziehung – und damit der Kippung – für die Herleitung des Limes-Dissipationspotenzials auf. Die Modelle mit Wasserstein-Dissipation zeigen, dass das Limes-Dissipationspotenzial nicht der naive Grenzwert ist. Insbesondere können klassische Gradientensysteme mit quadratischer Dissipation zu verallgemeinerten Gradientensysteme konvergieren. / In the realm of generalized gradient systems and metric gradient systems we study a notion of convergence suited for gradient flows which depend on a small parameter. This notion is called EDP-convergence. In order to understand the convergence of gradient systems we need an algorithm to derive the limiting energy as well as the limiting dissipation potential. The fundamental question of evolutionary Γ-convergence is how to compute the limit dissipation potential. The aim of this thesis is to show that EDP-convergence connects the microscopic dissipation potential with the macroscopic, i.e. limiting, dissipation potential in a meaningful and unique way. As a proof of concept 3 different examples are presented: (i) the diffusion equation on a thin sandwich-like domain, (ii) the porous medium equation with a thin interface and (iii) a wiggly energy model. We show how the gradient flow concept that is used in this thesis can be used to obtain also gradient flows with respect to the Wasserstein metric. We motivate the definition of relaxed EDP-convergence and EDP- convergence with tilting. EDP-convergence is based upon the principle that there is an energy-dissipation-balance involving the total dissipation functional and the energy difference – the energy-dissipation-principle (EDP). The limit passage, in both the energy and the total dissipation functional, is performed in terms of Γ-convergence. By perturbing the flow as well as the driving force, the dissipation-landscape is explored and a kinetic relation for the limit system can be established. The wiggly energy model demonstrates the importance of the kinetic relation for the construction of the limiting dissipation potential and thus the introduction of tilts. The models with a Wasserstein dissipation show that the limiting dissipation potential is not the naive limit. In particular, classical gradient systems with a quadratic dissipation potential converge to a generalized gradient systems.

Die vorliegende Arbeit beschäftigt sich mit einer Klasse von Variationsproblemen, die im Kontext des Ferromagnetismus entstehen. Es soll hierbei ein numerischer und analytischer Hintergrund zur Behandlung von harten magnetischen dünnen Filmen mit senkrechter Anisotropie gegeben werden. Bei magnetischen dünnen Filmen handelt es sich um Schichten von magnetischen Materialien mit Dicken von wenigen Mikrometern bis hin zu einigen Nanometern. Ausgangspunkt der Betrachtungen ist ein Modell von Landau und Lifshitz, das die Grundzustände der Magnetisierung in einem dreidimensionalen Körpers mit den Minimierer eines nichtkonvexen und nichtlokalen Energiefunktionals, der sogenannten mikromagnetischen Energie, verbindet. Unter der Annahme sehr kleiner Filmdicken wird aus dem betrachteten Modell ein zwei-dimensionales Modell hergeleitet. Anschließend wird mit Hilfe der Gamma-Konvergenz die Konvergenz zu einem Sharp-Interface-Modell gezeigt. Das resultierende Energiefunktional besteht aus konkurrierenden Interface- und Dipolenergieanteilen. Der zweite Teil der Arbeit beschäftigt sich mit der Analyse einer numerischen Methode, die die Lösungen des vorher hergeleiteten Modells approximiert. Hierbei stützen sich die Betrachtungen auf ein relaxiertes Modell, in dem der Interfaceenergiebeitrag durch seine Modica-Mortola Approximation ersetzt und dann der entsprechende L^2 Gradientenfluß betrachtet wird. Die daraus resultierende nichtlineare und nichtlokale parabolische Gleichung wird anschließend durch ein Crank-Nicolson-Verfahren in der Zeitvariablen und einem Fourieransatz für die Raumvariablen diskretisiert. Wir beweisen die Existenz und Eindeutigkeit von Lösungen des numerischen Verfahrens, sowie deren Konvergenz zu Lösungen des anfänglich betrachteten stetigen Modells. Ferner werden auch a priori Fehlerabschätzungen für die numerische Methode hergeleitet. Abschließend werden die analytischen Resultate anhand numerischer Experimente illustriert. / This thesis is concerned with the study of a class of variational problems arising in the context of ferromagnetism. More precisely, it aims at providing a numerical and analytical background to the study of hard magnetic thin films with perpendicular anisotropy. Magnetic thin films are sheets of magnetic materials with thicknesses of a few micrometers down to a few nanometers used mainly in electronic industry, for example as magnetic data storage media for computers. Our initial considerations are based on a model of Landau and Lifshitz that associates the ground states of the magnetization within a three-dimensional body to the minimizers of a nonconvex and nonlocal energy functional, the so-called micromagnetic energy. Under film thickness considerations (thin film regime), we first reduce the aforementioned model to two dimensions and then carry out a Gamma-limit for a sharp-interface model. The resulting energy functional features a competition between an interfacial and a dipolar energy contribution. The second part of the thesis is concerned with the analysis of a numerical method to approximate solutions of the previously derived sharp-interface model. We base our considerations on a relaxed model in which we replace the interfacial energy contribution by its Modica-Mortola approximation, and then study the associated L^2 gradient flow. The resulting evolution equation, a nonlinear and nonlocal parabolic equation, is discretized by a Crank-Nicolson approximation for the time variable and a Fourier collocation method for the space variable. We prove the existence and uniqueness of the solutions of the numerical scheme, the convergence of these solutions towards solutions of the initial continuous model and also derive a-priori error estimates for the numerical method. Finally, we illustrate the analytical results by a series of numerical experiments.

We consider thin bodies made from elastomers and nematic liquid crystal elastomers. Starting from a nonlinear 3d hyperelastic model, and using the Gamma-convergence method, we derive lower dimensional models for 2d and 1d. The limit models describe the interplay between free liquid crystal orientations and bending deformations.:1 Introduction 1.1 Main results and structure of the text 1.2 Survey of the literature 1.2.1 Dimension reduction in nonlinear elasticity 1.2.2 Relation to other bending regime results in detail 1.2.3 Relation to other Gamma-convergence results of LCEs 2 Liquid crystal elastomers 2.1 Properties 2.2 Modeling 3 Rods 3.1 Setup and statement of analytical main results 3.1.1 The 3d-model and assumptions 3.1.2 The effective 1d-model 3.1.3 The Gamma-convergence result without boundary conditions 3.1.4 Boundary conditions for y 3.1.5 Weak and strong anchoring of n 3.1.6 Definition and properties of the effective coefficients 3.2 Numerical 1d-model exploration 3.3 Dimensional analysis and scalings 3.3.1 Non-dimensionalization and rescaling 3.3.2 Scaling assumptions 3.3.3 Dimensional analysis and applicability of the 1d-model 3.4 Smooth approximation of framed curves 3.5 Proofs 3.5.1 Compactness: proofs of Theorem 3.1.3 (a) and Proposition 3.1.4 (a) 3.5.2 Lower bound: proof of Theorem 3.1.3 (b) . . . . . . . . . . . . 68 3.5.3 Upper bound: proofs of Theorem 3.1.3 (c) and Proposition 3.1.4 (b) 3.5.4 Anchoring: proof of Proposition 3.1.5 3.5.5 Properties of the effective coefficients 4 Plates 4.1 Setup and statement of analytical main results 4.1.1 The 3d-model and assumptions 4.1.2 The effective 2d-model 4.1.3 The Gamma-convergence result without boundary conditions 4.1.4 Definition and properties of the effective coefficients 4.1.5 Boundary conditions for y 4.1.6 Weak and strong anchoring of n 4.2 Analytical and numerical 2d-model exploration 4.2.1 Analytical 2d-model exploration 4.2.2 Numerical 2d-model exploration 4.3 Dimensional analysis and scalings 4.3.1 Non-dimensionalization and rescaling 4.3.2 Scaling assumptions 4.3.3 Dimensional analysis and applicability 4.4 Geometry and approximation of bending deformations 4.4.1 Proofs of the geometric properties in the smooth case 4.4.2 Proof for the smooth approximations 4.5 Proofs 4.5.1 Compactness: proofs of Theorems 4.1.1 (a) and 4.1.8 (a) 4.5.2 Lower bound: proof of Theorem 4.1.1 (b) 4.5.3 Upper bound: proofs of Theorem 4.1.1 (c) and Theorem 4.1.8 (b) 4.5.4 Properties of the effective coefficients 4.5.5 Anchorings 4.5.6 Approximation of nonlinear strains: proof of Proposition 4.5.4 5 Conclusions and outlooks Bibliography

Der Grenzwert kleiner Deformationen in Anwesenheit eines gegebenen Risses wird in drei verschiedenen kontinuumsmechanischen Modellen betrachtet. Erstens wird für rein statische Elastizität mit finiter Spannung im Grenzwert kleiner Belastung bewiesen, dass die Nebenbedingung globaler Injektivität im Sinne der Gamma-Konvergenz eine lokale Nichtdurchdringungsbedingung auf dem Riss ergibt. Zweitens wird Deformationsplastizität mit finiten Spannungen und multiplikativer Zerlegung des Spannungstensors behandelt und die Gamma-Konvergenz zu linearisierter Deformationsplastizität mit Rissbedingungen gezeigt. Drittens wird die ratenunabhängige Evolution der Elastoplastizität betrachtet mit einer allgemeineren Klasse globaler Injektivitätsbedingungen für den finiten Fall. Hierbei wird einerseits die evolutionäre Gamma-Konvergenz unter Vernachlässigung der Nebenbedinung gezeigt, andererseits eine Vermutung aufgestellt, unter deren Voraussetzung die evolutionäre Gamma-Konvergenz auch mit Rissbedingungen gilt. / The small-deformation limit in presence of a given crack is considered in three distinct continuummechanical models. First, a purely static finite-strain elasticity model is considered in the limit of small loading, where the constraint of global injectivity is shown to converge in the sense of Gamma-convergence to a local constraint of non-interpenetration along the crack. Second, finitestrain deformation plasticity based on the multiplicative decomposition of the strain tensor is shown to Gamma-converge to linearized deformation elastoplasticity with crack conditions. Third, the rate-independent evolution of elastoplasticity is considered with a generalized class of global injectivity constraints for the finite-strain model. On the one hand, neglecting the constraints the evolutionary Gamma-converge to linearized elastoplasticity is proven. On the other hand, a conjecture is made, subject to which the evolutionary Gamma-convergence with constraints still holds.

Diese Arbeit beschäftigt sich mit Coarse-Graining (dt. ``Vergröberung", ``Zusammenfassung von Zuständen") für Gradientensysteme und Markov-Prozesse. Coarse-Graining ist ein etabliertes Verfahren in der Mathematik und in den Naturwissenschaften und hat das Ziel, die Komplexität eines physikalischen Systems zu reduzieren und effektive Modelle herzuleiten. Die mathematischen Probleme in dieser Arbeit stammen aus der Theorie der Systeme interagierender Teilchen. Hierbei werden zwei Ziele verfolgt: Erstens, Coarse-Graining mathematisch rigoros zu beweisen, zweitens, mathematisch äquivalente Beschreibungen für die effektiven Modelle zu formulieren. Die ersten drei Teile der Arbeit befassen sich mit dem Grenzwert schneller Reaktionen für Reaktionssysteme und Reaktions-Diffusions-Systeme. Um effektive Modelle herzuleiten, werden nicht nur die zugehörigen Reaktionsratengleichungen betrachtet, sondern auch die zugrunde liegende Gradientenstruktur. Für Gradientensysteme wurde in den letzten Jahren eine strukturelle Konvergenz, die sogenannte ``EDP-Konvergenz", entwickelt. Dieses Coarse-Graining-Verfahren wird in der vorliegenden Arbeit auf folgende Systeme mit langsamen und schnellen Reaktionen angewandt: lineare Reaktionssysteme (bzw. Markov-Prozesse auf endlichem Zustandsraum), nichtlineare Reaktionssysteme, die das Massenwirkungsgesetz erfüllen, und lineare Reaktions-Diffusions-Systeme. Für den Grenzwert schneller Reaktionen wird eine mathematisch rigorose und strukturerhaltende Vergröberung auf dem Level des Gradientensystems inform von EDP-Konvergenz bewiesen. Im vierten Teil wird der Zusammenhang zwischen Gleichungen mit Gedächtnis und Markov-Prozessen untersucht. Für Gleichungen mit Gedächtnisintegralen wird explizit ein größer Markov-Prozess konstruiert, der die Gleichung mit Gedächtnis als Teilsystem enthält. Der letzte Teil beschäftigt sich mit verschieden Diskretisierungen für den Fokker-Planck-Operator. Dazu werden numerische und analytische Eigenschaften untersucht. / This thesis deals with coarse-graining for gradient systems and Markov processes. Coarse-graining is a well-established tool in mathematical and natural sciences for reducing the complexity of a physical system and for deriving effective models. The mathematical problems in this work originate from interacting particle systems. The aim is twofold: first, providing mathematically rigorous results for physical coarse-graining, and secondly, formulating mathematically equivalent descriptions for the effective models. The first three parts of the thesis deal with fast-reaction limits for reaction systems and reaction-diffusion systems. Instead of deriving effective models by solely investigating the associated reaction-rate equation, we derive effective models using the underlying gradient structure of the evolution equation. For gradient systems a structural convergence, the so-called ``EDP-convergence", has been derived in recent years. In this thesis, this coarse-graining procedure has been applied to the following systems with slow and fast reactions: linear reaction systems (or Markov process on finite state space), nonlinear reaction systems of mass-action type, and linear reaction-diffusion systems. For the fast-reaction limit, we perform rigorous and structural coarse-graining on the level of the gradient system by proving EDP-convergence. In the fourth part, the connection between memory equations and Markov processes is investigated. Considering linear memory equations, which can be motivated from spatial homogenization, we explicitly construct a larger Markov process that includes the memory equation as a subsystem. The last part deals with different discretization schemes for the Fokker–Planck operator and investigates their analytical and numerical properties.