Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userIndividual Polymer Chain Dynamics in an Entangled Polymeric Liquid Using a Stochastic Tube Model
74
0
0.0
(
0
)
Authors
Behrouz Behdani
Ask ChatGPT about the research
No Arabic abstract
This study focuses on comparing the individual polymer chain dynamics in an entangled polymeric liquid under different shear and extension rates. Polymer chains under various shear rates and extension rates were simulated using a stochastic-tube model [J. Rheol. 56: 1057 (2012)]. We developed a Matlab code to visualize and analyze the simulated configurations from the stochastic-tube model. We introduced new variables to determine how the extent of linearity changes with time for different shear rates, which is more useful than a typical end-to-end distance analysis. We identified whether the polymer chains undergo a tumbling rotation (slight elongation not accompanying contraction) or flipping rotation (elongation accompanying contraction). The simulation results indicate that the polymer chains exhibit a significant tendency to elongate at higher shear rates and occasionally experience flipping, while lower shear rates tend to exhibit very frequent tumbling. Furthermore, no rotations were observed under extensional flows. These results help clarifying uncertainty of previously proposed polymer deformation mechanisms of the convective constraint release and the configuration-dependent friction coefficient.
rate research
Read More
Many atomic liquids can form transient covalent bonds reminiscent of those in the corresponding solid states. These directional interactions dictate many important properties of the liquid state, necessitating a quantitative, atomic-scale understanding of bonding in these complex systems. A prototypical example is liquid silicon, wherein transient covalent bonds give rise to local tetrahedral order and consequent non-trivial effects on liquid state thermodynamics and dynamics. To further understand covalent bonding in liquid silicon, and similar liquids, we present an ab initio simulation-based approach for quantifying the structure and dynamics of covalent bonds in condensed phases. Through the examination of structural correlations among silicon nuclei and maximally localized Wannier function centers, we develop a geometric criterion for covalent bonds in liquid Si. We use this to monitor the dynamics of transient covalent bonding in the liquid state and estimate a covalent bond lifetime. We compare covalent bond dynamics to other processes in liquid Si and similar liquids and suggest experiments to measure the covalent bond lifetime.
Development of reliable interatomic potentials is crucial for theoretical studies of relationship between chemical composition, structure and observable properties in glass-forming metallic alloys. Due to ambiguity of potential parametrization procedure, some crucial properties of the system, such as crystallization stability or symmetry of the ground state crystal phase, may not be correctly reproduced in computer simulations. Here we address this issue for ${rm Cu_{64.5}Zr_{35.5}}$ alloy described by two modifications of embedded atom model potential as well as by textit{ab initio} molecular dynamics. We observe that, at low supercooling, both models provide very similar liquid structure, which agrees with that obtained by textit{ab initio} simulations. Hoverer, deeply supercooled liquids demonstrate essentially different local structure and so different crystallization stability. The system, which demonstrate more pronounced icosahedral sort-range order, is more stable to crystallization that is in agreement with Frank hypothesis.
Although common in nature, the self-assembly of small molecules at sold-liquid interfaces is difficult to control in artificial systems. The high mobility of dissolved small molecules limits their residence at the interface, typically restricting the self-assembly to systems under confinement or with mobile tethers between the molecules and the surface. Small hydrogen-bonding molecules can overcome these issues by exploiting group-effect stabilization to achieve non-tethered self-assembly at hydrophobic interfaces. Significantly, the weak molecular interactions with the solid makes it possible to influence the interfacial hydrogen bond network, potentially creating a wide variety of supramolecular structures. Here we investigate the nanoscale details of water and alcohols mixtures self-assembling at the interface with graphite through group effect. We explore the interplay between inter-molecular and surface interactions by adding small amounts of foreign molecules able to interfere with the hydrogen bond network and systematically varying the length of the alcohol hydrocarbon chain. The resulting supramolecular structures forming at room temperature are then examined using atomic force microscopy with insights from computer simulations. We show that the group-based self-assembly approach investigated here is general and can be reproduced on other substrates such as molybdenum disulphide and graphene oxide, potentially making it relevant for a wide variety of systems.
Two of the most successful methods that are presently available for simulating the quantum dynamics of condensed phase systems are centroid molecular dynamics (CMD) and ring polymer molecular dynamics (RPMD). Despite their conceptual differences, practical implementations of these methods differ in just two respects: the choice of the Parrinello-Rahman mass matrix and whether or not a thermostat is applied to the internal modes of the ring polymer during the dynamics. Here we explore a method which is halfway between the two approximations: we keep the path integral bead masses equal to the physical particle masses but attach a Langevin thermostat to the internal modes of the ring polymer during the dynamics. We justify this by showing analytically that the inclusion of an internal mode thermostat does not affect any of the desirable features of RPMD: thermostatted RPMD (TRPMD) is equally valid with respect to everything that has actually been proven about the method as RPMD itself. In particular, because of the choice of bead masses, the resulting method is still optimum in the short-time limit, and the transition state approximation to its reaction rate theory remains closely related to the semiclassical instanton approximation in the deep quantum tunneling regime. In effect, there is a continuous family of methods with these properties, parameterised by the strength of the Langevin friction. Here we explore numerically how the approximation to quantum dynamics depends on this friction, with a particular emphasis on vibrational spectroscopy. We find that a broad range of frictions approaching optimal damping give similar results, and that these results are immune to both the resonance problem of RPMD and the curvature problem of CMD.
Molecular dynamics (MD) simulation is a powerful computational tool to study the behavior of macromolecular systems. But many simulations of this field are limited in spatial or temporal scale by the available computational resource. In recent years, graphics processing unit (GPU) provides unprecedented computational power for scientific applications. Many MD algorithms suit with the multithread nature of GPU. In this paper, MD algorithms for macromolecular systems that run entirely on GPU are presented. Compared to the MD simulation with free software GROMACS on a single CPU core, our codes achieve about 10 times speed-up on a single GPU. For validation, we have performed MD simulations of polymer crystallization on GPU, and the results observed perfectly agree with computations on CPU. Therefore, our single GPU codes have already provided an inexpensive alternative for macromolecular simulations on traditional CPU clusters and they can also be used as a basis to develop parallel GPU programs to further speedup the computations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
