Topological phases and curvature-driven pattern formation in cholesteric shells

We study the phase behaviour of cholesteric liquid crystal shells with different geometries. We compare the cases of tangential anchoring and no anchoring at the surface, focussing on the former case, which leads to a competition between the intrinsic tendency of the cholesteric to twist and the anchoring free energy which suppresses it. We then characterise the topological phases arising close to the isotropic-cholesteric transition. These typically consist of quasi-crystalline or amorphous tessellations of the surface by half-skyrmions, which are stable at lower and larger shell sizes, respectively. For ellipsoidal shells, defects in the tessellation couple to a local curvature, and according to the shell size, they either migrate to the poles or distribute uniformly on the surface. For toroidal shells, the variations in the local curvature of the surface stabilise heterogeneous phases where cholesteric or isotropic patterns coexist with hexagonal lattices of half-skyrmions.

Circular Polycatenanes: Supramolecular Structures with Topologically Tunable Properties

Polycatenanes, macrochains of topologically interlocked rings with unique physical properties have recently gained considerable interest in supramolecular chemistry, biology, and soft matter. Most of the work has been, so far, focused on linear chains and on their variety of conformational properties compared to standard polymers. Here we go beyond the linear case and show that, by circularizing such macrochains, one can exploit the topology of the local interlockings to store twist in the system, significantly altering its metric and local properties. Moreover, by properly defining the twist (Tw) and writhe (Wr) of these macrorings we show the validity of a relation equivalent to the Călugăreanu-White-Fuller theorem Tw+Wr=const, originally proved for ribbonlike structures such as double stranded DNA. Our results suggest that circular polycatenanes with storable and tunable twist can form a new category of highly designable multiscale structures with potential applications in supramolecular chemistry and material science.

Cholesteric Shells: Two-Dimensional Blue Fog and Finite Quasicrystals

We study the phase behavior of a quasi-two-dimensional cholesteric liquid crystal shell. We characterize the topological phases arising close to the isotropic-cholesteric transition and show that they differ in a fundamental way from those observed on a flat geometry. For spherical shells, we discover two types of quasi-two-dimensional topological phases: finite quasicrystals and amorphous structures, both made up of mixtures of polygonal tessellations of half-skyrmions. These structures generically emerge instead of regular double twist lattices because of geometric frustration, which disallows a regular hexagonal tiling of curved space. For toroidal shells, the variations in the local curvature of the surface stabilizes heterogeneous phases where cholesteric patterns coexist with hexagonal lattices of half-skyrmions. Quasicrystals and amorphous and heterogeneous structures could be sought experimentally by self-assembling cholesteric shells on the surface of emulsion droplets.

Phase diagrams of confined square lattice linked polygons

The phase diagrams of two models of two confined and dense two-dimensional ring polymers are examined numerically. The ring polymers are modeled by square lattice polygons in a square cavity and are placed to be either unlinked or linked in the plane. The phase diagrams of the two models are found to be a function of the placement of the ring polymers and include multicritical points where first-order and continuous phase boundaries meet. We estimate numerically the critical exponents associated with the phase boundaries and the multicritical points.

Investigating site-selection mechanisms of retroviral integration in supercoiled DNA braids

We theoretically study the integration of short viral DNA in a DNA braid made up by two entwined double-stranded DNA molecules. We show that the statistics of single integration events substantially differ in the straight and buckled, or plectonemic, phase of the braid and are more likely in the latter. We further discover that integration is most likely close to plectoneme tips, where the larger bending energy helps overcome the associated energy barrier and that successive integration events are spatio-temporally correlated, suggesting a potential mechanistic explanation of clustered integration sites in host genomes. The braid geometry we consider provides a novel experimental set-up to quantify integration in a supercoiled substrate in vitro, and to better understand the role of double-stranded DNA topology during this process.

Lamellar ordering, droplet formation and phase inversion in exotic active emulsions

We study numerically the behaviour of a two-dimensional mixture of a passive isotropic fluid and an active polar gel, in the presence of a surfactant favouring emulsification. Focussing on parameters for which the underlying free energy favours the lamellar phase in the passive limit, we show that the interplay between nonequilibrium and thermodynamic forces creates a range of multifarious exotic emulsions. When the active component is contractile (e.g., an actomyosin solution), moderate activity enhances the efficiency of lamellar ordering, whereas strong activity favours the creation of passive droplets within an active matrix. For extensile activity (occurring, e.g., in microtubule-motor suspensions), instead, we observe an emulsion of spontaneously rotating droplets of different size. By tuning the overall composition, we can create high internal phase emulsions, which undergo sudden phase inversion when activity is switched off. Therefore, we find that activity provides a single control parameter to design composite materials with a strikingly rich range of morphologies.

