Wave optics of Liesegang rings

Periodic Precipitation in a Confined Liquid Layer

Pattern formation is a ubiquitous phenomenon in animate and inanimate systems generated by mass transport and reaction of chemical species. The Liesegang phenomenon is a self-organized periodic precipitation pattern always studied in porous media such as hydrogels and aerogels for over a century. The primary consideration of applying the porous media is to prevent the disintegration of the precipitation structures due to the sedimentation of the precipitate and induced fluid flow. Here, we show that the periodic precipitation patterns can be engineered using a Hele-Shaw cell in a confined liquid phase, restricting hydrodynamic instability. The patterns generated in several precipitation reaction systems exhibit spatiotemporal properties consistent with patterns obtained in solid hydrogels. Furthermore, analysis considering the Rayleigh-Darcy number emphasizes the crucial role of fluidity in generating periodic precipitation structures in a thin liquid film. This exploration promises breakthroughs at the intersection of fundamental understanding and practical applications.

Diffusion-controlled annihilation A+B→0: Coalescence, fragmentation, and collapse of nonidentical A-particle islands submerged in the B-particle sea

We present a systematic analysis of diffusion-controlled interaction and collapse of two nonidentical spatially separated d-dimensional A-particle islands in the B-particle sea at propagation of the sharp reaction front A+B→0 at equal species diffusivities. We show that at a sufficiently large initial distance between the centers of islands 2ℓ and a relatively large initial ratio of island-to-sea concentrations, the evolution dynamics of the island-sea-island system demonstrates remarkable universality and, depending on the system dimension, is determined unambiguously by two dimensionless parameters Λ=N_{0}^{+}/N_{Ω} and q=N_{0}^{-}/N_{0}^{+}, where N_{0}^{+} and N_{0}^{-} are the initial particle numbers in the larger and smaller of the islands, respectively, and N_{Ω} is the initial number of sea particles in the volume Ω=(2ℓ)^{d}. We find that at each fixed 0<q≤1, there are threshold values Λ_{★}(q) and Λ_{s}(q)≥Λ_{★}(q) that depend on the dimension and separate the domains of individual death of each of the islands Λ<Λ_{★}(q), coalescence and subsequent fragmentation (division) of a two-centered island Λ_{★}(q)<Λ<Λ_{s}(q), and collapse of a single-centered island formed by coalescence Λ>Λ_{s}(q). We demonstrate that regardless of d, the trajectories of the island centers are determined unambiguously by the parameter q, and we reveal a detailed picture of the evolution of islands and front trajectories with an increase in Λ, focusing on the scaling laws of evolution at the final collapse stage and in the vicinity of starting coalescence and fragmentation points.

Chemical Tracking of Temperature by Concurrent Periodic Precipitation Pattern Formation in Polyacrylamide Gels

In nature, nonequilibrium systems reflect environmental changes, and these changes are often "recorded" in their solid body as they develop. Periodic precipitation patterns, aka Liesegang patterns (LPs), are visual sums of complex events in nonequilibrium reaction-diffusion processes. Here we aim to achieve an artificial system that "records" the temperature changes in the environment with the concurrent LP formation. We first illustrate the differences in 1-D LPs developing at different temperatures in terms of band spacings, which can demonstrate the time, ramp steepness, and extent of a temperature change. These results are discussed and augmented by a mathematical model. Using scanning electron microscopy, we show that the average size of the CuCrO₄ precipitate also reflects the temperature changes. Finally, we show that these changes can also be "recorded" in the 2-D and 3-D LPs, which can have applications in long-term temperature tracking and complex soft material design.

Modification of the Matalon-Packter Law for Self-Organized Periodic Precipitation Patterns by Incorporating Time-Dependent Diffusion Flux

Spontaneous pattern formation is common in both inanimate and living systems. Although the Liesegang pattern (LP) is a well-studied chemical model for precipitation patterns, various recent LP systems based on artificial control could not be easily evaluated using classical tools. The Matalon-Packter (MP) law describes the effect of the initial electrolyte concentration, which governs the diffusion flux (Fdiff), on the spatial distribution of LP. Note that the classical MP law only considers Fdiff through the initial concentration of electrolytes, even though it should also depend on the volume of the reservoir used for the outer electrolyte because of the temporal change in the concentration therein due to diffusion. However, there has been no report on the relationship between the MP law, the reservoir volume, and Fdiff. Here, we experimentally demonstrated and evaluated the effect of the reservoir volume on LP periodicity according to the classical MP law. Numerical simulations revealed that the reservoir volume affects the temporal modulation of Fdiff. By expressing the MP law as a function of estimated Fdiff after a certain period of time, we provide a uniform description of the changes in periodicity for both small and large reservoir volumes. Such modification should make the MP law a more robust tool for studying LP systems.

