Spatial Organization of Slit-Confined Melts of Ring Polymers with Nonconserved Topology: A Lattice Monte Carlo Study
Ontology highlight
ABSTRACT:
SUBMITTER:
PROVIDER: S-EPMC10569094 | biostudies-literature | 2023 Sep
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Abnormally slower diffusional processes than its internal structure relaxation have been observed in ring polymeric melt systems recently. A key structural feature in ring polymer melts is topological constraints which allow rings to assume a threading configuration in the melt phase. In this work, we constructed a lattice model under the assumption of asymmetric diffusivity between two threading rings, and investigated a link between the structural correlation and its dynamic behavior via Monte Carlo simulations. We discovered that the hierarchical threading configurations render the whole system to exhibit abnormally slow dynamics. By analyzing statistical distributions of timescales of threading configurations, we found that the decoupling between internal structure relaxation and diffusion is crucial to understand the threading effects on the dynamics of a ring melt. In particular, in the limit of small but threaded rings, scaling exponents of the diffusion coefficient D and timescale τ diff with respect to the degree of polymerization N agree well with that of the annealed tree model as well as our mean-field analysis. As N increases, however, the ring diffusion abruptly slows down to the glassy behavior, which is supported by a breakdown of the Stokes⁻Einstein relation.
| S-EPMC6473489 | biostudies-literature
Project description:Ionic transport in solid electrolytes can often be approximated as ions performing a sequence of hops between distinct lattice sites. If these hops are uncorrelated, quantitative relationships can be derived that connect microscopic hopping rates to macroscopic transport coefficients; i.e. tracer diffusion coefficients and ionic conductivities. In real materials, hops are uncorrelated only in the dilute limit. At non-dilute concentrations, the relationships between hopping frequency, diffusion coefficient and ionic conductivity deviate from the random walk case, with this deviation quantified by single-particle and collective correlation factors, f and fI, respectively. These factors vary between materials, and depend on the concentration of mobile particles, the nature of the interactions, and the host lattice geometry. Here, we study these correlation effects for the garnet lattice using lattice-gas Monte Carlo simulations. We find that, for non-interacting particles (volume exclusion only), single-particle correlation effects are more significant than for any previously studied three-dimensional lattice. This is attributed to the presence of two-coordinate lattice sites, which causes correlation effects intermediate between typical three-dimensional and one-dimensional lattices. Including nearest-neighbour repulsion and on-site energies produces more complex single-particle correlations and introduces collective correlations. We predict particularly strong correlation effects at xLi=3 (from site energies) and xLi=6 (from nearest-neighbour repulsion), where xLi=9 corresponds to a fully occupied lithium sublattice. Both effects are consequences of ordering of the mobile particles. Using these simulation data, we consider tuning the mobile-ion stoichiometry to maximize the ionic conductivity, and show that the 'optimal' composition is highly sensitive to the precise nature and strength of the microscopic interactions. Finally, we discuss the practical implications of these results in the context of lithium garnets and other solid electrolytes.
| S-EPMC5717647 | biostudies-literature
Project description:The motion of nanoparticles (NPs) in entangled melts of linear polymers and nonconcatenated ring polymers are compared by large-scale molecular dynamics simulations. The comparison provides a paradigm for the effects of polymer architecture on the dynamical coupling between NPs and polymers in nanocomposites. Strongly suppressed motion of NPs with diameter d larger than the entanglement spacing a is observed in a melt of linear polymers before the onset of Fickian NP diffusion. This strong suppression of NP motion occurs progressively as d exceeds a and is related to the hopping diffusion of NPs in the entanglement network. In contrast to the NP motion in linear polymers, the motion of NPs with d > a in ring polymers is not as strongly suppressed prior to Fickian diffusion. The diffusion coefficient D decreases with increasing d much slower in entangled rings than in entangled linear chains. NP motion in entangled nonconcatenated ring polymers is understood through a scaling analysis of the coupling between NP motion and the self-similar entangled dynamics of ring polymers.
| S-EPMC5379250 | biostudies-literature
Project description:Connecting the viscoelastic behavior of stressed ring melts to the various forms of entanglement that can emerge in such systems is still an open challenge. Here, we consider active ring melts, where stress is generated internally, and introduce a topology-based method to detect and track consequential forms of ring entanglements, namely, deadlocks. We demonstrate that, as stress accumulates, more and more rings are co-opted in a growing web of deadlocks that entrap many other rings by threading, bringing the system to a standstill. The method ought to help the study of topological aging in more general polymer contexts.
| S-EPMC10883035 | biostudies-literature
Project description:By incorporating predicted secondary and tertiary restraints derived from multiple sequence alignments into ab initio folding simulations, it has been possible to assemble native-like tertiary structures for a test set of 19 nonhomologous proteins ranging from 29 to 100 residues in length and representing all secondary structural classes. Secondary structural restraints are provided by the PHD secondary structure prediction algorithm that incorporates multiple sequence information. Multiple sequence alignments also provide predicted tertiary restraints via a two-step process: First, seed side chain contacts are selected from a correlated mutation analysis, and then an inverse folding algorithm expands these seed contacts. The predicted secondary and tertiary restraints are incorporated into a lattice-based, reduced protein model for structure assembly and refinement. The resulting native-like topologies exhibit a coordinate root-mean-square deviation from native for the whole chain between 3.1 and 6.7 A, with values ranging from 2.6 to 4.1 A over approximately 80% of the structure. Overall, this study suggests that the use of restraints derived from multiple sequence alignments combined with a fold assembly algorithm is a promising approach to the prediction of the global topology of small proteins.
| S-EPMC18658 | biostudies-literature
Project description:A scaling model of self-similar conformations and dynamics of nonconcatenated entangled ring polymers is developed. Topological constraints force these ring polymers into compact conformations with fractal dimension df = 3 that we call fractal loopy globules (FLGs). This result is based on the conjecture that the overlap parameter of subsections of rings on all length scales is the same and equal to the Kavassalis-Noolandi number OKN ≈ 10-20. The dynamics of entangled rings is self-similar and proceeds as loops of increasing sizes are rearranged progressively at their respective diffusion times. The topological constraints associated with smaller rearranged loops affect the dynamics of larger loops through increasing the effective friction coefficient but have no influence on the entanglement tubes confining larger loops. As a result, the tube diameter defined as the average spacing between relevant topological constraints increases with time t, leading to "tube dilation". Analysis of the primitive paths in molecular dynamics simulations suggests a complete tube dilation with the tube diameter on the order of the time-dependent characteristic loop size. A characteristic loop at time t is defined as a ring section that has diffused a distance equal to its size during time t. We derive dynamic scaling exponents in terms of fractal dimensions of an entangled ring and the underlying primitive path and a parameter characterizing the extent of tube dilation. The results reproduce the predictions of different dynamic models of a single nonconcatenated entangled ring. We demonstrate that traditional generalization of single-ring models to multi-ring dynamics is not self-consistent and develop a FLG model with self-consistent multi-ring dynamics and complete tube dilation. This selfconsistent FLG model predicts that the longest relaxation time of nonconcatenated entangled ring polymers scales with their degree of polymerization N as τrelax ~ N7/3, while the diffusion coefficient of these rings scales as D3d ~ N-5/3. For the entangled solutions and melts of rings, we predict power law stress relaxation function G(t) ~ t-3/7 at t < τrelax without a rubbery plateau and the corresponding viscosity scaling with the degree of polymerization N as η ~ N4/3. These theoretical predictions are in good agreement with recent computer simulations and are consistent with experiments of melts of nonconcatenated entangled rings.
| S-EPMC4819263 | biostudies-literature
Project description:Droplet nucleation and evaporation are ubiquitous in nature and many technological applications, such as phase-change cooling and boiling heat transfer. So far, the description of these phenomena at the molecular scale has posed challenges for modelling with most of the models being implemented on a lattice. Here, we propose an off-lattice Monte-Carlo approach combined with a grid that can be used for the investigation of droplet formation and evaporation. We provide the details of the model, its implementation as Python code, and results illustrating its dependence on various parameters. The method can be easily extended for any force-field (e.g., coarse-grained, all-atom models, and external fields, such as gravity and electric field). Thus, we anticipate that the proposed model will offer opportunities for a wide range of studies in various research areas involving droplet formation and evaporation and will also form the basis for further method developments for the molecular modelling of such phenomena.
| S-EPMC8122685 | biostudies-literature
Project description:Due to the complexity of biological membrane, computer simulation of transmembrane protein's folding is challenging. In this paper, an implicit biological membrane environment has been constructed in lattice space, in which the lipid chains and water molecules were represented by the unoccupied lattice sites. The biological membrane was characterized with three features: stronger hydrogen bonding interaction, membrane lateral pressure, and lipophobicity index for the amino acid residues. In addition to the hydrocarbon core spanning region and the water solution, the lipid interface has also been represented in this implicit membrane environment, which was proved to be effective for the transmembrane protein's folding. The associated Monte Carlo simulations have been performed for SARS-CoV E protein and M2 protein segment (residues 18-60) of influenza A virus. It was found that the coil-helix transition of the transmembrane segment occurred earlier than the coil-globule transition of the two terminal domains. The folding process and final orientation of the amphipathic helical block in water solution are obviously influenced by its corresponding hydrophobicity/lipophobicity. Therefore, this implicit membrane environment, though in lattice space, can make an elaborate balance between different driving forces for the membrane protein's folding, thus offering a potential means for the simulation of transmembrane protein oligomers in feasible time.
| S-EPMC7117040 | biostudies-literature
Project description:We developed a hybrid Monte Carlo self-consistent field technique to model physical gels composed of ABA triblock copolymers and gain insight into the structure and interactions in such gels. The associative A blocks of the polymers are confined to small volumes called nodes, while the B block can move freely as long as it is connected to the A blocks. A Monte Carlo algorithm is used to sample the node configurations on a lattice, and Scheutjens-Fleer self-consistent field (SF-SCF) equations are used to determine the change in free energy. The advantage of this approach over more coarse grained methods is that we do not need to predefine an interaction potential between the nodes. Using this MC-SCF hybrid simulation, we determined the radial distribution functions of the nodes and structure factors and osmotic compressibilities of the gels. For a high number of polymers per node and a solvent-B Flory-Huggins interaction parameter of 0.5, phase separation is predicted. Because of limitations in the simulation volume, we did however not establish the full phase diagram. For comparison, we performed some coarse-grained MC simulations in which the nodes are modeled as single particles with pair potentials extracted from SF-SCF calculations. At intermediate concentrations, these simulations gave qualitatively similar results as the MC-SCF hybrid. However, at relatively low and high polymer volume fractions, the structure of the coarse-grained gels is significantly different because higher-order interactions between the nodes are not accounted for. Finally, we compare the predictions of the MC-SCF simulations with experimental and modeling data on telechelic polymer networks from literature.
| S-EPMC6328284 | biostudies-literature
Project description:The residual deposits usually left near the contact line after pinned sessile colloidal droplet evaporation are commonly known as a "coffee-ring" effect. However, there were scarce attempts to simulate the effect, and the realistic fully three-dimensional (3D) model is lacking since the complex drying process seems to limit the further investigation. Here we develop a stochastic method to model the particle deposition in evaporating a pinned sessile colloidal droplet. The 3D Monte Carlo model is developed in the spherical-cap-shaped droplet. In the algorithm, the analytical equations of fluid flow are used to calculate the probability distributions for the biased random walk, associated with the drift-diffusion equations. We obtain the 3D coffee-ring structures as the final results of the simulation and analyze the dependence of the ring profile on the particle volumetric concentration and sticking probability.
| S-EPMC3945784 | biostudies-literature