Physical principles of retroviral integration in the human genome

Certain retroviruses, including HIV, insert their DNA in a non-random fraction of the host genome via poorly understood selection mechanisms. Here, we develop a biophysical model for retroviral integration as stochastic and quasi-equilibrium topological reconnections between polymers. We discover that physical effects, such as DNA accessibility and elasticity, play important and universal roles in this process. Our simulations predict that integration is favoured within nucleosomal and flexible DNA, in line with experiments, and that these biases arise due to competing energy barriers associated with DNA deformations. By considering a long chromosomal region in human T-cells during interphase, we discover that at these larger scales integration sites are predominantly determined by chromatin accessibility. Finally, we propose and solve a reaction-diffusion problem that recapitulates the distribution of HIV hot-spots within T-cells. With few generic assumptions, our model can rationalise experimental observations and identifies previously unappreciated physical contributions to retroviral integration site selection.

Entropic elasticity and dynamics of the bacterial chromosome: A simulation study

We study the compression and extension dynamics of a DNA-like polymer interacting with non-DNA binding and DNA-binding proteins, by means of computer simulations. The geometry we consider is inspired by recent experiments probing the compressional elasticity of the bacterial nucleoid (DNA plus associated proteins), where DNA is confined into a cylindrical container and subjected to the action of a "piston"-a spherical bead to which an external force is applied. We quantify the effect of steric interactions (excluded volume) on the force-extension curves as the polymer is compressed. We find that non-DNA-binding proteins, even at low densities, exert an osmotic force which can be a lot larger than the entropic force exerted by the compressed DNA. The trends we observe are qualitatively robust with respect to changes in protein sizes and are similar for neutral and charged proteins (and DNA). We also quantify the dynamics of DNA expansion following removal of the "piston": while the expansion is well fitted by power laws, the apparent exponent depends on protein concentration and protein-DNA interaction in a significant way. We further highlight an interesting kinetic process which we observe during the expansion of DNA interacting with DNA-binding proteins when the interaction strength is intermediate: the proteins bind while the DNA is packaged by the compression force, but they "pop-off" one-by-one as the force is removed, leading to a slow unzipping kinetics. Finally, we quantify the importance of supercoiling, which is an important feature of bacterial DNA in vivo.

Epigenetic Transitions and Knotted Solitons in Stretched Chromatin

The spreading and regulation of epigenetic marks on chromosomes is crucial to establish and maintain cellular identity. Nonetheless, the dynamic mechanism leading to the establishment and maintenance of tissue-specific, epigenetic pattern is still poorly understood. In this work we propose, and investigate in silico, a possible experimental strategy to illuminate the interplay between 3D chromatin structure and epigenetic dynamics. We consider a set-up where a reconstituted chromatin fibre is stretched at its two ends (e.g., by laser tweezers), while epigenetic enzymes (writers) and chromatin-binding proteins (readers) are flooded into the system. We show that, by tuning the stretching force and the binding affinity of the readers for chromatin, the fibre undergoes a sharp transition between a stretched, epigenetically disordered, state and a crumpled, epigenetically coherent, one. We further investigate the case in which a knot is tied along the chromatin fibre, and find that the knotted segment enhances local epigenetic order, giving rise to "epigenetic solitons" which travel and diffuse along chromatin. Our results point to an intriguing coupling between 3D chromatin topology and epigenetic dynamics, which may be investigated via single molecule experiments.

Shear dynamics of an inverted nematic emulsion

Here we study theoretically the dynamics of a 2D and a 3D isotropic droplet in a nematic liquid crystal under a shear flow. We find a large repertoire of possible nonequilibrium steady states as a function of the shear rate and of the anchoring of the nematic director field at the droplet surface. We first discuss homeotropic anchoring. For weak anchoring, we recover the typical behaviour of a sheared isotropic droplet in a binary fluid, which rotates, stretches and can be broken by the applied flow. For intermediate anchoring, new possibilities arise due to elastic effects in the nematic fluid. We find that in this regime the 2D droplet can tilt and move in the flow, or tumble incessantly at the centre of the channel. For sufficiently strong anchoring, finally, one or both of the topological defects which form close to the surface of the isotropic droplet in equilibrium detach from it and get dragged deep into the nematic state by the flow. In 3D, instead, the Saturn ring associated with the normal anchoring disclination line can be deformed and shifted downstream by the flow, but remains always localized in the proximity of the droplet, at least for the parameter range we explored. Tangential anchoring in 2D leads to a different dynamic response, as the boojum defects characteristic of this situation can unbind from the droplet under a weaker shear with respect to the normal anchoring case. Our results should stimulate further experiments with inverted liquid crystal emulsions under shear, as most of the predictions can be testable in principle by monitoring the evolution of liquid crystalline orientation patterns or by tracking the position and shape of the droplet over time.

