Minimal knotted polygons in cubic lattices

An implementation of BFACF-style algorithms on knotted polygons in the simple cubic, face centered cubic and body centered cubic lattice is used to estimate the statistics and writhe of minimal length knotted polygons in each of the lattices. Data are collected and analysed on minimal length knotted polygons, their entropy, and their lattice curvature and writhe

Partially directed paths in a wedge

The enumeration of lattice paths in wedges poses unique mathematical challenges. These models are not translationally invariant, and the absence of this symmetry complicates both the derivation of a functional recurrence for the generating function, and solving for it. In this paper we consider a model of partially directed walks from the origin in the square lattice confined to both a symmetric wedge defined by $Y = \pm pX$, and an asymmetric wedge defined by the lines $Y= pX$ and Y=0, where $p > 0$ is an integer. We prove that the growth constant for all these models is equal to $1+\sqrt{2}$, independent of the angle of the wedge. We derive functional recursions for both models, and obtain explicit expressions for the generating functions when $p=1$. From these we find asymptotic formulas for the number of partially directed paths of length $n$ in a wedge when $p=1$. The functional recurrences are solved by a variation of the kernel method, which we call the ``iterated kernel method''. This method appears to be similar to the obstinate kernel method used by Bousquet-Melou. This method requires us to consider iterated compositions of the roots of the kernel. These compositions turn out to be surprisingly tractable, and we are able to find simple explicit expressions for them. However, in spite of this, the generating functions turn out to be similar in form to Jacobi $\theta$-functions, and have natural boundaries on the unit circle.Comment: 26 pages, 5 figures. Submitted to JCT

Forcing Adsorption of a Tethered Polymer by Pulling

We present an analysis of a partially directed walk model of a polymer which at one end is tethered to a sticky surface and at the other end is subjected to a pulling force at fixed angle away from the point of tethering. Using the kernel method, we derive the full generating function for this model in two and three dimensions and obtain the respective phase diagrams. We observe adsorbed and desorbed phases with a thermodynamic phase transition in between. In the absence of a pulling force this model has a second-order thermal desorption transition which merely gets shifted by the presence of a lateral pulling force. On the other hand, if the pulling force contains a non-zero vertical component this transition becomes first-order. Strikingly, we find that if the angle between the pulling force and the surface is beneath a critical value, a sufficiently strong force will induce polymer adsorption, no matter how large the temperature of the system. Our findings are similar in two and three dimensions, an additional feature in three dimensions being the occurrence of a reentrance transition at constant pulling force for small temperature, which has been observed previously for this model in the presence of pure vertical pulling. Interestingly, the reentrance phenomenon vanishes under certain pulling angles, with details depending on how the three-dimensional polymer is modeled

Collapsing lattice animals and lattice trees in two dimensions

We present high statistics simulations of weighted lattice bond animals and lattice trees on the square lattice, with fugacities for each non-bonded contact and for each bond between two neighbouring monomers. The simulations are performed using a newly developed sequential sampling method with resampling, very similar to the pruned-enriched Rosenbluth method (PERM) used for linear chain polymers. We determine with high precision the line of second order transitions from an extended to a collapsed phase in the resulting 2-dimensional phase diagram. This line includes critical bond percolation as a multicritical point, and we verify that this point divides the line into two different universality classes. One of them corresponds to the collapse driven by contacts and includes the collapse of (weakly embeddable) trees, but the other is {\it not yet} bond driven and does not contain the Derrida-Herrmann model as special point. Instead it ends at a multicritical point $P^*$ where a transition line between two collapsed phases (one bond-driven and the other contact-driven) sparks off. The Derrida-Herrmann model is representative for the bond driven collapse, which then forms the fourth universality class on the transition line (collapsing trees, critical percolation, intermediate regime, and Derrida-Herrmann). We obtain very precise estimates for all critical exponents for collapsing trees. It is already harder to estimate the critical exponents for the intermediate regime. Finally, it is very difficult to obtain with our method good estimates of the critical parameters of the Derrida-Herrmann universality class. As regards the bond-driven to contact-driven transition in the collapsed phase, we have some evidence for its existence and rough location, but no precise estimates of critical exponents.Comment: 11 pages, 16 figures, 1 tabl

Simulations of lattice animals and trees

