Source author record

S. Mineshige

S. Mineshige 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

21works
4topics
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

21 published item(s)

preprint2019arXiv

Similar shot profile morphology of fast variability in cataclysmic variable, X-ray binary and blazar; the MV Lyr case

The cataclysmic variable MV Lyr was present in the Kepler field yielding a light curve with the duration of almost 1500 days with 60 second cadence. Such high quality data of this nova-like system with obvious fast optical variability show multicomponent power density spectra. Our goal is to study the light curve from different point of view, and perform a shot profile analysis. We search for characteristics not discovered with standard power density spectrum based methods. The shot profile method identifies individual shots in the light curve, and averages them in order to get all substructures with typical time scales. We also tested the robustness of our analysis using simple shot noise model. We obtained mean profiles with multicomponent features. The shot profile method distinguishes substructures with similar time scales which appear as a single degenerate feature in power density spectra. Furthermore, this method yields the identification of another high frequency component in the power density spectra of Kepler and XMM-Newton data not detected so far. Moreover, we found side-lobes accompanied with the central spike, making the profile very similar to another Kepler data of blazar W2R 1926+42, and Ginga data of Cyg X-1. All three objects show similar time scale ratios of the rising vs. declining part of the central spikes, while the two binaries have also similar rising profiles of the shots described by a power-law function. The similarity of both binary shot profiles suggests that the shots originate from the same origin, e.g. aperiodic mass accretion in the accretion disc. Moreover, the similarity with the blazar may imply that the ejection fluctuations in the blazar jet are connected to accretion fluctuations driving the variability in binaries. This points out to connection between jet and the accretion disc.

preprint2015arXiv

Shot Analysis of Kepler Blazar W2R 1926+42

Blazars show rapid and violent variabilities, which timescale are often less than a day. We studied intraday variations by applying a shot analysis technique to Kepler monitoring of blazar W2R 1926+42 in Quarter 14. We obtained a mean profile calculated from 195 rapid variations. The mean profile shows three components; one is a sharp structure distributed within $\pm$0.1 day of the peak, and two slow-varying components. This spiky-peak component reflects features of rapid variations directly. The profile of peak component shows an exponential rise and decay of which timescales are different, 0.0416 and 0.0588 day respectively. This component is too sharp to represent a standard function which is often used to express blazar variations. This asymmetric profile at the peak is difficult to be explained by a simple variation of the Doppler factor by changing a geometry of the emitting region. This result indicates that intraday variations arise from a production of high-energy accelerated particles in the jet.

preprint2014arXiv

Radiation Hydrodynamic Instability in Plane-Parallel, Super-Eddington Atmosphere: A Mechanism for Clump Formation

In order to understand the physical processes underlying clump formation in outflow from supercritical accretion flow, we performed two-dimensional radiation hydrodynamic (RHD) simulations. We focus our discussion on the nature of RHD instability in marginally optically thick, plane-parallel, super-Eddington atmosphere. Initially we set two-layered atmosphere with a density contrast of 100 exposed to strong, upward continuum-radiation force; the lower layer is denser than the upper one, condition for an RHD instability. We assume non-zero but negligible gravitational force, compared with the radiation force. We find that short wavelength perturbations first grow, followed by growth of longer wavelength patterns, which lead to the formation of clumpy structure. The typical size of clumps (clouds) corresponds to about one optical depth. An anti-correlation between the radiation pressure and the gas pressure is confirmed: this anti-correlation provides a damping mechanism of longer wavelength perturbations than the typical clump size. Matter and radiation energy densities are correlated. These features are exactly what we found in the radiation-magnetohydrodynamic (radiation-MHD) simulations of supercritical outflow.

preprint2014arXiv

Rms-flux relation and fast optical variability simulations of the nova-like system MV Lyr