Electric field controlled columnar and planar patterning of cholesteric colloids

We study how dispersions of colloidal particles in a cholesteric liquid crystal behave under a time-dependent electric field. By controlling the amplitude and shape of the applied field wave, we show that the system can be reproducibly driven out of equilibrium through different kinetic pathways and navigated through a glassylike free energy landscape encompassing many competing metastable equilibria. Such states range from simple Saturn rings to complex structures featuring amorphous defect networks, or stacks of disclination loops. A nonequilibrium electric field can also trigger the alignment of particles into columnar arrays, through defect-mediated force impulses, or their repositioning within a plane. Our results are promising in terms of providing new avenues towards controlled patterning and self-assembly of soft colloid-liquid crystal composite materials.

Interacting elastic lattice polymers: a study of the free energy of globular rings

We introduce and implement a Monte Carlo scheme to study the equilibrium statistics of polymers in the globular phase. It is based on a model of "interacting elastic lattice polymers" and allows a sufficiently good sampling of long and compact configurations, an essential prerequisite to study the scaling behavior of free energies. By simulating interacting self-avoiding rings at several temperatures in the collapsed phase, we estimate both the bulk and the surface free energy. Moreover from the corresponding estimate of the entropic exponent α-2 we provide evidence that, unlike for swollen and Θ-point rings, the hyperscaling relation is not satisfied for globular rings.

Hospitalization rates of complicated pneumococcal community-acquired pneumonia is increasing in Tuscan children

To provide epidemiological data on community-acquired pneumonia (CAP) and complicated CAP, a retrospective study was conducted on a partially vaccinated paediatric population. Data from children hospitalized for CAP in Tuscan hospitals between January 1st, 1999 and December 31st, 2009 were analysed. A total of 5,450 children with CAP were hospitalized. Annual hospitalization rates for CAP did not change significantly over the study period (X2 for trend= 0.652; p=0.419). The total annual hospitalization rate for pneumococcal CAP varied according to age (28.04 per 100,000 children aged less than 5 years, 10.06 per 100,000 children aged 6-12 years and 0.98 per 100,000 children aged greater than13years). Hospitalization rates for pneumococcal CAP increased from12.84 (95 percent CI:7.35-18.34) in 2001 to 45.4 (95 percent CI:35.93-54.90) per 100,000 children aged less than 5 years in 2009 (p less than 0.0001). In addition, a significant increase of hospitalization rates for complicated CAP (from 6.07 in 1999 to 13.66 in 2009 per 100,000 children; P less than 0.0001) and pneumococcal complicated CAP (from 0.19 in 1999 to 3.41 in 2009 per 100,000 children) over the study period were highlighted. Our epidemiological data confirm the decision to introduce the PCV13 vaccine, to satisfy the need to prevent a wider group of pneumococcal serotypes.

TTR fibril formation inhibitors: is there a SAR?

Transthyretin is a homotetrameric protein that carries thyroxine and retinol binding protein in plasma and is associated with a variety of amyloid diseases. One approach to the potential treatment of TTR amyloidosis is the stabilization of the native tetramer, over the dissociative transition state, through the binding of small molecules; this increases the kinetic barrier for tetramer dissociation and prevents protein misfolding. Several molecules discovered through focused screening, or created utilizing the structure-based design, were studied to identify the structural features that could make up for a good candidate drug. In this review, we examine several different chemical classes of TTR fibril formation inhibitors, highlighting the structural modifications that have led to an improvement or to a decrease of their potency and/or selectivity.

Universal properties of knotted polymer rings

