Diffusion of rigid rodlike polymer in isotropic solutions studied by dissipative particle dynamics simulation

Numerical Modeling of Anisotropic Particle Diffusion through a Cylindrical Channel

The transport of molecules and particles through single pores is the basis of biological processes, including DNA and protein sequencing. As individual objects pass through a pore, they cause a transient change in the current that can be correlated with the object size, surface charge, and even chemical properties. The majority of experiments and modeling have been performed with spherical objects, while much less is known about the transport characteristics of aspherical particles, which would act as a model system, for example, for proteins and bacteria. The transport kinetics of aspherical objects is an especially important, yet understudied, problem in nanopore analytics. Here, using the Wiener process, we present a simplified model of the diffusion of rod-shaped particles through a cylindrical pore, and apply it to understand the translation and rotation of the particles as they pass through the pore. Specifically, we analyze the influence of the particles' geometrical characteristics on the effective diffusion type, the first passage time distribution, and the particles' orientation in the pore. Our model shows that thicker particles pass through the channel slower than thinner ones, while their lengths do not affect the passage time. We also demonstrate that both spherical and rod-shaped particles undergo normal diffusion, and the first passage time distribution follows an exponential asymptotics. The model provides guidance on how the shape of the particle can be modified to achieve an optimal passage time.

Percolation analysis of the electrical conductive network in a polymer nanocomposite by nanorod functionalization

Chemical functionalization of nanofillers is an effective strategy to benefit the formation of the conductive network in the matrix which can enhance the electrical conductivity of polymer nanocomposites (PNCs). In this work, we adopted a coarse-grained molecular dynamics simulation to investigate the effect of the nanorod (NR) functionalization on the conductive probability of PNCs under the quiescent state or under a shear field. It is found that the direct aggregation structure of NRs is gradually broken down with increasing the NR functionalization degree λ_A, which improves their dispersion state. Moreover, a local bridging structure of NRs sandwiched via one polymer layer is formed at high λ_A. Corresponding to it, the percolation threshold of PNCs first quickly decreases, then increases and last slightly decreases again with the increase of λ_A, which exhibits an anti N-type under the quiescent state. Meanwhile, it shows a non-monotonic dependence on the interaction between polymer and the functionalized beads which reaches the lowest value at the moderate interaction. However, the percolation threshold is nearly independent of λ_A under the shear field. Compared with in the quiescent state, the decrease or the increase of the percolation threshold can be tuned by λ_A under the shear field. The significant change in the percolation threshold is attributed to the orientation and the dispersion state of NRs under the shear field, which affects the conductive network. Especially, we found that the dispersion state of NRs is different for different λ_A under the shear field. However, the percolation threshold is similar which indicates that the dispersion state of NRs is not completely correlated to the conductive network. In summary, this work presents some further understanding of how the NR functionalization affects the electrical conductivity of PNCs.

Chemical functionalization of nanofillers is an effective strategy to benefit the formation of the conductive network in the matrix which can enhance the electrical conductivity of polymer nanocomposites (PNCs).

Dynamically crowded solutions of infinitely thin Brownian needles

We study the dynamics of solutions of infinitely thin needles up to densities deep in the semidilute regime by Brownian dynamics simulations. For high densities, these solutions become strongly entangled and the motion of a needle is essentially restricted to a one-dimensional sliding in a confining tube composed of neighboring needles. From the density-dependent behavior of the orientational and translational diffusion, we extract the long-time transport coefficients and the geometry of the confining tube. The sliding motion within the tube becomes visible in the non-Gaussian parameter of the translational motion as an extended plateau at intermediate times and in the intermediate scattering function as an algebraic decay. This transient dynamic arrest is also corroborated by the local exponent of the mean-square displacements perpendicular to the needle axis. Moreover, the probability distribution of the displacements perpendicular to the needle becomes strongly non-Gaussian; rather, it displays an exponential distribution for large displacements. On the other hand, based on the analysis of higher-order correlations of the orientation we find that the rotational motion becomes diffusive again for strong confinement. At coarse-grained time and length scales, the spatiotemporal dynamics of the needle for the high entanglement is captured by a single freely diffusing phantom needle with long-time transport coefficients obtained from the needle in solution. The time-dependent dynamics of the phantom needle is also assessed analytically in terms of spheroidal wave functions. The dynamic behavior of the needle in solution is found to be identical to needle Lorentz systems, where a tracer needle explores a quenched disordered array of other needles.

Molecular dynamics simulation of the electrical conductive network formation of polymer nanocomposites with polymer-grafted nanorods