The stochastic variability (flickering) of the nova-like system (subclass of cataclysmic variable) MV Lyr yields a complicated power density spectrum with four break frequencies. Scaringi et al. (2012) analysed high-cadence Kepler data of MV Lyr, taken almost continuously over 600 days, giving the unique opportunity to study multicomponent Power Density Spectra (PDS) over a wide frequency range. We modelled this variability with our statistical model based on disc angular momentum transport via discrete turbulent bodies with an exponential distribution of the dimension scale. Two different models were used, a full disc (developed from the white dwarf to the outer radius of ~10^10 cm) and a radially thin disc (a ring at a distance of ~10^10 cm from the white dwarf) that imitates an outer disc rim. We succeed in explaining the two lowest observed break frequencies assuming typical values for a disc radius of 0.5 and 0.9 times the primary Roche lobe and an alpha parameter of 0.1-0.4. The highest observed break frequency was also modelled, but with a rather small accretion disc with radius of 0.3 times the primary Roche lobe and a high alpha value of 0.9 consistent with previous findings by Scaringi (2014). Furthermore, the simulated light curves exhibit the typical linear rms-flux proportionality linear relation and the typical log-normal flux distribution. As the turbulent process is generating fluctuations in mass accretion that propagate through the disc, this confirms the general knowledge that the typical rms-flux relation is mainly generated by these fluctuations. In general a higher rms is generated by a larger amount of superposed flares which is compatible with a higher mass accretion rate expressed by a larger flux.

preprint2013arXiv

Clumpy Outflow from Supercritical Accretion Flows

Significant fraction of matter in supercritical (or super-Eddington) accretion flow is blown away by radiation force, thus forming outflows, however, the properties of such radiation-driven outflows have been poorly understood. We have performed global two-dimensional radiaion-magnetohydrodynamic simulations of supercritical accretion flow onto a black hole with 10 or 10^8 solar masses in a large simulation box of 514 r_S x 514 r_S (with r_S being the Schwarzschild radius). We confirm that uncollimated outflows with velocities of 10 percents of the speed of light emerge from the innermost part of the accretion flow over wide angles of 10 - 50 degree from the disk rotation axis. Importantly, the outflows exhibit clumpy structure above heights of ~ 250 r_S. The typical size of the clumps is ~ 10 r_S, which corresponds to one optical depth, and their shapes are slightly elongated along the outflow direction. Since clumps start to form in the layer above which (upward) radiation force overcomes (downward) gravity force, Rayleigh-Taylor instability seems to be of primary cause. In addition, a radiation hydrodynamic instability, which arises when radiation funnels through radiation-pressure supported atmosphere, may also help forming clumps of one optical depth. Magnetic photon bubble instability seems not to be essential, since similar clumpy outflow structure is obtained in non-magnetic radiation-hydrodynamic simulations. Since the spatial covering factor of the clumps is estimated to be ~ 0.3 and since they are marginally optically thick, they will explain at least some of rapid light variations of active galactic nuclei. We further discuss a possibility of producing broad-line clouds by the clumpy outflow.

preprint2013arXiv

Highly Ionized Fe-K Absorption Line from Cygnus X-1 in the High/Soft State Observed with Suzaku

We present observations of a transient He-like Fe K alpha absorption line in Suzaku observations of the black hole binary Cygnus X-1 on 2011 October 5 near superior conjunction during the high/soft state, which enable us to map the full evolution from the start and the end of the episodic accretion phenomena or dips for the first time. We model the X-ray spectra during the event and trace their evolution. The absorption line is rather weak in the first half of the observation, but instantly deepens for ~10 ks, and weakens thereafter. The overall change in equivalent width is a factor of ~3, peaking at an orbital phase of ~0.08. This is evidence that the companion stellar wind feeding the black hole is clumpy. By analyzing the line with a Voigt profile, it is found to be consistent with a slightly redshifted Fe XXV transition, or possibly a mixture of several species less ionized than Fe XXV. The data may be explained by a clump located at a distance of ~10^(10-12) cm with a density of ~10^((-13)-(-11)) g cm^-3, which accretes onto and/or transits the line-of-sight to the black hole, causing an instant decrease in the observed degree of the ionization and/or an increase in density of the accreting matter. Continued monitoring for individual events with future X-ray calorimeter missions such as ASTRO-H and AXSIO will allow us to map out the accretion environment in detail and how it changes between the various accretion states.

preprint2013arXiv

Inner Disc Obscuration in GRS 1915+105 Based on Relativistic Slim Disc Model

