Theory of wormlike polymer chains in confinement

Forced extension of a wormlike chain in the Gibbs and Helmholtz ensembles

A semiflexible polymer can be stretched by either applying a force to it or by fixing the positions of its endpoints. The two approaches generally yield different results and correspond to experiments performed in either the Gibbs or Helmholtz statistical ensembles. Here, we derive the Helmholtz force-extension relationship for the commonly used wormlike-chain model in the strongly stretched regime. By analyzing it in comparison with the Gibbs ensemble result, we show that equivalence between the two relationships is achieved only in the long-chain thermodynamic limit.

Scaling regimes for wormlike chains confined to cylindrical surfaces under tension

We compute the free energy of confinement [Formula: see text] for a wormlike chain (WLC), with persistence length [Formula: see text], that is confined to the surface of a cylinder of radius R under an external tension f using a mean field variational approach. For long chains, we analytically determine the behavior of the chain in a variety of regimes, which are demarcated by the interplay of [Formula: see text], the Odijk deflection length ([Formula: see text]), and the Pincus length ([Formula: see text], with [Formula: see text] being the thermal energy). The theory accurately reproduces the Odijk scaling for strongly confined chains at [Formula: see text], with [Formula: see text]. For moderate values of f, the Odijk scaling is discernible only when [Formula: see text] for strongly confined chains. Confinement does not significantly alter the scaling of the mean extension for sufficiently high tension. The theory is used to estimate unwrapping forces for DNA from nucleosomes.

Intrinsically disordered regions in TRPV2 mediate protein-protein interactions

Transient receptor potential (TRP) ion channels are gated by diverse intra- and extracellular stimuli leading to cation inflow (Na+, Ca2+) regulating many cellular processes and initiating organismic somatosensation. Structures of most TRP channels have been solved. However, structural and sequence analysis showed that ~30% of the TRP channel sequences, mainly the N- and C-termini, are intrinsically disordered regions (IDRs). Unfortunately, very little is known about IDR 'structure', dynamics and function, though it has been shown that they are essential for native channel function. Here, we imaged TRPV2 channels in membranes using high-speed atomic force microscopy (HS-AFM). The dynamic single molecule imaging capability of HS-AFM allowed us to visualize IDRs and revealed that N-terminal IDRs were involved in intermolecular interactions. Our work provides evidence about the 'structure' of the TRPV2 IDRs, and that the IDRs may mediate protein-protein interactions.

Macroscopic analogue to entangled polymers

The entangled structure of polymeric materials is often described as resembling a bowl of spaghetti, swarms of earthworms, or snakes. These analogies not only illustrate the concept, but form the foundation of polymer physics. However, the similarity between these macroscopic, athermal systems and polymers in terms of topology remains uncertain. To better understand this relationship, we conducted an experiment using X-ray tomography to study the structure of arrays of linear rubber bands. We found that, similar to linear polymers, the average number of entanglements increases linearly with the length of the ribbons. Additionally, we observed that entanglements are less frequent near the surface of the container, where there are also more ends, similar to what has been seen in trapped polymers. These findings provide the first experimental evidence supporting the visualization of polymer structures using macroscopic, athermal analogues, confirming the initial intuitive insights of the pioneers of polymer physics.

Depletion interaction mediated by semiflexible polymers

We present a simple mean-field theory to describe the polymer-mediated depletion attraction between colloidal particles that accounts for the polymer's chain stiffness. We find that for fixed polymer radius of gyration and volume fraction, the strength of this attraction increases with increasing chain stiffness in both dilute and semidilute concentration regimes. In contrast, the range of attraction monotonically decreases with chain stiffness in the dilute regime, while it attains a maximum in the semidilute regime. The obtained analytical expressions for the depletion interaction were compared with numerical self-consistent field lattice computations and shown to be in quantitative agreement. From the interaction potential between two spheres, we calculated the second osmotic virial coefficient B₂, which appears to be a convex function of chain stiffness. A minimum of B₂ as a function of chain stiffness was observed both in the numerical self-consistent field computations and the analytical theory. These findings help explain the general observation that semiflexible polymers are more effective depletants than flexible polymers and give insight into the phase behavior of mixtures containing spherical colloids and semiflexible polymers.

