Next Article in Journal
Low-Temperature Deposition of Transparent Conducting Films Applied to Flexible Electrochromic Devices
Previous Article in Journal
A Theoretical Model for Debonding Prediction in the RC Beams Externally Strengthened with Steel Strip and Inorganic Matrix
Previous Article in Special Issue
Modeling of Nonlinear Optical Phenomena in Host-Guest Systems Using Bond Fluctuation Monte Carlo Model: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of Monte Carlo Dynamic/Kinetic Properties of Local Structure in Bond Fluctuation Model of Polymer System

Department of Theoretical Physics, Wroclaw University of Science and Technology, Wybrzeze Wyspianskiego 27, 50-370 Wroclaw, Poland
*
Author to whom correspondence should be addressed.
Materials 2021, 14(17), 4962; https://doi.org/10.3390/ma14174962
Submission received: 21 June 2021 / Revised: 19 August 2021 / Accepted: 25 August 2021 / Published: 31 August 2021
(This article belongs to the Special Issue Monte Carlo Simulations of Polymeric Materials)

Abstract

:
We report the results of the characterization of local Monte Carlo (MC) dynamics of an equilibrium bond fluctuation model polymer matrix (BFM), in time interval typical for MC simulations of non-linear optical phenomena in host-guest systems. The study contributes to the physical picture of the dynamical aspects of quasi-binary mosaic states characterized previously in the static regime. The polymer dynamics was studied at three temperatures (below, above and close to the glass transition), using time-dependent generalization of the static parameters which characterize local free volume and local mobility of the matrix. Those parameters play the central role in the kinetic MC model of host-guest systems. The analysis was done in terms of the probability distributions of instantaneous and time-averaged local parameters. The main result is the characterization of time scales characteristic of various local structural processes. Slowing down effects close to the glass transition are clearly marked. The approach yields an elegant geometric criterion for the glass transition temperature. A simplified quantitative physical picture of the dynamics of guest molecules dispersed in BFM matrix at low temperatures offers a starting point for stochastic modeling of host-guest systems.

1. Introduction

1.1. Bond Fluctuation Model and Simulations of Physical Phenomena in Polymer Systems

Polymer systems are successfully used in various areas of science and technology. While a significant progress has been made in the theoretical statistical physics-based description of polymer physics [1], there still remain essential challenges related, in particular, to the characterization of static and dynamic features of local structure of the polymer matrix. Also, the origin of complex dynamics of physical phenomena in polymer-based systems close to the glass transition temperature remains unclear.
Computer simulations offer complementary theoretical tools for studies of polymer systems. Various simulation techniques are used, depending on the temporal and spatial scales relevant for the specific system [2,3,4].
In those cases when a large number of polymer chains is indispensable, Monte Carlo simulations use simplified models which neglect the detailed atomistic features of polymers. Amongst those models a special role is played by lattice bond fluctuation model (BFM) [5,6,7]. This non-specific model of polymer dynamics has given proof of its capability to predict/simulate a large variety of static and kinetic/dynamic effects in polymer materials like, e.g., properties of linear chains [8,9,10,11,12,13,14,15,16,17,18,19] and polymer rings [20,21], polymer blends and interfaces [22,23], gels and networks [24], localization of glass transition [25,26,27], polymer blends [28], (co)polymers at surfaces [29], polymer brushes in good solvents [30,31,32,33,34,35,36], polymer thin films [37,38,39], equilibrium polymers [17,40,41], general self-assembly [42,43,44], networks and gel point [45,46,47,48], olympic gels [49,50,51], hyperbranched polymers [52,53], dendrimers [54,55,56,57], lipid membranes [58,59,60,61], see also review papers [62,63,64]. Studies of dynamics of BFM chains included, among others, such topics as reptation predictions and scaling of relaxation times [7], Rouse and reptation dynamics at finite temperatures [65], crossover from Rouse to Reptation Dynamics [66], dynamical properties of the slithering-snake algorithm [67], or diffusion of single long polymers [68].
The generalized BFM model, which takes into account the interaction of external electromagnetic fields with guest molecules in a polymer matrix (host-guest system), has offered a theoretical description of a multitude of nonlinear optical phenomena, summarized in a review paper of this Special Issue [69]. This model is based on two concepts: local free volume in polymer matrix, quantified in terms of time-dependent local void parameter V ( r , t ) , and local mobility of the matrix, characterized by time-dependent local mobility parameter C ( r , t ) , see Section 2.2 and Section 2.3 for details. Those parameters are parts of the expressions for transition rates of kinetic/dynamic processes in host-guest systems like, e.g., photoisomerization and reorientation of guest molecules, promoted by light-matter interactions and thermal fluctuations [69], and, correspondingly, play an important role in Monte Carlo simulations of BFM host-guest systems. A detailed probabilistic characterization of random fields V ( r , t ) and C ( r , t ) in terms of adequate probability distributions and correlation functions constitutes one of challenges in the physics of model polymer systems.

1.2. Local Free Volume—Static Characterization of Mosaic-Like States

Figure 1 shows two-dimensional (2D) maps of parameter V ( r , t ) — instantaneous sections of 3D system—for three temperatures: below, close to and above glass transition temperature T g . A typical map is a complex intertwining mosaic (originally proposed a long time ago [70]) of local areas (cells) occupied (red) and unoccupied (green) by monomers (see Section 2.2 for technical details).
The characterization of a typical instantaneous configuration of local free volume V ( r ) in terms of static probability distributions and correlation functions, carried out on large, intermediate and local spatial scales, results in a simplified physical picture of quasi-binary mosaic-like states [69,71]. Approximately 50% of the cells are blocked in the following sense: there is a non-negligible amount of monomers in their closest neighbourhood. On the other hand, about one third cells are free or nearly free of monomers in their close neighbourhood. The latter are not distributed randomly in the space but display a tendency to group into scale-free (fractal) clusters with fractal dimension dependent on the size of a cluster.
An eye-inspection of the three configurations in Figure 1 does not reveal any differences between them—actually, all statistical distributions and correlation functions are practically independent on the temperature. Thus, the structure of mosaic-like states depends, first of all, on the density of the polymer matrix.
We point out that the mosaic-like spatial organization seems to be an inherent feature of dense condensed matter systems, which display some kind of local order in the absence of long-range correlations [72,73].

1.3. Heterogeneous Dynamics and Complexity

The dynamics of guest molecules immersed in the matrix or attached to the polymer chains, depends, among others, on steric hindrances imposed by the polymer matrix, characterized by field V ( r ) . The cells marked red in Figure 1 are blocked and the dynamics of guest molecules which occupy those cells is strongly hindered. This kind of restricted dynamics is referred to as a trapping event. On the contrary, the dynamics of guest molecules in cells marked green depends solely on their interactions with external fields. The resulting dynamics becomes heterogeneous in space, reflecting the heterogeneity of the distribution of local free volume.
The dynamics of the polymer matrix results in the time-dependence of the local free volume: V = V ( r , t ) . It modifies the dynamics of trapping events which become random variables characterized by the spectrum of trapping-time periods. Correspondingly, the dynamics of guest molecules becomes heterogeneous both in space and time. This effect can promote complex behaviour of the host-guest system. Namely, theoretical studies, not related to the polymer physics, lead to the general conclusion that the long-tailed distribution (power law) of trapping-time periods promotes stretched exponential complexity [74,75,76,77,78]. This scenario was demonstrated for the 2D Lennard-Jones liquid close to the melting temperature: the long-tailed distribution of the intervals of trapping events has promoted the stretched exponential behaviour of correlation functions which characterize topological features of mosaic-like states [72,73].
Complex dynamics was reported in recent BFM modeling of all-optical poling [79], Surface Relief Grating inscription [80,81], and dynamics of model polymer chains functionalized with azo-dye molecules in two [82] and three [80] dimensions. It was tentatively ascribed to the physical picture presented above, in which the broad spectrum of temporal scales characterizing the times of life of clusters of blocked and free cells determines an overall dynamics of host-guest systems.

1.4. Objective

