Researcher profile

Hiroshi Koibuchi

Hiroshi Koibuchi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
19works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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 graph slice

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

Published work

19 published item(s)

preprint2022arXiv

Langevin Navier-Stokes simulation of protoplasmic streaming by 2D MAC method

We study protoplasmic streaming in plant cells such as chara brauni by simplifying the flow field to a two-dimensional Couette flow with Brownian random motion inside parallel plates. Protoplasmic streaming is receiving a lot of attention in many areas, such as agriculture-technology and biotechnology. The plant size depends on the velocity of streaming and the driving force originating in molecular motors. Therefore, it is interesting to study detailed information on the velocity of streaming. Recently, experimentally observed peaks in the velocity distribution have been simulated by a 2D Langevin Navier-Stokes (LNS) equation for vortex and flow function. However, to simulate actual 3D flows, we have to use the NS equation for velocity, which, in the case of 2D flows, is not always equivalent to that for vorticity and stream function. In this paper, we report that a 2D LNS equation for velocity and pressure successfully simulates protoplasmic streaming by comparing the results with the experimental data and those obtained by 2D LNS simulations for vortex and flow function. Moreover, a dimensional analysis clarifies the dependence of numerical results on the strength $D$ of Brownian random force and physical parameters such as kinematic viscosity and cell size. We find from this analysis how the peak position in normalized velocity distribution moves depending on these parameters.

preprint2022arXiv

The stability of 3D skyrmions under mechanical stress studied via Monte Carlo calculations

Using Monte Carlo (MC) simulations, we study the skyrmion stability/instability as a response to uniaxial mechanical stresses. Skyrmions emerge in chiral magnetic materials as a stable spin configuration under external magnetic field $\vec{B}$ with the competition of ferromagnetic interaction and Dzyaloshinskii-Moriya interaction (DMI) at low temperature $T$. Skyrmion configurations are also known to be stable (unstable) under a compressive stress applied parallel (perpendicular) to $\vec{B}$. To understand the origin of such experimentally confirmed stability/instability, we use the Finsler geometry modeling technique with a new degree of freedom for strains, which plays an essential role in DMI being anisotropic. We find from MC data that the area of the skyrmion state on the $B$-$T$ phase diagram increases (decreases) depending on the direction of applied stresses, in agreement with reported experimental results. This change in the area of the skyrmion state indicates that skyrmions become more (less) stable if the tensile strain direction is parallel (perpendicular) to $\vec{B}$. From the numerical data in this paper, we find that the so-called magneto-elastic effect is suitably implemented in the effective DMI theory with the strain degree of freedom without complex magneto-elastic coupling terms for chiral magnetic materials. This result confirms that experimentally-observed skyrmion stability and instability are caused by DMI anisotropy.

preprint2020arXiv

Stochastic Fluid Dynamics Simulations of the Velocity Distribution in Protoplasmic Streaming

Protoplasmic streaming in plant cells is directly visible in the cases of \textit{Chara corallina} and \textit{Nitella flexilis}, and this streaming is understood to play a role in the transport of biological materials. For this reason, related studies have focused on molecular transportation from a fluid mechanics viewpoint. However, the experimentally observed distribution of the velocity along the flow direction $x$, which exhibits two peaks at $V_x\!=\!0$ and at a finite $V_x(\not=\!0)$, remains to be studied. In this paper, we numerically study whether this behavior of the flow field can be simulated by a 2D stochastic Navier-Stokes (NS) equation for Couette flow, in which random Brownian force is assumed. We present the first numerical evidence that these peaks are reproduced by the stochastic NS equation, which implies that the Brownian motion of the fluid particles plays an essential role in the emergence of these peaks in the velocity distribution. We also find that the position of the peak at $V_x(\not=\!0)$ moves with the variation in the strength $D$ of the random Brownian force, which also changes depending on physical parameters such as the kinematic viscosity, boundary velocity and diameter of the plant cells.

preprint2016arXiv

Dependence of the surface tension on the shape of surface boundary

We numerically check that the surface tension of membranes is independent of the shape of surface boundary. The surface tension is calculated by means of the Monte Carlo simulation technique on two types of cylinders made of rubans of size $L_1$ and $L_2$, where the rubans are the same for the projected area and different in the ratio $L_1/L_2$. The difference of the surface tension disappears in the thermodynamic limit in both models of Helfrich-Polyakov and Landau-Ginzburg.

preprint2016arXiv

Finsler geometry modeling of phase separation in multi-component membranes

