Physically consistent simulation of mesoscale chemical kinetics: The non-negative FIS-alpha method

Biochemical pathways involving chemical kinetics in medium concentrations (i.e., at mesoscale) of the reacting molecules can be approximated as chemical Langevin equations (CLE) systems. We address the physically consistent non-negative simulation of the CLE sample paths as well as the issue of non-Lipschitz diffusion coefficients when a species approaches depletion and any stiffness due to faster reactions. The non-negative Fully Implicit Stochastic alpha (FIS alpha) method in which stopped reaction channels due to depleted reactants are deleted until a reactant concentration rises again, for non-negativity preservation and in which a positive definite Jacobian is maintained to deal with possible stiffness, is proposed and analysed. The method is illustrated with the computation of active Protein Kinase C response in the Protein Kinase C pathway. (C) 2011 Elsevier Inc. All rights reserved.

Structural characterisation of reaction zone for friction stir welded aluminium-stainless steel joint

In the present investigation, commercially pure Al has been joined with 304 stainless steel (SS) by friction stir welding. The assembly finds widespread application in the field of cryogenics, nuclear, structural industries and domestic appliances. Microstructural characterisation was carried out using scanning and transmission electron microscopes. It has been found that diffusion of Fe, Cr and Ni is substantial within Al; however, diffusion of Al within 304SS is limited. Owing to interdiffusion of chemical species across the bondline, discrete islands of Fe3Al intermetallic form within the reaction zone. The rubbing action of tool over the butting edge of 304SS removed fine particles from 304SS, which were embedded in the stirring zone of Al matrix. Subsequently, austenite underwent phase transformation to ferrite due to large strain within this grain. Fracture path mainly moves through stirring zone of Al alloy under tensile loading; however, in some places, presence of Fe3Al compound has been also found.

Growth Kinetics of Nanoclusters in Solution

We study the growth kinetics of nanoclusters in solution. There are two generic factors that drive growth: (a) reactions that produce the nanomaterial; and (b) diffusion of the nanomaterial due to chemical-potential gradients. We model the growth kinetics of ZnO nanoparticles via coupled dynamical equations for the relevant order parameters, We study this model both analytically and numerically. We find that there is a crossover in thenanocluster growth law: from L(t) similar to t(1/2) in the reaction-controlled regime to L(t) t(1/3) in the diffusion-controlled regime.

The diffusion and relaxation of Gaussian chains in narrow rectangular slits

The confinement of a polymer to volumes whose characteristic linear dimensions are comparable to or smaller than its bulk radius of gyration R-G,R-bulk can produce significant changes in its static and dynamic properties, with important implications for the understanding of single-molecule processes in biology and chemistry. In this paper, we present calculations of the effects of a narrow rectangular slit of thickness d on the scaling behavior of the diffusivity D and relaxation time tau(r) of a Gaussian chain of polymerization index N and persistence length l(0). The calculations are based on the Rouse-Zimm model of chain dynamics, with the pre-averaged hydrodynamic interaction being obtained from the solutions to Stokes equations for an incompressible fluid in a parallel plate geometry in the limit of small d. They go beyond de Gennes' purely phenomenological analysis of the problem based on blobs, which has so far been the only analytical route to the determination of chain scaling behavior for this particular geometry. The present model predicts that D similar to dN(-1) ln(N/d(2)) and tau(r) similar to N(2)d(-1) ln(N/d(2))(-1) in the regime of moderate confinement, where l(0) << d < R-G,R-bulk. The corresponding results for the blob model have exactly the same power law behavior, but contain no logarithmic corrections; the difference suggests that segments within a blob may actually be partially draining and not non-draining as generally assumed.

The diffusion and relaxation of Gaussian chains in narrow rectangular slits

The confinement of a polymer to volumes whose characteristic linear dimensions are comparable to or smaller than its bulk radius of gyration R-G,R-bulk can produce significant changes in its static and dynamic properties, with important implications for the understanding of single-molecule processes in biology and chemistry. In this paper, we present calculations of the effects of a narrow rectangular slit of thickness d on the scaling behavior of the diffusivity D and relaxation time tau(r) of a Gaussian chain of polymerization index N and persistence length l(0). The calculations are based on the Rouse-Zimm model of chain dynamics, with the pre-averaged hydrodynamic interaction being obtained from the solutions to Stokes equations for an incompressible fluid in a parallel plate geometry in the limit of small d. They go beyond de Gennes' purely phenomenological analysis of the problem based on blobs, which has so far been the only analytical route to the determination of chain scaling behavior for this particular geometry. The present model predicts that D similar to dN(-1) ln(N/d(2)) and tau(r) similar to N(2)d(-1) ln(N/d(2))(-1) in the regime of moderate confinement, where l(0) << d < R-G,R-bulk. The corresponding results for the blob model have exactly the same power law behavior, but contain no logarithmic corrections; the difference suggests that segments within a blob may actually be partially draining and not non-draining as generally assumed. (C) 2013 AIP Publishing LLC.