The qualitative physical picture of a heterogeneous spatio-temporal dynamics presented above needs quantitative confirmation. So far, only the static spatial heterogeneity was analyzed and, at least partially, confirmed [69,71].
The goal of the paper is to make the first steps towards a detailed probabilistic characterization of the dynamics of local free volume, hypothetically responsible for complex dynamics of physical phenomena in model polymer-based systems. To this end we develop methods and tools which generalize their counterparts used for the characterization of static configurations. Particular emphasis is placed on characterization of the field V ( r , t ) ; other fields, in particular the local mobility field C ( r , t ) , are analyzed to a lesser extent.
The paper is organized as follows. The next Section summarizes briefly the main features of bond fluctuation Monte Carlo model and introduces static and dynamic parameters which characterize locally model polymer matrix. Third Section is devoted to the probabilistic analysis of chosen aspects of Monte Carlo dynamics of local structure of BFM matrix. In the last Section we discuss briefly an emerging physical picture of local dynamics in the polymer matrix.

2. Materials and Methods

The methods and concepts used for the simulations and for an analysis of polymer’s local structure were worked out in detail in the review paper of this Special Issue [69]; in what follows we discuss only their milestones.

2.1. Monte Carlo Bond Fluctuation Model

Bond fluctuation model of a polymer chain stems from the bead-spring model [1]. The chain is composed of effective units called monomers located on a simple cubic lattice, connected by bonds with variable lengths equal (in lattice constants) 2 , 5 , 6 , 3 , 3 , and 10 [5,9,64]. The corresponding potential energies read E = E 0 for the three shortest bonds and E = 0 for the longer ones. Parameter E 0 sets the energy scale and defines the reduced temperature T * = k B T E 0 , where k B and T denote, respectively, Boltzmann constant and an absolute temperature. In what follows T * is referred to as T.
The polymer matrix consisted of N = 24 , 000 polymer chains, each with M = 20 monomers, see Figure 2b. The size of the simulation box read (in lattice constants) V p = 200 3 [2,9,71]; the reduced density ρ = 8 N M / V p = 0.48 characterizes a concentrated polymer solution [9]. Periodic boundary conditions were used.
The polymer matrix was simulated using standard Metropolis Monte Carlo method [2,83]. In a single Monte Carlo Step (MCS) each of the monomers performed a trial move of a unit length along randomly chosen direction x , y or z. The move was accepted under the following conditions [5,69]: (i) the lengths of emerging bonds belonged to the fixed set of lengths, (ii) the move did not violate any steric limitations and (iii) the move was accepted by Metropolis acceptance rule [83]. The configurations of the polymer matrix were prepared in the same way as in Ref. [71]. After the equilibration period the data for further analysis were sampled in the time interval Δ t = 2 · 10 5 MCS. This specific choice of the simulation interval was dictated by the fact that it is sufficient to observe chosen nonlinear optical phenomena in BFM host-guest systems like, e.g., surface relief gratings formation [80,82] ( Δ t 4 · 10 4 MCS), all-optical poling [79] ( Δ t = 2 · 10 5 MCS) and others [71,84,85,86,87,88]. It is important to stress that a different choice of the time of the simulation would influence the results. The simulations were done for three representative temperatures T = 0.1 , 0.25 and 0.5 corresponding, respectively, to glassy phase, close neighbourhood of glass transition and liquid phase. The glass transition temperature was estimated as T g 0.23 0.26 [71,79].
MC method, basically applied for the simulations of equilibrium systems, can be also used to study the dynamics or non-equilibrium systems. Namely, the MC-time evolution of coarse-grained systems does not contradict real dynamics under the condition that the trial moves are limited to local conformational changes [2,89]. Thus, we interpret the results of MC simulations in terms of dynamical parameters.

2.2. Local Void Parameter V ( r , t )

The local void parameter V ( r , t ) characterizes the degree of local inhomogeneity of the distribution of monomers around lattice point r at time t. The quantification of this concept is based on the occupancy of a 3 × 3 × 3 cube with center at lattice point r , see Figure 2a. The central cross consists of the cell at lattice point r and its 6 nearest neighbours. If a monomer occupies a position inside the cross then the central cell is considered as blocked and the value of local void parameter is set to zero: V ( r , t ) = 0 . Otherwise V ( r , t ) = V m a x k , where 0 k V m a x is the number of monomers inside the cube. Thus, parameter V characterizes the strength of steric hindrances—they become important for low values of V and, on the other hand, negligible for its high values. We have used V m a x = 7 , which is the maximal occupancy of a 3 × 3 × 3 cube by monomers found in the simulations [69,71,90]. The cube used in the modeling corresponds approximately to the volume 10 nm 3 for PMMA chains [69,71].

2.3. Local Mobility Parameter C ( r , t )

Local mobility parameter C ( r , t ) characterizes the local dynamics of a polymer matrix around lattice point r at time t in terms of an evolution of the occupancies of lattice sites by monomers inside the cube with center at r , in a single MC step: [69,85,90]
C ( r , t ) = i | n i ( t + 1 ) n i ( t ) | ,
where n i ( t ) denotes the occupancy of i-th cell at time t: n i = 1 when it is occupied by a monomer, otherwise n i = 0 . The sum extends over cells in a 9 × 9 × 9 cube—the large size of the cube is dictated by the fact that monomers in the system are relatively sparse.

2.4. Parameter D ( r , t )

Finally, we introduce the parameter
D r , t = V r , t · C r , t
which combines two features of a polymer matrix: free volume and the dynamics around lattice site r . It characterizes the transition rate for rotational diffusion of guest molecules driven by thermal mobility of the polymer chains [69,79,85,90,91].

3. Results

3.1. Evolution in Time of Local Free Volume: Parameter V ( r , t )

