Mesoscopic Modeling in the Kinetic Theory of Adsorbates

Structure formation in a conserved mass model of a set of individuals interacting with attractive and repulsive forces

We study a set of interacting individuals that conserve their total mass. In order to describe its dynamics we resort to mesoscopic equations of reaction diffusion including currents driven by attractive and repulsive forces. For the mass conservation we consider a linear response parameter that maintains the mass in the vicinity of a optimal value which is determined by the set. We use the reach and intensity of repulsive forces as control parameters. When sweeping a wide range of parameter space we find a great diversity of localized structures, stationary as well as other ones with cyclical and chaotic dynamics. We compare our results with real situations.

The Generalized Stochastic Smoluchowski Equation

We study the dynamics of a system of overdamped Brownian particles governed by the generalized stochastic Smoluchowski equation associated with a generalized form of entropy and involving a long-range potential of interaction [P.H. Chavanis, Entropy 17, 3205 (2015)]. We first neglect fluctuations and provide a macroscopic description of the system based on the deterministic mean field Smoluchowski equation. We then take fluctuations into account and provide a mesoscopic description of the system based on the stochastic mean field Smoluchowski equation. We establish the main properties of this equation and derive the Kramers escape rate formula, giving the lifetime of a metastable state, from the theory of instantons. We relate the properties of the generalized stochastic Smoluchowski equation to a principle of maximum dissipation of free energy. We also discuss the connection with the dynamical density functional theory of simple liquids.

Nano-sized Adsorbate Structure Formation in Anisotropic Multilayer System

In this article, we study dynamics of adsorbate island formation in a model plasma-condensate system numerically. We derive the generalized reaction-diffusion model for adsorptive multilayer system by taking into account anisotropy in transfer of adatoms between neighbor layers induced by electric field. It will be found that with an increase in the electric field strength, a structural transformation from nano-holes inside adsorbate matrix toward separated nano-sized adsorbate islands on a substrate is realized. Dynamics of adsorbate island sizes and corresponding distributions are analyzed in detail. This study provides an insight into details of self-organization of adatoms into nano-sized adsorbate islands in anisotropic multilayer plasma-condensate systems.

1/f 2 noise in bistable electrocatalytic reactions on mesoscale electrodes

The formation of a self-organized spatial domain during current-controlled CO oxidation, a kinetically bistable reaction, is investigated experimentally and by deterministic simulations as a function of the electrode size and of the supporting electrolyte concentration. Decreasing the microelectrode size leads to the suppression of the spatial instability at the electrode and thus stabilizes the S-NDR branch of the reaction. The critical microelectrode size capable of supporting sustained domain formation is shown to be strongly affected by the sulfuric acid concentration, the characteristic time of the positive feedback loop increasing with the sulfate concentration. Furthermore, we demonstrate that for microelectrode diameters close to the instability threshold, small amplitude electrochemical potential fluctuations appear in the system. These potential fluctuations cannot be captured by deterministic mathematical models and are attributed to a strong enhancement of molecular fluctuations or intrinsic noise in the vicinity of the spatial instability. Analysis of the electrochemical noise revealed a 1/f ² frequency dependence and several common features with neuronal shot noise.

Calculation of displacement correlation tensor indicating vortical cooperative motion in two-dimensional colloidal liquids

As an indicator of cooperative motion in a system of Brownian particles that models two-dimensional colloidal liquids, a displacement correlation tensor is calculated analytically and compared with numerical results. The key idea for the analytical calculation is to relate the displacement correlation tensor, which is a kind of four-point space-time correlation, to the Lagrangian two-time correlation of the deformation gradient tensor. Tensorial treatment of the statistical quantities, including the displacement correlation itself, allows capturing the vortical structure of the cooperative motion. The calculated displacement correlation also implies a negative long-time tail in the velocity autocorrelation, which is a manifestation of the cage effect. Both the longitudinal and transverse components of the displacement correlation are found to be expressible in terms of a similarity variable, suggesting that the cages are nested to form a self-similar structure in the space-time.

