Wetting Behavior of Kerogen Surfaces: Insights from Molecular Dynamics.

In this study, the wettability of a kerogen surface, a key component of shale reservoirs, is investigated by using molecular dynamics simulations. Specifically, we examined the impact of droplet size and morphology as well as surface roughness on the water contact angles. The findings highlighted that the contact angle dependency on the droplet size intensifies with increased rigidity of the surface. Conversely, as the surface becomes more flexible and rougher, it gains hydrophilicity. The higher hydrophilicity stems from the ability of water molecules to penetrate the kerogen corrugations and form more hydrogen bonds with heteroatoms, particularly oxygen. Notably, the contact angle of kerogen hovers between 65 and 75°, thereby crossing the transition from an underoil hydrophilic to an underoil hydrophobic state. Consequently, minor alterations in the kerogen nanostructure can dramatically alter the wetting preference between water and oil. This insight is of paramount significance for refining strategies in managing fluid interactions in shale reservoirs such as geological carbon storage or oil extraction.

Water Structuring Induces Nonuniversal Hydration Repulsion between Polar Surfaces: Quantitative Comparison between Molecular Simulations, Theory, and Experiments.

Polar surfaces in water typically repel each other at close separations, even if they are charge-neutral. This so-called hydration repulsion balances the van der Waals attraction and gives rise to a stable nanometric water layer between the polar surfaces. The resulting hydration water layer is crucial for the properties of concentrated suspensions of lipid membranes and hydrophilic particles in biology and technology, but its origin is unclear. It has been suggested that surface-induced molecular water structuring is responsible for the hydration repulsion, but a quantitative proof of this water-structuring hypothesis is missing. To gain an understanding of the mechanism causing hydration repulsion, we perform molecular simulations of different planar polar surfaces in water. Our simulated hydration forces between phospholipid bilayers agree perfectly with experiments, validating the simulation model and methods. For the comparison with theory, it is important to split the simulated total surface interaction force into a direct contribution from surface-surface molecular interactions and an indirect water-mediated contribution. We find the indirect hydration force and the structural water-ordering profiles from the simulations to be in perfect agreement with the predictions from theoretical models that account for the surface-induced water ordering, which strongly supports the water-structuring hypothesis for the hydration force. However, the comparison between the simulations for polar surfaces with different headgroup architectures reveals significantly different decay lengths of the indirect water-mediated hydration-force, which for laterally homogeneous water structuring would imply different bulk-water properties. We conclude that laterally inhomogeneous water ordering, induced by laterally inhomogeneous surface structures, shapes the hydration repulsion between polar surfaces in a decisive manner. Thus, the indirect water-mediated part of the hydration repulsion is caused by surface-induced water structuring but is surface-specific and thus nonuniversal.

Understanding the "Berg limit": the 65° contact angle as the universal adhesion threshold of biomatter.

Surface phenomena in aqueous environments such as long-range hydrophobic attraction, macromolecular adhesion, and even biofouling are predominantly influenced by a fundamental parameter-the water contact angle. The minimal contact angle required for these and related phenomena to occur has been repeatedly reported to be around 65° and is commonly referred to as the "Berg limit." However, the universality of this specific threshold across diverse contexts has remained puzzling. In this perspective article, we aim to rationalize the reoccurrence of this enigmatic contact angle. We show that the relevant scenarios can be effectively conceptualized as three-phase problems involving the surface of interest, water, and a generic oil-like material that is representative of the nonpolar constituents within interacting entities. Our analysis reveals that attraction and adhesion emerge when substrates display an underwater oleophilic character, corresponding to a "hydrophobicity under oil", which occurs for contact angles above approximately 65°. This streamlined view provides valuable insights into macromolecular interactions and holds implications for technological applications.

Cavitation in lipid bilayers poses strict negative pressure stability limit in biological liquids.

Biological and technological processes that involve liquids under negative pressure are vulnerable to the formation of cavities. Maximal negative pressures found in plants are around -100 bar, even though cavitation in pure bulk water only occurs at much more negative pressures on the relevant timescales. Here, we investigate the influence of small solutes and lipid bilayers, both constituents of all biological liquids, on the formation of cavities under negative pressures. By combining molecular dynamics simulations with kinetic modeling, we quantify cavitation rates on biologically relevant length scales and timescales. We find that lipid bilayers, in contrast to small solutes, increase the rate of cavitation, which remains unproblematically low at the pressures found in most plants. Only when the negative pressures approach -100 bar does cavitation occur on biologically relevant timescales. Our results suggest that bilayer-based cavitation is what generally limits the magnitude of negative pressures in liquids that contain lipid bilayers.