The macroscopic dynamics of guest molecules is determined by two main factors related to the polymer matrix: the distribution of local free volume in space and its temporal evolution. The latter is characterized by the times of life t L ( V , r ) , i.e., the intervals of time in which parameter V remains constant for the cell at r . In particular, t L ( V = 6 ) characterizes the periods of local free dynamics and t L ( V = 0 ) —the periods of trapping events.
The dynamics of local free volume displays a variety of patterns depending on the temperature. Figure 3 shows typical trajectories V ( t ) for MC-time interval of 5 × 10 3 MCS. In spite of the fact that it constitutes only a few percent of the typical time interval in which non-linear optical effects were observed in the simulations [69], it is sufficient for the demonstration of various types of local dynamics of the matrix.
For T = 0.15 , well below T g , most of local free volume is frozen (panel (a), V = 7 and V = 0 , red and green lines, respectively). A low concentration of cells display rapid transitions between locally blocked ( V = 0 ) and (nearly) free regimes ( V = 5 , 6 ). So, there are two typical patterns of local dynamics, corresponding to long (dominating) and short temporal scales (times of life).
Close to the glass transition temperature the dynamics displays a large degree of diversity (panel (b)). Long-lived local free voids ( V = 6 , 7 —black line) are present, as well as long-lived blocked cells ( V = 0 , green line) which change their state twice, for a very short period, to a nearly free one ( V = 6 , 7 ). A very different pattern of the dynamics corresponds to the oscillations between blocked and (nearly) free states (red line). In this case the times of life do not display any typical (regular) pattern. Other kinds of the dynamics (not shown) were also found. We conclude, rather tentatively, that a wide spectrum of temporal scales (times of life), from the shortest to the longest ones, is present.
At higher temperatures ( T = 0.5 ) the dynamical patterns which were found both at low temperature and close to the glass transition are present: oscillations between blocked and free states with very short as well as with long times of life and, on the other hand, long-lasting trapping events with rapid changes of parameter V. We believe that a wide spectrum of times of life is also present in this case. On the other hand, both cases are qualitatively different. This topic is briefly discussed in Section 4.
To quantify the aforementioned qualitative conclusions two approaches were used. The first one concerns the characterization of times of life t L ( V ) of local voids with fixed value of parameter V. To this end, the times of life of such cells were monitored during the whole simulation ( Δ t = 2 · 10 5 MCS) and the empirical histograms, representing the probability distribution ρ ( t L , V ) , were constructed. This kind of analysis provides an information about the maximal as well as typical time scales related to local dynamics of free volume.
Figure 4 shows the plots of ρ ( t L , V ) against MC time for V = 0 and V = 6 . This choice is dictated by their exceptional properties [71,79]—the corresponding cells form the two main components of quasi-binary mosaic-like states. Consider first the case V = 6 (Figure 4a). Such cells are nearly free of monomers in their closest neighbourhood and group into scale-free (fractal) clusters [69,71]. At low temperature ( T = 0.1 ) a smooth decrease of the probability distribution is followed, starting from t L 10 3 MCS, by an irregular behaviour corresponding to longer times of life, up to 2 · 10 5 MCS (i.e., the time of the simulation). We point out that the plot is incomplete—its missing part corresponds to a large concentration (around 44%) of V = 6 cells with t L > 2 · 10 5 MCS. This situation is typical for MC simulations at low temperatures, and is caused by weak thermal fluctuations of the polymer matrix. For a high temperature ( T = 0.5 ) a smooth decay of the probability distribution is present, with a maximal value of t L around 4 · 10 3 MCS. Close to T g ( T = 0.25 ) much longer times of life are present, up to t L 10 5 MCS. The plot displays a new feature: it has two inflection points which mark two different regimes. In the first one the plot matches the plot for T = 0.1 . In the second regime some features characteristic for a power law (over 2 but not 3 decades) are present:
ρ ( t L , T T g ) t L α ,
with α 2.38 . The origin of the power law can be ascribed to the slowing down mechanism.
The probability distributions ρ ( t L , 0 ) for the cells with V = 0 , which constitute the blocked part of the matrix in the quasi-binary mosaic-like states, bear a strong resemblance to their counterparts for V = 6 , with one significant exception: the cells with V = 0 display longer temporal scales than the cells with V = 6 . Namely, the times of life t L reach the values up to 6 · 10 4 MCS for T = 0.5 and up to 2 · 10 5 MCS close to the glass transition temperature. As was the case with V = 6 , at low temperature a vast majority (around 78%) of V = 0 cells are of long-lived ( t L > 2 · 10 5 MCS).
The typical time scales correspond to the times of life t L of the majority of the cells. They can be inferred from the cumulative distribution function
F ( V , t ) = 0 t ρ ( V , t L ) d t L ,
shown in Figure 5 for V = 6 and V = 0 . In both cases the time of life for the majority of the cells does not exceeds 10 3 MCS. For example, at high temperature ( T = 0.5 ) 95% of the cells with V = 6 have t L < 70 MCS and for 99% of the cells t L < 200 MCS. For V = 0 one observes a slight increase of the corresponding times of life—200 MCS and 300 MCS, respectively. Close to the glass transition temperature the times of life are approximately three times larger.
The second method of the characterization of local dynamics of parameter V concerns its evolution in a fixed cell r , in the full simulation period. To this end we define the time average of parameter V ( r , t ) :
V r ¯ = 1 Δ t k = 1 Δ t V r , t k .
Figure 6 shows normalized empirical histograms which represent the probability density ρ ( V r ¯ ) , calculated from all the cells in the system. Clearly, three distinct regimes are present. At low temperatures ( T = 0.1 , Figure 6a), the distribution displays peaks for integer values of V, with negligibly small bars otherwise. The heights of the peaks reproduce the probability distribution ρ S ( V ) of random variable in a static configuration V [69,71]. Correspondingly, the vast majority of the cells are frozen in the time interval Δ t , reflecting, in particular, the presence of a very long tail of times of life for V = 6 discussed above. Close to the glass transition ( T = 0.25 , Figure 6b) a new feature appears: the probability distribution becomes smooth, with the exception of the peak at V = 0 . The remaining part displays a flat maximum in the interval V = 4 5 . Finally, in the high temperature regime ( T = 0.5 , Figure 6c) the probability density becomes nearly gaussian, with a maximum at V 3.2 and standard deviation σ 0.47 .
The latter result offers a semi-quantitative estimation based on the central limit theorem [92], of a typical period of time in which the parameter V remains constant in a cell at high temperature ( T = 0.5 ). Namely, if the realizations of a random variable V in subsequent MCS were independent, the probability density ρ ( V r ¯ ) would be gaussian with the standard deviation σ S / Δ t 0.01 , where σ S = 2.7 denotes the standard deviation for the static probability distribution ρ S ( V ) . The actual value σ 0.47 is much larger, as an effect of the correlations in time between the values of V at a given lattice point. Let us assume, for simplicity, that V remains constant in time interval τ , and that the random variables V in two such intervals are statistically independent. This simple (rather naive) model yields the following estimation of parameter τ . The number N G of independent random variables which, when added together, yield the gaussian probability density, satisfies the equation σ = σ S / N G . One finds N G 34 , which yields an estimation of the block length in which V remains constant: τ = Δ t / N G 7500 MCS.
A short discussion of the time scales reported above and their relation to the characteristic times in which non-linear effects are observed in the MC simulations of host-guest systems can be found in Section 4.

3.2. Evolution in Time of Local Mobility Parameter C ( r , t )

Figure 7 shows exemplary instantaneous 2-dimensional maps of parameter C for three temperatures: T = 0.1 , 0.25 and 0.5 . At low temperature the vast majority of the cells, referred to as immobile, are not influenced by the thermal motions of the chains. The mobile cells, i.e., those with a higher mobility of the chains in their neighbourhood, are rather sporadic and scattered throughout the system. Close to the glass transition the map matches the low-temperature map with the only difference that the concentration of mobile cells is slightly higher and some traces of spatial correlations between them are present. At high temperature mobile cells dominate the system and form a complex spatially correlated structure.
Those observations are fully supported by empirical probability distributions ρ ( C ) calculated from a single realization of the configuration C ( r , t ) , shown in Figure 8. At low and intermediate temperatures, the dominant mobility corresponds to C = 2 , while at a higher temperature a wider spectrum of values C (up to C = 15 ) is present.
The local dynamics of parameter C in a fixed cell r , in the full simulation period, was studied in the same way as for parameter V. Namely, we define the time average of parameter C ( r , t ) :
C r ¯ = 1 Δ t k = 1 Δ t C r , t k .
Figure 9 shows the normalized empirical probability distributions of C r ¯ . Unlike the case with parameter V, the low-temperature probability distribution is rather smooth and does not reproduce the probability mass function ρ ( C ) . The difference between the two cases stems from different sizes of cubes used in the definitions for V and C. Namely, since the latter is much larger, a change of parameter C is not obligatorily accompanied by the change of local free volume V. Close to the glass transition temperature an interesting effect appears—the probability distribution clearly differs from the low-temperature distribution, in spite of the fact that the corresponding probability distributions ρ ( C ) are very similar (Figure 8). It reflects the presence of different time scales in both cases.
A semi-quantitative estimation of characteristic time scales was done using the simple approach based on central limit theorem, introduced in previous Section. Close to the glass transition the probability distribution is reasonably well fitted using gaussian function (solid line in Figure 9, middle panel). The number of independent random variables reads N G 15 , which yields an estimation of the block length in which C remains constant: τ = Δ t / N G 1.7 · 10 4 MCS. At high temperature the probability distribution displays a sharp peak accompanied by a low-probability shoulder. The gaussian fitting procedure of the peak (solid line in Figure 9, right panel) yields N G 140 and τ 1800 MCS. It is worth making two comments. Firstly, the characteristic time scale at T = 0.5 is of an order of magnitude shorter than close ( T = 0.25 ) to the glass transition. Two phenomena bring this effect on: a higher mobility at higher temperature and critical slowing down at glass transition, the latter observed in MC simulations of an all-optical poling [79]. Secondly, the time scale at T = 0.5 is of the same order of magnitude (actually it is four times shorter) than the time scale characteristic for parameter V at the same temperature.

3.3. Parameter D ( r , t )

