Search

Friction is a phenomenon that manifests across all spatial and temporal scales, from the molecular to the macroscopic scale. It describes the dissipation of energy from the motion of particles or abstract reaction coordinates and arises in the transition from a detailed molecular-level description to a simplified, coarse-grained model. It has long been understood that time-dependent (non-Markovian) friction effects are critical for describing the dynamics of many systems, but that they are notoriously difficult to evaluate for complex physical, chemical, and biological systems. In recent years, the development of advanced numerical friction extraction techniques and methods to simulate the generalized Langevin equation have enabled exploration of the role of time-dependent friction across all scales. We discuss recent applications of these friction extraction techniques and the growing understanding of the role of friction in complex equilibrium and non-equilibrium dynamic many-body systems., Comment: 26 pages, 6 figures

Water at negative pressures can remain in a metastable state for a surprisingly long time before it reaches equilibrium by cavitation, i.e. by the formation of vapor bubbles. The wide spread of experimentally measured cavitation pressures depending on water purity, surface contact angle and surface quality implicates the relevance of water cavitation in bulk, at surfaces and at surface defects for different systems. We formulate a kinetic model that includes all three different cavitation pathways and determine the needed nucleation attempt frequencies in bulk, at surfaces and at defects from atomistic molecular dynamics simulations. Our model reveals that cavitation occurs in pure bulk water only for defect-free hydrophilic surfaces with wetting contact angles below 50{\deg} to 60{\deg} and at pressures of the order of $-$100 MPa, depending only slightly on system size and observation time. Cavitation on defect-free surfaces occurs only for higher contact angles, with the typical cavitation pressure rising to about $-$30 MPa for very hydrophobic surfaces. Nanoscopic hydrophobic surface defects act as very efficient cavitation nuclei and can dominate the cavitation kinetics in a macroscopic system. In fact, a nanoscopic defect that hosts a pre-existing vapor bubble can raise the critical cavitation pressure much further. Our results explain the wide variation of experimentally observed cavitation pressures in synthetic and biological systems and highlight the importance of surface and defect mechanisms for the nucleation of metastable systems.

The air-water interface is a highly prevalent phase boundary with a far-reaching impact on natural and industrial processes. Water molecules behave differently at the interface compared to the bulk, exhibiting anisotropic orientational distributions, reduced intermolecular connectivity in the hydrogen bond network, and significantly slower dynamics. Despite many decades of research, the thickness of the structural anisotropy in the interfacial layer remains controversial, with a direct experimental measurement being absent. In this study, we utilise an advancement in non-linear vibrational spectroscopy to gain access to this important parameter. Combining phase-resolved sum- and difference-frequency generation (SFG and DFG) responses, we directly measure the decay in structural anisotropy of the air-water interface. We find a decay length of ~6-8\r{A}, in excellent agreement with depth-resolved SFG spectra calculated from ab initio parameterised molecular dynamics (MD) simulations. The result reveals surprisingly short anisotropic orientational correlations from the interfacial layer that are even shorter than in the bulk. Furthermore, the recorded SFG and DFG responses are decomposed into a vibrationally resonant and non-resonant contribution through isotopic exchange measurements. Through their separate analysis, we show that the resonant response is a sensitive probe of the structural anisotropy at the interface whereas the non-resonant contribution contains a significant isotropic contribution from the bulk and therefore only partially reports on the interfacial structure. This finding places stringent restrictions on the insight available through both purely non-resonant and second-order intensity studies., Comment: 18 pages, 4 figures, 1 table, supplementary information

The mean first-passage time (MFPT) is one standard measure for the reaction time in thermally activated barrier-crossing processes. While the relationship between MFPTs and phenomenological rate coefficients is known for systems that satisfy Markovian dynamics, it is not clear how to interpret MFPTs for experimental and simulation time-series data generated by non-Markovian systems. Here, we simulate a one-dimensional generalized Langevin equation (GLE) in a bistable potential and compare two related numerical methods for evaluating MFPTs: one that only incorporates information about first arrivals between subsequent states and is equivalent to calculating the waiting time, or dwell time, and one that incorporates information about all first-passages associated with a given barrier-crossing event and is therefore typically employed to enhance numerical sampling. In the Markovian limit, the two methods are equivalent. However, for significant memory times, the two methods suggest dramatically different reaction kinetics. By focusing on first-passage distributions, we systematically reveal the influence of memory-induced rapid state-recrossing on the MFPTs, which we compare to various other numerical or theoretical descriptions of reaction times. Overall, we demonstrate that it is necessary to consider full first-passage distributions, rather than just the mean barrier-crossing kinetics when analyzing non-Markovian time series data., Comment: 7 pages, 5 figures