By performing Monte Carlo sampling of N-steps self-avoiding polygons embedded on different Bravais lattices we explore the robustness of universality in the entropic, metric, and geometrical properties of knotted polymer rings. In particular, by simulating polygons with N up to 10(5) we furnish a sharp estimate of the asymptotic values of the knot probability ratios and show their independence on the lattice type. This universal feature was previously suggested, although with different estimates of the asymptotic values. In addition, we show that the scaling behavior of the mean-squared radius of gyration of polygons depends on their knot type only through its correction to scaling. Finally, as a measure of the geometrical self-entanglement of the self-avoiding polygons we consider the standard deviation of the writhe distribution and estimate its power-law behavior in the large N limit. The estimates of the power exponent do depend neither on the lattice nor on the knot type, strongly supporting an extension of the universality property to some features of the geometrical entanglement.

Facilitated diffusion on confined DNA

In living cells, proteins combine three-dimensional bulk diffusion and one-dimensional sliding along the DNA to reach a target faster. This process is known as facilitated diffusion and we investigate its dynamics in the physiologically relevant case of confined DNA. The confining geometry and DNA elasticity are key parameters: We find that facilitated diffusion is most efficient inside an isotropic volume and on a flexible polymer. By considering the typical copy numbers of proteins in vivo, we show that the speedup due to sliding becomes insensitive to fine tuning of parameters, rendering facilitated diffusion a robust mechanism to speed up intracellular diffusion-limited reactions. The parameter range we focus on is relevant for in vitro systems and for facilitated diffusion on yeast chromatin.

Bistable defect structures in blue phase devices

Blue phases are liquid crystals made up by networks of defects, or disclination lines. While existing phase diagrams show a striking variety of competing metastable topologies for these networks, very little is known as to how to kinetically reach a target structure, or how to switch from one to the other, which is of paramount importance for devices. We theoretically identify two confined blue phase I systems in which by applying an appropriate series of electric field it is possible to select one of two bistable defect patterns. Our results may be used to realize new generation and fast switching energy-saving bistable devices in ultrathin surface treated blue phase I wafers.

Shearing self-propelled suspensions: arrest of coarsening and suppression of giant density fluctuations

We study the effect of a linear shear flow on a collection of interacting active, self-propelled particles modeled via the Vicsek model. The imposed flow has a dramatic effect on the behavior of the model. We find that in the presence of shear there is no order-disorder transition, and that coarsening of the domains is arrested. Shear also suppresses the so-called giant density fluctuations that are observed in the quiescent limit.

Topological signatures of globular polymers

Simulations in which a globular ring polymer with delocalized knots is separated in two interacting loops by a slipping link, or in two noninteracting globuli by a wall with a hole, show how the minimal crossing number of the knots controls the equilibrium statistics. With slipping link the ring length is divided between the loops according to a simple law, but with unexpectedly large fluctuations. These are suppressed only for unknotted loops, whose length distribution always shows a fast power-law decay. We also discover and explain a topological effect interfering with that of surface tension in the globule translocation through a membrane nanopore.

Biopolymer organization upon confinement

Biopolymers in vivo are typically subject to spatial restraints, either as a result of molecular crowding in the cellular medium or of direct spatial confinement. DNA in living organisms provides a prototypical example of a confined biopolymer. Confinement prompts a number of biophysics questions. For instance, how can the high level of packing be compatible with the necessity to access and process the genomic material? What mechanisms can be adopted in vivo to avoid the excessive geometrical and topological entanglement of dense phases of biopolymers? These and other fundamental questions have been addressed in recent years by both experimental and theoretical means. A review of the results, particularly of those obtained by numerical studies, is presented here. The review is mostly devoted to DNA packaging inside bacteriophages, which is the best studied example both experimentally and theoretically. Recent selected biophysical studies of the bacterial genome organization and of chromosome segregation in eukaryotes are also covered.

Different pulling modes in DNA overstretching: a theoretical analysis

We study the thermally driven denaturation of a double-stranded polymer in the presence of a stretching force via Monte-Carlo simulations. When one strand only is stretched, the denaturation transition is first order, while when both strands are stretched, melting is second order. By revisiting the Poland-Scheraga model for DNA melting, we show that at room temperature, the most likely scenario is that DNA melts as it overstretches. Our results are in general agreement with the most recent experiments and suggest how varying temperature and stretching mode may help settle the question whether S-DNA exists or not.

Supercoil formation in DNA denaturation

We generalize the Poland-Scheraga model to the case of a circular DNA, taking into account the twisting of the two strains around each other. Guided by recent single-molecule experiments on DNA strands, we assume that the torsional stress induced by denaturation enforces the formation of supercoils whose writhe absorbs the linking number expelled by the loops. Our model predicts that when the entropy parameter of a loop satisfies c<or=2, denaturation transition does not take place. On the other hand, for c>2, a first-order denaturation transition is consistent with our model and may take place in the actual system, as in the case with no supercoils. These results are in contrast with other treatments of circular DNA melting where denaturation is assumed to be accompanied by an increase in twist rather than writhe on the bound segments.