We study the observational signatures of the relativistic slim disc of 10 M_sun black hole, in a wide range of mass accretion rate, mdot, dimensionless spin parameter, a_ast, and viewing angle, i. In general, the innermost temperature, T_in increases with the increase of i for a fixed value of mdot and a_ast, due to the Doppler effect. However, for i > 50 and mdot > mdot_turn, T_in starts to decrease with the increase of mdot. This is a result of self-obscuration -- the radiation from the innermost hot part of the disc is blocked by the surrounding cooler part. The value of mdot_turn and the corresponding luminosities depend on a_ast and i. Such obscuration effects cause an interesting behavior on the disc luminosity (L_disc) -- T_in plane for high inclinations. In addition to the standard-disc branch which appears below mdot_turn and which obeys L_disc propto T_in^4 -relation, another branch above mdot_turn, which is nearly horizontal, may be observed at luminosities close to the Eddington luminosity. We show that these features are likely observed in a Galactic X-ray source, GRS 1915+105. We support a high spin parameter (a_ast > 0.9) for GRS 1915+105 since otherwise the high value of T_in and small size of the emitting region (r_in < 1r_S) cannot be explained.

preprint2013arXiv

Rapid Spectral Changes of Cygnus X-1 in the Low/Hard State with Suzaku

Rapid spectral changes in the hard X-ray on a time scale down to ~0.1 s are studied by applying "shot analysis" technique to the Suzaku observations of the black hole binary Cygnus X-1, performed on 2008 April 18 during the low/hard state. We successfully obtained the shot profiles covering 10--200 keV with the Suzaku HXD-PIN and HXD-GSO detector. It is notable that the 100-200 keV shot profile is acquired for the first time owing to the HXD-GSO detector. The intensity changes in a time-symmetric way, though the hardness does in a time-asymmetric way. When the shot-phase-resolved spectra are quantified with the Compton model, the Compton y-parameter and the electron temperature are found to decrease gradually through the rising phase of the shot, while the optical depth appears to increase. All the parameters return to their time-averaged values immediately within 0.1 s past the shot peak. We have not only confirmed this feature previously found in energies below ~60 keV, but also found that the spectral change is more prominent in energies above ~100 keV, implying the existence of some instant mechanism for direct entropy production. We discuss possible interpretations on the rapid spectral changes in the hard X-ray band.

preprint2013arXiv

Resolving different sources of fast X-ray variability of the dwarf nova RU Peg in quiescence

We analysed an X-ray light curve of the dwarf nova RU Peg taken by XMM-Newton with a duration of 46300 s. The power density spectrum has a complicated shape with two red noise and two white noise components, indicating the presence of two turbulent regions. We developed a statistical "toy model" to study light curves with variability produced by an unstable turbulent accretion flow from the inner disc. Our results are consistent with a disc truncation radius maximally 0.8 x 10^9 cm. We found that any fluctuation of the viscous mass accretion at the inner disc are visible as UV and X-ray variations with the same break frequency in the power density spectrum. This process is generating low frequency variability. A second break suggests the presence of a faster X-ray variability component which must be generated by another process likely localised between the inner disc and the white dwarf.

preprint2012arXiv

The Structure and Spectral Features of a Thin Disk and Evaporation-Fed Corona in High-Luminosity AGNs

We investigate the accretion process in high-luminosity AGNs (HLAGNs) in the scenario of the disk evaporation model. Based on this model, the thin disk can extend down to the innermost stable circular orbit (ISCO) at accretion rates higher than $0.02\dot{M}_{\rm Edd}$; while the corona is weak since part of the coronal gas is cooled by strong inverse Compton scattering of the disk photons. This implies that the corona cannot produce as strong X-ray radiation as observed in HLAGNs with large Eddington ratio. In addition to the viscous heating, other heating to the corona is necessary to interpret HLAGN. In this paper, we assume that a part of accretion energy released in the disk is transported into the corona, heating up the electrons and thereby radiated away. We for the first time, compute the corona structure with additional heating, taking fully into account the mass supply to the corona and find that the corona could indeed survive at higher accretion rates and its radiation power increases. The spectra composed of bremsstrahlung and Compton radiation are also calculated. Our calculations show that the Compton dominated spectrum becomes harder with the increase of energy fraction ($f$) liberating in the corona, and the photon index for hard X-ray($2-10 \rm keV$) is $2.2 < Γ< 2.7 $. We discuss possible heating mechanisms for the corona. Combining the energy fraction transported to the corona with the accretion rate by magnetic heating, we find that the hard X-ray spectrum becomes steeper at larger accretion rate and the bolometric correction factor ($L_{\rm bol}/L_{\rm 2-10keV}$) increases with increasing accretion rate for $f<8/35$, which is roughly consistent with the observational results.

preprint2010arXiv

A Novel Jet Model: Magnetically Collimated, Radiation-Pressure Driven Jet

