2D materials, characterized by their extensive surface area and customizable chemical and electronic properties, offer compelling advantages as advanced materials. These unique attributes pave the way for the development of next-generation electronics and optoelectronics, photo- and electro-catalysis, energy storage and conversion devices, and sensors. The most prominent and commonly available 2D transition metal dichalcogenide, molybdenum disulfide (MoS₂), has already shown its potential for advanced applications. However, its relatively unfavorable electronic structure and limited intrinsic conductivity lower its suitability for applications that require high conductivity, such as electrocatalysts. One way to enhance its conductivity is by electrochemically intercalating alkali metal ions, e.g., Na⁺ and K⁺, into its layered structure, potentially adjusting its electronic structure. Here, we present a comprehensive investigation into the atomic-scale intercalation mechanism using molecular dynamics simulations, complemented by experimental analysis of structural and electronic properties at the macro scale through various characterization techniques. It is demonstrated that the hydration shell of ions serves as an energy barrier to intercalation as it undergoes a structural change during the intercalation. When alkali metal ions are intercalated into MoS₂, they introduce more defects and enhance conductivity. Notably, these effects are more pronounced for potassium than for sodium.

Nanofluidic membranes (NFMs) are gaining prominence as alternative ion-exchange membranes, because of their distinct selectivity mechanism, which does not rely on functional groups on a polymeric backbone but rather on charged nanopores that allow straight ion-conductive pathways for efficient ion transport. We measured the conductivity of commercial anodized aluminum oxide membranes with different pore sizes under different current densities and electrolyte concentrations. We also simulated a nanopore channel with charged walls between two electrolyte reservoirs. Our findings indicate that electrolyte concentration is the main parameter that determines NFM conductivity, with a linear dependence at least up to 1 M. Our study shows that the optimal pore length is between 0.5 and 5 μm considering the trade-off between selectivity and conductance. On the other hand, the conductance is not sensitive to the pore diameter. Conical nanopores are a way to increase conductance, but according to our results, this increase comes at the expense of selectivity. Our findings suggest that NFMs can outperform polymeric ion-exchange membranes in certain electrochemical applications, such as reverse electrodialysis, but not in applications that use low electrolyte concentrations on both sides of the membrane.

The local conditions inside a gas diffusion electrode (GDE) pore, especially in the electrical double layer (EDL) region, influence the charge transfer reactions and the selectivity of desired CO₂ER products. Most GDE computational models ignore the EDL or are limited in their applicability at high potentials. In this work, we present a continuum model to describe the local environment inside a catalytic pore at varying potentials, electrolyte concentrations and pore diameters. The systems studied in this work are based on an Ag catalyst in contact with KHCO₃ solution. Our study shows that steric effects dominate the local environment at high cathodic potentials (≪−25 mV vs pzc at the OHP), leading to a radial drop of CO₂ concentration. We also observe a drop in pH value within 1 nm of the reaction plane due to electrostatic repulsion and attraction of OH⁻ and H⁺ ions, respectively. We studied the influence of pore radii (1-10 nm) on electric field and concentrations. Pores with a radius smaller than 5 nm show a higher mean potential, which lowers the mean CO₂ concentration. Pores with a favourable local environment can be designed by regulating the ratio between the pore radius and Debye length.

Metal-organic frameworks (MOF) are a class of crystalline, porous materials possessing well-defined channels that have widespread applications across the sustainable landscape. Analogous to zeolites, these materials are well-suited for adsorption processes targeting environmental contaminants. Herein, a zirconium MOF, UiO-66, was functionalized with ferrocene for the selective removal of oxyanion contaminants, specifically NO₃^-, SO₄^2-, and PO₄^3-. Electrochemical oxidation of the embedded ferrocene pendants induces preferential adsorption of these oxyanions, even in the presence of Cl^- in a 10-fold excess. Anion selectivity strongly favoring PO₄^3- (S_oxy/comp = 3.80) was observed following an adsorption trend of PO₄^3- > SO₄^2- > NO₃^- > (10-fold)Cl^- in multi-ion solution mixtures. The underlying mechanisms responsible for ion selectivity were elucidated by performing ex situ X-ray photoelectron spectroscopy (XPS) on the heterogeneous electrode surface postadsorption and by calculating the electronic structure of various adsorption configurations. It was eventually shown that oxyanion selectivity stemmed from strong ion association with a positively charged pore interior due to the spatial distribution of charge by oxygen constituents. While ferrocenium provided the impetus for ion migration-diffusion, it was the formation of stable complexes with zirconium nodes that ultimately contributed to selective adsorption of oxyanions.

