Tracer dynamics in polymer networks: Generalized Langevin description.

Tracer diffusion in polymer networks and hydrogels is relevant in biology and technology, while it also constitutes an interesting model process for the dynamics of molecules in fluctuating, heterogeneous soft matter. Here, we systematically study the time-dependent dynamics and (non-Markovian) memory effects of tracers in polymer networks based on (Markovian) implicit-solvent Langevin simulations. In particular, we consider spherical tracer solutes at high dilution in regular, tetrafunctional bead-spring polymer networks and control the tracer-network Lennard-Jones (LJ) interactions and the polymer density. Based on the analysis of the memory (friction) kernels, we recover the expected long-time transport coefficients and demonstrate how the short-time tracer dynamics, polymer fluctuations, and the viscoelastic response are interlinked. Furthermore, we fit the characteristic memory modes of the tracers with damped harmonic oscillations and identify LJ contributions, bond vibrations, and slow network relaxations. Tuned by the LJ interaction parameter, these modes enter the kernel with an approximately linear to quadratic scaling, which we incorporate into a reduced functional form for convenient tracer memory interpolation and extrapolation. This eventually leads to highly efficient simulations utilizing the generalized Langevin equation, in which the polymer network acts as an additional thermal bath with a tunable intensity.

Photo-excited charge transfer from adamantane to electronic bound states in water.

Aqueous nanodiamonds illuminated by UV light produce free solvated electrons, which may drive high-energy reduction reactions in water. However, the influence of water conformations on the excited-state electron-transfer mechanism are still under debate. In this work, we offer a theoretical study of charge-transfer states in adamantane-water structures obtained by linear-response time-dependent density-functional theory. Small water clusters with broken hydrogen bonds are found to efficiently bind the electron from adamantane. A distinction is made with respect to the nature of the water clusters: some bind the electron in a water cavity, others along a strong permanent total dipole. These two types of bound states are more strongly binding, the higher their electron affinity and their positive electrostatic potential, the latter being dominated by the energy of the lowest unoccupied molecular orbital of the isolated water clusters. Structural sampling in a thermal equilibrium at room temperature via molecular dynamics snapshots confirms under which conditions the underlying waters clusters can occur and verifies that broken hydrogen bonds in the water network close to adamantane can create traps for the solvated electron.

Unveiling the Borohydride Ion through Force-Field Development.

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.

External field-driven property localization in liquids of responsive macromolecules.

We explore theoretically the effects of external potentials on the spatial distribution of particle properties in a liquid of explicitly responsive macromolecules. In particular, we focus on the bistable particle size as a coarse-grained internal degree of freedom (DoF, or "property"), σ, that moves in a bimodal energy landscape, in order to model the response of a state-switching (big-to-small) macromolecular liquid to external stimuli. We employ a mean-field density functional theory (DFT) that provides the full inhomogeneous equilibrium distributions of a one-component model system of responsive colloids (RCs) interacting with a Gaussian pair potential. For systems confined between two parallel hard walls, we observe and rationalize a significant localization of the big particle state close to the walls, with pressures described by an exact RC wall theorem. Application of more complex external potentials, such as linear (gravitational), osmotic, and Hamaker potentials, promotes even stronger particle size segregation, in which macromolecules of different size are localized in different spatial regions. Importantly, we demonstrate how the degree of responsiveness of the particle size and its coupling to the external potential tune the position-dependent size distribution. The DFT predictions are corroborated by Brownian dynamics simulations. Our study highlights the fact that particle responsiveness can be used to localize liquid properties and therefore helps to control the property- and position-dependent function of macromolecules, e.g., in biomedical applications.

Salt Effects on Caffeine across Concentration Regimes.

Salts affect the solvation thermodynamics of molecules of all sizes; the Hofmeister series is a prime example in which different ions lead to salting-in or salting-out of aqueous proteins. Early work of Tanford led to the discovery that the solvation of molecular surface motifs is proportional to the solvent accessible surface area (SASA), and later studies have shown that the proportionality constant varies with the salt concentration and type. Using multiscale computer simulations combined with vapor-pressure osmometry on caffeine-salt solutions, we reveal that this SASA description captures a rich set of molecular driving forces in tertiary solutions at changing solute and osmolyte concentrations. Central to the theoretical work is a new potential energy function that depends on the instantaneous surface area, salt type, and concentration. Used in, e.g., Monte Carlo simulations, this allows for a highly efficient exploration of many-body interactions and the resulting thermodynamics at elevated solute and salt concentrations.