Fine tuning of pattern selection in the cadmium-hydroxide-system

Controlling self-organization in precipitation reactions has received growing attention in the efforts of engineering highly ordered spatial structures. Experiments have been successful in regulating the band patterns of the Liesegang phenomenon on various scales. Herein, we show that by adjusting the composition of the hydrogel medium, we can switch the final pattern between the classical band structure and the rare precipitate spots with hexagonal symmetry. The accompanying modeling study reveals that besides the modification of gel property, tuning of the time scale of diffusional spreading of hydroxide ions with respect to that of the phase separation drives the mode selection between one-dimensional band and two-dimensional spot patterns.

Diffusion-controlled coalescence, fragmentation, and collapse of d-dimensional A-particle islands in the B-particle sea

We present a systematic analysis of diffusion-controlled interaction and collapse of two identical spatially separated d-dimensional A-particle islands in the B-particle sea at propagation of the sharp reaction front A+B→0 at equal species diffusivities. We show that at a sufficiently large initial distance between the centers of islands 2ℓ compared to their characteristic initial size and a relatively large initial ratio of island to sea concentrations, the evolution dynamics of the island-sea-island system is determined unambiguously by the dimensionless parameter Λ=N_{0}/N_{Ω}, where N_{0} is the initial particle number in the island and N_{Ω} is the initial number of sea particles in the volume Ω=(2ℓ)^{d}. It is established that (a) there is a d-dependent critical value Λ_{★} above which island coalescence occurs; (b) regardless of d the centers of each of the islands move toward each other along a universal trajectory merging in a united center at the d-dependent critical value Λ_{s}≥Λ_{★}; (c) in one-dimensional systems Λ_{★}=Λ_{s}, therefore, at Λ<Λ_{★} each of the islands dies individually, whereas at Λ>Λ_{★} coalescence is completed by collapse of a single-centered island in the system center; (d) in two- and three-dimensional systems in the range Λ_{★}<Λ<Λ_{s} coalescence is accompanied by subsequent fragmentation of a two-centered island and is completed by individual collapse of each of the islands. We discuss a detailed picture of coalescence, fragmentation, and collapse of islands focusing on evolution of their shape and on behavior of the relative width of the reaction front at the final collapse stage and in the vicinity of starting coalescence and fragmentation points. We demonstrate that in a wide range of parameters, the front remains sharp up to a narrow vicinity of the coalescence, fragmentation, and collapse points.

Regularization of the Ostwald supersaturation model for Liesegang bands

In a previous paper, we analysed the Keller-Rubinow formulation of Ostwald's supersaturation theory for the formation of Liesegang rings or Liesegang bands, and found that the model is ill-posed, in the sense that after the termination of the first crystal front growth, secondary bands form, as in the experiment, but these are numerically found to be a single grid space wide, and thus an artefact of the numerical method. This ill-posedness is due to the discontinuity in the crystal growth rate, which itself reflects the supersaturation threshold inherent in the theory. Here we show that the ill-posedness can be resolved by the inclusion of a relaxation mechanism describing an impurity coverage fraction, which physically enables the transition in heterogeneous nucleation from precipitate-free impurity to precipitate-covered impurity.

Role of Electrolyte in Liesegang Pattern Formation