Relativistic jets from compact objects are ubiquitous phenomena in the Unvierse, but their driving mechanism has been an enigmatic issue over many decades. Two basic models have been extensively discussed: magnetohydrodynamic (MHD) jets and radiation-hydrodynamic (RHD) jets. Currently, the former is more widely accepted, since magnetic field is expected to provide both the acceleration and collimation mechanisms, whereas radiation field cannot collimate outflow. Here, we propose a new type of jets, radiation-magnetohydrodynamic (RMHD) jets, based on our global RMHD simulation of luminous accretion flow onto a black hole shining above the Eddington luminosity. The RMHD jet can be accelerated up to the relativistic speed by the radiation-pressure force and is collimated by the Lorentz force of a magnetic tower, inflated magnetic structure made by toroidal magnetic field lines accumulated around the black hole, though radiation energy greatly dominates over magnetic energy. This magnetic tower is collimated by a geometrically thick accretion flow supported by radiation-pressure force. This type of jet may explain relativistic jets from Galactic microquasars, appearing at high luminosities.

preprint2010arXiv

Probing the Peculiar Behavior of GRS 1915+105 at Near-Eddington Luminosity

To understand the nature of supercritical accretion, we systematically analyze the {\it RXTE}/PCA data of GRS 1915+105 in its quasi-steady states, by choosing data with small variability during 1999 -- 2000. We apply a multicolor disk plus a thermal Comptonization model and take into consideration accurate interstellar absorption, a reflection component, and absorption features from the disk wind self-consistently. There is a strong correlation between the inner disk temperature and the fraction of the disk component. Most of the Comptonization-dominated spectra show $T_{\rm in} \sim 1$ keV with a high electron temperature of $>10$ keV, which may correspond to the very high state in canonical black hole X-ray binaries (BHBs). By contrast, the disk-dominated spectra have $T_{\rm in} \sim 2$ keV with a low temperature ($<10$ keV) and optically thick Comptonization, and show two separate branches in the ($L$ -- $T_{\rm in}$) diagram. The lower branch clearly follows the $L \propto T_{\rm in}^4$-track. Furthermore, applying the extended disk blackbody model, we find that 9 out of 12 datasets with disk luminosity above $0.3 L_{\rm E}$ prefer a flatter temperature gradient than that in the standard disk ($p < 0.7$). We interpret that, in the lower branch, the disk extends down to the innermost stable circular orbit, and is most probably in the slim disk state. A rapidly spinning black hole can explain both the lack of the $L \propto T_{\rm in}^2$-track and a high value of spectral hardening factor ($\sim 4$) that would be required for a non-rotating black hole. The spectra in the upper branch are consistent with the picture of a truncated disk with low temperature Comptonization. This state is uniquely observed from GRS 1915+105 among BHBs, which may be present at near-Eddington luminosity.

preprint2009arXiv

Viscous propagation of mass flow variability in accretion discs

We study mass flow rate through a disc resulting from a varying mass supply rate. Variable mass supply rate occurs, e.g., during disc state transitions, and in interacting eccentric binaries. It is, however, damped by the viscosity of the disc. Here, we calculate this damping in detail. We derive an analytical description of the propagation of the flow rate using the solution of Lynden-Bell & Pringle, in which the disc is assumed to extend to infinity. In particular, we derive the accretion-rate Green's function, and its Fourier transform, which gives the fractional damping at a given variability frequency. We then compare this model to that of a finite disc with the mass supply at its outer edge. We find significant differences with respect to the infinite disc solution, which we find to overestimate the viscous damping. In particular, the asymptotic form of the Green's function is power-law for the infinite disc and exponential for the finite one. We then find a simple fitting form for the latter, and also calculate its Fourier transform. In general, the damping becomes very strong when the viscous time at the outer edge of the disc becomes longer than the modulation time scale. We apply our results to a number of astrophysical systems. We find the effect is much stronger in low-mass X-ray binaries, where the disc size is comparable to that of the Roche lobe, than in high-mass binaries, where the wind-fed disc can have a much smaller size.

preprint2009arXiv

X-ray spectral variability in the ultraluminous X-ray source Holmberg IX X-1