Synergistic chemomechanical dynamics of feedback-controlled microreactors.

The experimental control of synergistic chemomechanical dynamics of catalytically active microgels (microreactors) is a key prerequisite for the design of adaptive and biomimetic materials. Here, we report a minimalistic model of feedback-controlled microreactors based on the coupling between the hysteretic polymer volume phase transition and a volume-controlled permeability for the internal chemical conversion. We categorize regimes of mono- and bistability, excitability, damped oscillations, as well as sustained oscillatory states with tunable amplitude, as indicated by experiments and representable by the FitzHugh-Nagumo dynamics for neurons. We summarize the features of such a "colloidal neuron" in bifurcation diagrams with respect to microgel design parameters, such as permeability and relaxation times, as a guide for experimental synthesis.

Mean-field models for the chemical fueling of transient soft matter states.

The chemical fueling of transient states (CFTS) is a powerful process to control the nonequilibrium structuring and the homeostatic function of adaptive soft matter systems. Here, we introduce a simple mean-field model of CFTS based on the activation of metastable equilibrium states in a tilted 'Landau' bistable energy landscape along a coarse-grained reaction coordinate (or 'order parameter') triggered by a nonmonotonic two-step chemical fueling reaction. Evaluation of the model in the quasi-static (QS) limit-valid for fast system relaxation-allows us to extract useful analytical laws for the critical activation concentration and duration of the transient states in dependence of physical parameters, such as rate constants, fuel concentrations, and the system's distance to its equilibrium transition point. We apply our model in the QS limit explicitly to recent experiments of CFTS of collapsing responsive microgels and find a very good performance with only a few global and physically interpretable fitting parameters, which can be employed for programmable material design. Moreover, our model framework also allows a thermodynamic analysis of the energy and performed work in the system. Finally, we go beyond the QS limit, where the system's response is slow and retarded versus the chemical reaction, using an overdamped Smoluchowski approach. The latter demonstrates how internal system time scales can be used to tune the time-dependent behavior and programmed delay of the transient states in full nonequilibrium.

Intramolecular structural heterogeneity altered by long-range contacts in an intrinsically disordered protein.

Short-range interactions and long-range contacts drive the 3D folding of structured proteins. The proteins' structure has a direct impact on their biological function. However, nearly 40% of the eukaryotes proteome is composed of intrinsically disordered proteins (IDPs) and protein regions that fluctuate between ensembles of numerous conformations. Therefore, to understand their biological function, it is critical to depict how the structural ensemble statistics correlate to the IDPs' amino acid sequence. Here, using small-angle X-ray scattering and time-resolved Förster resonance energy transfer (trFRET), we study the intramolecular structural heterogeneity of the neurofilament low intrinsically disordered tail domain (NFLt). Using theoretical results of polymer physics, we find that the Flory scaling exponent of NFLt subsegments correlates linearly with their net charge, ranging from statistics of ideal to self-avoiding chains. Surprisingly, measuring the same segments in the context of the whole NFLt protein, we find that regardless of the peptide sequence, the segments' structural statistics are more expanded than when measured independently. Our findings show that while polymer physics can, to some level, relate the IDP's sequence to its ensemble conformations, long-range contacts between distant amino acids play a crucial role in determining intramolecular structures. This emphasizes the necessity of advanced polymer theories to fully describe IDPs ensembles with the hope that it will allow us to model their biological function.

How Regiochemistry Influences Aggregation Behavior and Charge Transport in Conjugated Organosulfur Polymer Cathodes for Lithium-Sulfur Batteries.

