Source author record

Takahiro Sakaue

Takahiro Sakaue 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

27works
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

27 published item(s)

preprint2022arXiv

Mechanical properties of Nucleic Acids and the non-local Twistable Wormlike Chain model

Mechanical properties of nucleic acids play an important role in many biological processes which often involve physical deformations of these molecules. At sufficiently long length scales (say above $\sim 20-30$ base pairs) the mechanics of DNA and RNA double helices is described by a homogeneous Twistable Wormlike Chain (TWLC), a semiflexible polymer model characterized by twist and bending stiffnesses. At shorter scales this model breaks down for two reasons: the elastic properties become sequence-dependent and the mechanical deformations at distal sites gets coupled. We discuss in this paper the origin of the latter effect using the framework of a non-local Twistable Wormlike Chain (nlTWLC). We show, by comparing all-atom simulations data for DNA and RNA double helices, that the non-local couplings are of very similar nature in these two molecules: couplings between distal sites are strong for tilt and twist degrees of freedom and weak for roll. We introduce and analyze a simple double-stranded polymer model which clarifies the origin of this universal distal couplings behavior. In this model, referred to as the ladder model, a nlTWLC description emerges from the coarsening of local (atomic) degrees of freedom into angular variables which describe the twist and bending of the molecule. Differently from its local counterpart, the nlTWLC is characterized by a length-scale-dependent elasticity. Our analysis predicts that nucleic acids are mechanically softer at the scale of a few base pairs and are asymptotically stiffer at longer length scales, a behavior which matches experimental data.

preprint2021arXiv

Quantifying the mobility of chromatin during embryogenesis: Nuclear size matters

Chromatin moves dynamically inside the cell nucleus, and its motion is often correlated with gene functions such as DNA recombination and transcription. A recent study has shown that during early embryogenesis of the nematode, Caenorhabiditis elegans, the chromatin motion markedly decreases. However, the underlying mechanism for this transition has yet to be elucidated. We systematically investigated the impact of nuclear size to demonstrate that it is indeed a decisive factor in chromatin mobility. To this end, we established a method to quantify chromatin motion inside the nucleus, while excluding the contribution of the movement of the nucleus itself, which allowed us to extract the intrinsic mean-squared displacement (iMSD) of individual chromosomal loci in moving nuclei from the correlated motion of two loci. We show that a simple theoretical description, which takes into account the topological constraints of chromatin polymers, can quantitatively describe the relationship between the nucleus size and the chromatin motion in vivo. Our results emphasize a regulatory role of nuclear size in restricting chromatin motion, and a generic polymer physics model plays a guiding role in capturing this essential feature.

preprint2020arXiv

Dynamical Entanglement and Cooperative Dynamics in Entangled Solutions of Ring and Linear Polymers

Understanding how entanglements affect the behaviour of polymeric complex fluids is an open challenge in many fields. To elucidate the nature and consequence of entanglements in dense polymer solutions, we propose a novel method: a "dynamical entanglement analysis" (DEA) to extract spatio-temporal entanglement structures from the pair-wise displacement correlation of entangled chains. By applying this method to large-scale Molecular Dynamics simulations of linear and unknotted, nonconcatenated ring polymers, we find a strong and unexpected cooperative dynamics: the footprint of mutual entrainment between entangled chains. We show that DEA is a powerful and sensitive probe that reveals previously unnoticed, and architecture-dependent, spatio-temporal structures of dynamical entanglement in polymeric solutions. We also propose a mean-field approximation of our analysis which provides previously under-appreciated physical insights into the dynamics of generic entangled polymers. We envisage DEA will be useful to analyse the dynamical evolution of entanglements in generic polymeric systems such as blends and composites.

preprint2020arXiv

Stabilization of straight longitudinal dune under bimodal wind with large directional variation

It has been observed that the direction in which a sand dune extends its crest line depends on seasonal variation of wind direction; when the variation is small, the crest line develops more or less perpendicularly to the mean wind direction to form a transverse dune with some undulation. In the case of bimodal wind with a large relative angle, however, the dune extends its crest along the mean wind direction and evolves into an almost straight longitudinal dune. Motivated by these observations, we investigate the dynamical stability of isolated dunes using the crest line model, where the dune dynamics is represented by its crest line motion. First, we extend the previous linear stability analysis under the unidirectional wind to the case with non-zero slant angle between the wind direction and the normal direction of the crest line, and show that the stability diagram does not depend on the slant angle. Secondly, we examine how the linear stability is affected by the seasonal changes of wind direction in the case of bimodal wind with equal strength and duration. For the transverse dune, we find that the stability is virtually the same with that for the unidirectional wind as long as the dune evolution during a season is small. On the other hand, in the case of the longitudinal dune, the dispersions of the growth rates for the perturbation are drastically different from those of the unidirectional wind, and we find that the largest growth rate is always located at $k = 0$. This is because the growth of the perturbation with $k \ne 0$ is canceled by the alternating wind from opposite sides of the crest line even though it grows during each duration period of the bimodal wind.

preprint2016arXiv

Miscibility Phase Diagram of Ring Polymer Blends: A Topological Effect

The miscibility of polymer blends, a classical problem in polymer science, may be altered, if one or both of the component do not have chain ends. Based on the idea of {\it topological volume}, we propose a mean-field theory to clarify how the topological constraints in ring polymers affect the phase behavior of the blends. While the large enhancement of the miscibility is expected for ring-linear polymer blends, the opposite trend toward demixing, albeit comparatively weak, is predicted for ring-ring polymer blends. Scaling formulas for the shift of critical point for both cases are derived. We discuss the valid range of the present theory, and the crossover to the linear polymer blends behaviors, which is expected for short chains. These analysis put forward a view that the topological constraints could be represented as an effective excluded-volume effects, in which the topological length plays a role of the screening factor.

preprint2016arXiv

Torque-Induced Rotational Dynamics in Polymers: Torsional Blobs and Thinning

By using the blob theory and computer simulations, we investigate the properties of a linear polymer performing a stationary rotational motion around a long impenetrable rod. In particular, in the simulations the rotation is induced by a torque applied to the end of the polymer that is tethered to the rod. Three different regimes are found, in close analogy with the case of polymers pulled by a constant force at one end. For low torques the polymer rotates maintaining its equilibrium conformation. At intermediate torques the polymer assumes a trumpet shape, being composed by blobs of increasing size. At even larger torques the polymer is partially wrapped around the rod. We derive several scaling relations between various quantities as angular velocity, elongation and torque. The analytical predictions match the simulation data well. Interestingly, we find a "thinning" regime where the torque has a very weak (logarithmic) dependence on the angular velocity. We discuss the origin of this behavior, which has no counterpart in polymers pulled by an applied force.

preprint2015arXiv

Driven Anomalous Diffusion: An example from Polymer Stretching

The way tension propagates along a chain is a key to govern many of anomalous dynamics in macromolecular systems. After introducing the weak and the strong force regimes of the tension propagation, we focus on the latter, in which the dynamical fluctuations of a segment in a long polymer during its stretching process is investigated. We show that the response, i.e., average drift, is anomalous, which is characterized by the nonlinear memory kernel, and its relation to the fluctuation is nontrivial. These features are discussed on the basis of the generalized Langevin equation, in which the role of the temporal change in spring constant due to the stress hardening is pinpointed. We carried out the molecular dynamics simulation, which supports our theory.

preprint2014arXiv

Dynamical Scaling of Polymerized Membranes

Monte Carlo simulations have been performed to analyze the sub-diffusion dynamics of a tagged monomer in self-avoiding polymerized membranes in the flat phase. By decomposing the mean square displacement into the out-of-plane ($\parallel$) and the in-plane ($\perp$) components, we obtain good data collapse with two distinctive diffusion exponents $2 α_{\parallel} = 0.36 \pm 0.01$ and $2 α_{\perp} = 0.21 \pm 0.01$, and the roughness exponents $ζ_{\parallel} = 0.6 \pm 0.05$ and $ζ_{\perp} = 0.25 \pm 0.05 $, respectively for each component. Their values are consistent with the relation from the rotational symmetry. We derive the generalized Langevin equations to describe the sub-diffusional behaviors of a tagged monomer in the intermediate time regime where the collective effect of internal modes in the membrane dominate the dynamics to produce negative memory kernels with a power-law. We also briefly discuss how the long-range hydrodynamic interactions alter the exponents.

preprint2014arXiv

Dynamics and electro-rheology of sheared immiscible fluid mixtures

