Master Equation: Biological Applications and Thermodynamic Description

It is well known that many realistic mathematical models of biological systems, such as cell growth, cellular development and differentiation, gene expression, gene regulatory networks, enzyme cascades, synaptic plasticity, aging and population growth need to include stochasticity. These systems are not isolated, but rather subject to intrinsic and extrinsic fluctuations, which leads to a quasi equilibrium state (homeostasis). The natural framework is provided by Markov processes and the Master equation (ME) describes the temporal evolution of the probability of each state, specified by the number of units of each species. The ME is a relevant tool for modeling realistic biological systems and allow also to explore the behavior of open systems. These systems may exhibit not only the classical thermodynamic equilibrium states but also the nonequilibrium steady states (NESS). This thesis deals with biological problems that can be treat with the Master equation and also with its thermodynamic consequences. It is organized into six chapters with four new scientific works, which are grouped in two parts: (1) Biological applications of the Master equation: deals with the stochastic properties of a toggle switch, involving a protein compound and a miRNA cluster, known to control the eukaryotic cell cycle and possibly involved in oncogenesis and with the propose of a one parameter family of master equations for the evolution of a population having the logistic equation as mean field limit. (2) Nonequilibrium thermodynamics in terms of the Master equation: where we study the dynamical role of chemical fluxes that characterize the NESS of a chemical network and we propose a one parameter parametrization of BCM learning, that was originally proposed to describe plasticity processes, to study the differences between systems in DB and NESS.

Time-lapse pressure tomography for characterizing CO2 plume evolution in a deep saline aquifer

A time-lapse pressure tomography inversion approach is applied to characterize the CO2 plume development in a virtual deep saline aquifer. Deep CO2 injection leads to flow properties of the mixed-phase, which vary depending on the CO2 saturation. Analogous to the crossed ray paths of a seismic tomographic experiment, pressure tomography creates streamline patterns by injecting brine prior to CO2 injection or by injecting small amounts of CO2 into the two-phase (brine and CO2) system at different depths. In a first step, the introduced pressure responses at observation locations are utilized for a computationally rapid and efficient eikonal equation based inversion to reconstruct the heterogeneity of the subsurface with diffusivity (D) tomograms. Information about the plume shape can be derived by comparing D-tomograms of the aquifer at different times. In a second step, the aquifer is subdivided into two zones of constant values of hydraulic conductivity (K) and specific storage (Ss) through a clustering approach. For the CO2 plume, mixed-phase K and Ss values are estimated by minimizing the difference between calculated and “true” pressure responses using a single-phase flow simulator to reduce the computing complexity. Finally, the estimated flow property is converted to gas saturation by a single-phase proxy, which represents an integrated value of the plume. This novel approach is tested first with a doublet well configuration, and it reveals a great potential of pressure tomography based concepts for characterizing and monitoring deep aquifers, as well as the evolution of a CO2 plume. Still, field-testing will be required for better assessing the applicability of this approach.

A large strain anisotropic elastoplastic continuum theory for nonlinear kinematic hardening and texture evolution

In this paper we present a continuum theory for large strain anisotropic elastoplasticity based on a decomposition of the modified plastic velocity gradient into energetic and dissipative parts. The theory includes the Armstrong and Frederick hardening rule as well as multilayer models as special cases even for large strain anisotropic elastoplasticity. Texture evolution may also be modelled by the formulation, which allows for a meaningful interpretation of the terms of the dissipation equation

Decomposition driven interface evolution for layers of binary mixtures: I. Model dirivation and base states

A dynamical model is proposed to describe the coupled decomposition and profile evolution of a free surfacefilm of a binary mixture. An example is a thin film of a polymer blend on a solid substrate undergoing simultaneous phase separation and dewetting. The model is based on model-H describing the coupled transport of the mass of one component (convective Cahn-Hilliard equation) and momentum (Navier-Stokes-Korteweg equations) supplemented by appropriate boundary conditions at the solid substrate and the free surface. General transport equations are derived using phenomenological nonequilibrium thermodynamics for a general nonisothermal setting taking into account Soret and Dufour effects and interfacial viscosity for the internal diffuse interface between the two components. Focusing on an isothermal setting the resulting model is compared to literature results and its base states corresponding to homogeneous or vertically stratified flat layers are analyzed.

