Search

Recent advances in machine learning and their applications have lead to the development of diverse structure-property relationship models for crucial chemical properties, and the solvation free energy is one of them. Here, we introduce a novel ML-based solvation model, which calculates the solvation energy from pairwise atomistic interactions. The novelty of the proposed model consists of a simple architecture: two encoding functions extract atomic feature vectors from the given chemical structure, while the inner product between two atomistic features calculates their interactions. The results on 6,493 experimental measurements achieve outstanding performance and transferability for enlarging training data due to its solvent-non-specific nature. Analysis of the interaction map shows there is a great potential that our model reproduces group contributions on the solvation energy, which makes us believe that the model not only provides the predicted target property but also gives us more detailed physicochemical insights., Comment: 20 pages, 5 figures

Polymer brushes normally swell in a good solvent and collapse in a poor solvent. An abnormal response of polymer brushes, so-called co-nonsolvency, is the phenomenon where the brush counter-intuitively collapses in a good solvent mixture. In this work, we employed molecular dynamics simulations to investigate the structural properties of the grafted polymers on the occurrence of co-nonsolvency. Brushes with various grafting densities were considered to study the effect of topological excluded volumes on the co-nonsolvency. We found that the brush height follows a novel scaling behavior of the grafting density $h \sim \sigma_{\text g}^{0.71}$ in the co-nonsolvent mixture. Using the scaling exponent and Alexander-de Gennes theory, an analytic function that predicts the monomer density was obtained. The many-chain effects in the co-nonsolvent lead to the formation of both intermolecular and intramolecular bridging structures. Increasing the grafting density entails lower looping events occuring because of the intermolcular bridging, causing diverse structural properties. We report how the average thickness, the polymer orientation, and the looping probability varies as the grafting density increases. Based on these observations, we constructed a phase diagram of the polymer brush system using the average thickness and orientation as order parameters. Our simulations and analytical results reveal the nature of co-nonsolvency in polymer brushes in an explicit way and will help to provide practical guidelines for applications such as drug delivery and sensor devices.

We present a coarse-graining strategy for reducing the number of particle species in mixtures to achieve a simpler system with higher diffusion while preserving the total particle number and characteristic dynamic features. As a system of application, we chose the bidisperse Lennard-Jones-like mixture discovered by Kob and Andersen [Phys. Rev. Lett. 73, 1376 (1994)] possessing a slow dynamics due to the fluid's multi-component character with its apparently unconventional choice for the pair potential of the type-A-type-B arrangement. We further established in a so-formed coarse-grained and temperature-independent monodisperse system an equilibrium structure with a radial distribution function resembling its mixture counterpart. This one-component system further possesses similar dynamic features like glass transition temperature and critical exponents while subjected to Newtonian mechanics. This strategy may finally lead to the manufacturing of new nanoparticle/colloidal fluids by experimentally modeling only the outcoming effective pair potential(s) and no other macroscopic quantity., Comment: 19 pages, 16 figures

We numerically investigate the dynamic heterogeneity and its length scale found in the coarse-grained ionic liquid model systems. In our ionic liquid model systems, cations are modeled as dimers with positive charge, while anions are modeled as monomers with negative charge, respectively. To study the effect of the charge distributions on the cations, two ionic liquid models with different charge distributions are used and the model with neutral charge is also considered as a counterpart. To reveal the heterogeneous dynamics in the model systems, we examine spatial distributions of displacement and time distributions of exchange and persistence times. All the models show significant increase of the dynamic heterogeneity as the temperature is lowered. The dynamic heterogeneity is quantified via the well-known four-point susceptibility which measures the fluctuation of a time correlation function. The dynamic correlation length is calculated by fitting the dynamic structure factor, S4(k,t), with Ornstein-Zernike form at the time scale at which the dynamic heterogeneity reaches the maximum value. Obtained time and length scales exhibit a power law relation at the low temperatures, similar to various supercooled liquid models. Especially, the charged model systems show unusual crossover behaviors which are not observed in the uncharged model system. We ascribe the crossover behavior to the enhanced cage effect caused by charges on the particles.

Kinetically constrained models (KCMs) have gained much interest as models that assign the origins of interesting dynamic properties of supercooled liquids to dynamical facilitation mechanisms that have been revealed in many expreiments and numerical simulations. In this work, we investigate the dynamic heterogeneity in the fragile-to-strong liquid via Monte Carlo method using the model that linearly interpolates between the strong-liquid like behavior and the fragile-liquid like behavior by an asymmetry parameter b. When the asymmetry parameter is sufficiently small, smooth fragile-to-strong transition is observed both in the relaxation time and the diffusion constant. Using these physical quantities, we investigate fractional Stokes-Einstein relations observed in this model. When b is fixed, the system shows constant power law exponent under the temperature change, and the exponent has the value between that of the FA model and the East model. Furthermore, we investigate the dynamic length scale of our systems and also find the crossover relation between the relaxation time. We ascribe the competition between energetically favored symmetric relaxation mechanism and entropically favored asymmetric relaxation mechanism to the fragile-to-strong crossover behavior.