Correction to: Scientific Reportshttps://doi.org/10.1038/s41598-024-68971-x, published online 13 August 2024 In the original version of this Article a previous rendition of Figure 2B, Figure 4 and Figure 5D was published. The original Figure 2, 4 and 5 and accompanying legends appear below. (Figure presented.) (Figure presented.) (Figure presented.) (A) Representative bubble dynamics for different channel geometries. (B) Universal motion of bubbles within channels with different size, shape and length. The dashed line represents the developed theory, Eq. (2). The marker colors represent the hydraulic diameters (d_h), the shapes represent the cross-section and the facecolor represent the lengths (L). The graphical marker symbols and colors established here are followed throughout this article. The black arrow represents the region of deviation(s) from the expected dynamics. The threshold laser energy absorbed for bubble formation estimated from experiments (E_th,exp) against theory (E_th,theory) presented in Eq. (5). (A,B) Representative dynamic bubble size curves illustrating the emergence of instabilities. The zones of the instabilities are highlighted using a shaded rectangular area. The arrows represent if the instabilities occur before or after X_max. (A) Illustrates the experimental data for different d_h with similar oscillation time. The instabilities emerge with increasing d_h. (B) Illustrates the data for d_h = 200 µm with increasing laser energies. The instabilities disappear with increasing E_abs. (C) Flow stability diagram with the transition line at W_o = 734. The markers represent the experiments and the lines represent the analytical estimate. The numbers correspond to the channel hydraulic diameters (in µm) with the dashed and solid lines representing the channel lengths L = 25 and 50 mm, respectively. (D) The dimensionless convective timescale against the L/d_h aspect ratio. The partition line is a linear relation between the x and y axes with 45 × 10⁻⁶ as the slope and the origin as the intercept. The original Article has been corrected.

Oscillatory flow in confined spaces is central to understanding physiological flows and rational design of synthetic periodic-actuation based micromachines. Using theory and experiments on oscillating flows generated through a laser-induced cavitation bubble, we associate the dynamic bubble size (fluid velocity) and bubble lifetime to the laser energy supplied—a control parameter in experiments. Employing different channel cross-section shapes, sizes and lengths, we demonstrate the characteristic scales for velocity, time and energy to depend solely on the channel geometry. Contrary to the generally assumed absence of instability in low Reynolds number flows (<1000), we report a momentary flow distortion that originates due to the boundary layer separation near channel walls during flow deceleration. The emergence of distorted laminar states is characterized using two stages. First the conditions for the onset of instabilities is analyzed using the Reynolds number and Womersley number for oscillating flows. Second the growth and the ability of an instability to prevail is analyzed using the convective time scale of the flow. Our findings inform rational design of microsystems leveraging pulsatile flows via cavitation-powered microactuation.

Electrochemical reduction of CO2 heavily depends on the reaction conditions found near the electrode surface. These local conditions are affected by phenomena such as electric double layer formation and steric effects of the solution species, which in turn impact the passage of CO2 molecules to the catalytic surface. Most models for CO2 reduction ignore these effects, leading to an incomplete understanding of the local electrode environment. In this work, we present a modeling approach consisting of a set of size-modified Poisson–Nernst–Planck equations and the Frumkin interpretation of Tafel kinetics. We introduce a modification to the steric effects inside the transport equations which results in more realistic concentration profiles. We also show how the modification lends the model numerical stability without adopting any separate stabilization technique. The model can replicate experimental current densities and faradaic efficiencies till −1.5 vs. SHE/V of applied electrode potential. We also show the utility of this approach for systems operating at elevated CO2 pressures. Using Frumkin-corrected kinetics gels well with the theoretical understanding of the double layer. Hence, this work provides a sound mechanistic understanding of the CO2 reduction process, from which new insights on key performance controlling parameters can be obtained.

Crystallization abounds in nature and industrial practice. A plethora of indispensable products ranging from agrochemicals and pharmaceuticals to battery materials are produced in crystalline form in industrial practice. Yet, our control over the crystallization process across scales, from molecular to macroscopic, is far from complete. This bottleneck not only hinders our ability to engineer the properties of crystalline products essential for maintaining our quality of life but also hampers progress toward a sustainable circular economy in resource recovery. In recent years, approaches leveraging light fields have emerged as promising alternatives to manipulate crystallization. In this review article, we classify laser-induced crystallization approaches where light-material interactions are utilized to influence crystallization phenomena according to proposed underlying mechanisms and experimental setups. We discuss nonphotochemical laser-induced nucleation, high-intensity laser-induced nucleation, laser trapping-induced crystallization, and indirect methods in detail. Throughout the review, we highlight connections among these separately evolving subfields to encourage the interdisciplinary exchange of ideas.

We demonstrate that a cavitation bubble initiated by a Nd:YAG laser pulse below breakdown threshold induces crystallization from supersaturated aqueous solutions with supersaturation and laser-energy-dependent nucleation kinetics. Combining high-speed video microscopy and simulations, we argue that a competition between the dissipation of absorbed laser energy as latent and sensible heat dictates the solvent evaporation rate and creates a momentary supersaturation peak at the vapor-liquid interface. The number and morphology of crystals correlate to the characteristics of the simulated supersaturation peak.

Optofluidic devices have revolutionized the manipulation and transportation of fluid at smaller length scales ranging from micrometers to millimeters. We describe a dedicated optical setup for studying laser-induced cavitation inside a microchannel. In a typical experiment, we use a tightly focused laser beam to locally evaporate the solution laced with a dye resulting in the formation of a microbubble. The evolving bubble interface is tracked using high-speed microscopy and digital image analysis. Furthermore, we extend this system to analyze fluid flow through fluorescence-Particle Image Velocimetry (PIV) technique with minimal adaptations. In addition, we demonstrate the protocols for the in-house fabrication of a microchannel tailored to function as a sample holder in this optical setup. In essence, we present a complete guide for constructing a fluorescence microscope from scratch using standard optical components with flexibility in the design and at a lower cost compared to its commercial analogues.

One of the important parameters in water management of proton exchange membranes is the electro-osmotic drag (EOD) coefficient of water. The value of the EOD coefficient is difficult to justify, and available literature data on this for Nafion membranes show scattering from in experiments and simulations. Here, we use a classical all-atom model to compute the EOD coefficient and thermodynamic properties of water from molecular dynamics simulations for temperatures between 330 and 420 K, and for different water contents between λ = 5 and λ = 20. λ is the ratio between the moles of water molecules to the moles of sulfonic acid sites. This classical model does not capture the Grotthuss mechanism; however, it is shown that it can predict the EOD coefficient within the range of experimental values for λ = 5 where the vehicular mechanism dominates proton transfer. For λ > 5, the Grotthuss mechanism becomes dominant. To obtain the EOD coefficient, average velocities of water and ions are computed by imposing different electric fields to the system. Our results show that the velocities of water and hydronium scale linearly with the electric field, resulting in a constant ratio of ca. 0.4 within the error bars. We find that the EOD coefficient of water linearly increases from 2 at λ = 5 to 8 at λ = 20 and the results are not sensitive to temperature. The EOD coefficient at λ = 5 is within the range of experimental values, confirming that the model can capture the vehicular transport of protons well. At λ = 20, due to the absence of proton hopping in the model, the EOD coefficient is overestimated by a factor of 3 compared to experimental values. To analyze the interactions between water and Nafion, the partial molar enthalpies and partial molar volumes of water are computed from molecular dynamics simulations. At different water uptakes, multiple linear regression is used on raw simulation data within a narrow composition range of water inside the Nafion membrane. The partial molar volumes and partial molar excess enthalpies of water asymptotically approach the molar volumes and molar excess enthalpies of pure water for water uptakes above 5. This confirms the model can capture the bulklike behavior of water in the Nafion at high water uptakes.

Surface conductivity in the electrical double layer (EDL) is known to be affected by proton hopping and diffusion at solid-liquid interfaces. Yet, the role of surface protolysis and its kinetics on the thermodynamic and transport properties of the EDL are usually ignored as physical models consider static surfaces. Here, using a novel molecular dynamics method mimicking surface protolysis, we unveil the impact of such chemical events on the system’s response. Protolysis is found to strongly affect the EDL and electrokinetic aspects with major changes in Formula Presented potential and electro-osmotic flow.

In heterogeneous catalysis, reactivity and selectivity are not only influenced by chemical processes occurring on catalytic surfaces but also by physical transport phenomena in the bulk fluid and fluid near the reactive surfaces. Because these processes take place at a large range of time and length scales, it is a challenge to model catalytic reactors, especially when dealing with complex surface reactions that cannot be reduced to simple mean-field boundary conditions. As a particle-based mesoscale method, Stochastic Rotation Dynamics (SRD) is well suited for studying problems that include both microscale effects on surfaces and transport phenomena in fluids. In this work, we demonstrate how to simulate heterogeneous catalytic reactors by coupling an SRD fluid with a catalytic surface on which complex surface reactions are explicitly modeled. We provide a theoretical background for modeling different stages of heterogeneous surface reactions. After validating the simulation method for surface reactions with mean-field assumptions, we apply the method to non-mean-field reactions in which surface species interact with each other through a Monte Carlo scheme, leading to island formation on the catalytic surface. We show the potential of the method by simulating a more complex three-step reaction mechanism with reactant dissociation.

Biological nanopores such as aquaporins combine the opposing functions of high water permeation and total ion exclusion in part by the virtue of their hourglass shape. Here, we perform molecular dynamics simulations to examine water and ion conduction through hourglass shaped nanopores created from carbon nanotubes (CNTs) of chirality (6,6), (8,8), and (10,10) in combination with carbon nanocones of half cone angles 41.8°, 30.0°, 19.45°, 9.6° and 0.0°. We observe large variations in flow through the nanopores with change in half cone angles and tube diameters. By computing the pore-water interactions we find a correlated change between the flux and the density profiles of water inside the nanopores. Further, from the orientation, and the hydrogen bonding characteristics of water, we uncover some unexplored facets of flow through hourglass shaped nanopores. The results are insightful for devising novel separation membranes based on nanopores that mimic the shape of biological nanochannels.

Many physical and chemical processes involve energy change with rates that depend sensitively on local temperature. Important examples include heterogeneously catalyzed reactions and activated desorption. Because of the multiscale nature of such systems, it is desirable to connect the macroscopic world of continuous hydrodynamic and temperature fields to mesoscopic particle-based simulations with discrete particle events. In this work we show how to achieve real-time measurement of the local temperature in stochastic rotation dynamics (SRD), a mesoscale method particularly well suited for problems involving hydrodynamic flows with thermal fluctuations. We employ ensemble averaging to achieve local temperature measurement in dynamically changing environments. After validation by heat diffusion between two isothermal plates, heating of walls by a hot strip, and by temperature programed desorption, we apply the method to a case of a model flow reactor with temperature-sensitive heterogeneously catalyzed reactions on solid spherical catalysts. In this model, adsorption, chemical reactions, and desorption are explicitly tracked on the catalyst surface. This work opens the door for future projects where SRD is used to couple hydrodynamic flows and thermal fluctuations to solids with complex temperature-dependent surface mechanisms.

Molecular Dynamics (MD) simulations are uniquely suitable for providing molecular-level insights into the Electric Double Layer (EDL) that forms when a charged surface is in contact with an aqueous solution. However, simulations are only as accurate in predicting EDL properties as permitted by the atomic interaction models. Experimental ζ-potential values and surface charges could provide a potentially suitable reference to validate and tune the interaction models, if not for the fact that they themselves are a product of imperfect models used to interpret the raw measurement data. Here, we present an approach to tune an interaction model by comparing Electro-Osmotic Flow (EOF) MD simulations against experimental Streaming Current (SC) measurements while minimizing potential modeling errors arising from both approaches. The point that is least susceptible to interpretation and modeling errors is argued to be at the concentration for which zero flow velocity is observed in EOF simulations and a net zero electric current is measured in SC experiments. At this concentration, the ζ-potential is also zero. We were able to match the experimental concentration at which ζ = 0 in MD simulations for a CaCl2 solution at pH 7.5 in contact with fused silica by tuning the ion-surface Lennard-Jones cross interactions. These interactions were found to greatly affect the ion distribution within the EDL and particularly the formation of inner-sphere surface-complexes, which, in turn, affects the electrokinetic flow. With the ion distribution determined explicitly, a series of properties can be calculated unambiguously, such as the capacitance needed for surface complexation models.