Modelling relativistic solitary wave interactions in over-dense plasmas: a perturbed nonlinear Schröndinger equation framework

We investigate the dynamics of localized solutions of the relativistic cold-fluid plasma model in the small but finite amplitude limit, for slightly overcritical plasma density. Adopting a multiple scale analysis, we derive a perturbed nonlinear Schrödinger equation that describes the evolution of the envelope of circularly polarized electromagnetic field. Retaining terms up to fifth order in the small perturbation parameter, we derive a self-consistent framework for the description of the plasma response in the presence of localized electromagnetic field. The formalism is applied to standing electromagnetic soliton interactions and the results are validated by simulations of the full cold-fluid model. To lowest order, a cubic nonlinear Schrödinger equation with a focusing nonlinearity is recovered. Classical quasiparticle theory is used to obtain analytical estimates for the collision time and minimum distance of approach between solitons. For larger soliton amplitudes the inclusion of the fifth-order terms is essential for a qualitatively correct description of soliton interactions. The defocusing quintic nonlinearity leads to inelastic soliton collisions, while bound states of solitons do not persist under perturbations in the initial phase or amplitude

Nonlinear analysis of a simple model of temperature evolution in a satellite

We analyze a simple model of the heat transfer to and from a small satellite orbiting round a solar system planet. Our approach considers the satellite isothermal, with external heat input from the environment and from internal energy dissipation, and output to the environment as black-body radiation. The resulting nonlinear ordinary differential equation for the satellite’s temperature is analyzed by qualitative, perturbation and numerical methods, which prove that the temperature approaches a periodic pattern (attracting limit cycle). This approach can occur in two ways, according to the values of the parameters: (i) a slow decay towards the limit cycle over a time longer than the period, or (ii) a fast decay towards the limit cycle over a time shorter than the period. In the first case, an exactly soluble average equation is valid. We discuss the consequences of our model for the thermal stability of satellites.

Self-similar intermediate asymptotics for a degenerate parabolic filtration-absorption equation

The equation ∂tu = u∂xx2u − (c − 1)(∂xu)2 is known in literature as a qualitative mathematical model of some biological phenomena. Here this equation is derived as a model of the groundwater flow in a water-absorbing fissurized porous rock; therefore, we refer to this equation as a filtration-absorption equation. A family of self-similar solutions to this equation is constructed. Numerical investigation of the evolution of non-self-similar solutions to the Cauchy problems having compactly supported initial conditions is performed. Numerical experiments indicate that the self-similar solutions obtained represent intermediate asymptotics of a wider class of solutions when the influence of details of the initial conditions disappears but the solution is still far from the ultimate state: identical zero. An open problem caused by the nonuniqueness of the solution of the Cauchy problem is discussed.

A new code for the Hall-driven magnetic evolution of neutron stars

Over the past decade, the numerical modeling of the magnetic field evolution in astrophysical scenarios has become an increasingly important field. In the crystallized crust of neutron stars the evolution of the magnetic field is governed by the Hall induction equation. In this equation the relative contribution of the two terms (Hall term and Ohmic dissipation) varies depending on the local conditions of temperature and magnetic field strength. This results in the transition from the purely parabolic character of the equations to the hyperbolic regime as the magnetic Reynolds number increases, which presents severe numerical problems. Up to now, most attempts to study this problem were based on spectral methods, but they failed in representing the transition to large magnetic Reynolds numbers. We present a new code based on upwind finite differences techniques that can handle situations with arbitrary low magnetic diffusivity and it is suitable for studying the formation of sharp current sheets during the evolution. The code is thoroughly tested in different limits and used to illustrate the evolution of the crustal magnetic field in a neutron star in some representative cases. Our code, coupled to cooling codes, can be used to perform long-term simulations of the magneto-thermal evolution of neutron stars.

Semisolid microstructural evolution of AlSi7Mg alloy during partial remelting

Different as-cast microstructures of an AlSi7Mg alloy were produced by controlling the solidification conditions. The as-cast grain size ranged from 1.4 mm to 160 mum and the morphology varied from dendritic to rosette-like to globular. The as-cast materials were then partially remelted and isothermally held at 580degreesC for microstructure evolution. The final microstructure depended on the initial as-cast microstructure and the isothermal holding time. After partial remelting and isothermal holding, coarse-grained dendritic structures were not able to evolve to a globular structure, while structures with medium sized dendritic grains evolved to a globular structure with a relatively large particle size after a long isothermal holding time. Fine-grained structures evolved to well-rounded globular grains within times ranging front 10 min to 5 min as the dendritic nature of the starting structure diminished. An empirical equation has been established to describe the relationship between the evolved microstructure and the as-cast microstructure. (C) 2003 Elsevier B.V. All rights reserved.