Equilibrium microphase separation in the two-leaflet model of lipid membranes

Because of the coupling between local lipid composition and the thickness of the membrane, microphase separation in two-component lipid membranes can take place; such effects may underlie the formation of equilibrium nanoscale rafts. Using a kinetic description, this phenomenon is analytically and numerically investigated. The phase diagram is constructed through the stability analysis for linearized kinetic equations, and conditions for microphase separation are discussed. Simulations of the full kinetic model reveal the development of equilibrium membrane nanostructures with various morphologies from the initial uniform state.

Effects of molecular noise on bistable protein distributions in rod-shaped bacteria

The distributions of many proteins in rod-shaped bacteria are far from homogeneous. Often they accumulate at the cell poles or in the cell centre. At the same time, the copy number of proteins in a single cell is relatively small making the patterns noisy. To explore limits to protein patterns due to molecular noise, we studied a generic mechanism for spontaneous polar protein assemblies in rod-shaped bacteria, which are based on cooperative binding of proteins to the cytoplasmic membrane. For mono-polar assemblies, we find that the switching time between the two poles increases exponentially with the cell length and with the protein number. This feature could be beneficial to organelle maintenance in ageing bacteria.

Thermodynamics of the one-dimensional parallel Kawasaki model: exact solution and mean-field approximations

The adsorption isotherm for the recently proposed parallel Kawasaki (PK) lattice-gas model [Phys. Rev. E 88, 062144 (2013)] is calculated exactly in one dimension. To do so, a third-order difference equation for the grand-canonical partition function is derived and solved analytically. In the present version of the PK model, the attraction and repulsion effects between two neighboring particles and between a particle and a neighboring empty site are ruled, respectively, by the dimensionless parameters ϕ and θ. We discuss the inflections induced in the isotherms by situations of high repulsion, the role played by finite lattice sizes in the emergence of substeps, and the adequacy of the two most widely used mean-field approximations in lattice gases, namely, the Bragg-Williams and the Bethe-Peierls approximations.

Conciliating synchronicity with spatial discretization, exclusion, interactions, and detailed balance

The construction of a discrete stochastic system of interacting particles that evolves through a fully synchronous evolution rule while satisfying detailed balance is a highly demanding task. As a consequence, the presence of nontrivial interaction fields can make synchronicity and thermodynamic equilibrium look as two conflicting counterparts. We show that, with the proper prescriptions, the process of migration of particles in a lattice of mutually exclusive nodes can be simulated with a fully synchronous algorithm, which we call parallel Kawasaki dynamics (PKD), that incorporates site exclusion, local interactions, and detailed balance without the need of system partitioning schemes. We show that the underlying pseudo-Hamiltonian (which is derived from the PKD dynamics instead of being assumed a priori as usual in a sequential Monte Carlo scheme) is temperature dependent and causes the resulting equilibrium properties to differ substantially from the conventional hopping model when the system is near critical conditions.

Synchronous equilibrium model for the diffusion of mutually exclusive particles in a heterogeneous lattice of adsorption sites

Through straight synchronization and proper manipulation of a sequential Monte Carlo glass-forming rule introduced by Fröbose and Jäckle [J. Stat. Phys. 42, 551 (1986)], we constructed a synchronous, non-glass-forming rule for diffusion of mutually exclusive particles in a lattice of adsorption sites. The rule satisfies detailed balance in the presence of both homogeneous and heterogeneous adsorption energies. Our model differs from the usual lattice-gas cellular automata diffusion rules in that the mutual exclusion holds on the lattice sites rather than on the channels which connect neighboring sites, and from the mass-conserving cellular automata rules in the use of a no-partitioning scheme. The first aim of this work is to show that, although some prescriptions in the synchronous rule are introduced just to allow that both detailed balance and mutual exclusion can coexist with synchronicity, the diffusion process produced by the rule is not anomalous so that the rule can be regarded as a diffusion model. We then compare the diffusion isotherms of several test systems with the ones obtained by means of sequential Monte Carlo simulations of Arrhenius jumps of particles on a lattice. Finally, we apply the rule to the case of a (100) fcc model surface and estimate the amount of time correlation in the migration process, and show that the synchronous rule produces higher correlations and slightly lower diffusivity than the sequential Monte Carlo counterpart.

