Source author record

Wen-Sheng Xu

Wen-Sheng Xu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
1topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2022arXiv

Thermodynamic-Dynamic Interrelations in Glass-Forming Polymer Fluids

We utilize the lattice cluster theory (LCT) of polymer fluids to show that the configurational entropy, enthalpy, and internal energy are all closely interrelated, as suggested by recent measurements by Caruthers and Medvedev, so that the generalized entropy theory (GET) of glass formation, a combination of the LCT and Adam-Gibbs model, can be recast in terms of any of these thermodynamic properties as a matter of convenience. Thermodynamic scaling, a form of density-temperature scaling exhibited by dynamic and some thermodynamic properties, is used to assess which thermodynamic properties are most naturally linked to dynamics, and we explore the origin of this scaling by both direct calculations based on the GET and molecular dynamics simulations of a coarse-grained polymer model. Through a combination of our comprehensive modeling of thermodynamic properties using the LCT and the highly predictive GET model for how the fluid thermodynamics relate to its dynamics, along with simulation results confirming these theoretical frameworks, we obtain insights into thermodynamic aspects of collective motion and the slow $β$-relaxation processes of glass-forming liquids.

preprint2021arXiv

Polymer Glass Formation: Role of Activation Free Energy, Configurational Entropy, and Collective Motion

