Search

The transport of water and protons in the cathode catalyst layer (CCL) of proton exchange membrane (PEM) fuel cells is critical for cell performance, but the underlying mechanism is still unclear. Herein, the ionomer structure and the distribution/transport characteristics of water and protons in CCLs are investigated via all-atom molecular dynamics simulations. The results show that at low water contents, isolated water clusters form in ionomer pores, while proton transport is mainly via the charged sites of the ionomer side chains and the Grotthuss mechanism. Moreover, with increasing water content, water clusters are interconnected to form continuous water channels, which provide effective paths for proton transfer via the vehicular and Grotthuss mechanisms. Increasing the ionomer mass content can enhance the dense arrangement of the ionomer, which in turn increases the density of charge sites and improves the proton transport efficiency. When the ionomer mass content is high, the clustering effect reduces the space for water diffusion, increases the proton transport path, and finally decreases the proton transport efficiency. By providing physics insights into the proton transport mechanism, this study is helpful for the structural design and performance improvement of CCLs of PEM fuel cells., Comment: 14 pages, 15 figures

We investigate the coalescence of surfactant-laden water droplets by using several different surfactant types and a wide range of concentrations by means of a coarse-grained model obtained by the statistical associating fluid theory. Our results demonstrate in detail a universal mass transport mechanism of surfactant across many concentrations and several surfactant types during the process. Coalescence initiation is seen to occur via a single pinch due to aggregation of surface surfactant, and its remnants tend to become engulfed in part inside the forming bridge. Across the board we confirm the existence of an initial thermal regime with constant bridge width followed by a later inertial regime with bridge width scaling roughly as the square root of time, but see no evidence of an intermediate viscous regime. Coalescence becomes slower as surfactant concentration grows, and we see evidence of the appearance of a further slowdown of a different nature for several times the critical concentration. We anticipate that our results provide further insights in the mechanisms of coalescence of surfactant-laden droplets., Comment: 13 pages, 11 figures

Durotaxis motion is a spectacular phenomenon manifesting itself by the autonomous motion of a nano-object between parts of a substrate with different stiffness. This motion usually takes place along a stiffness gradient from softer to stiffer parts of the substrate. Here, we propose a new design of a polymer brush substrate that demonstrates antidurotaxis droplet motion, that is, droplet motion from stiffer to softer parts of the substrate. By carrying out extensive molecular dynamics simulation of a coarse-grained model, we find that antidurotaxis is solely controlled by the gradient in the grafting density of the brush and is favorable for fluids with a strong attraction to the substrate (low surface energy). The driving force of the antidurotaxial motion is the minimization of the droplet-substrate interfacial energy, which is attributed to the penetration of the droplet into the brush. Thus, we anticipate that the proposed substrate design offers a new understanding and possibilities in the area of autonomous motion of droplets for applications in microfluidics, energy conservation, and biology., Comment: 12 pages, 9 figures

MARTINI is a popular coarse-grained force-field that is mainly used in molecular dynamics (MD) simulations. It is based on the ``Lego'' approach where intermolecular interactions between coarse-grained beads representing chemical units of different polarity are obtained through water--octanol partition coefficients. This enables the simulation of a wide range of molecules by only using a finite number of parametrized coarse-grained beads, similar to the Lego game, where a finite number of bricks are used to create larger structures. Moreover, the MARTINI force-field is based on the Lennard-Jones potential with the shortest possible cutoff including attractions, thus rendering it very efficient for MD simulations. However, MD simulation is in general a computationally expensive method. Here, we demonstrate that using the MARTINI ``Lego'' approach is suitable for many-body dissipative particle (MDPD) dynamics, a method that can simulate multi-component and multi-phase soft matter systems in a much faster time (about 4--7 times) than MD. In this study, a DPPC lipid bilayer is chosen to provide evidence for the validity of this approach and various properties are compared to highlight the potential of the method. Thus, we anticipate that our study opens new possibilities for faster simulations of a wide range of soft matter systems by using the MDPD method., Comment: 9 pages, 5 figures