Anisotropy of water droplets on single rectangular posts

We report results of extensive experimental and numerical studies of the anisotropy of water drops deposited on single rectangular posts of mesoscopic size sculpted on different materials. Drops of different volume deposited on the top face of the posts assume an elongated shape along the post direction. Systematic investigations show that while the angle measured along the direction parallel to the post does not change, the one measured across them increases monotonically with the drop volume. The difference in these two angles is found to be proportional to the contact line eccentricity even for very elongated drops, regardless of the post size and material. Results obtained with the lattice Boltzmann method are consistent with these observations and indicate useful trends on the evolution of the drop shape with the system main parameters. We argue that drops deposited on single posts having a very sharp profile represent an ideal model system to investigate anisotropic wetting.

Simulations of knotting in confined circular DNA

The packing of DNA inside bacteriophages arguably yields the simplest example of genome organization in living organisms. As an assay of packing geometry, the DNA knot spectrum produced upon release of viral DNA from the P4 phage capsid has been analyzed, and compared to results of simulation of knots in confined volumes. We present new results from extensive stochastic sampling of confined self-avoiding and semiflexible circular chains with volume exclusion. The physical parameters of the chains (contour length, cross section, and bending rigidity) have been set to match those of P4 bacteriophage DNA. By using advanced sampling techniques, involving multiple Markov chain pressure-driven confinement combined with a thermodynamic reweighting technique, we establish the knot spectrum of the circular chains for increasing confinement up to the highest densities for which available algorithms can exactly classify the knots. Compactified configurations have an enclosing hull diameter approximately 2.5 times larger than the P4 caliper size. The results are discussed in relation to the recent experiments on DNA knotting inside the capsid of a P4 tailless mutant. Our investigation indicates that confinement favors chiral knots over achiral ones, as found in the experiments. However, no significant bias of torus over twist knots is found, contrary to the P4 results. The result poses a crucial question for future studies of DNA packaging in P4: is the discrepancy due to the insufficient confinement of the equilibrium simulation or does it indicate that out-of-equilibrium mechanisms (such as rotation by packaging motors) affect the genome organization, hence its knot spectrum in P4?

Shearing active gels close to the isotropic-nematic transition

We study numerically the rheological properties of a slab of active gel close to the isotropic-nematic transition. The flow behavior shows a strong dependence on the sample size, boundary conditions, and on the bulk constitutive curve, which, on entering the nematic phase, acquires an activity-induced discontinuity at the origin. The precursor of this within the metastable isotropic phase for contractile systems (e.g., actomyosin gels) gives a viscosity divergence; its counterpart for extensile suspensions admits instead a shear-banded flow with zero apparent viscosity.

Knot localization in adsorbing polymer rings

We study by Monte Carlo simulations a model of a knotted polymer ring adsorbing onto an impenetrable, attractive wall. The polymer is described by a self-avoiding polygon on the cubic lattice. We find that the adsorption transition temperature, the crossover exponent phi, and the metric exponent nu are the same as in the model where the topology of the ring is unrestricted. By measuring the average length of the knotted portion of the ring, we are able to show that adsorbed knots are localized. This knot localization transition is triggered by the adsorption transition but is accompanied by a less sharp variation of the exponent related to the degree of localization. Indeed, for a whole interval below the adsorption transition, one can not exclude a continuous variation with temperature of this exponent. Deep into the adsorbed phase we are able to verify that knot localization is strong and well described in terms of the flat knot model.

Steady-state hydrodynamic instabilities of active liquid crystals: hybrid lattice Boltzmann simulations

We report hybrid lattice Boltzmann (HLB) simulations of the hydrodynamics of an active nematic liquid crystal sandwiched between confining walls with various anchoring conditions. We confirm the existence of a transition between a passive phase and an active phase, in which there is spontaneous flow in the steady state. This transition is attained for sufficiently "extensile" rods, in the case of flow-aligning liquid crystals, and for sufficiently "contractile" ones for flow-tumbling materials. In a quasi-one-dimensional geometry, deep in the active phase of flow-aligning materials, our simulations give evidence of hysteresis and history-dependent steady states, as well as of spontaneous banded flow. Flow-tumbling materials, in contrast, rearrange themselves so that only the two boundary layers flow in steady state. Two-dimensional simulations, with periodic boundary conditions, show additional instabilities, with the spontaneous flow appearing as patterns made up of "convection rolls." These results demonstrate a remarkable richness (including dependence on anchoring conditions) in the steady-state phase behavior of active materials, even in the absence of external forcing; they have no counterpart for passive nematics. Our HLB methodology, which combines lattice Boltzmann for momentum transport with a finite difference scheme for the order parameter dynamics, offers a robust and efficient method for probing the complex hydrodynamic behavior of active nematics.

