Coordinate and time-dependent diffusion dynamics in protein folding

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Particular solution for anomalous diffusion equation with source term

In literature the phenomenon of diffusion has been widely studied, however for nonextensive systems which are governed by a nonlinear stochastic dynamic, there are a few soluble models. The purpose of this study is to present the solution of the nonlinear Fokker-Planck equation for a model of potential with barrier considering a term of absorption. Systems of this nature can be observed in various chemical or biological processes and their solution enriches the studies of existing nonextensive systems.

Comparison of SRIM, MCNPX and GEANT simulations with experimental data for thick Al absorbers

Proton computerized tomography deals with relatively thick targets like the human head or trunk. In this case precise analytical calculation of the proton final energy is a rather complicated task, thus the Monte Carlo simulation stands out as a solution. We used the GEANT4.8.2 code to calculate the proton final energy spectra after passing a thick Al absorber and compared it with the same conditions of the experimental data. The ICRU49, Ziegler85 and Ziegler2000 models from the low energy extension pack were used. The results were also compared with the SRIM2008 and MCNPX2.4 simulations, and with solutions of the Boltzmann transport equation in the Fokker-Planck approximation. (C) 2009 Elsevier Ltd. All rights reserved.

Relating network rigidity, time scale hierarchies, and expression noise in gene networks

Fluctuation-dissipation theorems can be used to predict characteristics of noise from characteristics of the macroscopic response of a system. In the case of gene networks, feedback control determines the "network rigidity," defined as resistance to slow external changes. We propose an effective Fokker-Planck equation that relates gene expression noise to topology and to time scales of the gene network. We distinguish between two situations referred to as normal and inverted time hierarchies. The noise can be buffered by network feedback in the first situation, whereas it can be topology independent in the latter.

Robustness of the non-Markovian Alzheimer walk under stochastic perturbation

The elephant walk model originally proposed by Schutz and Trimper to investigate non-Markovian processes led to the investigation of a series of other random-walk models. Of these, the best known is the Alzheimer walk model, because it was the first model shown to have amnestically induced persistence-i.e. superdiffusion caused by loss of memory. Here we study the robustness of the Alzheimer walk by adding a memoryless stochastic perturbation. Surprisingly, the solution of the perturbed model can be formally reduced to the solutions of the unperturbed model. Specifically, we give an exact solution of the perturbed model by finding a surjective mapping to the unperturbed model. Copyright (C) EPLA, 2012

Non-Markovian stochastic processes and their applications: from anomalous diffusion to time series analysis

This work provides a forward step in the study and comprehension of the relationships between stochastic processes and a certain class of integral-partial differential equation, which can be used in order to model anomalous diffusion and transport in statistical physics. In the first part, we brought the reader through the fundamental notions of probability and stochastic processes, stochastic integration and stochastic differential equations as well. In particular, within the study of H-sssi processes, we focused on fractional Brownian motion (fBm) and its discrete-time increment process, the fractional Gaussian noise (fGn), which provide examples of non-Markovian Gaussian processes. The fGn, together with stationary FARIMA processes, is widely used in the modeling and estimation of long-memory, or long-range dependence (LRD). Time series manifesting long-range dependence, are often observed in nature especially in physics, meteorology, climatology, but also in hydrology, geophysics, economy and many others. We deepely studied LRD, giving many real data examples, providing statistical analysis and introducing parametric methods of estimation. Then, we introduced the theory of fractional integrals and derivatives, which indeed turns out to be very appropriate for studying and modeling systems with long-memory properties. After having introduced the basics concepts, we provided many examples and applications. For instance, we investigated the relaxation equation with distributed order time-fractional derivatives, which describes models characterized by a strong memory component and can be used to model relaxation in complex systems, which deviates from the classical exponential Debye pattern. Then, we focused in the study of generalizations of the standard diffusion equation, by passing through the preliminary study of the fractional forward drift equation. Such generalizations have been obtained by using fractional integrals and derivatives of distributed orders. In order to find a connection between the anomalous diffusion described by these equations and the long-range dependence, we introduced and studied the generalized grey Brownian motion (ggBm), which is actually a parametric class of H-sssi processes, which have indeed marginal probability density function evolving in time according to a partial integro-differential equation of fractional type. The ggBm is of course Non-Markovian. All around the work, we have remarked many times that, starting from a master equation of a probability density function f(x,t), it is always possible to define an equivalence class of stochastic processes with the same marginal density function f(x,t). All these processes provide suitable stochastic models for the starting equation. Studying the ggBm, we just focused on a subclass made up of processes with stationary increments. The ggBm has been defined canonically in the so called grey noise space. However, we have been able to provide a characterization notwithstanding the underline probability space. We also pointed out that that the generalized grey Brownian motion is a direct generalization of a Gaussian process and in particular it generalizes Brownain motion and fractional Brownain motion as well. Finally, we introduced and analyzed a more general class of diffusion type equations related to certain non-Markovian stochastic processes. We started from the forward drift equation, which have been made non-local in time by the introduction of a suitable chosen memory kernel K(t). The resulting non-Markovian equation has been interpreted in a natural way as the evolution equation of the marginal density function of a random time process l(t). We then consider the subordinated process Y(t)=X(l(t)) where X(t) is a Markovian diffusion. The corresponding time-evolution of the marginal density function of Y(t) is governed by a non-Markovian Fokker-Planck equation which involves the same memory kernel K(t). We developed several applications and derived the exact solutions. Moreover, we considered different stochastic models for the given equations, providing path simulations.

