Researcher profile

Hendrick W. de Haan

Hendrick W. de Haan contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
9works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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 map preview

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

Published work

9 published item(s)

preprint2022arXiv

Studying First Passage Problems using Neural Networks: A Case Study in the Slit-Well Microfluidic Device

This study presents deep neural network solutions to a time-integrated Smoluchowski equation modeling the mean first passage time of nanoparticles traversing the slit-well microfluidic device. This physical scenario is representative of a broader class of parameterized first passage problems in which key output metrics are dictated by a complicated interplay of problem parameters and system geometry. Specifically, whereas these types of problems are commonly studied using particle simulations of stochastic differential equation models, here the corresponding partial differential equation model is solved using a method based on deep neural networks. The results illustrate that the neural network method is synergistic with the time-integrated Smoluchowski model: together, these are used to construct continuous mappings from key physical inputs (applied voltage and particle diameter) to key output metrics (mean first passage time and effective mobility). In particular, this capability is a unique advantage of the time-integrated Smoluchowski model over the corresponding stochastic differential equation models. Furthermore, the neural network method is demonstrated to easily and reliably handle geometry-modifying parameters, which is generally difficult to accomplish using other methods.

preprint2020arXiv

Neural Network Solutions to Differential Equations in Non-Convex Domains: Solving the Electric Field in the Slit-Well Microfluidic Device

The neural network method of solving differential equations is used to approximate the electric potential and corresponding electric field in the slit-well microfluidic device. The device's geometry is non-convex, making this a challenging problem to solve using the neural network method. To validate the method, the neural network solutions are compared to a reference solution obtained using the finite element method. Additional metrics are presented that measure how well the neural networks recover important physical invariants that are not explicitly enforced during training: spatial symmetries and conservation of electric flux. Finally, as an application-specific test of validity, neural network electric fields are incorporated into particle simulations. Conveniently, the same loss functional used to train the neural networks also seems to provide a reliable estimator of the networks' true errors, as measured by any of the metrics considered here. In all metrics, deep neural networks significantly outperform shallow neural networks, even when normalized by computational cost. Altogether, the results suggest that the neural network method can reliably produce solutions of acceptable accuracy for use in subsequent physical computations, such as particle simulations.

preprint2016arXiv

Waves of DNA: Propagating Excitations in Extended Nanoconfined Polymers

We use a nanofluidic system to investigate the emergence of thermally driven collective phenomena along a single polymer chain. In our approach, a single DNA molecule is confined in a nanofluidic slit etched with arrays of embedded nanocavities; the cavity lattice is designed so that a single chain occupies multiple cavities. Fluorescent video-microscopy data shows that waves of excess fluorescence propagate across the cavity-straddling molecule, corresponding to propagating fluctuations of contour overdensity in the cavities. The waves are quantified by examining the correlation in intensity fluctuations between neighbouring cavities. Correlations grow from an anti-correlated minimum to a correlated maximum before decaying, corresponding to a transfer of contour between neighbouring cavities at a fixed transfer time-scale. The observed dynamics can be modelled using Langevin dynamics simulations and a minimal lattice model of coupled diffusion. This study shows how confinement-based sculpting of the polymer equilibrium configuration, by renormalizing the physical system into a series of discrete cavity states, can lead to new types of dynamic collective phenomena.

preprint2014arXiv

Coarse-Grained Molecular Dynamics Simulations of Depletion-Induced Interactions for Soft Matter Systems

Given the ubiquity of depletion effects in biological and other soft matter systems, it is desirable to have coarse-grained Molecular Dynamics simulation approaches appropriate for the study of complex systems. This paper examines the use of two common truncated Lennard-Jones (WCA) potentials to describe a pair of colloidal particles in a thermal bath of depletants. The shifted-WCA model is the steeper of the two repulsive potentials considered, while the combinatorial-WCA model is the softer. It is found that the depletion-induced well depth for the combinatorial-WCA model is significantly deeper than the shifted-WCA model because the resulting overlap of the colloids yields extra accessible volume for depletants. For both shifted- and combinatorial-WCA simulations, the second virial coefficients and pair potentials between colloids are demonstrated to be well approximated by the Morphometric Thermodynamics (MT) model. This agreement suggests that the presence of depletants can be accurately modelled in MD simulations by implicitly including them through simple, analytical MT forms for depletion-induced interactions. Although both WCA potentials are found to be effective generic coarse-grained simulation approaches for studying depletion effects in complicated soft matter systems, combinatorial-WCA is the more efficient approach as depletion effects are enhanced at lower depletant densities. The findings indicate that for soft matter systems that are better modelled by potentials with some compressibility, predictions from hard-sphere systems could greatly underestimate the magnitude of depletion effects at a given depletant density.

preprint2014arXiv

Force-Extension for DNA in a Nanoslit: Using an Effective Dimensionality to Map between the 3D and 2D Limits

The force-extension relation for a semi-flexible polymer such as DNA confined in a nanoslit is investigated and it is found that both the effective persistence length and the form of the force-extension relation change as the chain goes from 3D (very large slit heights) to 2D (very tight confinement). Generalizations of the Marko-Siggia relation appropriate for polymers in nanoconfinement are presented. The forms for both strong and weak confinement regimes are characterized by an \textit{effective dimensionality}. At low forces, the effective dimensionality is given by the correlations along the polymer in the plane of the confining walls. At high forces, the theoretical force must account for reduced conformation space. Together the interpolations give good agreement for all slit heights at all forces. As DNA and other semi-flexible biopolymers are commonly confined \textit{in situ} to various degrees, both the idea of an effective dimensionality and the associated generalized Marko-Siggia interpolations are useful for qualitatively understanding and quantitatively modeling polymers in nanoconfinement.

