Ring closure dynamics for a chemically active polymer.

The principles that underlie the motion of colloidal particles in concentration gradients and the propulsion of chemically-powered synthetic nanomotors are used to design active polymer chains. The active chains contain catalytic and noncatalytic monomers, or beads, at the ends or elsewhere along the polymer chain. A chemical reaction at the catalytic bead produces a self-generated concentration gradient and the noncatalytic bead responds to this gradient by a diffusiophoretic mechanism that causes these two beads to move towards each other. Because of this chemotactic response, the dynamical properties of these active polymer chains are very different from their inactive counterparts. In particular, we show that ring closure and loop formation are much more rapid than those for inactive chains, which rely primarily on diffusion to bring distant portions of the chain in close proximity. The mechanism presented in this paper can be extended to other chemical systems which rely on diffusion to bring reagents into contact for reactions to occur. This study suggests the possibility that synthetic systems could make use of chemically-powered active motion or chemotaxis to effectively carry out complex transport tasks in reaction dynamics, much like those that molecular motors perform in biological systems.

Batroxobin mobilizes circulating endothelial progenitor cells in patients with deep vein thrombosis.

Batroxobin, a thrombin-like enzyme from Bothrops atrox moojeni venom, is associated with the reduction of fibrinogen levels in plasma and the enhancement of anticoagulation and fibrinolysis. In this study, 15 patients with deep vein thrombosis (DVT) achieved successful limb salvage after the administration of batroxobin. We found that the levels of CD34+, CD31+, CD34+/CD31+, and vascular endothelial cadherin (VE-cadherin+) cells had increased in the peripheral blood of patients at 7 days and 14 days after treatment. At 0 day, 7 days, and 14 days, the percentages of CD34+ cells, which are assumed to be hematopoietic stem cells, are 0.39% ± 0.43%, 0.71% ± 0.50%, and 1.11% ± 0.66%, respectively. The levels of CD34+ cells at 14 days are significantly higher than the levels on the first day (P = .004). The levels of CD31+ cells and VE-cadherin+ cells, which represent mature endothelial cells, at 7 days (34.15% ± 11.32%, P = .013; 1.25% ± 1.39%, P = .014) and 14 days (35.21% ± 7.66%, P = .071; 1.85% ± 2.60%, P = .117) were slightly elevated compared with those at 0 day (27.55% ± 8.65%; 0.25 ± 0.39%). The double positive of CD34 and CD31 cells are assumed to be endothelial progenitor cells (EPCs). The levels of CD34+/CD31+ cells at 7 days (0.69% ± 0.50%, P = .001) and 14 days (1.07% ± 0.66%, P = .006) are significantly higher than that on the initial day (0.28% ± 0.30%). The number of CD34+/CD31+ cells significantly increased, indicating that in addition to its role in anticoagulation and fibrinolysis, treatment with batroxobin might simultaneously activate circulating EPCs that might promote the recanalization of the damaged vessel wall.

Catalytic nanomotors: self-propelled sphere dimers.

Experimental and theoretical studies of the self-propelled motional dynamics of a new genre of catalytic sphere dimer, which comprises a non-catalytic silica sphere connected to a catalytic platinum sphere, are reported for the first time. Using aqueous hydrogen peroxide as the fuel to effect catalytic propulsion of the sphere dimers, both quasi-linear and quasi-circular trajectories are observed in the solution phase and analyzed for different dimensions of the platinum component. In addition, well-defined rotational motion of these sphere dimers is observed at the solution-substrate interface. The nature of the interaction between the sphere dimer and the substrate in the aqueous hydrogen peroxide phase is discussed. In computer simulations of the sphere dimer in solution and the solution-substrate interface, sphere-dimer dynamics are simulated using molecular-dynamics methods and solvent dynamics are modeled by mesoscopic multiparticle collision methods taking hydrodynamic interactions into account. The rotational and translational dynamics of the sphere dimer are found to be in good accord with the predictions of computer simulations.

Dynamics of chemically powered nanodimer motors subject to an external force.

The chemically powered self-propelled directed motions of nanodimer motors confined in a rectangular channel and subject to an applied external conservative force are investigated using hybrid molecular dynamics/multiparticle collision dynamics. The influence of factors, such as dimer sizes, chemical reaction type, and the nature of the interaction potentials between dimer monomers and solvent molecules, on the propulsion force and friction constant are examined. The stall force, for which the nanodimer has zero net velocity, and the thermodynamic efficiency of the motor are calculated. Both irreversible and reversible chemical reactions are considered. The simulation results are compared to theoretical predictions which are able to capture the major features of the self-propelled motion.

Self-propelled polymer nanomotors.