Beiträge zur Modellierung technischer dynamischer Systeme

Die Nichtlineare Dynamik verallgemeinert Aussagen über dynamische Systeme durch Abstraktion von konkreten Systemen. In der Technik sind Maschinen dagegen sehr konkret und die Behandlung auftretender Probleme mit Methoden der theoretischen Physik ist nicht trivial. Diese Arbeit versucht einige der Schwierigkeiten einer technischen Anwendung der nichtlinearen Theorie zu lokalisieren. Am Beispiel von vier Klassen von Modellansätzen, werden Anwendungsschnittstellen beleuchtet und systematisiert. Die Anwendung von Modellen, die explizit auf bekannten physikalischen Gesetzmäßigkeiten aufbauen, findet Grenzen in der Anzahl der Freiheitsgrade und den Nebenbedingungen konkreter Systeme. Solche Modelle liefern jedoch wichtige Hinweise auf die Vielfalt der nichtlinearen Phänomene und tragen zu ihrem Verständnis bei. Daher sind sie für die Konstruktionspraxis wichtig. Es werden typisch nichtlineare Phänomene und ihre zugrundeliegenden Mechanismen vorgestellt und klassifiziert, sowie grundsätzliche Probleme der Berechenbarkeit analytisch formulierter Modelle betrachtet. Eine zweite Schnittstelle bieten die Darstellungen des Systemverhaltens als überlagerung spezieller Funktionen, diez.B. Symmetrieeigenschaften des betrachteten Systems besonders deutlich widerspiegeln. Gegenüber der klassischen Fourierzerlegung nach Frequenz und Phase bringt die Analyse nach Detaillierungsgrad und Position von Waveletfunktionen wichtige Vorteile für die nichtlineare zustandsraumbasierte Datenanalyse. Viele Verfahren der Nichtlinearen Datenanalyse beruhen auf metrischen Eigenschaften der dynamischen Systeme. Als dritte Gruppe werden demgegenüber topologische Methoden beleuchtet. Die Konstruktion von Simplexen aus Zeitreihen mittels der Zeitversatzmethode ist die Grundlage für eine Triangulation der Zustandsräume. Die Methoden, z.B. Templateverfahren, die auf der Einbettung von eindimensionalen Trajektorien in den R^3 basieren, lassen sich hingegen nicht einfach auf hochdimensionale Zustandsmannigfaltigkeiten anwenden. Schließlich werden stochastische Aspekte behandelt. Schwankungen des Systemverhaltens können auf Schwankungen der Anfangswerte und/oder auf Schwankungen der eigentlichen Systemdynamik beruhen. Die Einordnung des konkreten Anwendungsfalles setzt jedoch ein sicheres Verständnis stochastischer Prozesse voraus. Am Beispiel der Rekonstruktion der stochastischen Dynamik über eine eindimensionale Fokker-Planck-Gleichung zeigen sich deutlich die praktischen Grenzen solcher Ansätze.

