Researcher profile

Joachim Dzubiella

Joachim Dzubiella contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

Nonequilibrium free energy during polymer chain growth

During fast diffusion-influenced polymerization, nonequilibrium behavior of the polymer chains and the surrounding reactive monomers has been reported recently. Based on the laws of thermodynamics, the emerging nonequilibrium structures should be characterisable by some "extra free energy" (excess over the equilibrium Helmholtz free energy). Here, we study the nonequilibrium thermodynamics of chain-growth polymerization of ideal chains in a dispersion of free reactive monomers, using off-lattice, reactive Brownian Dynamics (R-BD) computer simulations in conjunction with approximative statistical mechanics and relative entropy (Gibbs-Shannon and Kullback-Leibler) concepts. In case of fast growing polymers, we indeed report increased nonequilibrium free energies of several kBT compared to equilibrium and near-equilibrium, slowly growing chains. Interestingly, the extra free energy is a non-monotonic function of the degree of polymerization and thus also of time. Our decomposition of the thermodynamic contributions shows that the initial dominant extra free energy is stored in the nonequilibrium inhomogeneous density profiles of the free monomer gas (showing density depletion and wakes) in the vicinity of the active center at the propagating polymer end. At later stages of the polymerization process, we report significant extra contributions stored in the nonequilibrium polymer conformations. Finally, our study implies a nontrivial relaxation kinetics and restoring of the extra free energy during the equilibration process after polymerization.

preprint2021arXiv

Active interaction switching controls the dynamic heterogeneity of soft colloidal dispersions

We employ Reactive Dynamical Density Functional Theory (R-DDFT) and Reactive Brownian Dynamics (R-BD) simulations to investigate the dynamics of a suspension of active soft Gaussian colloids with binary interaction switching, i.e., a one-component colloidal system in which every particle stochastically switches at predefined rates between two interaction states with different mobility. Using R-DDFT we extend a theory previously developed to access the dynamics of inhomogeneous liquids to study the influence of the switching activity on the self and distinct part of the Van Hove function in bulk solution, and determine the corresponding mean squared displacement of the switching particles. Our results demonstrate that, even though the average diffusion coefficient is not affected by the switching activity, it significantly modifies the non-equilibrium dynamics and diffusion coefficients of the individual particles, leading to a crossover from short to long times, with a regime for intermediate times showing anomalous diffusion. In addition, the self-part of the van Hove function has a Gaussian form at short and long times, but becomes non-Gaussian at intermediates ones, having a crossover between short and large displacements. The corresponding self-intermediate scattering function shows the two-step relaxation patters typically observed in soft materials with heterogeneous dynamics such as glasses and gels. We also introduce a phenomenological Continuous Time Random Walk (CTRW) theory to understand the heterogeneous diffusion of this system. R-DDFT results are in excellent agreement with R-BD simulations and the analytical predictions of CTRW theory, thus confirming that R-DDFT constitutes a powerful method to investigate not only the structure and phase behavior, but also the dynamical properties of non-equilibrium active switching colloidal suspensions.

preprint2021arXiv

Nanochannels and nanodroplets in polymer membranes controlling ionic transport

Polymer materials with low water uptake exhibit a highly heterogeneous interior, characterized by water clusters in the form of nanodroplets and nanochannels. Here, based on our recent insights from computer simulations, we argue that water cluster structure has large implications for ionic transport and selective permeability in polymer membranes. Importantly, we demonstrate that the two key quantities for transport, the ion diffusion and the solvation free energy inside the polymer, are extremely sensitive to molecular details of the water clusters. In particular, we highlight the significance of water droplet interface potentials and the nature of hopping diffusion through transient water channels. These mechanisms can be harvested and fine-tuned to optimize selectivity in ionic transport in a wide range of applications.

preprint2021arXiv

Transport in polymer membranes beyond linear response: Controlling permselectivity by the driving force