We use XMM-Newton and Swift data to study spectral variability in the ultraluminous X-ray source (ULX), Holmberg IX X-1. The source luminosity varies by a factor 3-4, giving rise to corresponding spectral changes which are significant, but subtle, and not well tracked by a simple hardness ratio. Instead, we co-add the Swift data in intensity bins and do full spectral fitting with disc plus thermal Comptonisation models. All the data are well-fitted by a low temperature, optically thick Comptonising corona, and the variability can be roughly characterised by decreasing temperature and increasing optical depth as the source becomes brighter, as expected if the corona is becoming progressively mass loaded by material blown off the super-Eddington inner disc. This variability behaviour is seen in other ULX which have similar spectra, but is opposite to the trend seen in ULX with much softer spectra. This supports the idea that there are two distinct physical regimes in ULXs, where the spectra go from being dominated by a disc-corona to being dominated by a wind.

preprint2006arXiv

Black hole accretion: theoretical limits and observational implications

Recently, the issue of the role of the Eddington limit in accretion discs became a matter of debate. While the classical (spherical) Eddington limit is certainly an over-simplification, it is not really clear how to treat it in a flattened structure like an accretion disc. We calculate the critical accretion rates and resulting disc luminosities for various disc models corresponding to the classical Eddington limit by equating the attractive and repulsive forces locally. We also discuss the observational appearance of such highly accreting systems by analyzing their spectral energy distributions. Our calculations indicate that the allowed mass accretion rates differ considerably from what one expects by applying the Eddington limit in its classical form, while the luminosities only weakly exceed their classical equivalent. Depending on the orientation of the disc relative to the observer, mild relativistic beaming turns out to have an important influence on the disc spectra. Thus, possible super-Eddington accretion, combined with mild relativistic beaming, supports the idea that ultraluminous X-ray sources host stellar mass black holes and accounts partially for the observed high temperatures of these objects.

preprint2002arXiv

Gravitational-Wave Radiation from Magnetized Accretion Disks

The detectability of gravitational wave (GW) radiation from accretion disks is discussed based on various astrophysical contexts. In order to emit GW radiation, the disk shape should lose axial symmetry. We point out that a significant deformation is plausible in non-radiative hot accretion disks because of enhanced magnetic activity, whereas it is unlikely for standard-type cool disks. We have analyzed the 3D magnetohydrodynamical (MHD) simulation data of magnetized accretion flow, finding non-axisymmetric density patterns. The corresponding ellipticity is $ε\sim 0.01$. The expected time variations of GW radiation are overall chaotic, but there is a hint of quasi-periodicity. GW radiation has no interesting consequence, however, in the case of close binaries, because of very tiny disk masses. GW radiation is not significant, either, for AGN because of very slow rotation velocities. The most promising case can be found in gamma-ray bursts or supernovae, in which a massive torus (or disk) with a solar mass or so may be formed around a stellar-mass compact object as the result of a merger of compact objects, or by the fallback of exploded material towards the center in a supernova. Although much more intense GW radiation is expected before the formation of the torus, the detection of GW radiation in the subsequent accretion phase is of great importance, since it will provide a good probe to investigating their central engines.

preprint2002arXiv

Log-Normal Distributions in Cygnus X-1: Possible Physical Link with Gamma-Ray Bursts and Blazars

Prompted by recent discoveries of log-normal distributions in gamma-ray/X-ray temporal variabilities of gamma-ray bursts (GRBs) and blazars, we re-examine the X-ray variability of Cygnus X-1 in the hard/low state using Ginga data in 1990. It was previously reported that the distributions of the time intervals between X-ray shots (flares) deviated from the Poisson distributions at short time intervals: the occurrence of shots tended to be suppressed for several seconds before and/or after a shot event. Detailed analyses show that this deviation is larger for shots with larger peaks, and that the time-interval distributions for large shots approach log-normal distributions with a peak at the interval of 7--8 s. Furthermore, we also show that the peak-intensity distribution for the shots is consistent with the log-normal distribution, though no typical peak intensity can be seen. This might indicate the presence of a physical link connecting the physics of black hole accretion flow and that of jet/GRB formation.

preprint2001arXiv

Theoretical Models of Multi-waveband QSO Luminosity Functions