Vortex dynamics in superconductors in the presence of anisotropic pinning

Vortex dynamics in two different classes of superconductors with anisotropic unidirected pinning sites was experimentally investigated by magnetoresistivity measurements: YBCO−films with unidirected twins and Nb-films deposited on faceted $mathrm Al_2O_3$ substrate surfaces. For the interpretation of the experimental results a theoretical model based on the Fokker-Planck equation was used. It was proved by X-ray measurements that YBCO films prepared on (001) $mathrm NdGaO_3$ substrates exhibit only one twin orientation in contrast to YBCO films grown on (100) $mathrm SrTiO_$3 substrates. The magnetoresistivity measurements of the YBCO films with unidirected twin boundaries revealed the existence of two new magnetoresistivity components, which is a characteristic feature of a guided vortex motion: an odd longitudinal component with respect to the magnetic field sign reversal and an even transversal component. However, due to the small coherence length in YBCO and the higher density of point-like defects comparing to high-quality YBCO single crystals, the strength of the isotropic point pinning was comparable with the strength of the pinning produced by twins. This smeared out all effects caused by the pinning anisotropy. The behaviour of the odd longitudinal component was found to be independent of the transport current direction with respect to the twin planes. The magnetoresistivity measurements of faceted Nb films demonstrated the appearance of an odd longitudinal and even transversal component of the magnetoresistivity. The temperature and magnetic field dependences of all relevant magnetoresistivity components were measured. The angles between the average vortex velocity vector and the transport current direction calculated from the experimental data for the different transport current orientations with respect to the facet ridges showed that the vortices moved indeed along the facet ridges. An anomalous Hall effect, i.e. a sign change of the odd transversal magnetoresistivity, has been found in the temperature and magnetic field dependences of the Hall resisitivity of the samples. The theory developed by V.~A.~Shklovskij was used for the explanation of the experimental data. It shows very good agreement with the experiment. The temperature dependence of the even longitudinal magnetoresistivity component of the samples could be very well fitted within the theoretical approach, using for the isotropic and anisotropic pinning potential simple potential with a symmetric triangular potential wells whose depths were estimated from the experimental data.

Analysis of nonlinear diffusion equations of second and fourth order

Wegen der fortschreitenden Miniaturisierung von Halbleiterbauteilen spielen Quanteneffekte eine immer wichtigere Rolle. Quantenphänomene werden gewöhnlich durch kinetische Gleichungen beschrieben, aber manchmal hat eine fluid-dynamische Beschreibung Vorteile: die bessere Nutzbarkeit für numerische Simulationen und die einfachere Vorgabe von Randbedingungen. In dieser Arbeit werden drei Diffusionsgleichungen zweiter und vierter Ordnung untersucht. Der erste Teil behandelt die implizite Zeitdiskretisierung und das Langzeitverhalten einer degenerierten Fokker-Planck-Gleichung. Der zweite Teil der Arbeit besteht aus der Untersuchung des viskosen Quantenhydrodynamischen Modells in einer Raumdimension und dessen Langzeitverhaltens. Im letzten Teil wird die Existenz von Lösungen einer parabolischen Gleichung vierter Ordnung in einer Raumdimension bewiesen, und deren Langzeitverhalten studiert.

Quantum hydrodynamic models for semiconductors with and without collisions