Confinement free energy for a polymer chain: Corrections to scaling

Spatial confinement of a polymer chain results in a reduction of conformational entropy. For confinement of a flexible N-mer chain in a planar slit or cylindrical pore (confining dimension D), a blob model analysis predicts the asymptotic scaling behavior ΔF/N ∼ D^-γ with γ ≈ 1.70, where ΔF is the free energy increase due to confinement. Here, we extend this scaling analysis to include the variation of local monomer density upon confinement giving ΔF/N ∼ D^-γ(1 - h(N, D)), where the correction-to-scaling term has the form h ∼ D^y/N^Δ with exponents y = 3 - γ ≈ 1.30 and Δ = 3/γ - 1 ≈ 0.76. To test these scaling predictions, we carry out Wang-Landau simulations of confined and unconfined tangent-hard-sphere chains (bead diameter σ) in the presence of a square-well trapping potential. The fully trapped chain provides a common reference state, allowing for an absolute determination of the confinement free energy. Our simulation results for 32 ≤ N ≤ 1024 and 3 ≤ D/σ ≤ 14 are well-described by the extended scaling relation giving exponents of γ = 1.69(1), y = 1.25(2), and Δ = 0.75(6).

Surface enrichment and interdiffusion in blends of semiflexible polymers of different stiffness

A model for a mixture of two kinds of semiflexible polymers (A and B) with the same chain length (N_A = N_B = 32), but different persistence lengths, confined between parallel planar repulsive walls in a common good solvent is studied by molecular dynamics simulations. In the isotropic phase at low polymer concentrations, both polymers are repelled by the walls, and the system is anisotropic near the walls over a range controlled by the polymer linear dimensions. Close to the concentrations where in the bulk nematic order sets in, precursors of thick nematic layers at the walls are observed, strongly enriched by a stiffer component, which hence is depleted in the center of the slit pore. At larger concentrations, where in the bulk a uniformly mixed nematic phase occurs, the enrichment of B-chains at the walls is rather minor, extending over the scale of the transverse correlation length of concentration fluctuations, which is of the order of a few monomeric diameters only for the present model. In this ordered phase, both self-diffusion and interdiffusion of chains (in the direction perpendicular to the director) are found to be significantly slowed down in comparison to dilute solutions. Since equilibration times scale with the square of the slit thickness, incomplete equilibration is predicted when polymeric coatings on substrate containing polymers differing in stiffness are produced.

Defect patterns of two-dimensional nematic liquid crystals in confinement

A two-dimensional or quasi-two-dimensional nematic liquid crystal refers to a surface-confined system. When such a system is further confined by external line boundaries or excluded from internal line boundaries, the nematic directors form a deformed texture that may display defect points or defect lines, for which winding numbers can be clearly defined. Here, a particular attention is paid to the case when the liquid crystal molecules prefer to form a boundary nematic texture in parallel to the wall surface (i.e., following the homogeneous boundary condition). A general theory, based on geometric argument, is presented for the relationship between the sum of all winding numbers in the system (the total winding number) and the type of confinement angles and curved segments. The conclusion is validated by comparing the theoretical defect rule with existing nematic textures observed experimentally and theoretically in recent years.

Crowding-induced polymer trapping in a channel

In this work, we report an intriguing phenomenon: crowding-induced polymer trapping in a channel. Using Langevin dynamics simulations and analytical calculations, we find that for a polymer confined in a channel, crowding particles can push a polymer into the channel corner through inducing an effective polymer-corner attraction due to the depletion effect. This phenomenon is referred to as polymer trapping. The occurrence of polymer trapping requires a minimum volume fraction of crowders, ϕ^{*}, which scales as ϕ^{*}∼(a_{c}/L_{p})^{1/3} for a_{c}≫a_{m} and ϕ^{*}∼(a_{c}/L_{p})^{1/3}(a_{c}/a_{m})^{1/2} for a_{c}≪a_{m}, where a_{c} is the crowder diameter, a_{m} is the monomer diameter, and L_{p} is the polymer persistence length. For DNA, ϕ^{*} is estimated to be around 0.25 for crowders with a_{c}=2nm. We find that ϕ^{*} also strongly depends on the shape of the channel cross section, and ϕ^{*} is much smaller for a triangle channel than a square channel. The polymer trapping leads to a nearly fully stretched polymer conformation along a channel corner, which may have practical applications, such as full stretching of DNA for the nanochannel-based genome mapping technology.

