Kinematic reduction of reaction-diffusion fronts with multiplicative noise. Derivation of stochastic sharp-interface equations

We study the dynamics of generic reaction-diffusion fronts, including pulses and chemical waves, in the presence of multiplicative noise. We discuss the connection between the reaction-diffusion Langevin-like field equations and the kinematic (eikonal) description in terms of a stochastic moving-boundary or sharp-interface approximation. We find that the effective noise is additive and we relate its strength to the noise parameters in the original field equations, to first order in noise strength, but including a partial resummation to all orders which captures the singular dependence on the microscopic cutoff associated with the spatial correlation of the noise. This dependence is essential for a quantitative and qualitative understanding of fluctuating fronts, affecting both scaling properties and nonuniversal quantities. Our results predict phenomena such as the shift of the transition point between the pushed and pulled regimes of front propagation, in terms of the noise parameters, and the corresponding transition to a non-Kardar-Parisi-Zhang universality class. We assess the quantitative validity of the results in several examples including equilibrium fluctuations and kinetic roughening. We also predict and observe a noise-induced pushed-pulled transition. The analytical predictions are successfully tested against rigorous results and show excellent agreement with numerical simulations of reaction-diffusion field equations with multiplicative noise.

Relevance of dynamic wetting in viscous fingering patterns

We demonstrate that wetting effects at moving contact lines have a strong impact in viscous fingering patterns. Experiments in a rotating Hele-Shaw (HS) cell, dry or prewetted, show consistent morphological differences. When the wetting fluid invades a dry region, contact angle dynamics yield a kinetic contribution to the interface pressure drop that scales with capillary number as Ca2¿3 but is significantly larger than the Park-Homsy kinetic correction. Numerical results are in very good agreement with experiments and show that standard HS equations work best for prewetted cells.

Time-dependent couplings and crossover length scales in nonequilibrium surface roughening

We show that time-dependent couplings may lead to nontrivial scaling properties of the surface fluctuations of the asymptotic regime in nonequilibrium kinetic roughening models. Three typical situations are studied. In the case of a crossover between two different rough regimes, the time-dependent coupling may result in anomalous scaling for scales above the crossover length. In a different setting, for a crossover from a rough to either a flat or damping regime, the time-dependent crossover length may conspire to produce a rough surface, although the most relevant term tends to flatten the surface. In addition, our analysis sheds light into an existing debate in the problem of spontaneous imbibition, where time-dependent couplings naturally arise in theoretical models and experiments.

Distorted-wave calculation of cross sections for inner-shell ionization by electron and positron impact

The relativistic distorted-wave Born approximation is used to calculate differential and total cross sections for inner shell ionization of neutral atoms by electron and positron impact. The target atom is described within the independent-electron approximation using the self-consistent Dirac-Fock-Slater potential. The distorting potential for the projectile is also set equal to the Dirac-Fock-Slater potential. For electrons, this guarantees orthogonality of all the orbitals involved and simplifies the calculation of exchange T-matrix elements. The interaction between the projectile and the target electrons is assumed to reduce to the instantaneous Coulomb interaction. The adopted numerical algorithm allows the calculation of differential and total cross sections for projectiles with kinetic energies ranging from the ionization threshold up to about ten times this value. Algorithm accuracy and stability are demonstrated by comparing differential cross sections calculated by our code with the distorting potential set to zero with equivalent results generated by a more robust code that uses the conventional plane-wave Born approximation. Sample calculation results are presented for ionization of K- and L-shells of various elements and compared with the available experimental data.

Magnetic hardening mechanisms in Nd-Fe-B nanocrystalline material