Diffusion of long ring polymers in a melt is much slower than the reorganization of their internal structures. While direct evidences for entanglements have not been observed in the long ring polymers unlike linear polymer melts, threading between the rings is suspected to be the main reason for slowing down of ring polymer diffusion. It is, however, difficult to define the threading configuration between two rings because the rings have no chain end. In this work, evidences for threading dynamics of ring polymers are presented by using molecular dynamics simulation and applying a novel analysis method. The simulation results are analyzed in terms of the statistics of persistence and exchange times that have proved useful in studying heterogeneous dynamics of glassy systems. We find that the threading time of ring polymer melts increases more rapidly with the degree of polymerization than that of linear polymer melts. This indicates that threaded ring polymers cannot diffuse until unthreading event occurs, which results in the slowing down of ring polymer diffusion.

We study structural properties of a ring polymeric melt confined in a film in comparison to a linear counterpart using molecular dynamics simulations. Local structure orderings of ring and linear polymers in the vicinity of the surface are similar to each other because the length scale of surface-monomer excluded volume interaction is smaller than the size of an ideal blob of the ring. In a long length scale, while the Silberberg hypothesis can be used to provide a physical origin of confined linear polymer results, it no longer holds for a ring polymer case. We also present different structural properties of ring and linear polymers in a melt, including the size of polymers, an adsorbed amount, and the coordination number of a polymer. Our observation reveals that a confined ring in a melt adopts highly segregated conformation due to a topological excluded volume repulsion, which may provide a new perspective to understand the nature of biological processes, such as territorial segregation of chromosomes in eukaryotic nuclei.

Dynamics of a coarse-grained model for the room-temperature ionic liquid, 1-ethyl-3-methylimidazolium hexafluorophosphate, couched in the united-atom site representation are studied via molecular dynamics simulations. The dynamically heterogeneous behavior of the model resembles that of fragile supercooled liquids. At or close to room temperature, the model ionic liquid exhibits slow dynamics, characterized by nonexponential structural relaxation and subdiffusive behavior. The structural relaxation time, closely related to the viscosity, shows a super-Arrhenius behavior. Local excitations, defined as displacement of an ion exceeding a threshold distance, are found to be mainly responsible for structural relaxation in the alternating structure of cations and anions. As the temperature is lowered, excitations become progressively more correlated. This results in the decoupling of exchange and persistence times, reflecting a violation of the Stokes-Einstein relation., Comment: Published on the Phys. Chem. Chem. Phys. website

To understand the non-exponential relaxation associated with solvation dynamics in the ionic liquid 1-ethyl-3-methylimidazolium hexafluorophosphate, we study power spectra of the fluctuating Franck-Condon energy gap of a diatomic probe solute via molecular dynamics simulations. Results show 1/f dependence in a wide frequency range over 2 to 3 decades, indicating distributed relaxation times. We analyze the memory function and solvation time in the framework of the generalized Langevin equation using a simple model description for the power spectrum. It is found that the crossover frequency toward the white noise plateau is directly related to the time scale for the memory function and thus the solvation time. Specifically, the low crossover frequency observed in the ionic liquid leads to a slowly-decaying tail in its memory function and long solvation time. By contrast, acetonitrile characterized by a high crossover frequency and (near) absence of 1/f behavior in its power spectra shows fast relaxation of the memory function and single-exponential decay of solvation dynamics in the long-time regime., Comment: 10 pages, 4 figures

By applying the concept of dynamical facilitation and analyzing the excitation lines that result from this facilitation, we investigate the origin of decoupling of transport coefficients in supercooled liquids. We illustrate our approach with two classes of models. One depicts diffusion in a strong glass former, and the other in a fragile glass former. At low temperatures, both models exhibit violation of the Stokes-Einstein relation, $D\sim\tau^{-1}$, where $D$ is the self diffusion constant and $\tau$ is the structural relaxation time. In the strong case, the violation is sensitive to dimensionality $d$, going as $D\sim\tau^{-2/3}$ for $d=1$, and as $D\sim \tau^{-0.95}$ for $d=3$. In the fragile case, however, we argue that dimensionality dependence is weak, and show that for $d=1$, $D \sim \tau^{-0.73}$. This scaling for the fragile case compares favorably with the results of a recent experimental study for a three-dimensional fragile glass former., Comment: 7 pages, 7 figures, submitted to Phys. Rev. E