Time-delayed fronts from biased random walks

We generalize a previous model of time-delayed reaction–diffusion fronts (Fort and Méndez 1999 Phys. Rev. Lett. 82 867) to allow for a bias in the microscopic random walk of particles or individuals. We also present a second model which takes the time order of events (diffusion and reproduction) into account. As an example, we apply them to the human invasion front across the USA in the 19th century. The corrections relative to the previous model are substantial. Our results are relevant to physical and biological systems with anisotropic fronts, including particle diffusion in disordered lattices, population invasions, the spread of epidemics, etc

Ion-exchange resins as stationary phase in reactive chromatography

Dynamic behavior of bothisothermal and non-isothermal single-column chromatographic reactors with an ion-exchange resin as the stationary phase was investigated. The reactor performance was interpreted by using results obtained when studying the effect of the resin properties on the equilibrium and kinetic phenomena occurring simultaneously in the reactor. Mathematical models were derived for each phenomenon and combined to simulate the chromatographic reactor. The phenomena studied includes phase equilibria in multicomponent liquid mixture¿ion-exchange resin systems, chemicalequilibrium in the presence of a resin catalyst, diffusion of liquids in gel-type and macroporous resins, and chemical reaction kinetics. Above all, attention was paid to the swelling behavior of the resins and how it affects the kinetic phenomena. Several poly(styrene-co-divinylbenzene) resins with different cross-link densities and internal porosities were used. Esterification of acetic acid with ethanol to produce ethyl acetate and water was used as a model reaction system. Choosing an ion-exchange resin with a low cross-link density is beneficial inthe case of the present reaction system: the amount of ethyl acetate as well the ethyl acetate to water mole ratio in the effluent stream increase with decreasing cross-link density. The enhanced performance of the reactor is mainly attributed to increasing reaction rate, which in turn originates from the phase equilibrium behavior of the system. Also mass transfer considerations favor the use ofresins with low cross-link density. The diffusion coefficients of liquids in the gel-type ion-exchange resins were found to fall rapidly when the extent of swelling became low. Glass transition of the polymer was not found to significantlyretard the diffusion in sulfonated PS¿DVB ion-exchange resins. It was also shown that non-isothermal operation of a chromatographic reactor could be used to significantly enhance the reactor performance. In the case of the exothermic modelreaction system and a near-adiabatic column, a positive thermal wave (higher temperature than in the initial state) was found to travel together with the reactive front. This further increased the conversion of the reactants. Diffusion-induced volume changes of the ion-exchange resins were studied in a flow-through cell. It was shown that describing the swelling and shrinking kinetics of the particles calls for a mass transfer model that explicitly includes the limited expansibility of the polymer network. A good description of the process was obtained by combining the generalized Maxwell-Stefan approach and an activity model that was derived from the thermodynamics of polymer solutions and gels. The swelling pressure in the resin phase was evaluated by using a non-Gaussian expression forthe polymer chain length distribution. Dimensional changes of the resin particles necessitate the use of non-standard mathematical tools for dynamic simulations. A transformed coordinate system, where the mass of the polymer was used as a spatial variable, was applied when simulating the chromatographic reactor columns as well as the swelling and shrinking kinetics of the resin particles. Shrinking of the particles in a column leads to formation of dead volume on top of the resin bed. In ordinary Eulerian coordinates, this results in a moving discontinuity that in turn causes numerical difficulties in the solution of the PDE system. The motion of the discontinuity was eliminated by spanning two calculation grids in the column that overlapped at the top of the resin bed. The reactive and non-reactive phase equilibrium data were correlated with a model derived from thethermodynamics of polymer solution and gels. The thermodynamic approach used inthis work is best suited at high degrees of swelling because the polymer matrixmay be in the glassy state when the extent of swelling is low.

Spatial pattern of landslides in Swiss Rhone valley