Nanosize pattern formation in overdamped stochastic reaction-diffusion systems with interacting adsorbate

We study overdamped stochastic model describing adsorption or desorption processes with nonequilibrium chemical reactions on the surface. It is shown that internal noise satisfying the fluctuation-dissipation relation at small intensities governs transitions between ordered thermodynamical dense and diluted phases. These phase transitions are characterized by an increase of fluctuations of the coverage filed and correlation radius of spatial modulation. At large noise intensity a transition towards disordered phase with chaotic spatial configuration is realized. We have shown that organized stationary patterns are of nanometer range. We define that both period of stationary structures and corresponding correlation radius depend on the noise intensity.

Nonequilibrium pattern formation by subdominant attractive forces: a simple model and a stabilization strategy by means of noise

We introduce a simple model describing a mechanism for transient pattern formation driven by subdominant attractive forces. The patterns can be stabilized if they are confined by means of a particular multiplicative noise into the region where such mechanism is active. The scope of the results appears to transcend the original application context.

Metastability in Schloegl's second model for autocatalysis: Lattice-gas realization with particle diffusion

We analyze metastability associated with a discontinuous nonequilibrium phase transition in a stochastic lattice-gas realization of Schloegl's second model for autocatalysis. This model realization involves spontaneous annihilation, autocatalytic creation, and diffusion of particles on a square lattice, where creation at empty sites requires an adjacent diagonal pair of particles. This model, also known as the quadratic contact process, exhibits discontinuous transition between a populated active state and a particle-free vacuum or "poisoned" state, as well as generic two-phase coexistence. The poisoned state exists for all particle annihilation rates p>0 and hop rates h≥0 and is an absorbing state in the sense of Markovian processes. The active or reactive steady state exists only for p below a critical value, p{e}=p{e}(h) , but a metastable extension appears for a range of higher p up to an effective upper spinodal point, p{s+}=p{s+}(h) (i.e., p{s+}>p{e} ). For selected h , we assess the location of p{s+}(h) by characterizing both the poisoning kinetics and the propagation of interfaces separating vacuum and active states as a function of p .

Reversible formation of nanodomains in monolayers of DPPC studied by kinetic modeling

Dipalmitoylphosphatidylcholine (DPPC) is the most abundant component in pulmonary surfactants and is believed to be responsible for maintaining low surface tension in alveoli during breathing. In this work, a kinetic model is introduced that describes the phase separation in DPPC films that produces the liquid-condensed (LC) and liquid-expanded (LE) fractions, which differ according to the area density of DPPC. The phase separation in an initially homogeneous film has been investigated numerically. Furthermore, explicit simulations of periodic compression-expansion cycles are reported. In this process, a moderate change of the surface area resulted in a dramatic change in the total amount of LC fraction, as well as in the surface morphology. Depending on the extent of the film's compression, the simulated surface morphologies comprised individual nanosized LC domains embedded in the LE fraction, interconnected networks of such domains, or continuous LC films with nanopores. Equilibration of the total area of the LC nanodomains occurred over a few milliseconds, indicating that the rate of the LE-LC phase transformation is sufficient for maintaining low surface tension during breathing, and that nanoscale LC domains are likely to play a major role in this process. Unlike the total content of the LC fraction, which stabilized quickly, the average size of LC nanodomains showed a tendency to increase slowly, at a rate determined by the diffusivity of DPPC. The computed average domain size seems to be compatible with published experiments for DPPC films. The numeric results also elucidate the distinction between thermodynamically determined and kinetically limited structural properties during phase separation in the major structure-forming component of pulmonary surfactants.