We provide a perspective on polymer glass formation, with an emphasis on models in which the fluid entropy and collective particle motion dominate the theoretical description and data analysis. We first discuss the dynamics of liquids in the high temperature Arrhenius regime, where transition state theory is formally applicable. We then summarize the evolution of the entropy theory from a qualitative framework for organizing and interpreting temperature-dependent viscosity data by Kauzmann to the formulation of a hypothetical `ideal thermodynamic glass transition' by Gibbs and DiMarzio, followed by seminal measurements linking entropy and relaxation by Bestul and Chang and the Adam-Gibbs (AG) model of glass formation rationalizing the observations of Bestul and Chang. These developments laid the groundwork for the generalized entropy theory (GET), which merges an improved lattice model of polymer thermodynamics accounting for molecular structural details and enabling the analytic calculation of the configurational entropy with the AG model, giving rise to a highly predictive model of the segmental structural relaxation time of polymeric glass-forming liquids. The development of the GET has occurred in parallel with the string model of glass formation in which concrete realizations of the cooperatively rearranging regions are identified and quantified for a wide range of polymeric and other glass-forming materials. The string model has shown that many of the assumptions of AG are well supported by simulations, while others are certainly not, giving rise to an entropy theory of glass formation that is largely in accord with the GET. As the GET and string models continue to be refined, these models progressively grow into a more unified framework, and this Perspective reviews the present status of development of this promising approach to the dynamics of polymeric glass-forming liquids.

preprint2016arXiv

Generalized entropy theory of glass-formation in fully flexible polymer melts

The generalized entropy theory (GET) offers many insights into how molecular parameters influence polymer glass-formation. Given the fact that chain rigidity often plays a critical role in understanding the glass-formation of polymer materials, the GET was originally developed based on models of semiflexible chains. Consequently, all previous calculations within the GET considered polymers with some degree of chain rigidity. Motivated by unexpected results from computer simulations of fully flexible polymer melts concerning the dependence of thermodynamic and dynamic properties on the cohesive interaction strength ($ε$), the present paper employs the GET to explore the influence of $ε$ on glass-formation in models of polymer melts with a vanishing bending rigidity, i.e., fully flexible polymer melts. In accord with simulations, the GET for fully flexible polymer melts predicts that basic dimensionless thermodynamic properties (such as the thermal expansion coefficient and isothermal compressibility) are universal functions of the temperature scaled by $ε$ in the regime of low pressures. Similar scaling behavior is also found for the configurational entropy density in the GET for fully flexible polymer melts. Moreover, we find that the characteristic temperatures of glass-formation increase linearly with $ε$ and that the fragility is independent of $ε$ in fully flexible polymer melts, predictions that are again consistent with simulations of glass-forming polymer melts composed of fully flexible chains. Beyond an explanation of these general trends observed in simulations, the GET for fully flexible polymer melts predicts the presence of a positive residual configurational entropy at low temperatures, indicating a return to Arrhenius relaxation in the low temperature glassy state.

preprint2016arXiv

Polymer Glass-Formation in Variable Dimension

We explore the nature of glass-formation in variable spatial dimensionality ($d$) based on the generalized entropy theory, a synthesis of the Adam-Gibbs model with direct computation of the configurational entropy of polymer fluids using an established statistical mechanical model. We find that structural relaxation in the fluid state asymptotically becomes Arrhenius in the $d\rightarrow\infty$ limit and that the fluid transforms upon sufficient cooling above a critical dimension near $d=8$ into a dense amorphous state with a finite positive residual configurational entropy. Direct computations of the isothermal compressibility and thermal expansion coefficient, taken to be physical measures of packing frustration, demonstrate that these fluid properties strongly correlate with the fragility of glass-formation.

preprint2016arXiv

Self-assembly and glass-formation in a lattice model of telechelic polymer melts: Influence of stiffness of the sticky bonds

The lattice cluster theory (LCT) for strongly interacting, self-assembling telechelic polymers provides a theoretical tool that enables establishing the connections between important microscopic molecular details of self-assembling polymers and their bulk thermodynamics. The original LCT for self-assembly of telechelic polymers considers a model of fully flexible linear chains [J. Dudowicz and K. F. Freed, J. Chem. Phys. \textbf{136}, 064902 (2012)], while our recent work introduces a significant improvement to the LCT by including a description of chain semiflexibility for the bonds within each individual telecheic chain [W.-S. Xu and K. F. Freed, J. Chem. Phys. \textbf{143}, 024901 (2015)], but the physically associative (or called "sticky") bonds between the ends of the telechelics are left as fully flexible. Motivated by the ubiquitous presence of steric constraints on the association of real telechelic polymers that impart an additional degree of bond stiffness (or rigidity), the present paper further extends the LCT to permit the sticky bonds to be semiflexible but to have a stiffness differing from that within each telechelic chain. An analytical expression for the Helmholtz free energy is provided for this model of linear telechelic polymer melts, and illustrative calculations demonstrate the significant influence of the stiffness of the sticky bonds on the self-assembly and thermodynamics of telechelic polymers. A brief discussion is also provided for the impact of self-assembly on glass-formation by combining the LCT description for this extended model of telechelic polymers with the Adam-Gibbs relation between the structural relaxation time and the configurational entropy.

preprint2015arXiv

Glass formation in a mixture of hard disks and hard ellipses

We present an event-driven molecular dynamics study of glass formation in two-dimensional binary mixtures composed of hard disks and hard ellipses, where both types of particles have the same area. We demonstrate that characteristic glass-formation behavior appears upon compression under appropriate conditions in such systems. In particular, while a rotational glass transition occurs only for the ellipses, both types of particles undergo a kinetic arrest in the translational degrees of freedom at a single density. The translational dynamics for the ellipses is found to be faster than that for the disks within the same system, indicating that shape anisotropy promotes the translational motion of particles. We further examine the influence of mixture's composition and aspect ratio on the glass formation. For the mixtures with an ellipse aspect ratio of $k=2$, both translational and rotational glass transition densities decrease with increasing the disk concentration at a similar rate and hence, the two glass transitions remain close to each other at all concentrations investigated. By elevating $k$, however, the rotational glass transition density diminishes at a faster rate than the translational one, leading to the formation of an orientational glass for the ellipses between the two transitions. Our simulations imply that mixtures of particles with different shapes emerge as a promising model for probing the role of particle shape in determining the properties of glass-forming liquids. Furthermore, our work illustrates the potential of using knowledge concerning the dependence of glass-formation properties on mixture's composition and particle shape to assist in the rational design of amorphous materials.

preprint2015arXiv

Lattice model of linear telechelic polymer melts. I. Inclusion of chain semiflexibility in the lattice cluster theory

The lattice cluster theory (LCT) for the thermodynamics of polymer systems has recently been reformulated to treat strongly interacting self-assembling polymers composed of fully flexible linear telechelic chains [J. Dudowicz and K. F. Freed, J. Chem. Phys. \textbf{136}, 064902 (2012)]. Here, we further extend the LCT for linear telechelic polymer melts to include a description of chain semiflexibility, which is treated by introducing a bending energy penalty whenever a pair of consecutive bonds from a single chain lies along orthogonal directions. An analytical expression for the Helmholtz free energy is derived for the model of semiflexible linear telechelic polymer melts. The extension provides a theoretical tool for investigating the influence of chain stiffness on the thermodynamics of self-assembling telechelic polymers, and for further exploring the influence of self-assembly on glass formation in such systems.

preprint2015arXiv

Lattice model of linear telechelic polymer melts. II. Influence of chain stiffness on basic thermodynamic properties

The lattice cluster theory (LCT) for semiflexible linear telechelic melts, developed in paper I, is applied to examine the influence of chain stiffness on the average degree of self-assembly and the basic thermodynamic properties of linear telechelic polymer melts. Our calculations imply that chain stiffness promotes self-assembly of linear telechelic polymer melts that assemble on cooling when either polymer volume fraction $ϕ$ or temperature $T$ is high, but opposes self-assembly when both $ϕ$ and $T$ are sufficiently low. This allows us to identify a boundary line in the $ϕ$-$T$ plane that separates two regions of qualitatively different influence of chain stiffness on self-assembly. The enthalpy and entropy of self-assembly are usually treated as adjustable parameters in classical Flory-Huggins type theories for the equilibrium self-assembly of polymers, but they are demonstrated here to strongly depend on chain stiffness. Moreover, illustrative calculations for the dependence of the entropy density of linear telechelic polymer melts on chain stiffness demonstrate the importance of including semiflexibility within the LCT when exploring the nature of glass formation in models of linear telechelic polymer melts.

preprint2014arXiv

Influence of Cohesive Energy and Chain Stiffness on Polymer Glass Formation

The generalized entropy theory is applied to assess the joint influence of the microscopic cohesive energy and chain stiffness on glass formation in polymer melts using a minimal model containing a single bending energy and a single (monomer averaged) nearest neighbor van der Waals energy. The analysis focuses on the combined impact of the microscopic cohesive energy and chain stiffness on the magnitudes of the isobaric fragility parameter $m_P$ and the glass transition temperature $T_g$. The computations imply that polymers with rigid structures and weak nearest neighbor interactions are the most fragile, while $T_g$ becomes larger when the chains are stiffer and/or nearest neighbor interactions are stronger. Two simple fitting formulas summarize the computations describing the dependence of $m_P$ and $T_g$ on the microscopic cohesive and bending energies. The consideration of the combined influence of the microscopic cohesive and bending energies leads to the identification of some important design concepts, such as iso-fragility and iso-$T_g$ lines, where, for instance, iso-fragility lines are contours with constant $m_P$ but variable $T_g$. Several thermodynamic properties are found to remain invariant along the iso-fragility lines, while no special characteristics are detected along the iso-$T_g$ lines. Our analysis supports the widely held view that fragility provides more fundamental insight for the description of glass formation than $T_g$.

preprint2014arXiv

Lattice cluster theory for polymer melts with specific interactions

Despite the long-recognized fact that chemical structure and specific interactions greatly influence the thermodynamic properties of polymer systems, a predictive molecular theory that enables systematically addressing the role of chemical structure and specific interactions has been slow to develop even for polymer melts. While the lattice cluster theory (LCT) provides a powerful vehicle for understanding the influence of various molecular factors, such as monomer structure, on the thermodynamic properties of polymer melts and blends, the application of the LCT has heretofore been limited to the use of the simplest polymer model in which all united atom groups within the monomers of a species interact with a common monomer averaged van der Waals energy. Thus, the description of a compressible polymer melt involves a single van der Waals energy. As a first step towards developing more realistic descriptions to aid in the analysis of experimental data and the design of new materials, the LCT is extended here to treat models of polymer melts in which the backbone and side groups have different interaction strengths, so three energy parameters are present, namely, backbone-backbone, side group-side group, and backbone-side group interaction energies. Because of the great algebraic complexity of this extension, we retain maximal simplicity within this class of models by further specializing this initial study to models of polymer melts comprising chains with poly($n$-$α$-olefin) structures where only the end segments on the side chains may have different, specific van der Waals interaction energies with the other united atom groups. An analytical expression for the LCT Helmholtz free energy is derived for the new model. Illustrative calculations are presented to demonstrate the degree to which the thermodynamic properties of polymer melts can be controlled by specific interactions.

preprint2014arXiv

Relaxation dynamics in a binary hard-ellipse liquid

Structural relaxation in binary hard spherical particles has been shown recently to exhibit a wealth of remarkable features when size disparity or mixture's composition is varied. In this paper, we test whether or not similar dynamical phenomena occur in glassy systems composed of binary hard ellipses. We demonstrate via event-driven molecular dynamics simulation that a binary hard-ellipse mixture with an aspect ratio of two and moderate size disparity displays characteristic glassy dynamics upon increasing density in both the translational and the rotational degrees of freedom. The rotational glass transition density is found to be close to the translational one for the binary mixtures investigated. More importantly, we assess the influence of size disparity and mixture's composition on the relaxation dynamics. We find that an increase of size disparity leads, both translationally and rotationally, to a speed up of the long-time dynamics in the supercooled regime so that both the translational and the rotational glass transition shift to higher densities. By increasing the number concentration of the small particles, the time evolution of both translational and rotational relaxation dynamics at high densities displays two qualitatively different scenarios, i.e., both the initial and the final part of the structural relaxation slow down for small size disparity, while the short-time dynamics still slows down but the final decay speeds up in the binary mixture with large size disparity. These findings are reminiscent of those observed in binary hard spherical particles. Therefore, our results suggest a universal mechanism for the influence of size disparity and mixture's composition on the structural relaxation in both isotropic and anisotropic particle systems.

preprint2013arXiv

Hard ellipses: Equation of state, structure and self-diffusion

We present an event-driven molecular dynamics study for hard ellipses and assess the effects of aspect ratio and area fraction on their physical properties. For state points in the plane of aspect ratio (k=1-9) and area fraction (phi=0.01-0.8), we identify three different phases, including isotropic, plastic and nematic states. The equation of state (EOS) is shown for a wide range of aspect ratios and is compared with the scaled particle theory (SPT) for the isotropic states. We find that SPT provides a good description of the EOS for the isotropic phase of hard ellipses. At large fixed phi, the reduced pressure p increases with k in both the isotropic and the plastic phases, and interestingly, its dependence on k is rather weak in the nematic phase. We rationalize the thermodynamics of hard ellipses in terms of particle motions. The plastic crystal is shown to form for aspect ratios up to k=1.4, while appearance of the stable nematic phase starts approximately at k=3. We quantitatively determine the locations of the isotropic-plastic (I-P) transition and the isotropic-nematic (I-N) transition by analyzing the bond-orientation correlations and the angular correlations, respectively. As expected, the I-P transition point is found to increase with k, while a larger k leads to a smaller area fraction where the I-N transition takes place. Moreover, our simulations strongly support that the two-dimensional nematic phase in hard ellipses has only quasi-long-range orientational order. The self-diffusion of hard ellipses is further explored and connections are revealed between the structure and the self-diffusion. We discuss the relevance of our results to the glass transition in hard ellipses. Finally, the results of the isodiffusivity lines are evaluated for hard ellipses and we discuss the effect of spatial dimension on the diffusive dynamics of hard ellipsoidal particles.

preprint2013arXiv

Thermodynamic scaling of dynamics in polymer melts: Predictions from the generalized entropy theory

Many glass-forming fluids exhibit a remarkable thermodynamic scaling in which dynamic properties, such as the viscosity, the relaxation time, and the diffusion constant, can be described under different thermodynamic conditions in terms of a unique scaling function of the ratio rho^gamma/T, where rho is the density, T is the temperature, and gamma is a material dependent constant. Given the successes of the generalized entropy theory in elucidating the influence of molecular details on the universal properties of glass-forming polymers, this theory is extended here to investigate the thermodynamic scaling in polymer melts. The predictions of theory are in accord with the appearance of thermodynamic scaling for pressures not in excess of about 50 MPa. (The failure at higher pressures arises due to inherent limitations of a lattice model.) In line with arguments relating the magnitude of gamma to the steepness of the repulsive part of the intermolecular potential, the abrupt, square-well nature of the lattice model interactions lead, as expected, to much larger values of the scaling exponent. Nevertheless, the theory is employed to study how individual molecular parameters affect the scaling exponent in order to extract a molecular understanding of the information content contained in the exponent. The chain rigidity, cohesive energy, chain length, and the side group length are all found to significantly affect the magnitude of the scaling exponent, and the computed trends agree well with available experiments. The variations of gamma with these molecular parameters are explained by establishing a correlation between the computed molecular dependence of the scaling exponent and the fragility. Thus, the efficiency of packing the polymers is established as the universal physical mechanism determining both the fragility and the scaling exponent gamma.

preprint2012arXiv

Dynamics and correlation length scales of a glass-forming liquid in quiescent and sheared conditions

We numerically study dynamics and correlation length scales of a colloidal liquid in both quiescent and sheared conditions to further understand the origin of slow dynamics and dynamic heterogeneity in glass-forming systems. The simulation is performed in a weakly frustrated two-dimensional liquid, where locally preferred order is allowed to develop with increasing density. The four-point density correlations and bond-orientation correlations, which have been frequently used to capture dynamic and static length scales $ξ$ in a quiescent condition, can be readily extended to a system under steady shear in this case. In the absence of shear, we confirmed the previous findings that the dynamic slowing down accompanies the development of dynamic heterogeneity. The dynamic and static length scales increase with $α$-relaxation time $τ_α$ as power-law $ξ\simτ_α^μ$ with $μ>0$. In the presence of shear, both viscosity and $τ_α$ have power-law dependence on shear rate in the marked shear thinning regime. However, dependence of correlation lengths cannot be described by power laws in the same regime. Furthermore, the relation $ξ\simτ_α^μ$ between length scales and dynamics holds for not too strong shear where thermal fluctuations and external forces are both important in determining the properties of dense liquids. Thus, our results demonstrate a link between slow dynamics and structure in glass-forming liquids even under nonequilibrium conditions.

preprint2012arXiv

Effect of attractions on correlation length scales in a glass-forming liquid

There is growing evidence that slow dynamics and dynamic heterogeneity possess structural signatures in glass-forming liquids. However, even in the weakly frustrated glass-forming liquids, whether or not the dynamic heterogeneity has a structural origin is a matter of debate. Via molecular dynamics simulation, we present a study of examining the connection between dynamic heterogeneity and bond orientational order in a weakly frustrated glass-forming liquid in two dimensions by taking advantage of assessing the effect of attractions on the correlation length scales. We find that attractions can strongly affect relaxation dynamics, dynamic heterogeneity and the associated dynamic correlation length of the liquid, but their influence on bond orientational order and the associated static correlation length shows a manner reminiscent of the effect of attractions on the thermodynamics of liquids. This implies that the growth of bond orientational order and static correlation length scale might be merely a manifestation of favoring the configurational entropy in weakly frustrated glass-forming liquids. Thus, our results lead strong evidence that bond orientational order cannot provide a complete description of dynamic heterogeneity even in weakly frustrated glass-forming systems.

preprint2012arXiv

Structure, compressibility factor and dynamics of highly size-asymmetric binary hard-disk liquids

By using event-driven molecular dynamics simulation, we investigate effects of varying the area fraction of the smaller component on structure, compressibility factor and dynamics of the highly size-asymmetric binary hard-disk liquids. We find that the static pair correlations of the large disks are only weakly perturbed by adding small disks. The higher-order static correlations of the large disks, by contrast, can be strongly affected. The compressibility factor of the system first decreases and then increases upon increasing the area fraction of the small disks and separating different contributions to it allows to rationalize this non-monotonic phenomenon. Furthermore, adding small disks can influence dynamics of the system in quantitative and qualitative ways. For the large disks, the structural relaxation time increases monotonically with increasing the area fraction of the small disks at low and moderate area fractions of the large disks. In particular, "reentrant" behavior appears at sufficiently high area fractions of the large disks, strongly resembling the reentrant glass transition in short-ranged attractive colloids and the inverted glass transition in binary hard spheres with large size disparity. By tuning the area fraction of the small disks, relaxation process for the small disks shows concave-to-convex crossover and logarithmic decay behavior, as found in other binary mixtures with large size disparity. Moreover, diffusion of both species is suppressed by adding small disks. Long-time diffusion for the small disks shows power-law-like behavior at sufficiently high area fractions of the small disks, which implies precursors of a glass transition for the large disks and a localization transition for the small disks. Therefore, our results demonstrate the generic dynamic features in highly size-asymmetric binary mixtures.

preprint2011arXiv

Diffusive redistribution of small spheres in crystallization of highly asymmetric binary hard-sphere mixtures

We report a molecular dynamics study of crystallization in highly asymmetric binary hard-sphere mixtures, in which the large spheres can form a crystal phase while the small ones remain disordered during the crystallization process of the large spheres. By taking advantage of assisting crystal nucleation with a patterned substrate, direct evidence is presented that there is a close link between the diffusive redistribution of the small spheres and the crystal formation of the large spheres. Although the addition of a second component with large size disparity will not alter the crystal structure formed by the large spheres, the density profile of the small spheres displays corresponding changes at different crystallization stages and closely relates to the crystal growth, suggesting possible effect of small spheres on the crystallization kinetics.