Finsler geometric surface model is studied as a coarse-grained model for membranes of three-component such as DOPC, DPPC and Cholesterol. To understand the phase separation of liquid ordered (DPPC rich) $L_o$ and the liquid disordered (DOPC rich) $L_d$, we introduce a variable $σ(\in \{1,-1\})$ in the triangulated surface model. We numerically find that there appear two circulars and stripe domains on the surface and that these two morphologies are separated by a phase transition. The morphological change from the one to the other with respect to the variation of the area fraction of $L_o$ is consistent with existing experimental results. This gives us a clear understanding of the origin of the line tension energy, which has been used to understand those morphological changes in the three-component membranes. In addition to these two circulars and stripe domains, raft-like domain and budding domain are also observed, and the corresponding several phase diagrams are obtained. Technical details of the Finsler geometry modeling are also shown.

preprint2015arXiv

Internal phase transition induced by external forces in Finsler geometric model for membranes

We numerically study an anisotropic shape transformation of membranes under external forces for two-dimensional triangulated surfaces on the basis of Finsler geometry. The Finsler metric is defined by using a vector field, which is the tangential component of a three dimensional unit vector $σ$ corresponding to the tilt or some external macromolecules on the surface of disk topology. The sigma model Hamiltonian is assumed for the tangential component of $σ$ with the interaction coefficient $λ$. For large (small) $λ$, the surface becomes oblong (collapsed) at relatively small bending rigidity. For the intermediate $λ$, the surface becomes planar. Conversely, fixing the surface with the boundary of area $A$ or with the two point boundaries of distance $L$, we find that the variable $σ$ changes from random to aligned state with increasing of $A$ or $L$ for the intermediate region of $λ$. This implies that an internal phase transition for $σ$ is triggered not only by the thermal fluctuations but also by external mechanical forces. We also find that the frame (string) tension shows the expected scaling behavior with respect to $A/N$ ($L/N$) at the intermediate region of $A$ ($L$) where the $σ$ configuration changes between the disordered and ordered phases. Moreover, we find that the string tension $γ$ at sufficiently large $λ$ is considerably smaller than that at small $λ$. This phenomenon resembles the so-called soft-elasticity in the liquid crystal elastomer, which is deformed by small external tensile forces.

preprint2015arXiv

Parallel tempering Monte Carlo simulations of spherical fixed-connectivity model for polymerized membranes

We study the first order phase transition of the fixed-connectivity triangulated surface model using the Parallel Tempering Monte Carlo (PTMC) technique on relatively large lattices. From the PTMC results, we find that the transition is considerably stronger than the reported ones predicted by the conventional Metropolis MC (MMC) technique and the flat histogram MC technique. We also confirm that the results of the PTMC on relatively smaller lattices are in good agreement with those known results. This implies that the PTMC is successfully used to simulate the first order phase transitions. The parallel computation in the PTMC is implemented by OpenMP, where the speed of the PTMC on multi-core CPUs is considerably faster than that on the single-core CPUs.

preprint2015arXiv

Surface tension and Laplace pressure in triangulated surface models for membranes without fixed boundary

A Monte Carlo (MC) study is performed to evaluate the surface tension $γ$ of spherical membranes that may be regarded as the models of the lipid layers. We use the canonical surface model defined on the self-avoiding triangulated lattices. The surface tension $γ$ is calculated by keeping the total surface area $A$ constant during the MC simulations. In the evaluation of $γ$, we use $A$ instead of the projected area $A_p$, which is unknown due to the fluctuation of the spherical surface without boundary. The pressure difference ${\itΔ}p $ between the inner and the outer sides of the surface is also calculated by maintaining the enclosed volume constant. Using ${\itΔ}p $ and the Laplace formula, we obtain the tension, which is considered to be equal to the frame tension $τ$ conjugate to $A_p$, and check whether or not $γ$ is consistent with $τ$. We find reasonable consistency between $γ$ and $τ$ in the region of sufficiently large bending rigidity $κ$ or sufficiently large $A/N$. It is also found that $τ$ becomes constant in the limit of $A/N\to \infty$ both in the tethered and fluid surfaces.

preprint2014arXiv

Branched-polymer to inflated transition of self-avoiding fluid surfaces

We study phase transition of self-avoiding fluid surface model on dynamically triangulated lattices using the Monte Carlo simulation technique. We report the continuous transition between the branched polymer and inflated phases at ${\it Δ}p \!=\!0$, where ${\it Δ}p (=\!p_{\rm in}\!-\!p_{\rm out})$ is the pressure difference between the inner and outer sides of the surface. This transition is characterized by almost discontinuous change of the enclosed volume versus the variations of the bending rigidity $κ$ and the pressure difference ${\it Δ}p$. No surface fluctuation transition accompanies this transition up to the surface with the number of vertices $N\!=\!2562$.

