ABSTRACT
The efficient transport of small molecules through dense hydrogel networks is crucial for various applications, including drug delivery, biosensing, catalysis, nanofiltration, water purification, and desalination. In dense polymer matrices, such as collapsed microgels, molecular transport follows the solution-diffusion principle: Molecules dissolve in the polymeric matrix and subsequently diffuse due to a concentration gradient. Employing dynamical density functional theory (DDFT), we investigate the nonequilibrium release kinetics of nonionic subnanometer-sized molecules from a microgel particle, using parameters derived from prior molecular simulations of a thermoresponsive hydrogel. The kinetics is primarily governed by the microgel radius and two intensive parameters: the diffusion coefficient and solvation free energy of the molecule. Our results reveal two limiting regimes: a diffusion-limited regime for large, slowly diffusing, and poorly soluble molecules within the hydrogel; and a reaction-limited regime for small, rapidly diffusing, and highly soluble molecules. These principles allow us to derive an analytical equation for release time, demonstrating excellent quantitative agreement with the DDFT results-a valuable and straightforward tool for predicting release kinetics from microgels.
ABSTRACT
Adsorption of surfactants to fluid interfaces occurs in numerous technological and daily-life contexts. The coverage at the interface and other properties of the formed adsorption layers determine the performance of a surfactant with regard to the desired application. Given the importance of these applications, there is a great demand for the comprehensive characterization and understanding of surfactant adsorption layers. In this review, we provide an overview of suitable experimental and simulation-based techniques and review the literature in which they were used for the investigation of surfactant adsorption layers. We come to the conclusion that, while these techniques have been successfully applied to investigate Langmuir monolayers of water-insoluble surfactants, their application to the study of Gibbs adsorption layers of water-soluble surfactants has not been fully exploited. Finally, we emphasize the great potential of these methods in providing a deeper understanding of the behavior of soluble surfactants at interfaces, which is crucial for optimizing their performance in various applications.
ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
The borohydride ion, BH4-, is an essential reducing agent in many technological processes, yet its full understanding has been elusive, because of at least two significant challenges. One challenge arises from its marginal stability in aqueous solutions outside of basic pH conditions, which considerably limits the experimental thermodynamic data. The other challenge comes from its unique and atypical hydration shell, stemming from the negative excess charge on its hydrogen atoms, which complicates the accurate modeling in classical atomistic simulations. In this study, we combine experimental and computer simulation techniques to devise a classical force field for NaBH4 and deepen our understanding of its characteristics. We report the first measurement of the ion's activity coefficient and extrapolate it to neutral pH conditions. Given the difficulties in directly measuring its solvation free energies, owing to its instability, we resort to quantum chemistry calculations. This combined strategy allows us to derive a set of nonpolarizable force-field parameters for the borohydride ion for classical molecular dynamics simulations. The derived force field simultaneously captures the solvation free energy, the hydration structure, as well as the activity coefficient of NaBH4 salt across a broad concentration range. The obtained insights into the hydration shell of the BH4- ion are crucial for accurately modeling and understanding its interactions with other molecules, ions, materials, and interfaces.
ABSTRACT
Surfactants play essential roles in many commonplace applications and industrial processes. Although significant progress has been made over the past decades with regard to model-based predictions of the behavior of surfactants, important challenges have remained. Notably, the characteristic time scales of surfactant exchange among micelles, interfaces, and the bulk solution typically exceed the time scales currently accessible with atomistic molecular dynamics (MD) simulations. Here, we circumvent this problem by introducing a framework that combines the general thermodynamic principles of self-assembly and interfacial adsorption with atomistic MD simulations. This approach provides a full thermodynamic description based on equal chemical potentials and connects the surfactant bulk concentration, the experimental control parameter, with the surfactant surface density, the suitable control parameter in MD simulations. Self-consistency is demonstrated for the nonionic surfactant C12EO6 (hexaethylene glycol monododecyl ether) at an alkane/water interface, for which the adsorption and pressure isotherms are computed. The agreement between the simulation results and experiments is semiquantitative. A detailed analysis reveals that the used atomistic model captures well the interactions between surfactants at the interface but less so their adsorption affinities to the interface and incorporation into micelles. Based on a comparison with other recent studies that pursued similar modeling challenges, we conclude that the current atomistic models systematically overestimate the surfactant affinities to aggregates, which calls for improved models in the future.
ABSTRACT
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.
ABSTRACT
Lipid monolayers are ubiquitous in biological systems and have multiple roles in biotechnological applications, such as lipid coatings that enhance colloidal stability or prevent surface fouling. Despite the great technological importance of surface-adsorbed lipid monolayers, the connection between their formation and the chemical characteristics of the underlying surfaces has remained poorly understood. Here, we elucidate the conditions required for stable lipid monolayers nonspecifically adsorbed on solid surfaces in aqueous solutions and water/alcohol mixtures. We use a framework that combines the general thermodynamic principles of monolayer adsorption with fully atomistic molecular dynamics simulations. We find that, very universally, the chief descriptor of adsorption free energy is the wetting contact angle of the solvent on the surface. It turns out that monolayers can form and remain thermodynamically stable only on substrates with contact angles above the adsorption contact angle, θads. Our analysis establishes that θads falls into a narrow range of around 60∘-70∘ in aqueous media and is only weakly dependent on the surface chemistry. Moreover, to a good approximation, θads is roughly determined by the ratio between the surface tensions of hydrocarbons and the solvent. Adding small amounts of alcohol to the aqueous medium lowers θads and thereby facilitates monolayer formation on hydrophilic solid surfaces. At the same time, alcohol addition weakens the adsorption strength on hydrophobic surfaces and results in a slowdown of the adsorption kinetics, which can be useful for the preparation of defect-free monolayers.
ABSTRACT
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.
ABSTRACT
Adsorption of small amphiphilic molecules occurs in various biological and technological processes, sometimes desired while other times unwanted (e.g., contamination). Surface-active molecules preferentially bind to interfaces and affect their wetting properties. We use molecular dynamics simulations to study the adsorption of short-chained alcohols (simple surfactants) to the water-vapor interface and solid surfaces of various polarities. With a theoretical analysis, we derive an equation for the adsorption coefficient, which scales exponentially with the molecular surface area and the surface wetting coefficient and is in good agreement with the simulation results. We apply the outcomes to aqueous sessile droplets containing surfactants, where the competition of surfactant adsorptions to both interfaces alters the contact angle in a nontrivial way. The influence of surfactants is the strongest on very hydrophilic and hydrophobic surfaces, whereas droplets on moderately hydrophilic surfaces are less affected.
Subject(s)
Surface-Active Agents , Water , Adsorption , Hydrophobic and Hydrophilic Interactions , Surface-Active Agents/chemistry , Water/chemistry , WettabilityABSTRACT
RNA is a functionally rich molecule with multilevel, hierarchical structures whose role in the adsorption to molecular substrates is only beginning to be elucidated. Here, we introduce a multiscale simulation approach that combines a tractable coarse-grained RNA structural model with an interaction potential of a structureless flat adsorbing substrate. Within this approach, we study the specific role of stem-hairpin and multibranch RNA secondary structure motifs on its adsorption phenomenology. Our findings identify a dual regime of adsorption for short RNA fragments with and without the secondary structure and underline the adsorption efficiency in both cases as a function of the surface interaction strength. The observed behavior results from an interplay between the number of contacts formed at the surface and the conformational entropy of the RNA molecule. The adsorption phenomenology of RNA seems to persist also for much longer RNAs as qualitatively observed by comparing the trends of our simulations with a theoretical approach based on an ideal semiflexible polymer chain.
ABSTRACT
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.
ABSTRACT
The performance of gold nanoparticles (NPs) in applications depends critically on the structure of the NP-solvent interface, at which the electrostatic surface polarization is one of the key characteristics that affects hydration, ionic adsorption, and electrochemical reactions. Here, we demonstrate significant effects of explicit metal polarizability on the solvation and electrostatic properties of bare gold NPs in aqueous electrolyte solutions of sodium salts of various anions (Cl-, BF4-, PF6-, nitrophenolate, and 3- and 4-valent hexacyanoferrate), using classical molecular dynamics simulations with a polarizable core-shell model for the gold atoms. We find considerable spatial heterogeneity of the polarization and electrostatic potentials on the NP surface, mediated by a highly facet-dependent structuring of the interfacial water molecules. Moreover, ion-specific, facet-dependent ion adsorption leads to considerable alterations of the interfacial polarization. Compared to nonpolarizable NPs, surface polarization modifies water local dipole densities only slightly but has substantial effects on the electrostatic surface potentials and leads to significant lateral redistributions of ions on the NP surface. Besides, interfacial polarization effects cancel out in the far field for monovalent ions but not for polyvalent ions, as anticipated from continuum "image-charge" concepts. Far-field effective Debye-Hückel surface potentials change accordingly in a valence-specific fashion. Hence, the explicit charge response of metal NPs is crucial for the accurate description and interpretation of interfacial electrostatics (e.g., for charge transfer and interfacial polarization in catalysis and electrochemistry).
ABSTRACT
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.
ABSTRACT
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.
Subject(s)
Humidity , AerosolsABSTRACT
HYPOTHESIS: In previous studies we looked at the foam stability of various surfactants with C12 alkyl chains but different head groups and found that stable foams are only generated if the head groups are capable of forming hydrogen bonds with each other. Despite the consistency of the experimental data with the conclusions drawn from it we had no direct proof for our hypothesis that H-bonds are formed between surfactant head groups. EXPERIMENTS: To fill this gap, i.e. to demonstrate intersurfactant H-bond formation, we chose the non-ionic sugar surfactant n-dodecyl-ß-d-maltoside (ß-C12G2) and used molecular dynamics (MD) simulations as well as grazing-incidence X-ray (GIX) scattering and diffraction to study the surfactant-loaded air-water interface. FINDINGS: (1) In a densely packed monolayer, close to the critical micelle concentration (cmc), each head group of the sugar surfactant is involved in â¼5 intersurfactant H-bonds with other head groups and in â¼5 H-bonds with water molecules. (2) The number of intersurfactant H-bonds decreases, while the number of surfactant-water H-bonds increases with increasing distance between the head groups (below the cmc). (3) Even at very large distances (well below the cmc) there are still intersurfactant H-bonds, which we ascribe to the formation of clusters at the surface. (4) GIX scattering revealed that a homogeneous surfactant monolayer is formed at full coverage (around the cmc), i.e. cluster formation only happens below the cmc.
ABSTRACT
The permeability of hydrogels for the selective transport of molecular penetrants (drugs, toxins, reactants, etc.) is a central property in the design of soft functional materials, for instance in biomedical, pharmaceutical, and nanocatalysis applications. However, the permeation of dense and hydrated polymer membranes is a complex multifaceted molecular-level phenomenon, and our understanding of the underlying physicochemical principles is still very limited. Here, we uncover the molecular principles of permeability and selectivity in hydrogel permeation. We combine the solution-diffusion model for permeability with comprehensive atomistic simulations of molecules of various shapes and polarities in a responsive hydrogel in different hydration states. We find in particular that dense collapsed states are extremely selective, owing to a delicate balance between the partitioning and diffusivity of the penetrants. These properties are sensitively tuned by the penetrant size, shape, and chemistry, leading to vast cancellation effects, which nontrivially contribute to the permeability. The gained insights enable us to formulate semiempirical rules to quantify and extrapolate the permeability categorized by classes of molecules. They can be used as approximate guiding ("rule-of-thumb") principles to optimize penetrant or membrane physicochemical properties for a desired permeability and membrane functionality.
Subject(s)
Hydrogels , Polymers , Diffusion , PermeabilityABSTRACT
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.
ABSTRACT
We study the nonequilibrium diffusive release of electroneutral molecular cargo encapsulated inside hollow hydrogel nanoparticles. We propose a theoretical model that includes osmotic, steric, and short-range polymer-cargo attractions to determine the effective cargo-hydrogel interaction, ueff*, and the effective diffusion coefficient of the cargo inside the polymer network, Deff*. Using dynamical density functional theory (DDFT), we investigate the scaling of the characteristic release time, τ1/2, with the key parameters involved in the process, namely, ueff*, Deff*, and the swelling ratio. This effort represents a full study of the problem, covering a broad range of cargo sizes and providing predictions for repulsive and attractive polymer shells. Our calculations show that the release time through repulsive polymer networks scales with q2eßueff*/Deff* for ßueff* â« 1. In this case, the cargo molecules are excluded from the shell of the hydrogel. For attractive shells, the polymer retains the cargo molecules on its internal surface and its interior, and the release time grows exponentially with the attraction strength. The DDFT calculations are compared to an analytical model for the mean first passage time, which provides an excellent quantitative description of the kinetics for both repulsive and attractive shells without fitting parameters. Finally, we apply the method to reproduce experimental results on the release of paclitaxel from hollow poly(4-vinylpyridine) nanoparticles and find that the slow release of the drug can be explained in terms of the strong binding attraction between the drug and the polymer.
ABSTRACT
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.