The stress deformation response influenced by the chain rigidity for mesostructures in diblock copolymers

A self-consistent field theory formalism based on the wormlike chain model is developed to investigate the stress-strain relation for mesostructures in diblock copolymers under the influence of chain rigidity, involving the adjustable simulation cell in the non-orthogonal coordinates by means of optimization of free energy. We elucidate the effect of the chain persistency broadly spanning from the Gaussian chain to the rigid rodlike chain on the elastic response of mesophases that deviate from the initial equilibrium structures. We analytically and numerically demonstrate that our current approach in the long chain limit recovers to the Gaussian-chain-based theory. Being ascribed to the distinct conformational behaviors for flexible chains and rigid rodlike chains, the tensile and compressive stresses applied to lamellae exhibit asymmetric deformation behaviors and the shear stress applied to the initial equilibrium hexagonal cylinders results in noticeable deviations in the shape and spatial arrangement of cylindroids for various chain rigidity values. For the zero stress, in addition, our approach can be straightforwardly utilized to explore the optimal size and shape of the simulation cell in order to achieve a stress free configuration of systems.

Adsorption of Semiflexible Polymers in Cylindrical Tubes

Conformations of wormlike chains in cylindrical pores with attractive walls are explored for varying pore radius and strength of the attractive wall potential by molecular dynamics simulations of a coarse-grained model. Local quantities such as the fraction of monomeric units bound to the surface and the bond-orientational order parameter as well as the radial density distribution are studied, as well as the global chain extensions parallel to the cylinder axis and perpendicular to the cylinder surface. A nonmonotonic convergence of these properties to their counterparts for adsorption on a planar substrate is observed due to the conflict between pore surface curvature and chain stiffness. Also the interpretation of partially adsorbed chains in terms of trains, loops, and tails is discussed.

Cylindrical confinement of solutions containing semiflexible macromolecules: surface-induced nematic order versus phase separation

Solutions of semiflexible polymers confined in cylindrical pores with repulsive walls are studied by Molecular Dynamics simulations for a wide range of polymer concentrations. Both the case where both lengths are of the same order and the case when the persistence length by far exceeds the contour length are considered, and the enhancement of nematic order along the cylinder axis is characterized. With increasing density the character of the surface effect changes from depletion to the formation of a layered structure. For binary 50 : 50 mixtures of the two types of polymers an interplay between surface enrichment of the stiffer component and the isotropic-nematic transition is found, and a phase separated structure with cylindrical symmetry occurs, with the isotropic phase located around the cylinder axis. For melt densities the mixed nematic phase forms at the wall a layer with a screw-like structure of a tilted smectic phase. The observed behavior is tentatively interpreted in terms of the competition of the chain orientational entropy with entropy of mixing and excluded volume due to the wall.

Elasticity of connected semiflexible quadrilaterals

Using the positional-orientational propagator of a semiflexible filament in the weakly bending regime, we analytically calculate the probability densities associated with the fluctuating tip and the corners of a grafted system of connected quadrilaterals. We calculate closed analytic expressions for the probability densities within the framework of the worm-like chain model, which are valid in the weakly bending regime. The probability densities give the physical quantities related to the elasticity of the system such as the force-extension relation in the fixed extension ensemble, the Poisson's ratio and the average of the force exerted to a confining stiff planar wall by the fluctuating tip of the system. Our analysis reveals that the force-extension relations depend on the contour length of the system (material content), the bending stiffness (chemical nature), the geometrical angle and the number of the quadrilaterals, while the Poisson's ratio depends only on the geometrical angle and the number of the quadrilaterals, and is thus a purely geometric property of the system.

The influence of side-chain conformations on the phase behavior of bottlebrush block polymers