The present study analyses the spatial pattern of quaternary gravitational slope deformations (GSD) and historical/present-day instabilities (HPI) inventoried in the Swiss Rhone Valley. The main objective is to test if these events are clustered (spatial attraction) or randomly distributed (spatial independency). Moreover, analogies with the cluster behaviour of earthquakes inventoried in the same area were examined. The Ripley's K-function was applied to measure and test for randomness. This indicator allows describing the spatial pattern of a point process at increasing distance values. To account for the non-constant intensity of the geological phenomena, a modification of the K-function for inhomogeneous point processes was adopted. The specific goal is to explore the spatial attraction (i.e. cluster behaviour) among landslide events and between gravitational slope deformations and earthquakes. To discover if the two classes of instabilities (GSD and HPI) are spatially independently distributed, the cross K-function was computed. The results show that all the geological events under study are spatially clustered at a well-defined distance range. GSD and HPI show a similar pattern distribution with clusters in the range 0.75?9 km. The cross K-function reveals an attraction between the two classes of instabilities in the range 0?4 km confirming that HPI are more prone to occur within large-scale slope deformations. The K-function computed for GSD and earthquakes indicates that both present a cluster tendency in the range 0?10 km, suggesting that earthquakes could represent a potential predisposing factor which could influence the GSD distribution.

Propagation through fractal media: The Sierpinski gasket and the Koch curve

We present new analytical tools able to predict the averaged behavior of fronts spreading through self-similar spatial systems starting from reaction-diffusion equations. The averaged speed for these fronts is predicted and compared with the predictions from a more general equation (proposed in a previous work of ours) and simulations. We focus here on two fractals, the Sierpinski gasket (SG) and the Koch curve (KC), for two reasons, i.e. i) they are widely known structures and ii) they are deterministic fractals, so the analytical study of them turns out to be more intuitive. These structures, despite their simplicity, let us observe several characteristics of fractal fronts. Finally, we discuss the usefulness and limitations of our approa

Description of diffusive and propagative behavior on fractals

The known properties of diffusion on fractals are reviewed in order to give a general outlook of these dynamic processes. After that, we propose a description developed in the context of the intrinsic metric of fractals, which leads us to a differential equation able to describe diffusion in real fractals in the asymptotic regime. We show that our approach has a stronger physical justification than previous works on this field. The most important result we present is the introduction of a dependence on time and space for the conductivity in fractals, which is deduced by scaling arguments and supported by computer simulations. Finally, the diffusion equation is used to introduce the possibility of reaction-diffusion processes on fractals and analyze their properties. Specifically, an analytic expression for the speed of the corresponding travelling fronts, which can be of great interest for application purposes, is derived

Cohabitation effect on the slowdown of the Neolithic expansion

We introduce the effect of cohabitation between generations to a previous model on the slowdown of the Neolithic transition in Europe. This effect consists on the fact that human beings do not leave their children alone when they migrate, but on the contrary they cohabit until their children reach adulthood. We also use archaeological data to estimate the variation of the Mesolithic population density with distance, and use this information to predict the slowdown of the Neolithic front speed. The new equation leads to a substantial correction, up to 37%, relative to previous results. The new model is able to provide a satisfactory explanation not only to the relative speed but also to the absolute speed of the Neolithic front obtained from archaeological data

The Role of the Delay Time in the Modeling of Biological Range Expansions

The time interval between successive migrations of biological species causes a delay time in the reaction-diffusion equations describing their space-time dynamics. This lowers the predicted speed of the waves of advance, as compared to classical models. It has been shown that this delay-time effect improves the modeling of human range expansions. Here, we demonstrate that it can also be important for other species. We present two new examples where the predictions of the time-delayed and the classical (Fisher) approaches are compared to experimental data. No free or adjustable parameters are used. We show that the importance of the delay effect depends on the dimensionless product of the initial growth rate and the delay time. We argue that the delay effect should be taken into account in the modeling of range expansions for biological species

COMPREENSÃO MATEMÁTICA DA DIFUSÃO NO CONTEXTO ELETROQUÍMICO

The aim of this study was to explain in detail the mathematical methods used to deal with diffusion equations, mainly for students and researchers interested in electrochemistry and related areas. Emphasis was placed on the deduction and resolution of diffusion equations, as well as addressing cartesian, spherical and cylindrical coordinates. Different aspects of mass transfer processes were discussed including the importance of the resolution of Fick's laws equations to understand and derive parameters of the electroactive species (e.g., diffusion coefficients, formal electrode potentials) from the electrochemical techniques. As an example, the resolution of diffusion equations for a reversible reduction process of soluble oxidized species was presented for the chronopotentiometry technique. This study is envisaged to broaden the understanding of these frequently used methods, in which mathematical deductions are not always completely understood.