Cosmological evolution of the QSO luminosity functions (LFs) at NIR/optical/X-ray bands for 1.3 < z < 3.5 is investigated based on the realistic QSO spectra. The accretion-disk theory predicts that although QSO luminosities only depend on mass-accretion rate, \Mdot, QSO spectra have a dependence on black-hole mass, M_{BH}, as well. The smaller M_{BH} is and/or the larger \Mdot is, the harder becomes the QSO NIR/optical/UV spectrum. We model disk spectra which can reproduce these features and calculated LFs for redshift z ~ 3 with the assumption of new-born QSOs being shining at the Eddington luminosity. The main results are: (i) the observed LFs at optical and X-rays can be simultaneously reproduced. (ii) LFs at optical and X-ray bands are not sensitive to M_{BH}, while LFs at NIR bands are; about one order of magnitude difference is expected in volume number densities at L_{I, J} ~ 10^{46} erg s^{-1} between the case that all QSOs would have the same spectral shape as that of M_{BH} = 10^{9} M_sun and the case with M_{BH} = 10^{11} M_sun. (iii) The resultant LFs at NIR are dominated by 10^{7} M_sun black-holes at L_{I, J} ~ 10^{44} erg s^{-1}, and by 10^{11} M_sun black-holes at L_{I, J} \~ 10^{46} erg s^{-1}. Future infrared observations from space(e.g.NGST) will probe cosmological evolution of black hole masses. For redshift z < 3, on the other hand, the observed optical/X-ray LFs can be fitted, if the initial QSO luminosity L_0 is below the Eddington luminosity. Interestingly, the best fitting values of l = L_0/L_{Edd} are different in B- and X-ray bands; l_B ~ 2.5 l_X. The reason for this discrepancy is briefly discussed.

preprint1999arXiv

Temporal 1/f^αFluctuations from Fractal Magnetic Fields in Black Hole Accretion Flow

Rapid fluctuation with a frequency dependence of $1/f^α$ (with $α\simeq 1 - 2$) is characteristic of radiation from black-hole objects. Its origin remains poorly understood. We examine the three-dimensional magnetohydrodynamical (MHD) simulation data, finding that a magnetized accretion disk exhibits both $1/f^α$ fluctuation (with $α\simeq 2$) and a fractal magnetic structure (with the fractal dimension of $D \sim 1.9$). The fractal field configuration leads reconnection events with a variety of released energy and of duration, thereby producing $1/f^α$ fluctuations.

preprint1997arXiv

The influence of Galactic wind upon the star formation histories of Local Group galaxies

We examine the possibility that ram pressure exerted by the galactic wind from the Galaxy could have stripped gas from the Local Group dwarf galaxies, thereby affecting their star formation histories. Whether gas stripping occurs or not depends on the relative magnitudes of two counteracting forces acting on gas in a dwarf galaxy: ram pressure force by the wind and the gravitational binding force by the dwarf galaxy itself. We suggest that the galactic wind could have stripped gas in a dwarf galaxy located within the distance of $R_{c}\simeq 120(r_{s}/1 {kpc})^{3/2} ({\cal E}_{b}/10^{50} {erg})^{-1/2}$ kpc (where $r_{s}$ is the surface radius and ${\cal E}_{b}$ is the total binding energy of the dwarf galaxy, respectively) from the Galaxy within a timescale of Gyr, thereby preventing star formation there. Our result based on this Galactic wind model explains the recent observation that dwarfs located close to the Galaxy experienced star formation only in the early phase of their lifetimes, whereas distant dwarfs are still undergoing star formation. The present star formation in the Large Magellanic Cloud can also be explained through our Galactic wind model.

preprint1995arXiv

Subcluster merging in clusters of galaxies and the cosmological density parameter

Recent X-ray observations have established that collisions between subclusters of galaxies are rather common phenomena. Prompted by such observations, we have performed N-body simulations of two equal-mass subclusters of galaxies, which are going to merge. We first have confirmed that only a part of kinetic energy associated with the relative motion of two subclusters is converted to internal energy of each subcluster at the first encounter. We performed simulations for a variety of dark-matter (DM) distributions, and find that the time scale for washing out the double peak structures after the first encounter strongly depends on the distribution of the density and velocity distributions of the DM. It takes longer when the DM is spatially extended. According to our calculation it takes more than $4 imes 10^9$yr after the first encounter until the density contour shows only a single peak. To explain the high fraction of the clusters with substructures among nearby clusters, Richstone,Loeb Turner(1992; hereafter RLT) required recent($$ $10^9 h^{-1}$yr) cluster formation. However, our results show that the timescale for subcluster merging is still uncertain and possibly much longer than the time scale assumed by RLT. Caution should be exercised when concluding that the density of the universe is high by using RLT's method.