Polymer conformations: The picture shows instantaneous conformations of product molecules produced in a catalytic reaction at the head of a polymer chain. The extended (left) and collapsed (right) polymers correspond to chains in good and poor solvent conditions, respectively. The product molecule concentration gradient leads to a directed force on the polymer, so that it functions as a self-propelled nanomotor.

Design of chemically propelled nanodimer motors.

The self-propelled motion of nanodimers fueled by a chemical reaction taking place under nonequilibrium steady state conditions is investigated. The nanodimer consists of a pair of catalytic and chemically inactive spheres, in general with different sizes, with a fixed internuclear separation. The solvent in which the dimer moves is treated at a particle-based mesoscopic level using multiparticle collision dynamics. The directed motion of the dimer can be controlled by adjusting the interaction potentials between the solvent molecules and the dimer spheres, the internuclear separation, and sphere sizes. Dimers can be designed so that the directed motion along the internuclear axis occurs in either direction and is much larger than the thermal velocity fluctuations, a condition needed for such nanodimers to perform tasks involving targeted dynamics.

Multiparticle collision dynamics modeling of viscoelastic fluids.

In order to investigate the rheological properties of viscoelastic fluids by mesoscopic hydrodynamics methods, we develop a multiparticle collision (MPC) dynamics model for a fluid of harmonic dumbbells. The algorithm consists of alternating streaming and collision steps. The advantage of the harmonic interactions is that the integration of the equations of motion in the streaming step can be performed analytically. Therefore, the algorithm is computationally as efficient as the original MPC algorithm for Newtonian fluids. The collision step is the same as in the original MPC method. All particles are confined between two solid walls moving oppositely, so that both steady and oscillatory shear flows can be investigated. Attractive wall potentials are applied to obtain a nearly uniform density everywhere in the simulation box. We find that both in steady and oscillatory shear flows, a boundary layer develops near the wall, with a higher velocity gradient than in the bulk. The thickness of this layer is proportional to the average dumbbell size. We determine the zero-shear viscosities as a function of the spring constant of the dumbbells and the mean free path. For very high shear rates, a very weak "shear thickening" behavior is observed. Moreover, storage and loss moduli are calculated in oscillatory shear, which show that the viscoelastic properties at low and moderate frequencies are consistent with a Maxwell fluid behavior. We compare our results with a kinetic theory of dumbbells in solution, and generally find good agreement.

Periodic orientational motions of rigid liquid-crystalline polymers in shear flow.

The collective periodic motions of liquid-crystalline polymers in a nematic phase in shear flow have, for the first time, been simulated at the particle level by Brownian dynamics simulations. A wide range of parameter space has been scanned by varying the aspect ratio L/D between 10 and 60 at three different scaled volume fractions Lphi/D and an extensive series of shear rates. The influence of the start configuration of the box on the final motion has also been studied. Depending on these parameters, the motion of the director is either characterized as tumbling, kayaking, log-rolling, wagging, or flow-aligning. The periods of kayaking and wagging motions are given by T=4.2(Lphi/D)gamma(-1) for high aspect ratios. Our simulation results are in agreement with theoretical predictions and recent shear experiments on fd viruses in solution. These calculations of elongated rigid rods have become feasible with a newly developed event-driven Brownian dynamics algorithm.

Isotropic-nematic spinodals of rigid long thin rodlike colloids by event-driven Brownian dynamics simulations.

The isotropic-nematic spinodals of solutions of rigid spherocylindrical colloids with various shape anisotropies L/D in a wide range from 10 to 60 are investigated by means of Brownian dynamics simulations. To make these simulations feasible, we developed a new event-driven algorithm that takes the excluded volume interactions between particles into account as instantaneous collisions, but neglects the hydrodynamic interactions. This algorithm is applied to dense systems of highly elongated rods and proves to be efficient. The calculated isotropic-nematic spinodals lie between the previously established binodals in the phase diagram and extrapolate for infinitely long rods to Onsager's [Ann. N. Y. Acad. Sci. 51, 627 (1949)] theoretical predictions. Moreover, we investigate the shear induced shifts of the spinodals, qualitatively confirming the theoretical prediction of the critical shear rate at which the two spinodals merge and the isotropic-nematic phase transition ceases to exist.

Kayaking and wagging of rods in shear flow.

For the first time, we have simulated the periodic collective orientational motions performed by rigid liquid-crystalline polymers with large aspect ratio in the nematic state in shear flow. In order to be able to do so, we developed a new, event-driven Brownian dynamics technique. We present the results of simulations of rods with aspect ratios L/d ranging from 20 to 60 at volume fractions phi given by Lphi/d = 3.5 and 4.5. By studying the path of the director, i.e., the average direction of the rods, we observe kayaking, wagging, flow aligning, and log-rolling type of orbits, depending on the parameters of the simulation and the initial orientation. We find that the tumbling periods depend on Lphi/d and the shear rate but not on the type of motion. Our simulation results qualitatively confirm theoretical predictions and are in good agreement with the experimental measurements of tumbling times of fd viruses.