Phonon spectra and thermal properties of some fcc metals using embedded-atom potentials /

By employing the embedded-atom potentials of Mei et ai.[l], we have calculated the dynamical matrices and phonon dispersion curves for six fee metals (Cu,Ag,Au,Ni,Pd and Pt). We have also investigated, within the quasiharmonic approximation, some other thermal properties of these metals which depend on the phonon density of states, such as the temperature dependence of lattice constant, coefficient of linear thermal expansion, isothermal and adiabatic bulk moduli, heat capacities at constant volume and constant pressure, Griineisen parameter and Debye temperature. The computed results are compared with the experimental findings wherever possible. The comparison shows a generally good agreement between the theoretical values and experimental data for all properties except the discrepancies of phonon frequencies and Debye temperature for Pd, Pt and Au. Further, we modify the parameters of this model for Pd and Pt and obtain the phonon dispersion curves which is in good agreement with experimental data.

Modélisation d'arythmies auriculaires modulées par le système nerveux autonome

La fibrillation auriculaire (FA) est la forme d’arythmie la plus fréquente et représente environ un tiers des hospitalisations attribuables aux troubles du rythme cardiaque. Les mécanismes d’initiation et de maintenance de la FA sont complexes et multiples. Parmi ceux-ci, une contribution du système nerveux autonome a été identifiée mais son rôle exact demeure mal compris. Ce travail cible l’étude de la modulation induite par l’acétylcholine (ACh) sur l’initiation et le maintien de la FA, en utilisant un modèle de tissu bidimensionnel. La propagation de l’influx électrique sur ce tissu est décrite par une équation réaction-diffusion non-linéaire résolue sur un maillage rectangulaire avec une méthode de différences finies, et la cinétique d'ACh suit une évolution temporelle prédéfinie qui correspond à l’activation du système parasympathique. Plus de 4400 simulations ont été réalisées sur la base de 4 épisodes d’arythmies, 5 tailles différentes de région modulée par l’ACh, 10 concentrations d’ACh et 22 constantes de temps de libération et de dégradation d’ACh. La complexité de la dynamique des réentrées est décrite en fonction de la constante de temps qui représente le taux de variation d’ACh. Les résultats obtenus suggèrent que la stimulation vagale peut mener soit à une dynamique plus complexe des réentrées soit à l’arrêt de la FA en fonction des quatre paramètres étudiés. Ils démontrent qu’une décharge vagale rapide, représentée par des constantes de temps faibles combinées à une quantité suffisamment grande d’ACh, a une forte probabilité de briser la réentrée primaire provoquant une activité fibrillatoire. Cette activité est caractérisée par la création de plusieurs ondelettes à partir d’un rotor primaire sous l’effet de l’hétérogénéité du gradient de repolarisation causé par l’activité autonomique.

Utilisation du motif imidazole en transport d'anions et en catalyse organométallique

Le motif imidazole, un hétérocycle à 5 atomes contenant 2 atomes d’azote et trois atomes de carbone, présente des propriétés physico-chimiques intéressantes qui en font un composé de choix pour plusieurs applications. Parmi ces propriétés, la fonctionnalisation simple des deux atomes d’azote pour former un sel d’imidazolium est très intéressante. Ces sels sont d’excellents précurseurs de carbènes N-hétérocycliques (NHC) et sont couramment utilisés pour synthétiser des ligands en vue d’une utilisation en catalyse organométallique. D’autre part, cette famille de composés possède des propriétés anionophores permettant une utilisation en transport anionique. Le présent travail contient les résultats de travaux concernant ces deux domaines, soit la catalyse et le transport anionique. Dans un premier temps, les propriétés de dérivés de l’imidazole sont exploitées pour former un catalyseur de type palladium-NHC qui est utilisé pour catalyser la réaction de Suzuki-Miyaura en milieu aqueux. L’efficacité de ce catalyseur a été démontrée en utilisant aussi peu que 0,001 mol% pour un rendement quantitatif. Il s’agit de la première occurrence d’un processus hétérogène et recyclable dans l’eau, utilisant un catalyseur de type Pd-NHC et qui ne nécessite aucun additif ou co-solvant. Le recyclage a été prouvé jusqu’à 10 cycles sans diminution apparente de l’activité du catalyseur. Dans un second temps, plusieurs sels d’imidazolium ont été testés en tant que transporteurs transmembranaires d’anions chlorures. Les propriétés intrinsèques des sels utilisés qui en font des transporteurs efficaces ont été élucidées. Ainsi, les paramètres qui semblent affecter le plus le transport anionique sont le changement du contre-anion du sel d’imidazolium de même que la propension de ce dernier à s’auto-assembler via une succession d’empilements-π. De plus, les propriétés du transport ont été élucidées, montrant la formation de canaux transmembranaires qui permettent non-seulement la diffusion d’ions Cl-, mais aussi le transport de protons et d’ions Ca2+. L’intérêt de cette recherche repose d’abord dans le traitement de diverses pathologies voyant leur origine dans le dysfonctionnement du transport anionique. Cependant, les propriétés bactéricides des sels d’imidazolium utilisés ont été identifiées lors des dernières expériences.