Parameter D ( r , t ) , which characterizes a cumulative effect of local free volume and local mobility on the orientational diffusion of guest molecules at point r at time t, inherits main features of its both components. Low values indicate that the diffusion at point r is hindered, because of a small amount of local free volume or low level of local mobility. On the contrary, strong diffusion is promoted by high values of both parameters. A moderate local diffusion corresponds to the cases when one of the parameters has a high value and the second—a low value, see below for more details.
Figure 10 shows the evolution in time of parameter D ( r , t ) for three representative cells, at three temperatures. In what follows we interpret the results in terms of a rotational diffusion promoted by the temporal evolution of this field. At a low temperature the maximal value of D corresponds to V = 7 and C 5 , thus D m a x 35 . The first cell represents the case without an orientational diffusion. In the second cell a weak diffusion is present in a half of the short simulation period, with D 5 . Third cell displays a moderate diffusion interrupted by short periods of a high activity, with D 20 . Close to the glass transition parameter D m a x reads approximately 7 · 6 = 42 (see Figure 8). A large variety of patterns of temporal evolution are present, including lack of diffusion, moderate/strong diffusion in the simulation period, as well as periods of moderate diffusion separated by periods of strongly hindered diffusion. The case of a high temperature ( D m a x 100 ) is similar, with one exception—no alternating periods of activity/inactivity were observed.
The field D ( r , t ) is strongly inhomogeneous in space, see Figure 11, which shows its instantaneous 2D sections/maps. As expected, we find a striking qualitative similarity with the maps of the field C ( r , t ) shown in Figure 7. The inhomogeneity is characterized, as in the case of parameter C, in terms of the normalized empirical probability distributions ρ ( D ) , shown in Figure 12. They are dominated by the peaks at D = 0 , which originate from the peaks at V = 0 [69,71] and C = 0 (Figure 8). At low temperature and close to the glass transition around 75 80 % of the cells do not support orientational diffusion. At high temperature the D = 0 peak is mainly due to the V = 0 peak which corresponds to approximately 40 % of blocked cells. The probabilities of non-zero values of D are small; the probability distribution displays some features of a uniform distribution in the interval 0 < D < D m a x .
The local dynamics of parameter D in a fixed cell r , in the full simulation period, was studied using the time average of parameter D ( r , t ) :
D r ¯ = 1 Δ t k = 1 Δ t D r , t k .
Figure 13 shows the normalized empirical probability distributions ρ ( D r ¯ ) . As expected, they inherit some characteristic features of histograms from Figure 6 and Figure 9: a sharp peak at D r ¯ = 0 and an asymmetry of the plots. The former plays an important role in the localization of the glass transition temperature discussed in the next Section.

3.4. Localization of Glass Transition

The MC dynamics of local free volume marks a narrow interval of temperatures around glass transition temperature. Two approaches were used. In the first one we study a global parameter D being an average . . . over the whole system of local parameter D r ¯ :
D = D r ¯ = 1 L 3 i = 1 L 3 D r i ¯ .
Figure 14 shows the plot of temperature dependence of parameter D. Two linear regimes, for low and high temperatures, are separated by a non-linear regime in the interval T = 0.2 0.3 . The linear fits intersect at T = 0.25 , clearly marking the change of the regimes, and providing an estimation of T g . It is in a good agreement with previous estimations of T g based on MC simulations of geometric parameters like square radius of gyration ( T g = 0.26 ) or mean free volume ( T g = 0.25 ), physical parameters ( T g = 0.23 for mean bond energy) and non-linear optical processes (angular hole burning yields T g = 0.25 0.29 , all-optical poling— T g = 0.22 0.25 ) [71,79].
The second approach, based on the probability densities ρ ( V r ¯ ) (Figure 6) and ρ ( D r ¯ ) (Figure 13), offers an elegant, geometric-based criterion for the glass transition temperature. Namely, both plots display sharp peaks located at V r ¯ = D r ¯ = 0 , which corresponds to the blocked ( V = 0 ) part of the matrix. The amplitude of those peaks decreases as the temperature increases, see Figure 15. Both plots, which display a change of convexity at some temperature, were reasonably well fitted using Boltzmann function (Origin 9.5):
f ( T ) = A 2 + A 1 A 2 1 + exp T T 0 d T ,
where A 1 , A 2 , d T , T 0 denote the parameters of the fit. The points of inflection are the solutions of the equation:
d 2 d T 2 f ( T ) = 0 ,
and read T 0 = 0.23 for V r ¯ and T 0 = 0.22 for D r ¯ . Those temperatures are interpreted as the estimations of T g . They are, again, in a good agreement with previous MC-based estimations of T g .

4. Discussion

The study was oriented onto the characterization of local Monte Carlo dynamics/kinetics of an equilibrium BFM polymer matrix in MC-time interval ( 5 25 ) · 10 4 MCS, typical for MC simulations of various non-linear optical phenomena [69]. It constitutes a follow-up of the static characterization recently reported in Refs. [69,71].
The polymer dynamics at point r and time t was characterized at three temperatures (below, above and close to the glass transition), using MC-time dependent parameters V ( r , t ) , C ( r , t ) and D ( r , t ) corresponding, respectively, to local free volume, local mobility, and their product. Special emphasis was put on the dynamics of local volume V ( r , t ) , because it plays the main role in the kinetic model used in MC simulations [69]. To this end, basic probabilistic concepts were used, including probability distributions of various instantaneous and time-averaged parameters. Additional information, of a rather qualitative character, followed from 2D instantaneous maps and short time series.
An extensive discussion of the results was already presented in the main text. In what follows, we discuss them in a wider context.
The paper offers a semi-quantitative support to the physical picture of a spatio-temporal heterogeneity of the polymer (BFM) matrix. The spatial heterogeneity is characterized by instantaneous maps of local free volume (Figure 1), local mobility (Figure 7) and their product (Figure 11), as well as by their probability distributions (Figure 8 and Figure 12). The temporal heterogeneities are illustrated by a highly irregular short time series (Figure 3 and Figure 10).
The main result of the study is the estimation of both maximal and typical time scales characteristic of various structural processes. The general conclusion is that the typical time scales are of an order of magnitude 10 3 MCS, up or down, around two orders of magnitude smaller than the period of the simulations of non-linear optical phenomena. Namely, at a high temperature ( T = 0.5 ) the time of life of a vast majority ( 99 % ) of the V = 6 and V = 0 cells does not exceed a few hundred MCS. In the glass transition region, it becomes approximately three times larger. Similar increase, ascribed to the slowing-down phenomenon, was found in the simulations of an all-optical poling effect [79]. Another approach, based on Central Limit Theorem, provides a rough estimation of the time of life of an arbitrary value of parameters V ( r ) or C ( r ) at cell r . The corresponding time scales read a few thousand MC steps. At a low temperature ( T = 0.1 ) the local dynamics is substantially frozen. 60 % of the cells preserve the values of V during the simulation period, among them 78 % of the population of blocked cells ( V = 0 ). Close to the glass transition the concentration of frozen cells drops down to a few percent. The maximal time-scales correspond to a negligible concentration of the cells, and are substantially larger. The times of life of V = 6 and V = 0 cells constitute, away from the glass transition region, a few thousand and a few ten of thousand MCS, respectively. Close to the glass transition they increase, at least by one order of magnitude, and become comparable to the full period of the simulations.
Another issue of interest is the joint effect of a small amount of local free volume and of a low local mobility at point r , characterized by parameter D ( r , t ) , which promotes or hinders orientational diffusion of guest molecules in MC modeling. We have found that orientational diffusion is fully blocked ( D = 0 ) for 75 80 % of the cells at temperatures up to T g ; at a high temperature half of the cells do not support any orientational diffusion.
The fields V ( r , t ) and D ( r , t ) offer an elegant geometric criterion (inflection point in Figure 15) for the localization of T g , based on those parts of the probability distributions which characterize the absence of rotational diffusion in MC modeling.
It is generally accepted that a polymer matrix is, to some extent, frozen at low temperatures. This study offers a simplified, but nevertheless quantitative physical picture of the dynamics of guest molecules dispersed in the BFM matrix at low temperature ( T = 0.1 ). It is based on the observation that approximately 60 % of the cells preserve their local volume V during the simulation (cases V = 0 , 6 were discussed in Section 3.1). Thus, around 40 % of the cells change local volume V, among them fully blocked cells ( V = 0 ) which, as a rule, after a short time return to fully blocked state. They constitute approximately 15 % of all the cells and, correspondingly, around 75 % of the cells support local dynamics in a fixed local environment in the course of the simulation. Moreover, since around 80 % of the cells do not support any orientational diffusion, its impact on an overall dynamics becomes rather limited. All this leads to the quasi-frozen configuration physical picture in which an orientational diffusion is neglected and the dynamics takes place in time-independent local environment V ( r ) , set by an initial configuration. The active dynamics can be tentatively ascribed to large clusters of cells with V = 6 . This physical picture offers a plausible starting point for a stochastic modeling of the dynamics of guest molecules at low temperatures.
This study yields some preliminary contributions to the physical picture of the dynamical aspects of quasi-binary mosaic states studied previously in the static regime [69,71]. The next steps toward the characterization of random fields (stochastic processes) V ( r , t ) , C ( r , t ) and D ( r , t ) require addressing two main issues. We are of an opinion that the most important one is the study of the cross-correlations of those fields in space and time. The former includes the study of the correlations between times of life of a cell and a large cluster membership, in the context of a hypothetical but plausible clusterization of long-lived areas. The temporal correlations concern, first of all, those at fixed point r . This approach was successfully used in the analysis of structural correlations in two-dimensional model liquids [72]. The second issue is the concept of typical time scales, which should be analyzed in more detail, using advanced mathematical methods for time-series analysis. It includes, in particular, the characterization of times of life of the clusters of cells, especially those with V = 6 . Some of those topics are at progress now.
Finally, let us point out that the time scales, in particular the probability distributions of relaxation times, play one of the central roles in modeling of the complexity which accompany relaxation processes [74]. This study has introduced some parameters and methods which can be useful for a validation of model probabilistic concepts [74] in physical systems. In particular, we have found some traces of temporal heterogeneity in BFM matrix.
This study does not focus on any specific type of polymer system. However, the model might be expanded onto various materials and problems [93,94].