Apparent line tension induced by surface-active impurities.

Line tension in wetting processes is of high scientific and technological relevance, but its understanding remains vague, mainly because it is difficult to determine. A widely used method to extract line tension relies on the variation of a droplet's contact angle with the droplet's size. Such an approach yields the apparent line tension, which is an effective parameter that factors in numerous contributions to the finite-size dependence, thus masking the actual line tension in terms of the excess free energy of the three-phase contact line. Based on our recent computer simulation study, we investigate how small amounts of nonionic surfactants, such as surface-active impurities, contribute to the apparent line tension in aqueous droplets. When depositing polydisperse droplets, their different surface area-to-volume ratios can result in different final bulk concentrations of surfactants, different excess adsorptions to the interfaces, and, consequently, different contact angles. We show that already trace amounts of longer-chained surfactants in a pre-contaminated liquid are enough to affect measurements of the apparent line tension. Our analysis quantifies to what extent "background" impurities, inevitably present in all kinds of experimental settings, limit the resolution of line tension measurements, which is crucial for avoiding data misinterpretation.

Hydrophobicity of Self-Assembled Monolayers of Alkanes: Fluorination, Density, Roughness, and Lennard-Jones Cutoffs.

The interplay of fluorination and structure of alkane self-assembled monolayers and how these affect hydrophobicity are explored via molecular dynamics simulations, contact angle goniometry, and surface-enhanced infrared absorption spectroscopy. Wetting coefficients are found to grow linearly in the monolayer density for both alkane and perfluoroalkane monolayers. The larger contact angles of monolayers of perfluorinated alkanes are shown to be primarily caused by their larger molecular volume, which leads to a larger nearest-neighbor grafting distance and smaller tilt angle. Increasing the Lennard-Jones force cutoff in simulations is found to increase hydrophilicity. Specifically, wetting coefficients scale like the inverse square of the cutoff, and when extrapolated to the infinite cutoff limit, they yield contact angles that compare favorably to experimental values. Nanoscale roughness is also found to reliably increase monolayer hydrophobicity, mostly via the reduction of the entropic part of the work of adhesion. Analysis of depletion lengths shows that droplets on nanorough surfaces partially penetrate the surface, intermediate between Wenzel and Cassie-Baxter states.

Tuning the permeability of regular polymeric networks by the cross-link ratio.

The amount of cross-linking in the design of polymer materials is a key parameter for the modification of numerous physical properties, importantly, the permeability to molecular solutes. We consider networks with a diamond-like architecture and different cross-link ratios, concurring with a wide range of the polymer volume fraction. We particularly focus on the effect and the competition of two independent component-specific solute-polymer interactions, i.e., we distinguish between chain-monomers and cross-linkers, which individually act on the solutes and are altered to cover attractive and repulsive regimes. For this purpose, we employ coarse-grained, Langevin computer simulations to study how the cross-link ratio of polymer networks controls the solute partitioning, diffusion, and permeability. We observe different qualitative behaviors as a function of the cross-link ratio and interaction strengths. The permeability can be tuned ranging over two orders of magnitude relative to the reference bulk permeability. Finally, we provide scaling theories for the partitioning and diffusion that explicitly account for the component-specific interactions as well as the cross-link ratio and the polymer volume fraction. These are in overall good agreement with the simulation results and grant insight into the underlying physics, rationalizing how the cross-link ratio can be exploited to tune the solute permeability of polymeric networks.

Relative humidity in droplet and airborne transmission of disease.

A large number of infectious diseases are transmitted by respiratory droplets. How long these droplets persist in the air, how far they can travel, and how long the pathogens they might carry survive are all decisive factors for the spread of droplet-borne diseases. The subject is extremely multifaceted and its aspects range across different disciplines, yet most of them have only seldom been considered in the physics community. In this review, we discuss the physical principles that govern the fate of respiratory droplets and any viruses trapped inside them, with a focus on the role of relative humidity. Importantly, low relative humidity-as encountered, for instance, indoors during winter and inside aircraft-facilitates evaporation and keeps even initially large droplets suspended in air as aerosol for extended periods of time. What is more, relative humidity affects the stability of viruses in aerosol through several physical mechanisms such as efflorescence and inactivation at the air-water interface, whose role in virus inactivation nonetheless remains poorly understood. Elucidating the role of relative humidity in the droplet spread of disease would permit us to design preventive measures that could aid in reducing the chance of transmission, particularly in indoor environment.

Ion-Specific Adsorption on Bare Gold (Au) Nanoparticles in Aqueous Solutions: Double-Layer Structure and Surface Potentials.