Brownian dynamics simulations of the self- and collective rotational diffusion coefficients of rigid long thin rods.

Recently a microscopic theory for the dynamics of suspensions of long thin rigid rods was presented, confirming and expanding the well-known theory by Doi and Edwards [The Theory of Polymer Dynamics (Clarendon, Oxford, 1986)] and Kuzuu [J. Phys. Soc. Jpn. 52, 3486 (1983)]. Here this theory is put to the test by comparing it against computer simulations. A Brownian dynamics simulation program was developed to follow the dynamics of the rods, with a length over a diameter ratio of 60, on the Smoluchowski time scale. The model accounts for excluded volume interactions between rods, but neglects hydrodynamic interactions. The self-rotational diffusion coefficients D(r)(phi) of the rods were calculated by standard methods and by a new, more efficient method based on calculating average restoring torques. Collective decay of orientational order was calculated by means of equilibrium and nonequilibrium simulations. Our results show that, for the currently accessible volume fractions, the decay times in both cases are virtually identical. Moreover, the observed decay of diffusion coefficients with volume fraction is much quicker than predicted by the theory, which is attributed to an oversimplification of dynamic correlations in the theory.

Theoretical study on the mechanism of the (3)CH(2) + NO(2) reaction.

The complex doublet potential energy surface of the CH(2)NO(2) system is investigated at the B3LYP/6-31G(d,p) and QCISD(T)/6-311G(d,p) (single-point) levels to explore the possible reaction mechanism of the triplet CH(2) radical with NO(2). Forty minimum isomers and 92 transition states are located. For the most relevant reaction pathways, the high-level QCISD(T)/6-311 + G(2df,2p) calculations are performed at the B3LYP/6-31G(d,p) geometries to accurately determine the energetics. It is found that the top attack of the (3)CH(2) radical at the N-atom of NO(2) first forms the branched open-chain H(2)CNO(2) a with no barrier followed by ring closure to give the three-membered ring isomer cC(H(2))ON-O b that will almost barrierlessly dissociate to product P(1) H(2)CO + NO. The lesser followed competitive channel is the 1,3-H-shift of a to isomer HCN(O)OH c, which will take subsequent cis-trans conversion and dissociation to P(2) OH + HCNO. The direct O-extrusion of a to product P(3) (3)O + H(2)CNO is even much less feasible. Because the intermediates and transition states involved in the above three channels are all lower than the reactants in energy, the title reaction is expected to be rapid, as is consistent with the measured large rate constant at room temperature. Formation of the other very low-lying dissociation products such as NH(2) + CO(2), OH + HNCO and H(2)O + NCO seems unlikely due to kinetic hindrance. Moreover, the (3)CH(2) attack at the end-O of NO(2) is a barrier-consumed process, and thus may only be of significance at very high temperatures. The reaction of the singlet CH(2) with NO(2) is also briefly discussed. Our calculated results may assist in future laboratory identification of the products of the title reaction.

Theoretical study on the mechanism of the 1CHCl + NO reaction.

The complex doublet potential energy surface of the CHClNO system, including 31 minimum isomers and 84 transition states, is investigated at the QCISD(T)/6-311G(d, p)//B3LYP/6-31G(d, p) level in order to explore the possible reaction mechanism of the singlet CHCl with NO. Various possible isomerization and dissociation channels are probed. The initial association between 1CHCl and NO at the terminal N-site can almost barrierlessly lead to the chainlike adducts HClCNO a (a1, a2) followed by the direct Cl-extrusion to product P9 Cl + HCNO, which is the most feasible channel. Much less competitively, a (a1, a2) undergoes a ring-closure leading to the cyclic isomer c-C(HCl)NO d followed by a concerted Cl-shift and N-O cleavage of d to form the branched isomers ClNC(H)O f (f1, f2). Eventually, f (f1, f2) may take a direct H-extrusion to produce P7 H + ClNCO or a concerted 1,2-H-shift and Cl-extrusion to form P1 Cl + HNCO. The low-lying products P2 HCl + NCO, P3 Cl + HOCN, P14 HCO + 3NCl, P6 ClO + HCN, and P13 ClNC + OH may have the lowest yields observed. Our calculations show that the product distributions of the title reaction are quite different from those of the analogous 1CHF + NO reaction, yet are similar to those of another analogous 3CH2 + NO reaction. The similarities and discrepancies among the three reactions are discussed in terms of the substitution effect. The present article may assist in future experimental identification of the product distributions for the title reaction and may be helpful for understanding the halogenated carbene chemistry.