Researcher profile

Masaharu Isobe

Masaharu Isobe contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
7works
0followers
3topics
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 graph slice

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

Published work

7 published item(s)

preprint2015arXiv

Hard-sphere melting and crystallization with event-chain Monte Carlo

We simulate crystallization and melting with local Monte Carlo (LMC), event-chain Monte Carlo (ECMC), and with event-driven molecular dynamics (EDMD) in systems with up to one million three-dimensional hard spheres. We illustrate that our implementations of the three algorithms rigorously coincide in their equilibrium properties. We then study nucleation in the NVE ensemble from the fcc crystal into the homogeneous liquid phase and from the liquid into the homogeneous crystal. ECMC and EDMD both approach equilibrium orders of magnitude faster than LMC. ECMC is also notably faster than EDMD, especially for the equilibration into a crystal from a disordered initial condition at high density. ECMC can be trivially implemented for hard-sphere and for soft-sphere potentials, and we suggest possible applications of this algorithm for studying jamming and the physics of glasses, as well as disordered systems.

preprint2013arXiv

Hard-disk equation of state: First-order liquid-hexatic transition in two dimensions with three simulation methods

We report large-scale computer simulations of the hard-disk system at high densities in the region of the melting transition. Our simulations reproduce the equation of state, previously obtained using the event-chain Monte Carlo algorithm, with a massively parallel implementation of the local Monte Carlo method and with event-driven molecular dynamics. We analyze the relative performance of these simulation methods to sample configuration space and approach equilibrium. Our results confirm the first-order nature of the melting phase transition in hard disks. Phase coexistence is visualized for individual configurations via the orientational order parameter field. The analysis of positional order confirms the existence of the hexatic phase.

preprint2012arXiv

Fluctuation-Dissipation Relations for Motions of Center of Mass in Driven Granular Fluids under Gravity

We investigated the validity of fluctuation-dissipation relations in the nonequilibrium stationary state of fluidized granular media under gravity by two independent approaches, based on theory and numerical simulations. A phenomenological Langevin-type theory describing the fluctuation of center of mass height, which was originally constructed for a one-dimensional granular gas on a vibrating bottom plate, was generalized to any dimensionality, even for the case in which the vibrating bottom plate is replaced by a thermal wall. The theory predicts a fluctuation-dissipation relation known to be satisfied at equilibrium, with a modification that replaces the equilibrium temperature by an effective temperature defined by the center of mass kinetic energy. To test the validity of the fluctuation-dissipation relation, we performed extensive and accurate event-driven molecular dynamics simulations for the model system with a thermal wall at the bottom. The power spectrum and response function of the center of mass height were measured and closely compared with theoretical predictions. It is shown that the fluctuation-dissipation relation for the granular system is satisfied, especially in the high-frequency (short time) region, for a wide range of system parameters. Finally, we describe the relationship between systematic deviations in the low-frequency (long time) region and the time scales of the driven granular system.

preprint2012arXiv

Generalized Bond Order Parameters to Characterize Transient Crystals

Higher order parameters in the hard disk fluid are computed to investigate the number, the life time and size of transient crystal nuclei in the pre-freezing phase. The methodology introduces further neighbor shells bond orientational order parameters and coarse-grains the correlation functions needed for the evaluation of the stress autocorrelation function for the viscosity. We successfully reproduce results by the previous collision method for the pair orientational correlation function, but some two orders of magnitude faster. This speed-up allows calculating the time dependent four body orientational correlation between two different pairs of particles as a function of their separation, needed to characterize the size of the transient crystals. The result is that the slow decay of the stress autocorrelation function near freezing is due to a large number of rather small crystal nuclei lasting long enough to lead to the molasses tail.

preprint2010arXiv

A New Approach based on Langevin type Equation for Driven Granular Gas under Gravity

We propose a novel approach based on a Langevin equation for fluctuating motion of the center of mass of granular media fluidized by energy injection from a bottom plate. In this framework, the analytical solution of the Langevin equation is used to derive analytic expressions for several macroscopic quantities and the power spectrum for the center of mass. In order to test our theory, we performed event-driven molecular dynamics simulations for one- and two-dimensional systems. Energy is injected from a vibrating bottom plate in the one-dimensional case and from a thermal wall at the bottom in the two-dimensional case. We found that the theoretical predictions are in good agreement with the results of those simulations under the assumption that the fluctuation-dissipation relation holds in the case of nearly elastic collisions between particles. However, as the inelasticity of the interparticle collisions increases, the power spectrum for the center of mass obtained by the simulations gradually deviates from the prediction of theoretical curve. Connection between this deviation and violation of the fluctuation-dissipation relation is discussed.

preprint2010arXiv

Study of Transient Nuclei near Freezing

The molasses tail in dense hard core fluids is investigated by extensive event-driven molecular dynamics simulation through the orientational autocorrelation functions. Near the fluid-solid phase transition, there exist three regimes in the relaxation of the pair orientational autocorrelation function, namely the kinetic, molasses (stretched exponential), and diffusional power decay. The density dependence of both the molasses and diffusional power regimes are evaluated and the latter compares with theoretical predictions in three dimensions. The largest cluster at the freezing density of only a few sphere diameter in size persist for only about 30 picoseconds (~ 2.8 x 10^{-11}[s]). The most striking observation through the bond orientatinal order parameter is the dramatic increase of the cluster size as the freezing density is approached.

preprint2008arXiv

A Langevin Approach to One-Dimensional Granular Media Fluidized by Vibrations

We present a Langevin approach to describe the steady-state dynamics of one-dimensional granular media fluidized by a vibrating bottom plate. We adopt a linear Langevin equation to describe the motion of the center of mass. Within this framework, we derive analytical expressions for several macroscopic quantities. We also predict the power spectrum for the height of the center of mass. We find good agreement between our theoretical predictions and extensive event-driven molecular dynamics simulations.