Analysis of Fluctuations in the Early Universe Using Fokker-Planck Formalism

The thesis begins with a review of basic elements of general theory of relativity (GTR) which forms the basis for the theoretical interpretation of the observations in cosmology. The first chapter also discusses the standard model in cosmology, namely the Friedmann model, its predictions and problems. We have also made a brief discussion on fractals and inflation of the early universe in the first chapter. In the second chapter we discuss the formulation of a new approach to cosmology namely a stochastic approach. In this model, the dynam ics of the early universe is described by a set of non-deterministic, Langevin type equations and we derive the solutions using the Fokker—Planck formalism. Here we demonstrate how the problems with the standard model, can be eliminated by introducing the idea of stochastic fluctuations in the early universe. Many recent observations indicate that the present universe may be approximated by a many component fluid and we assume that only the total energy density is conserved. This, in turn, leads to energy transfer between different components of the cosmic fluid and fluctuations in such energy transfer can certainly induce fluctuations in the mean to factor in the equation of state p = wp, resulting in a fluctuating expansion rate for the universe. The third chapter discusses the stochastic evolution of the cosmological parameters in the early universe, using the new approach. The penultimate chapter is about the refinements to be made in the present model, by means of a new deterministic model The concluding chapter presents a discussion on other problems with the conventional cosmology, like fractal correlation of galactic distribution. The author attempts an explanation for this problem using the stochastic approach.

Analysis of Stochastic Volatility Sequences Generated by Product Autoregressive Models

The classical methods of analysing time series by Box-Jenkins approach assume that the observed series uctuates around changing levels with constant variance. That is, the time series is assumed to be of homoscedastic nature. However, the nancial time series exhibits the presence of heteroscedasticity in the sense that, it possesses non-constant conditional variance given the past observations. So, the analysis of nancial time series, requires the modelling of such variances, which may depend on some time dependent factors or its own past values. This lead to introduction of several classes of models to study the behaviour of nancial time series. See Taylor (1986), Tsay (2005), Rachev et al. (2007). The class of models, used to describe the evolution of conditional variances is referred to as stochastic volatility modelsThe stochastic models available to analyse the conditional variances, are based on either normal or log-normal distributions. One of the objectives of the present study is to explore the possibility of employing some non-Gaussian distributions to model the volatility sequences and then study the behaviour of the resulting return series. This lead us to work on the related problem of statistical inference, which is the main contribution of the thesis

Time-delayed fronts from biased random walks

We generalize a previous model of time-delayed reaction–diffusion fronts (Fort and Méndez 1999 Phys. Rev. Lett. 82 867) to allow for a bias in the microscopic random walk of particles or individuals. We also present a second model which takes the time order of events (diffusion and reproduction) into account. As an example, we apply them to the human invasion front across the USA in the 19th century. The corrections relative to the previous model are substantial. Our results are relevant to physical and biological systems with anisotropic fronts, including particle diffusion in disordered lattices, population invasions, the spread of epidemics, etc

Anisotropic dispersion, space competition, and the slowdown of the Neolithic transition

The front speed of the Neolithic (farmer) spread in Europe decreased as it reached Northern latitudes, where the Mesolithic (huntergatherer) population density was higher. Here, we describe a reaction diffusion model with (i) an anisotropic dispersion kernel depending on the Mesolithic population density gradient and (ii) a modified population growth equation. Both effects are related to the space available for the Neolithic population. The model is able to explain the slowdown of the Neolithic front as observed from archaeological data