We analyze the electro-rheological effect in immiscible fluid mixtures with dielectric mismatch. By taking the electric field effect into account, which couples to the dynamics of domain morphology under flow, we propose a set of electro-rheological constitutive equations valid under the condition where the relative magnitude of the flow field is stronger than that of the electric field. Through the comparison with recent experiment, we point out a unique dynamical stress response inherent in situations, where the cross-coupling between different fields is essential.

preprint2014arXiv

Nano-pore based characterization of branched polymers

We propose a novel characterization method of randomly branched polymers based on the geometrical property of such objects in confined spaces. The central idea is that randomly branched polymers exhibit passing/clogging transition across the nano-channel as a function of the channel size. This critical channel size depends on the degree of the branching, whereby allowing the extraction of the branching information of the molecule.

preprint2013arXiv

Cis-Trans Dynamical Asymmetry in Driven Polymer Translocation

During polymer translocation driven by e.g. voltage drop across a nanopore, the segments in the cis-side is incessantly pulled into the pore, which are then pushed out of it into the trans-side. This pulling and pushing polymer segments are described in the continuum level by nonlinear transport processes known, respectively, as fast and slow diffusions. By matching solutions of both sides through the mass conservation across the pore, we provide a physical basis for the cis and trans dynamical asymmetry, a feature repeatedly reported in recent numerical simulations. We then predict how the total driving force is dynamically allocated between cis (pulling) and trans (pushing) sides, demonstrating that the trans-side event adds a finite-chain length effect to the dynamical scaling, which may become substantial for weak force and/or high pore friction cases.

preprint2013arXiv

Hamilton-Jacobi method for molecular distribution function in a chemical oscillator

Using the Hamilton-Jacobi method, we solve chemical Fokker-Planck equations within the Gaussian approximation and obtain a simple and compact formula for a conditional probability distribution. The formula holds in general transient situations, and can be applied not only for a steady state but also for a oscillatory state. By analyzing the long time behavior of the solution in the oscillatory case, we obtain the phase diffusion constant along the periodic orbit and the steady distribution perpendicular to it. A simple method for numerical evaluation of these formulas are devised, and they are compared with Monte Carlo simulations in the case of Brusselator as an example. Some results are shown to be identical to previously obtained expressions.

preprint2013arXiv

Inelastic collapse in one-dimensional driven systems under gravity

We study the inelastic collapse in the one-dimensional $N$-particle systems in the situation where the system is driven from below under the gravity. We investigate the hard-sphere limit of the inelastic soft-sphere systems by numerical simulations to find how the collision rate per particle $n_{coll}$ increases as a function of the elastic constant of the sphere $k$ when the restitution coefficient $e$ is kept constant. For the systems with large enough $N \agt 20$, we find three regimes in $e$ depending on the behavior of $n_{coll}$ in the hard-sphere limit: (i) uncollapsing regime for $1 \ge e > e_{c1}$, where $n_{coll}$ converges to a finite value, (ii) logarithmically collapsing regime for $e_{c1} > e > e_{c2}$, where $n_{coll}$ diverges as $n_{coll} \sim \log k$, and (iii) power-law collapsing regime for $e_{c2} > e > 0$, where $n_{coll}$ diverges as $n_{coll} \sim k^α$ with an exponent $α$ that depends on $N$. The power-law collapsing regime shrinks as $N$ decreases and seems not to exist for the system with N=3 while, for large $N$, the size of the uncollapsing and the logarithmically collapsing regime decreases as $e_{c1} \simeq 1-2.6/N$ and $e_{c2} \simeq 1-3.0/N$. We demonstrate that this difference between large and small systems exists already in the inelastic collapse without the external drive and the gravity.

preprint2012arXiv

Dragging a polymer in a viscous fluid: steady-state and transient

We study the conformation and dynamics of a single polymer chain that is pulled by a constant force applied at its one end with the other end free. Such a situation is relevant to the growing technology of manipulating individual macromolecules, which offers a paradigm research for probing far-from-equilibrium responses of long flexible biological polymers. We first analyze the Rouse model for the Gaussian chains for which the exact analytical results can be obtained. More realistic features such as the finite extensibility, the excluded volume and the hydrodynamic interactions are taken into account with the help of the scaling argument, which leads to various nontrivial predictions such as the stretching-force-dependent friction constants. We elucidate (i) generalized dynamical equations of state describing extension/friction laws in steady-state and (ii) the tension propagation laws in the transient process. We point out that the time evolutions of the dynamic friction in the transient process crucially depend on the experimental protocol, i.e., either constant force or velocity ensemble, which might be detectable in experiments using giant DNAs and chromosomes.

