h andan arkın 1,2 , w olfhard janke 1
DESCRIPTION
Thermodynamics of a p olymer c hain in an a ttractive spherical cage. H andan Arkın 1,2 , W olfhard Janke 1. 1 Institut f ü r Theoretische Physik, Universität Leipzig, Postfach 100 920, D-04009 Leipzig, Germany 2 Department of Physics Engineering, Faculty of Engineering, - PowerPoint PPT PresentationTRANSCRIPT
![Page 1: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/1.jpg)
1
Handan Arkın1,2, Wolfhard Janke1
1 Institut für Theoretische Physik, Universität Leipzig,Postfach 100 920, D-04009 Leipzig, Germany
2 Department of Physics Engineering, Faculty of Engineering,Ankara University Tandogan, 06100 Ankara, Turkey
![Page 2: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/2.jpg)
2
Outline
■ Motivation
■ Model: Off-lattice Polymer inside an Attractive Sphere
■ Multicanonical Monte Carlo Method
■ Some Simulation results
![Page 3: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/3.jpg)
3
Motivation
Investigating basis structure formation mechanisms of biomolecules at different interfaces is one of the major challenges of modern interdisciplinary research and possible application in nanotechnology
Applications: Polymer adhesion to metals, semiconductors biomedical implants and biosensors, smart drugs etc.
Due to the complexity introduced by the huge amount of possible substrate structures and sequence variations this problem is not trivial.
Therefore the theoretical treatment of the adsorption of macromolecules within the framework of minimalistic coarse-grained polymer models in statistical mechanics has been a longstanding problem.
![Page 4: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/4.jpg)
4
The polymer Model
2
1 2612
2
1
114cos1
4
1 N
i
N
ij ijij
N
kip
rrE
Bending energyLennard-Jones
potential
Interaction energy:
■ Coarse-grained off-lattice
■ Semi-flexible
■ Three-dimensional
Adjacent monomers are connected by rigid covalent bonds and the distance between them is fixed and set to unityThe position vector of the ith monomer is ir
![Page 5: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/5.jpg)
5
Polymer Chain inside an Attractice Sphere Potential
■ The interaction of polymer monomers and the attractive sphere is of van der Walls type, modeled by LJ 12-6
We integrate this potential over the entire sphere inner surface
Where the surface element in spherical coordinates is
![Page 6: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/6.jpg)
6
Polymer Chain inside an Attractice Sphere Potential
icicicicr
rrRrRrRrR
RVi
ccs
441010
25
14)(
cR The radius of the sphere
ir The distance of a monomer to the origin
c Set to unity
Varied during simulation
![Page 7: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/7.jpg)
7
Polymer Chain inside an Attractice Sphere Potential
2
1 2612
2
1
114cos1
4
1 N
i
N
ij ijij
N
kip
rrE
icicicic rRrRrRrRV
N
i iccs rRr
441010
125
114)(
20N
20cR4.11.0
![Page 8: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/8.jpg)
Simulations in Generalized Ensembles
• Idea: choose ensemble that allows better sampling
• Example: Multicanonical Ensemble • General Procedure: Determine the weight factors Large scale simulation Calculate expectation values for desired temperatures
• Problem: Find good estimators for the weights
• Advantages: Any energy barrier can be crossed. The probability of finding the global minimum is enhanced. Thermodynamic quantities for a range of temperatures
![Page 9: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/9.jpg)
9
Multicanonical Monte Carlo Method
• The canonical ensemble samples with the Boltzmann probability density
where x labels the configuration. The probability of the
energy E is
where n(E) is the density of states.
• The Muca ensemble is based on a probability function in which the different energies are equally probable:
where w(E) are multicanonical weight factors.
ZTkExP Bx
B //exp
ZTkEEnEP B
B
T //exp
.constEwEnEPMU
![Page 10: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/10.jpg)
10
Observables
Specific Heat:2
22
)(NT
TCEE
V
Radius of gyration:
Mean number of adsorbed monomers to the inner wall of the sphere:
We define a monomer i is being adsorbed if
and this can be expressed as .
N
i
N
j jicmi
N
ig NrrNrrR1 1
222
1
2 2/)(/)(
2.1 cic rrR
N
i ics rrN1
)(
![Page 11: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/11.jpg)
11
Observables
■ Gyration Tensor
Transformation to principal axis system diagonalizes S
Where the eigenvalues are sorted in descending order
![Page 12: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/12.jpg)
12
Observables
■ The first invariant of Gyration Tensor
The second invariant shape descriptor is Shape Anisotropy (reflects both symmetry and dimensionality)
Where
The last shape desciptor is the asphericity parameter
![Page 13: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/13.jpg)
13
Results: Pseudo Phase Diagram
D : Desorbed
A: Adsorbed
E: Extended
G: Globular
C:Compact
![Page 14: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/14.jpg)
14
Results: Specific -Heat
![Page 15: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/15.jpg)
15
Results: Radius of Gyration
![Page 16: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/16.jpg)
16
Results: The mean number of adsorbed monomers
![Page 17: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/17.jpg)
17
Results: Relative Shape Anisotropy
![Page 18: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/18.jpg)
18
Results: The eigenvalues of the Gyration tensor
![Page 19: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/19.jpg)
19
Results: Fluctuations of the Observables
![Page 20: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/20.jpg)
20
Results: Low Energy Conformations
1.0 4.0 7.0 0.1
![Page 21: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/21.jpg)
21
Acknowlegments
Thanks to ...
CQT Team, Institute of Theoretical Physics, Leipzig
![Page 22: H andan Arkın 1,2 , W olfhard Janke 1](https://reader035.vdocument.in/reader035/viewer/2022062309/56814918550346895db651e5/html5/thumbnails/22.jpg)
22
Results: The ratio of the greatest eigenvalue to the smallest