Dynamics of Diblock Copolymers in Dilute Solutions

We consider the dynamics of freely translating and rotating diblock (A-B), Gaussian copolymers, in dilute solutions. Using the multiple scattering technique, we have computed the diffusion and the friction coefficients D_AB and Zeta_AB, and the change Eta_AB in the viscosity of the solution as functions of x = N_A/N and t = l_B/l_A, where N_A, N are the number of segments of the A block and of the whole copolymer, respectively, and l_A, l_B are the Kuhn lengths of the A and B blocks. Specific regimes that maximize the efficiency of separation of copolymers with distinct "t" values, have been identified.Comment: 20 pages Revtex, 7 eps figures, needs epsf.tex and amssymb.sty, submitted to Macromolecule

Efficient Set Sharing Using ZBDDs

Set sharing is an abstract domain in which each concrete object is represented by the set of local variables from which it might be reachable. It is a useful abstraction to detect parallelism opportunities, since it contains definite information about which variables do not share in memory, i.e., about when the memory regions reachable from those variables are disjoint. Set sharing is a more precise alternative to pair sharing, in which each domain element is a set of all pairs of local variables from which a common object may be reachable. However, the exponential complexity of some set sharing operations has limited its wider application. This work introduces an efficient implementation of the set sharing domain using Zero-suppressed Binary Decision Diagrams (ZBDDs). Because ZBDDs were designed to represent sets of combinations (i.e., sets of sets), they naturally represent elements of the set sharing domain. We show how to synthesize the operations needed in the set sharing transfer functions from basic ZBDD operations. For some of the operations, we devise custom ZBDD algorithms that perform better in practice. We also compare our implementation of the abstract domain with an efficient, compact, bit set-based alternative, and show that the ZBDD version scales better in terms of both memory usage and running time

Microphase separation in polyelectrolytic diblock copolymer melt : weak segregation limit

We present a generalized theory of microphase separation for charged-neutral diblock copolymer melt. Stability limit of the disordered phase for salt-free melt has been calculated using Random Phase Approximation (RPA) and self-consistent field theory (SCFT). Explicit analytical free energy expressions for different classical ordered microstructures (lamellar, cylinder and sphere) are presented. We demonstrate that chemical mismatch required for the onset of microphase separation ($\chi^{\star} N$) in charged-neutral diblock melt is higher and the period of ordered microstructures is lower than those for the corresponding neutral-neutral diblock system. Theoretical predictions on the period of ordered structures in terms of Coulomb electrostatic interaction strength, chain length, block length, and the chemical mismatch between blocks are presented. SCFT has been used to go beyond the stability limit, where electrostatic potential and charge distribution are calculated self-consistently. Stability limits calculated using RPA are in perfect agreement with the corresponding SCFT calculations. Limiting laws for stability limit and the period of ordered structures are presented and comparisons are made with an earlier theory. Also, transition boundaries between different morphologies have been investigated

Anomalous Dynamics of Translocation

We study the dynamics of the passage of a polymer through a membrane pore (translocation), focusing on the scaling properties with the number of monomers $N$. The natural coordinate for translocation is the number of monomers on one side of the hole at a given time. Commonly used models which assume Brownian dynamics for this variable predict a mean (unforced) passage time $\tau$ that scales as $N^2$, even in the presence of an entropic barrier. However, the time it takes for a free polymer to diffuse a distance of the order of its radius by Rouse dynamics scales with an exponent larger than 2, and this should provide a lower bound to the translocation time. To resolve this discrepancy, we perform numerical simulations with Rouse dynamics for both phantom (in space dimensions $d=1$ and 2), and self-avoiding (in $d=2$) chains. The results indicate that for large $N$, translocation times scale in the same manner as diffusion times, but with a larger prefactor that depends on the size of the hole. Such scaling implies anomalous dynamics for the translocation process. In particular, the fluctuations in the monomer number at the hole are predicted to be non-diffusive at short times, while the average pulling velocity of the polymer in the presence of a chemical potential difference is predicted to depend on $N$.Comment: 9 pages, 9 figures. Submitted to Physical Review

Anomalous Dynamics of Forced Translocation