Reaction dynamics under confinement: an exact path integral treatment of a two-stage model of stochastic gene expression

Gene expression in living systems is inherently stochastic, and tends to produce varying numbers of proteins over repeated cycles of transcription and translation. In this paper, an expression is derived for the steady-state protein number distribution starting from a two-stage kinetic model of the gene expression process involving p proteins and r mRNAs. The derivation is based on an exact path integral evaluation of the joint distribution, P(p, r, t), of p and r at time t, which can be expressed in terms of the coupled Langevin equations for p and r that represent the two-stage model in continuum form. The steady-state distribution of p alone, P(p), is obtained from P(p, r, t) (a bivariate Gaussian) by integrating out the r degrees of freedom and taking the limit t -> infinity. P(p) is found to be proportional to the product of a Gaussian and a complementary error function. It provides a generally satisfactory fit to simulation data on the same two-stage process when the translational efficiency (a measure of intrinsic noise levels in the system) is relatively low; it is less successful as a model of the data when the translational efficiency (and noise levels) are high.

Mechanistic Insights into a Non-Classical Diffusion Pathway for the Formation of Hollow Intermetallics: A Route to Multicomponent Hollow Structures

Hollow nanostructures are used for various applications including catalysis, sensing, and drug delivery. Methods based on the Kirkendall effect have been the most successful for obtaining hollow nanostructures of various multicomponent systems. The classical Kirkendall effect relies on the presence of a faster diffusing species in the core; the resultant imbalance in flux results in the formation of hollow structures. Here, an alternate non-Kirkendall mechanism that is operative for the formation of hollow single crystalline particles of intermetallic PtBi is demonstrated. The synthesis method involves sequential reduction of Pt and Bi salts in ethylene glycol under microwave irradiation. Detailed analysis of the reaction at various stages indicates that the formation of the intermetallic PtBi hollow nanoparticles occurs in steps. The mechanistic details are elucidated using control experiments. The use of microwave results in a very rapid synthesis of intermetallics PtBi that exhibits excellent electrocatalytic activity for formic acid oxidation reaction. The method presented can be extended to various multicomponent systems and is independent of the intrinsic diffusivities of the species involved.

Dynamics of the reaction between the free end of a tethered self-avoiding polymer and a flat penetrable surface: A renormalization group study

The average time tau(r) for one end of a long, self-avoiding polymer to interact for the first time with a flat penetrable surface to which it is attached at the other end is shown here to scale essentially as the square of the chain's contour length N. This result is obtained within the framework of the Wilemski-Fixman approximation to diffusion-limited reactions, in which the reaction time is expressed as a time correlation function of a ``sink'' term. In the present work, this sink-sink correlation function is calculated using perturbation expansions in the excluded volume and the polymer-surface interactions, with renormalization group methods being used to resum the expansion into a power law form. The quadratic dependence of tau(r) on N mirrors the behavior of the average time tau(c) of a free random walk to cyclize, but contrasts with the cyclization time of a free self-avoiding walk (SAW), for which tau(r) similar to N-2.2. A simulation study by Cheng and Makarov J. Phys. Chem. B 114, 3321 (2010)] of the chain-end reaction time of an SAW on a flat impenetrable surface leads to the same N-2.2 behavior, which is surprising given the reduced conformational space a tethered polymer has to explore in order to react. (C) 2014 AIP Publishing LLC.

Digestive ripening facilitated atomic diffusion at nanosize regime: Case of AuIn2 and Ag3In intermetallic nanoparticles

Monodisperse colloidal gold-indium (AuIn2) intermetallic nanoparticles have been synthesized from Au and In colloids using the digestive ripening process. Formation of the intermetallic proceeds via digestive ripening facilitated atomic diffusion of Au and In atoms from the Au and In nanoparticles followed simultaneously by their growth in the solution. Optimization of the reaction temperature was found to be crucial for the formation of AuIn2 intermetallic from gold and indium nanoparticles. Transmission electron microscopy revealed the presence of nearly monodisperse nanoparticles of Au and AuIn2 with particle size distribution of 3.7 +/- 1.0 nm and 5.0 +/- 1.6 nm, respectively. UV-visible spectral studies brought out the absence of SPR band in pure AuIn2 intermetallic nanoparticles. Optical study and electron microscopy, in combination with powder X-ray diffraction established phase pure AuIn2 intermetallic nanoparticles unambiguously. The potential of such an unprecedented approach has been further exploited in the synthesis of Ag3In intermetallic nanoparticles with the dimension of less than 10 nm. (C) 2014 Elsevier B.V. All rights reserved.

Tunability of monodispersed intermetallic AuCu nanoparticles through understanding of reaction pathways