preprint2014arXiv

Monte Carlo simulations of Landau-Ginzburg model for membranes

The Landau-Ginzburg (LG) model for membranes is numerically studied on triangulated spheres in ${\bf R}^3$. The LG model is in sharp contrast to the model of Helfrich-Polyakov (HP). The reason for this difference is that the curvature energy of the LG (HP) Hamiltonian is defined by means of the tangential (normal) vector of the surface. For this reason the curvature energy of the LG model includes the in-plane bending or shear energy component, which is not included in the curvature energy of the HP model. From the simulation data, we find that the LG model undergoes a first-order collapse transition. The results of the LG model in the higher dimensional spaces ${\bf R}^d (d>3)$ and on the self-avoiding surfaces in ${\bf R}^3$ are presented and discussed. We also study the David-Guitter (DG) model, which is a variant of the LG model, and find that the DG model undergoes a first-order transition. It is also found that the transition can be observed only on the homogeneous surfaces, which are composed of almost uniform triangles according to the condition that the induced metric ${\partial_a {\bf r}}\cdot{\partial_b {\bf r}}$ is close to $δ_{ab}$.

preprint2013arXiv

Crumpled-to-tubule transition and shape transformations of a model of self-avoiding spherical meshwork

This paper analyzes a new self-avoiding (SA) meshwork model using the canonical Monte Carlo simulation technique on lattices that consist of connection-fixed triangles. The Hamiltonian of this model includes a self-avoiding potential and a pressure term. The model identifies a crumpled-to-tubule (CT) transition between the crumpled and tubular phases. This is a second-order transition, which occurs when the pressure difference between the inner and outer sides of the surface is close to zero. We obtain the Flory swelling exponents $ν_{{\rm R}^2}(=\!D_f/2)$ and $\barν_{\rm v}$ corresponding to the mean square radius of gyration $R_g^2$ and enclosed volume $V$, where $D_f$ is the fractal dimension. The analysis shows that $\barν_{\rm v}$ at the transition is almost identical to the one of the smooth phase of previously reported SA model which has no crumpled phase.

preprint2013arXiv

Monte Carlo studies of a Finsler geometric surface model

This paper presents a new type of surface models constructed on the basis of Finsler geometry. A Finsler metric is defined on the surface by using an underlying vector field, which is an in-plane tilt order. According to the orientation of the vector field, the Finsler length becomes dependent on both position and direction on the surface, and for this reason the parameters such as the surface tension and bending rigidity become anisotropic. To confirm that the model is well-defined, we perform Monte Carlo simulations under several isotropic conditions such as those given by random vector fields. The results are comparable to those of previous simulations of the conventional model. It is also found that a tubular phase appears when the vector field is constant. Moreover, we find that the tilts form the Kosterlitz-Thouless and low temperature configurations, which correspond to two different anisotropic phases such as disk and tubular, in the model in which the tilt variable is assumed to be a dynamical variable. This confirms that the model in this paper may be used as an anisotropic model for membranes.

preprint2012arXiv

In-plane deformation of a triangulated surface model with metric degrees of freedom

Using the canonical Monte Carlo simulation technique, we study a Regge calculus model on triangulated spherical surfaces. The discrete model is statistical mechanically defined with the variables $X$, $g$ and $ρ$, which denote the surface position in ${\bf R}^3$, the metric on a two-dimensional surface $M$ and the surface density of $M$, respectively. The metric $g$ is defined only by using the deficit angle of the triangles in {$M$}. This is in sharp contrast to the conventional Regge calculus model, where {$g$} depends only on the edge length of the triangles. We find that the discrete model in this paper undergoes a phase transition between the smooth spherical phase at $b to infty$ and the crumpled phase at $b to 0$, where $b$ is the bending rigidity. The transition is of first-order and identified with the one observed in the conventional model without the variables $g$ and $ρ$. This implies that the shape transformation transition is not influenced by the metric degrees of freedom. It is also found that the model undergoes a continuous transition of in-plane deformation. This continuous transition is reflected in almost discontinuous changes of the surface area of $M$ and that of $X(M)$, where the surface area of $M$ is conjugate to the density variable $ρ$.

preprint2011arXiv

First-order phase transition of triangulated surfaces on a spherical core

