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.

Explicit and implicit modeling of nanobubbles in hydrophobic confinement

Water at normal conditions is a fluid thermodynamically close to the liquid-vapor phase coexistence and features a large surface tension. This combination can lead to interesting capillary phenomena on microscopic scales. Explicit water molecular dynamics (MD) computer simulations of hydrophobic solutes, for instance, give evidence of capillary evaporation on nanometer scales, i.e., the formation of nanometer-sized vapor bubbles (nanobubbles) between confining hydrophobic surfaces. This phenomenon has been exemplified for solutes with varying complexity, e.g., paraffin plates, coarse-grained homopolymers, biological and solid-state channels, and atomistically resolved proteins. It has been argued that nanobubbles strongly impact interactions in nanofluidic devices, translocation processes, and even in protein stability, function, and folding. As large-scale MD simulations are computationally expensive, the efficient multiscale modeling of nanobubbles and the prediction of their stability poses a formidable task to the'nanophysical' community. Recently, we have presented a conceptually novel and versatile implicit solvent model, namely, the variational implicit solvent model (VISM), which is based on a geometric energy functional. As reviewed here, first solvation studies of simple hydrophobic solutes using VISM coupled with the numerical level-set scheme show promising results, and, in particular, capture nanobubble formation and its subtle competition to local energetic potentials in hydrophobic confinement.

Água em condições normais consiste de um fluido termodinamicamente próximo à fase líquida-vapor exibindo alta tensão superficial. Esta combinação conduz a fenômenos capilares interessantes na escala microscópica. Simulações computacionais baseadas em técnicas de Dinâmica Molecular em solutos hidrofóbicos por exemplo fornecem evidências do fenômeno de evaporação capilar em escalas nanométricas dando origem à formação de bolhas nanométricas confinadas entre superfícies hidrofóbicas. Este fenômeno tem sido exemplificado em solutos de complexidade variável, ex placas de parafina, polímeros, canais biológicos e de estado sólido bem como proteínas tratadas atomicamente. Tem sido arguido quebolhas nanométricas possuem impacto significativo nos mecanismos fluídicos nanométricos e na estabilidade e dobramento de proteínas. Dado que simulações baseadas em dinâmica molecular são custosas computacionalmente, o desenvolvimento de uma modelagem multiescala eficiente sob o ponto de vista computacional impõe uma tarefa formidável à comunidadenano-física. Recentemente apresentamos um novo e versátil modelo de solvente denominado modelo implícito variacional de solvente (VISM) o qual é baseado em um funcional deenergia geométrica. Tal como apresentado aqui os primeiros estudos de solvatação de solutos hidrofóbicos simples usando VISM acoplados com esquemas numéricos de conjunto de níveis mostraram resultados promissores e em particular capturaram a formação de nano-bolhas e a sua competição com potenciais localmente energéticos em condições de confinamento hidrofóbico.