Next Article in Journal
Evaluation of Glycerylphytate Crosslinked Semi- and Interpenetrated Polymer Membranes of Hyaluronic Acid and Chitosan for Tissue Engineering
Previous Article in Journal
Dielectric Elastomers UV-Cured from Poly(dimethylsiloxane) Solution in Vinyl Acetate
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Topological Constraints on Penetration Structures of Semi-Flexible Ring Polymers

1
Department of Physics, Zhejiang University, Hangzhou 310027, China
2
Department of Physics, Wenzhou University, Wenzhou 325035, China
*
Authors to whom correspondence should be addressed.
Submission received: 7 September 2020 / Revised: 3 November 2020 / Accepted: 8 November 2020 / Published: 11 November 2020
(This article belongs to the Special Issue Topology and Dynamics of Ring Polymers)

Abstract

:
The effects of topological constraints on penetration structures of semi-flexible ring polymers in a melt are investigated using molecular dynamics simulations, considering simultaneously the effects of the chain stiffness. Three topology types of rings are considered: 01-knot (the unknotted), 31-knot and 61-knot ring polymers, respectively. With the improved algorithm to detect and quantify the inter-ring penetration (or inter-ring threading), the degree of ring threading does not increase monotonously with the chain stiffness, existing a peak value at the intermediate stiffness. It indicates that rings interpenetrate most at intermediate stiffness where there is a balance between coil expansion (favoring penetrations) and stiffness (inhibiting penetrations). Meanwhile, the inter-ring penetration would be suppressed with the knot complexity of the rings. The analysis of effective potential between the rings provides a better understanding for this non-monotonous behavior in inter-ring penetration.

Graphical Abstract

1. Introduction

A ring polymer is formed by the simple process of joining the ends of a linear polymer chain, corresponding to the so called 01-knot or trivial knot, that is, an unknotted ring polymer. For the semi-flexible rings in the range of the crossover from the rod-like limit to the random coil one, the conformations show a strong anisotropic character relying on the stiffness. As is well known, DNA often exhibits a knotted ring conformation in bacteria and phages [1]. Now, synthetic ring polymers display very complex knotted configurations, such as 31-knot, 41-knot, 51-knot, 61-knot and so forth in Conway’s notation (e.g., Nk-knot denotes the kth type of knotted ring with N crossings) [2]. The presence of knots in a ring polymer can be viewed as a self-entanglement phenomenon [3]. The topological constraints of rings has a dramatic effect on the physical properties with respect to their linear counterparts due to the change in the conformational degrees of freedom [4,5]. The most prominent examples are their different scaling behaviors [6,7,8,9,10] and rheological properties [11,12]. First, it has been demonstrated that the size of isolated rings scales with the polymerization degree (N) as Dg~NVF, where Dg is the diameter of gyration and VF ≈ 0.588 is the Flory exponent, respectively. However, the size of linear chain scales as Dg~N0.5 [7]. Second, no free ends implies that rings do not relax via common reptation which is primarily responsible for linear chains [13]. Unlike the linear chains, the ring polymer melts exhibit self-similar dynamics, thereby yielding a power-law stress relaxation instead of the rubber-like plateau of linear melts [14]. Additionally, the topological constraints dramatically reduces its attainable states in the phase space, resulting in a repulsive effective interaction between two ring chains, while it vanishes between two linear chains or a linear chain and a ring chain [15,16,17].
Absent in systems of linear chains, inter-ring penetration (or inter-ring threading) is a unique feature of ring polymers but the effects of threading is still controversial until now. Earlier studies [18,19] neglect the threading and assume rings to adopt a double-folded annealed branched structure. This so-called annealed tree model [20,21] shows a great agreement with simulations whether on static or dynamic properties so far. However, studies by Michieletto et al. [22,23] and Lee et al. [24,25] believe in the profound effects of the ring threading, responsible for the observed “topologically-driven glassy” behavior. To address this controversy, it is extremely necessary to have a comprehensive understanding of penetration structures of semi-flexible rings. It is quite essential for theoretical simulation and also for experiment of bio-packing [26,27], such as observed in circular DNA of bacteriophage viruses.
In this work, we focus on the penetration structures of three knot types of rings, that is, the 01-knot, the 31-knot and the 61-knot rings [2]. We define the degree of penetration and monitor the internal structure of penetrated rings in detail by changing the complexity of the knot, the bending energy of rings and the number density of system. Meanwhile, we also derive the effective potential between penetrated rings to support our findings about the effects of the topological constraints and chain stiffness on inter-ring penetration. Our results can provide a deeper understanding about the penetration structure of rings, especially for the case of semi-flexible rings.

2. Simulation Model and Method

2.1. Model