Droplet coalescence is ubiquitous in nature and the same time key to various technologies, such as inkjet printing. Here, we report on the coalescence of polymer droplets with different chain lengths coalescing on substrates of different wettability. By means of molecular dynamics simulations of a coarse-grained model, it is found that the rate of bridge growth is higher in the case of droplets with smaller contact angles (more wettable substrates) and decreases with the increase of the chain length of the polymers. Different behavior has also been identified in the dynamics of the approach of the two droplets during coalescence with the substrate wettability playing a more important role compared to the chain length of the polymers. While the dynamics of the droplet are greatly affected by the latter parameters, the density profile and flow patterns remain the same for the different cases. Thus, we anticipate that our work provides further insights into the coalescence of liquid polymer droplets on solid substrates with implications for relevant technologies., Comment: 10 pages, 5 figures

Computer simulation is an important tool for scientific progress, especially when lab experiments are either extremely costly and difficult or lack the required resolution. However, all of the simulation methods come with limitations. In molecular dynamics (MD) simulation, the length and time scales that can be captured are limited, while computational fluid dynamics (CFD) methods are built on a range of assumptions, from the continuum hypothesis itself, to a variety of closure assumptions. To address these issues, the coupling of different methodologies provides a way to retain the best of both methods. Here, we provide a perspective on multiscale simulation based on the coupling of MD and CFD with each a distinct part of the simulation domain. This style of coupling allows molecular detail to be present only where it is needed, so CFD can model larger scales than possible with MD alone. We present a unified perspective of the literature, showing the links between state and flux coupling and discuss the various assumptions required for both. A unique challenge in such coupled simulation is obtaining averages and constraining local parts of a molecular simulation. We highlight that incorrect localisation has resulted in an error in the literature for both pressure tensor and coupling constraints. We then finish with some applications, focused on the simulation of fluids. Thus, we hope to motivate further research in this exciting area with applications across the spectrum of scientific disciplines., Comment: 21 pages, 8 figures perspective paper submitted to Physical Chemistry Chemical Physics (PCCP)

The breakup of liquid threads into smaller droplets is a fundamental problem in fluid dynamics. In this study, we estimate the characteristic wavelength of the breakup process by means of many-body dissipative particle dynamics. This wavelength shows a power-law dependence on the Ohnesorge number in line with results from stability analysis. We also discover that the number of satellite droplets exhibits a power-law decay with exponent $0.72 \pm 0.04$ in the product of the Ohnesorge and thermal capillary numbers, while the overall size of main droplets is larger than that based on the characteristic wavelength thanks to the asynchronous breakup of the thread. Finally, we show that the formation of satellite droplets is the result of the advection of pinching points towards the main droplets in a remaining thinning neck, when the velocity gradient of the fluid exhibits two symmetric maxima., Comment: 18 pages, 9 figures

Nano-patterned substrates offer possibilities for controlling the motion of fluids without external energy supply in novel technologies in microfluidics, coatings, etc. Here, we report on the rugotaxial motion of droplets on wrinkled substrates with gradient in the wavelength of the wrinkles by exploring a broad range of parameters, such as amplitude of the wrinkles, substrate wettability, droplet size and wavelength gradient. Adopting a theoretical and molecular dynamics approach, we determine the Cassie-Baxter and Wenzel states of the droplets, investigate the efficiency of rugotaxis as a function of different parameters, and discuss additional effects, such as pinning. We find that shallow wrinkles characterised by small wavelength gradients, and moderate adhesion of the droplet to the substrate favour the rugotaxis motion with growing droplet size, when pinning is avoided. We also find that the driving force in rugotaxis is the gain in interfacial energy between the droplet and the substrate as the droplet enters regions of denser wrinkles (smaller wavelengths of the wrinkles)., Comment: 7 pages, 5 figures