Synthesis of size selective monodispersed nanoparticles particularly intermetallic with well-defined compositions represents a challenge. This paper presents a way for the synthesis of intermetallic AuCu nanoparticles as a model system. We show that reduction of Au and Cu precursors is sensitive to the ratio of total molar concentrations of surfactant to metal precursors. A careful design of experiments to understand the kinetics of the reduction process reveals initial formation of seed nanoparticles of pure Au. Reduction of Cu occurs on the surface of the seed followed by diffusion to yield AuCu. This understanding allows us to develop a two step synthesis where the precise size controlled seed of Au nanoparticles produced in the first step is used in the second step reaction mixture as an Au precursor to allow deposition and interdiffusion of Cu that yields size selected AuCu intermetallics of sub 10 nm sizes.

Uninterrupted galvanic reaction for scalable and rapid synthesis of metallic and bimetallic sponges/dendrites as efficient catalysts for 4-nitrophenol reduction

Here, we demonstrate an uninterrupted galvanic replacement reaction (GRR) for the synthesis of metallic (Ag, Cu and Sn) and bimetallic (Cu M, M=Ag, Au, Pt and Pd) sponges/dendrites by sacrificing the low reduction potential metals (Mg in our case) in acidic medium. The acidic medium prevents the oxide formation on Mg surface and facilitates the uninterrupted reaction. The morphology of dendritic/spongy structures is controlled by the volume of acid used for this reaction. The growth mechanism of the spongy/dendritic microstructures is explained by diffusion-limited aggregate model (DLA), which is also largely affected by the volume of acid. The significance of this method is that the yield can be easily predicted, which is a major challenge for the commercialization of the products. Furthermore, the synthesis is complete in 1-2 minutes at room temperature. We show that the sponges/dendrites efficiently act as catalysts to reduce 4-nitrophenol (4-NP) to 4-aminophenol (4-AP) using NaBH4-a widely studied conversion process.

A new dynamic model for study of dislocation pattern formation

Based on the principle given in nonlinear diffusion-reaction dynamics, a new dynamic model for dislocation patterning is proposed by introducing a relaxation time to the relation between dislocation density and dislocation flux. The so-called chemical potential like quantities, which appear in the model can be derived from variation principle for free energy functional of dislocated media, where the free energy density function is expressed in terms of not only the dislocation density itself but also their spatial gradients. The Linear stability analysis on the governing equations of a simple dislocation density shows that there exists an intrinsic wave number leading to bifurcation of space structure of dislocation density. At the same time, the numerical results also demonstrate the coexistence and transition between different dislocation patterns.

A perturbational h4 exponential finite-difference scheme for the convective diffusion equation

A perturbational h4 compact exponential finite difference scheme with diagonally dominant coefficient matrix and upwind effect is developed for the convective diffusion equation. Perturbations of second order are exerted on the convective coefficients and source term of an h2 exponential finite difference scheme proposed in this paper based on a transformation to eliminate the upwind effect of the convective diffusion equation. Four numerical examples including one- to three-dimensional model equations of fluid flow and a problem of natural convective heat transfer are given to illustrate the excellent behavior of the present exponential schemes, the h4 accuracy of the perturbational scheme is verified using double precision arithmetic.

Some theoretical aspects of the simplified Navier-Stokes(SNS) equations

This study deals with the formulation, mathematical property and physical meaning of the simplified Navier-Stokes (SNS) equations. The tensorial SNS equations proposed is the simplest in form and is applicable to flow fields with arbitrary body boundaries. The zones of influence and dependence of the SNS equations, which are of primary importance to numerical solutions, are expounded for the first time from the viewpoint of subcharacteristics. Besides, a detailed analysis of the diffusion process in flow fields shows that the diffusion effect has an influence zone globally windward and an upwind propagation greatly depressed by convection. The maximum upwind influential distance of the viscous effect and the relative importance of the viscous effect in the flow direction to that in the direction normal to the flow are represented by the Reynolds number, which illustrates the conversion of the complete Navier-Stokes (NS) equations to the SNS equations for flows with large Reynolds number.

Characteristics of a CW HF chemical-laser calculated from a simplified two-dimensional model

A two-dimensional simplified model of an HF chemical laser is introduced. Using an implicit finite difference scheme, the solution of two adjacent parallel streams with diffusion mixing and chemical reaction is generated. A contour of mixing and reaction boundary is obtained without presupposition. The distribution of the HF(v) concentrations, gas temperature and the optical small signal gain (alpha sub V, J) on the flowing plane (X, Y) are presented. Compared with the solution solved directly from a set of Navier-Stokes equations, the results of these two methods agree with each other qualitatively. The influences of the different velocity, temperature (T sub 0) and composition of the two streams on the small signal gain after the nozzle exit are investigated. It is interesting that for larger J with a fixed v, the peaks of alpha sub v-T sub 0 profiles move towards higher T sub 0. The computing method is simple and only a short computing time is needed.