Source author record

Bernhard Mehlig

Bernhard Mehlig appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2024arXiv

Bayesian posterior approximation with stochastic ensembles

We introduce ensembles of stochastic neural networks to approximate the Bayesian posterior, combining stochastic methods such as dropout with deep ensembles. The stochastic ensembles are formulated as families of distributions and trained to approximate the Bayesian posterior with variational inference. We implement stochastic ensembles based on Monte Carlo dropout, DropConnect and a novel non-parametric version of dropout and evaluate them on a toy problem and CIFAR image classification. For both tasks, we test the quality of the posteriors directly against Hamiltonian Monte Carlo simulations. Our results show that stochastic ensembles provide more accurate posterior estimates than other popular baselines for Bayesian inference.

preprint2022arXiv

Collision fluctuations of lucky droplets with superdroplets

It was previously shown that the superdroplet algorithm for modeling the collision-coalescence process can faithfully represent mean droplet growth in turbulent clouds. But an open question is how accurately the superdroplet algorithm accounts for fluctuations in the collisional aggregation process. Such fluctuations are particularly important in dilute suspensions. Even in the absence of turbulence, Poisson fluctuations of collision times in dilute suspensions may result in substantial variations in the growth process, resulting in a broad distribution of growth times to reach a certain droplet size. We quantify the accuracy of the superdroplet algorithm in describing the fluctuating growth history of a larger droplet that settles under the effect of gravity in a quiescent fluid and collides with a dilute suspension of smaller droplets that were initially randomly distributed in space ('lucky droplet model'). We assess the effect of fluctuations upon the growth history of the lucky droplet and compute the distribution of cumulative collision times. The latter is shown to be sensitive enough to detect the subtle increase of fluctuations associated with collisions between multiple lucky droplets. The superdroplet algorithm incorporates fluctuations in two distinct ways: through the random spatial distribution of superdroplets and through the Monte Carlo collision algorithm involved. Using specifically designed numerical experiments, we show that both on their own give an accurate representation of fluctuations. We conclude that the superdroplet algorithm can faithfully represent fluctuations in the coagulation of droplets driven by gravity.

preprint2020arXiv

Condensational and collisional growth of cloud droplets in a turbulent environment

We investigate the effect of turbulence on the combined condensational and collisional growth of cloud droplets by means of high resolution direct numerical simulations of turbulence and a superparticle approximation for droplet dynamics and collisions. The droplets are subject to turbulence as well as gravity, and their collision and coalescence efficiencies are taken to be unity. We solve the thermodynamic equations governing temperature, water-vapor mixing ratio, and the resulting supersaturation fields together with the Navier-Stokes equation. We find that the droplet-size distribution broadens with increasing Reynolds number and/or mean energy dissipation rate. Turbulence affects the condensational growth directly through supersaturation fluctuations, and it influences collisional growth indirectly through condensation. Our simulations show for the first time that, in the absence of the mean updraft cooling, supersaturation fluctuation-induced broadening of droplet-size distributions enhances the collisional growth. This is contrary to classical (non-turbulent) condensational growth, which leads to a growing mean droplet size, but a narrower droplet-size distribution. Our findings, instead, show that condensational growth facilitates collisional growth by broadening the size distribution in the tails at an early stage of rain formation. With increasing Reynolds numbers, evaporation becomes stronger. This counteracts the broadening effect due to condensation at late stages of rain formation. Our conclusions are consistent with results of laboratory experiments and field observations, and show that supersaturation fluctuations are important for precipitation.

preprint2020arXiv

Intercellular communication induces glycolytic synchronization waves between individually oscillating cells

Metabolic oscillations in single cells underlie the mechanisms behind cell synchronization and cell-cell communication. For example, glycolytic oscillations mediated by biochemical communication between cells may synchronize the pulsatile insulin secretion by pancreatic tissue, and a link between glycolytic synchronization anomalies and type-2 diabetes has been hypotesized. Cultures of yeast cells have provided an ideal model system to study synchronization and propagation waves of glycolytic oscillations in large populations. However, the mechanism by which synchronization occurs at individual cell-cell level and overcome local chemical concentrations and heterogenic biological clocks, is still an open question because of experimental limitations in sensitive and specific handling of single cells. Here, we show how the coupling of intercellular diffusion with the phase regulation of individual oscillating cells induce glycolytic synchronization waves. We directly measure the single-cell metabolic responses from yeast cells in a microfluidic environment and characterize a discretized cell-cell communication using graph theory. We corroborate our findings with simulations based on a kinetic detailed model for individual yeast cells. These findings can provide insight into the roles cellular synchronization play in biomedical applications, such as insulin secretion regulation at the cellular level.

preprint2020arXiv

Stochastic unfolding of nanoconfined DNA: experiments, model and Bayesian analysis

