-
Surface viscosity in simple liquids
Authors:
Paolo Malgaretti,
Ubaldo Bafile,
Renzo Vallauri,
Pál Jedlovszky,
Marcello Sega
Abstract:
The response of Newtonian liquids to small perturbations is usually considered to be fully described by homogeneous transport coefficients like shear and dilatational viscosity. However, the presence of strong density gradients at the liquid/vapor boundary of fluids hints at the possible existence of an inhomogeneous viscosity. Here, we show that a surface viscosity emerges from the collective dyn…
▽ More
The response of Newtonian liquids to small perturbations is usually considered to be fully described by homogeneous transport coefficients like shear and dilatational viscosity. However, the presence of strong density gradients at the liquid/vapor boundary of fluids hints at the possible existence of an inhomogeneous viscosity. Here, we show that a surface viscosity emerges from the collective dynamics of interfacial layers in molecular simulations of simple liquids. We estimate the surface viscosity to be 8-16 times smaller than that of the bulk fluid at the thermodynamic point considered. This result can have important implications for reactions at liquid surfaces in atmospheric chemistry and catalysis.
△ Less
Submitted 6 March, 2023;
originally announced March 2023.
-
Suspensions of viscoelastic capsules: effect of membrane viscosity on transient dynamics
Authors:
Fabio Guglietta,
Francesca Pelusi,
Marcello Sega,
Othmane Aouane,
Jens Harting
Abstract:
Membrane viscosity is known to play a central role in the transient dynamics of isolated viscoelastic capsules by decreasing their deformation, inducing shape oscillations and reducing the loading time, that is, the time required to reach the steady-state deformation. However, for dense suspensions of capsules, our understanding of the influence of the membrane viscosity is minimal. In this work,…
▽ More
Membrane viscosity is known to play a central role in the transient dynamics of isolated viscoelastic capsules by decreasing their deformation, inducing shape oscillations and reducing the loading time, that is, the time required to reach the steady-state deformation. However, for dense suspensions of capsules, our understanding of the influence of the membrane viscosity is minimal. In this work, we perform a systematic numerical investigation based on coupled immersed boundary -- lattice Boltzmann (IB-LB) simulations of viscoelastic spherical capsule suspensions in the non-inertial regime. We show the effect of the membrane viscosity on the transient dynamics as a function of volume fraction and capillary number. Our results indicate that the influence of membrane viscosity on both deformation and loading time strongly depends on the volume fraction in a non-trivial manner: dense suspensions with large surface viscosity are more resistant to deformation but attain loading times that are characteristic of capsules with no surface viscosity, thus opening the possibility to obtain richer combinations of mechanical features.
△ Less
Submitted 3 July, 2023; v1 submitted 7 February, 2023;
originally announced February 2023.
-
Viscous to inertial coalescence of liquid lenses: a lattice Boltzmann investigation
Authors:
Thomas Scheel,
Qingguang Xie,
Marcello Sega,
Jens Harting
Abstract:
Liquid lens coalescence is an important mechanism involved in many industrial and scientific applications. It has been investigated both theoretically and experimentally, yet it is numerically very challenging to obtain consistent results over the wide ranges of surface tension and viscosity values that are necessary to capture the asymptotic temporal behavior in the viscous and inertial limits. W…
▽ More
Liquid lens coalescence is an important mechanism involved in many industrial and scientific applications. It has been investigated both theoretically and experimentally, yet it is numerically very challenging to obtain consistent results over the wide ranges of surface tension and viscosity values that are necessary to capture the asymptotic temporal behavior in the viscous and inertial limits. We report results of massively parallel simulations based on the color gradient lattice Boltzmann method, which overcome these limitations, and investigate the scaling laws of both regimes. For the two-dimensional case we find good agreement with the similarity solution of the thin-sheet equation, where in the viscous regime the connecting bridge grows linearly with time and in the inertial regime proportionally to $t^{2/3}$. In three dimensions, the viscous growth of the bridge also exhibits a linear time dependence, while in the inertial regime the growth of both the bridge height and the bridge width is proportional to $t^{1/2}$.
△ Less
Submitted 14 June, 2023; v1 submitted 13 January, 2023;
originally announced January 2023.
-
High-density liquid (HDL) adsorption at the supercooled water/vapor interface and its possible relation to the second surface tension inflection point
Authors:
Alexander Gorfer,
Christoph Dellago,
Marcello Sega
Abstract:
We investigate the properties of water along the liquid/vapor coexistence line in the supercooled regime down to the no-man's land. Extensive molecular dynamics simulations of the TIP4P/2005 liquid/vapor interface in the range 198 -- 348 K allow us to locate the second surface tension inflection point with high accuracy at 283$\pm$5 K, close to the temperature of maximum density (TMD). This temper…
▽ More
We investigate the properties of water along the liquid/vapor coexistence line in the supercooled regime down to the no-man's land. Extensive molecular dynamics simulations of the TIP4P/2005 liquid/vapor interface in the range 198 -- 348 K allow us to locate the second surface tension inflection point with high accuracy at 283$\pm$5 K, close to the temperature of maximum density (TMD). This temperature also coincides with the appearance of a density anomaly at the interface known as the apophysis. We relate the emergence of the apophysis to the observation of HDL water adsorption in the proximity of the liquid/vapor interface.
△ Less
Submitted 16 November, 2022;
originally announced November 2022.
-
Resolving the microscopic hydrodynamics at the moving contact line
Authors:
Amal K. Giri,
Paolo Malgaretti,
Dirk Peschka,
Marcello Sega
Abstract:
The molecular structure of moving contact lines (MCLs) and the emergence of a corresponding macroscopic dissipation have made the MCL a paradigm of fluid dynamics. Through novel averaging techniques that remove capillary waves smearing we achieve an unprecedented resolution in molecular dynamics (MD) simulations and find that they match with the continuum description obtained by finite element met…
▽ More
The molecular structure of moving contact lines (MCLs) and the emergence of a corresponding macroscopic dissipation have made the MCL a paradigm of fluid dynamics. Through novel averaging techniques that remove capillary waves smearing we achieve an unprecedented resolution in molecular dynamics (MD) simulations and find that they match with the continuum description obtained by finite element method (FEM) down to molecular scales. This allows us to distinguish dissipation at the liquid-solid interface (Navier-slip) and at the contact line, the latter being negligible for the rather smooth substrate considered.
△ Less
Submitted 1 November, 2022;
originally announced November 2022.
-
Structure and scaling laws of liquid/vapor interfaces close to the critical point
Authors:
György Hantal,
Pál Jedlovszky,
Marcello Sega
Abstract:
To reach a deeper understanding of fluid interfaces it is necessary to identify a meaningful coarse-graining length that separates intrinsic fluctuations from capillary ones, given the lack of a proper statistical mechanical definition of the latter. Here, with the help of unsupervised learning techniques, we introduce a new length scale based on the local density of the fluid. This length scale f…
▽ More
To reach a deeper understanding of fluid interfaces it is necessary to identify a meaningful coarse-graining length that separates intrinsic fluctuations from capillary ones, given the lack of a proper statistical mechanical definition of the latter. Here, with the help of unsupervised learning techniques, we introduce a new length scale based on the local density of the fluid. This length scale follows a scaling law that diverges more mildly than the bulk correlation length upon approaching the critical point. This allows to distinguish regimes of correlated and uncorrelated capillary waves from that of intrinsic fluctuations.
△ Less
Submitted 3 January, 2022;
originally announced January 2022.
-
On the accuracy and performance of the lattice Boltzmann method with 64-bit, 32-bit and novel 16-bit number formats
Authors:
Moritz Lehmann,
Mathias J. Krause,
Giorgio Amati,
Marcello Sega,
Jens Harting,
Stephan Gekle
Abstract:
Fluid dynamics simulations with the lattice Boltzmann method (LBM) are very memory-intensive. Alongside reduction in memory footprint, significant performance benefits can be achieved by using FP32 (single) precision compared to FP64 (double) precision, especially on GPUs. Here, we evaluate the possibility to use even FP16 and Posit16 (half) precision for storing fluid populations, while still car…
▽ More
Fluid dynamics simulations with the lattice Boltzmann method (LBM) are very memory-intensive. Alongside reduction in memory footprint, significant performance benefits can be achieved by using FP32 (single) precision compared to FP64 (double) precision, especially on GPUs. Here, we evaluate the possibility to use even FP16 and Posit16 (half) precision for storing fluid populations, while still carrying arithmetic operations in FP32. For this, we first show that the commonly occurring number range in the LBM is a lot smaller than the FP16 number range. Based on this observation, we develop novel 16-bit formats - based on a modified IEEE-754 and on a modified Posit standard - that are specifically tailored to the needs of the LBM. We then carry out an in-depth characterization of LBM accuracy for six different test systems with increasing complexity: Poiseuille flow, Taylor-Green vortices, Karman vortex streets, lid-driven cavity, a microcapsule in shear flow (utilizing the immersed-boundary method) and finally the impact of a raindrop (based on a Volume-of-Fluid approach). We find that the difference in accuracy between FP64 and FP32 is negligible in almost all cases, and that for a large number of cases even 16-bit is sufficient. Finally, we provide a detailed performance analysis of all precision levels on a large number of hardware microarchitectures and show that significant speedup is achieved with mixed FP32/16-bit.
△ Less
Submitted 31 January, 2022; v1 submitted 16 December, 2021;
originally announced December 2021.
-
Monolayer Structure of Supramolecular Antagonistic Salt Aggregates
Authors:
David Jung,
Jens Harting,
Marcello Sega
Abstract:
The speculated presence of monomolecular lamellae of antagonistic salts in oil-water mixtures has left several open questions besides their hypothetical existence, including their microscopic structure and stabilization mechanism. Here, we simulate the spontaneous formation of supramolecular aggregates of the antagonistic salt sodium tetraphenylborate (NaBPh$_4$) in water and 3-methylpyridine (3-M…
▽ More
The speculated presence of monomolecular lamellae of antagonistic salts in oil-water mixtures has left several open questions besides their hypothetical existence, including their microscopic structure and stabilization mechanism. Here, we simulate the spontaneous formation of supramolecular aggregates of the antagonistic salt sodium tetraphenylborate (NaBPh$_4$) in water and 3-methylpyridine (3-MP) at the atomistic level. We show that, indeed, the lamellae are formed by a monomolecular layer of the anion, enveloped by 3-MP and hydrated sodium counterions. To understand which thermodynamic forces drive the aggregation, we compare the full-atomistic model with a simplified one for the salt and show that the strong hydrophobic effect granted by the large excluded volume of the anion, together with electrostatic repulsion, suffice to explain the stability of the monomolecular lamellae.
△ Less
Submitted 18 March, 2021;
originally announced March 2021.
-
Capillary-bridge Forces Between Solid Particles: Insights from Lattice Boltzmann Simulations
Authors:
Lei Yang,
Marcello Sega,
Jens Harting
Abstract:
Liquid capillary-bridge formation between solid particles has a critical influence on the rheological properties of granular materials and, in particular, on the efficiency of fluidized bed reactors. The available analytical and semi-analytical methods have inherent limitations, and often do not cover important aspects, like the presence of non-axisymmetric bridges. Here, we conduct numerical simu…
▽ More
Liquid capillary-bridge formation between solid particles has a critical influence on the rheological properties of granular materials and, in particular, on the efficiency of fluidized bed reactors. The available analytical and semi-analytical methods have inherent limitations, and often do not cover important aspects, like the presence of non-axisymmetric bridges. Here, we conduct numerical simulations of the capillary bridge formation between equally and unequally-sized solid particles using the lattice Boltzmann method, and provide an assessment of the accuracy of different families of analytical models. We find that some of the models taken into account are shown to perform better than others. However, all of them fail to predict the capillary force for contact angles larger than $π/2$, where a repulsive capillary force attempts to push the solid particle outwards to minimize the surface energy, especially at a small separation distance.
△ Less
Submitted 16 June, 2021; v1 submitted 23 November, 2020;
originally announced November 2020.
-
Lattice Boltzmann simulations of drying suspensions of soft particles
Authors:
Maarten Wouters,
Othmane Aouane,
Marcello Sega,
Jens Harting
Abstract:
The ordering of particles in the drying process of a colloidal suspension is crucial in determining the properties of the resulting film. For example, microscopic inhomogeneities can lead to the formation of cracks and defects that can deteriorate the quality of the film considerably. This type of problem is inherently multiscale and here we study it numerically, using our recently developed metho…
▽ More
The ordering of particles in the drying process of a colloidal suspension is crucial in determining the properties of the resulting film. For example, microscopic inhomogeneities can lead to the formation of cracks and defects that can deteriorate the quality of the film considerably. This type of problem is inherently multiscale and here we study it numerically, using our recently developed method for the simulation of soft polymeric capsules in multicomponent fluids. We focus on the effect of the particle softness on the film microstructure during the drying phase and how it relates to the formation of defects. We quantify the order of the particles by measuring both the Voronoi entropy and the isotropic order parameter. Surprisingly, both observables exhibit a non-monotonic behaviour when the softness of the particles is increased. We further investigate the correlation between the interparticle interaction and the change in the microstructure during the evaporation phase. We observe that the rigid particles form chain-like structures that tend to scatter into small clusters when the particle softness is increased.
△ Less
Submitted 9 June, 2021; v1 submitted 20 October, 2020;
originally announced October 2020.
-
On the weak scaling of the contact distance between two fluctuating interfaces with system size
Authors:
Clemens Moritz,
Marcello Sega,
Max Innerbichler,
Phillip L. Geissler,
Christoph Dellago
Abstract:
A pair of flat parallel surfaces, each freely diffusing along the direction of their separation, will eventually come into contact. If the shapes of these surfaces also fluctuate, then contact will occur when their centers of mass remain separated by a nonzero distance $\ell$. Here we examine the statistics of $\ell$ at the time of first contact for surfaces that evolve in time according to the Ed…
▽ More
A pair of flat parallel surfaces, each freely diffusing along the direction of their separation, will eventually come into contact. If the shapes of these surfaces also fluctuate, then contact will occur when their centers of mass remain separated by a nonzero distance $\ell$. Here we examine the statistics of $\ell$ at the time of first contact for surfaces that evolve in time according to the Edwards-Wilkinson equation. We present a general approach to calculate its probability distribution and determine how its most likely value $\ell^*$ depends on the surfaces' lateral size $L$. We are motivated by an interest in the motion of interfaces between two phases at conditions of thermodynamic coexistence, and in particular the annihilation of domain wall pairs under periodic boundary conditions. Computer simulations of this scenario verify the predicted scaling behavior in two and three dimensions. In the latter case, slow growth where $\ell^\ast$ is an algebraic function of $\log L$ implies that slab-shaped domains remain topologically intact until $\ell$ becomes very small, contradicting expectations from equilibrium thermodynamics.
△ Less
Submitted 20 August, 2020;
originally announced August 2020.
-
Capillary interactions between soft capsules protruding through thin fluid films
Authors:
Maarten Wouters,
Othmane Aouane,
Marcello Sega,
Jens Harting
Abstract:
When a suspension dries, the suspending fluid evaporates, leaving behind a dry film composed of the suspended particles. During the final stages of drying, the height of the fluid film on the substrate drops below the particle size, inducing local interface deformations that lead to strong capillary interactions among the particles. Although capillary interactions between rigid particles are well…
▽ More
When a suspension dries, the suspending fluid evaporates, leaving behind a dry film composed of the suspended particles. During the final stages of drying, the height of the fluid film on the substrate drops below the particle size, inducing local interface deformations that lead to strong capillary interactions among the particles. Although capillary interactions between rigid particles are well studied, much is still to be understood about the behaviour of soft particles and the role of their softness during the final stages of film drying. Here, we use our recently-introduced numerical method that couples a fluid described using the lattice Boltzmann approach to a finite element description of deformable objects to investigate the drying process of a film with suspended soft particles. Our measured menisci deformations and lateral capillary forces, which agree well with previous theoretical and experimental works in case of rigid particles, show the deformations become smaller with increasing particles softness, resulting in weaker lateral interaction forces. At large interparticle distances, the force approaches that of rigid particles. Finally, we investigate the time dependent formation of particle clusters at the late stages of the film drying.
△ Less
Submitted 19 October, 2020; v1 submitted 30 July, 2020;
originally announced July 2020.
-
Ab-initio Structure and Thermodynamics of the RPBE-D3 Water/Vapor Interface by Neural-Network Molecular Dynamics
Authors:
Oliver Wohlfahrt,
Christoph Dellago,
Marcello Sega
Abstract:
Aided by a neural network representation of the density functional theory (DFT) potential energy landscape of water in the RPBE approximation corrected for dispersion, we calculate several structural and thermodynamic properties of its liquid/vapor interface. The neural network speed allows us to bridge the size and time scale gaps required to sample the properties of water along its liquid/vapor…
▽ More
Aided by a neural network representation of the density functional theory (DFT) potential energy landscape of water in the RPBE approximation corrected for dispersion, we calculate several structural and thermodynamic properties of its liquid/vapor interface. The neural network speed allows us to bridge the size and time scale gaps required to sample the properties of water along its liquid/vapor coexistence line with unprecedented precision.
△ Less
Submitted 23 October, 2020; v1 submitted 20 July, 2020;
originally announced July 2020.
-
Atomistic simulation of PDADMAC/PSS oligoelectrolyte multilayers: overall comparison of tri- and tetra-layer systems
Authors:
Pedro A. Sánchez,
Martin Vögele,
Jens Smiatek,
Baofu Qiao,
Marcello Sega,
Christian Holm
Abstract:
By employing large-scale molecular dynamics simulations of atomistically resolved oligoelectrolytes in aqueous solutions, we study in detail the first four layer-by-layer deposition cycles of an oligoelectrolyte multilayer made of poly(diallyl dimethyl ammonium chloride)/poly(styrene sulfonate sodium salt) (PDADMAC/PSS). The multilayers are grown on a silica substrate in 0.1M {NaCl} electrolyte so…
▽ More
By employing large-scale molecular dynamics simulations of atomistically resolved oligoelectrolytes in aqueous solutions, we study in detail the first four layer-by-layer deposition cycles of an oligoelectrolyte multilayer made of poly(diallyl dimethyl ammonium chloride)/poly(styrene sulfonate sodium salt) (PDADMAC/PSS). The multilayers are grown on a silica substrate in 0.1M {NaCl} electrolyte solutions and the swollen structures are then subsequently exposed to varying added salt concentration. We investigated the microscopic properties of the films, analyzing in detail the differences between three- and four-layers systems. Our simulations provide insights on the early stages of growth of a multilayer, which are particularly challenging for experimental observations. We found a rather strong entanglement of the oligoelectrolytes, with a fuzzy layering of the film structure. The main charge compensation mechanism is for all cases intrinsic, whereas extrinsic compensation is relatively ehanced for the layer of the last deposition cycle. In addition, we quantified other fundamental observables of these systems, as the film thickness, water uptake, and overcharge fractions for each deposition layer.
△ Less
Submitted 31 October, 2019;
originally announced October 2019.
-
A simple approximation for the distribution of ions between charged plates in the weak coupling regime
Authors:
Andrei Bazarenko,
Marcello Sega
Abstract:
The solution of the Poisson--Boltzmann equation for counterions confined between two charged plates is known analytically up to a constant, namely, the ion density in the middle of the channel. This quantity is relevant also because it gives access, through the contact theorem, to the osmotic pressure of the system. Here we compare the values of the ion density obtained by numerical and simulation…
▽ More
The solution of the Poisson--Boltzmann equation for counterions confined between two charged plates is known analytically up to a constant, namely, the ion density in the middle of the channel. This quantity is relevant also because it gives access, through the contact theorem, to the osmotic pressure of the system. Here we compare the values of the ion density obtained by numerical and simulation approaches, and report a useful analytic approximation for the weak coupling regime in the absence of added salt, which predicts the value of the ion density in the worst case within 5%. The inclusion of higher order terms in a Laurent expansion can further improve the accuracy, at the expense of simplicity.
△ Less
Submitted 5 September, 2018;
originally announced September 2018.
-
Pressure Profile Calculation with Mesh Ewald Methods
Authors:
Marcello Sega,
Balázs Fábián,
Pál Jedlovszky
Abstract:
The importance of calculating pressure profiles across liquid interfaces is increasingly gaining recognition, and efficient methods for the calculation of long-range contributions are fundamental in addressing systems with a large number of charges. Here, we show how to compute the local pressure contribution for mesh-based Ewald methods, retaining the typical N log N scaling as a function of the…
▽ More
The importance of calculating pressure profiles across liquid interfaces is increasingly gaining recognition, and efficient methods for the calculation of long-range contributions are fundamental in addressing systems with a large number of charges. Here, we show how to compute the local pressure contribution for mesh-based Ewald methods, retaining the typical N log N scaling as a function of the lattice nodes N. This is a considerable improvement on existing methods, which include approximating the electrostatic contribution using a large cut-off and the, much slower, Ewald calculation. As an application, we calculate the contribution to the pressure profile across the water/vapour interface, coming from different molecular layers, both including and removing the effect of thermal capillary waves. We compare the total pressure profile with the one obtained using the cutoff approximation for the calculation of the stresses, showing that the stress distribution obtained by the Harasima and Irving-Kirkwood are quite similar and shifted with respect to each other at most 0.05~nm.
△ Less
Submitted 1 June, 2017;
originally announced June 2017.
-
Long-range dispersion effects on the water/vapor interface simulated using the most common models
Authors:
Marcello Sega,
Christoph Dellago
Abstract:
The long-range contribution to dispersion forces is known to have a major impact on the properties of inhomogeneous fluids, and its correct treatment is increasingly recognized as being a necessary requirement to avoid cutoff-related artefacts. Although analytical corrections for quantities like the surface tension are known, these can not take into account the structural changes induced by the lo…
▽ More
The long-range contribution to dispersion forces is known to have a major impact on the properties of inhomogeneous fluids, and its correct treatment is increasingly recognized as being a necessary requirement to avoid cutoff-related artefacts. Although analytical corrections for quantities like the surface tension are known, these can not take into account the structural changes induced by the long-range contributions. Here, we analyze the interfacial properties of seven popular water models, comparing the results with the cut-off version of the dispersion potential. The differences in surface tension estimates are in all cases found to be less than 2 mN/m.
△ Less
Submitted 10 December, 2016;
originally announced December 2016.
-
Extended friction elucidates the breakdown of fast water transport in graphene oxide membranes
Authors:
A. Montessori,
C. A. Amadei,
G. Falcucci,
M. Sega,
C. D. Vecitis,
S. Succi
Abstract:
The understanding of water transport in graphene oxide (GO) membranes stands out as a major theoretical problem in graphene research. Notwithstanding the intense efforts devoted to the subject in the recent years, a consolidated picture of water transport in GO membranes is yet to emerge. By performing mesoscale simulations of water transport in ultrathin GO membranes, we show that even small amou…
▽ More
The understanding of water transport in graphene oxide (GO) membranes stands out as a major theoretical problem in graphene research. Notwithstanding the intense efforts devoted to the subject in the recent years, a consolidated picture of water transport in GO membranes is yet to emerge. By performing mesoscale simulations of water transport in ultrathin GO membranes, we show that even small amounts of oxygen functionalities can lead to a dramatic drop of the GO permeability, in line with experimental findings. The coexistence of bulk viscous dissipation and spatially extended molecular friction results in a major decrease of both slip and bulk flow, thereby suppressing the fast water transport regime observed in pristine graphene nanochannels. Inspection of the flow structure reveals an inverted curvature in the near-wall region, which connects smoothly with a parabolic profile in the bulk region. Such inverted curvature is a distinctive signature of the coexistence between single-particle Langevin friction and collective hydrodynamics. The present mesoscopic model with spatially extended friction may offer a computationally efficient tool for future simulations of water transport in nanomaterials.
△ Less
Submitted 6 August, 2016;
originally announced August 2016.
-
The importance of chemical potential in the determination of water slip in nanochannels
Authors:
Marcello Sega,
Mauro Sbragaglia,
Luca Biferale,
Sauro Succi
Abstract:
We investigate the slip properties of water confined in graphite-like nano-channels by non-equilibrium molecular dynamics simulations, with the aim of identifying and analyze separately the influence of different physical quantities on the slip length. In a system under confinement but connected to a reservoir of fluid, the chemical potential is the natural control parameter: we show that two nano…
▽ More
We investigate the slip properties of water confined in graphite-like nano-channels by non-equilibrium molecular dynamics simulations, with the aim of identifying and analyze separately the influence of different physical quantities on the slip length. In a system under confinement but connected to a reservoir of fluid, the chemical potential is the natural control parameter: we show that two nanochannels characterized by the same macroscopic contact angle -- but a different microscopic surface potential -- do not exhibit the same slip length unless the chemical potential of water in the two channels is matched. Some methodological issues related to the preparation of samples for the comparative analysis in confined geometries are also discussed.
△ Less
Submitted 28 July, 2015;
originally announced July 2015.
-
Layer-by-layer formation of oligoelectrolyte multilayers: a combined experimental and computational study
Authors:
Samantha Micciulla,
Pedro A. Sánchez,
Jens Smiatek,
Baofu Qiao,
Marcello Sega,
André Laschewsky,
Christian Holm,
Regine von Klitzing
Abstract:
For the first time, the combination of experimental preparation and results of fully atomistic simulations of an oligoelectrolyte multilayer (OEM) made of poly(diallyl dimethyl ammonium chloride)/poly(styrene sulfonate sodium salt) (PDADMAC/PSS) is presented. The layer-by-layer growth was carried out by dipping silica substrates in oligoelectrolyte solutions and was modeled by means of atomistic m…
▽ More
For the first time, the combination of experimental preparation and results of fully atomistic simulations of an oligoelectrolyte multilayer (OEM) made of poly(diallyl dimethyl ammonium chloride)/poly(styrene sulfonate sodium salt) (PDADMAC/PSS) is presented. The layer-by-layer growth was carried out by dipping silica substrates in oligoelectrolyte solutions and was modeled by means of atomistic molecular dynamics simulations with a protocol that mimics the experimental procedure up to the assembly of four layers. Measurements of OEM thickness, surface roughness and amount of adsorbed oligoelectrolyte chains obtained from both approaches are compared. A good agreement between simulated and experimental results was found, with some deviations due to intrinsic limitations of both methods. However, the combination of information extracted from simulations to support the analysis of experimental data can overcome such restrictions and improve the interpretation of experimental results. On the other hand, processes dominated by slower kinetics, like the destabilization of adsorbed layers upon equilibration with the surrounding environment, are out of reach for the simulation modeling approach, but they can be investigated by monitoring in situ the oligoelectrolyte adsorption during the assembly process. This demonstrates how the synergistic use of simulation and experiments improves the knowledge of OEM properties down to the molecular scale.
△ Less
Submitted 23 May, 2014;
originally announced May 2014.
-
Kinetic and ion pairing contributions in the dielectric spectra of electrolyte aqueous solutions
Authors:
Marcello Sega,
Sofia Kantorovich,
Christian Holm,
Axel Arnold
Abstract:
Understanding dielectric spectra can reveal important information about the dynamics of solvents and solutes from the dipolar relaxation times down to electronic ones. In the late 1970s, Hubbard and Onsager predicted that adding salt ions to a polar solution would result in a reduced dielectric permittivity that arises from the unexpected tendency of solvent dipoles to align opposite to the applie…
▽ More
Understanding dielectric spectra can reveal important information about the dynamics of solvents and solutes from the dipolar relaxation times down to electronic ones. In the late 1970s, Hubbard and Onsager predicted that adding salt ions to a polar solution would result in a reduced dielectric permittivity that arises from the unexpected tendency of solvent dipoles to align opposite to the applied field. So far, this effect has escaped an experimental verification, mainly because of the concomitant appearance of dielectric saturation from which the Hubbard-Onsager decrement cannot be easily separated. Here we develop a novel non-equilibrium molecular dynamics simulation approach to determine this decrement accurately for the first time. Using a thermodynamic consistent all-atom force field we show that for an aqueous solution containing sodium chloride around 4.8 Mol/l, this effect accounts for 12\% of the total dielectric permittivity. The dielectric decrement can be strikingly different if a less accurate force field for the ions is used. Using the widespread GROMOS parameters, we observe in fact an {\it increment} of the dielectric permittivity rather than a decrement. We can show that this increment is caused by ion pairing, introduced by a too low dispersion force, and clarify the microscopic connection between long-living ion pairs and the appearance of specific features in the dielectric spectrum of the solution.
△ Less
Submitted 3 April, 2014;
originally announced April 2014.
-
Kinetic dielectric decrement revisited: phenomenology of finite ion concentrations
Authors:
Marcello Sega,
Sofia Kantorovich,
Axel Arnold
Abstract:
With the help of a recently developed non-equilibrium approach, we investigate the ionic strength dependence of the Hubbard--Onsager dielectric decrement. We compute the depolarization of water molecules caused by the motion of ions in sodium chloride solutions from the dilute regime (0.035 M) up close to the saturation concentration (4.24 M), and find that the kinetic decrement displays a strong…
▽ More
With the help of a recently developed non-equilibrium approach, we investigate the ionic strength dependence of the Hubbard--Onsager dielectric decrement. We compute the depolarization of water molecules caused by the motion of ions in sodium chloride solutions from the dilute regime (0.035 M) up close to the saturation concentration (4.24 M), and find that the kinetic decrement displays a strong nonmonotonic behavior, in contrast to the prediction of available models. We introduce a phenomenological modification of the Hubbard--Onsager continuum theory, that takes into account the screening due to the ionic cloud at mean field level, and, is able to describe the kinetic decrement at high concentrations including the presence of a pronounced minimum.
△ Less
Submitted 16 July, 2014; v1 submitted 3 April, 2014;
originally announced April 2014.
-
Mesoscale Structures at Complex Fluid-Fluid Interfaces: a Novel Lattice Boltzmann / Molecular Dynamics Coupling
Authors:
Marcello Sega,
Mauro Sbragaglia,
Sofia Sergeevna Kantorovich,
Alexey Olegovich Ivanov
Abstract:
Complex fluid-fluid interfaces featuring mesoscale structures with adsorbed particles are key components of newly designed materials which are continuously enriching the field of soft matter. Simulation tools which are able to cope with the different scales characterizing these systems are fundamental requirements for efficient theoretical investigations. In this paper we present a novel simulatio…
▽ More
Complex fluid-fluid interfaces featuring mesoscale structures with adsorbed particles are key components of newly designed materials which are continuously enriching the field of soft matter. Simulation tools which are able to cope with the different scales characterizing these systems are fundamental requirements for efficient theoretical investigations. In this paper we present a novel simulation method, based on the approach of Ahlrichs and Dünweg [Ahlrichs and Dünweg, Int. J. Mod. Phys. C, 1998, 9, 1429], that couples the "Shan-Chen" multicomponent Lattice Boltzmann technique to off-lattice molecular dynamics to simulate efficiently complex fluid-fluid interfaces.
We demonstrate how this approach can be used to study a wide class of challenging problems. Several examples are given, with an accent on bicontinuous phases formation in polyelectrolyte solutions and ferrofluid emulsions. We also show that the introduction of solvation free energies in the particle-fluid interaction unveils the hidden, multiscale nature of the particle-fluid coupling, allowing to treat symmetrically (and interchangeably) the on-lattice and off-lattice components of the system.
△ Less
Submitted 19 February, 2014; v1 submitted 4 June, 2013;
originally announced June 2013.
-
Regularization of the slip length divergence in water nanoflows by inhomogeneities at the Angstrom scale
Authors:
Marcello Sega,
Mauro Sbragaglia,
Luca Biferale,
Sauro Succi
Abstract:
We performed non-equilibrium Molecular Dynamics simulations of water flow in nano-channels with the aim of discriminating {\it static} from {\it dynamic} contributions of the solid surface to the slip length of the molecular flow. We show that the regularization of the slip length divergence at high shear rates, formerly attributed to the wall dynamics, is controlled instead by its static properti…
▽ More
We performed non-equilibrium Molecular Dynamics simulations of water flow in nano-channels with the aim of discriminating {\it static} from {\it dynamic} contributions of the solid surface to the slip length of the molecular flow. We show that the regularization of the slip length divergence at high shear rates, formerly attributed to the wall dynamics, is controlled instead by its static properties. Surprisingly, we find that atomic displacements at the Angstrom scale are sufficient to remove the divergence of the slip length and realize the no-slip condition. Since surface thermal fluctuations at room temperature are enough to generate these displacements, we argue that the no-slip condition for water can be achieved also for ideal surfaces, which do not present any surface roughness.
△ Less
Submitted 19 February, 2014; v1 submitted 30 April, 2013;
originally announced April 2013.
-
Identifying interfacial molecules in nonplanar interfaces: the generalized ITIM algorithm
Authors:
Marcello Sega,
Sofia Kantorovich,
Pál Jedlovszky,
Miguel Jorge
Abstract:
We present a generalized version of the ITIM algorithm for the identification of interfacial molecules, which is able to treat arbitrarily shaped interfaces. The algorithm exploits the similarities between the concept of probe sphere used in ITIM and the circumsphere criterion used in the alpha-shapes approach, and can be regarded either as a reference-frame independent version of the former, or a…
▽ More
We present a generalized version of the ITIM algorithm for the identification of interfacial molecules, which is able to treat arbitrarily shaped interfaces. The algorithm exploits the similarities between the concept of probe sphere used in ITIM and the circumsphere criterion used in the alpha-shapes approach, and can be regarded either as a reference-frame independent version of the former, or as an extended version of the latter that includes the atomic excluded volume. The new algorithm is applied to compute the intrinsic orientational order parameters of water around a DPC and a cholic acid micelle in aqueous environment, and to the identification of solvent-reachable sites in four model structures for soot. The additional algorithm introduced for the calculation of intrinsic density profiles in arbitrary geometries proved to be extremely useful also for planar interfaces, as it allows to solve the paradox of smeared intrinsic profiles far from the interface
△ Less
Submitted 10 October, 2012;
originally announced October 2012.
-
Dominant folding pathways of a peptide chain, from ab-initio quantum-mechanical simulations
Authors:
S. a Beccara,
P. Faccioli,
G. Garberoglio,
M. Sega,
F. Pederiva,
H. Orland
Abstract:
Using the Dominant Reaction Pathways method, we perform an ab-initio quantum-mechanical simulation of a conformational transition of a peptide chain. The method we propose makes it possible to investigate the out-of-equilibrium dynamics of these systems, without resorting to an empirical representation of the molecular force field. It also allows to study rare transitions involving rearrangements…
▽ More
Using the Dominant Reaction Pathways method, we perform an ab-initio quantum-mechanical simulation of a conformational transition of a peptide chain. The method we propose makes it possible to investigate the out-of-equilibrium dynamics of these systems, without resorting to an empirical representation of the molecular force field. It also allows to study rare transitions involving rearrangements in the electronic structure. By comparing the results of the ab-initio simulation with those obtained employing a standard force field, we discuss its capability to describe the non-equilibrium dynamics of conformational transitions.
△ Less
Submitted 29 July, 2010;
originally announced July 2010.
-
The ICC* Algorithm: A fast way to include dielectric boundary effects into molecular dynamics simulations
Authors:
Stefan Kesselheim,
Marcello Sega,
Christian Holm
Abstract:
We employ a fast and accurate algorithm to treat dielectric interfaces within molecular dynamics simulations and demonstrate the importance of dielectric boundary forces (DBFs) in two systems of interests in soft-condensed matter science. We investigate a salt solution confined to a slit pore, and a model of a DNA fragment translocating thorugh a narrow pore.
We employ a fast and accurate algorithm to treat dielectric interfaces within molecular dynamics simulations and demonstrate the importance of dielectric boundary forces (DBFs) in two systems of interests in soft-condensed matter science. We investigate a salt solution confined to a slit pore, and a model of a DNA fragment translocating thorugh a narrow pore.
△ Less
Submitted 5 March, 2010;
originally announced March 2010.
-
Influence of pore dielectric boundaries on the translocation barrier of DNA
Authors:
Stefan Kesselheim,
Marcello Sega,
Christian Holm
Abstract:
We investigate the impact of dielectric boundary forces on the translocation process of charged rigid DNA segments through solid neutral nanopores. We assess the electrostatic contribution to the translocation free energy barrier of a model DNA segment by evaluating the potential of mean force in absence and presence of polarization effects by means of coarse-grained molecular dynamics simulatio…
▽ More
We investigate the impact of dielectric boundary forces on the translocation process of charged rigid DNA segments through solid neutral nanopores. We assess the electrostatic contribution to the translocation free energy barrier of a model DNA segment by evaluating the potential of mean force in absence and presence of polarization effects by means of coarse-grained molecular dynamics simulations. The effect of induced polarization charges has been taken into account by employing ICC*, a recently developed algorithm that can efficiently compute induced polarization charges induced on suitably discretized dielectric boundaries. Since water has a higher dielectric constant than the pore walls, polarization effects repel charged objects in the vicinity of the interface, with the effect of significantly increasing the free energy barrier. Another investigated side effect is the change of the counterion distribution around the charged polymer in presence of the induced pore charges. Furthermore we investigate the influence of adding salt to the solution.
△ Less
Submitted 14 February, 2010;
originally announced February 2010.
-
Mesoscopic simulations of the counterion-induced electroosmotic flow - a comparative study
Authors:
Jens Smiatek,
Marcello Sega,
Christian Holm,
Ulf D. Schiller,
Friederike Schmid
Abstract:
We present mesoscopic simulations of the counterion-induced electroosmotic flow in different electrostatic coupling regimes. Two simulation methods are compared, Dissipative Particle Dynamics (DPD) and coupled Lattice-Boltzmann/Molecular Dynamics (LB/MD). A general mapping scheme to match DPD to LB/MD is developed. For the weak-coupling regime, analytic expressions for the flow profiles in the p…
▽ More
We present mesoscopic simulations of the counterion-induced electroosmotic flow in different electrostatic coupling regimes. Two simulation methods are compared, Dissipative Particle Dynamics (DPD) and coupled Lattice-Boltzmann/Molecular Dynamics (LB/MD). A general mapping scheme to match DPD to LB/MD is developed. For the weak-coupling regime, analytic expressions for the flow profiles in the presence of partial-slip as well as no-slip boundary conditions are derived from the Poisson-Boltzmann and Stokes equations, which are in good agreement with the numerical results. The influence of electrofriction and partial slip on the flow profiles is discussed.
△ Less
Submitted 19 May, 2009; v1 submitted 22 January, 2009;
originally announced January 2009.
-
Ensemble Inequivalence in Single Molecule Experiments
Authors:
M. Süzen,
M. Sega,
C. Holm
Abstract:
In bulk systems the calculation of the main thermodynamic quantities leads to the same expectation values in the thermodynamic limit, regardless of the choice of the statistical ensemble. Single linear molecules can be still regarded as statistical systems, where the thermodynamic limit is represented by infinitely long chains. The question of equivalence between different ensembles is not at al…
▽ More
In bulk systems the calculation of the main thermodynamic quantities leads to the same expectation values in the thermodynamic limit, regardless of the choice of the statistical ensemble. Single linear molecules can be still regarded as statistical systems, where the thermodynamic limit is represented by infinitely long chains. The question of equivalence between different ensembles is not at all obvious and has been addressed in the literature, with sometimes contradicting conclusions. We address this problem by studying the scaling properties of the ensemble difference for two different chain models, as a function of the degree of polymerization. By characterizing the scaling behavior of the difference between the isotensional (Gibbs) and isometric (Helmholtz) ensembles in the transition from the low-stretching to the high-stretching regime, we show that ensemble equivalence cannot be reached for macroscopic chains in the low force regime, and we characterize the transition from the inequivalence to the equivalence regime.
△ Less
Submitted 1 May, 2009; v1 submitted 19 October, 2008;
originally announced October 2008.
-
Dominant Reaction Pathways in High Dimensional Systems
Authors:
E. Autieri,
P. Faccioli,
M. Sega,
F. Pederiva,
H. Orland
Abstract:
This paper is devoted to the development of a theoretical and computational framework to efficiently sample the statistically significant thermally activated reaction pathways, in multi-dimensional systems obeying Langevin dynamics. We show how to obtain the set of most probable reaction pathways and compute the corrections due to quadratic thermal fluctuations around such trajectories. We discu…
▽ More
This paper is devoted to the development of a theoretical and computational framework to efficiently sample the statistically significant thermally activated reaction pathways, in multi-dimensional systems obeying Langevin dynamics. We show how to obtain the set of most probable reaction pathways and compute the corrections due to quadratic thermal fluctuations around such trajectories. We discuss how to obtain predictions for the evolution of arbitrary observables and how to generate conformations which are representative of the transition state ensemble. We present an illustrative implementation of our method by studying the diffusion of a point particle in a 2-dimensional funneled external potential.
△ Less
Submitted 27 June, 2008; v1 submitted 2 June, 2008;
originally announced June 2008.
-
Quantitative Protein Dynamics from Dominant Folding Pathways
Authors:
M. Sega,
P. Faccioli,
F. Pederiva,
G. Garberoglio,
H. Orland
Abstract:
We develop a theoretical approach to the protein folding problem based on out-of-equilibrium stochastic dynamics. Within this framework, the computational difficulties related to the existence of large time scale gaps in the protein folding problem are removed and simulating the entire reaction in atomistic details using existing computers becomes feasible. In addition, this formalism provides a…
▽ More
We develop a theoretical approach to the protein folding problem based on out-of-equilibrium stochastic dynamics. Within this framework, the computational difficulties related to the existence of large time scale gaps in the protein folding problem are removed and simulating the entire reaction in atomistic details using existing computers becomes feasible. In addition, this formalism provides a natural framework to investigate the relationships between thermodynamical and kinetic aspects of the folding. For example, it is possible to show that, in order to have a large probability to remain unchanged under Langevin diffusion, the native state has to be characterized by a small conformational entropy. We discuss how to determine the most probable folding pathway, to identify configurations representative of the transition state and to compute the most probable transition time. We perform an illustrative application of these ideas, studying the conformational evolution of alanine di-peptide, within an all-atom model based on the empiric GROMOS96 force field.
△ Less
Submitted 18 January, 2007;
originally announced January 2007.
-
Dominant Pathways in Protein Folding
Authors:
P. Faccioli,
M. Sega,
F. Pederiva,
H. Orland
Abstract:
We present a method to investigate the kinetics of protein folding on a long time-scale and the dynamics underlying the formation of secondary and tertiary structures during the entire reaction. The approach is based on the formal analogy between thermal and quantum diffusion: by writing the solution of the Fokker-Planck equation for the time-evolution of a protein in a viscous heat-bath in term…
▽ More
We present a method to investigate the kinetics of protein folding on a long time-scale and the dynamics underlying the formation of secondary and tertiary structures during the entire reaction. The approach is based on the formal analogy between thermal and quantum diffusion: by writing the solution of the Fokker-Planck equation for the time-evolution of a protein in a viscous heat-bath in terms of a path integral, we derive a Hamilton-Jacobi variational principle from which we are able to compute the most probable pathway of folding. The method is applied to the folding of the Villin Headpiece Subdomain, in the framework of a Go-model. We have found that, in this model, the transition occurs through an initial collapsing phase driven by the starting coil configuration and a later rearrangement phase, in which secondary structures are formed and all computed paths display strong similarities. This method is completely general, does not require the prior knowledge of any reaction coordinate and represents an efficient tool to perfom ab-initio simulations of the entire folding process with available computers.
△ Less
Submitted 27 July, 2006; v1 submitted 24 October, 2005;
originally announced October 2005.