1
|
Grigas AT, Fisher A, Shattuck MD, O'Hern CS. Connecting polymer collapse and the onset of jamming. Phys Rev E 2024; 109:034406. [PMID: 38632799 DOI: 10.1103/physreve.109.034406] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2023] [Accepted: 02/13/2024] [Indexed: 04/19/2024]
Abstract
Previous studies have shown that the interiors of proteins are densely packed, reaching packing fractions that are as large as those found for static packings of individual amino-acid-shaped particles. How can the interiors of proteins take on such high packing fractions given that amino acids are connected by peptide bonds and many amino acids are hydrophobic with attractive interactions? We investigate this question by comparing the structural and mechanical properties of collapsed attractive disk-shaped bead-spring polymers to those of three reference systems: static packings of repulsive disks, of attractive disks, and of repulsive disk-shaped bead-spring polymers. We show that the attractive systems quenched to temperatures below the glass transition T≪T_{g} and static packings of both repulsive disks and bead-spring polymers possess similar interior packing fractions. Previous studies have shown that static packings of repulsive disks are isostatic at jamming onset, i.e., the number of interparticle contacts N_{c} matches the number of degrees of freedom, which strongly influences their mechanical properties. We find that repulsive polymer packings are hypostatic at jamming onset (i.e., with fewer contacts than degrees of freedom) but are effectively isostatic when including stabilizing quartic modes, which give rise to quartic scaling of the potential energy with displacements along these modes. While attractive disk and polymer packings are often considered hyperstatic with excess contacts over the isostatic number, we identify a definition for interparticle contacts for which they can also be considered as effectively isostatic. As a result, we show that the mechanical properties (e.g., scaling of the potential energy with excess contact number and low-frequency contribution to the density of vibrational modes) of weakly attractive disk and polymer packings are similar to those of isostatic repulsive disk and polymer packings. Our results demonstrate that static packings generated via attractive collapse or compression of repulsive particles possess similar structural and mechanical properties.
Collapse
Affiliation(s)
- Alex T Grigas
- Graduate Program in Computational Biology and Bioinformatics, Yale University, New Haven, Connecticut 06520, USA
- Integrated Graduate Program in Physical and Engineering Biology, Yale University, New Haven, Connecticut 06520, USA
| | - Aliza Fisher
- Department of Mechanical Engineering and Materials Science, Yale University, New Haven, Connecticut 06520, USA
| | - Mark D Shattuck
- Benjamin Levich Institute and Physics Department, The City College of New York, New York, New York 10031, USA
| | - Corey S O'Hern
- Graduate Program in Computational Biology and Bioinformatics, Yale University, New Haven, Connecticut 06520, USA
- Integrated Graduate Program in Physical and Engineering Biology, Yale University, New Haven, Connecticut 06520, USA
- Department of Mechanical Engineering and Materials Science, Yale University, New Haven, Connecticut 06520, USA
- Department of Physics, Yale University, New Haven, Connecticut 06520, USA
- Department of Applied Physics, Yale University, New Haven, Connecticut 06520, USA
| |
Collapse
|
2
|
Ramos PM, Herranz M, Martínez-Fernández D, Foteinopoulou K, Laso M, Karayiannis NC. Crystallization of Flexible Chains of Tangent Hard Spheres under Full Confinement. J Phys Chem B 2022; 126:5931-5947. [PMID: 35904560 DOI: 10.1021/acs.jpcb.2c03424] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
We present results from extensive Monte Carlo simulations on the crystallization of athermal polymers under full confinement. Polymers are represented as freely jointed chains of tangent hard spheres of uniform size. Confinement is applied through the presence of flat, parallel, and impenetrable walls in all dimensions. We analyze crystallization as the summation of two contributions: one that occurs in the bulk volume of the system (bulk crystallization), and one on the wall surfaces (surface crystallization). Depending on volume fraction initially amorphous (disordered) hard-sphere chain packings transit to the stable crystal phase. The established ordered morphologies consist primarily of hexagonal close-packed (HCP) crystals in the bulk volume and of triangular (TRI) crystals on the surface. As in the case of athermal packings in the bulk (without confinement), a structural competition is observed between the 5-fold local symmetry and the formation of close-packed crystallites. Effectively, the full confinement inside a cube favors the growth of the HCP crystal, as the FCC one is quite incompatible with the imposed spatial constraints. Consequently, we observe the formation of noncompact ordered motifs which grow from the surface to the inner volume of the simulation cell. We further compare the 2D and 3D crystals formed by monomeric hard spheres under the same simulation conditions. Significant differences are observed at low densities that tend to diminish as concentration increases.
Collapse
Affiliation(s)
- Pablo Miguel Ramos
- Institute for Optoelectronic Systems and Microtechnology (ISOM) and Escuela Técnica Superior de Ingenieros Industriales (ETSII), Universidad Politécnica de Madrid (UPM), José Gutiérrez Abascal 2, 28006 Madrid, Spain
| | - Miguel Herranz
- Institute for Optoelectronic Systems and Microtechnology (ISOM) and Escuela Técnica Superior de Ingenieros Industriales (ETSII), Universidad Politécnica de Madrid (UPM), José Gutiérrez Abascal 2, 28006 Madrid, Spain
| | - Daniel Martínez-Fernández
- Institute for Optoelectronic Systems and Microtechnology (ISOM) and Escuela Técnica Superior de Ingenieros Industriales (ETSII), Universidad Politécnica de Madrid (UPM), José Gutiérrez Abascal 2, 28006 Madrid, Spain
| | - Katerina Foteinopoulou
- Institute for Optoelectronic Systems and Microtechnology (ISOM) and Escuela Técnica Superior de Ingenieros Industriales (ETSII), Universidad Politécnica de Madrid (UPM), José Gutiérrez Abascal 2, 28006 Madrid, Spain
| | - Manuel Laso
- Institute for Optoelectronic Systems and Microtechnology (ISOM) and Escuela Técnica Superior de Ingenieros Industriales (ETSII), Universidad Politécnica de Madrid (UPM), José Gutiérrez Abascal 2, 28006 Madrid, Spain
| | - Nikos Ch Karayiannis
- Institute for Optoelectronic Systems and Microtechnology (ISOM) and Escuela Técnica Superior de Ingenieros Industriales (ETSII), Universidad Politécnica de Madrid (UPM), José Gutiérrez Abascal 2, 28006 Madrid, Spain
| |
Collapse
|
3
|
Ramos PM, Herranz M, Foteinopoulou K, Karayiannis NC, Laso M. Entropy-Driven Heterogeneous Crystallization of Hard-Sphere Chains under Unidimensional Confinement. Polymers (Basel) 2021; 13:1352. [PMID: 33919100 PMCID: PMC8122411 DOI: 10.3390/polym13091352] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2021] [Revised: 04/15/2021] [Accepted: 04/19/2021] [Indexed: 11/16/2022] Open
Abstract
We investigate, through Monte Carlo simulations, the heterogeneous crystallization of linear chains of tangent hard spheres under confinement in one dimension. Confinement is realized through flat, impenetrable, and parallel walls. A wide range of systems is studied with respect to their average chain lengths (N = 12 to 100) and packing densities (ϕ = 0.50 to 0.61). The local structure is quantified through the Characteristic Crystallographic Element (CCE) norm descriptor. Here, we split the phenomenon into the bulk crystallization, far from the walls, and the projected surface crystallization in layers adjacent to the confining surfaces. Once a critical volume fraction is met, the chains show a phase transition, starting from regions near the hard walls. The established crystal morphologies consist of alternating hexagonal close-packed or face-centered cubic layers with a stacking direction perpendicular to the confining walls. Crystal layer perfection is observed with an increasing concentration. As in the case of the unconstrained phase transition of athermal polymers at high densities, crystal nucleation and growth compete with the formation of sites of a fivefold local symmetry. While surface crystallites show perfection with a predominantly triangular character, the morphologies of square crystals or of a mixed type are also formed. The simulation results show that the rate of perfection of the surface crystallization is not significantly faster than that of the bulk crystallization.
Collapse
Affiliation(s)
| | | | | | - Nikos Ch. Karayiannis
- Institute for Optoelectronic Systems and Microtechnology (ISOM) and Escuela Técnica Superior de Ingenieros Industriales (ETSII), Universidad Politécnica de Madrid (UPM), José Gutierrez Abascal 2, 28006 Madrid, Spain; (P.M.R.); (M.H.); (K.F.); (M.L.)
| | | |
Collapse
|
4
|
Kwon T, Sung BJ. Confinement effects on the mechanical heterogeneity of polymer fiber glasses. Phys Rev E 2020; 102:052501. [PMID: 33327119 DOI: 10.1103/physreve.102.052501] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2020] [Accepted: 10/13/2020] [Indexed: 11/07/2022]
Abstract
Both polymer fiber glasses and bulk polymer glasses exhibit nonlinear mechanical responses under uniaxial deformation. In polymer fibers, however, polymer chains are confined strongly and the surface area is relatively large compared to their volume. The confinement and the surface may lead to the spatially heterogeneous relaxation of chains in polymer fibers. In this work we perform molecular dynamics simulations and investigate the relation between the heterogeneous dynamics and the nonlinear mechanical responses at a molecular level. Our molecular simulations capture successfully not only the nonlinear mechanical response but also the dependence of mechanical properties on the strain rate of typical polymer glasses as in experiments. We find that the local elastic modulus and the nonaffine displacement are spatially heterogeneous in the pre-yield regime, which results in a lower elastic modulus for polymer fibers than bulk polymer glasses. In the post-yield regime, those mechanical properties become relatively homogeneous. Monomers with large nonaffine displacement are localized mainly at the interfacial region in the pre-yield regime while highly nonaffine monomers are distributed throughout the fibers in the post-yield regime. We show that the nonaffine displacement during deformation relates closely to the mechanical response of the polymer fibers. We also find that in the strain-hardening regime there is a significant difference in the energetic contribution to the stress between polymer fibers and bulk polymers, for which the modulus of the strain-hardening regime of the polymer fibers is smaller than that of bulk polymers.
Collapse
Affiliation(s)
- Taejin Kwon
- Department of Chemistry, Sogang University, Seoul 04107, Republic of Korea
| | - Bong June Sung
- Department of Chemistry, Sogang University, Seoul 04107, Republic of Korea
| |
Collapse
|