Author Contributions

Conceptualization, A.C.M., W.R. and G.P.; investigation, W.R. and A.C.M.; methodology, A.C.M., G.P. and W.R.; software, W.R. and G.P.; writing—original draft preparation, A.C.M. and W.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Doi, M.; Edwards, S.F. The Theory of Polymer Dynaics; Oxford Univesity Press: New York, NY, USA, 1986. [Google Scholar]
  2. Binder, K. Monte Carlo and Molecular Dynamics Simulations in Polymer Science; Oxford University Press: New York, NY, USA, 1995. [Google Scholar]
  3. Kremer, K. Computer Simulations for Macromolecular Science. Macromol. Chem. Phys. 2003, 204, 257–264. [Google Scholar] [CrossRef]
  4. Binder, K.; Ciccotti, G. Monte Carlo and Molecular Dynamics of Condensed Matter Systems; Italian Physical Society: Bologna, Italy, 1996. [Google Scholar]
  5. Carmesin, I.; Kremer, K. The Bond Fluctuation Method: A New Effective Algorithm for the Dynamics of Polymers in All Spatial Dimensions. Macromolecules 1988, 21, 2819–2823. [Google Scholar] [CrossRef]
  6. Wittmann, H.-P.; Kremer, K. Vectorized version of the bond fluctuation method for lattice polymers. Comp. Phys. Commun. 1990, 61, 309–330. [Google Scholar] [CrossRef]
  7. Paul, W.; Binder, K.; Heermann, D.; Kremer, K. Dynamics of polymer solutions and melts. Reptation predictions and scaling of relaxation Times. J. Chem. Phys. 1991, 95, 7726. [Google Scholar] [CrossRef]
  8. Deutsch, H.; Dickman, R. Equation of state for athermal lattice chains in a 3d fluctuating bond model. J. Chem. Phys. 1990, 93, 8983. [Google Scholar] [CrossRef]
  9. Deutsch, H.-P.; Binder, K. Interdiffusion and self-diffusion in polymer mixtures: A Monte Carlo study. J. Chem. Phys. 1991, 94, 2294–2304. [Google Scholar] [CrossRef]
  10. Paul, W.; Binder, K.; Heermann, D.; Kremer, K. Crossover scaling in semidilute polymer solutions: A Monte Carlo test. J. Phys. II 1991, 1, 37. [Google Scholar] [CrossRef]
  11. Muller, M.; Paul, W. Measuring the chemical potential of polymer solutions and melts in computer simulations. J. Chem. Phys. 1994, 100, 719. [Google Scholar] [CrossRef]
  12. Wilding, N.; Muller, M. Accurate measurements of the chemical potential of polymeric systems by Monte Carlo simulation. J. Chem. Phys. 1994, 101, 4324. [Google Scholar] [CrossRef] [Green Version]
  13. Muller, M.; Binder, K.; Schafer, L. Intra- and Interchain Correlations in Semidilute Polymer Solutions: Monte Carlo Simulations and Renormalization Group Results. Macromolecules 2000, 33, 4568. [Google Scholar] [CrossRef]
  14. Stukan, M.; Ivanov, V.; Muller, M.; Paul, W.; Binder, K. Finite size effects in pressure measurements for Monte Carlo simulations of lattice polymer models. J. Chem. Phys. 2002, 117, 9934. [Google Scholar] [CrossRef]
  15. Wittmer, J.P.; Meyer, H.; Baschnagel, J.; Johner, A.; Obukhov, S.P.; Mattioni, L.; Muller, M.; Semenov, A.N. Long Range Bond-Bond Correlations in Dense Polymer Solutions. Phys. Rev. Lett. 2004, 93, 147801. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Wittmer, J.P.; Beckrich, P.; Johner, A.; Semenov, A.N.; Obukhov, S.P.; Meyer, H.; Baschnagel, J. Why polymer chains in a melt are not random walks. Europhys. Lett. 2007, 77, 56003. [Google Scholar] [CrossRef] [Green Version]
  17. Beckrich, P.; Johner, A.; Semenov, A.N.; Obukhov, S.P.; Benoit, H.; Wittmer, J.P. Intramolecular Form Factor in Dense Polymer Systems: Systematic Deviations from the Debye Formula. Macromolecules 2007, 40, 3805. [Google Scholar] [CrossRef] [Green Version]
  18. Wittmer, J.P.; Beckrich, P.; Meyer, H.; Cavallo, A.; Johner, A.; Baschnagel, J.; Johner, A.; Semenov, A.N.; Obukhov, S.P.; Benoit, H. Intramolecular long-range correlations in polymer melts: The segmental size distribution and its moments. Phys. Rev. E 2007, 76, 011803. [Google Scholar] [CrossRef] [Green Version]
  19. Meyer, H.; Wittmer, J.P.; Kreer, T.; Beckrich, P.; Johner, A.; Farago, J.; Baschnagel, J. Static Rouse modes and related quantities: Corrections to chain ideality in polymer melts. Eur. Phys. J. E 2008, 26, 25. [Google Scholar] [CrossRef] [PubMed]
  20. Muller, M.; Wittmer, J.P.; Cates, M.E. Topological effects in ring polymers: A computer simulation study. Phys. Rev. E 1996, 53, 5063. [Google Scholar] [CrossRef] [Green Version]
  21. Muller, M.; Wittmer, J.P.; Barrat, J.-L. On two intrinsic length scales in polymer physics: Topological constraints vs. entanglement length. Europhys. Lett. 2000, 52, 406. [Google Scholar] [CrossRef]
  22. Muller, M. Miscibility behavior and single chain properties in polymer blends: A bond fluctuation model study. Macromol. Theory Simul. 1999, 8, 343. [Google Scholar] [CrossRef]
  23. Cavallo, A.; Muller, M.; Binder, K. Anomalous scaling of the critical temperature of unmixing with chain length for two-dimensional polymer blends. Europhys. Lett. 2003, 61, 214. [Google Scholar] [CrossRef]
  24. Sommer, J.-U.; Saalwachter, K. Segmental order in end-linked polymer networks: A Monte Carlo study. Eur. Phys. J. E 2005, 18, 167. [Google Scholar] [CrossRef]
  25. Binder, K.; Baschnagel, J.; Paul, W. Glass transition of polymer melts: Test of theoretical concepts by computer simulation. Prog. Polym. Sci. 2003, 28, 115. [Google Scholar] [CrossRef]
  26. Baschnagel, J.; Binder, K.; Wittmann, H.-P. The influence of the cooling rate on the glass transition and the glassy state in three-dimensional dense polymer melts: A Monte Carlo study. J. Phys. Condens. Matter 1993, 5, 1597. [Google Scholar] [CrossRef]
  27. Wittmann, H.-P.; Kremer, K.; Binder, K. Glass transition of polymer melts: A two-dimensional Monte Carlo study in the framework of the bond fluctuation method. J. Chem. Phys. 1996, 96, 6291. [Google Scholar] [CrossRef]
  28. Deutsch, H.-P.; Binder, K. Critical Behavior and Crossover Scaling in Symmetric Polyme Mixtures: A Monte Carlo Investigation. Macromolecules 1992, 25, 6214–6230. [Google Scholar] [CrossRef]
  29. Werner, A.; Schmid, F.; Muller, M. Monte Carlo simulations of copolymers at homopolymer interfaces: Interfacial structure as a function of the copolymer density. J. Chem. Phys. 1999, 110, 5370. [Google Scholar] [CrossRef] [Green Version]
  30. Wittmer, J.; Johner, A.; Joanny, J.F.; Binder, K. Chain desorption from a semidilute polymer brush: A Monte Carlo simulation. J. Chem. Phys. 1994, 101, 4379. [Google Scholar] [CrossRef] [Green Version]
  31. Kopf, A.; Baschnagel, J.; Wittmer, J.; Binder, K. On the Adsorption Process in Polymer Brushes: A Monte Carlo Study. Macromolecules 1996, 29, 1433. [Google Scholar] [CrossRef]
  32. Wittmer, J.P.; Cates, M.E.; Johner, A.; Turner, M.S. Diffusive growth of a polymer layer by in situ polymerization. Europhys. Lett. 1996, 33, 397. [Google Scholar] [CrossRef] [Green Version]
  33. Lang, M.; Werner, M.; Dockhorn, R.; Kreer, T. Arm Retraction Dynamics in Dense Polymer Brushes. Macromolecules 2016, 49, 5190–5201. [Google Scholar] [CrossRef]
  34. Lang, M.; Hoffmann, M.; Dockhorn, R.; Werner, M.; Sommer, J.-U. Fluctuation driven height reduction of crosslinked polymer brushes: A Monte Carlo study. J. Chem. Phys. 2013, 139, 164903. [Google Scholar] [CrossRef] [PubMed]
  35. Lai, P.-Y.; Binder, K. Structure and dynamics of grafted polymer layers: A Monte Carlo simulation. J. Chem. Phys. 1991, 95, 9288. [Google Scholar] [CrossRef]
  36. Lai, P.-Y. Grafted polymer layers with chain exchange: A Monte Carlo simulation. J. Chem. Phys. 1993, 98, 669. [Google Scholar] [CrossRef]
  37. Khalatur, P.; Khokhlov, A.; Prokhorova, S.; Sheiko, S.; Moller, M.; Reineker, P.; Shirvanyanz, D. Unusual conformation of molecular cylindrical brushes strongly adsorbed on a flat solid Surface. Eur. Phys. J. E 2000, 1, 99. [Google Scholar] [CrossRef]
  38. Mischler, C.; Baschnagel, J.; Binder, K. Polymer films in the normal-liquid and supercooled state: A review of recent Monte Carlo simulation results. Adv. Colloid Interface Sci. 2001, 94, 197. [Google Scholar] [CrossRef] [Green Version]
  39. Cavallo, A.; Müller, M.; Wittmer, J.P.; Johner, A. Single chain structure in thin polymer films: Corrections to Flory’s and Silberberg’s hypotheses. J. Phys. Condens. Matter. 2005, 17, S1697. [Google Scholar] [CrossRef]
  40. Wittmer, J.P.; Milchev, A.; Cates, M.E. Dynamical Monte Carlo study of equilibrium polymers: Static properties. J. Chem. Phys. 1998, 109, 834. [Google Scholar] [CrossRef] [Green Version]
  41. Wittmer, J.P.; Beckrich, P.; Crevel, F.; Huang, C.C.; Cavallo, A.; Kreer, T.; Meyer, H. Are polymer melts “ideal”? Comput. Phys. Commun. 2007, 177, 146. [Google Scholar] [CrossRef] [Green Version]
  42. Cavallo, A.; Muller, M.; Binder, K. Formation of Micelles in Homopolymer-Copolymer Mixtures: Quantitative Comparison between Simulations of Long Chains and Self- Consistent Field Calculations. Macromolecules 2006, 39, 9539. [Google Scholar] [CrossRef]
  43. Cavallo, A.; Muller, M.; Binder, K. Monte Carlo Simulation of a Homopolymer-Copolymer Mixture Interacting with a Surface: Bulk versus Surface Micelles and Brush Formation. Macromolecules 2008, 41, 4937. [Google Scholar] [CrossRef]
  44. Wengenmayr, M.; Dockhorn, R.; Sommer, J.-U. Multicore Unimolecular Structure Formation in Single Dendritic-Linear Copolymers under Selective Solvent Conditions. Macromolecules 2016, 49, 9215–9227. [Google Scholar] [CrossRef]
  45. Lang, M.; Müller, T. Analysis of the Gel Point of Polymer Model Networks by Computer Simulations. Macromolecules 2020, 53, 498–512. [Google Scholar] [CrossRef]
  46. Muller, T.; Sommer, J.-U.; Lang, M. Tendomers—Force sensitive bis-rotaxanes with jump-like deformation behavior. Soft Matter 2019, 15, 3671–3679. [Google Scholar] [CrossRef] [PubMed]
  47. Rabbel, H.; Breier, P.; Sommer, J.-U. Swelling Behavior of Single-Chain Polymer Nanoparticles: Theory and Simulation. Macromolecules 2017, 50, 7410–7418. [Google Scholar] [CrossRef]
  48. Lang, M.; Schwenke, K.; Sommer, J.-U. Short Cyclic Structures in Polymer Model Networks: A Test of Mean Field Approximation by Monte Carlo Simulations. Macromolecules 2012, 45, 4886–4895. [Google Scholar] [CrossRef]
  49. Lang, M.; Fischer, J.; Werner, M.; Sommer, J.-U. Olympic Gels: Concatenation and Swelling. Macromol. Symp. 2015, 358, 140–147. [Google Scholar] [CrossRef]
  50. Fischer, J.; Lang, M.; Sommer, J.-U. The formation and structure of Olympic gels. J. Chem. Phys. 2015, 143, 243114. [Google Scholar] [CrossRef]
  51. Lang, M.; Fischer, J.; Werner, M.; Sommer, J.-U. Swelling of Olympic Gels. PRL 2014, 112, 238001. [Google Scholar] [CrossRef]
  52. Dockhorn, R.; Pluuschke, L.; Geisler, M.; Zessin, J.; Lindner, P.; Mundil, R.; Merna, J.; Sommer, J.-U.; Lederer, A. Polyolefins Formed by Chain Walking Catalysis—A Matter of Branching Density Only? J. Am. Chem. Soc. 2019, 141, 15586–15596. [Google Scholar] [CrossRef]
  53. Jurjiu, A.; Dockhorn, R.; Mironova, O.; Sommer, J.-U. Two universality classes for random hyperbranched polymers. Soft Matter 2014, 10, 4935–4946. [Google Scholar] [CrossRef] [Green Version]
  54. Wengenmayr, M.; Dockhorn, R.; Sommer, J.-U. Dendrimers in Solution of Linear Polymers: Crowding Effects. Macromolecules 2019, 52, 2616–2626. [Google Scholar] [CrossRef]
  55. Klos, J.S.; Sommer, J.-U. Dendrimer solutions: A Monte Carlo study. Soft Matter 2016, 12, 9007–9013. [Google Scholar] [CrossRef] [Green Version]
  56. Sommer, J.-U.; Klos, J.S.; Mironova, O. Adsorption of branched and dendritic polymers onto flat surfaces: A Monte Carlo study. J. Chem. Phys. 2013, 139, 244903. [Google Scholar] [CrossRef]
  57. Klos, J.S.; Sommer, J.-U. Simulations of Terminally Charged Dendrimers with Flexible Spacer Chains and Explicit Counterions. Macromolecules 2010, 43, 4418–4427. [Google Scholar] [CrossRef]
  58. Checkervarty, A.; Werner, M.; Sommer, J.-U. Formation and stabilization of pores in bilayermembranes by peptide-like amphiphilic polymers. Soft Matter 2018, 14, 2526–2534. [Google Scholar] [CrossRef] [PubMed]
  59. Rabbel, H.; Werner, M.; Sommer, J.-U. Interactions of Amphiphilic Triblock Copolymers with Lipid Membranes: Modes of Interaction and Effect on Permeability Examined by Generic Monte Carlo Simulations. Macromolecules 2015, 48, 4724–4732. [Google Scholar] [CrossRef]
  60. Werner, M.; Sommer, J.-U. Translocation and Induced Permeability of Random AmphiphilicCopolymers Interacting with Lipid Bilayer Membranes. Biomacromolecules 2015, 16, 125–135. [Google Scholar] [CrossRef] [PubMed]
  61. Sommer, J.-U.; Werner, M.; Baulin, V.A. Critical adsorption controls translocation of polymer chains through lipid bilayers and permeation of solvent. EPL 2012, 98, 18003. [Google Scholar] [CrossRef]
  62. Baschnagel, J.; Wittmer, J.P.; Meyer, H. From Synthetic Polymers to Proteins. Comput. Soft Matter 2004, 23, 83. [Google Scholar]
  63. Muller, M. Handbook of Materials Modeling; Springer: New Yorkg, NY, USA, 2005. [Google Scholar]
  64. Wittmer, J.P.; Cavallo, A.; Kreer, T.; Baschnagel, J.; Johner, A. A finite excluded volume bond-fluctuation model: Static properties of dense polimer melts revisited. J. Chem. Phys. 2009, 131, 064901. [Google Scholar] [CrossRef] [Green Version]
  65. Wittmer, J.P.; Paul, W.; Binder, K. Rouse and reptation dynamics at finite temperatures: A Monte Carlo simulation. Macromolecules 1992, 25, 7211. [Google Scholar] [CrossRef]
  66. Kreer, T.; Baschnagel, J.; Muller, M.; Binder, K. Monte Carlo Simulation of Long Chain Polymer Melts: Crossover from Rouse to Reptation Dynamics. Macromolecules 2001, 34, 1105. [Google Scholar] [CrossRef] [Green Version]
  67. Mattioni, L.; Wittmer, J.P.; Baschnagel, J.; Barrat, J.-L.; Luijten, E. Dynamical properties of the slithering-snake algorithm: A numerical test of the activated-reptation hypothesis. Eur. Phys. J. E 2003, 10, 369. [Google Scholar] [CrossRef]
  68. Azuma, R.; Takayama, H. Diffusion of single long polymers in fixed and low density matrix of obstacles confined to two dimensions. J. Chem. Phys. 1999, 111, 8666. [Google Scholar] [CrossRef] [Green Version]
  69. Mitus, A.C.; Saphiannikova, M.; Radosz, W.; Toshchevikov, V.; Pawlik, G. Modeling of Nonlinear Optical Phenomena in Host-Guest Systems Using Bond Fluctuation Monte Carlo Model: A Review. Materials 2021, 14, 1454. [Google Scholar] [CrossRef]
  70. Adam, G.; Gibbs, J.H. On the temperature dependence of cooperative relaxation properties in glass-forming liquids. J. Chem. Phys. 1965, 43, 139–146. [Google Scholar] [CrossRef] [Green Version]
  71. Radosz, W.; Orlik, R.; Pawlik, G.; Mitus, A.C. On complex structure of local free volume in bond fluctuation model of polymer matrix. Polymer 2019, 177, 1–9. [Google Scholar] [CrossRef]
  72. Patashinski, A.Z.; Ratner, M.A.; Grzybowski, B.A.; Orlik, R.; Mitus, A.C. Heterogeneous structure, heterogeneous dynamics, and complex behavior in two-dimensional liquids. J. Phys. Chem. Lett. 2012, 3, 2431–2435. [Google Scholar] [CrossRef] [PubMed]
  73. Patashinski, A.Z.; Ratner, M.A.; Orlik, R.; Mitus, A.C. Nanofluidic Manifestations of Structure in Liquids: A Toy Model. J. Phys. Chem. C 2019, 123, 16787–16795. [Google Scholar] [CrossRef]
  74. Stanislavsky, A.; Weron, K. Stochastic tools hidden behind the empirical dielectric relaxation laws. Rep. Prog. Phys. 2017, 80, 036001. [Google Scholar] [CrossRef] [Green Version]
  75. Metzler, R.; Klafter, J. From stretched exponential to inverse power-law: Franctional dynamics, Cole-Cole relaxation processes, and beyond. J. Non-Cryst. Solids 2002, 305, 81–87. [Google Scholar] [CrossRef]
  76. Vainstein, M.H.; Costa, I.V.; Morgado, R.; Oliveira, F.A. Non-exponential relaxation for anomalous diffusion. EPL 2006, 73, 726–732. [Google Scholar] [CrossRef]
  77. Bohmer, R.; Ngai, K.; Angell, C.; Plazek, D. Nonexponential relaxations in strong and fragile glass formers. J. Chem. Phys. 1993, 99, 4201–4209. [Google Scholar] [CrossRef]
  78. Shlesinger, M.F. Fractal time in condensed matter. Annu. Rev. Phys. Chem. 1988, 39, 269–290. [Google Scholar] [CrossRef]
  79. Radosz, W.; Pawlik, G.; Mitus, A.C. Complex Dynamics of Photo-Switchable Guest Molecules in All-Optical Poling Close to the Glass Transition: Kinetic Monte Carlo Modeling. J. Phys. Chem. B 2018, 122, 1756–1765. [Google Scholar] [CrossRef] [PubMed]
  80. Pawlik, G.; Wysoczanski, T.; Mitus, A.C. Complex Dynamics of Photoinduced Mass Transport and Surface Relief Gratings Formation. Nanomaterials 2019, 9, 352. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Pawlik, G.; Sobolewska, A.; Miniewicz, A.; Mitus, A.C. Generic stochastic Monte Carlo model of the photoinduced mass transport in azo-polymers and fine structure of surface relief gratings. Europhys. Lett. 2014, 105, 26002. [Google Scholar] [CrossRef]
  82. Pawlik, G.; Mitus, A.C. Photoinduced Mass Transport in Azo-Polymers in 2D: Monte Carlo Study of Polarization Effects. Materials 2020, 13, 4724. [Google Scholar] [CrossRef]
  83. Metropolis, N.; Rosenbluth, A.W.; Rosenbluth, M.N.; Teller, A.N.; Teller, E. Equation of state calculations by fast computing machines. J. Chem. Phys. 1953, 21, 1087–1092. [Google Scholar] [CrossRef] [Green Version]
  84. Pawlik, G.; Orlik, R.; Radosz, W.; Mitus, A.C.; Kuzyk, M.G. Towards understanding the photomechanical effect in polymeric fibers: Analysis of free volume in a model polymeric matrix. Proc. SPIE 2008, 2012, 84740A1. [Google Scholar]
  85. Pawlik, G.; Mitus, A.C.; Miniewicz, A.; Sobolewska, A.; Kajzar, F. Temperature dependence of the kinetics of diffraction gratings formation in a polymer matrix containing azobenzene chromophores: Monte Carlo simulations and experiment. Mol. Cryst. Liq. Cryst. 2005, 426, 243–252. [Google Scholar] [CrossRef]
  86. Pawlik, G.; Rau, I.; Kajzar, F.; Mitus, A.C. Second-harmonic generation in poled polymers: Pre-poling history paradigm. Opt. Express 2010, 18, 18793–18804. [Google Scholar] [CrossRef] [Green Version]
  87. Pawlik, G.; Mitus, A.C.; Mysliwiec, J.; Miniewicz, A.; Grote, J.G. Photochromic dye semi-intercalation into DNA-based polymeric matrix: Computer modeling and experiment. Chem. Phys. Lett. 2010, 484, 321–323. [Google Scholar] [CrossRef]
  88. Pawlik, G.; Radosz, W.; Mitus, A.C.; Mysliwiec, J.; Miniewicz, A.; Kajzar, F.; Rau, I. Holographic grating inscription in DR1: DNA-CTMA thin films: The puzzle of time scales. Cent. Eur. J. Chem. 2014, 12, 886–892. [Google Scholar] [CrossRef]
  89. Kawakatsu, K. Statistical Physics of Polymers; Springer: Berlin/Heidelberg, Germany, 2004. [Google Scholar]
  90. Pawlik, G.; Mitus, A.C.; Miniewicz, A.; Kajzar, F. Monte Carlo simulations of temperature dependence of the kinetics of diffraction gratings formation in a polymer matrix containing azobenzene chromophores. J. Non. Opt. Phys. Mat. 2004, 13, 481–489. [Google Scholar] [CrossRef]
  91. Pawlik, G.; Mitus, A.C.; Rau, I.; Kajzar, F. Monte Carlo modelling of chosen non-linear optical effects for systems of guest molecules in polymeric and liquid-crystal matrices. J. Nonlinear Opt. Quantum Opt. 2009, 38, 227–244. [Google Scholar]
  92. Feller, W. An Introduction to Probability Theory and Its Applications; John Wiley and Sons: Hoboken, NJ, USA, 1968. [Google Scholar]
  93. Krauklis, A.E.; Akulichev, A.G.; Gagani, A.I.; Echtermeyer, A.T. Time–Temperature–Plasticization Superposition Principle: Predicting Creep of a Plasticized Epoxy. Polymers 2019, 11, 1848. [Google Scholar] [CrossRef] [Green Version]
  94. Lalwani, S.M.; Eneh, C.I.; Lutkenhaus, J.L. Emerging trends in the dynamics of polyelectrolyte complexes. Phys. Chem. Chem. Phys. 2020, 22, 24157. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Instantaneous 2-dimensional maps of parameter V: T = 0.1 (left), T = 0.25 (close to glass transition, (middle)) and T = 0.5 (right).