Deviations between macrorheological and particle-based microrheological measurements are often considered a nuisance and neglected. We study aqueous poly(ethylene oxide) (PEO) hydrogels for varying PEO concentrations and chain lengths that contain microscopic tracer particles and show that these deviations in fact reveal the nanoscopic viscoelastic properties of the particle-hydrogel interface. Based on the transient Stokes equation, we first demonstrate that the deviations are not due to finite particle radius, compressibility or surface-slip effects. Small-angle neutron scattering rules out hydrogel heterogeneities. Instead, we show that a generalized Stokes-Einstein relation, accounting for a nanoscopic interfacial shell around tracers with viscoelastic properties that significantly deviate from bulk, consistently explains our macrorheological and microrheological measurements. The extracted shell diameter is comparable with the PEO end-to-end distance, indicating the importance of dangling chain ends. Our methodology reveals the nanoscopic interfacial rheology of hydrogels and is generally applicable to different kinds of viscoelastic fluids and particles., Comment: 25 pages, 4 figures, 52 pages of Supporting Information

Molecular isomerization kinetics in liquid solvents are determined by a complex interplay between the friction acting on a rotating dihedral due to interactions with the solvent, internal dissipation effects (also known as internal friction), the viscosity of the solvent, and the free energy profile over which a dihedral rotates. Currently, it is not understood how these quantities are related at the molecular scale. Here, we combine molecular dynamics simulations of isomerizing n-alkane chains and dipeptide molecules in mixed water-glycerol solvents with memory-kernel extraction techniques to directly evaluate the frequency-dependent friction acting on a rotating dihedral. We extract the friction and isomerization times over a range of glycerol concentrations and accurately evaluate the relationships between solvent viscosity, isomerization kinetics, and dihedral friction. We show that the total friction acting on a rotating dihedral does not scale linearly with solvent viscosity, thus violating the Stokes-Einstein relation. Additionally, we demonstrate that the kinetics of isomerization are significantly faster compared to the Kramers prediction in the overdamped limit. We suggest that isomerization kinetics are determined by the multi-time-scale friction coupling between a rotating dihedral and its solvent environment, which results in non-Markovian kinetic speed-up effects., Comment: 9 pages, 4 figures

Based on a Hamiltonian that incorporates the elastic coupling between a tracer and active particles, we derive a generalized Langevin model for the non-equilibrium mechanical response of active viscoelastic biomatter. Our model accounts for the power-law viscoelastic response of the embedding polymeric network as well as for the non-equilibrium energy transfer between active and tracer particles. Our analytical expressions for the frequency-dependent response function and the positional autocorrelation function agree nicely with experimental data for red blood cells and actomyosin networks with and without ATP. The fitted effective active-particle temperature, elastic constants and effective friction coefficients of our model allow straightforward physical interpretation., Comment: 6 pages, 2 figures; supplemental material (16 pages, 5 figures) is appended

We introduce a generalized Langevin model system for different non-Markovian effects in the well and barrier regions of a potential, and use it to numerically study the dependence of the barrier-crossing time. In the appropriate limits, our model interpolates between the theoretical barrier-crossing-time predictions by Grote and Hynes (GH), as well as by Pollak et al., which for a single barrier memory time can differ by several orders of magnitude. Our model furthermore allows to test an analytic rate theory for space-inhomogeneous memory, which disagrees with our numerical results in the long well-memory regime. In this regime, we find that short barrier memory decreases the barrier-crossing time as compared to long barrier memory. This is in contrast with the short well-memory regime, where both our numerical results and GH theory predict an acceleration of the barrier crossing time with increasing barrier memory time. Both effects, the `Markovian-barrier acceleration' and GH `non-Markovian-barrier acceleration' can be understood from a committor analysis. Our model combines finite relaxation times of orthogonal degrees of freedom with a space-inhomogeneous coupling to such degrees, and represents a step towards more realistic modeling of physical reaction coordinates.

Signatures of solvated excess protons in infrared difference absorption spectra, such as the continuum band between the water bend and stretch bands, have been experimentally known for a long time, but the theoretical basis for linking spectral signatures with the microscopic proton-transfer mechanism so far relied on normal-mode analysis. We analyze the excess-proton dynamics in ab initio molecular-dynamics simulations of aqueous hydrochloric acid solutions by trajectory-decomposition techniques. The continuum band in the 2000 - 3000 cm$^{-1 }$ range is shown to be due to normal-mode oscillations of temporary H$_3$O$^+$ complexes. An additional prominent peak at 400 cm$^{-1}$ reports on the coupling of excess-proton motion to the relative vibrations of the two flanking water molecules. The actual proton transfer between two water molecules, which for large water separations involves crossing of a barrier and thus is not a normal mode, is characterized by two characteristic time scales: Firstly, the waiting time for transfer to occur in the range of 200 - 300 fs, which leads to a broad weak shoulder around ~100 cm$^{-1}$, consistent with our experimental THz spectra. Secondly, the mean duration of a transfer event of about 14 fs, which produces a rather well-defined spectral contribution around 1200 cm$^{-1}$ and agrees in location and width with previous experimental mid-infrared spectra.

We introduce a hybrid projection scheme that combines linear Mori projection and conditional Zwanzig projection techniques and use it to derive a Generalized Langevin Equation (GLE) for a general interacting many-body system. The resulting GLE includes i) explicitly the potential of mean force (PMF) that describes the equilibrium distribution of the system in the chosen space of reaction coordinates, ii) a random force term that explicitly depends on the initial state of the system, and iii) a memory friction contribution that splits into two parts: a part that is linear in the past reaction-coordinate velocity and a part that is in general non-linear in the past reaction coordinates but does not depend on velocities. Our hybrid scheme thus combines all desirable properties of the Zwanzig and Mori projection schemes. The non-linear memory friction contribution is shown to be related to correlations between the reaction-coordinate velocity and the random force. We present a numerical method to compute all parameters of our GLE, in particular, the non-linear memory friction function and the random force distribution, from a trajectory in reaction coordinate space. We apply our method to the dihedral-angle dynamics of a butane molecule in water obtained from atomistic molecular dynamics simulations. For this example, we demonstrate that non-linear memory friction is present and that the random force exhibits significant non-Gaussian corrections. We also present the derivation of the GLE for multidimensional reaction coordinates that are general functions of all positions in the phase space of the underlying many-body system; this corresponds to a systematic coarse-graining procedure that preserves not only the correct equilibrium behavior but also the correct dynamics of the coarse-grained system.

From ab initio simulations of liquid water, the time-dependent friction functions and time-averaged non-linear effective bond potentials for the OH stretch and HOH bend vibrations are extracted. The obtained friction exhibits adiabatic contributions at and below the vibrational time scales, but also much slower non-adiabatic contributions, reflecting homogeneous and inhomogeneous line broadening mechanisms, respectively. Compared to the gas phase, hydration softens both stretch and bend potentials, which by itself would lead to a red-shift of the corresponding vibrational bands. In contrast, non-adiabatic friction contributions cause a spectral blue shift. For the stretch mode, the potential effect dominates and thus a significant red shift when going from gas to the liquid phase results. For the bend mode, potential and non-adiabatic friction effects are of comparable magnitude, so that a slight blue shift results, in agreement with well-known but puzzling experimental findings. The observed line broadening is shown to be roughly equally caused by adiabatic and non-adiabatic friction contributions for both, the stretch and bend modes in liquid water. Thus, the understanding of infrared vibrational frequencies and line shapes is considerably advanced by the quantitative analysis of the time-dependent friction that acts on vibrational modes in liquids.

A stochastic theory is developed to predict the spectral signature of proton-transfer processes and applied to infrared spectra computed from ab initio molecular-dynamics simulations of a single H$_5$O$_2{}^{+}$ cation. By constraining the oxygen atoms to a fixed distance, this system serves as a tunable model for general proton-transfer processes with variable barrier height. Three spectral contributions at distinct frequencies are identified and analytically predicted: the quasi-harmonic motion around the most probable configuration, amenable to normal-mode analysis, the contribution due to transfer paths when the proton moves over the barrier and a shoulder for low frequencies stemming from the stochastic transfer-waiting-time distribution; the latter two contributions are not captured by normal-mode analysis but exclusively report on the proton-transfer kinetics. In accordance with reaction kinetic theory, the transfer-waiting-contribution frequency depends inverse exponentially on the barrier height, whereas the transfer-path-contribution frequency is rather insensitive to the barrier height.

Friction in liquids arises from conservative forces between molecules and atoms. Although the hydrodynamics at the nanoscale is subject of intense research and despite the enormous interest in the non-Markovian dynamics of single molecules and solutes, the onset of friction from the atomistic scale so far could not be demonstrated. Here, we fill this gap based on frequency-resolved friction data from high-precision simulations of three prototypical liquids, including water. Combining with rigorous theoretical arguments, we show that friction in liquids emerges abruptly at a characteristic frequency, beyond which viscous liquids appear as non-dissipative, elastic solids; as a consequence, its origin is non-local in time. Concomitantly, the molecules experience Brownian forces that display persistent correlations and long-lasting memory. A critical test of the generalised Stokes-Einstein relation, mapping the friction of single molecules to the viscoelastic response of the macroscopic sample, disproves the relation for Newtonian fluids, but substantiates it exemplarily for water and a moderately supercooled liquid. The employed approach is suitable to yield novel insights into vitrification mechanisms and the intriguing mechanical properties of soft materials.

Energy dissipation in water is very fast and more efficient than in many other liquids. This behavior is commonly attributed to the intermolecular interactions associated with hydrogen bonding. Here, we investigate the dynamic energy flow in the hydrogen-bond network of liquid water by a pump-probe experiment. We resonantly excite intermolecular degrees of freedom with ultrashort single-cycle terahertz pulses and monitor its Raman response. By using ultrathin sample-cell windows, a background-free bipolar signal whose tail relaxes mono-exponentially is obtained. The relaxation is attributed to the molecular translational motions, using complementary experiments, force-field and ab initio molecular dynamics simulations. They reveal an initial coupling of the terahertz electric field to the molecular rotational degrees of freedom whose energy is rapidly transferred, within the excitation pulse duration, to the restricted-translational motion of neighboring molecules. This rapid energy transfer may be rationalized by the strong anharmonicity of the intermolecular interactions.

Using classical as well as ab-initio molecular dynamics simulations, we calculate the frequency-dependent shear viscosity of pure water and water-glycerol mixtures. In agreement with recent experiments, we find deviations from Newtonian-fluid behavior in the THz regime. Based on an extension of the Maxwell model, we introduce a viscoelastic model to describe the observed viscosity spectrum of pure water. We find four relaxation modes in the spectrum which we attribute to i) hydrogen-bond network topology changes, ii) hydrogen-bond stretch vibrations of water pairs, iii) collective vibrations of water molecule triplets, and iv) librational excitations of individual water molecules. Our model quantitatively describes the viscoelastic response of liquid water on short timescales, where the hydrodynamic description via a Newtonian-fluid model breaks down., Comment: Supplemental material is provided as ancillary file

Absorption spectra of liquid water at 300 K are calculated from both classical and density functional theory molecular dynamics simulation data, which together span from 1 MHz to hundreds of THz, agreeing well with experimental data qualitativley and quantitavely over the entire range, including the IR modes, the microwave peak, and the intermediate THz bands. The spectra are decomposed into single-molecular and collective components as well as into components due to molecular reorientations and changes in induced intramolecular dipole moments. These decompositions shed new light on the motions underlying the librational and translational (hydrogen-bond stretch) bands at 20 and 5 THz respectively: interactions between donor protons and acceptor lone pair electrons are shown to be important for the line shape in both librational and translational regimes, and in- and out-of-phase librational dimer modes are observed and explored., Comment: Main text: 7 pages with 4 figures. Supplement: 8 pages with 5 figures. Submitted to Nature Communications on February 29, 2020

We investigate non-Markovian barrier-crossing kinetics of a massive particle in one dimension in the presence of a memory function that is the sum of two exponentials with different memory times $\tau_1$ and $\tau_2$. Our Langevin simulations for the special case where both exponentials contribute equally to the total friction show that the barrier crossing time becomes independent of the longer memory time if at least one of the two memory times is larger than the intrinsic diffusion time. When we associate memory effects with coupled degrees of freedom that are orthogonal to a one-dimensional reaction coordinate, this counterintuitive result shows that the faster orthogonal degrees of freedom dominate barrier-crossing kinetics in the non-Markovian limit and that the slower orthogonal degrees become negligible, quite contrary to the standard time-scale separation assumption and with important consequences for the proper setup of coarse-graining procedures in the non-Markovian case. By asymptotic matching and symmetry arguments, we construct a crossover formula for the barrier crossing time that is valid for general multi-exponential memory kernels. This formula can be used to estimate barrier-crossing times for general memory functions for high friction, i.e. in the overdamped regime, as well as for low friction, i.e. in the inertial regime. Typical examples where our results are important include protein folding in the high-friction limit and chemical reactions such as proton-transfer reactions in the low-friction limit.

Based on the trajectories of the separation between water molecule pairs from MD simulations, we investigate the bond breakage dynamics in bulk water. From the spectrum of mean first-passage times, the Fokker-Planck equation allows us to derive the diffusivity profile along the separation coordinate and thus to unambiguously disentangle the effects of free-energy and local friction on the separation kinetics. For tightly coordinated water the friction is six times higher than in bulk, which can be interpreted in terms of a dominant reaction path that involves additional orthogonal coordinates., Comment: 11 pages, 10 figures

The concept of a protein diffusing in its free energy folding landscape has been fruitful for both theory and experiment. Yet the choice of the reaction coordinate (RC) introduces an undesirable degree of arbitrariness into the problem. We analyze extensive simulation data of an alpha-helix in explicit water solvent as it stochastically folds and unfolds. The free energy profiles for different RCs exhibit significant variation, some having an activation barrier, others not. We show that this variation has little effect on the predicted folding kinetics if the diffusivity profiles are properly taken into account. This kinetic quasi-universality is rationalized by an RC rescaling, which, due to the reparameterization invariance of the Fokker-Planck equation, allows the combination of free energy and diffusivity effects into a single function, the rescaled free energy profile. This rescaled free energy indeed shows less variation among different RCs than the bare free energy and diffusivity profiles separately do, if we properly distinguish between RCs that contain knowledge of the native state and those that are purely geometric in nature. Our method for extracting diffusivity profiles is easily applied to experimental single molecule time series data and might help to reconcile conflicts that arise when comparing results from different experimental probes for the same protein., Comment: 22 pages, 9 figures + supplementary material (8 pages, 5 figures); to appear in J. Chem. Phys.

Molecular dynamics simulations of ionic solutions depend sensitively on the force fields employed for the ions. To resolve the fine differences between ions of the same valence and roughly similar size and in particular to correctly describe ion-specific effects, it is clear that accurate force fields are necessary. In the past, optimization strategies for ionic force fields either considered single-ion properties (such as the solvation free energy at infinite dilution or the ion-water structure) or ion-pair properties (in the form of ion-ion distribution functions). In this paper we investigate strategies to optimize ionic force fields based on single-ion and ion-pair properties simultaneously. To that end, we simulate five different salt solutions, namely CsCl, KCl, NaI, KF, and CsI, at finite ion concentration. The force fields of these ions are systematically varied under the constraint that the single-ion solvation free energy matches the experimental value, which reduces the two-dimensional $\{\sigma,\epsilon\}$ parameter space of the Lennard Jones interaction to a one dimensional line for each ion. From the finite-concentration simulations, the pair-potential is extracted and the osmotic coefficient is calculated, which is compared to experimental data. We find a strong dependence of the osmotic coefficient on the force field, which is remarkable as the single-ion solvation free energy and the ion-water structure remain invariant under the parameter variation. Optimization of the force field is achieved for the cations Cs$^+$ and K$^+$, while for the anions I$^-$ and F$^-$ the experimental osmotic coefficient cannot be reached. This suggests that in the long run, additional parameters might have to be introduced into the modeling, for example by modified mixing rules., Comment: 29 pages, 6 figures, to appear in J. Chem. Phys

We investigate diffusion-limited reactions between a diffusing particle and a target site on a semiflexible polymer, a key factor determining the kinetics of DNA-protein binding and polymerization of cytoskeletal filaments. Our theory focuses on two competing effects: polymer shape fluctuations, which speed up association, and the hydrodynamic coupling between the diffusing particle and the chain, which slows down association. Polymer bending fluctuations are described using a mean field dynamical theory, while the hydrodynamic coupling between polymer and particle is incorporated through a simple heuristic approximation. Both of these we validate through comparison with Brownian dynamics simulations. Neither of the effects has been fully considered before in the biophysical context, and we show they are necessary to form accurate estimates of reaction processes. The association rate depends on the stiffness of the polymer and the particle size, exhibiting a maximum for intermediate persistence length and a minimum for intermediate particle radius. In the parameter range relevant to DNA-protein binding, the rate increase is up to 100% compared to the Smoluchowski result for simple center-of-mass motion. The quantitative predictions made by the theory can be tested experimentally., Comment: 21 pages, 11 figures, 1 table