We study an intrinsic curvature model defined on fixed-connectivity triangulated lattices enclosing a spherical core by using the canonical Monte Carlo simulation technique. We find that the model undergoes a discontinuous transition of shape transformation between the smooth state and a collapsed state even when the core radius $R$ is sufficiently large; the transition depends on $R$. The origin of the multitude of transitions is considered to be a degeneracy of the collapsed states. We also find that the Gaussian bond potential $S_1/N$, which is the sum of bond length squares, discontinuously changes at the transition. The discontinuity in $S_1/N$ implies a possibility of large fluctuations of the distance between lipids, or the density of lipids, in biological membranes such as giant vesicles or liposomes enclosing some materials.

preprint2010arXiv

Flat histogram Monte Carlo simulations of triangulated fixed-connectivity surface models

Using the Wang-Landau flat histogram Monte Carlo (FHMC) simulation technique, we were able to study two types of triangulated spherical surface models in which the two-dimensional extrinsic curvature energy is assumed in the Hamiltonian. The Gaussian bond potential is also included in the Hamiltonian of the first model, but it is replaced by a hard-wall potential in the second model. The results presented in this paper are in good agreement with the results previously reported by our group. The transition of surface fluctuations and collapsing transition were studied using the canonical Metropolis Monte Carlo simulation technique and were found to be of the first-order. The results obtained in this paper also show that the FHMC technique can be successfully applied to triangulated surface models. It is non-trivial whether the technique is applicable or not to surface models because the simulations are performed on relatively large surfaces.

preprint2010arXiv

Monte Carlo studies of triangulated spherical surfaces in the two-dimensional space

We numerically study a triangulated surface model in R^2 by taking into account a viewpoint of string model. The models are defined by a mapping X from a two-dimensional surface M to R^2, where the mapping X and the metric g of M are the dynamical variables. The sum over g in the partition function is simulated by the sum over bond lengths and deficit angles by using the Regge calculus technique, and the sum over g is defined to be performed independently of the sum over X. We find that the model undergoes a first-order transition of surface fluctuations, which accompanies a collapsing transition, and that the transitions are reflected in the internal geometry of surface. Fluid surface models are also studied on dynamically triangulated surfaces, and the transitions are found to be of second order. The order of the transition remains unchanged from that of the conventional model defined only by the variable X both in the fixed-connectivity and the fluid models.

preprint2010arXiv

Shape transformations of a model of self-avoiding triangulated surfaces of sphere topology

We study a surface model with a self-avoiding (SA) interaction using the canonical Monte Carlo simulation technique on fixed-connectivity (FC) triangulated lattices of sphere topology. The model is defined by an area energy, a deficit angle energy, and the SA potential. A pressure term is also included in the Hamiltonian. The volume enclosed by the surface is well defined because of the self-avoidance. We focus on whether or not the interaction influences the phase structure of the FC model under two different conditions of pressure ${\it Δ} p$; zero and small negative. The results are compared with the previous results of the self-intersecting model, which has a rich variety of phases; the smooth spherical phase, the tubular phase, the linear phase, and the collapsed phase. We find that the influence of the SA interaction on the multitude of phases is almost negligible except for the evidence that no crumpled surface appears under ${\it Δ} p\=\0$ at least even in the limit of zero bending rigidity $α\to \0$. The Hausdorff dimension is obtained in the limit of $α\to \0$ and compared with previous results of SA models, which are different from the one in this paper.

preprint2009arXiv

Spherical surface models with directors

A triangulated spherical surface model is numerically studied, and it is shown that the model undergoes phase transitions between the smooth phase and the collapsed phase. The model is defined by using a director field, which is assumed to have an interaction with a normal of the surface. The interaction between the directors and the surface maintains the surface shape. The director field is not defined within the two-dimensional differential geometry, and this is in sharp contrast to the conventional surface models, where the surface shape is maintained only by the curvature energies. We also show that the interaction makes the Nambu-Goto model well-defined, where the bond potential is given by the area of triangles; the Nambu-Goto model is well-known as an ill-defined one even when the conventional two-dimensional bending energy is included in the Hamiltonian.

preprint2007arXiv

Phase structure of a spherical surface model on fixed connectivity meshes

An elastic surface model is investigated by using the canonical Monte Carlo simulation technique on triangulated spherical meshes. The model undergoes a first-order collapsing transition and a continuous surface fluctuation transition. The shape of surfaces is maintained by a one-dimensional bending energy, which is defined on the mesh, and no two-dimensional bending energy is included in the Hamiltonian.