Time-evolution of highly localized positive-energy states of the free Dirac electron

Highly localized positive-energy states of the free Dirac electron are constructed and shown to evolve in a simple way under the action of Dirac's equation. When the initial uncertainty in position is small on the scale of the Compton wavelength, there is an associated uncertainty in the mean energy that is large compared with the rest mass of the electron. However, this does not lead to any breakdown of the one-particle description, associated with the possibility of pair-production, but rather leads to a rapid expansion of the probability density outwards from the point of localization, at speeds close to the speed of light.

The effect of Cu addition and milling contaminations on the microstructure evolution of ball milled Al-Pb alloy during sintering

Al-10 wt.%Pb and Al-10 wt.%Pb-x wt.%Cu (x = 0-7.0) bulk alloys were prepared by sintering the mechanically alloyed powders at various temperatures. The microstructure changes of the as consolidated powders in the course of sintering were analyzed by differential scanning calorimetry, scanning electron microscopy, X-ray diffraction analysis and transmission electron microscopy. It has been found that, with respect to the Al-10 wt.%Pb-x wt.%Cu alloy, CuAl2 and Cu9Al4 phases formed in the milling process, and the amount of CuAl2 phase increased while the Cu9Al4 phase disappeared gradually in the sintering process. In both Al-10 wt.%Pb and Al-10 wt.%Pb-x wt.%Cu alloys, the sintering process results in the coarsening of Pb phase and the growth rate of Pb phase fulfills the Lifshitz-Slyozov-Wagner equation even though the size of the Pb phase was in nanometer range. The Pb particle exhibits cuboctahedral morphology and has a cubic to cubic orientation relationship with the Al matrix. The addition of Cu strongly depressed the growth rate of Pb. Contamination induced by milling has apparent influence on the microstructure of the sintered alloys. Al7Cu2Fe and aluminium oxide phases were identified in the sintered alloys. The cuboctahedral morphology of Pb particles was broken up by the presence of the oxide phase. (c) 2006 Elsevier B.V. All rights reserved.

Impact of third-order fibre dispersion on the evolution of parabolic optical pulses

We develop a perturbation analysis that describes the effect of third-order dispersion on the similariton pulse solution of the nonlinear Schrodinger equation in a fibre gain medium. The theoretical model predicts with sufficient accuracy the pulse structural changes induced, which are observed through direct numerical simulations.

Impact of third-order dispersion on the evolution of parabolic pulses

We present a perturbation analysis that describes the effect of third-order dispersion on the similariton pulse solution of the nonlinear Schrödinger equation in a fibre gain medium. The theoretical model predicts with sufficient accuracy the pulse structural changes induced, which are observed through direct numerical simulations.

Breakup of a multisoliton state of the linearly damped nonlinear Schrödinger equation

We address the breakup (splitting) of multisoliton solutions of the nonlinear Schrödinger equation (NLSE), occurring due to linear loss. Two different approaches are used for the study of the splitting process. The first one is based on the direct numerical solution of the linearly damped NLSE and the subsequent analysis of the eigenvalue drift for the associated Zakharov-Shabat spectral problem. The second one involves the multisoliton adiabatic perturbation theory applied for studying the evolution of the solution parameters, with the linear loss taken as a small perturbation. We demonstrate that in the case of strong nonadiabatic loss the evolution of the Zakharov-Shabat eigenvalues can be quite nontrivial. We also demonstrate that the multisoliton breakup can be correctly described within the framework of the adiabatic perturbation theory and can take place even due to small linear loss. Eventually we elucidate the occurrence of the splitting and its dependence on the phase mismatch between the solitons forming a two-soliton bound state. © 2007 The American Physical Society.

Impact of third-order dispersion on the evolution of parabolic pulses

We present a perturbation analysis that describes the effect of third-order dispersion on the similariton pulse solution of the nonlinear Schrödinger equation in a fibre gain medium. The theoretical model predicts with sufficient accuracy the pulse structural changes induced, which are observed through direct numerical simulations.