Pattern formation based on the Liesegang phenomenon is considered one of the useful models for gaining a mechanistic understanding of spontaneous spatiotemporal pattern formations in nature. However, for more than a century, the Liesegang phenomenon in chemical systems has been investigated by using electrolytes as both the reaction substrate and aggregation promoter, which has obfuscated the role of the electrolyte. Here, we distinguish the electrolyte (Na₂SO₄) from the reaction substrates (Ag⁺ ion and citrate), where Na₂SO₄ does not participate in the reaction step and acts as an aggregation promoter. The addition of Na₂SO₄ in Ag⁺-citrate-type Liesegang rings gave well-resolvable clear bands with a larger spacing coefficient. The observed changes were discussed by using the classical DLVO (Derjaguin-Landau-Verwey-Overbeek) theory, where the role of the electrolyte is to shield the electrostatic repulsive interaction among the reaction products. Furthermore, the numerical simulation of the reaction-diffusion equation with different aggregation thresholds reproduced the salt-dependent change in the spacing coefficient. We expect that an understanding of the exact role of the electrolyte as the aggregation promoter reported here will offer novel insight into how nature spontaneously forms beautiful spatiotemporal patterns.

On the Keller-Rubinow model for Liesegang ring formation

We study the model of Keller & Rubinow (Keller & Rubinow 1981 J. Chem. Phys74, 5000-5007. (doi:10.1063/1.441752)) describing the formation of Liesegang rings due to Ostwald's supersaturation mechanism. Keller and Rubinow provided an approximate solution both for the growth and equilibration of the first band, and also for the formation of secondary bands, based on a presumed asymptotic limit. However, they did not provide a parametric basis for the assumptions in their solution, nor did they provide any numerical corroboration, particularly of the secondary band formation. Here, we provide a different asymptotic solution, based on a specific parametric limit, and we show that the growth and subsequent cessation of the first band can be explained. We also show that the model is unable to explain the formation of finite width secondary bands, and we confirm this result by numerical computation. We conclude that the model is not fully posed, lacking a transition variable which can describe the hysteretic switch across the nucleation threshold.

Diffusion-controlled formation and collapse of a d-dimensional A-particle island in the B-particle sea

We consider diffusion-controlled evolution of a d-dimensional A-particle island in the B-particle sea at propagation of the sharp reaction front A+B→0 at equal species diffusivities. The A-particle island is formed by a localized (point) A-source with a strength λ that acts for a finite time T. We reveal the conditions under which the island collapse time t_{c} becomes much longer than the injection period T (long-living island) and demonstrate that regardless of d the evolution of the long-living island radius r_{f}(t) is described by the universal law ζ_{f}=r_{f}/r_{f}^{M}=sqrt[eτ|lnτ|], where τ=t/t_{c} and r_{f}^{M} is the maximal island expansion radius at the front turning point t_{M}=t_{c}/e. We find that in the long-living island regime the ratio t_{c}/T changes with the increase of the injection period T by the law ∝(λ^{2}T^{2-d})^{1/d}, i.e., increases with the increase of T in the one-dimensional (1D) case, does not change with the increase of T in the 2D case and decreases with the increase of T in the 3D case. We derive the scaling laws for particles death in the long-living island and determine the limits of their applicability. We demonstrate also that these laws describe asymptotically the evolution of the d-dimensional spherical island with a uniform initial particle distribution generalizing the results obtained earlier for the quasi-one-dimensional geometry. As striking results, we present a systematic analysis of the front relative width evolution for fluctuation, logarithmically modified, and mean-field regimes, and we demonstrate that in a wide range of parameters the front remains sharp up to a narrow vicinity of the collapse point.

Self-organization in precipitation reactions far from the equilibrium

Far from the thermodynamic equilibrium, many precipitation reactions create complex product structures with fascinating features caused by their unusual origins. Unlike the dissipative patterns in other self-organizing reactions, these features can be permanent, suggesting potential applications in materials science and engineering. We review four distinct classes of precipitation reactions, describe similarities and differences, and discuss related challenges for theoretical studies. These classes are hollow micro- and macrotubes in chemical gardens, polycrystalline silica carbonate aggregates (biomorphs), Liesegang bands, and propagating precipitation-dissolution fronts. In many cases, these systems show intricate structural hierarchies that span from the nanometer scale into the macroscopic world. We summarize recent experimental progress that often involves growth under tightly regulated conditions by means of wet stamping, holographic heating, and controlled electric, magnetic, or pH perturbations. In this research field, progress requires mechanistic insights that cannot be derived from experiments alone. We discuss how mesoscopic aspects of the product structures can be modeled by reaction-transport equations and suggest important targets for future studies that should also include materials features at the nanoscale.

Mesoscopic study of the effects of gel concentration and materials on the formation of precipitation patterns