preprint2014arXiv

Simulating the Entropic Collapse of Coarse-Grained Chromosomes

Depletion forces play a role in the compaction and de-compation of chromosomal material in simple cells but it remains debatable whether they are sufficient to account for chromosomal collapse. We present coarse-grained molecular dynamics simulations, which reveal that depletion-induced attraction is sufficient to cause the collapse of a flexible chain of large structural monomers immersed in a bath of smaller depletants. These simulations use an explicit coarse-grained computational model that treats both the supercoiled DNA structural monomers and the smaller protein crowding agents as combinatorial, truncated Lennard-Jones spheres. By presenting a simple theoretical model, we quantitatively cast the action of depletants on supercoiled bacterial DNA as an effective solvent quality. The rapid collapse of the simulated flexible chromosome at the predicted volume fraction of depletants is a continous phase transition. Additional physical effects to such simple chromosome models, such as enthalpic interactions between structural monomers or chain rigidity, are required if the collapse is to be a first-order phase transition.

preprint2014arXiv

Using a Péclet Number for the Translocation of Polymer through a Nanopore to Tune Coarse-Grained Simulations to Experimental Conditions

Coarse-grained simulations are often employed to study the translocation of DNA through a nanopore. The majority of these studies investigate the translocation process in a relatively generic sense and do not endeavour to match any particular set of experimental conditions. In this manuscript, we use the concept of a Péclet number for translocation, $P_t$, to compare the drift-diffusion balance in a typical experiment vs a typical simulation. We find that the standard coarse-grained approach over-estimates diffusion effects by anywhere from a factor of 5 to 50 compared to experimental conditions using dsDNA. By defining a coarse-graining parameter, $λ$, we are able to correct this and tune the simulations to replicate the experimental $P_t$ (for dsDNA and other scenarios). To show the effect that a particular $P_t$ can have on the dynamics of translocation, we perform simulations across a wide range of $P_t$ values for two different types of driving forces: a force applied in the pore and a pulling force applied to the end of the polymer. As $P_t$ brings the system from a diffusion dominated to a drift dominated regime, a variety of effects are observed including a non-monotonic dependence of the translocation time $τ$ on $P_t$ and a steep rise in the probability of translocating. Comparing the two force cases illustrates the impact of the crowding effects that occur on the trans side: a non-monotonic dependence of the width of the $τ$ distributions is obtained for the in-pore force but not for the pulling force.

preprint2012arXiv

Monte Carlo Approaches for Simulating a Particle at a Diffusivity Interface and the "Ito--Stratonovich Dilemma"

The possibility of different interpretations of the stochastic term (or calculi) in the overdamped Langevin equation for the motion of a particle in an inhomogeneous medium is often referred to as the "Ito--Stratonovich dilemma," although there is, in fact, a continuum of choices. We introduce two Monte Carlo (MC) simulation approaches for studying such systems, with both approaches giving the choice between different calculi (in particular, Ito, Stratonovich, and "isothermal"). To demonstrate these approaches, we study the diffusion on a 1D interval of a particle released at an interface (in the middle of the system) between two media where this particle has different diffusivities (e.g., two fluids with different viscosities). For reflecting boundary conditions at the ends of the interval, a discontinuity at the interface in the stationary-state particle distribution is found, except for the isothermal case, as expected. We also study the first-passage problem using absorbing boundary conditions. Good agreement is found when comparing the MC approaches to theory as well as Brownian and Langevin dynamics simulations. Also, the results themselves turn out to be interesting. For instance: 1) for some calculi, there are more particles on the low-viscosity side (LVS) at earlier times and then fewer particles at later times; 2) there is no preference to end up on a particular wall for the Ito variant, but a bias towards the wall on the LVS in all other cases; 3) the mean first-passage time to the wall on the LVS grows as the viscosity on the other side is increased, except for the isothermal case where it approaches a constant; 4) when the viscosity ratio is high, the first-passage-time distribution for the wall on the LVS is much broader than for the other wall, with a power-law dependence of the former in a certain time interval whose exponent depends on the calculus.

preprint2012arXiv

Translocation of a Polymer through a Nanopore across a Viscosity Gradient

The translocation of a polymer through a pore in a membrane separating fluids of different viscosities is studied via several computational approaches. Starting with the polymer halfway, we find that as a viscosity difference across the pore is introduced, translocation will predominately occur towards one side of the membrane. These results suggest an intrinsic pumping mechanism for translocation across cell walls which could arise whenever the fluid across the membrane is inhomogeneous. Somewhat surprisingly, the sign of the preferred direction of translocation is found to be strongly dependent on the simulation algorithm: for Langevin Dynamics (LD) simulations, a bias towards the low viscosity side is found while for Brownian Dynamics (BD), a bias towards the high viscosity is found. Examining the translocation dynamics in detail across a wide range of viscosity gradients and developing a simple force model to estimate the magnitude of the bias, the LD results are demonstrated to be more physically realistic. The LD results are also compared to those generated from a simple, one dimensional random walk model of translocation to investigate the role of the internal degrees of freedom of the polymer and the entropic barrier. To conclude, the scaling of the results across different polymer lengths demonstrates the saturation of the preferential direction with polymer length and the non-trivial location of the maximum in the exponent corresponding to the scaling of the translocation time with polymer length.