In the popular solution-diffusion picture, the membrane permeability is defined as the product of the partition ratio and the diffusivity of penetrating solutes inside the membrane in the linear response regime, i.e., in equilibrium. However, of practical importance is the penetrants' flux across the membrane driven by external forces. Here, we study nonequilibrium membrane permeation orchestrated by a uniform external driving field using molecular computer simulations and continuum (Smoluchowski) theory in the stationary state. In the simulations, we explicitly resolve the penetrants' transport across a finite monomer-resolved polymer network, addressing one-component penetrant systems and mixtures. We introduce and discuss possible definitions of nonequilibrium, force-dependent permeability, representing `system' and `membrane' permeability. In particular, we present for the first time a definition of the differential permeability response to the force. We demonstrate that the latter turns out to be significantly nonlinear for low-permeable systems, leading to a high amount of selectiveness in permeability, called `permselectivity', and is tunable by the driving force. Our continuum-level analytical solutions exhibit remarkable qualitative agreement with the penetrant- and polymer-resolved simulations, thereby allowing us to characterize the underlying mechanism of permeabilities and steady-state transport beyond the linear response level.

preprint2020arXiv

Affinity, Kinetics, and Pathways of Anisotropic Ligands Binding to Hydrophobic Model Pockets

Using explicit-water molecular dynamics (MD) simulations of a generic pocket-ligand model we investigate how chemical and shape anisotropy of small ligands influences the affinities, kinetic rates and pathways for their association to hydrophobic binding sites. In particular, we investigate aromatic compounds, all of similar molecular size, but distinct by various hydrophilic or hydrophobic residues. We demonstrate that the most hydrophobic sections are in general desolvated primarily upon binding to the cavity, suggesting that specific hydration of the different chemical units can steer the orientation pathways via a `hydrophobic torque'. Moreover, we find that ligands with bimodal orientation fluctuations have significantly increased kinetic barriers for binding compared to the kinetic barriers previously observed for spherical ligands due to translational fluctuations. We exemplify that these kinetic barriers, which are ligand specific, impact both binding and unbinding times for which we observe considerable differences between our studied ligands.

preprint2020arXiv

Aqueous Nanoclusters Govern Ion Partitioning in Dense Polymer Membranes

The uptake and sorption of charged molecules by responsive polymer membranes and hydrogels in aqueous solutions is of key importance for the development of soft functional materials. Here we investigate the partitioning of simple monoatomic (Na$^+$, K$^+$, Cs$^+$, Cl$^-$, I$^-$) and one molecular ion (4-nitrophenolate; NP$^-$) within a dense, electroneutral poly($N$-isopropylacrylamide) membrane using explicit-water molecular dynamics simulations. Inside the predominantly hydrophobic environment water distributes in a network of polydisperse water nanoclusters. The average cluster size determines the mean electrostatic self-energy of the simple ions, which preferably reside deeply inside them; we therefore find substantially larger partition ratios $K\simeq\>$10$^{-1}$ than expected from a simple Born picture using a uniform dielectric constant. Despite their irregular shapes we observe that the water clusters possess a universal negative electrostatic potential with respect to their surrounding, as is known for aqueous liquid-vapor interfaces. This potential, which we find concealed in cases of symmetric monoatomic salts, can dramatically impact the transfer free energies of larger charged molecules because of their weak hydration and increased affinity to interfaces. Consequently, and in stark contrast to the simple ions, the molecular ion NP$^-$ can have a partition ratio much larger than unity, $K\simeq\>$10-30 (depending on the cation type) or even $10^3$ in excess of monovalent salt, which explains recent observations of enhanced reaction kinetics of NP$^-$ reduction catalyzed within dense polymer networks. These results also suggest that ionizing a molecule can even enhance the partitioning in a collapsed, rather hydrophobic gel, which strongly challenges the traditional simplistic reasoning.

preprint2020arXiv

Competitive sorption of mono- versus di-valent 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 mono- versus di-valent 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 as well as Poisson-Boltzmann approaches, to systematically study the competitive uptake of mono- and di-valent 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 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).

preprint2020arXiv

Controlling the microstructure and phase behavior of confined soft colloids by active interaction switching

We explore the microstructure and phase behavior of confined soft colloids which can actively switch their interactions at a predefined kinetic rate. For this, we employ a reaction-diffusion approach based on a reactive dynamical density-functional theory (R-DDFT) and study the effect of a binary (two-state) switching of the size of colloids interacting with a Gaussian pair potential. The rate of switching interpolates between a near-equilibrium binary Gaussian mixture at low rates and an effective monodisperse liquid for large rates, hence strongly affecting the one-body density profiles, adsorption, and pressure at confining walls. Importantly, we demonstrate that sufficiently fast switching impedes the phase separation of an (in equilibrium) unstable liquid, allowing the dynamic control of the degree of mixing/condensation and local micro-structuring in a cellular confinement by tuning the switching rate.

preprint2020arXiv

Enhanced Catalytic Activity of Gold@Polydopamine Nanoreactors with Multi-compartment Structure Under NIR Irradiation

Photothermal conversion (PTC) nanostructures have great potential for applications in many fields, and therefore, they have attracted tremendous attention. However, the construction of a PTC nanoreactor with multi-compartment structure to achieve the combination of unique chemical properties and structural feature is still challenging due to the synthetic difficulties. Herein, we designed and synthesized a catalytically active, PTC gold (Au)@polydopamine (PDA) nanoreactor driven by infrared irradiation using assembled PS-b-P2VP nanosphere as soft template. The particles exhibit multi-compartment structure which is revealed by 3D electron tomography characterization technique. They feature permeable shells with tunable shell thickness. Full kinetics for the reduction reaction of 4-nitrophenol has been investigated using these particles as nanoreactors and compared with other reported systems. Notably, a remarkable acceleration of the catalytic reaction upon near-infrared irradiation is demonstrated, which reveals for the first time the importance of the synergistic effect of photothermal conversion and complex inner structure to the kinetics of the catalytic reduction. The ease of synthesis and fresh insights into catalysis will promote a new platform for novel nanoreactor studies.

preprint2020arXiv

Modeling of stimuli-responsive nanoreactors: rational rate control towards the design of colloidal enzymes

In modern applications of heterogeneous liquid-phase nanocatalysis, the catalysts (e.g., metal nanoparticles) need to be typically affixed to a colloidal carrier system for stability and easy handling. "Passive carriers" (e.g., simple polyelectrolytes) serve for a controlled synthesis of the nanoparticles and prevent coagulation during catalysis. Recently, however, hybrid conjugates of nanoparticles and synthetic thermosensitive polymers have been developed that enable to change the catalytic activity of the nanoparticles by external triggers. In particular, nanoparticles embedded in a stimuli-responsive network made from poly(N-isopropylacrylamide) (PNIPAM) have become the most-studied examples of such hybrids. It has been demonstrated that the permeability of the polymer network and thus the reactant flux can be switched and controlled by external stimuli. Such "active carriers" may thus be viewed as true nanoreactors that open up new design routes in nano-catalysis and elevate synthesis to create highly selective, programmable "colloidal enzymes". However, only a comprehensive understanding of these materials on all time and length scales can lead to a rational design of future, highly functional materials. Here we review the current state of the theoretical and multi-scale simulation approaches, aiming at a fundamental understanding of these nanoreactors. In particular, we summarize a theoretical approach for reaction rates of surface-catalyzed bimolecular reactions in responsive nanoreactors in terms of the key material parameters, the polymer shell permeability P and the reactant partition ratio K. We discuss recent computer simulation studies of both atomistic and coarse-grained polymer models in which these quantities have been characterized in some detail. We conclude with an outlook on selected open questions and future theoretical challenges in nanoreactor modeling.

preprint2020arXiv

Structural and Transport Properties of Li/S Battery Electrolytes: Role of the Polysulfide Species

Lithium--sulfur (Li/S) batteries are regarded as one of the most promising energy storage devices beyond lithium-ion batteries because of their high energy density of 2600 Wh/kg and an affordable cost of sulfur. Meanwhile, some challenges inherent to Li/S batteries remain to be tackled, for instance, the polysulfide (PS) shuttle effect, the irreversible solidification of Li$_2$S, and the volume expansion of the cathode material during discharge. On the molecular level, these issues originate from the structural and solubility behavior of the PS species in bulk and in the electrode confinement. In this study, we use classical molecular dynamics (MD) simulations to develop a working model for PS of different chain lengths in applied electrolyte solutions of lithium bistriflimide (LiTFSI) in 1,2-dimethoxyethane (DME) and 1,3-dioxolane (DOL) mixtures. We investigate conductivities, diffusion coefficients, solvation structures, and clustering behavior and verify our simulation model with experimental measurements available in literature and newly performed by us. Our results show that diffusion coefficients and conductivities are significantly influenced by the chain length of PS. The conductivity contribution of the short chains, like Li$_2$S$_4$, is lower than of longer PS chains, such as Li$_2$S$_6$ or Li$_2$S$_8$, despite the fact that the diffusion coefficient of Li$_2$S$_4$ is higher than for longer PS chains. The low conductivity of Li$_2$S$_4$ can be attributed to its low degree of dissociation and even to a formation of large clusters in the solution. It is also found that an addition of 1 M LiTFSI into PS solutions considerably reduces the clustering behavior. Our simulation model enables future systematic studies in various solvating and confining systems for the rational design of Li/S electrolytes.

preprint2020arXiv

Thermal compaction of disordered and elastin-like polypeptides: a temperature-dependent, sequence-specific coarse-grained simulation model

Elastin-like polypeptides (ELPs) undergo a sharp solubility transition from low temperature solvated phases to coacervates at elevated temperatures, driven by the increased strength of hydrophobic interactions at higher temperatures. The transition temperature, or 'cloud point', critically depends on sequence composition, sequence length, and concentration of the ELPs. In this work, we present a temperature-dependent, implicit solvent, sequence-specific coarse-grained (CG) simulation model that reproduces the transition temperatures as a function of sequence length and guest residue identity of various experimentally probed ELPs to appreciable accuracy. Our model builds upon the self-organized polymer model introduced recently for intrinsically disordered polypeptides (SOP-IDP), and introduces a semi-empirical functional form for the temperature-dependence of hydrophobic interactions. In addition to the fine performance for various ELPs, we demonstrate the ability of our model to capture the thermal compactions in dominantly hydrophobic intrinsically disordered polypeptides (IDPs), consistent with experimental scattering data. With the high computational efficiency afforded by the CG representation, we envisage that the model will be ideally suited for simulations of large-scale structures such as ELP networks and hydrogels, as well as agglomerates of IDPs.

preprint2019arXiv

Tuning the permeability of dense membranes by shaping nanoscale potentials

The 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 non-monotonic behavior and the occurrence of extreme values of the permeability is far from trivial and triggered by a strong anti-correlation and substantial (orders of magnitude) cancellation between penetrant partitioning and diffusivity, especially within dense and highly attractive membranes.

preprint2019arXiv

Variational implicit-solvent predictions of the dry-wet transition pathways for ligand-receptor binding and unbinding kinetics

Ligand-receptor binding and unbinding are fundamental biomolecular processes and particularly essential to drug efficacy. Environmental water fluctuations, however, impact the corresponding thermodynamics and kinetics and thereby challenge theoretical descriptions. Here, we devise a holistic, implicit-solvent, multi-method approach to predict the (un)binding kinetics for a generic ligand-pocket model. We use the variational implicit-solvent model (VISM) to calculate the solute-solvent interfacial structures and the corresponding free energies, and combine the VISM with the string method to obtain the minimum energy paths and transition states between the various metastable ('dry' and 'wet') hydration states. The resulting dry-wet transition rates are then used in a spatially-dependent multi-state continuous-time Markov chain Brownian dynamics simulations, and the related Fokker-Planck equation calculations, of the ligand stochastic motion, providing the mean first-passage times for binding and unbinding. We find the hydration transitions to significantly slow down the binding process, in semi-quantitative agreement with existing explicit-water simulations, but significantly accelerate the unbinding process. Moreover, our methods allow the characterization of non-equilibrium hydration states of pocket and ligand during the ligand movement, for which we find substantial memory and hysteresis effects for binding versus unbinding. Our study thus provides a significant step forward towards efficient, physics-based interpretation and predictions of the complex kinetics in realistic ligand-receptor systems.