Droplet coalescence is an important process in nature and various technologies (e.g. inkjet printing). Here, we unveil the surfactant mass-transport mechanism and report on several major differences in the coalescence of surfactant-laden droplets as compared to pure water droplets by means of molecular dynamics simulation of a coarse-grained model. Large scale changes to bridge growth dynamics are identified, such as the lack of multiple thermally excited precursors, attenuated collective excitations after contact, slowing down in the inertial regime due to aggregate-induced rigidity and reduced water flow, and a slowing down in the coalescence rate (deceleration) when surfactant concentration increases, while at the same time we also confirm the existence of an initial thermal, and a power-law, inertial, regime of the bridge growth dynamics in both the pure and the surfactant-laden droplets. Thus, we unveil the key mechanisms in one of the fundamental topological processes of liquid droplets containing surfactant, which is crucial in relevant technologies., Comment: 23 pages, 7 figures

Using extensive molecular dynamics simulation of a coarse-grained model, we demonstrate the possibility of sustained unidirectional motion (durotaxis) of droplets without external energy supply when placed on a polymer brush substrate with stiffness gradient in a certain direction. The governing key parameters for the specific substrate design studied, which determine the durotaxis efficiency, are found to be the grafting density of the brush and the droplet adhesion to the brush surface, whereas the strength of the stiffness gradient, the viscosity of the droplet, or the length of the polymer chains of the brush have only a minor effect on the process. It is shown that this durotaxial motion is driven by the steady increase of the interfacial energy between droplet and brush as the droplet moves from softer to stiffer parts of the substrate whereby the mean driving force gradually declines with decreasing roughness of the brush surface. We anticipate that our findings indicate further possibilities in the area of nanoscale motion without external energy supply.

As one of the most fascinating confined water/ice phenomena, two-dimensional square ice has been extensively studied and experimentally confirmed in recent years. Apart from the unidirectional homogeneous square icing patterns considered in previous studies, the multidirectional partitioned square icing patterns are discovered in this study and characterized by molecular dynamics (MD) simulations. Square icing parameters are proposed to quantitatively distinguish the partitioned patterns from the homogeneous patterns and the liquid water. The number of graphene monolayers n is varied in this study, and the results show that it is more energetically favorable to form partitioned square icing patterns when the water molecules are confined between graphite sheets (n >= 2) compared to graphene (n = 1). This phenomenon is insensitive to n as long as n >= 2, because of the short-range nature of the interaction between water molecules and the carbon substrate. Moreover, it is energetically unfavorable to form partitioned square icing patterns for a single layer of water molecules even for n >= 2, verifying that the interaction between layers of water molecules is another dominant factor in the formation of partitioned structures. The conversion from partitioned structure to homogenous square patterns is investigated by changing the pressure and the temperature. Based on the comprehensive MD simulations, this study unveils the formation mechanism of the partitioned square icing patterns., Comment: 10 pages, 9 figures

Atmospheric aerosols can consist of inorganic and organic substances, including surfactants at a significant concentration. Importantly, the latter can reduce the surface tension at the liquid-vapor surfaces, where they preferentially adsorb due to their amphiphilic structure. As a result, processes such as droplet coalescence, development of precipitation and ultimately cloud lifetime, may depend on the presence of surfactants in the aerosols. Here, we present a numerical model for cloud droplet formation, which is based on the Lagrangian particle-based microphysics-scheme super-droplet method and takes into account the presence of surfactant in the droplets. Our results show that surfactant facilitates cloud formation by increasing the number and size of activated droplets, which concentrate at the bottom of the cloud, while the largest droplets are concentrated at the top of the cloud. This indicates a circulation of droplets that involves activation and growth processes from the bottom of the cloud towards the top. Moreover, our conclusions are independent of the particular approach used for modeling the diffusion of Eulerian variables due to the subgrid-scale turbulence. We anticipate that our results will enrich our understanding of the role of surfactants in the behavior of atmospheric aerosols and, importantly, will pave the way for further developments in the numerical modeling of systems with surfactants at macroscopic scales., Comment: 25 pages, 9 figures