Source author record

Debabrata Panja

Debabrata Panja 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

20works
9topics
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

20 published item(s)

preprint2022arXiv

Hidden dependence of spreading vulnerability on topological complexity

Many dynamical phenomena in complex systems concern spreading that plays out on top of networks with changing architecture over time -- commonly known as temporal networks. A complex system's proneness to facilitate spreading phenomena, which we abbreviate as its `spreading vulnerability', is often surmised to be related to the topology of the temporal network featured by the system. Yet, cleanly extracting spreading vulnerability of a complex system directly from the topological information of the temporal network remains a challenge. Here, using data from a diverse set of real-world complex systems, we develop the `entropy of temporal entanglement' as a novel and insightful quantity to measure topological complexities of temporal networks. We show that this parameter-free quantity naturally allows for topological comparisons across vastly different complex systems. Importantly, by simulating three different types of stochastic dynamical processes playing out on top of temporal networks, we demonstrate that the entropy of temporal entanglement serves as a quantitative embodiment of the systems' spreading vulnerability, irrespective of the details of the processes. In being able to do so, i.e., in being able to quantitatively extract a complex system's proneness to facilitate spreading phenomena from topology, this entropic measure opens itself for applications in a wide variety of natural, social, biological and engineered systems.

preprint2022arXiv

Influence maximization under limited network information: Seeding high-degree neighbors

The diffusion of information, norms, and practices across a social network can be initiated by compelling a small number of seed individuals to adopt first. Strategies proposed in previous work either assume full network information or large degree of control over what information is collected. However, privacy settings on the Internet and high non-response in surveys often severely limit available connectivity information. Here we propose a seeding strategy for scenarios with limited network information: Only the degrees and connections of some random nodes are known. This new strategy is a modification of "random neighbor sampling" and seeds the highest-degree neighbors of randomly selected nodes. In simulations of a linear threshold model on a range of synthetic and real-world networks, we find that this new strategy outperforms other seeding strategies, including high-degree seeding and clustered seeding.

preprint2022arXiv

Quantifying agent impacts on contact sequences in social interactions

Human social behavior plays a crucial role in how pathogens like SARS-CoV-2 or fake news spread in a population. Social interactions determine the contact network among individuals, while spreading, requiring individual-to-individual transmission, takes place on top of the network. Studying the topological aspects of a contact network, therefore, not only has the potential of leading to valuable insights into how the behavior of individuals impacts spreading phenomena, but it may also open up possibilities for devising effective behavioral interventions. Because of the temporal nature of interactions - since the topology of the network, containing who is in contact with whom, when, for how long, and in which precise sequence, varies (rapidly) in time - analyzing them requires developing network methods and metrics that respect temporal variability, in contrast to those developed for static (i.e., time-invariant) networks. Here, by means of event mapping, we propose a method to quantify how quickly agents mingle by transforming temporal network data of agent contacts. We define a novel measure called 'contact sequence centrality', which quantifies the impact of an individual on the contact sequences, reflecting the individual's behavioral potential for spreading. Comparing contact sequence centrality across agents allows for ranking the impact of agents and identifying potential 'behavioral super-spreaders'. The method is applied to social interaction data collected at an art fair in Amsterdam. We relate the measure to the existing network metrics, both temporal and static, and find that (mostly at longer time scales) traditional metrics lose their resemblance to contact sequence centrality. Our work highlights the importance of accounting for the sequential nature of contacts when analyzing social interactions.

preprint2022arXiv

Structural dynamics of polycrystalline graphene

The exceptional properties of the two-dimensional material graphene make it attractive for multiple functional applications, whose large-area samples are typically polycrystalline. Here, we study the mechanical properties of graphene in computer simulations and connect these to the experimentally relevant mechanical properties. In particular, we study the fluctuations in the lateral dimensions of the periodic simulation cell. We show that over short time scales, both the area A and the aspect ratio B of the rectangular periodic box show diffusive behavior under zero external field during dynamical evolution, with diffusion coefficients DA and DB that are related to each other. At longer times, fluctuations in A are bounded, while those in B are not. This makes the direct determination of DB much more accurate, from which DA can then be derived indirectly. We then show that the dynamic behavior of polycrystalline graphene under external forces can also be derived from DA and DB via the Nernst-Einstein relation. Additionally, we study how the diffusion coefficients depend on structural properties of the polycrystalline graphene, in particular, the density of defects.

preprint2020arXiv

Super slowing down in the bond-diluted Ising model

In models in statistical physics, the dynamics often slows down tremendously near the critical point. Usually, the correlation time $τ$ at the critical point increases with system size $L$ in power-law fashion: $τ\sim L^z$, which defines the critical dynamical exponent $z$. We show that this also holds for the 2D bond-diluted Ising model in the regime $p>p_c$, where $p$ is the parameter denoting the bond concentration, but with a dynamical critical exponent $z(p)$ which shows a strong $p$-dependence. Moreover, we show numerically that $z(p)$, as obtained from the autocorrelation of the total magnetisation, diverges when the percolation threshold $p_c=1/2$ is approached: $z(p)-z(1) \sim (p-p_c)^{-2}$. We refer to this observed extremely fast increase of the correlation time with size as {\it super slowing down}. Independent measurement data from the mean-square deviation of the total magnetisation, which exhibits anomalous diffusion at the critical point, supports this result.

preprint2016arXiv

Dynamics of a double-stranded DNA segment in a shear flow

We study the dynamics of a double-stranded DNA (dsDNA) segment, as a semiflexible polymer, in a shear flow, the strength of which is customarily expressed in terms of the dimensionless Weissenberg number Wi. Polymer chains in shear flows are well-known to undergo tumbling motion. When the chain lengths are much smaller than the persistence length, one expects a (semiflexible) chain to tumble as a rigid rod. At low Wi, a polymer segment shorter than the persistence length does indeed tumble as a rigid rod. However, for higher Wi the chain does not tumble as a rigid rod, even if the polymer segment is shorter than the persistence length. In particular, from time to time the polymer segment may assume a buckled form, a phenomenon commonly known as Euler buckling. Using a bead-spring Hamiltonian model for extensible dsDNA fragments, we first analyze Euler buckling in terms of the oriented deterministic state (ODS), which is obtained as the steady-state solution of the dynamical equations by turning off the stochastic (thermal) forces at a fixed orientation of the chain. The ODS exhibits symmetry breaking at a critical Weissenberg number Wi$_{\text c}$, analogous to a pitchfork bifurcation in dynamical systems. We then follow up the analysis with simulations and demonstrate symmetry breaking in computer experiments, characterized by a unimodal to bimodal transformation of the probability distribution of the second Rouse mode with increasing Wi. Our simulations reveal that shear can cause strong deformation for a chain that is shorter than its persistence length, similar to recent experimental observations.

preprint2015arXiv

Efficient simulation of semiflexible polymers

Using a recently developed bead-spring model for semiflexible polymers that takes into account their natural extensibility, we report an efficient algorithm to simulate the dynamics for polymers like double-stranded DNA (dsDNA) in the absence of hydrodynamic interactions. The dsDNA is modelled with one bead-spring element per basepair, and the polymer dynamics is described by the Langevin equation. The key to efficiency is that we describe the equations of motion for the polymer in terms of the amplitudes of the polymer's fluctuation modes, as opposed to the use of the physical positions of the beads. We show that, within an accuracy tolerance level of $5\%$ of several key observables, the model allows for single Langevin time steps of $\approx1.6$, 8, 16 and 16 ps for a dsDNA model-chain consisting of 64, 128, 256 and 512 basepairs (i.e., chains of 0.55, 1.11, 2.24 and 4.48 persistence lengths) respectively. Correspondingly, in one hour, a standard desktop computer can simulate 0.23, 0.56, 0.56 and 0.26 ms of these dsDNA chains respectively. We compare our results to those obtained from other methods, in particular, the (inextensible discretised) WLC model. Importantly, we demonstrate that at the same level of discretisation, i.e., when each discretisation element is one basepair long, our algorithm gains about 5-6 orders of magnitude in the size of time steps over the inextensible WLC model. Further, we show that our model can be mapped one-on-one to a discretised version of the extensible WLC model; implying that the speed-up we achieve in our model must hold equally well for the latter. We also demonstrate the use of the method by simulating efficiently the tumbling behaviour of a dsDNA segment in a shear flow.

preprint2014arXiv

Semiflexible Polymer Dynamics with a Bead-Spring Model

We study the dynamical properties of semiflexible polymers with a recently introduced bead-spring model. We focus on double-stranded DNA. The two parameters of the model, $T^*$ and $ν$, are chosen to match its experimental force-extension curve. The bead-spring Hamiltonian is approximated in the first order by the Hessian that is quadratic in the bead positions. The eigenmodels of the Hessian provide the longitudinal (stretching) and transverse (bending) eigenmodes of the polymer, and the corresponding eigenvalues match well with the established phenomenology of semiflexible polymers. Using the longitudinal and transverse eigenmodes, we obtain analytical expressions of (i) the autocorrelation function of the end-to-end vector, (ii) the autocorrelation function of a bond (i.e., a spring, or a tangent) vector at the middle of the chain, and (iii) the mean-square displacement of a tagged bead in the middle of the chain, as sum over the contributions from the modes. We also perform simulations with the full dynamics of the model. The simulations yield numerical values of the correlation functions (i-iii) that agree very well with the analytical expressions for the linearized dynamics. We also study the mean-square displacement of the longitudinal component of the end-to-end vector that showcases strong nonlinear effects in the polymer dynamics, and we identify at least an effective $t^{7/8}$ power-law regime in its time-dependence. Nevertheless, in comparison to the full mean-square displacement of the end-to-end vector the nonlinear effects remain small at all times --- it is in this sense we state that our results demonstrate that the linearized dynamics suffices for dsDNA fragments that are shorter than or comparable to the persistence length. Our results are consistent with those of the wormlike chain (WLC) model, the commonly used descriptive tool of semiflexible polymers.

preprint2013arXiv

Dynamical Eigenmodes of Star and Tadpole Polymers

The dynamics of phantom bead-spring chains with the topology of a symmetric star with $f$ arms and tadpoles ($f=3$, a special case) is studied, in the overdamped limit. In the simplified case where the hydrodynamic radius of the central monomer is $f$ times as heavy as the other beads, we determine their dynamical eigenmodes exactly, along the lines of the Rouse modes for linear bead-spring chains. These eigenmodes allow full analytical calculations of virtually any dynamical quantity. As examples we determine the radius of gyration, the mean square displacement of a tagged monomer, and, for star polymers, the autocorrelation function of the vector that spans from the center of the star to a bead on one of the arms.

preprint2013arXiv

Role of osmotic and hydrostatic pressures in bacteriophage genome ejection

A critical step in the bacteriophage life cycle is genome ejection into host bacteria. The ejection process for double-stranded DNA phages has been studied thoroughly \textit{in vitro}, where after triggering with the cellular receptor the genome ejects into a buffer. The experimental data have been interpreted in terms of the decrease in free energy of the densely packed DNA associated with genome ejection. Here we detail a simple model of genome ejection in terms of the hydrostatic and osmotic pressures inside the phage, a bacterium, and a buffer solution/culture medium. We argue that the hydrodynamic flow associated with the water movement from the buffer solution into the phage capsid and further drainage into the bacterial cytoplasm, driven by the osmotic gradient between the bacterial cytoplasm and culture medium, provides an alternative mechanism for phage genome ejection \textit{in vivo}; the mechanism is perfectly consistent with phage genome ejection \textit{in vitro}.

preprint2013arXiv

Through the Eye of the Needle: Recent Advances in Understanding Biopolymer Translocation

In recent years polymer translocation, i.e., transport of polymeric molecules through nanometer-sized pores and channels embedded in membranes, has witnessed strong advances. It is now possible to observe single-molecule polymer dynamics during the motion through channels with unprecedented spatial and temporal resolution. These striking experimental studies have stimulated many theoretical developments. In this short theory-experiment review, we discuss recent progress in this field with a strong focus on non-equilibrium aspects of polymer dynamics during the translocation process.

preprint2011arXiv

DNA translocation through nanopores with salt gradients: The role of osmotic flow

Recent experiments of translocation of double stranded DNA through nanopores [M. Wanunu \textit{et al.} Nature Nanotech. {\bf 5}, 160 (2010)] reveal that the DNA capture rate can be significantly influenced by a salt gradient across the pore. We show that osmotic flow combined with electrophoretic effects can quantitatively explain the experimental data on the salt-gradient dependence of the capture rate.

preprint2011arXiv

Probabilistic Phase Space Trajectory Description for Anomalous Polymer Dynamics

It has been recently shown that the phase space trajectories for the anomalous dynamics of a tagged monomer of a polymer --- for single polymeric systems such as phantom Rouse, self-avoiding Rouse, Zimm, reptation, and translocation through a narrow pore in a membrane; as well as for many-polymeric system such as polymer melts in the entangled regime --- is robustly described by the Generalized Langevin Equation (GLE). Here I show that the probability distribution of phase space trajectories for all these classical anomalous dynamics for single polymers is that of a fractional Brownian motion (fBm), while the dynamics for polymer melts between the entangled regime and the eventual diffusive regime exhibits small, but systematic deviations from that of a fBm.

preprint2011arXiv

Structural modes of a polymer in the repton model

Using extensive computer simulations, the behavior of the structural modes --- more precisely, the eigenmodes of a phantom Rouse polymer --- are characterized for a polymer in the three-dimensional repton model, and are used to study the polymer's dynamics at time scales well before the tube renewal. Although these modes are not the eigenmodes for a polymer in the repton model, we show that numerically the modes maintain a high degree of statistical independence. The correlations in the mode amplitudes decay exponentially with $(p/N)^2A(t)$, in which $p$ is the mode number, $N$ is the polymer length and $A(t)$ is a single function shared by all modes. In time, the quantity $A(t)$ causes an exponential decay for the mode amplitude correlation functions for times $<1$; a stretched exponential with an exponent 1/2 between times 1 and $τ_R\sim N^2$, the time-scale for diffusion of tagged reptons along the contour of the polymer; and again an exponential decay for times $t>τ_R$. Having assumed statistical independence and the validity of a single function $A(t)$ for all modes, we compute the temporal behavior of three structural quantities: the vectorial distance between the positions of the middle monomer and the center-of-mass, the end-to-end vector, and the vector connecting two nearby reptons around the middle of the polymer. Furthermore, we study the mean-squared displacement of the center-of-mass and the middle repton, and their relation with the temporal behavior of the modes.

preprint2010arXiv

Anomalous Polymer Dynamics Is Non-Markovian: Memory Effects and The Generalized Langevin Equation Formulation

Any first course on polymer physics teaches that the dynamics of a tagged monomer of a polymer is anomalously subdiffusive, i.e., the mean-square displacement of a tagged monomer increases as $t^α$ for some $α<1$ until the terminal relaxation time $τ$ of the polymer. Beyond time $τ$ the motion of the tagged monomer becomes diffusive. Classical examples of anomalous dynamics in polymer physics are single polymeric systems, such as phantom Rouse, self-avoiding Rouse, self-avoiding Zimm, reptation, translocation through a narrow pore in a membrane, and many-polymeric systems such as polymer melts. In this pedagogical paper I report that all these instances of anomalous dynamics in polymeric systems are robustly characterized by power-law memory kernels within a {\it unified} Generalized Langevin Equation (GLE) scheme, and therefore, are non-Markovian. The exponents of the power-law memory kernels are related to the relaxation response of the polymers to local strains, and are derived from the equilibrium statistical physics of polymers. The anomalous dynamics of a tagged monomer of a polymer in these systems is then reproduced from the power-law memory kernels of the GLE via the fluctuation-dissipation theorem (FDT). Using this GLE formulation I further show that the characteristics of the drifts caused by a (weak) applied field on these polymeric systems are also obtained from the corresponding memory kernels.

preprint2010arXiv

Dynamics of Bacteriophage Genome Ejection In Vitro and In Vivo

Bacteriophages, phages for short, are viruses of bacteria. The majority of phages contain a double-stranded DNA genome packaged in a capsid at a density of ~500 mg/ml. This high density requires substantial compression of the normal B form helix, leading to the conjecture that DNA in mature phage virions is under significant pressure, and that pressure is used to eject the DNA during infection. A large number of theoretical, computer simulation and in vitro experimental studies surrounding this conjecture has revealed many --- though often isolated and/or contradictory --- aspects of packaged DNA. This prompts us to present a unified view of the statistical physics and thermodynamics of DNA packaged in phage capsids. We argue that the DNA in a mature phage is in a (meta)stable state, wherein electrostatic self-repulsion is balanced by curvature stress due to confinement in the capsid. We show that in addition to the osmotic pressure associated with the packaged DNA and its counterions, there are four different pressures within the capsid: pressure on the DNA, hydrostatic pressure, the pressure experienced by the capsid, and the pressure associated with the chemical potential of DNA ejection. Significantly, we analyze the mechanism of force transmission in the packaged DNA, and demonstrate that the pressure on DNA is not important for ejection. We derive equations showing a strong hydrostatic pressure difference across the capsid shell. We propose that when a phage is triggered to eject by interaction with its receptor in vitro, the (thermodynamic) incentive of water molecules to enter the phage capsid flushes the DNA out of the capsid. In vivo, the difference between the osmotic pressures in the bacterial cell cytoplasm and the culture medium similarly results in a water flow that drags the DNA out of the capsid and into the bacterial cell.

preprint2010arXiv

Generalized Langevin Equation Formulation for Anomalous Polymer Dynamics

For reproducing the anomalous -- i.e., sub- or super-diffusive -- behavior in some stochastic dynamical systems, the Generalized Langevin Equation (GLE) has gained considerable popularity in recent years. Motivated by the question whether or not a system with anomalous dynamics can have the GLE formulation, here I consider polymer physics, where sub-diffusive behavior is commonplace. I provide an exact derivation of the GLE for phantom Rouse polymers, andby identifying polymeric response to local strains, I argue the case for the GLE formulation for self-avoiding polymers and polymer translocation through a narrow pore in a membrane. The number of instances in polymer physics, where the anomalous dynamics corresponds to the GLE, thus seems to be fairly common.

preprint2010arXiv

Simulations of Two-Dimensional Unbiased Polymer Translocation Using the Bond Fluctuation Model

We use the Bond Fluctuation Model (BFM) to study the pore-blockade times of a translocating polymer of length $N$ in two dimensions, in the absence of external forces on the polymer (i.e., unbiased translocation) and hydrodynamic interactions (i.e., the polymer is a Rouse polymer), through a narrow pore. Earlier studies using the BFM concluded that the pore-blockade time scales with polymer length as $τ_d \sim N^β$, with $β=1+2ν$, whereas some recent studies with different polymer models produce results consistent with $β=2+ν$, originally predicted by us. Here $ν$ is the Flory exponent of the polymer; $ν=0.75$ in 2D. In this paper we show that for the BFM if the simulations are extended to longer polymers, the purported scaling $τ_d \sim N^{1+2ν}$ ceases to hold. We characterize the finite-size effects, and study the mobility of individual monomers in the BFM. In particular, we find that in the BFM, in the vicinity of the pore the individual monomeric mobilities are heavily suppressed in the direction perpendicular to the membrane. After a modification of the BFM which counters this suppression (but possibly introduces other artifacts in the dynamics), the apparent exponent $β$ increases significantly. Our conclusion is that BFM simulations do not rule out our theoretical prediction for unbiased translocation, namely $β=2+ν$.

preprint2005arXiv

Systematic Density Expansion of the Lyapunov Exponents for a Two-dimensional Random Lorentz Gas

We study the Lyapunov exponents of a two-dimensional, random Lorentz gas at low density. The positive Lyapunov exponent may be obtained either by a direct analysis of the dynamics, or by the use of kinetic theory methods. To leading orders in the density of scatterers it is of the form $A_{0}\tilde{n}\ln\tilde{n}+B_{0}\tilde{n}$, where $A_{0}$ and $B_{0}$ are known constants and $\tilde{n}$ is the number density of scatterers expressed in dimensionless units. In this paper, we find that through order $(\tilde{n}^{2})$, the positive Lyapunov exponent is of the form $A_{0}\tilde{n}\ln\tilde{n}+B_{0}\tilde{n}+A_{1}\tilde{n}^{2}\ln\tilde{n} +B_{1}\tilde{n}^{2}$. Explicit numerical values of the new constants $A_{1}$ and $B_{1}$ are obtained by means of a systematic analysis. This takes into account, up to $O(\tilde{n}^{2})$, the effects of {\it all\/} possible trajectories in two versions of the model; in one version overlapping scatterer configurations are allowed and in the other they are not.