Methods for Creation and Linear Elastic Response Analysis of Packings of Semi-flexible Soft Polymer Chains

R.C. Dennis

2024-03-15

Abstract:From understanding the sand on the beach to the foam on your beer, soft sphere simulations have been crucial to the study of the amorphous world around us. However, many of the materials we interact with on a daily basis aren't comprised of individual grains, but complex molecules and chains of polymers. By extending the soft sphere model to a model of linked spheres, we can learn more about the materials we interact with on a daily basis. In this methods paper, I show how one can find and study the physical properties of packings of flexible chains, rigid molecules, and everything in between. In addition to describing the energy landscape of these materials, these methods describe how to shear stabilize polymer packings, classify them in the jamming hierarchy, describe their elastic properties, and much more. This simple modification to soft sphere simulations has the potential to yield new discoveries surrounding glasses.

Soft Condensed Matter

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to create and analyze the linear elastic response of semi - flexible soft polymer chains and molecular packings. Specifically, the paper proposes a method to study the physical properties of materials composed of complex molecules and polymer chains by modifying the existing soft - sphere simulation techniques. These materials are very common in daily life, but their interactions are more complex than those of individual particles. The method in the paper not only describes the energy landscapes of these materials but also provides information on how to shear - stabilize polymer packings, classify their positions in the jamming hierarchy, and describe their elastic properties, etc. Simple modifications of this method may lead to new discoveries regarding glassy materials.

### Main problems
1. **Creating packings of semi - flexible soft polymer chains and molecules**: The paper proposes a method to generate these packings, which are in force balance at local energy minima and can be achieved through constraint conditions.
2. **Analyzing the linear elastic response**: The paper details how to perform shear - stabilization treatment on these packings and calculate their physical properties such as elastic modulus.
3. **Studying physical properties**: The paper explores how to find the normal modes, rattling clusters, and classify the positions in the jamming hierarchy of the packings through these methods.

### Method overview
- **Generating polymer chains**: Constraining particles by defining links to form specific chain - like structures.
- **Orthogonal - normalized constraint basis**: Finding the constraint basis through QR decomposition or the parallelized Gram - Schmidt process to ensure that the projection of forces does not violate the constraint conditions.
- **Lagrangian method**: Using the Lagrangian method to find displacements that satisfy the constraint conditions.
- **Higher - order corrections**: Periodically re - fixing the constraints through the Newton - Raphson method to prevent the constraints from failing due to non - linear contributions.
- **Shear - stabilization**: Minimizing energy to consider both position and strain degrees of freedom simultaneously, ensuring that the packing is in the lowest energy state in all strain degrees of freedom.
- **Minimization algorithm**: Applying methods such as gradient descent and FIRE algorithm to minimize energy.
- **Calculating the constrained Hessian matrix and stiffness matrix**: Finding normal modes and elastic modulus by calculating second derivatives.
- **Strict jamming test**: Using a linear programming algorithm to test whether the packing remains stable under all possible combinations of strain and perturbation.

### Formula examples
- **Constraint function**:
\[
c_{\langle ij \rangle} = s\sum_{\alpha}(x_{ij}^{\alpha})^2 - l_{\langle ij \rangle} = 0
\]
where \(l_{\langle ij \rangle}\) is the constant length of the link \(\langle ij \rangle\).

- **Constraint Jacobian matrix**:
\[
\Omega_{k\alpha}^{\langle ij \rangle}=\hat{\tau}_{ij}^{\alpha}(\delta_{jk}-\delta_{ik})
\]

- **Constraint force vector**:
\[
\vec{F}_{\text{cons}}=(1 - \Omega_o\Omega_o^T)\vec{F}
\]

- **Constrained Hessian matrix**:
\[
H_E = BH_0B^T
\]

- **Stress - strain relationship**:
\[
\vec{\varepsilon}=S\vec{\sigma}
\]

Through these methods and formulas, the paper provides a comprehensive framework for generating and analyzing packings of semi - flexible soft polymer chains and molecules, thereby further understanding the physical properties of these complex materials.