A self-consistent field theory based on the wormlike chain model is implemented in the investigation of the self-assembly behavior of bottlebrush block polymers in the formation of a lamellar phase. We utilize the model in which the semi-flexible side chains of two types A and B are grafted at the semi-flexible backbone of type C to mimic the bottlebrush molecule, particularly allowing for the extended chain conformation due to the high grafting density. We examine the positional and orientational probability distribution for the segments along the backbone and side chains as a function of the grafting density and chain flexibility for all blocks, covering a broad regime spanning from the flexible chain to rigid rod chain. This reveals that the persistence length of side chains λ_SC which intrinsically tunes the chain conformation of bottlebrush polymers plays a pivotal role in determining the manner of the local monomer packing in microphase segregation. As an important adjustable factor, λ_SC has a remarkable impact on the backbone extension and then realizes the effective manipulation of the characteristic structural size of self-assembled microstructures, such as the domain spacing and the interfacial width.

Rodlike molecules in extreme confinement

A unique feature of colloid particles and biopolymers is the molecule's intrinsic rigidity characterized by a molecular-level length scale. Under extreme confinement conditions at cellular scales or in nanodevices, these molecules can display orientational ordering accompanied by severe density depletion. Conventional liquid-crystal theories, such as the Oseen-Frank or Landau-de Gennes theories, cannot capture the essential molecular-level properties: the boundary effects, which extend to a distance of the rigidity length scale, and the drastic variations of the inhomogeneous molecular density. Here we show, based on a simple interpretation of the Onsager model, that rodlike molecules in extreme annular confinement produce unusual liquid-crystal defect structures that are independent phases from the patterns usually seen in a weaker confinement environment.

Nematic ordering of worm-like polymers near an interface

The phase behavior of semi-flexible polymers is integral to various contexts, from materials science to biophysics, many of which utilize or require specific confinement geometries as well as the orientational behavior of the polymers. Inspired by collagen assembly, we study the orientational ordering of semi-flexible polymers, modeled as Maier-Saupe worm-like chains, using self-consistent field theory. We first examine the bulk behavior of these polymers, locating the isotropic-nematic transition and delineating the limit of stability of the isotropic and nematic phases. This effort explains how nematic ordering emerges from the isotropic phase and offers insight into how different (e.g., mono-domain vs multi-domain) nematic phases form. We then clarify the influence of planar confinement on the nematic ordering of the polymers. We find that while the presence of a single confining wall does not shift the location of nematic transition, it aligns the polymers in parallel to the wall; above the onset of nematic transition, this preference tends to propagate into the bulk phase. Introducing a second, perpendicular, wall leads to a nematic phase that is parallel to both walls, allowing the ordering direction to be uniquely set by the geometry of the experimental setup. The advantage of wall-confinement is that it can be used to propagate mono-domain nematic phases into the bulk phase.

Statistical ensemble inequivalence for flexible polymers under confinement in various geometries

The problem of statistical ensemble inequivalence for single polymers has been the subject of intense research. In a recent publication, we show that even though the force-extension relation of a free Gaussian chain exhibits ensemble equivalence, confinement to half-space due to tethering to a planar substrate induces significant inequivalence [S. Dutta and P. Benetatos, Soft Matter, 2018, 14, 6857-6866]. In the present article, we extend that work to the conformational response to confining forces distributed over surfaces. We analyze in both the Helmholtz and the Gibbs ensemble the pressure-volume equation of state of a chain in rectangular, spherical, and cylindrical confinement. We especially consider the case of a directed polymer in a cylinder. We also analyze the case of a tethered chain inside a rectangular box, a sphere, and outside a sphere. In general, confinement causes significant ensemble inequivalence. Remarkably, we recover ensemble equivalence in the limit of squashing confinement. We trace the ensemble inequivalence to the persistence of strong fluctuations. Our work may be useful in the interpretation of single molecule experiments and caging phenomena.

Semiflexible Polymers Interacting With Planar Surfaces: Weak versus Strong Adsorption

