RESUMO
In this work, we compare the performance of two structural indicators based on the degree of translational order up to the second coordination shell in three water models: SPC/E, TIP4P/2005, and TIP5P. Beyond directly contrasting their distributions for different temperatures to evidence their usefulness in estimating the fraction of structured and unstructured molecules and, when possible, their classification capability, we also correlate them with an indirect measure of structural constraint: the dynamic propensity. Furthermore, this procedure enables us to show the existence of evident correlations between structural and dynamical information. More specifically, we find that locally structured molecules display a preference for low dynamic propensity values and, more conspicuously, that locally unstructured molecules are extremely subject to high dynamic propensity. This result is particularly relevant for the supercooled regime where the establishment of firm links between the structure and dynamics has remained rather elusive since the occurrence of dynamics that vary in orders of magnitude upon supercooling usually contrast with barely noticeable overall structural changes.
RESUMO
We study the evolution of dynamic fluctuations averaged over different space lengths and time scales to characterize spatially and temporally heterogeneous behavior of TIP4P/2005 water in liquid and supercooled states. Analyzing a 250 000 molecules simulated system, we provide evidence of the existence, upon supercooling, of a significant enhancement of spatially localized dynamic fluctuations stemming from regions of correlated mobile molecules. We show that both the magnitude of the departure from the value expected for the system-size dependence of an uncorrelated system and the system size at which such a trivial regime is finally recovered clearly increase upon supercooling. This provides a means to estimate an upper limit to the maximum length scale of influence of the regions of correlated mobile molecules. Notably, such an upper limit grows two orders of magnitude on cooling, reaching a value corresponding to a few thousand molecules at the lowest investigated temperature.
RESUMO
We study the potential of mean force for pairs of parallel flat surfaces with attractive electrostatic interactions by employing model systems functionalized with different charged, hydrophobic and hydrophilic groups. We study the way in which the local environment (hydrophobic or hydrophilic moieties) modulates the interaction between the attractive charged groups on the plates by removing or attracting nearby water and thus screening or not the electrostatic interaction. To explicitly account for the role of the solvent and the local hydrophobicity, we also perform studies in vacuo. Additionally, the results are compared to that for non-charged plates in order to single out and rationalize the non-additivity of the different non-covalent interactions. Our simulations demonstrate that the presence of neighboring hydrophobic groups promote water removal in the vicinity of the charged groups, thus enhancing charge attraction upon self-assembly. This role of the local hydrophobicity modulating electrostatic interactions is consistent with recent qualitative descriptions in the protein binding context.
RESUMO
By combining the local structure index with potential energy minimisations we study the local environment of the water molecules for a couple of water models, TIP5P-Ew and SPC/E, in order to characterise low- and high-density "species". Both models show a similar behaviour within the supercooled regime, with two clearly distinguishable populations of unstructured and structured molecules, the fraction of the latter increasing with supercooling. Additionally, for TIP5P-Ew, we find that the structured component vanishes quickly at the normal liquid regime (above the melting temperature). Thus, while SPC/E provides a fraction of structured molecules similar to that found in X-ray experiments, we show that TIP5P-Ew underestimates such value. Moreover, unlike SPC/E, we demonstrate that TIP5P-Ew does not follow the linear dependence of the logarithm of the structured fraction with inverse temperature, as predicted by the two-order parameter model. Finally, we link structure to dynamics by showing that there exists a strong correlation between structural fluctuation and dynamics in the supercooled state with spatial correlations in both static and dynamic quantities.
Assuntos
Temperatura , Água/química , Modelos QuímicosRESUMO
Graphene and the graphene-based materials like graphite, carbon nanotubes, and fullerenes are not only usually regarded as hydrophobic but also have been widely employed as paradigms for the investigation of the behavior of water under nonpolar confinement, a question of major concern for fields ranging from biology to materials design. However, some experimental and theoretical insights seem to contradict, at least partially, such a picture. In this work, we will provide firm evidence for a neat hydrophilic nature of graphene surfaces. Our molecular dynamics studies will demonstrate that parallel graphene sheets present a strong tendency to remain fully hydrated for moderately long times (even when the equilibrium state is indeed the collapse of the plates), and thus, they are less prone to self-assembly than model hydrophobic surfaces we shall employ as control which readily undergo a hydrophobic collapse. Potential of mean force calculations will indeed make evident that the solvent exerts a repulsive contribution on the self-assembly of graphene surfaces. Moreover, we shall also quantify graphene hydrophilicity by means of the calculation of water density at two pressures and water density fluctuations. This latter study has never been performed on graphene and represents a means both to confirm and to quantify its neat hydrophilic behavior. We shall also make evident the relevance of the mildly attractive water-carbon interactions, since their artificial weakening will be shown to revert from typically hydrophilic to typically hydrophobic behavior.
RESUMO
In this paper, we investigate how the dielectric constant, ϵ, of an electrolyte solvent influences the current rectification characteristics of bipolar nanopores. It is well recognized that bipolar nanopores with two oppositely charged regions rectify current when exposed to an alternating electric potential difference. Here, we consider dilute electrolytes with NaCl only and with a mixture of NaCl and charged nanoparticles. These systems are studied using two levels of description, all-atom explicit water molecular dynamics (MD) simulations and coarse-grained implicit solvent MD simulations. The charge density and electric potential profiles and current-voltage relationship predicted by the implicit solvent simulations with ϵ = 11.3 show good agreement with the predictions from the explicit water simulations. Under nonequilibrium conditions, the predictions of the implicit solvent simulations with a dielectric constant closer to the one of bulk water are significantly different from the predictions obtained with the explicit water model. These findings are closely aligned with experimental data on the dielectric constant of water when confined to nanometric spaces, which suggests that ϵ decreases significantly compared to its value in the bulk. Moreover, the largest electric current rectification is observed in systems containing nanoparticles when ϵ = 78.8. Using enhanced sampling, we have shown that this larger rectification arises from the presence of a significantly deeper minimum in the free energy of the system with a larger ϵ, and when a negative voltage bias is applied. Since implicit solvent models and mean-field continuum theories are often used to design Janus membranes based on bipolar nanopores, this work highlights the importance of properly accounting for the effects of confinement on the dielectric constant of the electrolyte solvent. The results presented here indicate that the dielectric constant in implicit solvent simulations may be used as an adjustable parameter to approximately account for the effects of nanometric confinement on aqueous electrolyte solvents.
RESUMO
The two-liquids scenario for liquid water assumes the existence of two competing preferential local molecular structural states characterized by either low or high local density. While the former is expected to present good local order thus involving privileged structures, the latter is usually regarded as conforming a high-entropy unstructured state. A main difference in the local arrangement of such "classes" of water molecules can be inferred from the degree of translational order between the first and second molecular shells. This is so, since the low-local-density molecules present a clear gap between the first two shells while in the case of the high-local-density ones, one or more molecules from the second shell have collapsed toward the first one, thus populating the intershell region. Some structural indicators, like the widely employed local structure index and the recently introduced ζ index, have been devised precisely on the basis of this observation, being successful in detecting well-structured low-local-density molecules. However, the nature of the high-local-density state has been mainly disregarded over the years. In this work we employ molecular dynamics simulations for two water models (the extended simple point charge model and the five-site model) at the liquid and supercooled regimes combined with the inherent dynamics approach (energy minimizations of the instantaneous configurations) in order to both rationalize the detailed structural and topological information that these indicators provide and to advance in our understanding of the high-density state.