For lithium-sulfur (Li-S) batteries to become competitive, they require high stability and energy density. Organosulfur polymer-based cathodes have recently shown promising performance due to their ability to overcome common limitations of Li-S batteries, such as the insulating nature of sulfur. In this study, we use a multiscale modeling approach to explore the influence of the regiochemistry of a conjugated poly(4-(thiophene-3-yl)benzenethiol) (PTBT) polymer on its aggregation behavior and charge transport. Classical molecular dynamics simulations of the self-assembly of polymer chains with different regioregularity show that a head-to-tail/head-to-tail regularity can form a well-ordered crystalline phase of planar chains allowing for fast charge transport. Our X-ray diffraction measurements, in conjunction with our predicted crystal structure, confirm the presence of crystalline phases in the electropolymerized PTBT polymer. We quantitatively describe the charge transport in the crystalline phase in a band-like regime. Our results give detailed insights into the interplay between microstructural and electrical properties of conjugated polymer cathode materials, highlighting the effect of polymer chain regioregularity on its charge transport properties.

Machine Learning Frontier Orbital Energies of Nanodiamonds.

Nanodiamonds have a wide range of applications including catalysis, sensing, tribology, and biomedicine. To leverage nanodiamond design via machine learning, we introduce the new data set ND5k, consisting of 5089 diamondoid and nanodiamond structures and their frontier orbital energies. ND5k structures are optimized via tight-binding density functional theory (DFTB) and their frontier orbital energies are computed using density functional theory (DFT) with the PBE0 hybrid functional. From this data set we derive a qualitative design suggestion for nanodiamonds in photocatalysis. We also compare recent machine learning models for predicting frontier orbital energies for similar structures as they have been trained on (interpolation on ND5k), and we test their abilities to extrapolate predictions to larger structures. For both the interpolation and extrapolation task, we find the best performance using the equivariant message passing neural network PaiNN. The second best results are achieved with a message passing neural network using a tailored set of atomic descriptors proposed here.

Liquid structure of bistable responsive macromolecules using mean-field density-functional theory.

Macromolecular crowding typically applies to biomolecular and polymer-based systems in which the individual particles often feature a two-state folded/unfolded or coil-to-globule transition, such as found for proteins and peptides, DNA and RNA, or supramolecular polymers. Here, we employ a mean-field density functional theory (DFT) of a model of soft and bistable responsive colloids (RCs) in which the size of the macromolecule is explicitly resolved as a degree of freedom living in a bimodal 'Landau' energy landscape (exhibiting big and small states), thus directly responding to the crowding environment. Using this RC-DFT we study the effects of self-crowding on the liquid bulk structure and thermodynamics for different energy barriers and softnesses of the bimodal energy landscape, in conditions close to the coil-to-globule transition. We find substantial crowding effects on the internal distributions, a complex polydispersity behavior, and quasi-universal compression curves for increasing (generalized) packing fractions. Moreover, we uncover distinct signatures of bimodal versus unimodal behavior in the particle compression. Finally, the analysis of the pair structure - derived from the test particle route - reveals that the microstructure of the liquid is quite inhomogeneous due to local depletion effects, tuneable by particle softness.

Feedback-controlled solute transport through chemo-responsive polymer membranes.

Polymer membranes are typically assumed to be inert and nonresponsive to the flux and density of the permeating particles in transport processes. Here, we theoretically study the consequences of membrane responsiveness and feedback on the steady-state force-flux relations and membrane permeability using a nonlinear-feedback solution-diffusion model of transport through a slab-like membrane. Therein, the solute concentration inside the membrane depends on the bulk concentration, c0, the driving force, f, and the polymer volume fraction, ϕ. In our model, the solute accumulation in the membrane causes a sigmoidal volume phase transition of the polymer, changing its permeability, which, in return, affects the membrane's solute uptake. This feedback leads to nonlinear force-flux relations, j(f), which we quantify in terms of the system's differential permeability, Psys Δ∝dj/df. We find that the membrane feedback can increase or decrease the solute flux by orders of magnitude, triggered by a small change in the driving force and largely tunable by attractive vs repulsive solute-membrane interactions. Moreover, controlling the inputs, c0 and f, can lead to the steady-state bistability of ϕ and hysteresis in the force-flux relations. This work advocates that the fine-tuning of the membrane's chemo-responsiveness will enhance the nonlinear transport control features, providing great potential for future (self-)regulating membrane devices.