Semiflexible polymers bound to planar substrates by a short-range surface potential are studied by Molecular Dynamics simulations to clarify the extent to which these chain molecules can be considered as strictly two-dimensional. Applying a coarse-grained bead-spring model, the chain length N and stiffness κ as well as the strength of the adsorption potential ϵwall are varied over a wide range. The excluded-volume (EV) interactions inherent in this model can also be “switched off” to provide a discretized version of the Kratky–Porod wormlike chain model. We study both local order parameters (fraction f of monomers within the range of the potential, bond-orientational order parameter η) and the mean square gyration radius parallel, 〈Rg2〉||, and perpendicular, 〈Rg2〉⊥, to the wall. While for strongly adsorbed chains EV has negligible effect on f and η, we find that 〈Rg2〉|| is strongly affected when the chain contour length exceeds the persistence length. Monomer coordinates in perpendicular (⊥) direction are correlated over the scale of the deflection length which is estimated. It is found that f,η, and 〈Rg2〉⊥ converge to their asymptotic values with 1/N corrections. For both weakly and strongly adsorbed chains, the distribution functions of “loops”, “trains”, and “tails” are analyzed. Some consequences pertaining to the analysis of experiments on adsorbed semiflexible polymers are pointed out.

Extension distribution for DNA confined in a nanochannel near the Odijk regime

DNA confinement in a nanochannel typically is understood via mapping to the confinement of an equivalent neutral polymer by hard walls. This model has proven to be effective for confinement in relatively large channels where hairpin formation is frequent. An analysis of existing experimental data for Escherichia coli DNA extension in channels smaller than the persistence length, combined with an additional dataset for λ-DNA confined in a 34 nm wide channel, reveals a breakdown in this approach as the channel size approaches the Odijk regime of strong confinement. In particular, the predicted extension distribution obtained from the asymptotic solution to the weakly correlated telegraph model for a confined wormlike chain deviates significantly from the experimental distribution obtained for DNA confinement in the 34 nm channel, and the discrepancy cannot be resolved by treating the alignment fluctuations or the effective channel size as fitting parameters. We posit that the DNA-wall electrostatic interactions, which are sensible throughout a significant fraction of the channel cross section in the Odijk regime, are the source of the disagreement between theory and experiment. Dimensional analysis of the wormlike chain propagator in channel confinement reveals the importance of a dimensionless parameter, reflecting the magnitude of the DNA-wall electrostatic interactions relative to thermal energy, which has not been considered explicitly in the prevailing theories for DNA confinement in a nanochannel.

Stochastic curvature of enclosed semiflexible polymers

The conformational states of a semiflexible polymer enclosed in a compact domain of typical size a are studied as stochastic realizations of paths defined by the Frenet equations under the assumption that stochastic "curvature" satisfies a white noise fluctuation theorem. This approach allows us to derive the Hermans-Ullman equation, where we exploit a multipolar decomposition that allows us to show that the positional probability density function is well described by a telegrapher's equation whenever 2a/ℓ_{p}>1, where ℓ_{p} is the persistence length. We also develop a Monte Carlo algorithm for use in computer simulations in order to study the conformational states in a compact domain. In addition, the case of a semiflexible polymer enclosed in a square domain of side a is presented as an explicit example of the formulated theory and algorithm. In this case, we show the existence of a polymer shape transition similar to the one found by Spakowitz and Wang [Phys. Rev. Lett. 91, 166102 (2003)PRLTAO0031-900710.1103/PhysRevLett.91.166102] where in this case the critical persistence length is ℓ_{p}^{*}≃a/8 such that the mean-square end-to-end distance exhibits an oscillating behavior for values ℓ_{p}>ℓ_{p}^{*}, whereas for ℓ_{p}<ℓ_{p}^{*} it behaves monotonically increasing.

Orientationally ordered states of a wormlike chain in spherical confinement

One of the basic characteristics of a linear dsDNA molecule is its persistence length, typically of order 50 nm. The DNA chain inflicts a large energy penalty if it is bent sharply at that length scale. Viruses of bacteria, known as bacteriophages, typically have a dimension of a few tens of nanometers. Yet, it is known that a bacteriophage actively packages viral DNA inside the capsid and ejects it afterwards. Here, adopting a commonly used polymer model known as the wormlike chain, we answer an idealized question: Placing a linear DNA molecule inside a spherical cavity, what ordered states can we derive from known tools in statistical physics? Solving the model in a rigorous field-theory framework, we report a universal phase diagram for four orientationally ordered and disordered states, in terms of two relevant physical parameters.

Statistical Behaviors of Semiflexible Polymer Chains Stretched in Rectangular Tubes

