No Arabic abstract
We use transition path sampling to study evaporation in the SPC/E model of liquid water. Based on thousands of evaporation trajectories, we characterize the members of the transition state ensemble (TSE), which exhibit a liquid-vapor interface with predominantly negative mean curvature at the site of evaporation. We also find that after evaporation is complete, the distributions of translational and angular momenta of the evaporated water are Maxwellian with a temperature equal to that of the liquid. To characterize the evaporation trajectories in their entirety, we find that it suffices to project them onto just two coordinates: the distance of the evaporating molecule to the instantaneous liquid-vapor interface, and the velocity of the water along the average interface normal. In this projected space, we find that the TSE is well-captured by a simple model of ballistic escape from a deep potential well, with no additional barrier to evaporation beyond the cohesive strength of the liquid. Equivalently, they are consistent with a near-unity probability for a water molecule impinging upon a liquid droplet to condense. These results agree with previous simulations and with some, but not all, recent experiments.
We characterize the water repartition within the partially saturated (two-phase) zone (PSZ) during evaporation out of mixed wettable porous media by controlling the wettability of glass beads, their sizes, and as well the surrounding relative humidity. Here, Capillary numbers are low and under these conditions, the percolating front is stabilized by gravity. Using experimental and numerical analyses, we find that the PSZ saturation decreases with the Bond number, where packing of smaller particles have higher saturation values than packing made of larger particles. Results also reveal that the extent (height) of the PSZ, as well as water saturation in the PSZ, both increase with wettability. We also numerically calculate the saturation exclusively contained in connected liquid films and results show that values are less than the expected PSZ saturation. These results strongly reflect that the two-phase zone is not solely made up of connected capillary networks, but also made of disconnected water clusters or pockets. Moreover, we also find that global saturation (PSZ + full wet zone) decreases with wettability, confirming that greater quantity of water is lost via evaporation with increasing hydrophilicity. These results show that connected liquid films are favored in more hydrophilic systems while disconnected water pockets are favored in less hydrophilic systems.
We study the dynamics of flexible, semiflexible, and self-avoiding polymer chains moving under a Kramers metastable potential. Due to thermal noise, the polymers, initially placed in the metastable well, can cross the potential barrier, but these events are extremely rare if the barrier is much larger than thermal energy. To speed up the slow rate processes in computer simulations, we extend the recently proposed path integral hyperdynamics method to the cases of polymers. We consider the cases where the polymers radii of gyration are comparable to the distance between the well bottom and the barrier top. We find that, for a flexible polymer, the crossing rate ($mathcal{R}$) monotonically decreases with chain contour length ($L$), but with the magnitude much larger than the Kramers rate in the globular limit. For a semiflexible polymer, the crossing rate decreases with $L$ but becomes nearly constant for large $L$. For a fixed $L$, the crossing rate becomes maximum at an intermediate bending stiffness. For a self-avoiding chain, the rate is a nonmonotonic function of $L$, first decreasing with $L$, and then, above certain length, increasing with $L$. These findings can be instrumental for efficient separation of biopolymers.
We introduce a lattice model for active nematic composed of self-propelled apolar particles,study its different ordering states in the density-temperature parameter space, and compare with the corresponding equilibrium model. The active particles interact with their neighbours within the framework of the Lebwohl-Lasher model, and move anisotropically along their orientation to an unoccupied nearest neighbour lattice site. An interplay of the activity, thermal fluctuations and density gives rise distinct states in the system. For a fixed temperature, the active nematic shows a disordered isotropic state, a locally ordered inhomogeneous mixed state, and bistability between the inhomogeneous mixed and a homogeneous globally ordered state in different density regime. In the low temperature regime, the isotropic to the inhomogeneous mixed state transition occurs with a jump in the order parameter at a density less than the corresponding equilibrium disorder-order transition density. Our analytical calculations justify the shift in the transition density and the jump in the order parameter. We construct the phase diagram of the active nematic in the density-temperature plane.
Glass transition process gets affected in ultrathin films having thickness comparable to the size of the molecules. We observe systematic broadening of glass transition temperature (Tg) as the thickness of the polymer film reduces below the radius of gyration but the change in the average Tg was found to be very small. Existence of reversible negative and positive thermal expansion below and above Tg increased the sensitivity of our thickness measurements performed using energy dispersive x-ray reflectivity. A simple model of Tg variation as a function of depth expected from sliding motion could explain the results. We observe clear glass transition even for 4 nm polystyrene film that was predicted to be absent from ellipsometry measurements of thicker films.
We study the evolution of an initially random distribution of particles on a square lattice, under certain rules for `growing and `culling of particles. In one version we allow the particles to move laterally along the surface (mobile layer) and in the other version this motion is not allowed (immobile case). In the former case both analytical and computer simulation results are presented, while in the latter only simulation is possible. We introduce growth and culling probabilities appropriate for condensation and evaporation on a two-dimensional surface, and compare results with existing models for this problem. Our results show very interesting behaviour, under certain conditions quite different from earlier models. We find a possibility of hysteresis not reported earlier for such models. ~