The scaling behaviour of randomly branched polymers in a good solvent is studied in two to nine dimensions, using as microscopic models lattice animals and lattice trees on simple hypercubic lattices. As a stochastic sampling method we use a biased sequential sampling algorithm with re-sampling, similar to the pruned-enriched Rosenbluth method (PERM) used extensively for linear polymers. Essentially we start simulating percolation clusters (either site or bond), re-weigh them according to the animal (tree) ensemble, and prune or branch the further growth according to a heuristic fitness function. In contrast to previous applications of PERM, this fitness function is {\it not} the weight with which the actual configuration would contribute to the partition sum, but is closely related to it. We obtain high statistics of animals with up to several thousand sites in all dimension 2 <= d <= 9. In addition to the partition sum (number of different animals) we estimate gyration radii and numbers of perimeter sites. In all dimensions we verify the Parisi-Sourlas prediction, and we verify all exactly known critical exponents in dimensions 2, 3, 4, and >= 8. In addition, we present the hitherto most precise estimates for growth constants in d >= 3. For clusters with one site attached to an attractive surface, we verify the superuniversality of the cross-over exponent at the adsorption transition predicted by Janssen and Lyssy. Finally, we discuss the collapse of animals and trees, arguing that our present version of the algorithm is also efficient for some of the models studied in this context, but showing that it is {\it not} very efficient for the `classical' model for collapsing animals.Comment: 17 pages RevTeX, 29 figures include

Task shifting and integration of HIV care into primary care in South Africa: The development and content of the streamlining tasks and roles to expand treatment and care for HIV (STRETCH) intervention

Background: Task shifting and the integration of human immunodeficiency virus (HIV) care into primary care services have been identified as possible strategies for improving access to antiretroviral treatment (ART). This paper describes the development and content of an intervention involving these two strategies, as part of the Streamlining Tasks and Roles to Expand Treatment and Care for HIV (STRETCH) pragmatic randomised controlled trial. Methods: Developing the intervention: The intervention was developed following discussions with senior management, clinicians, and clinic staff. These discussions revealed that the establishment of separate antiretroviral treatment services for HIV had resulted in problems in accessing care due to the large number of patients at ART clinics. The intervention developed therefore combined the shifting from doctors to nurses of prescriptions of antiretrovirals (ARVs) for uncomplicated patients and the stepwise integration of HIV care into primary care services. Results: Components of the intervention: The intervention consisted of regulatory changes, training, and guidelines to support nurse ART prescription, local management teams, an implementation toolkit, and a flexible, phased introduction. Nurse supervisors were equipped to train intervention clinic nurses in ART prescription using outreach education and an integrated primary care guideline. Management teams were set up and a STRETCH coordinator was appointed to oversee the implementation process. Discussion: Three important processes were used in developing and implementing this intervention: active participation of clinic staff and local and provincial management, educational outreach to train nurses in intervention sites, and an external facilitator to support all stages of the intervention rollout

Simulations of grafted polymers in a good solvent

We present improved simulations of three-dimensional self avoiding walks with one end attached to an impenetrable surface on the simple cubic lattice. This surface can either be a-thermal, having thus only an entropic effect, or attractive. In the latter case we concentrate on the adsorption transition, We find clear evidence for the cross-over exponent to be smaller than 1/2, in contrast to all previous simulations but in agreement with a re-summed field theoretic $\epsilon$-expansion. Since we use the pruned-enriched Rosenbluth method (PERM) which allows very precise estimates of the partition sum itself, we also obtain improved estimates for all entropic critical exponents.Comment: 5 pages with 9 figures included; minor change

Average Structures of a Single Knotted Ring Polymer

Two types of average structures of a single knotted ring polymer are studied by Brownian dynamics simulations. For a ring polymer with N segments, its structure is represented by a 3N -dimensional conformation vector consisting of the Cartesian coordinates of the segment positions relative to the center of mass of the ring polymer. The average structure is given by the average conformation vector, which is self-consistently defined as the average of the conformation vectors obtained from a simulation each of which is rotated to minimize its distance from the average conformation vector. From each conformation vector sampled in a simulation, 2N conformation vectors are generated by changing the numbering of the segments. Among the 2N conformation vectors, the one closest to the average conformation vector is used for one type of the average structure. The other type of the averages structure uses all the conformation vectors generated from those sampled in a simulation. In thecase of the former average structure, the knotted part of the average structure is delocalized for small N and becomes localized as N is increased. In the case of the latter average structure, the average structure changes from a double loop structure for small N to a single loop structure for large N, which indicates the localization-delocalization transition of the knotted part.Comment: 15 pages, 19 figures, uses jpsj2.cl