We investigated the statistical behaviors of semiflexible polymer chains, which were simultaneously subjected to force stretching and rectangular tube confinement. Based on the wormlike chain model and Odijk deflection theory, we derived a new deflection length, by using which new compact formulas were obtained for the confinement free energy and force⁻confinement⁻extension relations. These newly derived formulas were justified by numerical solutions of the eigenvalue problem associated with the Fokker⁻Planck governing equation and extensive Brownian dynamics simulations based on the so-called generalized bead-rod (GBR) model. We found that, compared to classical deflection theory, these new formulas were valid for a much more extended range of the confinement size/persistence length ratio and had no adjustable fitting parameters for sufficiently long semiflexible chains in the whole deflection regime.

The entropic force generated by intrinsically disordered segments tunes protein function

Protein structures are dynamic and can explore a large conformational landscape^1,2. Only some of these structural substates are important for protein function (such as ligand binding, catalysis and regulation)^3-5. How evolution shapes the structural ensemble to optimize a specific function is poorly understood^3,4. One of the constraints on the evolution of proteins is the stability of the folded 'native' state. Despite this, 44% of the human proteome contains intrinsically disordered peptide segments greater than 30 residues in length⁶, the majority of which have no known function^7-9. Here we show that the entropic force produced by an intrinsically disordered carboxy terminus (ID-tail) shifts the conformational ensemble of human UDP-α-D-glucose-6-dehydrogenase (UGDH) towards a substate with a high affinity for an allosteric inhibitor. The function of the ID-tail does not depend on its sequence or chemical composition. Instead, the affinity enhancement can be accurately predicted based on the length of the intrinsically disordered segment, and is consistent with the entropic force generated by an unstructured peptide attached to the protein surface^10-13. Our data show that the unfolded state of the ID-tail rectifies the dynamics and structure of UGDH to favour inhibitor binding. Because this entropic rectifier does not have any sequence or structural constraints, it is an easily acquired adaptation. This model implies that evolution selects for disordered segments to tune the energy landscape of proteins, which may explain the persistence of intrinsic disorder in the proteome.

Energy required to pinch a DNA plectoneme

DNA supercoiling plays an important role from a biological point of view. One of its consequences at the supramolecular level is the formation of DNA superhelices named plectonemes. Normally separated by a distance on the order of 10 nm, the two opposite double strands of a DNA plectoneme must be brought closer if a protein or protein complex implicated in genetic regulation is to be bound simultaneously to both strands, as if the plectoneme was locally pinched. We propose an analytic calculation of the energetic barrier, of elastic nature, required to bring closer the two loci situated on the opposed double strands. We examine how this energy barrier scales with the DNA supercoiling. For physically relevant values of elastic parameters and of supercoiling density, we show that the energy barrier is in the k_{B}T range under physiological conditions, thus demonstrating that the limiting step to loci encounter is more likely the preceding plectoneme slithering bringing the two loci side by side.

Self-Avoiding Wormlike Chain Confined in a Cylindrical Tube: Scaling Behavior

Within a confining tube section, the multithreads of a strongly confined, backfolding polymer exert the excluded-volume repulsions on each other and produce physical properties that are very different from those of a confined ideal chain. The conformational properties of a such confined wormlike chain are of fundamental interest and are also practically useful in understanding the DNA confinement problems. Here, the excluded-volume effects are added to the standard wormlike-chain model by a self-consistent field theory. The numerical solutions are examined in light of their scaling properties.

Intradomain phase transitions in flexible block copolymers with self-aligning segments

We study a model of flexible block copolymers (BCPs) in which there is an enlthalpic preference for orientational order, or local alignment, among like-block segments. We describe a generalization of the self-consistent field theory of flexible BCPs to include inter-segment orientational interactions via a Landau-de Gennes free energy associated with a polar or nematic order parameter for segments of one component of a diblock copolymer. We study the equilibrium states of this model numerically, using a pseudo-spectral approach to solve for chain conformation statistics in the presence of a self-consistent torque generated by inter-segment alignment forces. Applying this theory to the structure of lamellar domains composed of symmetric diblocks possessing a single block of "self-aligning" polar segments, we show the emergence of spatially complex segment order parameters (segment director fields) within a given lamellar domain. Because BCP phase separation gives rise to spatially inhomogeneous orientation order of segments even in the absence of explicit intra-segment aligning forces, the director fields of BCPs, as well as thermodynamics of lamellar domain formation, exhibit a highly non-linear dependence on both the inter-block segregation (χN) and the enthalpy of alignment (ε). Specifically, we predict the stability of new phases of lamellar order in which distinct regions of alignment coexist within the single mesodomain and spontaneously break the symmetries of the lamella (or smectic) pattern of composition in the melt via in-plane tilt of the director in the centers of the like-composition domains. We further show that, in analogy to Freedericksz transition confined nematics, the elastic costs to reorient segments within the domain, as described by the Frank elasticity of the director, increase the threshold value ε needed to induce this intra-domain phase transition.