Figure 1. Instantaneous 2-dimensional maps of parameter V: T = 0.1 (left), T = 0.25 (close to glass transition, (middle)) and T = 0.5 (right).
Materials 14 04962 g001
Figure 2. (a) Parameter V—a scheme of calculations. (b) Snapshot of the simulation: configuration of polymer system. 10 % of all polymer chains are shown.
Figure 2. (a) Parameter V—a scheme of calculations. (b) Snapshot of the simulation: configuration of polymer system. 10 % of all polymer chains are shown.
Materials 14 04962 g002
Figure 3. Time-dependence of local free volume V for three representative cells at T = 0.15 (a), 0.25 (b) and 0.5 (c), see text for more details.
Figure 3. Time-dependence of local free volume V for three representative cells at T = 0.15 (a), 0.25 (b) and 0.5 (c), see text for more details.
Materials 14 04962 g003
Figure 4. Double logarithmic plot of non-normalized probability distribution ρ ( t L , V ) for V = 6 (a) and V = 0 (b). T = 0.1 (black line), T = 0.25 (red line) and T = 0.5 (blue line). Green dashed line represents power law.
Figure 4. Double logarithmic plot of non-normalized probability distribution ρ ( t L , V ) for V = 6 (a) and V = 0 (b). T = 0.1 (black line), T = 0.25 (red line) and T = 0.5 (blue line). Green dashed line represents power law.
Materials 14 04962 g004
Figure 5. Empirical cumulative functions F ( V , t ) for V = 6 (left) and V = 0 (right). T = 0.25 (blue line) and T = 0.5 (red line). Horizontal lines mark the values 0.95 and 0.99.
Figure 5. Empirical cumulative functions F ( V , t ) for V = 6 (left) and V = 0 (right). T = 0.25 (blue line) and T = 0.5 (red line). Horizontal lines mark the values 0.95 and 0.99.
Materials 14 04962 g005
Figure 6. Normalized probability density ρ ( V r ¯ ) for T = 0.1 (a), T = 0.25 (b) and T = 0.5 (c). Solid line ( T = 0.5 ) represents gaussian fit.
Figure 6. Normalized probability density ρ ( V r ¯ ) for T = 0.1 (a), T = 0.25 (b) and T = 0.5 (c). Solid line ( T = 0.5 ) represents gaussian fit.
Materials 14 04962 g006
Figure 7. Instantaneous 2-dimensional maps of parameter C: T = 0.1 (a), T = 0.25 (b) and T = 0.5 (c).
Figure 7. Instantaneous 2-dimensional maps of parameter C: T = 0.1 (a), T = 0.25 (b) and T = 0.5 (c).
Materials 14 04962 g007
Figure 8. Empirical probability distributions ρ ( C ) : T = 0.1 (a), T = 0.25 (b) and T = 0.5 (c).
Figure 8. Empirical probability distributions ρ ( C ) : T = 0.1 (a), T = 0.25 (b) and T = 0.5 (c).
Materials 14 04962 g008
Figure 9. Normalized probability density ρ ( C r ¯ ) for T = 0.1 (a), T = 0.25 (b) and T = 0.5 (c). Solid lines represent gaussian fits (see text).
Figure 9. Normalized probability density ρ ( C r ¯ ) for T = 0.1 (a), T = 0.25 (b) and T = 0.5 (c). Solid lines represent gaussian fits (see text).
Materials 14 04962 g009
Figure 10. Time-dependence of parameter D in three representative cells for T = 0.15 (a), 0.25 (b) and 0.5 (c).
Figure 10. Time-dependence of parameter D in three representative cells for T = 0.15 (a), 0.25 (b) and 0.5 (c).
Materials 14 04962 g010
Figure 11. Instantaneous 2-dimensional maps of parameter D: T = 0.1 (a), T = 0.25 (b) and T = 0.5 (c).
Figure 11. Instantaneous 2-dimensional maps of parameter D: T = 0.1 (a), T = 0.25 (b) and T = 0.5 (c).
Materials 14 04962 g011
Figure 12. Empirical probability distributions ρ ( D ) : T = 0.1 (a), T = 0.25 (b) and T = 0.5 (c).
Figure 12. Empirical probability distributions ρ ( D ) : T = 0.1 (a), T = 0.25 (b) and T = 0.5 (c).
Materials 14 04962 g012
Figure 13. Normalized empirical probability distributions ρ ( D r ¯ ) for T = 0.1 (a), T = 0.25 (b) and T = 0.5 (c).
Figure 13. Normalized empirical probability distributions ρ ( D r ¯ ) for T = 0.1 (a), T = 0.25 (b) and T = 0.5 (c).
Materials 14 04962 g013
Figure 14. Temperature dependence of parameter D, Equation (8). Solid lines represent linear fits.
Figure 14. Temperature dependence of parameter D, Equation (8). Solid lines represent linear fits.
Materials 14 04962 g014
Figure 15. Plots of the amplitudes of the peaks for V r ¯ = 0 (black circles) and D r ¯ = 0 (blue triangles) against temperature. Red solid and dashed lines represent Boltzmann fits, Equation (9).
Figure 15. Plots of the amplitudes of the peaks for V r ¯ = 0 (black circles) and D r ¯ = 0 (blue triangles) against temperature. Red solid and dashed lines represent Boltzmann fits, Equation (9).
Materials 14 04962 g015
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Radosz, W.; Pawlik, G.; Mituś, A.C. Characterization of Monte Carlo Dynamic/Kinetic Properties of Local Structure in Bond Fluctuation Model of Polymer System. Materials 2021, 14, 4962. https://doi.org/10.3390/ma14174962

AMA Style

Radosz W, Pawlik G, Mituś AC. Characterization of Monte Carlo Dynamic/Kinetic Properties of Local Structure in Bond Fluctuation Model of Polymer System. Materials. 2021; 14(17):4962. https://doi.org/10.3390/ma14174962

Chicago/Turabian Style

Radosz, Wojciech, Grzegorz Pawlik, and Antoni C. Mituś. 2021. "Characterization of Monte Carlo Dynamic/Kinetic Properties of Local Structure in Bond Fluctuation Model of Polymer System" Materials 14, no. 17: 4962. https://doi.org/10.3390/ma14174962

APA Style

Radosz, W., Pawlik, G., & Mituś, A. C. (2021). Characterization of Monte Carlo Dynamic/Kinetic Properties of Local Structure in Bond Fluctuation Model of Polymer System. Materials, 14(17), 4962. https://doi.org/10.3390/ma14174962

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop