The advancement of microbial electrosynthesis systems (MES) towards industrialization is currently hindered by a limited understanding of the fundamental constraints affecting selective production of high-value chemicals. To address this challenge, we develop a comprehensive computational model that integrates microbial, electrochemical, and acid–base reactions with pore-scale transport processes within a three-dimensionally resolved biofilm. This study investigates the H₂-mediated CO₂ fixation pathway to acetate, butyrate, and caproate. The effect of applied cathode potential and biofilm thickness on macroscopic parameters, such as efficiency and selectivity, is analyzed based on local concentrations and electrochemical and biochemical fluxes. Among the limiting factors, the availability of CO₂ emerges as the main limitation for biochemical reactions due to its low solubility and high half-saturation constant. Additionally, hydrogen – serving as the electron mediator – limits the reaction rate at low current densities and reduces electron transfer efficiency at higher current densities. A key insight from our study is the identification of an optimal electrode potential for each biofilm thickness, balancing both H₂ transfer and CO₂ consumption efficiencies. Furthermore, carbon selectivity shifts with increasing biofilm thickness: net acetate production declines while caproate production increases. This trend is attributed to the prolonged residence time of metabolic intermediates within thicker biofilms, promoting chain elongation pathways. Thus, our work takes an important step towards a fundamental understanding of caproate selectivity across different biofilms, which can be used to optimize the electrode structure and operating conditions to control the local biofilm thickness.

Ultrasound-responsive liposomes represent a promising strategy for the targeted delivery of therapeutic agents to deep tissues, combining the clinically validated biocompatibility of liposomes with precise spatiotemporal control via ultrasound. While the incorporation of polyethylene glycol-conjugated lipids into liposomal formulations has long been known to enhance ultrasound-triggered release, the mechanistic basis of this enhancement remains poorly understood, hindering the rational design of formulations with enhanced ultrasound-triggered release. To address this, we investigated the proposed hypotheses in the literature, including the influence of the packing parameter, differences in thermal, mechanical, and fluidity properties, as well as structural changes such as micelle formation or bilayer thinning for stealth liposomes containing polyethylene glycol, polycarboxybetaine, or polysarcosine. Our results indicate that PEG-enhanced ultrasound-triggered release cannot be attributed to changes in bulk bilayer properties, including the packing parameter, membrane mechanics, thermal behavior, or fluidity. Moreover, the membrane structure remained unchanged after insonation, with no evidence of micelle ejection or membrane thinning. Instead, experiments and computational simulations suggest that polyethylene glycol-conjugated lipopolymers facilitate sonoporation by rearranging into micelle-like structures at the periphery of ultrasound-induced pores in the lipid bilayer, enabled by their low critical micelle concentration, thereby enhancing the release through these transient openings. Finally, we discuss how these insights could be extended beyond polyethylene glycol lipopolymers to guide the design of liposomes with enhanced ultrasound-triggered release, by shifting optimization from bulk bilayer properties to a sonoporation-oriented framework that integrates computational screening, AI-enabled formulation design, and ultrasound waveform considerations.

Molybdenum disulfide (MoS ₂) has emerged as a promising electrocatalyst for the electrochemical reduction of CO ₂, primarily yielding carbon monoxide. However, product selectivity is known to be highly sensitive to structural features such as edge termination and defect density. In this work, we report the formation of higher hydrocarbons (C ₂+ products) enabled by the presence of inherent sulfur vacancies in MoS ₂ when combined with various ionic liquids as co-catalysts. While MoS ₂ has traditionally shown limited hydrocarbon output, our findings demonstrate for the first time that native defect sites, interacting synergistically with the electrolyte environment, can facilitate the production of significant amounts of C ₂+ species. These results provide new insights into defect-mediated catalytic pathways and highlight the importance of electrolyte design in tuning product distribution during CO ₂ electroreduction.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.