Topological defects in two-dimensional liquid crystals confined by a box

When a spatially uniform system that displays a liquid-crystal ordering on a two-dimensional surface is confined inside a rectangular box, the liquid crystal direction field develops inhomogeneous textures accompanied by topological defects because of the geometric frustrations. We show that the rich variety of nematic textures and defect patterns found in recent experimental and theoretical studies can be classified by the solutions of the rather fundamental, extended Onsager model. This is critically examined based on the determined free energies of different defect states, as functions of a few relevant, dimensionless geometric parameters.

Mean-field theory of active electrolytes: Dynamic adsorption and overscreening

We investigate active electrolytes within the mean-field level of description. The focus is on how the double-layer structure of passive, thermalized charges is affected by active dynamics of constituting ions. One feature of active dynamics is that particles adhere to hard surfaces, regardless of chemical properties of a surface and specifically in complete absence of any chemisorption or physisorption. To carry out the mean-field analysis of the system that is out of equilibrium, we develop the "mean-field simulation" technique, where the simulated system consists of charged parallel sheets moving on a line and obeying active dynamics, with the interaction strength rescaled by the number of sheets. The mean-field limit becomes exact in the limit of an infinite number of movable sheets.

Bistable curvature potential at hyperbolic points of nematic shells

Nematic shells are colloidal particles coated with nematic liquid crystal molecules which may freely glide and rotate on the colloid's surface while keeping their long axis on the local tangent plane. We describe the nematic order on a shell by a unit director field on an orientable surface. Equilibrium fields can then be found by minimising the elastic energy, which in general is a function of the surface gradient of the director field. We learn how to extract systematically out of this energy a fossil component, related only to the surface and its curvatures, which expresses a curvature potential for the molecular torque. At hyperbolic points on the colloid's surface, and only there, the alignment preferred by the curvature potential may fail to be a direction of principal curvature. There the fossil energy becomes bistable.

Conformational Properties of a Back-Folding Wormlike Chain Confined in a Cylindrical Tube

When a semiflexible chain is confined in a narrow cylindrical tube, the formation of a polymer hairpin is a geometrical conformation that accompanies an exponentially large local free energy and, hence, is a relatively rare event. Numerical solutions of the hairpin distribution functions for persistence-length-to-tube-radius ratios over a wide range are obtained in high precision, by using the Green's function approach for the wormlike-chain model. The crossover region between the narrow and moderately narrow tubes is critically investigated in terms of the hairpin free energy, global persistence length, mean hairpin-tip distance from the tube axis, and hairpin-plane orientational properties. Accurate representations of the solutions by simple interpolation formulae are suggested.

Role of growth rate on the orientational alignment of Escherichia coli in a slit

We present experimental data on the nematic alignment of Escherichia coli bacteria confined in a slit, with an emphasis on the effect of growth rate and corresponding changes in cell aspect ratio. Global alignment with the channel walls arises from the combination of local nematic ordering of nearby cells, induced by cell division and the elongated shape of the cells, and the preferential orientation of cells proximate to the side walls of the slit. Decreasing the growth rate leads to a decrease in alignment with the walls, which is attributed primarily to effects of changing cell aspect ratio rather than changes in the variance in cell area. Decreasing confinement also reduces the degree of alignment by a similar amount as a decrease in the growth rate, but the distribution of the degree of alignment differs. The onset of alignment with the channel walls is coincident with the slits reaching their steady-state occupancy and connected to the re-orientation of locally aligned regions with respect to the walls during density fluctuations.

Confinement-driven organization of a histone-complexed DNA molecule in a dense array of nanoposts