In dieser Arbeit werden Quantum-Hydrodynamische (QHD) Modelle betrachtet, die ihren Einsatz besonders in der Modellierung von Halbleiterbauteilen finden. Das QHD Modell besteht aus den Erhaltungsgleichungen für die Teilchendichte, das Momentum und die Energiedichte, inklusive der Quanten-Korrekturen durch das Bohmsche Potential. Zu Beginn wird eine Übersicht über die bekannten Ergebnisse der QHD Modelle unter Vernachlässigung von Kollisionseffekten gegeben, die aus ein­em Schrödinger-System für den gemischten-Zustand oder aus der Wigner-Glei­chung hergeleitet werden können. Nach der Reformulierung der eindimensionalen QHD Gleichungen mit linearem Potential als stationäre Schrö­din­ger-Gleichung werden die semianalytischen Fassungen der QHD Gleichungen für die Gleichspannungs-Kurve betrachtet. Weiterhin werden die viskosen Stabilisierungen des QHD Modells be­rück­sich­tigt, sowie die von Gardner vorgeschlagene numerische Viskosität für das {sf upwind} Finite-Differenzen Schema berechnet. Im Weiteren wird das viskose QHD Modell aus der Wigner-Glei­chung mit Fokker-Planck Kollisions-Ope­ra­tor hergeleitet. Dieses Modell enthält die physikalische Viskosität, die durch den Kollision-Operator eingeführt wird. Die Existenz der Lösungen (mit strikt positiver Teilchendichte) für das isotherme, stationäre, eindimensionale, viskose Modell für allgemeine Daten und nichthomogene Randbedingungen wird gezeigt. Die dafür notwendigen Abschätzungen hängen von der Viskosität ab und erlauben daher den Grenzübergang zum nicht-viskosen Fall nicht. Numerische Simulationen der Resonanz-Tunneldiode modelliert mit dem nichtisothermen, stationären, eindimensionalen, viskosen QHD Modell zeigen den Einfluss der Viskosität auf die Lösung. Unter Verwendung des von Degond und Ringhofer entwickelten Quanten-Entropie-Minimierungs-Verfahren werden die allgemeinen QHD-Gleichungen aus der Wigner-Boltzmann-Gleichung mit dem BGK-Kollisions-Operator hergeleitet. Die Herleitung basiert auf der vorsichtige Entwicklung des Quanten-Max­well­ians in Potenzen der skalierten Plankschen Konstante. Das so erhaltene Modell enthält auch vertex-Terme und dispersive Terme für die Ge­schwin­dig­keit. Dadurch bleibt die Gleichspannungs-Kurve für die Re­so­nanz-Tunnel­diode unter Verwendung des allgemeinen QHD Modells in einer Dimension numerisch erhalten. Die Ergebnisse zeigen, dass der dispersive Ge­schwin­dig­keits-Term die Lösung des Systems stabilisiert.

Single polymers at surfaces: adsorption and detachment

This thesis is concerned with the adsorption and detachment of polymers at planar, rigid surfaces. We have carried out a systematic investigation of adsorption of polymers using analytical techniques as well as Monte Carlo simulations with a coarse grained off-lattice bead spring model. The investigation was carried out in three stages. In the first stage the adsorption of a single multiblock AB copolymer on a solid surface was investigated by means of simulations and scaling analysis. It was shown that the problem could be mapped onto an effective homopolymer problem. Our main result was the phase diagram of regular multiblock copolymers which shows an increase in the critical adsorption potential of the substrate with decreasing size of blocks. We also considered the adsorption of random copolymers which was found to be well described within the annealed disorder approximation. In the next phase, we studied the adsorption kinetics of a single polymer on a flat, structureless surface in the regime of strong physisorption. The idea of a ’stem-flower’ polymer conformation and the mechanism of ’zipping’ during the adsorption process were used to derive a Fokker-Planck equation with reflecting boundary conditions for the time dependent probability distribution function (PDF) of the number of adsorbed monomers. The numerical solution of the time-dependent PDF obtained from a discrete set of coupled differential equations were shown to be in perfect agreement with Monte Carlo simulation results. Finally we studied force induced desorption of a polymer chain adsorbed on an attractive surface. We approached the problem within the framework of two different statistical ensembles; (i) by keeping the pulling force fixed while measuring the position of the polymer chain end, and (ii) by measuring the force necessary to keep the chain end at fixed distance above the adsorbing plane. In the first case we treated the problem within the framework of the Grand Canonical Ensemble approach and derived analytic expressions for the various conformational building blocks, characterizing the structure of an adsorbed linear polymer chain, subject to pulling force of fixed strength. The main result was the phase diagram of a polymer chain under pulling. We demonstrated a novel first order phase transformation which is dichotomic i.e. phase coexistence is not possible. In the second case, we carried out our study in the “fixed height” statistical ensemble where one measures the fluctuating force, exerted by the chain on the last monomer when a chain end is kept fixed at height h over the solid plane at different adsorption strength ε. The phase diagram in the h − ε plane was calculated both analytically and by Monte Carlo simulations. We demonstrated that in the vicinity of the polymer desorption transition a number of properties like fluctuations and probability distribution of various quantities behave differently, if h rather than the force, f, is used as an independent control parameter.