preprint2012arXiv

Easier sieving through narrower pores: fluctuations and barrier crossing in flow-driven polymer translocation

We show that the injection of polymer chains into nanochannels becomes easier as the channel becomes narrower. This counter intuitive result arises because of a decrease in the diffusive time scale of the chains with increasing confinement. The results are obtained by extending the de Gennes blob model of confined polymers, and confirmed by hybrid molecular dynamics - lattice-Boltzmann simulations.

preprint2012arXiv

Localization and size distribution of a polymer knot confined in a channel

We have examined the behaviors of a knotted linear polymer in narrow tubes using Langevin dynamics simulation to investigate the knot localization property in one-dimensional (1D) geometry. We have found that the knot is strongly localized in such a geometry. By observing the distribution function of the size of localized knot, we found the scaling behavior of the fluctuation around the most probable size with radius of confinement. Based on the analysis of the probability distribution of the knot size, we show that the strong localization behavior and the fluctuation around the most probable size can be encompassed by a simple argument based on the virtual tubes composed of parallel strands and overlapping among them.

preprint2012arXiv

Maxwell stress in fluid mixtures

We examine the structure of Maxwell stress in binary fluid mixtures under an external electric field and discuss its consequence. In particular, we show that, in immiscible blends, it is intimately related to the statistics of domain structure. This leads to a compact formula, which may be useful in the investigation of electro-rheological effects in such systems. The stress tensor calculated in a phase separated fluid under a steady electric field is in a good agreement with recent experiments.

preprint2012arXiv

Process Time Distribution of Driven Polymer Transport

We discuss the temporal distribution of dynamic processes in driven polymer transport inherent to flexible chains due to stochastic tension propagation. The stochasticity originates from the disordered initial configuration of an equilibrium polymer coil, which results in random paths for tension propagation. We consider the process time for when translocation occurs across a fixed pore and when stretching occurs by pulling the chain end. A scaling argument for the mean and standard deviation of the process time is provided using the two-phase picture for stochastic propagation. The two cases are found to differ remarkably. The process time distribution of the translocation exhibits substantial spreading even in the long-chain limit, unlike that found for the dynamics of polymer stretching. In addition, the process time distribution in the driven translocation is shown to have a characteristic asymmetric shape.

preprint2012arXiv

Statistics and Geometrical Picture of Ring Polymer Melts and Solutions

We present a detailed account of a recently proposed phenomenological theory for noncatenated ring polymer melts (Phys. Rev. Lett. 106, 167802 (2011)). A basic assumption lies in the implementation of the noncatenation constraint via the effective excluded-volume effect, from which a geometrical picture of melts emerges. The result captures many of the salient features observed so far, including (i) the overall spacial size of rings, (ii) the coordinate number, i.e., the number of rings surrounding a given ring, (iii) the topological length scale as a function of the molecular weight and (iv) the effect of the chain stiffness and concentration. We also suggest a geometrical interpretation of the topological length scale, which may shed some light on the entanglement concept in polymeric systems.

preprint2012arXiv

Transcription fluctuation effects on biochemical oscillations

Biochemical oscillation systems often consist of negative feedback loops with repressive transcription regulation. Such systems have distinctive characteristics in comparison with ordinary chemical systems: i) the numbers of molecules involved are small, ii) there are typically only a couple of genes in a cell with a finite regulation time scale. Due to the fluctuations caused by these features, the system behavior can be quite different from the one obtained by rate equations, because the rate equations ignore molecular fluctuations and thus are exact only in the infinite molecular number limit. The molecular fluctuations on a free-running circadian system have been studied by Gonze et al. (2002) by introducing a scale parameter $Ω$ for the system size. They consider, however, only the first effect, assuming that the gene process is fast enough for the second effect to be ignored, but this has not been examined systematically yet. In this work, we study fluctuation effects due to the finite gene regulation time by introducing a new scale parameter $τ$, which we take as the unbinding time of a nuclear protein from the gene. We focus on the case where the fluctuations due to small molecular numbers can be ignored. In simulations on the same system studied by Gonze et al., we find the system is sensitive to the fluctuation in the transcription regulation; the period of oscillation fluctuates about 30 min even when the regulation time scale $τ$ is around 30 s, that is even smaller than 1/1000 of its circadian period. We also demonstrate that the distribution width for the oscillation period and the amplitude scales with $\sqrtτ$, and the correlation time of the oscillation scales with $1/τ$ in the small $τ$ regime. The relative fluctuations for the period are about half of that for the amplitude, namely, the periodicity is more stable than the amplitude.