Active responsive colloids coupled to different thermostats.

We introduce a model of active responsive colloids (ARCs) in which an internal degree of freedom (DoF) of a single colloidal particle is "activated" by coupling it to a different thermostat than for the translational DoFs. As for the responsive internal DoF, we consider specifically the size (diameter) of the spherical particles, which is confined by a harmonic parent potential being either entropic or energetic in nature. The ARCs interact via a repulsive Hertzian pair potential, appropriate to model hydrogels or elastic colloids, and are studied for various densities using Brownian dynamics simulations. We tune the internal activity in the nonequilibrium steady state by scanning through a wide range of internal temperatures, both smaller ("colder") and larger ("hotter") than the translational temperature. The results show a rich and intriguing behavior for the emergent property distributions, colloidal pair structure, and the diffusive translational dynamics controlled by the internal activity, substantially depending on whether the internal DoF moves in an entropic or energetic potential. We discuss theoretical thermal limits and phenomenological models which can explain some of the nonequilibrium trends qualitatively. Our study indicates that particle dynamical polydispersity as well as the structure and dynamics of dense macromolecular suspensions can be vastly tuned by internal activity in terms of internal "hot" or "cold" fluctuating states.

Constructing Binder- and Carbon Additive-Free Organosulfur Cathodes Based on Conducting Thiol-Polymers through Electropolymerization for Lithium-Sulfur Batteries.

Herein, the concept of constructing binder- and carbon additive-free organosulfur cathode was proved based on thiol-containing conducting polymer poly(4-(thiophene-3-yl) benzenethiol) (PTBT). The PTBT featured the polythiophene-structure main chain as a highly conducting framework and the benzenethiol side chain to copolymerize with sulfur and form a crosslinked organosulfur polymer (namely S/PTBT). Meanwhile, it could be in-situ deposited on the current collector by electro-polymerization, making it a binder-free and free-standing cathode for Li-S batteries. The S/PTBT cathode exhibited a reversible capacity of around 870 mAh g-1 at 0.1 C and improved cycling performance compared to the physically mixed cathode (namely S&PTBT). This multifunction cathode eliminated the influence of the additives (carbon/binder), making it suitable to be applied as a model electrode for operando analysis. Operando X-ray imaging revealed the remarkable effect in the suppression of polysulfides shuttle via introducing covalent bonds, paving the way for the study of the intrinsic mechanisms in Li-S batteries.

Effects of oxidative adsorbates and cluster formation on the electronic structure of nanodiamonds.

Nanodiamonds (NDs) are modern high-potential materials relevant for applications in biomedicine, photocatalysis, and various other fields. Their electronic surface properties, especially in the liquid phase, are key to their function in the applications, but we show that they are sensitively modified by their interactions with the environment. Two important interaction modes are those with oxidative aqueous adsorbates as well as ND self-aggregation towards the formation of ND clusters. For planar diamond surfaces it is known that the electron density migrates from the diamond towards oxidative adsorbates, which is known as transfer doping. Here, we quantify this effect for highly curved NDs of varying sizes (35-147 C atoms) and surface terminations (H, OH, F), focusing on their interactions with the most abundant aqueous oxidative adsorbates (H3 O+ , O2 , O3 ). We prove that the concept of transfer doping stays valid for the case of the high-curvature NDs and can be tuned via the ND's specific properties. Secondly, we investigate the electronic structures of clusters of NDs which are known to form in particular in aqueous dispersions. Upon cluster formation, we find that the optical gaps of the structures are significantly reduced, which explains why different experimental values were obtained for the optical gap of the same structures, and the cluster's LUMO shapes resemble atom-type orbitals, as in the case of isolated spherical NDs. Our findings have implications for ND applications as photocatalysts or electronic devices, where the specific electronic properties are key to the functionality of the ND material.

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 characterizable 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 computer simulations in conjunction with approximative statistical mechanics and relative entropy (Gibbs-Shannon and Kullback-Leibler) concepts. In the case of fast growing polymers, we indeed report increased nonequilibrium free energies ΔFneq of several kBT compared to equilibrium and near-equilibrium, slowly growing chains. Interestingly, ΔFneq 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.

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 [Archer et al., Phys. Rev. E: Stat., Nonlinear, Soft Matter Phys., 2007, 75, 040501] 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.