Nanochannels provide means for detailed experiments on the effect of confinement on biomacromolecules, such as DNA. We here introduce a model for the complete unfolding of DNA from the circular to linear configuration. Two main ingredients are the entropic unfolding force as well as the friction coefficient for the unfolding process, and we describe the associated dynamics by a non-linear Langevin equation. By analyzing experimental data where DNA molecules are photo-cut and unfolded inside a nanochannel, our model allows us to extract values for the unfolding force as well as the friction coefficient for the first time. In order to extract numerical values for these physical quantities, we employ a recently introduced Bayesian inference framework. We find that the determined unfolding force is in agreement with estimates from a simple Flory type argument. The estimated friction coefficient is in agreement with theoretical estimates for motion of a cylinder in a channel. We further validate the estimated friction constant by extracting this parameter from DNA's center-of-mass motion before and after unfolding, yielding decent agreement.

preprint2019arXiv

Importance of fluid inertia for the orientation of spheroids settling in turbulent flow

How non-spherical particles orient as they settle in a flow has important practical implications in a number of scientific and engineering problems. In a quiescent fluid, a slowly settling particle orients so that it settles with its broad side first. This is an effect of the torque due to convective inertia of the fluid set in motion by the settling particle, which maximises the drag experienced by the particle. Turbulent flows tend to randomise the particle orientation. Recently the settling of non-spherical particles in turbulence was analysed neglecting the effect of convective fluid inertia, but taking into account the effect of the turbulent fluid-velocity gradients on the particle orientation. These studies reached the opposite conclusion, namely that a rod settles preferentially with its tip first, wheras a disk settles with its edge first, therefore minimizing the drag on the particle. Here, we consider both effects, the convective inertial torque as well as the torque due to fluctuating velocity gradients, and ask under which circumstances either one or the other dominate. To this end we estimate the ratio of the magnitudes of the two torques. Our estimates suggest that the fluid-inertia torque prevails in high-Reynolds number flows. In this case non-spherical particles are expected to settle with a maximal drag. But when the Reynolds number is small then the torque due to fluid-velocity gradients may dominate, causing the particle to settle with its broad side first.

preprint2015arXiv

Analysis of the Correlation Dimension for Inertial Particles

We obtain an implicit equation for the correlation dimension which describes clustering of inertial particles in a complex flow onto a fractal measure. Our general equation involves a propagator of a nonlinear stochastic process in which the velocity gradient of the fluid appears as additive noise. When the long-time limit of the propagator is considered our equation reduces to an existing large-deviation formalism, from which it is difficult to extract concrete results. In the short-time limit, however, our equation reduces to a solvability condition on a partial differential equation. We show how this approach leads to a perturbative expansion of the correlation dimension, for which the coefficients can be obtained exactly and in principle to any order. We derive the perturbation series for the correlation dimension of inertial particles suspended in three-dimensional spatially smooth random flows with white-noise time correlations, obtaining the first 33 non-zero coefficients exactly.

preprint2015arXiv

How nanochannel confinement affects the DNA melting transition within the Poland-Scheraga model

When double-stranded DNA molecules are heated, or exposed to denaturing agents, the two strands get separated. The statistical physics of this process has a long history, and is commonly described in term of the Poland-Scheraga (PS) model. Crucial to this model is the configurational entropy for a melted region (compared to the entropy of an intact region of the same size), quantified by the loop factor. In this study we investigate how confinement affects the DNA melting transition, by using the loop factor for an ideal Gaussian chain. By subsequent numerical solutions of the PS model, we demonstrate that the melting temperature depends on the persistence lengths of single-stranded and double-stranded DNA. For realistic values of the persistence lengths the melting temperature is predicted to decrease with decreasing channel diameter. We also demonstrate that confinement broadens the melting transition. These general findings hold for the three scenarios investigated: namely 1. homo-DNA, i.e. identical basepairs along the DNA molecule; 2. random sequence DNA, and 3. "real" DNA, here T4 phage DNA. We show that cases 2 and 3 in general give rise to broader transitions than case 1. Case 3 exhibits a similar phase transition as case 2 provided the random sequence DNA has the same ratio of AT to GC basepairs. A simple analytical estimate for the shift in melting temperature is provided as a function of nanochannel diameter. For homo-DNA, we also present an analytical prediction of the melting probability as a function of temperature.

preprint2015arXiv

Measurement of particle motion in optical tweezers embedded in a Sagnac interferometer

We have constructed a counterpropagating optical tweezers setup embedded in a Sagnac interferometer in order to increase the sensitivity of position tracking for particles in the geometrical optics regime. Enhanced position determination using a Sagnac interferometer has previously been described theoretically by Taylor et al. [Journal of Optics 13, 044014 (2011)] for Rayleigh-regime particles trapped in an antinode of a standing wave. We have extended their theory to a case of arbitrarily-sized particles trapped with orthogonally-polarized counterpropagating beams. The working distance of the setup was sufficiently long to optically induce particle oscillations orthogonally to the axis of the tweezers with an auxiliary laser beam. Using these oscillations as a reference, we have experimentally shown that Sagnac-enhanced back focal plane interferometry is capable of providing an improvement of more than 5 times in the signal-to-background ratio, corresponding to a more than 30-fold improvement of the signal-to-noise ratio. The experimental results obtained are consistent with our theoretical predictions. In the experimental setup, we used a method of optical levitator-assisted liquid droplet delivery in air based on commercial inkjet technology, with a novel method to precisely control the size of droplets.