Complex chemical dynamics through engineering-like methods

Most of the problems in modern structural design can be described with a set of equation; solutions of these mathematical models can lead the engineer and designer to get info during the design stage. The same holds true for physical-chemistry; this branch of chemistry uses mathematics and physics in order to explain real chemical phenomena. In this work two extremely different chemical processes will be studied; the dynamic of an artificial molecular motor and the generation and propagation of the nervous signals between excitable cells and tissues like neurons and axons. These two processes, in spite of their chemical and physical differences, can be both described successfully by partial differential equations, that are, respectively the Fokker-Planck equation and the Hodgkin and Huxley model. With the aid of an advanced engineering software these two processes have been modeled and simulated in order to extract a lot of physical informations about them and to predict a lot of properties that can be, in future, extremely useful during the design stage of both molecular motors and devices which rely their actions on the nervous communications between active fibres.

Analisi di un modello matematico unificato per lo studio del processo angiogenico

In questa tesi si studia l'angiogenesi tumorale, dapprima descrivendo i fenomeni biologici alla base della dinamica cellulare, e successivamente, dopo aver introdotto gli strumenti matematici necessari, sviluppandone un modello seguendo la letteratura esistente basato sulle equazioni differenziali stocastiche e su quelle di Fokker-Planck. Ne vengono infine realizzate simulazioni numeriche.

Un Modello Stocastico per la Teoria dell'Evoluzione

Il testo contiene nozioni base di probabilità necessarie per introdurre i processi stocastici. Sono trattati infatti nel secondo capitolo i processi Gaussiani, di Markov e di Wiener, l'integrazione stocastica alla Ito, e le equazioni differenziali stocastiche. Nel terzo capitolo viene introdotto il rapporto tra la genetica e la matematica, dove si introduce l'evoluzione la selezione naturale, e altri fattori che portano al cambiamento di una popolazione; vengono anche formulate le leggi basilari per una modellizzazione dell’evoluzione fenotipica. Successivamente si entra più nel dettaglio, e si determina un modello stocastico per le mutazioni, cioè un modello che riesca ad approssimare gli effetti dei fattori di fluttuazione all'interno del processo evolutivo.

Non-Maxwellian electron distributions in time-dependent simulations of low-Z materials illuminated by a high-intensity X-ray laser

The interaction of high intensity X-ray lasers with matter is modeled. A collisional-radiative timedependent module is implemented to study radiation transport in matter from ultrashort and ultraintense X-ray bursts. Inverse bremsstrahlung absorption by free electrons, electron conduction or hydrodynamic effects are not considered. The collisional-radiative system is coupled with the electron distribution evolution treated with a Fokker-Planck approach with additional inelastic terms. The model includes spontaneous emission, resonant photoabsorption, collisional excitation and de-excitation, radiative recombination, photoionization, collisional ionization, three-body recombination, autoionization and dielectronic capture. It is found that for high densities, but still below solid, collisions play an important role and thermalization times are not short enough to ensure a thermal electron distribution. At these densities Maxwellian and non-Maxwellian electron distribution models yield substantial differences in collisional rates, modifying the atomic population dynamics.