We consider the passage of long polymers of length N through a hole in a membrane. If the process is slow, it is in principle possible to focus on the dynamics of the number of monomers s on one side of the membrane, assuming that the two segments are in equilibrium. The dynamics of s(t) in such a limit would be diffusive, with a mean translocation time scaling as N^2 in the absence of a force, and proportional to N when a force is applied. We demonstrate that the assumption of equilibrium must break down for sufficiently long polymers (more easily when forced), and provide lower bounds for the translocation time by comparison to unimpeded motion of the polymer. These lower bounds exceed the time scales calculated on the basis of equilibrium, and point to anomalous (sub-diffusive) character of translocation dynamics. This is explicitly verified by numerical simulations of the unforced translocation of a self-avoiding polymer. Forced translocation times are shown to strongly depend on the method by which the force is applied. In particular, pulling the polymer by the end leads to much longer times than when a chemical potential difference is applied across the membrane. The bounds in these cases grow as N^2 and N^{1+\nu}, respectively, where \nu is the exponent that relates the scaling of the radius of gyration to N. Our simulations demonstrate that the actual translocation times scale in the same manner as the bounds, although influenced by strong finite size effects which persist even for the longest polymers that we considered (N=512).Comment: 13 pages, RevTeX4, 16 eps figure

An overview of the ciao multiparadigm language and program development environment and its design philosophy

We describe some of the novel aspects and motivations behind the design and implementation of the Ciao multiparadigm programming system. An important aspect of Ciao is that it provides the programmer with a large number of useful features from different programming paradigms and styles, and that the use of each of these features can be turned on and off at will for each program module. Thus, a given module may be using e.g. higher order functions and constraints, while another module may be using objects, predicates, and concurrency. Furthermore, the language is designed to be extensible in a simple and modular way. Another important aspect of Ciao is its programming environment, which provides a powerful preprocessor (with an associated assertion language) capable of statically finding non-trivial bugs, verifying that programs comply with specifications, and performing many types of program optimizations. Such optimizations produce code that is highly competitive with other dynamic languages or, when the highest levéis of optimization are used, even that of static languages, all while retaining the interactive development environment of a dynamic language. The environment also includes a powerful auto-documenter. The paper provides an informal overview of the language and program development environment. It aims at illustrating the design philosophy rather than at being exhaustive, which would be impossible in the format of a paper, pointing instead to the existing literature on the system

What drives the translocation of stiff chains?

We study the dynamics of the passage of a stiff chain through a pore into a cell containing particles that bind reversibly to it. Using Brownian Molecular Dynamics simulations we investigate the mean-first-passage time as a function of the length of the chain inside, for different concentrations of binding particles. As a consequence of the interactions with these particles, the chain experiences a net force along its length whose calculated value from the simulations accounts for the velocity at which it enters the cell. This force can in turn be obtained from the solution of a generalized diffusion equation incorporating an effective Langmuir adsorption free energy for the chain plus binding particles. These results suggest a role of binding particles in the translocation process which is in general quite different from that of a Brownian ratchet. Furthermore, non-equilibrium effects contribute significantly to the dynamics, \emph{e.g.}, the chain often enters the cell faster than particle binding can be saturated, resulting in a force several times smaller than the equilibrium value.Comment: 7 pages, 4 figure

Electronic structure of strongly correlated d-wave superconductors

We study the electronic structure of a strongly correlated d-wave superconducting state. Combining a renormalized mean field theory with direct calculation of matrix elements, we obtain explicit analytical results for the nodal Fermi velocity, v_F, the Fermi wave vector, k_F, and the momentum distribution, n_k, as a function of hole doping in a Gutzwiller projected d-wave superconductor. We calculate the energy dispersion, E_k, and spectral weight of the Gutzwiller-Bogoliubov quasiparticles, and find that the spectral weight associated with the quasiparticle excitation at the antinodal point shows a non monotonic behavior as a function of doping. Results are compared to angle resolved photoemission spectroscopy (ARPES) of the high temperature superconductors.Comment: final version, comparison to experiments added, 4+ pages, 4 figure

Single Stranded DNA Translocation Through A Nanopore: A Master Equation Approach

We study voltage driven translocation of a single stranded (ss) DNA through a membrane channel. Our model, based on a master equation (ME) approach, investigates the probability density function (pdf) of the translocation times, and shows that it can be either double or mono-peaked, depending on the system parameters. We show that the most probable translocation time is proportional to the polymer length, and inversely proportional to the first or second power of the voltage, depending on the initial conditions. The model recovers experimental observations on hetro-polymers when using their properties inside the pore, such as stiffness and polymer-pore interaction.Comment: 7 pages submitted to PR