Ranking knots of random, globular polymer rings

An analysis of extensive simulations of interacting self-avoiding polygons on cubic lattice shows that the frequencies of different knots realized in a random, collapsed polymer ring decrease as a negative power of the ranking order, and suggests that the total number of different knots realized grows exponentially with the chain length. Relative frequencies of specific knots converge to definite values because the free energy per monomer, and its leading finite size corrections, do not depend on the ring topology, while a subleading correction only depends on the crossing number of the knots.

Size of knots in ring polymers

We give two different, statistically consistent definitions of the length l of a prime knot tied into a polymer ring. In the good solvent regime the polymer is modeled by a self avoiding polygon of N steps on cubic lattice and l is the number of steps over which the knot "spreads" in a given configuration. An analysis of extensive Monte Carlo data in equilibrium shows that the probability distribution of l as a function of N obeys a scaling of the form p(l,N) approximately l(-c)f(l/N(D)) , with c approximately equal to 1.25 and D approximately equal to 1. Both D and c could be independent of knot type. As a consequence, the knot is weakly localized, i.e., <l> approximately N(t) , with t=2-c approximately equal to 0.75 . For a ring with fixed knot type, weak localization implies the existence of a peculiar characteristic length l(nu) approximately N(tnu) . In the scaling approximately N(nu) (nu approximately equal to 0.58) of the radius of gyration of the whole ring, this length determines a leading power law correction which is much stronger than that found in the case of unrestricted topology. The existence of this correction is confirmed by an analysis of extensive Monte Carlo data for the radius of gyration. The collapsed regime is studied by introducing in the model sufficiently strong attractive interactions for nearest neighbor sites visited by the self-avoiding polygon. In this regime knot length determinations can be based on the entropic competition between two knotted loops separated by a slip link. These measurements enable us to conclude that each knot is delocalized (t approximately equal to 1) .

Hydrodynamics and rheology of active liquid crystals: a numerical investigation

We report numerical studies of the hydrodynamics and rheology of an active liquid crystal. We confirm the existence of a transition between a passive and an active phase, with spontaneous flow in steady state. We explore how the velocity profile changes with activity, and we point out the difference in behavior for flow-aligning and tumbling materials. We find that an active material can thicken or thin under a flow, or even exhibit both behaviors as the forcing changes.

Permeative flows in cholesterics: Shear and Poiseuille flows

By using a lattice Boltzmann scheme that solves the Beris-Edwards equations of motion describing liquid-crystal hydrodynamics, we study the response of cholesterics to shear and Poiseuille flows. The geometry we focus on is a flow along the direction of the helical axis, which is known to give rise to permeation. For both shear and Poiseuille flow we find that the boundary conditions on the director field are crucial in determining the rheological properties of the liquid crystal. For helices pinned at the boundaries, a small forcing leads to a large viscosity increase whereas a stronger forcing induces a sharp decrease towards the Newtonian value. This shear thinning behavior is in agreement with experiments and previous analytic results. If, on the other hand, the director is free to rotate at the walls, different behaviors are found depending on the symmetry of the steady-state primary flow. Some of the cases considered are compared to a similar imposed flow but with the helix lying perpendicular to the plates, for which no viscosity increase is observed.

Knotting of random ring polymers in confined spaces

Stochastic simulations are used to characterize the knotting distributions of random ring polymers confined in spheres of various radii. The approach is based on the use of multiple Markov chains and reweighting techniques, combined with effective strategies for simplifying the geometrical complexity of ring conformations without altering their knot type. By these means we extend previous studies and characterize in detail how the probability to form a given prime or composite knot behaves in terms of the number of ring segments N and confining radius R. For 50 < or =N < or =450 we show that the probability of forming a composite knot rises significantly with the confinement, while the occurrence probability of prime knots are, in general, nonmonotonic functions of 1R. The dependence of other geometrical indicators, such as writhe and chirality, in terms of R and N is also characterized. It is found that the writhe distribution broadens as the confining sphere narrows.

