@misc{Metzler2020, author = {Metzler, Ralf}, title = {Superstatistics and non-Gaussian diffusion}, series = {The European physical journal special topics}, volume = {229}, journal = {The European physical journal special topics}, number = {5}, publisher = {Springer}, address = {Heidelberg}, issn = {1951-6355}, doi = {10.1140/epjst/e2020-900210-x}, pages = {711 -- 728}, year = {2020}, abstract = {Brownian motion and viscoelastic anomalous diffusion in homogeneous environments are intrinsically Gaussian processes. In a growing number of systems, however, non-Gaussian displacement distributions of these processes are being reported. The physical cause of the non-Gaussianity is typically seen in different forms of disorder. These include, for instance, imperfect "ensembles" of tracer particles, the presence of local variations of the tracer mobility in heteroegenous environments, or cases in which the speed or persistence of moving nematodes or cells are distributed. From a theoretical point of view stochastic descriptions based on distributed ("superstatistical") transport coefficients as well as time-dependent generalisations based on stochastic transport parameters with built-in finite correlation time are invoked. After a brief review of the history of Brownian motion and the famed Gaussian displacement distribution, we here provide a brief introduction to the phenomenon of non-Gaussianity and the stochastic modelling in terms of superstatistical and diffusing-diffusivity approaches.}, language = {en} } @article{KumarHesseRaoetal.2020, author = {Kumar, Rohini and Hesse, Fabienne and Rao, P. Srinivasa and Musolff, Andreas and Jawitz, James and Sarrazin, Francois and Samaniego, Luis and Fleckenstein, Jan H. and Rakovec, Oldrich and Thober, S. and Attinger, Sabine}, title = {Strong hydroclimatic controls on vulnerability to subsurface nitrate contamination across Europe}, series = {Nature Communications}, volume = {11}, journal = {Nature Communications}, number = {1}, publisher = {Nature Publishing Group UK}, address = {London}, issn = {2041-1723}, doi = {10.1038/s41467-020-19955-8}, pages = {1 -- 10}, year = {2020}, abstract = {Subsurface contamination due to excessive nutrient surpluses is a persistent and widespread problem in agricultural areas across Europe. The vulnerability of a particular location to pollution from reactive solutes, such as nitrate, is determined by the interplay between hydrologic transport and biogeochemical transformations. Current studies on the controls of subsurface vulnerability do not consider the transient behaviour of transport dynamics in the root zone. Here, using state-of-the-art hydrologic simulations driven by observed hydroclimatic forcing, we demonstrate the strong spatiotemporal heterogeneity of hydrologic transport dynamics and reveal that these dynamics are primarily controlled by the hydroclimatic gradient of the aridity index across Europe. Contrasting the space-time dynamics of transport times with reactive timescales of denitrification in soil indicate that similar to 75\% of the cultivated areas across Europe are potentially vulnerable to nitrate leaching for at least onethird of the year. We find that neglecting the transient nature of transport and reaction timescale results in a great underestimation of the extent of vulnerable regions by almost 50\%. Therefore, future vulnerability and risk assessment studies must account for the transient behaviour of transport and biogeochemical transformation processes.}, language = {en} } @article{JeonChechkinMetzler2014, author = {Jeon, Jae-Hyung and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Scaled Brownian motion: a paradoxical process with a time dependent diffusivity for the description of anomalous diffusion}, series = {Physical chemistry, chemical physics : PCCP}, volume = {30}, journal = {Physical chemistry, chemical physics : PCCP}, number = {16}, publisher = {The Royal Society of Chemistry}, address = {Cambridge}, doi = {10.1039/C4CP02019G}, pages = {15811 -- 15817}, year = {2014}, abstract = {Anomalous diffusion is frequently described by scaled Brownian motion (SBM){,} a Gaussian process with a power-law time dependent diffusion coefficient. Its mean squared displacement is ?x2(t)? [similar{,} equals] 2K(t)t with K(t) [similar{,} equals] t[small alpha]-1 for 0 < [small alpha] < 2. SBM may provide a seemingly adequate description in the case of unbounded diffusion{,} for which its probability density function coincides with that of fractional Brownian motion. Here we show that free SBM is weakly non-ergodic but does not exhibit a significant amplitude scatter of the time averaged mean squared displacement. More severely{,} we demonstrate that under confinement{,} the dynamics encoded by SBM is fundamentally different from both fractional Brownian motion and continuous time random walks. SBM is highly non-stationary and cannot provide a physical description for particles in a thermalised stationary system. Our findings have direct impact on the modelling of single particle tracking experiments{,} in particular{,} under confinement inside cellular compartments or when optical tweezers tracking methods are used.}, language = {en} } @article{AndersSaitHorsley2022, author = {Anders, Janet and Sait, Connor R. J. and Horsley, Simon A. R.}, title = {Quantum Brownian motion for magnets}, series = {New journal of physics : the open-access journal for physics}, volume = {24}, journal = {New journal of physics : the open-access journal for physics}, number = {3}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/ac4ef2}, pages = {21}, year = {2022}, abstract = {Spin precession in magnetic materials is commonly modelled with the classical phenomenological Landau-Lifshitz-Gilbert (LLG) equation. Based on a quantized three-dimensional spin + environment Hamiltonian, we here derive a spin operator equation of motion that describes precession and includes a general form of damping that consistently accounts for memory, coloured noise and quantum statistics. The LLG equation is recovered as its classical, Ohmic approximation. We further introduce resonant Lorentzian system-reservoir couplings that allow a systematic comparison of dynamics between Ohmic and non-Ohmic regimes. Finally, we simulate the full non-Markovian dynamics of a spin in the semi-classical limit. At low temperatures, our numerical results demonstrate a characteristic reduction and flattening of the steady state spin alignment with an external field, caused by the quantum statistics of the environment. The results provide a powerful framework to explore general three-dimensional dissipation in quantum thermodynamics.}, language = {en} } @article{JayVazdaCruzEckertetal.2020, author = {Jay, Raphael M. and Vaz da Cruz, Vinicius and Eckert, Sebastian and Fondell, Mattis and Mitzner, Rolf and F{\"o}hlisch, Alexander}, title = {Probing solute-solvent interactions of transition metal complexes using L-edge absorption spectroscopy}, series = {The journal of physical chemistry : B, Condensed matter, materials, surfaces, interfaces \& biophysical chemistry}, volume = {124}, journal = {The journal of physical chemistry : B, Condensed matter, materials, surfaces, interfaces \& biophysical chemistry}, number = {27}, publisher = {American Chemical Society}, address = {Washington}, issn = {1520-6106}, doi = {10.1021/acs.jpcb.0c00638}, pages = {5636 -- 5645}, year = {2020}, abstract = {In order to tailor solution-phase chemical reactions involving transition metal complexes, it is critical to understand how their valence electronic charge distributions are affected by the solution environment. Here, solute-solvent interactions of a solvatochromic mixed-ligand iron complex were investigated using X-ray absorption spectroscopy at the transition metal L-2,L-3-edge. Due to the selectivity of the corresponding core excitations to the iron 3d orbitals, the method grants direct access to the valence electronic structure around the iron center and its response to interactions with the solvent environment. A linear increase of the total L-2,L-3-edge absorption cross section as a function of the solvent Lewis acidity is revealed. The effect is caused by relative changes in different metal-ligand-bonding channels, which preserve local charge densities while increasing the density of unoccupied states around the iron center. These conclusions are corroborated by a combination of molecular dynamics and spectrum simulations based on time-dependent density functional theory. The simulations reproduce the spectral trends observed in the X-ray but also optical absorption experiments. Our results underscore the importance of solute-solvent interactions when aiming for an accurate description of the valence electronic structure of solvated transition metal complexes and demonstrate how L-2,L-3-edge absorption spectroscopy can aid in understanding the impact of the solution environment on intramolecular covalency and the electronic charge distribution.}, language = {en} } @article{vanVelzen2020, author = {van Velzen, Ellen}, title = {Predator coexistence through emergent fitness equalization}, series = {Ecology : a publication of the Ecological Society of America}, volume = {101}, journal = {Ecology : a publication of the Ecological Society of America}, number = {5}, publisher = {Wiley}, address = {Hoboken}, issn = {0012-9658}, doi = {10.1002/ecy.2995}, pages = {10}, year = {2020}, abstract = {The competitive exclusion principle is one of the oldest ideas in ecology and states that without additional self-limitation two predators cannot coexist on a single prey. The search for mechanisms allowing coexistence despite this has identified niche differentiation between predators as crucial: without this, coexistence requires the predators to have exactly the same R* values, which is considered impossible. However, this reasoning misses a critical point: predators' R* values are not static properties, but affected by defensive traits of their prey, which in turn can adapt in response to changes in predator densities. Here I show that this feedback between defense and predator dynamics enables stable predator coexistence without ecological niche differentiation. Instead, the mechanism driving coexistence is that prey adaptation causes defense to converge to the value where both predators have equal R* values ("fitness equalization"). This result is highly general, independent of specific model details, and applies to both rapid defense evolution and inducible defenses. It demonstrates the importance of considering long-standing ecological questions from an eco-evolutionary viewpoint, and showcases how the effects of adaptation can cascade through communities, driving diversity on higher trophic levels. These insights offer an important new perspective on coexistence theory.}, language = {en} } @article{Omelʹchenko2020, author = {Omelʹchenko, Oleh E.}, title = {Nonstationary coherence-incoherence patterns in nonlocally coupled heterogeneous phase oscillators}, series = {Chaos : an interdisciplinary journal of nonlinear science}, volume = {30}, journal = {Chaos : an interdisciplinary journal of nonlinear science}, number = {4}, publisher = {American Institute of Physics}, address = {Melville}, issn = {1054-1500}, doi = {10.1063/1.5145259}, pages = {8}, year = {2020}, abstract = {We consider a large ring of nonlocally coupled phase oscillators and show that apart from stationary chimera states, this system also supports nonstationary coherence-incoherence patterns (CIPs). For identical oscillators, these CIPs behave as breathing chimera states and are found in a relatively small parameter region only. It turns out that the stability region of these states enlarges dramatically if a certain amount of spatially uniform heterogeneity (e.g., Lorentzian distribution of natural frequencies) is introduced in the system. In this case, nonstationary CIPs can be studied as stable quasiperiodic solutions of a corresponding mean-field equation, formally describing the infinite system limit. Carrying out direct numerical simulations of the mean-field equation, we find different types of nonstationary CIPs with pulsing and/or alternating chimera-like behavior. Moreover, we reveal a complex bifurcation scenario underlying the transformation of these CIPs into each other. These theoretical predictions are confirmed by numerical simulations of the original coupled oscillator system.}, language = {en} } @article{ShinCherstvyMetzler2014, author = {Shin, Jaeoh and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Kinetics of polymer looping with macromolecular crowding: effects of volume fraction and crowder size}, series = {Soft Matter}, journal = {Soft Matter}, editor = {Metzler, Ralf}, publisher = {The Royal Society of Chemistry}, address = {Cambridge}, issn = {1744-683X}, pages = {472 -- 488}, year = {2014}, abstract = {The looping of polymers such as DNA is a fundamental process in the molecular biology of living cells, whose interior is characterised by a high degree of molecular crowding. We here investigate in detail the looping dynamics of flexible polymer chains in the presence of different degrees of crowding. From the analysis of the looping-unlooping rates and the looping probabilities of the chain ends we show that the presence of small crowders typically slows down the chain dynamics but larger crowders may in fact facilitate the looping. We rationalise these non-trivial and often counterintuitive effects of the crowder size on the looping kinetics in terms of an effective solution viscosity and standard excluded volume. It is shown that for small crowders the effect of an increased viscosity dominates, while for big crowders we argue that confinement effects (caging) prevail. The tradeoff between both trends can thus result in the impediment or facilitation of polymer looping, depending on the crowder size. We also examine how the crowding volume fraction, chain length, and the attraction strength of the contact groups of the polymer chain affect the looping kinetics and hairpin formation dynamics. Our results are relevant for DNA looping in the absence and presence of protein mediation, DNA hairpin formation, RNA folding, and the folding of polypeptide chains under biologically relevant high-crowding conditions.}, language = {en} } @article{ScherlerSchwanghart2020, author = {Scherler, Dirk and Schwanghart, Wolfgang}, title = {Drainage divide networks}, series = {Earth surface dynamics}, volume = {8}, journal = {Earth surface dynamics}, number = {2}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {2196-6311}, doi = {10.5194/esurf-8-261-2020}, pages = {261 -- 274}, year = {2020}, abstract = {Drainage divides are organized into tree-like networks that may record information about drainage divide mobility. However, views diverge about how to best assess divide mobility. Here, we apply a new approach of automatically extracting and ordering drainage divide networks from digital elevation models to results from landscape evolution model experiments. We compared landscapes perturbed by strike-slip faulting and spatiotemporal variations in erodibility to a reference model to assess which topographic metrics (hillslope relief, flow distance, and chi) are diagnostic of divide mobility. Results show that divide segments that are a minimum distance of similar to 5 km from river confluences strive to attain constant values of hillslope relief and flow distance to the nearest stream. Disruptions of such patterns can be related to mobile divides that are lower than stable divides, closer to streams, and often asymmetric in shape. In general, we observe that drainage divides high up in the network, i.e., at great distances from river confluences, are more susceptible to disruptions than divides closer to these confluences and are thus more likely to record disturbance for a longer time period. We found that across-divide differences in hillslope relief proved more useful for assessing divide migration than other tested metrics. However, even stable drainage divide networks exhibit across-divide differences in any of the studied topographic metrics. Finally, we propose a new metric to quantify the connectivity of divide junctions.}, language = {en} } @article{EhrlichGaedke2020, author = {Ehrlich, Elias and Gaedke, Ursula}, title = {Coupled changes in traits and biomasses cascading through a tritrophic plankton food web}, series = {Limnology and oceanography}, volume = {65}, journal = {Limnology and oceanography}, number = {10}, publisher = {Wiley}, address = {Hoboken}, issn = {0024-3590}, doi = {10.1002/lno.11466}, pages = {2502 -- 2514}, year = {2020}, abstract = {Trait-based approaches have broadened our understanding of how the composition of ecological communities responds to environmental drivers. This research has mainly focussed on abiotic factors and competition determining the community trait distribution, while effects of trophic interactions on trait dynamics, if considered at all, have been studied for two trophic levels at maximum. However, natural food webs are typically at least tritrophic. This enables indirect interactions of traits and biomasses among multiple trophic levels leading to underexplored effects on food web dynamics. Here, we demonstrate the occurrence of mutual trait adjustment among three trophic levels in a natural plankton food web (Lake Constance) and in a corresponding mathematical model. We found highly recurrent seasonal biomass and trait dynamics, where herbivorous zooplankton increased its size, and thus its ability to counter phytoplankton defense, before phytoplankton defense actually increased. This is contrary to predictions from bitrophic systems where counter-defense of the consumer is a reaction to prey defense. In contrast, counter-defense of carnivores by size adjustment followed the defense of herbivores as expected. By combining observations and model simulations, we show how the reversed trait dynamics at the two lower trophic levels result from a "trophic biomass-trait cascade" driven by the carnivores. Trait adjustment between two trophic levels can therefore be altered by biomass or trait changes of adjacent trophic levels. Hence, analyses of only pairwise trait adjustment can be misleading in natural food webs, while multitrophic trait-based approaches capture indirect biomass-trait interactions among multiple trophic levels.}, language = {en} } @article{MarschallSkorovZakharovetal.2020, author = {Marschall, Raphael and Skorov, Yuri and Zakharov, Vladimir and Rezac, Ladislav and Gerig, Selina-Barbara and Christou, Chariton and Dadzie, S. Kokou and Migliorini, Alessandra and Rinaldi, Giovanna and Agarwal, Jessica and Vincent, Jean-Baptiste and Kappel, David}, title = {Cometary comae-surface links the physics of gas and dust from the surface to a spacecraft}, series = {Space science reviews}, volume = {216}, journal = {Space science reviews}, number = {8}, publisher = {Springer}, address = {Dordrecht}, issn = {0038-6308}, doi = {10.1007/s11214-020-00744-0}, pages = {53}, year = {2020}, abstract = {A comet is a highly dynamic object, undergoing a permanent state of change. These changes have to be carefully classified and considered according to their intrinsic temporal and spatial scales. The Rosetta mission has, through its contiguous in-situ and remote sensing coverage of comet 67P/Churyumov-Gerasimenko (hereafter 67P) over the time span of August 2014 to September 2016, monitored the emergence, culmination, and winding down of the gas and dust comae. This provided an unprecedented data set and has spurred a large effort to connect in-situ and remote sensing measurements to the surface. In this review, we address our current understanding of cometary activity and the challenges involved when linking comae data to the surface. We give the current state of research by describing what we know about the physical processes involved from the surface to a few tens of kilometres above it with respect to the gas and dust emission from cometary nuclei. Further, we describe how complex multidimensional cometary gas and dust models have developed from the Halley encounter of 1986 to today. This includes the study of inhomogeneous outgassing and determination of the gas and dust production rates. Additionally, the different approaches used and results obtained to link coma data to the surface will be discussed. We discuss forward and inversion models and we describe the limitations of the respective approaches. The current literature suggests that there does not seem to be a single uniform process behind cometary activity. Rather, activity seems to be the consequence of a variety of erosion processes, including the sublimation of both water ice and more volatile material, but possibly also more exotic processes such as fracture and cliff erosion under thermal and mechanical stress, sub-surface heat storage, and a complex interplay of these processes. Seasons and the nucleus shape are key factors for the distribution and temporal evolution of activity and imply that the heliocentric evolution of activity can be highly individual for every comet, and generalisations can be misleading.}, language = {en} } @article{BanerjeeLipowskySanter2020, author = {Banerjee, Pallavi and Lipowsky, Reinhard and Santer, Mark}, title = {Coarse-grained molecular model for the Glycosylphosphatidylinositol anchor with and without protein}, series = {Journal of Chemical Theory and Computation}, volume = {16}, journal = {Journal of Chemical Theory and Computation}, number = {6}, publisher = {ACS Publications}, address = {Washington DC}, issn = {1549-9626}, doi = {10.1021/acs.jctc.0c00056}, pages = {15}, year = {2020}, abstract = {Glycosylphosphatidylinositol (GPI) anchors are a unique class of complex glycolipids that anchor a great variety of proteins to the extracellular leaflet of plasma membranes of eukaryotic cells. These anchors can exist either with or without an attached protein called GPI-anchored protein (GPI-AP) both in vitro and in vivo. Although GPIs are known to participate in a broad range of cellular functions, it is to a large extent unknown how these are related to GPI structure and composition. Their conformational flexibility and microheterogeneity make it difficult to study them experimentally. Simplified atomistic models are amenable to all-atom computer simulations in small lipid bilayer patches but not suitable for studying their partitioning and trafficking in complex and heterogeneous membranes. Here, we present a coarse-grained model of the GPI anchor constructed with a modified version of the MARTINI force field that is suited for modeling carbohydrates, proteins, and lipids in an aqueous environment using MARTINI's polarizable water. The nonbonded interactions for sugars were reparametrized by calculating their partitioning free energies between polar and apolar phases. In addition, sugar-sugar interactions were optimized by adjusting the second virial coefficients of osmotic pressures for solutions of glucose, sucrose, and trehalose to match with experimental data. With respect to the conformational dynamics of GPI-anchored green fluorescent protein, the accessible time scales are now at least an order of magnitude larger than for the all-atom system. This is particularly important for fine-tuning the mutual interactions of lipids, carbohydrates, and amino acids when comparing to experimental results. We discuss the prospective use of the coarse-grained GPI model for studying protein-sorting and trafficking in membrane models.}, language = {en} } @article{JingKumarHesseetal.2020, author = {Jing, Miao and Kumar, Rohini and Heße, Falk and Thober, Stephan and Rakovec, Oldrich and Samaniego, Luis and Attinger, Sabine}, title = {Assessing the response of groundwater quantity and travel time distribution to 1.5, 2, and 3 °C global warming in a mesoscale central German basin}, series = {Hydrology and Earth System Sciences}, volume = {24}, journal = {Hydrology and Earth System Sciences}, number = {3}, publisher = {Copernicus Publ.}, address = {G{\"o}ttingen}, issn = {1607-7938}, doi = {10.5194/hess-24-1511-2020}, pages = {1511 -- 1526}, year = {2020}, abstract = {Groundwater is the biggest single source of high-quality freshwater worldwide, which is also continuously threatened by the changing climate. In this paper, we investigate the response of the regional groundwater system to climate change under three global warming levels (1.5, 2, and 3 ∘C) in a central German basin (N{\"a}gelstedt). This investigation is conducted by deploying an integrated modeling workflow that consists of a mesoscale hydrologic model (mHM) and a fully distributed groundwater model, OpenGeoSys (OGS). mHM is forced with climate simulations of five general circulation models under three representative concentration pathways. The diffuse recharges estimated by mHM are used as boundary forcings to the OGS groundwater model to compute changes in groundwater levels and travel time distributions. Simulation results indicate that groundwater recharges and levels are expected to increase slightly under future climate scenarios. Meanwhile, the mean travel time is expected to decrease compared to the historical average. However, the ensemble simulations do not all agree on the sign of relative change. Changes in mean travel time exhibit a larger variability than those in groundwater levels. The ensemble simulations do not show a systematic relationship between the projected change (in both groundwater levels and travel times) and the warming level, but they indicate an increased variability in projected changes with adjusting the enhanced warming level from 1.5 to 3 ∘C. Correspondingly, it is highly recommended to restrain the trend of global warming.}, language = {en} }