Active responsive colloids driven by intrinsic dichotomous noise.

We study the influence of intrinsic noise on the structure and dynamics of responsive colloids (RCs), which actively change their size and mutual interactions. The colloidal size is explicitly resolved in our RC model as an internal degree of freedom (DOF) in addition to the particle translation. A Hertzian pair potential between the RCs leads to repulsion and shrinking of the particles, resulting in an explicit responsiveness of the system to self-crowding. To render the colloids active, their size is internally driven by a dichotomous noise, randomly switching ("breathing") between growing and shrinking states with a predefined rate, as motivated by recent experiments on synthetic active colloids. The polydispersity of this dichotomous active responsive colloid (D-ARC) model can be tuned by the parameters of the noise. Utilizing stochastic computer simulations, we study crowding effects on the spatial distributions, relaxation times, and self-diffusion of dense suspensions of the D-ARCs. We find a substantial influence of the "built-in" intrinsic noise on the system's behavior, in particular, transitions from unimodal to bimodal size distributions for an increasing colloid density as well as intrinsic noise-modified diffusive translational dynamics. We conclude that controlling the noise of internal DOFs of a macromolecule or cell is a powerful tool for active colloidal materials to enable autonomous changes in the system's collective structure and dynamics towards the adaptation of macroscopic properties to external perturbations.

Combined first-principles statistical mechanics approach to sulfur structure in organic cathode hosts for polymer based lithium-sulfur (Li-S) batteries.

Polymer-based batteries that utilize organic electrode materials are considered viable candidates to overcome the common drawbacks of lithium-sulfur (Li-S) batteries. A promising cathode can be developed using a conductive, flexible, and free-standing polymer, poly(4-thiophen-3-yl)benzenethiol) (PTBT), as the sulfur host material. By a vulcanization process, sulfur is embedded into this polymer. Here, we present a combination of electronic structure theory and statistical mechanics to characterize the structure of the initial state of the charged cathode on an atomic level. We perform a stability analysis of differently sulfurized TBT dimers as the basic polymer unit calculated within density-functional theory (DFT) and combine this with a statistical binding model for the binding probability distributions of the vulcanization process. From this, we deduce sulfur chain length ("rank") distributions and calculate the average sulfur rank depending on the sulfur concentration and temperature. This multi-scale approach allows us to bridge the gap between the local description of the covalent bonding process and the derivation of the macroscopic properties of the cathode. Our calculations show that the main reaction of the vulcanization process leads to high-probability states of sulfur chains cross-linking TBT units belonging to different polymer backbones, with a dominant rank around n = 5. In contrast, the connection of adjacent TBT units of the same polymer backbone by a sulfur chain is the side reaction. These results are experimentally supported by Raman spectroscopy.

Controlling solvent quality by time: Self-avoiding sprints in nonequilibrium polymerization.

A fundamental paradigm in polymer physics is that macromolecular conformations in equilibrium can be described by universal scaling laws, being key for structure, dynamics, and function of soft (biological) matter and in the materials sciences. Here we reveal that during diffusion-influenced, nonequilibrium chain-growth polymerization, scaling laws change qualitatively, in particular, the growing polymers exhibit a surprising self-avoiding walk behavior in poor and θ solvents. Our analysis, based on monomer-resolved, off-lattice reaction-diffusion computer simulations, demonstrates that this phenomenon is a result of (i) nonequilibrium monomer density depletion correlations around the active polymerization site, leading to a locally directed and self-avoiding growth, in conjunction with (ii) chain (Rouse) relaxation times larger than the competing polymerization reaction time. These intrinsic nonequilibrium mechanisms are facilitated by fast and persistent reaction-driven diffusion ("sprints") of the active site, with analogies to pseudochemotactic active Brownian particles. Our findings have implications for time-controlled structure formation in polymer processing, as in, e.g., reactive self-assembly, photocrosslinking, and three-dimensional printing.