Surfactant molecules, known as organic friction modifiers (OFMs), are routinely added to lubricants to reduce friction and wear between sliding surfaces. In macroscale experiments, friction generally decreases as the coverage of OFM molecules on the sliding surfaces increases; however, recent nanoscale experiments with sharp atomic force microscopy (AFM) tips have shown increasing friction. To elucidate the origin of these opposite trends, we use nonequilibrium molecular dynamics (NEMD) simulations and study kinetic friction between OFM monolayers and an indenting nanoscale asperity. For this purpose, we investigate various coverages of stearamide OFMs on iron oxide surfaces and silica AFM tips with different radii of curvature. We show that the differences between the friction-coverage relations from macroscale and nanoscale experiments are due to molecular plowing in the latter. For our small tip radii, the friction coefficient and indentation depth both have a nonmonotonic dependence on OFM surface coverage, with maxima occurring at intermediate coverage. We rationalize the nonmonotonic relations through a competition of two effects (confinement and packing density) that varying the surface coverage has on the effective stiffness of the OFM monolayers. We also show that kinetic friction is not very sensitive to the sliding velocity in the range studied, indicating that it originates from instabilities. Indeed, we find that friction predominately originates from plowing of the monolayers by the leading edge of the tip, where gauche defects are created, while thermal dissipation is mostly localized in molecules toward the trailing edge of the tip, where the chains return to a more extended conformation.

This study aims to elucidate the impact of salinity on the interactions governing the adsorption of polar aromatic oil compounds onto calcite. To this end, molecular dynamics simulations were employed to assess adsorption of a model polar organic molecule (deprotonated benzoic acid, benzoate) on the calcite surface in NaCl brines of different concentration levels, namely, deionized water (DW), low-salinity water (LS, 5000 ppm), and sea water (SW; 45,000 ppm). Calcite was found to be completely covered by several well-ordered water layers. The top hydration layer is very compact and prevents direct adsorption of benzoates onto the substrate. Instead, Na+ ions form a distinct positively charged layer by adhering on the calcite substrate through inner-sphere complexion mode. Cl– ions mostly lodge on top of the adsorbed sodium cations, forming a negatively charged layer. The distribution of ions at the calcite/brine interface thus exhibits the features of an electrical double layer, composed of a Stern-like positive layer followed by a negative one. The positive charged layer attracts benzoates toward the surface. As such, the sodium ions attached onto the calcite can act as adsorption sites to connect benzoates to the surface. By increasing the salinity, more Na+ ions adsorb onto the calcite surface, and the density of benzoate molecules at the interface is enhanced as a result of more Na+ bridging ions. The monotonic salinity-dependent adsorption of benzoate molecules on calcite offers a mechanism driving additional oil recovery upon injection of diluted brine into subsurface carbonate reservoirs.

The intriguing mass transport properties of carbon nanotubes (CNTs) have received widespread attention, especially the rapid transport of water through CNTs due to their atomically smooth wall interiors. Extensive research has been dedicated to the comprehension of various aspects of water flow in contact with CNTs, the most prominent ones being the studies on slip and flow rates. Experimental and computational studies have confirmed an enhanced water flow rate through this graphitic nanoconfinement. However, a quantitative agreement has not yet been attained. These disparities coupled with incomplete knowledge of the mechanisms of water transport at nanoscale regimes are hindering the possibilities to integrate CNTs in numerous nanofluidic applications. In the present review, we focus on the slip and flow rates of water through CNTs and the factors influencing them. We discuss the key sources of discrepancies in water flow rate and suggest directions for future study.

Recently, β-cyclodextrin (βCD)-based polymers with enhanced adsorption kinetics and high removal capacity of organic micropollutants (OMPs) and uptake rates have been synthesized and tested experimentally. Although the exact physical-chemical mechanisms via which these polymers capture the various types of OMPs are not yet fully understood, it is suggested that the inclusion complex formation of OMPs with βCD is very important. In this study, the inclusion complex formation of OMPs with βCD in an aqueous solution is investigated by using the well-established attach-pull-release method in force field-based molecular dynamics simulations. A representative set of OMPs is selected based on the measured occurrences in surface and ground waters and the directives published by the European Union. To characterize the formation of the inclusion complex, the binding free energies, enthalpies, and entropies are computed and compared to experimental values. It is shown that computations using the q4md-CD/GAFF/Bind3P force field combination yield binding free energies that are in reasonable agreement with the experimental results for all OMPs studied. The binding enthalpies are decomposed into the main contributing interaction types. It is shown that, for all studied OMPs, the van der Waals interactions are favorable for the inclusion complexion and the hydrogen bond formation of the guest with the solvent and βCD plays a crucial role in the binding mechanism. Our findings show that MD simulations can adequately describe the inclusion complex formation of βCD with OMPs, which is the first step toward understanding the underlying mechanisms via which the βCD-based polymers capture OMPs.