preprint2012arXiv

Two phase picture in driven polymer translocation

Two phase picture is a simple and effective methodology to capture the nonequilibrium dynamics of polymer associated with tension propagation. When applying it to the driven translocation process, there is a point to be noted, as briefly discussed in our recent article [Phys. Rev. E 85, 061803 (2012)]. In this article, we address this issue in detail and modify our previous prediction [Euro. Phys. J. E 34, 135 (2011)] by adopting an alternative steady-state ansatz. The modified scaling prediction turns out to be the same as that of the iso-flux model recently proposed by Rowghanian and Grosberg [J. Phys. Chem. B 115, 14127-14135 (2011)].

preprint2011arXiv

Computing a Knot Invariant as a Constraint Satisfaction Problem

We point out the connection between mathematical knot theory and spin glass/search problem. In particular, we present a statistical mechanical formulation of the problem of computing a knot invariant; p-colorability problem, which provides an algorithm to find the solution. The method also allows one to get some deeper insight into the structural complexity of knots, which is expected to be related with the landscape structure of constraint satisfaction problem.

preprint2011arXiv

Dynamical Diagram and Scaling in Polymer Driven Translocation

By analyzing the real space nonequilibrium dynamics of polymers, we elucidate the physics of driven translocation and propose its dynamical scaling scenario analogous to that in the surface growth phenomena. We provide a detailed account of the previously proposed tension-propagation formulation and extend it to cover the broader parameter space relevant to real experiments; In addition to a near-equilibrium regime, we identify three distinct nonequilibrium regimes reflecting the steady-state property of a dragged polymer with the finite extensibility. Finite size effects are also pointed out. These elements are shown to be crucial for the appropriate comparison with experiments and simulations.

preprint2011arXiv

Length-dependent translocation of polymers through nanochannels

We consider the flow-driven translocation of single polymer chains through nanochannels. Using analytical calculations based on the de Gennes blob model and mesoscopic numerical simulations, we estimate the threshold flux for the translocation of chains of different number of monomers. The translocation of the chains is controlled by the competition between entropic and hydrodynamic effects, which set a critical penetration length for the chain before it can translocate through the channel. We demonstrate that the polymers show two different translocation regimes depending on how their length under confinement compares to the critical penetration length. For polymer chains longer than the threshold, the translocation process is insensitive to the number of monomers in the chain as predicted in Sakaue {\it et al.}, {\it Euro. Phys. Lett.}, {\bf 72} 83 (2005). However, for chains shorter than the critical length we show that the translocation process is strongly dependent on the length of the chain. We discuss the possible relevance of our results to biological transport.

preprint2011arXiv

Ring polymers in melts and solutions: scaling and crossover

We propose a simple mean-field theory for the structure of ring polymer melts. By combining the notion of topological volume fraction and a classical van der Waals theory of fluids, we take into account many body effects of topological origin in dense systems. We predict that although the compact statistics with the Flory exponent $ν=1/3$ is realized for very long chains, most practical cases fall into the crossover regime with the apparent exponent $ν= 2/5$ during which the system evolves toward a topological dense-packed limit.

preprint2010arXiv

How does torsional rigidity affect the wrapping transition of a semiflexible chain around a spherical core?

We investigated the effect of torsional rigidity of a semiflexible chain on the wrapping transition around a spherical core, as a model of nucleosome, the fundamental unit of chromatin. Through molecular dynamics simulation, we show that the torsional effect has a crucial effect on the chain wrapping around the core under the topological constraints. In particular, the torsional stress (i) induces the wrapping/unwrapping transition, and (ii) leads to a unique complex structure with an antagonistic wrapping direction which never appears without the topological constraints. We further examine the effect of the stretching stress for the nucleosome model, in relation to the unique characteristic effect of the torsional stress on the manner of wrapping.

preprint2005arXiv

Flow-injection of branched polymers inside nanopores

Flexible chains (linear or branched) can be forced to enter into a narrow capillary by using a hydrodynamic flow. Here, we correct our earlier description of this problem by considering the progressive nature of the suction process. We find that the critical current for penetration, $J\_c$, is controlled by the entry of a single blob of the capillary size, and that its scaling structure is the same for branched and linear chains.