The first step in the controlled storage of lengthy DNA molecules is to keep DNA molecules separated while integrated in micrometer-sized space. Herein, we present hybrid Monte Carlo simulations of a histone-complexed DNA (hcDNA) molecule confined in a dense array of nanoposts. Depending on the nanopost dimension, a single, 8.7 kilobase pair hcDNA molecule was either localized and elongated in a single inter-post space surrounded by four nanoposts or spread over several inter-post spaces through passages between two neighboring nanoposts. The conformational change of a hcDNA molecule is interpreted in terms of competitive effects of confinements in the inter-post and passage spaces. We propose that, by elaborately designing nanopost arrays, the competitive confinement effects can be adjusted such that each hcDNA molecule is localized in a single inter-post space, and thereby multiple hcDNA molecules can be physically separated from each other while stored together in the nanopost array.

Theory of Semiflexible Filaments and Networks

We briefly review the recent developments in the theory of individual semiflexible filaments, and of a crosslinked network of such filaments, both permanent and transient. Starting from the free energy of an individual semiflexible chain, models on its force-extension relation and other mechanical properties such as Euler buckling are discussed. For a permanently crosslinked network of filaments, theories on how the network responds to deformation are provided, with a focus on continuum approaches. Characteristic features of filament networks, such as nonlinear stress-strain relation, negative normal stress, tensegrity, and marginal stability are discussed. In the new area of transient filament network, where the crosslinks can be dynamically broken and re-formed, we show some recent attempts for understanding the dynamics of the crosslinks, and the related rheological properties, such as stress relaxation, yield stress and plasticity.

Thermodynamics of a Compressible Maier-Saupe Model Based on the Self-Consistent Field Theory of Wormlike Polymer

This paper presents a theoretical formalism for describing systems of semiflexible polymers, which can have density variations due to finite compressibility and exhibit an isotropic-nematic transition. The molecular architecture of the semiflexible polymers is described by a continuum wormlike-chain model. The non-bonded interactions are described through a functional of two collective variables, the local density and local segmental orientation tensor. In particular, the functional depends quadratically on local density-variations and includes a Maier–Saupe-type term to deal with the orientational ordering. The specified density-dependence stems from a free energy expansion, where the free energy of an isotropic and homogeneous homopolymer melt at some fixed density serves as a reference state. Using this framework, a self-consistent field theory is developed, which produces a Helmholtz free energy that can be used for the calculation of the thermodynamics of the system. The thermodynamic properties are analysed as functions of the compressibility of the model, for values of the compressibility realizable in mesoscopic simulations with soft interactions and in actual polymeric materials.

Maier-Saupe model of polymer nematics: Comparing free energies calculated with Self Consistent Field theory and Monte Carlo simulations

Self Consistent Field (SCF) theory serves as an efficient tool for studying mesoscale structure and thermodynamics of polymeric liquid crystals (LC). We investigate how some of the intrinsic approximations of SCF affect the description of the thermodynamics of polymeric LC, using a coarse-grained model. Polymer nematics are represented as discrete worm-like chains (WLC) where non-bonded interactions are defined combining an isotropic repulsive and an anisotropic attractive Maier-Saupe (MS) potential. The range of the potentials, σ, controls the strength of correlations due to non-bonded interactions. Increasing σ (which can be seen as an increase of coarse-graining) while preserving the integrated strength of the potentials reduces correlations. The model is studied with particle-based Monte Carlo (MC) simulations and SCF theory which uses partial enumeration to describe discrete WLC. In MC simulations the Helmholtz free energy is calculated as a function of strength of MS interactions to obtain reference thermodynamic data. To calculate the free energy of the nematic branch with respect to the disordered melt, we employ a special thermodynamic integration (TI) scheme invoking an external field to bypass the first-order isotropic-nematic transition. Methodological aspects which have not been discussed in earlier implementations of the TI to LC are considered. Special attention is given to the rotational Goldstone mode. The free-energy landscape in MC and SCF is directly compared. For moderate σ the differences highlight the importance of local non-bonded orientation correlations between segments, which SCF neglects. Simple renormalization of parameters in SCF cannot compensate the missing correlations. Increasing σ reduces correlations and SCF reproduces well the free energy in MC simulations.