Date: Thu, 9 May 2024 | Total: 4
While the linear Pearson correlation coefficient represents a well-established normalized measure to quantify the interrelation of two stochastic variables $X$ and $Y$, it fails for multidimensional variables such as Cartesian coordinates. Avoiding any assumption about the underlying data, the mutual information $I(X, Y)$ does account for multidimensional correlations. However, unlike the normalized Pearson correlation, it has no upper bound ($I \in [0, \infty)$), i.e., it is not clear if say, $I = 0.4$ corresponds to a low or a high correlation. Moreover, the mutual information (MI) involves the estimation of high-dimensional probability densities (e.g., six-dimensional for Cartesian coordinates), which requires a k-nearest neighbor algorithm, such as the estimator by Kraskov et al. [Phys. Rev. E 69, 066138 (2004)]. As existing methods to normalize the MI cannot be used in connection with this estimator, a new approach is presented, which uses an entropy estimation method that is invariant under variable transformations. The algorithm is numerically efficient and does not require more effort than the calculation of the (un-normalized) MI. After validating the method by applying it to various toy models, the normalized MI between the $C_{\alpha}$ -coordinates of T4 lysozyme is considered and compared to a correlation analysis of inter-residue contacts.
Dynamics of a polymer chain in solution gets significantly affected by the temperature and the frictional forces arising due to solvent viscosity. Here, using an explicit solvent framework for polymer simulation with the liberty to tune the solvent viscosity, we study the nonequilibrium dynamics of a flexible homopolymer when it is suddenly quenched from an extended coil state in good solvent to poor solvent conditions. Results from our extensive simulations reveal that depending on the temperature $T$ and solvent viscosity, one encounters long-lived sausage-like intermediates following the usual pearl-necklace intermediates. Use of shape factors of polymers allows us to disentangle these two distinct stages of the overall collapse process, and the corresponding relaxation times. The relaxation time $\tau_s$ of the sausage stage, which is the rate-limiting stage of the overall collapse process, follows an anti-Arrhenius behavior in the high-$T$ limit, and the Arrhenius behavior in the low-$T$ limit. Furthermore, the variation of $\tau_s$ with the solvent viscosity provides evidence of internal friction of the polymer, that modulates the overall collapse significantly, analogous to what is observed for relaxation rates of proteins during their folding. This suggests that the origin of internal friction in proteins is plausibly intrinsic to its polymeric backbone rather than other specifications.
In this work, the univariant two-phase coexistence line of the hydrate of tetrahydrofuran (THF) is determined from 100 to 1000 bar by molecular dynamics simulations. The study is carried out by putting in contact a THF hydrate phase with a stoichiometric aqueous solution phase. Following the direct coexistence technique, the pressure has been fixed, and the coexistence line has been determined by analyzing if the hydrate phase grows or melts at different values of temperature. The model of water used is the well-known TIP4P/Ice model. We have used two different models of THF based on the transferable parameters for phase equilibria-united atom approach (TraPPE-UA), the original (flexible) TraPPe-UA model as well as a rigid and planar version of it. Overall, at high pressures, small differences have been observed in the results obtained by both models. Also, large differences have been observed in the computational efforts required by the simulations performed using both models, being the rigid and planar version much faster than the original one. The effect of the unlike dispersive interactions between the water and THF molecules has been also analyzed at 250 bar using the rigid and planar THF model. In particular, we have modified the Berthelot combining rule by adding a factor ({\xi}O-THF) that modifies the unlike water-THF dispersive interactions and we have analyzed the effect on the dissociation temperature when {\xi}O-THF is modified from 1.0 (original Berthelot combining rule) to 1.4 (modified Berthelot combining rule). We have extended the study using {\xi}O-THF = 1.4 and the rigid THF model to the rest of the pressures considered in this work, finding an excellent agreement with the scarce experimental data taken from the literature.
We consider a flocking system confined transversally between two infinite reflecting parallel walls separated by a distance $L_\perp$. Infinite or periodic boundary conditions are assumed longitudinally to the direction of collective motion, defining a ring geometry typical of experimental realizations with flocking active colloids. Such a confinement selects a flocking state with its mean direction aligned parallel to the wall, thus breaking explicitly the rotational symmetry locally by a boundary effect. Finite size scaling analysis and numerical simulations show that confinement induces an effective mass term ${M_c} \sim L_\perp^{-\zeta}$ (with positive $\zeta$ being the dynamical scaling exponent of the free theory) suppressing scale free correlations at small wave-numbers. However, due to the finite system size in the transversal direction, this effect can only be detected for large enough longitudinal system sizes (i.e. narrow ring geometries). Furthermore, in the longitudinal direction, density correlations are characterized by an anomalous effective mass term. The effective mass term also enhances the global scalar order parameter and suppresses fluctuations of the mean flocking direction. These results suggest an equivalence between transversal confinement and driving by an homogeneous external field, which breaks the rotational symmetry at the global level.