Nonequilibrium microstructures in reactive monolayers as soft matter systems

Chemical systems provide classical examples of nonequilibrium pattern formation. Reactions in weak aqueous solutions, such as the extensively investigated Belousov-Zhabotinsky reaction, demonstrate a rich variety of patterns, ranging from travelling fronts to rotating spiral waves and chemical turbulence. Pattern formation in such systems is based on interplay between the reactions and diffusion. Intrinsically, this puts a restriction on the minimum length scale of the developing structures, which cannot be shorter than the diffusion length of the reactants. However, much smaller nonequilibrium structures, with characteristic lengths reaching down to nanoscales, are also possible. They are found in reactive soft matter, where energetic interactions between molecules are present as well. In these systems, chemical reactions and diffusion interfere with phase transitions, yielding active, stationary or dynamic microstructures. Nonequilibrium soft-matter microstructures are of fundamental importance for biological cells and may have interesting engineering applications. In this Minireview, we focus on the microstructures found in reactive soft-matter monolayers at solid surfaces or liquid-air interfaces.

Schloegl's second model for autocatalysis with particle diffusion: Lattice-gas realization exhibiting generic two-phase coexistence

We analyze a discontinuous nonequilibrium phase transition between an active (or reactive) state and a poisoned (or extinguished) state occurring in a stochastic lattice-gas realization of Schloegl's second model for autocatalysis. This realization, also known as the quadratic contact process, involves spontaneous annihilation, autocatalytic creation, and diffusion of particles on a square lattice, where creation at empty sites requires a suitable nearby pair of particles. The poisoned state exists for all annihilation rates p>0 and is an absorbing particle-free "vacuum" state. The populated active steady state exists only for p below a critical value, p(e). If p(f) denotes the critical value below which a finite population can survive, then we show that p(f)<p(e). This strict inequality contrasts a postulate of Durrett, and is a direct consequence of the occurrence of coexisting stable active and poisoned states for a finite range p(f)<or=p<or=p(e) (which shrinks with increasing diffusivity). This so-called generic two-phase coexistence markedly contrasts behavior in thermodynamic systems. However, one still finds metastability and nucleation phenomena similar to those in discontinuous equilibrium transitions.

Continuum mesoscopic framework for multiple interacting species and processes on multiple site types and/or crystallographic planes

While recently derived continuum mesoscopic equations successfully bridge the gap between microscopic and macroscopic physics, so far they have been derived only for simple lattice models. In this paper, general deterministic continuum mesoscopic equations are derived rigorously via nonequilibrium statistical mechanics to account for multiple interacting surface species and multiple processes on multiple site types and/or different crystallographic planes. Adsorption, desorption, reaction, and surface diffusion are modeled. It is demonstrated that contrary to conventional phenomenological continuum models, microscopic physics, such as the interaction potential, determines the final form of the mesoscopic equation. Models of single component diffusion and binary diffusion of interacting particles on single-type site lattice and of single component diffusion on complex microporous materials' lattices consisting of two types of sites are derived, as illustrations of the mesoscopic framework. Simplification of the diffusion mesoscopic model illustrates the relation to phenomenological models, such as the Fickian and Maxwell-Stefan transport models. It is demonstrated that the mesoscopic equations are in good agreement with lattice kinetic Monte Carlo simulations for several prototype examples studied.

Fluctuation-induced phase transition in a spatially extended model for catalytic CO oxidation

A reaction-diffusion master equation has been introduced in order to model the bistable CO oxidation on single crystal metal surfaces at high pressure where the diffusion length becomes small and local fluctuations are important. Analytical solutions can be found in a reduced one-component nonlinear master equation after applying the Weiss mean-field approximation together with the adiabatic elimination of oxygen. It is shown that the Weiss mean-field approximation predicts a symmetry-breaking bifurcation associated with a phase transition. The corresponding stationary solutions of the nonlinear master equation are supported by Gillespie-type Monte Carlo simulations.

Min-oscillations in Escherichia coli induced by interactions of membrane-bound proteins

During division it is of primary importance for a cell to correctly determine the site of cleavage. The bacterium Escherichia coli divides in the center, producing two daughter cells of equal size. Selection of the center as the correct division site is in part achieved by the Min-proteins. They oscillate between the two cell poles and thereby prevent division at these locations. Here, a phenomenological description of these oscillations is presented, where lateral interactions between proteins on the cell membrane play a key role. Solutions to the dynamic equations are compared to experimental findings. In particular, the temporal period of the oscillations is measured as a function of the cell length and found to be compatible with the theoretical prediction.

Spectral method study of domain coarsening

The self-organization of particles in a two phase system in the coexistence region through a diffusive mechanism is known as Ostwald ripening. The late stage of Ostwald ripening is studied here through the use of a mesoscopic model in conjunction with a special configuration that allows for the direct measurement of system characteristics such as domain size. The mesoscopic model is a stochastic partial integrodifferential equation and is studied through the use of recently developed and benchmarked spectral schemes. The size of the growing region is not observed to increase as a power law in this model during the late stages of self-organization in contrast to the power law growth observed in simulations of the earlier stages of self-organization. The results included here also demonstrate the effect of adjusting the interparticle interaction on the morphological evolution of the system.

Surface waves of Min-proteins

In the bacterium Escherichia coli, the Min-proteins show pronounced pole-to-pole oscillations. They are functional for suppressing cell division at the cell ends, leaving the center as the only possible site for division. Analyzing different models of Min-protein dynamics in a bacterial geometry, we find waves on the cytoplasmic membrane. Interestingly, the surface wave solutions of different models belong to different symmetry classes. We suggest that experiments on Min-protein surface waves in vitro are helpful in distinguishing between different classes of models of Min-protein dynamics.

Theoretical analysis of internal fluctuations and bistability in CO oxidation on nanoscale surfaces

The bistable CO oxidation on a nanoscale surface is characterized by a limited number of reacting molecules on the catalytic area. Internal fluctuations due to finite-size effects are studied by the master equation with a Langmuir-Hinshelwood mechanism for CO oxidation. Analytical solutions can be found in a reduced one-component model after the adiabatic elimination of one variable which in our case is the oxygen coverage. It is shown that near the critical point, with decreasing surface area, one cannot distinguish between two macroscopically stable stationary states. This is a consequence of the large fluctuations in the coverage which occur on a fast time scale. Under these conditions, the transition times between the macroscopic states also are no longer separated from the short-time scale of the coverage fluctuations as is the case for large surface areas and far away from the critical point. The corresponding stationary solutions of the probability distribution and the mean first passage times calculated in the reduced model are supported by numerics of the full two-component model.

Molecular sieve valves driven by adsorbate-adsorbate interactions: hysteresis in permeation of microporous membranes

A recently derived mesoscopic framework describing activated micropore diffusion is employed to explore system criticality in microporous membranes under nonequilibrium conditions. Rapid exploration of parameter space, possible with this continuum framework, elucidates a novel temperature-induced ignition and extinction of the molecular flux under a macroscopic gradient in pressure (chemical potential). Deviation from equilibrium like phase behavior (i.e., shifting and narrowing of phase envelopes and double hysteresis) derives from asymmetry of the coupled boundaries of the nonequilibrium membrane. We confirm this new phase behavior, akin to "opening" and "closing" of a molecular valve, via gradient kinetic Monte Carlo simulations of thin one-dimensional and three-dimensional systems. The heat of adsorption, strength of adsorbate-adsorbate intermolecular forces, and chemical potential gradient are all shown to control 'valve' actuation, suggesting potential implications in chemical sensing and novel diffusion control.