Microstructural and magnetic measurements of the evolution by heat treatment of initially amorphous Nd16Fe76B8 alloys prepared by melt spinning are presented. Evidence of magnetic hardening above a threshold temperature induced by magnetic isolation of the Nd2Fe14B grains is provided. A thermodynamic and kinetic explanation of local melting of the intergranular nanostructured Nd¿rich eutectic phase at temperatures below 900 K based on capillary effects is presented. A subsequent Ostwald ripening process moves Nd to wet intimately the hard magnetic grains, becoming, on cooling, a real paramagnetic isolating thin film (~2.5 nm). By using a simple analogy, it is shown that the switching magnetization field in a single¿domain crystal can be drastically affected through the exchange coupling to neighboring grains with different orientation of the easy axis. This effect should be important enough to reinforce the coercive field of polycrystalline hard magnetic materials and explains the observed enhancement from 0.9 to 1.9 T.

Elastic scattering of fast electrons and positrons by atoms

Elastic scattering of relativistic electrons and positrons by atoms is considered in the framework of the static field approximation. The scattering field is expressed as a sum of Yukawa terms to allow the use of various approximations. Accurate phase shifts have been computed by combining Bühring¿s power-series method with the WKB and Born approximations. This combined procedure allows the evaluation of differential cross sections for kinetic energies up to several tens of MeV. Numerical results are used to analyze the validity of several approximate methods, namely the first- and second-order Born approximations and the screened Mott formula, which are frequently adopted as the basis of multiple scattering theories and Monte Carlo simulations of electron and positron transport.

n=1/4 domain-growth universality class: Crossover to the n=1/2 class