A mesoscopic model based on the lattice Boltzmann method (LBM) is proposed to simulate the formation of Liesegang precipitation patterns and investigate the effects of gel on the morphology of the precipitates. In this model, nucleation is introduced on the basis of Ostwald's supersaturation theory, and subsequent crystal growth on the precipitate surface is simulated using a crystal growth model by taking into account the heterogeneous reaction on the surface of the precipitate. This model can capture the porous structures of the precipitates and can take into account the effects of the gel concentration and material by adjusting the gel porosity and nucleation threshold. The density of the precipitate nodes in the model is limited. A wide range of precipitation patterns is predicted under different gel concentrations, including regular bands, treelike patterns, and for the first time with numerical models, transition patterns between regular bands and treelike patterns. Both the spacing law and the width law are carefully investigated with different gel concentrations and gel materials, and the obtained α in the width law (w(n) ≈ x(n)(α), where x(n) and w(n) are the position and the width of the nth band, respectively) is in the range of 0.58-0.67, sandwiched by previously predicted ranges.

Band, target, and onion patterns in Co(OH)2 Liesegang systems

The study of morphology and shape development has gained considerable interest in certain sciences, notably biology and geology. Liesegang experiments producing Co(OH)2 stratification are performed here, in one, two, and three dimensions for comparison of the pattern morphologies. We obtain well-resolved bands in one dimension, target patterns (rings) in two dimensions, and onion patterns (spherical shells) in three dimensions. The morphological characteristics of the various patterns (spacing coefficients, rate of growth of ring spacing with distance) were measured. The spacing ratio of the strata in the different spatial dimensions was found to be anticorrelated with the surface-to-volume ratio of the gel domain. Some studies featuring the importance of morphology in Liesegang systems are briefly surveyed.

Reaction-diffusion systems in intracellular molecular transport and control

Chemical reactions make cells work only if the participating chemicals are delivered to desired locations in a timely and precise fashion. Most research to date has focused on active-transport mechanisms, although passive diffusion is often equally rapid and energetically less costly. Capitalizing on these advantages, cells have developed sophisticated reaction-diffusion (RD) systems that control a wide range of cellular functions-from chemotaxis and cell division, through signaling cascades and oscillations, to cell motility. These apparently diverse systems share many common features and are "wired" according to "generic" motifs such as nonlinear kinetics, autocatalysis, and feedback loops. Understanding the operation of these complex (bio)chemical systems requires the analysis of pertinent transport-kinetic equations or, at least on a qualitative level, of the characteristic times of the constituent subprocesses. Therefore, in reviewing the manifestations of cellular RD, we also describe basic theory of reaction-diffusion phenomena.

Self-assembled right handed helical ribbons of the bone mineral hydroxyapatite

A self-assembled right handed helical ribbon of bone mineral, hydroxyapatite (HAp) was crystallized in sodium meta silicate gel matrix at 27 degrees C and the physiological pH (7.4). At temperatures 37 and 47 degrees C, helical structures were followed by many Liesegang rings. The samples were characterized by FT-IR, XRD, SEM, ICP-OES and TG-DTA techniques. The helical ribbon consisted of platy Ca-deficient apatite crystals of size 2.8 microm. Liesegang ring had a continuous network of fibers with interconnected pores. The samples exhibited bioactivity when soaked in SBF.

Guiding fields for phase separation: controlling Liesegang patterns

Liesegang patterns emerge from precipitation processes and may be used to build bulk structures at submicrometer length scales. Thus they have significant potential for technological applications provided adequate methods of control can be devised. Here we describe a simple, physically realizable pattern control based on the notion of driven precipitation, meaning that the phase separation is governed by a guiding field such as, for example, a temperature or pH field. The phase separation is modeled through a nonautonomous Cahn-Hilliard equation whose spinodal is determined by the evolving guiding field. Control over the dynamics of the spinodal gives control over the velocity of the instability front that separates the stable and unstable regions of the system. Since the wavelength of the pattern is largely determined by this velocity, the distance between successive precipitation bands becomes controllable. We demonstrate the above ideas by numerical studies of a one-dimensional system with a diffusive guiding field. We find that the results can be accurately described by employing a linear stability analysis (pulled-front theory) for determining the velocity-local-wavelength relationship. From the perspective of the Liesegang theory, our results indicate that the so-called revert patterns may be naturally generated by diffusive guiding fields.