Effects of internal noise for rate oscillations during CO oxidation on platinum surfaces

We have studied the influence of internal noise on the reaction rate oscillation during carbon-monoxide oxidation on single crystal platinum surfaces using chemical Langevin equations. Considering that the surface is divided into small well-mixed cells, we have focused on the dynamic behavior inside a single cell. Internal noise can induce rate oscillations and the performance of the stochastic rate oscillations shows double maxima with the variation of the internal noise intensity, demonstrating the occurrence of internal noise coherent biresonance. The relationship between such a phenomenon with the deterministic bifurcation features of the system is also discussed.

Model for spreading of liquid monolayers

Manipulating fluids at the nanoscale within networks of channels or chemical lanes is a crucial challenge in developing small scale devices to be used in microreactors or chemical sensors. In this context, ultrathin (i.e., monolayer) films, experimentally observed in spreading of nanodroplets or upon extraction from reservoirs in capillary rise geometries, represent an extreme limit which is of physical and technological relevance since the dynamics is governed solely by capillary forces. In this work we use kinetic Monte Carlo (KMC) simulations to analyze in detail a simple, but realistic model proposed by Phys. Rev. Lett. 76, 86 (1996)]] for the two-dimensional spreading on homogeneous substrates of a fluid monolayer which is extracted from a reservoir. Our simulations confirm the previously predicted time dependence of the spreading, X ( t--> infinity ) =A square root of t, with X (t) as the average position of the advancing edge at time t, and they reveal a nontrivial dependence of the prefactor A on the strength U0 of interparticle attraction and on the fluid density C0 at the reservoir as well as an U0 -dependent spatial structure of the density profile of the monolayer. The asymptotic density profile at long time and large spatial scale is carefully analyzed within the continuum limit. We show that including the effect of correlations in an effective manner into the standard mean-field description leads to predictions both for the value of the threshold interaction above which phase segregation occurs and for the density profiles in excellent agreement with KMC simulation results.

Promoter-induced reactive phase separation in surface reactions

Promoters are adsorbed mobile species which do not directly participate in a catalytic surface reaction, but can influence its rate. Often, they are characterized by strong attractive interactions with one of the reactants. We show that these conditions lead to a Turing instability of the uniform state and to the formation of reaction-induced periodic concentration patterns. Experimentally such patterns are observed in catalytic water formation on a Rh(110) surface in the presence of coadsorbed potassium.

Nonequilibrium orientational patterns in two-component Langmuir monolayers

A model of a phase-separating two-component Langmuir monolayer in the presence of a photoinduced reaction interconverting two components is formulated. An interplay between phase separation, orientational ordering, and reaction is found to lead to a variety of nonequilibrium self-organized patterns, both stationary and traveling. Examples of the patterns, observed in numerical simulations, include flowing droplets, traveling stripes, wave sources, and vortex defects.

Lattice density functional theory of molecular diffusion

A density functional theory of diffusion is developed for lattice fluids with molecular flux as a functional of the density distribution. The formalism coincides exactly with the generalized Ono-Kondo density functional theory when there is no gradient of chemical potential, i.e., at equilibrium. Away from equilibrium, it gives Fick's first law in the absence of a potential energy gradient, and it departs from Fickian behavior consistently with the Maxwell-Stefan formulation. The theory is applied to model a nanopore, predicting nonequilibrium phase transitions and the role of surface diffusion in the transport of capillary condensate.

Coarse-grained stochastic processes for microscopic lattice systems