Rheology of cholesteric blue phases

Cholesteric blue phases are a spectacular example of disclination line networks. Here we numerically investigate their response to an imposed Poiseuille flow. We show that shear forces bend and twist and can unzip the disclination lines. Under gentle forcing the network opposes the flow and the apparent viscosity is large. With increased forcing we find strong shear thinning corresponding to the disruption of the network. As the viscosity starts to drop, the imposed flow sets the network into motion. Disclinations break up and re-form with their neighbors along the flow.

Biological qualification of blood units: considerations about the effects of sample's handling and storage on stability of nucleic acids

BACKGROUND

In transfusional setting introduction of nucleic amplification technique (NAT) for HBV-DNA, HCV-RNA and HIV-RNA in biological qualification of blood units suggest some problems. At first the opportunity to operate on mini-pool, at second the need to store the samples at +4 degrees C. The authors therefore have tried to estimate the impact of these conditions on the operativity of NAT testing in the transfusional setting.

METHODS

The following parameters has been estimated: distribution of viral-load in untreated subjects, stability of nucleic acids during storage at +4 degrees C, stability of nucleic acids after repeated cycles of freezing and defrosting, robustness of the test to the cross-contamination, definition of the detection-limit (95%). Quantitative tests has been performed by using the following kits: Cobas Amplicor HBV Monitor, Cobas Amplicor HCV Monitor, Cobas Amplicor HIV Monitor; the qualitative tests has been performed by using the following kits: Ampliscreen HBV, Ampliscreen HCV 2,0, Ampliscreen HIV 1,5 all supplied by Roche Molecular System (Brancburg, NJ).

RESULTS

Viral load in untreated subjects showed wide variation for HBV, HCV and HIV. HBV has been demonstrated much stable to the conservation +4 degrees C also until 168 h while for HCV and HIV a greater decrease of the viral-load was observed. For all and three virus the conservation to +4 degrees C until 72 h does not seem to involve meaningful fall in the viral-load. A remarkable reduction of the viral-load has been observed after five cycles of freezing and defrosting. All the tests showed a good robustness to cross-contamination. The detection-limit (95%) was 8 U/ml for HBV, 21 U/ml for HCV and 27 copy/ml for HIV.

CONCLUSIONS

Samples for NAT testing, can be stored until 72 h to +4 degrees C without appreciable lowering of the viral-load. Repeated cycles of changes of state should be avoided. The tests showed a good robustness to cross-contamination. NAT tests for biological qualification of blood units had a minimal sensibility around 50 (copy/unit/ml). In our experience the detection-limit (95%) was 21 U/ml for HCV, 27 copies/ml for HIV, 8 U/ml for HBV. The availability of NAT test for HBV-DNA, HCV-RNA e HIV-RNA, sensitive and reliable, together with epidemiological data, suggest the opportunity to place side by side, in the biological qualification of the blood units, to add the tests for HBV-DNA and HIV-RNA to the test for HCV-RNA mandatory by low, in Italy in the biological qualification of blood units.

Entangled polymers in condensed phases

We present Monte Carlo results on a model of polymers in a condensed phase, over a range of monomer densities. We imagine cutting a cube out of the system. This cube will typically have several polymer molecules running through its interior, and starting and ending on the boundary. These subchains will be mutually entangled and we present a way to assess the extent of entanglement complexity as a function of the monomer density and the number of subchains in the cube. The model is a set of k self-avoiding and mutually avoiding walks, properly embedded in the cube.

Lattice Boltzmann algorithm for three-dimensional liquid-crystal hydrodynamics

We describe a lattice Boltzmann algorithm to simulate liquid-crystal hydrodynamics in three dimensions. The equations of motion are written in terms of a tensor order parameter. This allows both the isotropic and the nematic phases to be considered. Backflow effects and the hydrodynamics of topological defects are naturally included in the simulations, as are viscoelastic effects such as shear-thinning and shear-banding. We describe the implementation of velocity boundary conditions and show that the algorithm can be used to describe optical bounce in twisted nematic devices and secondary flow in sheared nematics with an imposed twist.

Permeative flows in cholesteric liquid crystals

We use lattice Boltzmann simulations to solve the Beris-Edwards equations of motion for a cholesteric liquid crystal subjected to Poiseuille flow along the direction of the helical axis (permeative flow). The results allow us to clarify and extend the approximate analytic treatments currently available. We find that if the cholesteric helix is pinned at the boundaries there is an enormous viscosity increase. If, instead, the helix is free the velocity profile is flattened, but the viscosity is essentially unchanged. We highlight the importance of secondary flows, and, for higher flow velocities, we identify a flow-induced double twist structure in the director field--reminiscent of the texture characteristic of blue phases.