We study the solvation and electrostatic properties of bare gold (Au) nanoparticles (NPs) of 1-2 nm in size in aqueous electrolyte solutions of sodium salts of various anions with large physicochemical diversity (Cl-, BF4-, PF6-, Nip- (nitrophenolate), 3- and 4-valent hexacyanoferrate (HCF)) using nonpolarizable, classical molecular dynamics computer simulations. We find a substantial facet selectivity in the adsorption structure and spatial distribution of the ions at the AuNPs: while sodium and some of the anions (e.g., Cl-, HCF3-) adsorb more at the "edgy" (100) and (110) facets of the NPs, where the water hydration structure is more disordered, other ions (e.g., BF4-, PF6-, Nip-) prefer to adsorb strongly on the extended and rather flat (111) facets. In particular, Nip-, which features an aromatic ring in its chemical structure, adsorbs strongly and perturbs the first water monolayer structure on the NP (111) facets substantially. Moreover, we calculate adsorptions, radially resolved electrostatic potentials as well as the far-field effective electrostatic surface charges and potentials by mapping the long-range decay of the calculated electrostatic potential distribution onto the standard Debye-Hückel form. We show how the extrapolation of these values to other ionic strengths can be performed by an analytical Adsorption-Grahame relation between the effective surface charge and potential. We find for all salts negative effective surface potentials in the range from -10 mV for NaCl down to about -80 mV for NaNip, consistent with typical experimental ranges for the zeta potential. We discuss how these values depend on the surface definition and compare them to the explicitly calculated electrostatic potentials near the NP surface, which are highly oscillatory in the ±0.5 V range.

Tuning the selective permeability of polydisperse polymer networks.

We study the permeability and selectivity ('permselectivity') of model membranes made of polydisperse polymer networks for molecular penetrant transport, using coarse-grained, implicit-solvent computer simulations. In our work, permeability P is determined on the linear-response level using the solution-diffusion model, P = KDin, i.e., by calculating the equilibrium penetrant partition ratio K and penetrant diffusivity Din inside the membrane. We vary two key parameters, namely the network-network interaction, which controls the degree of swelling and collapse of the network, and the network-penetrant interaction, which tunes the selective penetrant uptake and microscopic energy landscape for diffusive transport. We find that the partitioning K covers four orders of magnitude and is a non-monotonic function of the parameters, well interpreted by a second-order virial expansion of the free energy of transferring one penetrant from a reservoir into the membrane. Moreover, we find that the penetrant diffusivity Din in the polydisperse networks, in contrast to highly ordered membrane structures, exhibits relatively simple exponential decays. We propose a semi-empirical scaling law for the penetrant diffusion that describes the simulation data for a wide range of densities and interaction parameters. The resulting permeability P turns out to follow the qualitative behavior (including maximization and minimization) of partitioning. However, partitioning and diffusion are typically anti-correlated, yielding large quantitative cancellations, controlled and fine-tuned by the network density and interactions, as rationalized by our scaling laws. We finally demonstrate that even small changes of network-penetrant interactions, e.g., by half a kBT, modify the permselectivity by almost one order of magnitude.

Competitive sorption of monovalent and divalent ions by highly charged globular macromolecules.

When a highly charged globular macromolecule, such as a dendritic polyelectrolyte or charged nanogel, is immersed into a physiological electrolyte solution, monovalent and divalent counterions from the solution bind to the macromolecule in a certain ratio and thereby almost completely electroneutralize it. For charged macromolecules in biological media, the number ratio of bound monovalent vs divalent ions is decisive for the desired function. A theoretical prediction of such a sorption ratio is challenging because of the competition of electrostatic (valency), ion-specific, and binding saturation effects. Here, we devise and discuss a few approximate models to predict such an equilibrium sorption ratio by extending and combining established electrostatic binding theories such as Donnan, Langmuir, Manning, and Poisson-Boltzmann approaches, to systematically study the competitive uptake of monovalent and divalent counterions by the macromolecule. We compare and fit our models to coarse-grained (implicit-solvent) computer simulation data of the globular polyelectrolyte dendritic polyglycerol sulfate (dPGS) in salt solutions of mixed valencies. The dPGS molecule has high potential to serve in macromolecular carrier applications in biological systems and at the same time constitutes a good model system for a highly charged macromolecule. We finally use the simulation-informed models to extrapolate and predict electrostatic features such as the effective charge as a function of the divalent ion concentration for a wide range of dPGS generations (sizes).

Tuning the Permeability of Dense Membranes by Shaping Nanoscale Potentials.