A bead-spring model by Kremer and Grest [28] is used to model the non-concatenated ring polymer, each consisting of N spherical monomers with a monomer diameter of σ and a mass of m. Bonded monomers interact via the finitely extensible nonlinear elastic (FENE) potential [29]
U F E N E ( r ) = K R 0 2 2 l n [ 1 ( r R 0 ) 2 ] ,
where r is the distance between two neighboring monomers. K = 30kBT/σ2 is a spring constant and R0 = 1.5σ is a finite extensibility to avoid chain crossing, where kB is the Boltzmann constant and T is the temperature. The chain stiffness is introduced by means of a bending potential between adjacent bonds given by
U b e n d i n g ( θ ) = K b ( 1 c o s ( θ θ 0 ) ) ,  
where θ is the angle between two consecutive bonds with θ 0 = (N-2)*π/N and Kb represents the bending energy. For any two monomers, they also interact through a shift and cutoff Lennard-Jones (LJ) potential
U L J ( r ) = {   4 ϵ [ ( σ r ) 12 ( σ r ) 6 + 1 4 ]             r < r C 0                           r > r C ,
where r c is the cut off distance and fixed at 2 1 / 6 σ and ε is chosen to be ε = 1 k B T .

2.2. Molecule Dynamics Simulation

We perform molecule dynamics (MD) simulations by using the large-scale atomic/molecular massively parallel simulator (LAMMPS) [30]. Reduced units ( ε = 1 , σ = 1 , m = 1 and τ0 = (2/kBT)1/2 = 1 are chosen to be the units of energy, length, mass and time, respectively) and the timestep is τ = 0.001 τ 0 . The simulation box size is fixed at 50σ × 50σ × 50σ and the periodic boundary conditions are performed in NVT ensemble (canonical ensemble). All ring polymers are placed randomly in our simulation box initially. The corresponding dynamics does not allow for chain crossings to make sure that the topology is preserved. To get a high concentration of rings, we set a scaling factor big enough. Thousands of steps are used to obtain the local equilibrated conformation after each step of scaling, making sure that the scaling does not affect the topology of rings. After that, an additional equilibration run (up to 9   × 108 steps) is performed to get the steady penetration structures of rings.
In our simulation, three topology types or knotted rings are considered: 01-knot (the unknotted or trivial knot), 31-knot (the trefoil knot) and 61-knot ring polymers, named in Conway’s notation [2]. The typical snapshots for latter two knots are shown in Figure 1. Here the bending energy Kb ranges from Kb = 0 to Kb = 100 in the units of ε, corresponding to flexible, semi-flexible and stiff ring polymers, respectively. Intermediate Kb values correspond the semi-flexible rings. The number density is defined as ρ = M*N/L3, where M is the number of rings and N is the number of monomers per ring. Firstly, ρ = 0.1 with N=128 is considered, that is, M = 98. Then the number density ρ varies in the range of 0.05~0.4 with N = 128, that is, M= 49~392.
The effective potential of mean force (PMF) between a pair of rings can be calculated directly from NVT simulations given by References [31,32,33,34]
β V e f f ( r ) = ln ρ ( r ) ,
where β = ( k B T ) 1 and ρ ( r ) is the probability distribution of the centers-of-mass of rings. Since the range of interaction is finite (of the order of D g , the average diameter of gyration), we do not need to calculate those distances which are longer than D g . While for the short distance ( r < D g ), due to the topology repulsion, there are large energy barriers which prevent an effective exploration of the configurational space within the available computer time. To overcome this, we use the umbrella sampling [35,36] to calculate ρ ( r ) where the biased potential is a simple harmonic spring
V b i a s ( r ) = k j 2 ( r r j ) 2 ,
where k j = 4.0 ε / σ 2 and r j stands for the free length of the spring at different windows. We perform the simulations starting from r j = 0 and increase it up to a maximum value which is much larger than D g for different rings to get the biased distributions ρ j ( b ) ( r ) . To eliminate the influence of the harmonic potential on total probability distribution ρ ( r ) , we have to employ the weighted histogram analysis method (WHAM) [33,34,35,36,37] to combine the biased distributions ρ j ( b ) ( r ) at each window. In order to achieve that, we need derive ρ ( r ) in terms of unbiased densities ρ j ( u ) ( r ) at each simulation window [34,35]
ρ ( r ) = j = 1 N p j ( r ) ρ j ( u ) ( r ) ,
with
ρ j ( u ) ( r ) = e β ( V j ( r ) f j ) ρ j ( b ) ( r ) ,
where V j is the biased potential and f j is free energy parameter arising from the addition of the biased potential. The normalization condition requires the sum of p j ( r ) is equal to 1 and they are chosen so as to minimize the statistical error. Then we can get
ρ ( r ) = j = 1 N c j ( r ) ρ j ( b ) ( r ) ,
with coefficient
c j ( r ) = 1 / j N e β ( V j ( r ) f j ) .
Up to now, we already get the total probability distribution ρ ( r ) except that f j are unknown. In fact, f j can be obtained self-consistently by
e f k = i = 1 N j = 1 n i e β V k ( r i , l ) j = 1 N n j e β ( V j ( r i , l ) f j ) ,
where n j is the number of sampled conformations at jth window. r i , l is the reaction coordinate value for conformation snap l in the simulation window i.

2.3. KMT Algorithm

To define whether a chain is penetrated by others, that is, a definition of threading, we introduce the improved KMT algorithm firstly suggested by Koniaris and Muthukumar [38], developed independently by Taylor [39] and referred as KTM reduction by Virnau [40]. As shown in Figure 2, assuming that the ring A is penetrated by another ring B. The KMT algorithm divides the ring A into a series of consecutive triangles of adjacent monomers, such as 123, 345 and 567. The cross denotes the intersection point of which ring A is penetrated by ring B in the case of reduction shown in Figure 2. After the first reduction, monomers 4 and 6 are removed since neither triangle 345 nor 567 is intersected by any part of ring B, while monomer 5 is retained. Through enough iterations of reduction, ring A is reduced into a simple triangle 123. The surface of the ring A is composed of triangle 123, while other triangles are removed at each reduction (e.g., 345 and 567 at the first reduction). Among all the possible surfaces by different routines of reduction, we can find the minimal surface [25,41] and the corresponding final triangle (e.g., the case of Figure 2b) to define the threading. With the definition of threading produced by a pair of rings, we call them one passive ring (A) and one active ring (B) respectively. The ring A is passively threaded by ring B, which is actively threading the ring A [21]. In the following, all of the average statistics goes through all rings.

3. Results and Discussion

We firstly focus on the degree of ring threading with different knot types and chain stiffness. An average threading number Pth is defined as the number of rings penetrated by a single chain. Here three knot types are considered: 01-knot, 31-knot and 61-knot rings respectively. The Pth varied by bending energy Kb for three knot types is displayed in Figure 3. It can be seen that Pth drops rapidly from 01-knot to 31-knot and even to 61-knot rings, due to the decrease of the free space that can be penetrated by other rings due to the internal topology. For the given ring topology structure, for example, 01-knot ring, as Kb increases, Pth increases rapidly first and then decreases slightly, exhibiting a maximum value Pth-max corresponding to an intermediate stiffness. A similar Kb dependence of Pth is presented for other two knots. For flexible rings, the knot has little influence on Pth as rings shrink into crumpled globular conformations, highly unfavorable for penetration. With the increase of bending energy Kb, rings can expand and adopt open-up configurations that facilitate inter-ring penetration, in great agreement with the results by Bernabei et.al [42]. However, the intrinsic mechanism on this non-monotonous dependence on bending energy Kb should be further explored in the following. As well known, the persistence length lp is an important parameter to describe the conformation of semiflexible polymer. The ratio of lp/L depends on the bending energy (Kb), chain length (N) and the chain topological structure is given in Figure S1 (see Figure S1 in Supporting Information).
In detail, we turn to detect the internal structure of penetrated rings to illustrate the expansion of the rings. The probability density distributions P(r) of the centers of mass between two penetrated rings for three knot types are shown in Figure 4. The size of the rings is changed with the topology thus we rescale the distance r by Dg, which is more comparable to discuss about the threading structures. For flexible rings (Kb = 5) shown in Figure 4a, the peak appears at a far distance range of r = 0.5~0.75 D g , which means a weak penetration depth owing to the shrinkage of flexible rings. For semi-flexible rings (Kb = 25) and stiff rings (Kb = 80) shown in Figure 4b,c, the peaks shift to a short distance, especially for the 01-knot rings (r* ≈ 0.20 for Kb = 25). The results also confirm the fact the chain stiffness can help rings swell and free up more space for other rings, favorable for inter-rings penetration. For different knot types, the more complex the topological constraints, the farther the peak position and the lower the penetration depth, which is quite consistence with the results of Figure 3.
For a better view, the typical penetration structures of 01-knot rings with different bending energies Kb are investigated in Figure 5. The blue ring presents the active ring B and the red is the passive one A, while the white ring is not penetrated by the blue ring. Flexible rings barely penetrate each other due to the crumpled and filled structures shown in Figure 5a. For semi-flexible rings in Figure 5b, it is the deepest penetration between rings, which are involved in several rings. The further increase of bending energy Kb leads to an improve of the flatness of the polymer contours in Figure 5c. The planar rigid rings tend to be arranged parallel to each other, thus weakening the penetration depth with respect to their semi-flexible counterparts, showing a strong anisotropic character, which is strengthened by the increase of the stiffness.
As discussed above, the effects of increasing stiffness have two aspects. One is the expansion of the internal space of rings. The other is the anisotropic character which seems to weaken the penetration. We will focus on the anisotropic character in the following. Our finding shows that the competition of these two effects results in the non-monotonous dependence of Pth on stiffness. Only focus on the case of 01-knot rings, the probability distributions P(θ) of threading angle for different bending energies Kb are shown in Figure 6a. As we have mentioned in the KMT algorithm, one ring can be reduced into a triangle penetrated by another ring satisfying the minimal surface principle. Therefore, the threading angle θ is defined as the angle between a line and a triangular plane (ABC) based on the KMT algorithm, as shown in the Figure 6b. For the flexible rings, the black line exhibits an approximate Gaussian distribution with a mean angle of ~45°. It indicates the conformation of flexible rings is isotropic like a collapsed ball, as displayed in Figure 5a, leading to a full angle distribution. Flexible rings penetrate each other slightly due to the strong intramolecular barriers [42,43,44]. With the increase of bending energy Kb, not only the peak value of P(θ) increases but also the position of the peak shifts to a smaller angle. As well known, the shapes of rings change from prolate, crumpled structures to planar, rigid rings caused by the increase of chain stiffness, which can dramatically destroy the isotropic of the rings. Rigid rings (Kb = 80) tend to be parallel with each other as a result of the topology potential [42,43,44], leading to a quite small threading angle and even close to 0°. In a word, there is a typical transition from isotropic character to anisotropic character as chain stiffness increases.
Corresponding to the three knot types, the probability distributions P(θ) of threading angle for different bending energies are shown in Figure 7. For the flexible case, there is little difference among 01-knot, 31-knot and 61-knot rings, because all of them interact as an isotropic ball, corresponding to a full angle distribution. With the increase of bending energy Kb, the difference among the three knot types is highlighted, especially for rigid rings (Kb = 80) shown in Figure 7c. The peaks of P(θ) are suppressed and shifted to a broad angle as the knot complexity of the rings increases, which actually strengthens the isotropic of rings. As a result, the more rigid of the rings, the stronger of the topology constraints.
The probability distributions of P(θ) indicate that anisotropic plays an important role on the penetration structure. To explain the effect of anisotropic character on threading degree, we subdivide the threading angle θ (0°~90°) into three angle ranges: small angle (0°~30°), medium angle (30°~60°) and large angles (60°~90°). Corresponding to the black line in Figure 3, the threading level Pth of 01-knot rings vs bending energy Kb for three angle ranges are plotted separately in Figure 8. In general, the proportion of small angle threading (black line) increases monotonously along with Kb. It indicates that the increase of chain stiffness strengthens the anisotropic penetration. Therefore, the stronger the rigidity, the more penetration of the small angle will be. For medium angles (red line) and large angles (blue line), as Kb increases, Pth both increase first because of the expansion of the free space, then decrease dramatically owing to the anisotropic character. Especially for the case of large angles (60°~90°), the proportion of penetration dropped sharply or even close to 0, which fits well with the fact that rigid rings tend to be parallel with each other.
We look into the influence of topological constraints on different threading angles as well. As shown in Figure 9, the threading level of 01-knot, 31-knot and 61-knot rings present a similar dependence on Kb for three angle ranges. In detail, the more complex the topological constraints, the smaller the change in threading level Pth, owing to the decrease of the free space for penetration. What’s more, the complex topology also weakens the anisotropic character, resulting in the domain of the medium threading angles, which is also obvious in Figure 7 as well. Look closely at Figure 9a, there is a little decrease for 01-knot rings at a high bending energy (after Kb = 80). Rigid rings penetrate others nearly in parallel, where a ring can penetrate only a small amount of other rings, as visually presented in Figure 5c. This is responsible for a slight decline in threading level corresponding to the small threading angle.
The effective potential of mean force (PMF) between a pair of rings Veff for three knot types [45], obtained by umbrella sampling MD, is calculated in Figure 10. For all cases, Veff has decayed to zero as r ≈ Dg. From Figure 10a–c, it can be seen that the amplitude of the effective potential at short distances increases with the complexity of the knot, that is, the effective potential of 01-knot rings is lower than 31-knot or 61-knot rings. Combined with the results of Figure 4, take Kb = 25 for example, the threading peaks of 01-knot rings approximately occur at a distance of 0.2Dg, while 0.4Dg and 0.6Dg for 31-knot rings and for 61-knot rings, respectively. Corresponding to the PMF in Figure 10b, at r = 0.2Dg, the potential of 01-knot rings is much lower, while at r = 0.4Dg or r = 0.6Dg, 31-knot and 61-knot rings are lower than 01-knot rings. It also confirms that 01-knot rings penetrate with each other at a much shorter distance, while rings with complex knots occur at a longer distance. Correspondingly, the effective interaction preferentially occurs at a short distance providing more free space available for inter-ring penetration and causing a deep penetration. Therefore, the threading level of 01-knot rings is much higher than rings with complex knots, shown in Figure 3. Additionally, the PMF of only 01-knot rings for different bending energies Kb is presented in Figure 10d. The effective potential decreases first and then increases with the increase of bending energy, exactly explaining non-monotonous dependence of the threading level Pth on chain stiffness. Narros, et. al., also proposed that the penetration conformation for 01-knot rings is quite common for r < 0.25Dg, resulting in the cost in free energy which can be approximately estimated by a Flory theory [3].
Figure 11. The peak position shifts to the small bending energy a little when chain length increases from N = 64 to N = 256. For example, the peak is located at Kb = 25 for N = 256 and at Kb = 20 for N = 128 or 64. Meanwhile, the longer of the ring, the larger of the free space for penetration resulting in an increase of Pth for longer rings. Average threading number Pth decreases with knot topology and has a peak at Kb = 20~25 for a fixed length N = 256 is given in Figure S2 (see Figure S2 in Supporting Information). As for the density, we focus on the case of 01-knot rings and ρ = M*N/L3, ranging from 0.05 to 0.4 by changing the number of rings M. The threading level Pth vs bending energy Kb for different number densities ρ are shown in Figure 12. The non-monotonous behavior also becomes more and more prominent with increasing number density. Meanwhile, as ρ increases, the amplitude Pth increases greatly, since more rings have more opportunities to contact, which is available for inter-ring penetration. Then, the effects of number density ρ on the probability distribution P(θ) of threading angles is also presented in Figure 13. The number density has little effect on flexible rings due to the isotropic shrinkage shown in Figure 13a. The anisotropic character of rings is highlighted by the increase of chain stiffness. For semi-flexible rings (Kb = 25) shown in Figure 13b, the peak of P(θ) slightly shifts to a smaller angle with the increasing of the number density ρ, resulting from the fact that the anisotropic penetration is strengthened by the excluded volume effect for a large number density ρ. For rigid rings (Kb = 80) shown in Figure 13c, the chain stiffness absolutely dominates the anisotropic penetration, unaffected by number density. The corresponding typical structures of rigid rings are displayed in the inset. The planar rigid rings tend to be arranged parallel to each other with a quite small threading angle of ~15 ° .

4. Conclusions

The penetration structures of the semi-flexible ring polymers have been studied by MD simulations. With the implementation of the improved KMT algorithm, the degree of penetration and the internal structures have been characterized by the average threading number Pth and probability distribution P(r) of penetrated rings, respectively. Our results show that the penetration degree of rings does not increase monotonously with the chain stiffness. It depends on the balance between coil expansion (favoring penetrations) and stiffness (inhibiting penetrations). As the increase of chain stiffness, the expansion of ring conformations results in more space for penetration, while a subsequent decline in penetration results from the anisotropic character of semi-flexible rings. The anisotropic character of semi-flexible rings, which can be characterized by the probability distribution P(θ) of threading angle and corresponding average threading number Pth for different angle ranges, is further confirmed not favoring the penetration. The analysis of effective potential between a pair of rings provides a better understanding for this non-monotonous behavior in inter-ring penetration. Furthermore, it is shown that the penetration would be suppressed with the knot complexity of the rings. Our results reveal the essential roles of topological constraints and chain stiffness on inter-ring penetration.

Supplementary Materials

Author Contributions

Conceptualization, L.Z.; methodology, F.G.; resources, F.G. and K.L.; data curation, F.G. and J.W.; funding acquisition, L.H. and L.Z.; writing, original draft preparation, F.G. and L.H.; writing, review and editing, L.H. and L.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the National Natural Science Foundation of China (Grant Nos. 21674082, 21873082 and 21674096) and Zhejiang Province (No.LY19B040006).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Guo, H.; Mao, R.; Block, T.M.; Guo, J.T. Production and function of the cytoplasmic deproteinized relaxed circular DNA of hepadnaviruses. J. Virol. 2010, 84, 387–396. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Conway, J.H. An enumeration of knots and links and some of their related properties. In Computational Problems in Abstract Algebra; Leech, J., Ed.; Pergamon Press: Oxford, UK, 1967; pp. 329–358. [Google Scholar]
  3. Narros, A.; Moreno, A.J.; Likos, C.N. Influence of topology on effective potentials: Coarse-graining ring polymers. Soft Matter 2010, 6, 2435–2441. [Google Scholar] [CrossRef]
  4. Daichi, I. Semiflexible ring polymers in dilute solutions. React. Funct. Polym. 2018, 130, 111–117. [Google Scholar]
  5. Hiroki, A.; Daichi, I. Mean-square radius of gyration and scattering function of semiflexible ring Polymers of the trefoil knot. Polymers 2016, 8, 271. [Google Scholar]
  6. Moore, N.T.; Grosberg, A.Y. Limits of analogy between self-avoidance and topology-driven swelling of polymer loops. Phys. Rev. E 2005, 72, 061803. [Google Scholar] [CrossRef] [Green Version]
  7. Halverson, J.D.; Lee, W.B.; Grest, G.S.; Grosberg, A.Y.; Kremer, K. Molecular dynamics simulation study of nonconcatenated ring polymers in a melt. I. Statics. J. Chem. Phys. 2011, 134, 204904. [Google Scholar] [CrossRef] [Green Version]
  8. Reigh, S.Y.; Yoon, D.Y. Concentration dependence of ring polymer conformations from Monte Carlo simulations. ACS Macro. Lett. 2013, 2, 296–300. [Google Scholar] [CrossRef]
  9. Zhou, X.L.; Jiang, Y.W.; Deng, Z.Y.; Zhang, L.X. Glassy dynamics of nanoparticles in semiflexible ring polymer nanocomposite melts. Sci. Rep. 2017, 7, 44325. [Google Scholar] [CrossRef] [Green Version]
  10. Zhou, X.L.; Jiang, Y.W.; Chen, J.M.; He, L.L.; Zhang, L.X. Size-dependent nanoparticle dynamics in semiflexible ring polymer nanocomposites. Polymer 2017, 131, 243–251. [Google Scholar] [CrossRef]
  11. Doi, Y.; Matsubara, K.; Ohta, Y.; Nakano, T.; Kawaguchi, D.; Takahashi, Y.; Takano, A.; Matsushita, Y. Melt rheology of ring polystyrenes with ultrahigh purity. Macromolecules 2015, 48, 3140–3147. [Google Scholar] [CrossRef]
  12. Yan, Z.C.; Costanzo, S.; Jeong, Y.; Chang, T.; Vlassopoulos, D. Linear and nonlinear shear rheology of a marginally entangled ring polymer. Macromolecules 2016, 49, 1444–1453. [Google Scholar] [CrossRef]
  13. Doi, M.; Edwards, S.F. The Theory of Polymer Dynamics; Oxford University Press: New York, NY, USA, 1986. [Google Scholar]
  14. Kapnistos, M.; Lang, M.; Vlassopoulos, D.; Pyckhout-Hintzen, W.; Richter, D.; Cho, D.; Chang, T.; Rubinstein, M. Unexpected power-law stress relaxation of entangled ring polymers. Nat. Mater. 2008, 7, 997–1002. [Google Scholar] [CrossRef] [PubMed]
  15. Narros, A.; Moreno, A.J.; Likos, C.N. Effects of knots on ring polymers in solvents of varying quality. Macromolecules 2013, 46, 9437–9445. [Google Scholar] [CrossRef]
  16. Narros, A.; Moreno, A.J.; Likos, C.N. Effective Interactions of knotted ring polymers. Biochem. Soc. Trans. 2013, 41, 630–634. [Google Scholar] [CrossRef] [Green Version]
  17. Zhou, X.L.; Guo, F.C.; Li, K.; He, L.L.; Zhang, L.X. Entropy-induced separation of binary semiflexible ring polymer mixture in spherical confinement. Polymers 2019, 11, 1992. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Khokhlov, A.R.; Nechaev, S.K. Polymer chain in an array of obstacles. Phys. Lett. A 1985, 112, 156–160. [Google Scholar] [CrossRef]
  19. Obukhov, S.P.; Rubinstein, M.; Duke, T. Dynamics of a ring polymer in a gel. Phys. Rev. Lett. 1994, 73, 1263–1266. [Google Scholar] [CrossRef]
  20. Grosberg, A.Y. Annealed lattice animal model and Flory theory for the melt of non-concatenated rings: Towards the physics of crumpling. Soft Matter 2014, 10, 560–565. [Google Scholar] [CrossRef] [Green Version]
  21. Smrek, J.; Grosberg, A.Y. Understanding the dynamics of rings in the melt in terms of the annealed tree model. J. Phys. Condens. Matter 2015, 27, 064117. [Google Scholar] [CrossRef]
  22. Michieletto, D.; Marenduzzo, D.; Orlandini, E.; Alexander, G.P.; Turner, M.S. Threading dynamics of ring polymers in a gel. ACS Macro. Lett. 2014, 3, 255–259. [Google Scholar] [CrossRef]
  23. Michieletto, D.; Nahali, N.; Rosa, A. Glassiness and heterogeneous dynamics in dense solutions of ring polymers. Phys. Rev. Lett. 2017, 119, 197801. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Lee, E.; Kim, S.; Jung, Y. Slowing down of ring polymer diffusion caused by inter-ring threading. Macromol. Rapid Commun. 2015, 36, 1115–1121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Smrek, J.; Grosberg, A.Y. Minimal surfaces on unconcatenated polymer rings in melt. ACS Macro. Lett. 2016, 5, 750–754. [Google Scholar] [CrossRef]
  26. Halverson, J.D.; Smrek, J.; Kremer, K.; Grosberg, A.Y. From a melt of rings to chromosome territories: The role of topological constraints in genome folding. Rep. Prog. Phys. 2014, 77, 022601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Marenduzzo, D.; Micheletti, C.; Orlandini, E. Biopolymer organization upon confinement. J. Phys. Condens. Matter 2010, 22, 283102. [Google Scholar] [CrossRef]
  28. Kremer, K.; Grest, G.S. Dynamics of entangled linear polymer melts: A molecular-dynamics simulation. J. Chem. Phys. 1990, 92, 5057. [Google Scholar] [CrossRef]
  29. Grest, G.S.; Kremer, K. Molecular dynamics simulation for polymers in the presence of a heat bath. Phys. Rev. A 1986, 33, 3628. [Google Scholar] [CrossRef]
  30. Plimpton, S.J. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 1995, 117, 1–19. [Google Scholar] [CrossRef] [Green Version]
  31. Torrie, G.M.; Valleau, J.P. Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling. J. Comput. Phys. 1977, 23, 187–199. [Google Scholar] [CrossRef]
  32. Jiang, Y.W.; Zhang, D.; He, L.L.; Zhang, L.X. Entropic interactions in semiflexible polymer nanocomposite melts. J. Phys. Chem. B 2016, 120, 572–582. [Google Scholar] [CrossRef]
  33. Souaille, M.; Roux, B. Extension to the weighted histogram analysis method: Combining umbrella sampling with free energy calculations. Comput. Phys. Commun. 2001, 135, 40–57. [Google Scholar] [CrossRef]
  34. Deng, Z.Y.; Jiang, Y.W.; He, L.L.; Zhang, L.X. Aggregation–dispersion transition for nanoparticles in semiflexible ring polymer nanocomposite melts. J. Phys. Chem. B 2016, 120, 11574–11581. [Google Scholar] [CrossRef] [PubMed]
  35. Chandler, D. Introduction to Modern Statistical Mechanics; Oxford University Press: Oxford, UK, 1987. [Google Scholar]
  36. Frenkel, D.; Smit, B. Understanding Molecular Simulation; Academic Press: San Diego, CA, USA, 1996. [Google Scholar]
  37. Ferrenberg, A.M.; Swendsen, R.H. Optimized Monte Carlo data analysis. Phys. Rev. Lett. 1989, 63, 1195. [Google Scholar] [CrossRef] [PubMed]
  38. Koniaris, K.; Muthukumar, M. Knottedness in ring polymers. Phys. Rev. Lett. 1991, 66, 2211–2214. [Google Scholar] [CrossRef] [PubMed]
  39. Taylor, W.R. A deeply knotted protein structure and how it might fold. Nature 2000, 406, 916–919. [Google Scholar] [CrossRef]
  40. Virnau, P.; Kantor, Y.; Kardar, M. Knots in globule and coil phases of a model polyethylene. J. Am. Chem. Soc. 2005, 127, 15102–15106. [Google Scholar] [CrossRef]
  41. Gómez, L.R.; García, N.A.; Pöschel, T. Packing structure of semiflexible rings. Proc. Natl. Acad. Sci. USA 2020, 117, 3382–3387. [Google Scholar] [CrossRef] [Green Version]
  42. Bernabei, M.; Bacova, P.; Moreno, A.J.; Narros, A.; Likos, C.N. Fluids of semiflexible ring polymers: Effective potentials and clustering. Soft Matter 2013, 9, 1287–1300. [Google Scholar] [CrossRef] [Green Version]
  43. Poier, P.; Likos, C.N.; Moreno, A.J.; Blaak, R. An Anisotropic effective model for the simulation of semiflexible ring polymers. Macromolecules 2015, 9, 4983–4997. [Google Scholar] [CrossRef] [Green Version]
  44. Heermann, B. Topological interactions between ring polymers: Implications for chromatin loops. J. Chem. Phys. 2010, 132, 044904. [Google Scholar]
  45. Mladek, P.B.M.; Kahl, G.; Likos, C.N. Computer assembly of cluster-forming amphiphilic dendrimers. Phys. Rev. Lett. 2008, 100, 028301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Typical snapshots for two topological knots (a) 31-knot and (b) 61-knot.
Figure 1. Typical snapshots for two topological knots (a) 31-knot and (b) 61-knot.
Polymers 12 02659 g001
Figure 2. (a) Illustration of the KMT algorithm [38,39,40]. Ring A is passively threaded by ring B, which is omitted for simplicity with only one dot marking a cross bond. Monomers 1~7 are part of the ring A. After the first reduction, monomers 4 and 6 are removed as triangle 345 and 567 are not intersected by Ring B. (b) If the ring A is threaded by ring B, ring A can be reduced finally into a triangle after many KMT reductions.
Figure 2. (a) Illustration of the KMT algorithm [38,39,40]. Ring A is passively threaded by ring B, which is omitted for simplicity with only one dot marking a cross bond. Monomers 1~7 are part of the ring A. After the first reduction, monomers 4 and 6 are removed as triangle 345 and 567 are not intersected by Ring B. (b) If the ring A is threaded by ring B, ring A can be reduced finally into a triangle after many KMT reductions.
Polymers 12 02659 g002
Figure 3. Average threading number Pth shows a non-monotonous dependence on bending energy Kb for 01-knot, 31-knot and 61-knot ring polymers. The peak appears at the intermediate stiffness corresponding to the semi-flexible case. Inset: Average threading number Pth as a function of lp/L, where lp denotes the persistence length and L is the contour length of the chain.
Figure 3. Average threading number Pth shows a non-monotonous dependence on bending energy Kb for 01-knot, 31-knot and 61-knot ring polymers. The peak appears at the intermediate stiffness corresponding to the semi-flexible case. Inset: Average threading number Pth as a function of lp/L, where lp denotes the persistence length and L is the contour length of the chain.
Polymers 12 02659 g003
Figure 4. Probability density distributions P(r) of penetrated rings for different bending energies Kb corresponding to the case of (a) flexible rings, (b)semi-flexible rings, and (c) rigid rings with three knot types. P(r) denotes the probability distribution for the centers of mass between two penetrated rings. The r is rescaled by Dg thus the normalization condition satisfies   P ( r ) d r = Dg.
Figure 4. Probability density distributions P(r) of penetrated rings for different bending energies Kb corresponding to the case of (a) flexible rings, (b)semi-flexible rings, and (c) rigid rings with three knot types. P(r) denotes the probability distribution for the centers of mass between two penetrated rings. The r is rescaled by Dg thus the normalization condition satisfies   P ( r ) d r = Dg.
Polymers 12 02659 g004
Figure 5. Illustration of typical threading structures of 01-knot rings with different bending energies. (a) Kb = 5, flexible rings adopt coiled conformation. (b) Kb = 25, the semi-flexible rings pierce others forming a deep penetration. (c) Kb = 80, rigid rings get parallel with each other resulting in a slight threading.
Figure 5. Illustration of typical threading structures of 01-knot rings with different bending energies. (a) Kb = 5, flexible rings adopt coiled conformation. (b) Kb = 25, the semi-flexible rings pierce others forming a deep penetration. (c) Kb = 80, rigid rings get parallel with each other resulting in a slight threading.
Polymers 12 02659 g005
Figure 6. (a) Probability distributions P(θ) of threading angle for 01-knot rings with different bending energies Kb shows an obvious transition from isotropic character to anisotropic character. (b) The upper panel shows the typical snapshot for two penetrated rings and the lower panel is the simplified case to show the definition of the threading angle as a triangle intersected a line based on the KMT algorithm.
Figure 6. (a) Probability distributions P(θ) of threading angle for 01-knot rings with different bending energies Kb shows an obvious transition from isotropic character to anisotropic character. (b) The upper panel shows the typical snapshot for two penetrated rings and the lower panel is the simplified case to show the definition of the threading angle as a triangle intersected a line based on the KMT algorithm.
Polymers 12 02659 g006
Figure 7. Probability distributions P(θ) of threading angle for different bending energies (a) flexible rings, (b) semi-flexible rings, and (c) rigid rings with three knot types. The topology constraints suppress the anisotropic character, especially at the case of rigid rings.
Figure 7. Probability distributions P(θ) of threading angle for different bending energies (a) flexible rings, (b) semi-flexible rings, and (c) rigid rings with three knot types. The topology constraints suppress the anisotropic character, especially at the case of rigid rings.
Polymers 12 02659 g007
Figure 8. Average threading number Pth of 01-knot rings within different angel ranges as a function of bending energy Kb shows a strong anisotropic character. For rigid rings, the penetration of medium and big angles are suppressed due to the parallel structures.
Figure 8. Average threading number Pth of 01-knot rings within different angel ranges as a function of bending energy Kb shows a strong anisotropic character. For rigid rings, the penetration of medium and big angles are suppressed due to the parallel structures.
Polymers 12 02659 g008
Figure 9. Average threading number Pth within angles of 0°~30° (a), 30°~60° (b) and 60°~90° (c) as a function of bending energy Kb for three knot types. The anisotropic character and expansion effects are both suppressed by the complex topology.
Figure 9. Average threading number Pth within angles of 0°~30° (a), 30°~60° (b) and 60°~90° (c) as a function of bending energy Kb for three knot types. The anisotropic character and expansion effects are both suppressed by the complex topology.
Polymers 12 02659 g009
Figure 10. Potential of mean force (PMF) of three knot types (ac) and 01-knot rings shows the strong repulsion due to the topology constraints at short distances, resulting in the shift of the peak of P(r) in Figure 3 (d) PMF for different bending energies Kb present a similar non-monotonous dependence on Kb.
Figure 10. Potential of mean force (PMF) of three knot types (ac) and 01-knot rings shows the strong repulsion due to the topology constraints at short distances, resulting in the shift of the peak of P(r) in Figure 3 (d) PMF for different bending energies Kb present a similar non-monotonous dependence on Kb.
Polymers 12 02659 g010
Figure 11. Average threading number Pth increases with the chain length (N) and has a peak at Kb = 20~25 for 01-knot ring polymers at ρ = 0.1.
Figure 11. Average threading number Pth increases with the chain length (N) and has a peak at Kb = 20~25 for 01-knot ring polymers at ρ = 0.1.
Polymers 12 02659 g011
Figure 12. Average threading number Pth of 01-knot rings as a function of bending energy Kb for different number densities ρ. Pth increases with the increase of the number density.
Figure 12. Average threading number Pth of 01-knot rings as a function of bending energy Kb for different number densities ρ. Pth increases with the increase of the number density.
Polymers 12 02659 g012
Figure 13. Probability distributions P(θ) of threading angle for different bending energies Kb corresponding to the case of (a) flexible rings, (b)semi-flexible rings, and (c) rigid rings with different number densities ρ. The anisotropic character is strengthened by the increase of the density at the intermediate stiffness due to the excluded volume effect. The inset in (c) shows the rigid rings tend to be arranged parallel to each other with the increasing of number density ρ.
Figure 13. Probability distributions P(θ) of threading angle for different bending energies Kb corresponding to the case of (a) flexible rings, (b)semi-flexible rings, and (c) rigid rings with different number densities ρ. The anisotropic character is strengthened by the increase of the density at the intermediate stiffness due to the excluded volume effect. The inset in (c) shows the rigid rings tend to be arranged parallel to each other with the increasing of number density ρ.
Polymers 12 02659 g013aPolymers 12 02659 g013b
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Guo, F.; Li, K.; Wu, J.; He, L.; Zhang, L. Effects of Topological Constraints on Penetration Structures of Semi-Flexible Ring Polymers. Polymers 2020, 12, 2659. https://0-doi-org.brum.beds.ac.uk/10.3390/polym12112659

AMA Style

Guo F, Li K, Wu J, He L, Zhang L. Effects of Topological Constraints on Penetration Structures of Semi-Flexible Ring Polymers. Polymers. 2020; 12(11):2659. https://0-doi-org.brum.beds.ac.uk/10.3390/polym12112659

Chicago/Turabian Style

Guo, Fuchen, Ke Li, Jiaxin Wu, Linli He, and Linxi Zhang. 2020. "Effects of Topological Constraints on Penetration Structures of Semi-Flexible Ring Polymers" Polymers 12, no. 11: 2659. https://0-doi-org.brum.beds.ac.uk/10.3390/polym12112659

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop