Researcher profile

Roland R. Netz

Roland R. Netz contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
12topics
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

20 published item(s)

preprint2026arXiv

Barrier-crossing and energy relaxation dynamics of non-Markovian inertial systems connected via analytical Green-Fokker-Planck approach

From numerical simulations it is known that the barrier-crossing time of a non-Markovian one-dimensional reaction coordinate with a single exponentially decaying memory function exhibits a memory-turnover: for intermediate values of the memory decay time the barrier-crossing time is reduced compared to the Markovian limit and for long memory times increases quadratically with the memory time when keeping the total integrated friction and the mass constant. The intermediate memory acceleration regime is accurately predicted by Grote-Hynes theory, for the asymptotic long-memory slow-down behavior no systematic analytically tractable theory is available. Starting from the Green function for a general inertial (i.e. finite-mass) non-Markovian Gaussian reaction coordinate in a harmonic well, we derive by an exact mapping a generalized Fokker-Planck equation with a time-dependent effective diffusion constant. To first order in a systematic cumulant expansion we derive an analytical Arrhenius expression for the barrier-crossing time with the pre-exponential factor given by the energy relaxation time, which can be used to robustly predict barrier-crossing times from simulation or experimental trajectory data of general non-Markovian inertial systems without the need to extract memory functions. For a single exponential memory kernel we give a closed-form expression for the barrier-crossing time, which reproduces the Kramers turnover between the high-friction and high-mass limits as well as the memory turnover from the intermediate memory acceleration to the asymptotic long-memory slow-down regime. We also show that non-Markovian systems are singular in the zero-mass limit, which suggests that the long-memory barrier-crossing slow-down reflects the interplay between mass and memory effects. Thus, physically sound models for non-Markovian systems have to include a finite mass.

preprint2026arXiv

Beyond the Electric Dipole Approximation: Electric and Magnetic Multipole Contributions Reveal Biaxial Water Structure from SFG Spectra at the Air-Water Interface

The interpretation of sum-frequency-generation (SFG) spectra has been severely limited by the absence of quantitative theoretical predictions of higher-order multipole contributions. Magnetic dipole and electric quadrupole contributions are determined by bulk properties but appear in all experimental SFG spectra, obscuring the connection between measured spectra and interfacial structure. We present the simulation-based framework to predict the full set of multipole spectral contributions. This framework also yields depth-resolved spectra, enabling the precise spatial localization of spectroscopic features. Applied to the air-water interface, our approach achieves quantitative agreement with experimental spectra for different polarization combinations in both the bending and stretching regions. Higher-order multipole contributions are crucial for correctly interpreting SFG spectra: in the bending band, the electric dipole and the magnetic dipole contributions have similar intensities, while the electric quadrupole contribution is significantly larger. In the OH-stretch region, the electric quadrupole contribution is found to be in large part responsible for the characteristic shoulder at 3600/cm. Crucially, subtracting the quadrupole and magnetic contributions isolates the second-order electric dipole susceptibility, which is a quantitative probe for interfacial molecular orientational anisotropy. This electric-dipole susceptibility reveals a pronounced biaxial ordering of water at the air-water interface. By resolving a fundamental limitation of the interpretation of SFG spectroscopy, our framework allows for the detailed extraction of interfacial water ordering from SFG spectra.

preprint2026arXiv

Combined effects of evaporation, sedimentation and solute crystallization on the dynamics of aerosol size distributions on multiple length and time scales

We investigate three aspects of aerosol-mediated air-borne viral infection mechanisms on different length and time scales. First, we address the evolution of the size distribution of a non-interacting ensemble of droplets that are subject to evaporation and sedimentation using a sharp droplet-air interface model. From the exact solution of the evolution equation we derive the viral load in the air and show that it depends sensitively on the relative humidity. Secondly, from Molecular Dynamics simulations we extract the molecular reflection coefficient of single water molecules from the air-water interface. This parameter determines the water condensation and evaporation rate at a liquid droplet surface and therefore the evaporation rate of aqueous droplets. We find the reflection of water to be negligible at room temperature but to rise significantly at elevated temperatures and for grazing incidence angles. Thirdly, we derive a thermodynamically consistent three-dimensional diffuse-interface model for solute-containing droplets that is formulated as a three-phase Cahn-Hilliard/Allen-Cahn system. By numerically solving the coupled system of equations, we explore representative scenarios that show that this model reproduces and generalizes features of the sharp-interface model. These interconnected studies on the dynamics of aerosol droplet evaporation are relevant in order to quantitatively assess the airborne infection risk under varying environmental conditions.

preprint2022arXiv

Barrier-crossing times for different non-Markovian friction in well and barrier $-$ A numerical study

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.

preprint2022arXiv

Generalized Langevin Equation with a Non-Linear Potential of Mean Force and Non-Linear Memory Friction From a Hybrid Projection Scheme

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.

preprint2022arXiv

Spectral signatures of excess-proton waiting and transfer-path dynamics in aqueous hydrochloric acid solutions

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.

preprint2022arXiv

Time-dependent friction effects on vibrational infrared frequencies and line shapes of liquid water

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.

preprint2020arXiv

A rigorous derivation and energetics of a wave equation with fractional damping

We consider a linear system that consists of a linear wave equation on a horizontal hypersurface and a parabolic equation in the half space below. The model describes longitudinal elastic waves in organic monolayers at the water-air interface, which is an experimental setup that is relevant for understanding wave propagation in biological membranes. We study the scaling regime where the relevant horizontal length scale is much larger than the vertical length scale and provide a rigorous limit leading to a fractionally-damped wave equation for the membrane. We provide the associated existence results via linear semigroup theory and show convergence of the solutions in the scaling limit. Moreover, based on the energy-dissipation structure for the full model, we derive a natural energy and a natural dissipation function for the fractionally-damped wave equation with a time derivative of order 3/2

preprint2020arXiv

Approach to equilibrium and non-equilibrium stationary distributions of interacting many-particle systems that are coupled to different heat baths

A Hamiltonian-based model of many harmonically interacting massive particles that are subject to linear friction and coupled to heat baths at different temperatures is used to study the dynamic approach to equilibrium and non-equilibrium stationary states. Based on the exactly calculated dynamic approach to the stationary distribution, the functional that governs this approach, which is called the free entropy, is constructed. For the stationary distribution the free entropy becomes maximal and its time derivative is minimal and vanishes. Thus, the free entropy characterizes equilibrium as well as non-equilibrium stationary distributions by their extremal and stability properties. For an equilibrium system, i.e. if all heat baths have the same temperature, the free entropy equals the negative free energy divided by temperature. Using a systematic perturbative scheme for calculating velocity and position correlations in the overdamped massless limit, explicit results for few particles are presented: For two particles localization in position and momentum space is demonstrated in the non-equilibrium stationary state, indicative of a tendency to phase separate. For three elastically interacting particles heat flows from a particle coupled to a cold reservoir to a particle coupled to a warm reservoir if the third reservoir is sufficiently hot. Active particle models can be described in the same general framework, which thereby allows to characterize their entropy production not only in the stationary state but also in the approach to the stationary non-equilibrium state. Finally, the connection to non-equilibrium thermodynamics formulations that include the reservoir entropy production is discussed.

preprint2020arXiv

Emergence of molecular friction in liquids: bridging between the atomistic and hydrodynamic pictures

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.

preprint2020arXiv

Energy Transfer within the Hydrogen Bonding Network of Water Following Resonant Terahertz Excitation

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.

preprint2020arXiv

Exploring the Absorption Spectrum of Simulated Water from MHz to the Infrared

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.

preprint2020arXiv

Particle Diffusivity and Free-Energy Profiles in Inhomogeneous Hydrogel Systems from Time-Resolved Penetration Profiles

A combined experimental/theoretical method to simultaneously determine diffusivity and free-energy profiles of particles that penetrate into inhomogeneous hydrogel systems is presented. As the only input, arbitrarily normalized concentration profiles from fluorescence intensity data of labeled tracer particles for different penetration times are needed. The method is applied to dextran molecules of varying size which penetrate into hydrogels of polyethylene-glycol (PEG) chains with different lengths that are covalently cross-linked by hyperbranched polyglycerol (hPG) hubs. Extracted dextran bulk diffusivities agree well with fluorescence correlation spectroscopy data obtained separately. Scaling laws for dextran diffusivities and free energies inside the hydrogel are identified as a function of the dextran mass. An elastic free-volume model that includes dextran as well as PEG linker flexibility describes the repulsive dextran-hydrogel interaction free energy, which is of steric origin, quantitatively and furthermore suggests that the hydrogel mesh-size distribution is rather broad and particle penetration is dominated by large hydrogel pores. Particle penetration into hydrogels is for steric particle-hydrogel interactions thus suggested to be governed by an elastic size-filtering mechanism that involves the tail of the hydrogel pore-size distribution.

preprint2020arXiv

Viral air load due to sedimenting and evaporating droplets produced by speaking

The effect of evaporation on droplet sedimentation times is crucial for estimating the risk of infection from virus-containing airborne droplets. For droplet radii in the range 100 nm < R < 60 μm, evaporation can be described in the stagnant-flow approximation and is diffusion limited. Analytical equations are presented for the droplet evaporation rate, the time-dependent droplet size and the sedimentation time, including the significant effect of evaporation cooling. Evaporation makes the time for large droplets to sediment much longer and thus significantly increases the viral air load. Using recent estimates for SARS-CoV-2 concentrations in sputum and droplet production rates while speaking, a single infected person that constantly speaks without a mouth cover produces a total air load of more than 10^4 virions. In a mid-size closed room, this leads to a viral inhalation frequency of at least 2.5 per minute. Low relative humidity, as encountered inside buildings in winter and in airliners, speeds evaporation and thus keeps initially larger droplets suspended in air.

preprint2019arXiv

Non-Markovian barrier crossing with two-time-scale memory is dominated by the faster memory component

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 $τ_1$ and $τ_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.

preprint2018arXiv

Cold denaturation of RNA secondary structures with loop entropy and quenched disorder

We study the folding of RNA secondary structures with quenched sequence randomness by means of the constrained annealing method. A thermodynamic phase transition is induced by including the conformational weight of loop structures. In addition to the expected melting at high temperature, a cold melting transition appears. Our results suggest that the cold denaturation of RNA found experimentally is, in fact, a continuous phase transition triggered by quenched sequence disorder. We calculate both hot and cold melting critical temperatures for the competing energy scenario between favorable and unfavorable base pairs and present a phase diagram as a function of the loop exponent and temperature.

preprint2011arXiv

Anisotropic Hydrodynamic Mean-Field Theory for Semiflexible Polymers under Tension

We introduce an anisotropic mean-field approach for the dynamics of semiflexible polymers under intermediate tension, the force range where a chain is partially extended but not in the asymptotic regime of a nearly straight contour. The theory is designed to exactly reproduce the lowest order equilibrium averages of a stretched polymer, and treats the full complexity of the problem: the resulting dynamics include the coupled effects of long-range hydrodynamic interactions, backbone stiffness, and large-scale polymer contour fluctuations. Validated by Brownian hydrodynamics simulations and comparison to optical tweezer measurements on stretched DNA, the theory is highly accurate in the intermediate tension regime over a broad dynamical range, without the need for additional dynamic fitting parameters.

preprint2011arXiv

Deconvolution of dynamic mechanical networks

Time-resolved single-molecule biophysical experiments yield data that contain a wealth of dynamic information, in addition to the equilibrium distributions derived from histograms of the time series. In typical force spectroscopic setups the molecule is connected via linkers to a read-out device, forming a mechanically coupled dynamic network. Deconvolution of equilibrium distributions, filtering out the influence of the linkers, is a straightforward and common practice. We have developed an analogous dynamic deconvolution theory for the more challenging task of extracting kinetic properties of individual components in networks of arbitrary complexity and topology. Our method determines the intrinsic linear response functions of a given molecule in the network, describing the power spectrum of conformational fluctuations. The practicality of our approach is demonstrated for the particular case of a protein linked via DNA handles to two optically trapped beads at constant stretching force, which we mimic through Brownian dynamics simulations. Each well in the protein free energy landscape (corresponding to folded, unfolded, or possibly intermediate states) will have its own characteristic equilibrium fluctuations. The associated linear response function is rich in physical content, since it depends both on the shape of the well and its diffusivity---a measure of the internal friction arising from such processes like the transient breaking and reformation of bonds in the protein structure. Starting from the autocorrelation functions of the equilibrium bead fluctuations measured in this force clamp setup, we show how an experimentalist can accurately extract the state-dependent protein diffusivity using a straightforward two-step procedure.

preprint2011arXiv

Secondary structure formation of homopolymeric single-stranded nucleic acids including force and loop entropy: implications for DNA hybridization

Loops are essential secondary structure elements in folded DNA and RNA molecules and proliferate close to the melting transition. Using a theory for nucleic acid secondary structures that accounts for the logarithmic entropy c ln m for a loop of length m, we study homopolymeric single-stranded nucleic acid chains under external force and varying temperature. In the thermodynamic limit of a long strand, the chain displays a phase transition between a low temperature / low force compact (folded) structure and a high temperature / high force molten (unfolded) structure. The influence of c on phase diagrams, critical exponents, melting, and force extension curves is derived analytically. For vanishing pulling force, only for the limited range of loop exponents 2 < c < 2.479 a melting transition is possible; for c <= 2 the chain is always in the folded phase and for 2.479 < c always in the unfolded phase. A force induced melting transition with singular behavior is possible for all loop exponents c < 2.479 and can be observed experimentally by single molecule force spectroscopy. These findings have implications for the hybridization or denaturation of double stranded nucleic acids. The Poland-Scheraga model for nucleic acid duplex melting does not allow base pairing between nucleotides on the same strand in denatured regions of the double strand. If the sequence allows these intra-strand base pairs, we show that for a realistic loop exponent c ~ 2.1 pronounced secondary structures appear inside the single strands. This leads to a lower melting temperature of the duplex than predicted by the Poland-Scheraga model. Further, these secondary structures renormalize the effective loop exponent c^, which characterizes the weight of a denatured region of the double strand, and thus affect universal aspects of the duplex melting transition.

preprint2010arXiv

DNA-Protein Binding Rates: Bending Fluctuation and Hydrodynamic Coupling Effects

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.