Permeability is one of the most fundamental transport properties in soft matter physics, material engineering, and nanofluidics. Here, we report by means of Langevin simulations of ideal penetrants in a nanoscale membrane made of a fixed lattice of attractive interaction sites, how the permeability can be massively tuned, even minimized or maximized, by tailoring the potential energy landscape for the diffusing penetrants, depending on the membrane attraction, topology, and density. Supported by limiting scaling theories we demonstrate that the observed nonmonotonic behavior and the occurrence of extreme values of the permeability is far from trivial and triggered by a strong anticorrelation and substantial (orders of magnitude) cancellation between penetrant partitioning and diffusivity, especially within dense and highly attractive membranes.

Cross-linker effect on solute adsorption in swollen thermoresponsive polymer networks.

The selective solute partitioning within a polymeric network is of key importance to applications in which controlled release or uptake of solutes in a responsive hydrogel is required. In this work we investigate the impact of cross-links on solute adsorption in a swollen polymer network by means of all-atom, explicit-water molecular dynamics simulations. We focus on a representative network subunit consisting of poly(N-isopropylacrylamide) (PNIPAM) and N,N'-methylenebisacrylamide (BIS/MBA) cross-linker types. Our studied system consists of one BIS-linker with four atactic PNIPAM chains attached in a tetrahedral geometry. The adsorption of several representative solutes of different polarity in the low concentration limit at the linker region is examined. We subdivide the solute adsorption regions and distinguish between contributions stemming from polymer chains and cross-link parts. In comparison to a single polymer chain, we observe that the adsorption of the solutes to the cross-link region can significantly differ, with details depending on the specific compounds' size and polarity. In particular, for solutes that have already a relatively large affinity to PNIPAM chains the dense cross-link region (where many-body attractions are at play) amplifies the local adsorption by an order of magnitude. We also find that the cross-link region can serve as a seed for the aggregation of mutually attractive solutes at higher solute concentrations. Utilizing the microscopic adsorption coefficients in a mean-field model of an idealized macroscopic polymer network, we extrapolate these results to the global solute partitioning in a swollen hydrogel and predict that these adsorption features may lead to non-monotonic partition ratios as a function of the cross-link density.

Charge and hydration structure of dendritic polyelectrolytes: molecular simulations of polyglycerol sulphate.

Macromolecules based on dendritic or hyperbranched polyelectrolytes have been emerging as high potential candidates for biomedical applications. Here we study the charge and solvation structure of dendritic polyglycerol sulphate (dPGS) of generations 0 to 3 in aqueous sodium chloride solution by explicit-solvent molecular dynamics computer simulations. We characterize dPGS by calculating several important properties such as relevant dPGS radii, molecular distributions, the solvent accessible surface area, and the partial molecular volume. In particular, as the dPGS exhibits high charge renormalization effects, we address the challenges of how to obtain a well-defined effective charge and surface potential of dPGS for practical applications. We compare implicit- and explicit-solvent approaches in our all-atom simulations with the coarse-grained simulations from our previous work. We find consistent values for the effective electrostatic size (i.e., the location of the effective charge of a Debye-Hückel sphere) within all the approaches, deviating at most by the size of a water molecule. Finally, the excess chemical potential of water insertion into dPGS and its thermodynamic signature are presented and rationalized.

Ionic structure around polarizable metal nanoparticles in aqueous electrolytes.

Metal nanoparticles are receiving increased scientific attention owing to their unique physical and chemical properties that make them suitable for a wide range of applications in diverse fields, such as electrochemistry, biochemistry, and nanomedicine. Their high metallic polarizability is a crucial determinant that defines their electrostatic character in various electrolyte solutions. Here, we introduce a continuum-based model of a metal nanoparticle with explicit polarizability in the presence of different kinds of electrolytes. We employ several, variously sophisticated, theoretical approaches, corroborated by Monte Carlo simulations in order to elucidate the basic electrostatics principles of the model. We investigate how different kinds of asymmetries between the ions result in non-trivial phenomena, such as charge separation and a build-up of a so-called zero surface-charge double layer.

From hydration repulsion to dry adhesion between asymmetric hydrophilic and hydrophobic surfaces.