preprint2013arXiv

Demography-adjusted tests of neutrality based on genome-wide SNP data

Tests of the neutral evolution hypothesis are usually built on the standard null model which assumes that mutations are neutral and population size remains constant over time. However, it is unclear how such tests are affected if the last assumption is dropped. Here, we extend the unifying framework for tests based on the site frequency spectrum, introduced by Achaz and Ferretti, to populations of varying size. A key ingredient is to specify the first two moments of the frequency spectrum. We show that these moments can be determined analytically if a population has experienced two instantaneous size changes in the past. We apply our method to data from ten human populations gathered in the 1000 genomes project, estimate their demographies and define demography-adjusted versions of Tajima's $D$, Fay & Wu's $H$, and Zeng's $E$. The adjusted test statistics facilitate the direct comparison between populations and they show that most of the differences among populations seen in the original tests can be explained by demography. We carried out whole genome screens for deviation from neutrality and identified candidate regions of recent positive selection. We provide track files with values of the adjusted and original tests for upload to the UCSC genome browser.

preprint2012arXiv

Universal anomalous diffusion of weakly damped particles

We show that anomalous diffusion arises in two different models for the motion of randomly forced and weakly damped particles: one is a generalisation of the Ornstein-Uhlenbeck process with a random force which depends on position as well as time, the other is a generalisation of the Chandrasekhar-Rosenbluth model of stellar dynamics, encompassing non-Coulombic potentials. We show that both models exhibit anomalous diffusion of position $x$ and momentum $p$ with the same exponents: $<x^2> \sim C_x t^2$ and $<p^2> \sim C_p t^{2/5}$. We are able to determine the prefactors $C_x$, $C_p$ analytically.

preprint2009arXiv

Correlation Dimension of Inertial Particles in Random Flows

We obtain an implicit equation for the correlation dimension of dynamical systems in terms of an integral over a propagator. We illustrate the utility of this approach by evaluating the correlation dimension for inertial particles suspended in a random flow. In the limit where the correlation time of the flow field approaches zero, taking the short-time limit of the propagator enables the correlation dimension to be determined from the solution of a partial differential equation. We develop the solution as a power series in a dimensionless parameter which represents the strength of inertial effects.

preprint2009arXiv

Poincare Indices of Rheoscopic Visualisations

Suspensions of small anisotropic particles, termed 'rheoscopic fluids', are used for flow visualisation. By illuminating the fluid with light of three different colours, it is possible to determine Poincare indices for vector fields formed by the longest axis of the particles. Because this vector field is non-oriented, half-integer Poincare indices are possible, and are observed experimentally. An exact solution for the direction vector appears to preclude the existence of topological singularities. However, we show that upon averaging over the random initial orientations of particles, singularities with half-integer Poincare index appear. We describe their normal forms.

preprint2008arXiv

Fingerprints of Random Flows?

We consider the patterns formed by small rod-like objects advected by a random flow in two dimensions. An exact solution indicates that their direction field is non-singular. However, we find from simulations that the direction field of the rods does appear to exhibit singularities. First, ` scar lines' emerge where the rods abruptly change direction by $π$. Later, these scar lines become so narrow that they ` heal over' and disappear, but their ends remain as point singularities, which are of the same type as those seen in fingerprints. We give a theoretical explanation for these observations.

preprint2008arXiv

Planet Formation by Concurrent Collapse

After reviewing the difficulties faced by the conventional theory of planet formation (based upon the aggregation of microscopic dust particles), we describe an alternative hypothesis. We propose that planets form by gravitational collapse at the same time as the star about which they orbit. This 'concurrent collapse' hypothesis avoids theoretical difficulties associated with the conventional model and suggests satisfying explanations for various poorly understood phenomena. We introduce new explanations for FU Orionis outbursts seen in young stars, the discovery of exoplanets with eccentric orbits and the existence of small rocky objects such as chondrules in the solar system.

preprint1999arXiv

Spectral Statistics in Chaotic Systems with Two Identical Connected Cells

Chaotic systems that decompose into two cells connected only by a narrow channel exhibit characteristic deviations of their quantum spectral statistics from the canonical random-matrix ensembles. The equilibration between the cells introduces an additional classical time scale that is manifest also in the spectral form factor. If the two cells are related by a spatial symmetry, the spectrum shows doublets, reflected in the form factor as a positive peak around the Heisenberg time. We combine a semiclassical analysis with an independent random-matrix approach to the doublet splittings to obtain the form factor on all time (energy) scales. Its only free parameter is the characteristic time of exchange between the cells in units of the Heisenberg time.