Diverse scientific disciplines ranging from materials science to catalysis to biomolecular dynamics to climate modeling involve nonlinear interactions across a large range of physically significant length scales. Here a class of coarse-grained stochastic processes and corresponding Monte Carlo simulation methods, describing computationally feasible mesoscopic length scales, are derived directly from microscopic lattice systems. It is demonstrated below that the coarse-grained stochastic models can capture large-scale structures while retaining significant microscopic information. The requirement of detailed balance is used as a systematic design principle to guarantee correct noise fluctuations for the coarse-grained model. The coarse-grained stochastic algorithms provide large computational savings without increasing programming complexity or computer time per executive event compared to microscopic Monte Carlo simulations.

Self-organized nanostructures in surface chemical reactions: Mechanisms and mesoscopic modeling

Nanoscale patterns can form in reactive adsorbates on catalytic surfaces as a result of attractive lateral interactions. These structures can be described within a mesoscopic theory that is derived by coarse graining the microscopic master equation thus providing a link between microscopic lattice models and reaction-diffusion equations. Such mesoscopic models allow to systematically investigate mechanisms responsible for the formation of nanoscale nonequilibrium patterns in reactive condensed matter. We have found that stationary and traveling nanostructures may result from the interplay of the attractive lateral interactions and nonequilibrium reactions. Besides reviewing these results, a detailed investigation of a single reactive adsorbate in the presence of attractive lateral interactions and global coupling through the gas phase is presented. Finally, it is outlined how a mesoscopic theory should be constructed for a particular scanning tunneling microscopy experiment [the oxidation of hydrogen on a Pt(111) surface] in order to overcome the failure of a corresponding reaction-diffusion model to quantitatively reproduce the experiments. (c) 2002 American Institute of Physics.

Microscopic aspects of pattern formation on surfaces

Recent scanning tunneling microscopy (STM) work gives insight into microscopic processes of surface reactions that play a role for spatio-temporal pattern formation. STM allows to resolve adsorbed particles, follow their surface motion, and monitor reactions with other particles on the atomic scale. The data reveal pronounced deviations from the implicite assumptions of the reaction-diffusion equations traditionally used to model spatio-temporal patterns. In contrast to these descriptions, particles are often not randomly distributed, but cluster in islands because of attractive interactions, and particle hopping can be highly correlated. It is shown that such phenomena can even affect the macroscopic kinetics. The article also discusses a case where the atomic processes inside propagating reaction fronts could be resolved. Here particular strong interaction effects were observed, caused by hydrogen bonds between the reacting species. (c) 2002 American Institute of Physics.

Stochastic and mesoscopic models for tropical convection

A new way to parametrize certain aspects of tropical convection through stochastic and mesoscopic models is developed here. The technical idea is to adapt tools from statistical physics and materials science to model important unresolved features of tropical convection. This new strategy consists of modeling the unresolved effects of convective inhibition in a coarse mesh mesoscopic parametrization through a "heat bath" model involving a stochastic spin flip model with very natural interaction rules for convective inhibition combined with a suitable external potential defined by the coarse mesh values. In turn, the values of the order parameter from this heat bath alter the vertical mass flux in regions of deep convection. Both stochastic and systematic deterministic mesoscopic parametrizations are developed here. The deterministic mesoscopic models derived in this fashion exhibit new phenomena such as multiple radiative equilibria in suitable parameter regimes. The simplest first numerical experiments reported here with the mesoscopic deterministic parametrization qualitatively reproduce several key features of the observational record regarding convectively coupled tropical waves. The systematic stochastic modeling strategy proposed here could also be very useful for capturing other features of tropical convection such as those involving cloud radiation feedbacks.

Spatiotemporal self-organization in a surface reaction: from the atomic to the mesoscopic scale

Scanning tunneling microscopy data revealed the atomic processes in propagating reaction fronts that occur in the catalytic oxidation of hydrogen on Pt(111). The fronts were also characterized on mesoscopic length scales with respect to their velocity and width. Simulations on the basis of a reaction-diffusion model reproduce the experimental findings qualitatively well. The quantitative comparison reveals the limitations of this traditional approach to modeling spatiotemporal pattern formation in nonlinear dynamics.