Using all-atom molecular dynamics (MD) simulations at constant water chemical potential in combination with basic theoretical arguments, we study hydration-induced interactions between two overall charge-neutral yet polar planar surfaces with different wetting properties. Whether the water film between the two surfaces becomes unstable below a threshold separation and cavitation gives rise to long-range attraction, depends on the sum of the two individual surface contact angles. Consequently, cavitation-induced attraction also occurs for a mildly hydrophilic surface interacting with a very hydrophobic surface. If both surfaces are very hydrophilic, hydration repulsion dominates at small separations and direct attractive force contribution can-if strong enough-give rise to wet adhesion in this case. In between the regimes of cavitation-induced attraction and hydration repulsion we find a narrow range of contact angle combinations where the surfaces adhere at contact in the absence of cavitation. This dry adhesion regime is driven by direct surface-surface interactions. We derive simple laws for the cavitation transition as well as for the transition between hydration repulsion and dry adhesion, which favorably compare with simulation results in a generic adhesion state diagram as a function of the two surface contact angles.

Interaction of Charged Patchy Protein Models with Like-Charged Polyelectrolyte Brushes.

We study the adsorption of charged patchy particle models (CPPMs) on a thin film of a like-charged and dense polyelectrolyte (PE) brush (of 50 monomers per chain) by means of implicit-solvent, explicit-salt Langevin dynamics computer simulations. Our previously introduced set of CPPMs embraces well-defined one- and two-patched spherical globules, each of the same net charge and (nanometer) size, with mono- and multipole moments comparable to those of small globular proteins. We focus on electrostatic effects on the adsorption far away from the isoelectric point of typical proteins, i.e., where charge regulation plays no role. Despite the same net charge of the brush and globule, we observe large binding affinities up to tens of the thermal energy, kBT, which are enhanced by decreasing salt concentration and increasing charge of the patch(es). Our analysis of the distance-resolved potentials of mean force together with a phenomenological description of all leading interaction contributions shows that the attraction is strongest at the brush surface, driven by multipolar, Born (self-energy), and counterion-release contributions, dominating locally over the monopolar and steric repulsions.

Interactions between charged particles with bathing multivalent counterions: experiments vs. dressed ion theory.

We compare the recent experimentally measured forces between charged colloidal particles, as well as their effective surface potentials (surface charge) in the presence of multivalent counterions in a bathing monovalent salt solution, with the predictions of the dressed ion theory of strongly charged colloidal systems. The benchmark for comparison is provided by the DLVO theory and the deviations from its predictions at small separations are taken as an indication of the additional non-DLVO attractions that can be fitted by an additional phenomenological exponential term. The parameters characterizing this non-DLVO exponential term as well as the dependencies of the effective potential on the counterion concentration and valency predicted by the dressed ion theory are well within the experimental values. This suggests that the deviations from the DLVO theory are probably caused by ion correlations as formalized within the dressed ion theory.

Selective solute adsorption and partitioning around single PNIPAM chains.

Thermoresponsive polymer architectures have become integral building blocks of 'smart' functional materials in modern applications. For a large range of developments, e.g. for drug delivery or nanocatalytic carrier systems, the selective adsorption and partitioning of molecules (ligands or reactants) inside the polymeric matrix are key processes that have to be controlled and tuned for the desired material function. In order to gain insights into the nanoscale structure and binding details in such systems, we here employ molecular dynamics simulations of the popular poly(N-isopropylacrylamide) (PNIPAM) polymer in explicit water in the presence of various representative solute types with a focus on aromatic model reactants. We study a single polymer chain and explore the influence of its elongation, stereochemistry, and temperature on the solute binding affinities. While we find that the excess adsorption generally increases with the size of the solute, the temperature-dependent affinity to the chain is highly solute specific and has a considerable dependence on the polymer elongation (i.e. polymer swelling state). We elucidate the molecular mechanisms of the selective binding in detail and eventually present how the results can be extrapolated to macroscopic partitioning of the solutes in swollen polymer architectures, such as hydrogels.

Going beyond the standard line tension: Size-dependent contact angles of water nanodroplets.

The dependence of the contact angle on the size of a nanoscopic droplet residing on a flat substrate is traditionally ascribed solely to line tension. Other contributions, stemming from the droplet geometry dependence of the surface tension and line tension, are typically ignored. Here, we perform molecular dynamics simulations of water droplets of cylindrical morphology on surfaces of a wide range of polarities. In the cylindrical geometry, where the line tension is not operative directly, we find significant contact angle dependence on the droplet size. The effect is most pronounced on hydrophilic surfaces, with the contact angle increase of up to 10° with a decreasing droplet size. On hydrophobic surfaces, the trend is reversed and considerably weaker. Our analysis suggests that these effects can be attributed to the Tolman correction due to the curved water-vapor interface and to a generalized line tension that possesses a contact angle dependence. The latter is operative also in the cylindrical geometry and yields a comparable contribution to the contact angle as the line tension itself in case of spherical droplets.