The kinetic domain-growth exponent is studied by Monte Carlo simulation as a function of temperature for a nonconserved order-parameter model. In the limit of zero temperature, the model belongs to the n=(1/4 slow-growth unversality class. This is indicative of a temporal pinning in the domain-boundary network of mixed-, zero-, and finite-curvature boundaries. At finite temperature the growth kinetics is found to cross over to the Allen-Cahn exponent n=(1/2. We obtain that the pinning time of the zero-curvature boundary decreases rapidly with increasing temperature.

Molecular dynamics study of the longitudinal modes in disparate-mass binaryliquid mixtures.

A series of molecular dynamics simulations of simple liquid binary mixtures of soft spheres with disparate-mass particles were carried out to investigate the origin of the marked differences between the dynamic structure factors of some liquid binary mixtures such as the Li0.7Mg0.3 and Li0.8Pb0.2 alloys. It is shown that the facility for observing peaks associated with fast-propagating modes in the partial Li-Li dynamic structure factor of Li0.8Pb0.2 should be mainly attributed to the structure of this alloy, which is characterized by an incipient ABAB ordering as found in molten salts. The longitudinal dispersion relations at intermediate wave vectors obtained from the longitudinal current spectra are very similar for the two alloys and reflect the existence of both fast-and slow-propagating modes of kinetic character associated with light and heavy particles, respectively. The influence of the hardness of the repulsive potential cores as well as the composition of the mixture on the longitudinal collective modes is also discussed.

Nonequilibrium patterns and shape fluctuations in reactive membranes

A simple kinetic model of a two-component deformable and reactive bilayer is presented. The two differently shaped components are interconverted by a nonequilibrium reaction, and a phenomenological coupling between local composition and curvature is proposed. When the two components are not miscible, linear stability analysis predicts, and numerical simulations show, the formation of stationary nonequilibrium composition/curvature patterns whose typical size is determined by the reactive process. For miscible components, a linearization of the dynamic equations is performed in order to evaluate the correlation function for shape fluctuations from which the behavior of these systems in micropipet aspiration experiments can be predicted.

Experiments on anisotropic radial viscous fingering

We observe dendritic patterns in fluid flow in an anisotropic Hele-Shaw cell and measure the tip shapes and trajectories of individual dendritic branches under conditions where the pattern growth appears to be dominated by surface tension anisotropy and also under conditions where kinetic effects appear dominant. In each case, the tip position depends on a power law in the time, but the exponent of this power law can vary significantly among flow realizations. Averaging many growth exponents a yields a =0.640.09 in the surface tension dominated regime and a =0.660.09 in the kinetic regime. Restricting the analysis to realizations when a is very close to 0.6 shows great regularity across pattern regimes in the coefficient of the temporal dependence of the tip trajectory.

Generation of dynamic structures in nonequilibrium reactive membranes

We present a nonequlibrium approach for the study of a flexible bilayer whose two components induce distinct curvatures. In turn, the two components are interconverted by an externally promoted reaction. Phase separation of the two species in the surface results in the growth of domains characterized by different local composition and curvature modulations. This domain growth is limited by the effective mixing due to the interconversion reaction, leading to a finite characteristic domain size. In addition to these effects, first introduced in our earlier work [ Phys. Rev. E 71 051906 (2005)], the important new feature is the assumption that the reactive process actively affects the local curvature of the bilayer. Specifically, we suggest that a force energetically activated by external sources causes a modification of the shape of the membrane at the reaction site. Our results show the appearance of a rich and robust dynamical phenomenology that includes the generation of traveling and/or oscillatory patterns. Linear stability analysis, amplitude equations, and numerical simulations of the model kinetic equations confirm the occurrence of these spatiotemporal behaviors in nonequilibrium reactive bilayers.

Traveling waves and nonequilibrium stationary patterns in two-component reactive Langmuir monolayers

A simple kinetic model of a two-component phase-separating Langmuir monolayer with a chemical reaction is proposed. Its analysis and numerical simulations show that nonequilibrium periodic stationary structures and patterns of traveling stripes can spontaneously develop. The nonequilibrium phase diagram of this system is constructed and the properties of the patterns are discussed.

Edelfosine-induced metabolic changes in cancer cells that precede the overproduction of reactive oxygen species and apoptosis

Background: Metabolic flux profiling based on the analysis of distribution of stable isotope tracer in metabolites is an important method widely used in cancer research to understand the regulation of cell metabolism and elaborate new therapeutic strategies. Recently, we developed software Isodyn, which extends the methodology of kinetic modeling to the analysis of isotopic isomer distribution for the evaluation of cellular metabolic flux profile under relevant conditions. This tool can be applied to reveal the metabolic effect of proapoptotic drug edelfosine in leukemia Jurkat cell line, uncovering the mechanisms of induction of apoptosis in cancer cells. Results: The study of 13C distribution of Jukat cells exposed to low edelfosine concentration, which induces apoptosis in ¿5% of cells, revealed metabolic changes previous to the development of apoptotic program. Specifically, it was found that low dose of edelfosine stimulates the TCA cycle. These metabolic perturbations were coupled with an increase of nucleic acid synthesis de novo, which indicates acceleration of biosynthetic and reparative processes. The further increase of the TCA cycle fluxes, when higher doses of drug applied, eventually enhance reactive oxygen species (ROS) production and trigger apoptotic program. Conclusion: The application of Isodyn to the analysis of mechanism of edelfosine-induced apoptosis revealed primary drug-induced metabolic changes, which are important for the subsequent initiation of apoptotic program. Initiation of such metabolic changes could be exploited in anticancer therapy.

Shear thinning and thixotropy of HMHEC and HEC water solutions

Steady state viscosity and thixotropy of hydrophobically modified hydroxyethyl cellulose HMHEC and nonassociative cellulose water solutions are studied. Although all the samples are shear thinning, only the HMHEC is thixotropic, since the migration of hydrophobes to micelles is controlled by diffusion. The Cross model fits steady state curves. The Mewis model, a phenomenological model that proposes that the rate of change of viscosity when the shear rate is suddenly changed is related to the difference between the steady state and current values of viscosity raised to an exponent, fits structure construction experiments when the exponent, n, is estimated to be around 2. The Newtonian assumption used by Mewis cannot be used here, however. This seems to be related to the fact that the thickening is due to bridged micelle formation, which is a slow process, and also to topological constraints and entanglements, which are rapid processes. The kinetic parameter was redefined to kn in order to make it independent of initial conditions. So, kn depends only on how the shear affects the structure. kn reaches a plateau at shear rates too low to produce structure destruction and decreases at higher shear rates.