Snell's law of refraction observed in thermal frontal polymerization

We demonstrate that Snell's law of refraction can be applied to thermal fronts propagating through a boundary between regions that support distinct frontal velocities. We use the free-radical frontal polymerization of a triacrylate with clay filler that allows for two domains containing two different concentrations of a peroxide initiator to be molded together. Because the polymerization reaction rates depend on the initiator concentration, the propagation speed is different in each domain. We study fronts propagating in two parallel strips in which the incident angle is 90 degrees. Our data fit Snell's law v(r)/v(i)=sin theta(r)/sin theta(i), where v(r) is the refracted velocity, v(i) is the incident velocity, theta(r) is the angle of refraction, and theta(i) is the incident angle. Further, we study circular fronts propagating radially from an initiation point in a high-velocity region into a low-velocity region (and vice versa). We demonstrate the close resemblance between the numerically simulated and experimentally observed thermal reaction fronts. By measuring the normal velocity and the angle of refraction of both simulated and experimental fronts, we establish that Snell's law holds for thermal frontal polymerization in our experimental system. Finally we discuss the regimes in which Snell's law may not be valid.

Pattern formation induced by ion-selective surfaces: models and simulations

Simple inorganic reactions in gels, such as NaOH + CuCl(2), NaOH + AgNO(3), and CuCl(2) + K(3)[Fe(CN)(6)], can yield to various precipitation patterns. The first compound penetrates in a hydrogel by diffusion, and reacts with the second compound homogenized in the gel. The precipitate patterns formed in these reactions have got two kinds of bordering surfaces. Recent experimental results suggested that one of these surfaces has an ion-selective (semipermeable) character: It restrains the diffusion of the reacting ion contained by the reactant that diffuses into the gel. In this paper, we built the above experimental observation into a reaction-diffusion cellular-automata model of the pattern formation. Computer simulations showed that the model is able to reproduce the basic building elements of the patterns.

Nano- and Microscopic Surface Wrinkles of Linearly Increasing Heights Prepared by Periodic Precipitation

Arrays of surface wrinkles of linearly increasing heights (from tens of nanometers to tens of micrometers) were prepared via a spontaneous reaction-diffusion process based on periodic precipitation. The slopes, dimensions, and positions of the precipitation bands could be controlled precisely by adjusting the concentrations of the participating chemicals as well as the material properties of patterned substrates. Additional control of periodic precipitation by localized UV irradiation allowed for the preparation of discontinuous and curvilinear structures. The nonbinary 3D surface topographies were replicated into poly(dimethylsiloxane), and the applications of replicas in microfluidics, microseparations, and cell biology have been suggested.

Formation of Liesegang patterns in the presence of an electric field

The effects of an external electric field on the formation of Liesegang patterns are investigated. The patterns are assumed to emerge from a phase separation process in the wake of a diffusive reaction front. The dynamics is described by a Cahn-Hilliard equation with a moving source term representing the reaction zone, and the electric field enters through its effects on the properties of the reaction zone. We employ our previous results [I. Bena, F. Coppex, M. Droz, and Z. Rácz, J. Chem. Phys. 122, 024512 (2005)] on how the electric field changes both the motion of the front, as well as the amount of reaction product left behind the front, and our main conclusion is that the number of precipitation bands becomes finite in a finite electric field. The reason for the finiteness in case when the electric field drives the reagents towards the reaction zone is that the width of consecutive bands increases so that, beyond a distance l(+), the precipitation is continuous (plug is formed). In case of an electric field of opposite polarity, the bands emerge in a finite interval l(-), since the reaction product decreases with time and the conditions for phase separation cease to exist. We give estimates of l(+/-) in terms of measurable quantities and thus present an experimentally verifiable prediction of the "Cahn-Hilliard equation with a moving source" description of Liesegang phenomena.

Regular Liesegang patterns and precipitation waves in an open system

We investigate the regular and moving Liesegang pattern formation phenomena in an open system. First, simulations have been performed at fixed coupling between the reactive medium and the reservoir, later this control parameter was varied during the simulations resulting in various phenomena. We predicted and monitored for the first time various--dynamically changing--precipitation structures and a spatial hysteresis phenomenon, which is beyond the scope of the Turing instability. The dynamics of the reaction is well detectable using specific quantities: the total amount of precipitate and its center of gravity.