Grafting chains on the surface of a filler is an effective strategy to tune and control the filler conductive network, which can be utilized to fabricate polymer nanocomposites (PNCs) with high electrical conductivity.

Tube Concept for Entangled Stiff Fibers Predicts Their Dynamics in Space and Time

We study dynamically crowded solutions of stiff fibers deep in the semidilute regime, where the motion of a single constituent becomes increasingly confined to a narrow tube. The spatiotemporal dynamics for wave numbers resolving the motion in the confining tube becomes accessible in Brownian dynamics simulations upon employing a geometry-adapted neighbor list. We demonstrate that in such crowded environments the intermediate scattering function, characterizing the motion in space and time, can be predicted quantitatively by simulating a single freely diffusing phantom needle only, yet with very unusual diffusion coefficients.

Destruction and recovery of a nanorod conductive network in polymer nanocomposites via molecular dynamics simulation

By adopting coarse-grained molecular dynamics simulation, we investigate the effects of end-functionalization and shear flow on the destruction and recovery of a nanorod conductive network in a functionalized polymer matrix. We find that the end-functionalization of polymeric chains can enhance the electrical conductivity of nanorod filled polymer nanocomposites, indicated by the decrease of the percolation threshold. However, there exists an optimal end-functionalization extent to reach the maximum electrical conductivity. In the case of steady shear flow, both homogeneous conductive probability and directional conductive probability perpendicular to the shear direction decrease with the shear rate, while the directional conductive probability parallel to the shear direction increases. Importantly, we develop a semi-empirical equation to describe the change of the homogeneous conductive probability as a function of the shear rate. Meanwhile, we obtain an empirical formula describing the relationship between the anisotropy of the conductive probability and the orientation of the nanorods. In addition, the conductivity stability increases with increasing nanorod volume fraction. During the recovery process of the nanorod conductive network, it can be fitted well by the model combining classical percolation theory and a time-dependent nanorod aggregation kinetic equation. The fitted recovery rate is similar for different nanorod volume fractions. In summary, this work provides some rational rules for fabricating polymer nanocomposites with excellent performance of electrical conductivity.

Translational and rotational diffusion of a single nanorod in unentangled polymer melts

Polymer nanocomposites have been an issue of both academic and industrial interest due to promising electrical, mechanical, optical, and magnetic properties. The dynamics of nanoparticles in polymer nanocomposites is a key to understanding those properties of polymer nanocomposites and is important for applications such as self-healing nanocomposites. In this article we investigate the translational and the rotational dynamics of a single nanorod in unentangled polymer melts by employing extensive molecular dynamics simulations. A nanorod and polymers are modeled as semiflexible tangent chains of spherical beads. The stiffness of a nanorod is tuned by changing the bending potential between chemical bonds. When polymers are sufficiently long and the nanorod is stiff, the nanorod translates in an anisotropic fashion along the nanorod axis within time scales of translational relaxation times even in unentangled polymer melts. The rotational diffusion is suppressed more significantly than the translational diffusion as the polymer chain length is increased, thus the translational and rotational diffusion of the nanorod are decoupled. We also estimate the winding numbers of polymers, i.e., how many times a polymer winds the nanorod. The winding number increases with longer polymers but is relatively insensitive to the nanorod stiffness.

Induced polar order in sedimentation equilibrium of rod-like nanoswimmers

The active diffusion and sedimentation equilibrium of rod-like nanoswimmers with length L are investigated by dissipative particle dynamics. In the absence of propulsion, the nanorod has the rotational correlation time τθ and mobility μ0, which vary with L. On the basis of the mean squared displacement, the diffusive behavior of rod-like nanoswimmers subject to active force Fa is found to follow (D - D0) ∝ (μ0Fa)(2)τθ, where D0 depicts the Brownian diffusivity. When the nanoswimmer suspension is under the external force Fe, the balance between the downward migration and upward active diffusion yields the sedimentation length δ = D/(μ0Fe), which no longer obeys δ0 = kBT/Fe obtained from the Einstein-Smoluchowski relationship, D0 = μ0kBT. Different from the suspension of passive rods, the polar order is clearly seen for active rods. The local polar order is essentially constant within the distance of about 2δ from the bottom wall but decays as the distance is further increased. In this work, the active Peclet number is small compared to unity and the maximum polar order grows linearly with Fe/Fa. The polar order arises because it is easier for the nanoswimmer with the swimming direction opposite to the external force to escape from the bottom wall.