Interplay between shear flow and elastic deformations in liquid crystals

We study shear flow in liquid crystal cells with elastic deformations using a lattice Boltzmann scheme that solves the full, three-dimensional Beris-Edwards equations of hydrodynamics. We consider first twisted and hybrid aligned nematic cells, in which the deformation is imposed by conflicting anchoring at the boundaries. We find that backflow renders the velocity profile non Newtonian, and that the director profile divides into two regions characterized by different director orientations. We next consider a cholesteric liquid crystal, in which a twist deformation is naturally present. We confirm the presence of secondary flow for small shear rates, and are able to follow the dynamical pathway of shear-induced unwinding, for higher shear rates. Finally, we analyze how the coupling between shear and elastic deformation can affect shear banding in an initially isotropic phase. We find that for a nematic liquid crystal, elastic distortions may cause an asymmetry in the dynamics of band formation, whereas for a cholesteric, shear can induce twist in an initially isotropic sample.

RNA denaturation: excluded volume, pseudoknots, and transition scenarios

A lattice model of RNA denaturation which fully accounts for the excluded volume effects among nucleotides is proposed. A numerical study shows that interactions forming pseudoknots must be included in order to get a sharp continuous transition. Otherwise a smooth crossover occurs from the swollen linear polymer behavior to highly ramified, almost compact conformations with secondary structures. In the latter scenario, which is appropriate when these structures are much more stable than pseudoknot links, probability distributions for the lengths of both loops and main branches obey scaling with nonclassical exponents.

Polymer theta-point as a knot delocalization transition

We study numerically the tightness of prime flat knots in a model of self-attracting polymers with excluded volume. We find that these knots are localized in the high temperature swollen regime, but become delocalized in the low temperature globular phase. Precisely at the collapse transition, the knots are weakly localized. Some of our results can be interpreted in terms of the theory of polymer networks, which allows one to conjecture exact exponents for the knot length probability distributions.

Response to a plaque control regimen on different levels of gingival inflammation

BACKGROUND

The purpose of the present study was to evaluate the effectiveness of an oral hygiene regimen in subjects presenting with substantially different severity of plaque-associated gingivitis.

METHODS

The study population was selected from among a large pool of subjects undergoing an experimental gingivitis trial. At completion of the 21-day plaque accumulation period, 2 sub-groups of subjects were identified on the basis of uppermost and lowest quartile for Gingival Index (GI), respectively classified as highly-inflamed (Hinf; n=17; GI: 1.07+/-0.10) and slightly-inflamed (Sinf; n=22; GI: 0.28+/-0.09) groups. An oral hygiene regimen, based on use of amine fluoride/stannous fluoride-containing toothpaste and mouthrinse, was then prescribed for 21 days.

RESULTS

Plaque Index (PI), GI, gingival crevicular fluid (GCF) and Angulated Bleeding Index (AngBI) significantly decreased after treatment in both HInf and SInf groups (p<0.001). However, PI (0.77+/-0.41 vs 0.43+/-0.33, p<0.01), GI (0.23+/-0.30 vs 0.08+/-0.11, p<0.05), GCF (15.23 +/-7.11 vs 7.66+/-2.93, p<0.0000) remained significantly greater in the Hinf group compared to the Sinf group.

CONCLUSIONS

These results suggest that 1) an oral hygiene regimen based on amine/stannous fluoride-containing toothpaste and mouthrinse is effective in reducing plaque-associated gingivitis, regardless of pre-existing severity of gingival inflammation; 2) the level of improvement in gingival status, however, is dependent on the pre-existing severity of the inflammatory condition.

Interstrand distance distribution of DNA near melting

The distance distribution between complementary base pairs of the two strands of a DNA molecule is studied near the melting transition. Scaling arguments are presented for a generalized Poland-Scheraga-type model that includes self-avoiding interactions. At the transition temperature and for a large distance r, the distribution decays as 1/r(kappa) with kappa=1+(c-2)/nu. Here nu is the self-avoiding walk correlation length exponent and c is the exponent associated with the entropy of an open loop in the chain. Results for the distribution function just below the melting point are also presented. Numerical simulations that fully take into account the self-avoiding interactions are in good agreement with the scaling approach.