Ionic Polymer Nanocomposites Subjected to Uniaxial Extension: A Nonequilibrium Molecular Dynamics Study.
Ontology highlight
ABSTRACT:
SUBMITTER: Moghimikheirabadi A
PROVIDER: S-EPMC8621629 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:We investigate nanoparticle (NP) dispersion, polymer conformations, entanglements and dynamics in ionic nanocomposites. To this end, we study nanocomposite systems with various spherical NP loadings, three different molecular weights, two different Bjerrum lengths, and two types of charge-sequenced polymers by means of molecular dynamics simulations. NP dispersion can be achieved in either oligomeric or entangled polymeric matrices due to the presence of electrostatic interactions. We show that the overall conformations of ionic oligomer chains, as characterized by their radii of gyration, are affected by the presence and the amount of charged NPs, while the dimensions of charged entangled polymers remain unperturbed. Both the dynamical behavior of polymers and NPs, and the lifetime and amount of temporary crosslinks, are found to depend on the ratio between the Bjerrum length and characteristic distance between charged monomers. Polymer-polymer entanglements start to decrease beyond a certain NP loading. The dynamics of ionic NPs and polymers is very different compared with their non-ionic counterparts. Specifically, ionic NP dynamics is getting enhanced in entangled matrices and also accelerates with the increase of NP loading.
| S-EPMC7694256 | biostudies-literature
Project description:Since the proposal of self-healing materials, numerous researchers have focused on exploring their potential applications in flexible sensors, bionic robots, satellites, etc. However, there have been few studies on the relationship between the morphology of the dynamic crosslink network and the comprehensive properties of self-healing polymer nanocomposites (PNCs). In this study, we designed a series of modified nanoparticles with different sphericity (η) to establish a supramolecular network, which provide the self-healing ability to PNCs. We analyzed the relationship between the morphology of the supramolecular network and the mechanical performance and self-healing behavior. We observed that as η increased, the distribution of the supramolecular network became more uniform in most cases. Examination of the segment dynamics of polymer chains showed that the completeness of the supramolecular network significantly hindered the mobility of polymer matrix chains. The mechanical performance and self-healing behavior of the PNCs showed that the supramolecular network mainly contributed to the mechanical performance, while the self-healing efficiency was dominated by the variation of η. We observed that appropriate grafting density is the proper way to effectively enhance the mechanical and self-healing performance of PNCs. This study provides a unique guideline for designing and fabricating self-healing PNCs with modified Nanoparticles (NPs).
| S-EPMC10422474 | biostudies-literature
Project description:Mechanical signals regulate functions of mechanosensitive proteins by inducing structural changes that are determinant for force-dependent interactions. Talin is a focal adhesion protein that is known to extend under mechanical load, and it has been shown to unfold via intermediate states. Here, we compared different nonequilibrium molecular dynamics (MD) simulations to study unfolding of the talin rod. We combined boxed MD (BXD), steered MD, and umbrella sampling (US) techniques and provide free energy profiles for unfolding of talin rod subdomains. We conducted BXD, steered MD, and US simulations at different detail levels and demonstrate how these different techniques can be used to study protein unfolding under tension. Unfolding free energy profiles determined by BXD suggest that the intermediate states in talin rod subdomains are stabilized by force during unfolding, and US confirmed these results.
| S-EPMC7474207 | biostudies-literature
Project description:We present results for the steady state shear rheology of non-concatenated, unentangled and marginally entangled ring poly(ethylene oxide) (PEO) melts from detailed, atomistic nonequilibrium molecular dynamics (NEMD) simulations, and compare them to the behavior of the corresponding linear melts. The applied flow field spans a wide range of shear rates, from the linear (Newtonian) to the highly non-linear (described by a power law) regime. For all melts studied, rings are found to exhibit shear thinning but to a lesser degree compared to linear counterparts, mostly due to their reduced deformability and stronger resistance to alignment in the direction of flow. These features are attributed to the more compact structure of ring molecules compared to linear chains; the latter are capable of adopting wider and more open conformations even under shear due to the freedom provided by the free ends. Similar to linear melts, rings also exhibit a first and a second normal stress coefficient; the latter is negative. The ratio of the magnitude of the two coefficients remains practically constant with shear rate and is systematically higher than the corresponding one for linear melts. Emphasis was also given to the statistics of terminal (re-orientational) relaxation times which we computed by analyzing all chains in the simulated systems one by one; it was demonstrated that long time dynamics are strongly heterogeneous both for rings and (especially) linears. Repeating the analysis under flow conditions, and as expected, we found that the applied flow field significantly suppresses dynamic heterogeneity, especially for high shear rates well beyond the Newtonian plateau. Finally, a detailed geometrical analysis revealed that the average population of ring-ring threading events in the longest melt studied here (the PEO-5k ring) remains practically unaffected by the imposed flow rate even at strong shear rates, except for multi-threadings which disappear. To further analyze this peculiar and rather unexpected effect, we computed the corresponding survival times and penetration lengths, and found that the overwhelming majority of threadings under shear are extremely weak constraints, as they are characterized by very small penetration lengths, thus also by short survival times. They are expected therefore to play only a minor (if any) role on chain dynamics.
| S-EPMC6680584 | biostudies-literature
Project description:Electrophoresis of a mixture of NaCl and CaCl2 in a lysozyme crystal is investigated using nonequilibrium molecular dynamics (MD) simulations. Upon exposure to an electric field, the stability of lysozyme is found to decrease slightly. This finding is demonstrated by increases in the root mean-square deviations of the heavy atoms of lysozyme, in the solvent-accessible surface area of hydrophobic residues, and in the number of hydrogen bonds between lysozyme and water. The solvent-accessible surface area of hydrophilic residues changes marginally, and the number of hydrogen bonds between lysozyme molecules decreases. Water molecules tend to align preferentially parallel to the electric field, and the dipole moment along the pore axis increases linearly with increasing field strength. Two pronounced layered structures are observed for Na+ and Ca2+ in the vicinity of protein surface, but only one enriched layer is observed for Cl-. The number distributions of all ions are nearly independent of the electric field. The water coordination numbers of all ions are smaller in the crystal than in aqueous bulk solution; however, the reverse is found for the Cl- coordination numbers of cations. Both the water and the Cl- coordination numbers are insensitive to the electric field. Ion diffusivities in the crystal are approximately 2 orders of magnitude smaller than those in aqueous bulk solution. The drift velocities of ions increase proportionally to the electric field, particularly at high strengths, and depend on ionic charge and coordination with oppositely charged ions. Electrical current exhibits a linear relationship with the field strength. The zero-field electrical conductivity is estimated to be 0.56 S/m, which is very close to 0.61 S/m as predicted by the Nernst-Einstein equation.
| S-EPMC2567938 | biostudies-literature
Project description:Most computer simulations of molecular dynamics take place under equilibrium conditions-in a closed, isolated system, or perhaps one held at constant temperature or pressure. Sometimes, extra tensions, shears, or temperature gradients are introduced to those simulations to probe one type of nonequilibrium response to external forces. Catalysts and molecular motors, however, function based on the nonequilibrium dynamics induced by a chemical reaction's thermodynamic driving force. In this scenario, simulations require chemostats capable of preserving the chemical concentrations of the nonequilibrium steady state. We develop such a dynamic scheme and use it to observe cycles of a particle-based classical model of a catenane-like molecular motor. Molecular motors are frequently modeled with detailed-balance-breaking Markov models, and we explicitly construct such a picture by coarse graining the microscopic dynamics of our simulations in order to extract rates. This work identifies inter-particle interactions that tune those rates to create a functional motor, thereby yielding a computational playground to investigate the interplay between directional bias, current generation, and coupling strength in molecular information ratchets.
| S-EPMC9033874 | biostudies-literature
Project description:The size of real nanoparticles (NPs) is polydisperse which can influence the electrical property of polymer nanocomposites (PNCs). Here, we explored the percolated network of mixed NPs with different sizes (small NPs and big NPs) by adopting a molecular dynamics simulation. The simulated results reveal that the big NPs are adverse to building the percolated network compared to the small NPs. Thus, the percolation threshold becomes higher along with increasing the mixing ratio, which denotes the concentration ratio of big NPs to the total NPs. For a better understanding of it, the dispersion state and the number and the size of clusters are employed to analyze the percolated network, which can explain the percolation threshold well. Furthermore, by adopting the Sun's theory (Macromolecules, 2009, 42, 459-463), small and big NPs exhibit a weak antagonistic effect in the simulation if their total concentration is fixed. On the one hand, the number of small NPs is larger than that of big NPs at the same concentration. In addition, one big NP can connect to more others than one small NP. These two contrast effects are responsible for it. Interestingly, the shear flow leads to more contact aggregation structure of NPs which is beneficial to build the new percolated networks. Especially, the big NPs play a more important role in forming the percolated network than small NPs. Consequently, the percolation threshold is reduced at a higher shear rate. In total, our research work provides a further understanding of how the mixed NPs with different sizes form the percolated network in polymer matrix.
| S-EPMC8232767 | biostudies-literature
Project description:A computational method is developed to carry out explicit solvent simulations of complex molecular systems under conditions of constant pH. In constant-pH simulations, preidentified ionizable sites are allowed to spontaneously protonate and deprotonate as a function of time in response to the environment and the imposed pH. The method, based on a hybrid scheme originally proposed by H. A. Stern (J. Chem. Phys. 2007, 126, 164112), consists of carrying out short nonequilibrium molecular dynamics (neMD) switching trajectories to generate physically plausible configurations with changed protonation states that are subsequently accepted or rejected according to a Metropolis Monte Carlo (MC) criterion. To ensure microscopic detailed balance arising from such nonequilibrium switches, the atomic momenta are altered according to the symmetric two-ends momentum reversal prescription. To achieve higher efficiency, the original neMD-MC scheme is separated into two steps, reducing the need for generating a large number of unproductive and costly nonequilibrium trajectories. In the first step, the protonation state of a site is randomly attributed via a Metropolis MC process on the basis of an intrinsic pKa; an attempted nonequilibrium switch is generated only if this change in protonation state is accepted. This hybrid two-step inherent pKa neMD-MC simulation method is tested with single amino acids in solution (Asp, Glu, and His) and then applied to turkey ovomucoid third domain and hen egg-white lysozyme. Because of the simple linear increase in the computational cost relative to the number of titratable sites, the present method is naturally able to treat extremely large systems.
| S-EPMC4535364 | biostudies-literature
Project description:Polymeric nanopores show unique transport properties and have attracted a great deal of scientific interest as a test system to study ionic and molecular transport at the nanoscale. By means of all-atom molecular dynamics, we simulated the ion dynamics inside polymeric polyethylene terephthalate nanopores. For this purpose, we established a protocol to assemble atomic models of polymeric material into which we sculpted a nanopore model with the key features of experimental devices, namely a conical geometry and a negative surface charge density. Molecular dynamics simulations of ion currents through the pore show that the protonation state of the carboxyl group of exposed residues have a considerable effect on ion selectivity, by affecting ionic densities and electrostatic potentials inside the nanopores. The role of high concentrations of Ca2+ ions was investigated in detail.
| S-EPMC2907245 | biostudies-literature
Project description:We investigate experimentally and theoretically the sequence of phases that occurs when a self-assembled monolayer of gold nanoparticles supported on a fluid is compressed uniaxially in a Langmuir trough. Uniaxial compression of the monolayer results in the appearance of lines that have been shown to be regions of trilayer. These lines exhibit complex patterns that depend on the extent of compression. We show that these patterns can be understood in terms of an equilibrium statistical mechanical theory, originally developed in the context of commensurate-incommensurate transitions in krypton monolayers adsorbed on graphite, in which there is an energy cost to line deformations and to line intersections. Even though line intersections are energetically costly, they lower the free energy because they cause the entropy of the system to increase when the density of lines is low enough. Our analytic and Monte Carlo analyses of the model demonstrate that the model exhibits two-phase coexistence. Our experimental observations are qualitatively consistent with the predictions of the model.
| S-EPMC3549100 | biostudies-literature