BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Masayuki Umemura

Masayuki Umemura appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-ph.IMastro-ph.HEastro-phphysics.comp-phastro-ph.EPastro-ph.SRBiological Physicsphysics.med-ph

Catalog footprint

What is connected

28works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

28 published item(s)

preprint2022arXiv

TRINITY: a three-dimensional time-dependent radiative transfer code for in-vivo near-infrared imaging

We develop a new three-dimensional time-dependent radiative transfer code, TRINITY (Time-dependent Radiative transfer In Near-Infrared TomographY), for in-vivo diffuse optical tomography (DOT). The simulation code is based on the design of long radiation rays connecting boundaries of a computational domain, which allows us to calculate light propagation with little numerical diffusion. We parallelize the code with Message Passing Interface (MPI) using the domain decomposition technique and confirm the high parallelization efficiency, so that simulations with a spatial resolution of $\sim 1$ millimeter can be performed in practical time. As a first application, we study the light propagation for a pulse collimated within $θ\sim 15^\circ$ in a phantom, which is a uniform medium made of polyurethane mimicking biological tissue. We show that the pulse spreads in all forward directions over $\sim$ a few millimeters due to the multiple scattering process of photons. Our simulations successfully reproduce the time-resolved signals measured with eight detectors for the phantom. We also introduce the effects of reflection and refraction at the boundary of medium with a different refractive index and demonstrate the faster propagation of photons in an air hole that is an analogue for the respiratory tract.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

High circular polarization of near infrared light induced by micron-size dust grains

We explore the induction of circular polarization (CP) of near-infrared light in star-forming regions using three-dimensional radiative transfer calculations. The simulations trace the change of Stokes parameters at each scattering/absorption process in a dusty gas slab composed of aligned grains. We find that the CP degree enlarges significantly according as the size of dust grains increases and exceeds $\sim 20$ percent for micron-size grains. Therefore, if micron-size grains are dominant in a dusty gas slab, the high CP observed around luminous young stellar objects can be accounted for. The distributions of CP show the asymmetric quadrupole patters regardless of the grain sizes. Also, we find that the CP degree depends on the relative position of a dusty gas slab. If a dusty gas slab is located behind a star-forming region, the CP reaches $\sim 60$ percent in the case of $1.0~{\rm μm}$ dust grains. Hence, we suggest that the observed variety of CP maps can be explained by different size distributions of dust grains and the configuration of aligned grains.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Discrimination of heavy elements originating from Pop III stars in z = 3 intergalactic medium

We investigate the distribution of metals in the cosmological volume at $z\sim3$, in particular, provided by massive population III (Pop III) stars using a cosmological $N$-body simulation in which a model of Pop III star formation is implemented. Owing to the simulation, we can choose minihaloes where Pop III star formation occurs at $z>10$ and obtain the spatial distribution of the metals at lower-redshifts. To evaluate the amount of heavy elements provided by Pop III stars, we consider metal yield of pair-instability or core-collapse supernovae (SNe) explosions of massive stars. By comparing our results to the Illustris-1 simulation, we find that heavy elements provided by Pop III stars often dominate those from galaxies in low density regions. The median value of the volume averaged metallicity is $Z\sim 10^{-4.5 - -2} Z_{\odot}$ at the regions. Spectroscopic observations with the next generation telescopes are expected to detect the metals imprinted on quasar spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

MAGI: many-component galaxy initialiser

Providing initial conditions is an essential procedure for numerical simulations of galaxies. The initial conditions for idealised individual galaxies in $N$-body simulations should resemble observed galaxies and be dynamically stable for time scales much longer than their characteristic dynamical times. However, generating a galaxy model ab initio as a system in dynamical equilibrium is a difficult task, since a galaxy contains several components, including a bulge, disc, and halo. Moreover, it is desirable that the initial-condition generator be fast and easy to use. We have now developed an initial-condition generator for galactic $N$-body simulations that satisfies these requirements. The developed generator adopts a distribution-function-based method, and it supports various kinds of density models, including custom-tabulated inputs and the presence of more than one disc. We tested the dynamical stability of systems generated by our code, representing early- and late-type galaxies, with $N=$~2,097,152 and 8,388,608 particles, respectively, and we found that the model galaxies maintain their initial distributions for at least 1~Gyr. The execution times required to generate the two models were $8.5$ and $221.7$ seconds, respectively, which is negligible compared to typical execution times for $N$-body simulations. The code is provided as open-source software and is publicly and freely available at \url{https://bitbucket.org/ymiki/magi}.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Formation of globular clusters induced by external ultraviolet radiation II: Three-dimensional radiation hydrodynamics simulations

We explore the possibility of the formation of globular clusters under ultraviolet (UV) background radiation. One-dimensional spherical symmetric radiation hydrodynamics (RHD) simulations by Hasegawa et al. have demonstrated that the collapse of low-mass (10^6-10^7 solar masses) gas clouds exposed to intense UV radiation can lead to the formation of compact star clusters like globular clusters (GCs) if gas clouds contract with supersonic infall velocities. However, three-dimensional effects, such as the anisotropy of background radiation and the inhomogeneity in gas clouds, have not been studied so far. In this paper, we perform three-dimensional RHD simulations in a semi-cosmological context, and reconsider the formation of compact star clusters in strong UV radiation fields. As a result, we find that although anisotropic radiation fields bring an elongated shadow of neutral gas, almost spherical compact star clusters can be procreated from a "supersonic infall" cloud, since photo-dissociating radiation suppresses the formation of hydrogen molecules in the shadowed regions and the regions are compressed by UV heated ambient gas. The properties of resultant star clusters match those of GCs. On the other hand, in weak UV radiation fields, dark matter-dominated star clusters with low stellar density form due to the self-shielding effect as well as the positive feedback by ionizing photons. Thus, we conclude that the "supersonic infall" under a strong UV background is a potential mechanism to form GCs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

GOTHIC: Gravitational oct-tree code accelerated by hierarchical time step controlling

The tree method is a widely implemented algorithm for collisionless $N$-body simulations in astrophysics well suited for GPU(s). Adopting hierarchical time stepping can accelerate $N$-body simulations; however, it is infrequently implemented and its potential remains untested in GPU implementations. We have developed a Gravitational Oct-Tree code accelerated by HIerarchical time step Controlling named \texttt{GOTHIC}, which adopts both the tree method and the hierarchical time step. The code adopts some adaptive optimizations by monitoring the execution time of each function on-the-fly and minimizes the time-to-solution by balancing the measured time of multiple functions. Results of performance measurements with realistic particle distribution performed on NVIDIA Tesla M2090, K20X, and GeForce GTX TITAN X, which are representative GPUs of the Fermi, Kepler, and Maxwell generation of GPUs, show that the hierarchical time step achieves a speedup by a factor of around 3--5 times compared to the shared time step. The measured elapsed time per step of \texttt{GOTHIC} is 0.30~s or 0.44~s on GTX TITAN X when the particle distribution represents the Andromeda galaxy or the NFW sphere, respectively, with $2^{24} =$~16,777,216 particles. The averaged performance of the code corresponds to 10--30\% of the theoretical single precision peak performance of the GPU.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Mergers of accreting stellar-mass black holes

We present post-Newtonian $N$-body simulations on mergers of accreting stellar-mass black holes (BHs), where such general relativistic effects as the pericenter shift and gravitational wave (GW) emission are taken into consideration. The attention is concentrated on the effects of the dynamical friction and the Hoyle-Lyttleton mass accretion by ambient gas. We consider a system composed of ten BHs with initial mass of $30~M_\odot$. As a result, we show that mergers of accreting stellar-mass BHs are classified into four types: a gas drag-driven, an interplay-driven, a three body-driven, or an accretion-driven merger. We find that BH mergers proceed before significant mass accretion, even if the accretion rate is $\sim10$ Eddington accretion rate, and then all BHs can merge into one heavy BH. Using the simulation results for a wide range of parameters, we derive a critical accretion rate ($\dot{m}_{\rm c}$), below which the BH growth is promoted faster by mergers. Also, it is found that the effect of the recoil by the GW emission can reduce $\dot{m}_{\rm c}$ especially in gas number density higher than $10^8~{\rm cm}^{-3}$, and enhance the escape probability of merged BHs. Very recently, a gravitational wave event, GW150914, as a result of the merger of a $\sim 30~M_\odot$ BH binary has been detected (Abbott et al. 2016). Based on the present simulations, the BH merger in GW150914 is likely to be driven by three-body encounters accompanied by a few $M_\odot$ of gas accretion, in high-density environments like dense interstellar clouds or galactic nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Subparsec-scale dynamics of a dusty gas disk exposed to anisotropic AGN radiation with frequency-dependent radiative transfer

We explore the gas dynamics near the dust sublimation radius of active galactic nucleus (AGN). For the purpose, we perform axisymmetric radiation hydrodynamic simulations of a dusty gas disk of radius $\approx 1\,\mathrm{pc}$ around a supermassive black hole of mass $10^{7}\,\mathrm{M_{\odot}}$ taking into account (1) anisotropic radiation of accretion disk, (2) X-ray heating by corona, (3) radiative transfer of infrared (IR) photons reemitted by dust, (4) frequency dependency of direct and IR radiations, and (5) separate temperatures for gas and dust. As a result, we find that for Eddington ratio $\approx 0.77$, a nearly neutral, dense ($\approx 10^{6\operatorname{-}8}\;\mathrm{cm^{-3}}$), geometrically-thin ($h/r<0.06$) disk forms with a high velocity ($\approx 200 \sim 3000\;\mathrm{km/s}$) dusty outflow launched from the disk surface. The disk temperature is determined by the balance between X-ray heating and various cooling, and the disk is almost supported by thermal pressure. Contrary to \citet{krolik07:_agn}, the radiation pressure by IR photons is not effective to thicken the disk, but rather compresses it. Thus, it seems difficult for a radiation-supported, geometrically-thick, obscuring torus to form near the dust sublimation radius as far as the Eddington ratio is high ($\sim 1$). The mass outflow rate is $0.05$-$0.1\;\mathrm{M_{\odot}}/\mathrm{yr}$ and the column density of the outflow is $N_{\mathrm{H}}\lesssim 10^{21}\;\mathrm{cm^{-2}}$. To explain observed type-II AGN fraction, it is required that outflow gas is extended to larger radii ($r\gtrsim 10\;\mathrm{pc}$) or that a denser dusty wind is launched from smaller radii ($r\sim 10^{4}\;R_{g}$).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Early Cosmic Merger of Multiple Black Holes

We perform numerical simulations on the merger of multiple black holes (BHs) in primordial gas at early cosmic epochs. We consider two cases of BH mass: $M_{BH} = 30 M_{\odot}$ and $M_{BH} = 10^4 M_{\odot}$. Attention is concentrated on the effect of the dynamical friction by gas in a host object. The simulations incorporate such general relativistic effects as the pericentre shift and gravitational wave emission. As a result, we find that multiple BHs are able to merge into one BH within 100 Myr in a wide range of BH density. The merger mechanism is revealed to be categorized into three types: gas-drag-driven merger (type A), interplay-driven merger (type B), and three-body-driven merger (type C). We find the relation between the merger mechanism and the ratio of the gas mass within the initial BH orbit ($M_{gas}$) to the total BH mass ($ΣM_{BH}$). Type A merger occurs if $M_{gas} \gtrsim 10^5 ΣM_{BH}$, type B if $M_{gas} \lesssim 10^5 ΣM_{BH}$, and type C if $M_{gas} \ll 10^5 ΣM_{BH}$. Supposing the gas and BH density based on the recent numerical simulations on first stars, all the BH remnants from first stars are likely to merge into one BH through the type B or C mechanism. Also, we find that multiple massive BHs ($M_{BH} = 10^4 M_{\odot}$) distributed over several parsec can merge into one BH through the type B mechanism, if the gas density is higher than $5\times 10^6$ cm$^{-3}$. The present results imply that the BH merger may contribute significantly to the formation of supermassive BHs at high redshift epochs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Statistical Properties of Diffuse Lyman-alpha Halos around Star-forming Galaxies at z~2

We present statistical properties of diffuse Lyman-alpha halos (LAHs) around high-$z$ star-forming galaxies with large Subaru samples of Lyman-alpha emitters (LAEs) at $z=2.2$. We make subsamples defined by the physical quantities of LAEs' central Lyman-alpha luminosities, UV magnitudes, Lyman-alpha equivalent widths, and UV slopes, and investigate LAHs' radial surface brightness (SB) profiles and scale lengths $r_n$ as a function of these physical quantities. We find that there exist prominent LAHs around LAEs with faint Lyman-alpha luminosities, bright UV luminosities, and small Lyman-alpha equivalent widths in cumulative radial Lyman-alpha SB profiles. We confirm this trend with the anti-correlation between $r_n$ and Lyman-alpha luminosities (equivalent widths) based on the Spearman's rank correlation coefficient that is $ρ=-0.9$ ($-0.7$) corresponding to the $96\%$ ($93\%$) confidence level, although the correlation between $r_n$ and UV magnitudes is not clearly found in the rank correlation coefficient. Our results suggest that LAEs with properties similar to typical Lyman-break galaxies (with faint Lyman-alpha luminosities and small equivalent widths) possess more prominent LAHs. We investigate scenarios for the major physical origins of LAHs with our results, and find that the cold stream scenario is not preferred, due to the relatively small equivalent widths up to $77$Å in LAHs that include LAEs' central components. There remain two possible scenarios of Lyman-alpha scattering in circum-galactic medium and satellite galaxies that cannot be tested with our observational data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Diffuse Lyman-alpha Halos around Galaxies at z=2.2-6.6: Implications for Galaxy Formation and Cosmic Reionization

We present diffuse Lyman-alpha halos (LAHs) identified in the composite Subaru narrowband images of 100-3600 Lyman-alpha emitters (LAEs) at z=2.2, 3.1, 3.7, 5.7, and 6.6. First, we carefully examine potential artifacts mimicking LAHs that include a large-scale point-spread function (PSF) made by instrumental and atmospheric effects. Based on our critical test with composite images of non-LAE samples whose narrowband-magnitude and source-size distributions are the same as our LAE samples, we confirm that no artifacts can produce a diffuse extended feature similar to our LAHs. After this test, we measure the scale lengths of exponential profile for the LAHs estimated from our z=2.2-6.6 LAE samples of L(Lyman-alpha) > 2 x 10^42 erg s^-1. We obtain the scale lengths of ~ 5-10 kpc at z=2.2-5.7, and find no evolution of scale lengths in this redshift range beyond our measurement uncertainties. Combining this result and the previously-known UV-continuum size evolution, we infer that the ratio of LAH to UV-continuum sizes is nearly constant at z=2.2-5.7. The scale length of our z=6.6 LAH is larger than 5-10 kpc just beyond the error bar, which is a hint that the scale lengths of LAHs would increase from z=5.7 to 6.6. If this increase is confirmed by future large surveys with significant improvements of statistical and systematical errors, this scale length change at z > 6 would be a signature of increasing fraction of neutral hydrogen scattering Lyman-alpha photons, due to cosmic reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Merger Criteria of Multiple Massive Black Holes and the Impact on the Host Galaxy

We perform N-body simulations on a multiple massive black hole (MBH) system in a host galaxy to derive the criteria for successive MBH merger. The calculations incorporate the dynamical friction by stars and general relativistic effects as pericentre shift and gravitational wave recoil. The orbits of MBHs are pursed down to ten Schwarzschild radii (~ 1 AU). As a result, it is shown that about a half of MBHs merge during 1 Gyr in a galaxy with mass $10^{11}M_{\odot}$ and stellar velocity dispersion 240 km/s, even if the recoil velocity is two times as high as the stellar velocity dispersion. The dynamical friction allows a binary MBH to interact frequently with other MBHs, and then the decay of the binary orbits leads to the merger through gravitational wave radiation, as shown by Tanikawa & Umemura (2011). We derive the MBH merger criteria for the masses, sizes, and luminosities of host galaxies. It is found that the successive MBH mergers are expected in bright galaxies, depending on redshifts. Furthermore, we find that the central stellar density is reduced by the sling-shot mechanism and that high-velocity stars with ~ 1000 km/s are generated intermittently in extremely radial orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the evolution of gas clouds exposed to AGN radiation. I. Three-dimensional radiation hydrodynamic simulations

We perform three-dimensional radiation hydrodynamic simulations of uniform dusty gas clouds irradiated by an active galactic nucleus (AGN) to investigate the dependence of evolution of clouds on the ionization parameter $\mathcal{U}$ and the Str{ö}mgren number $\mathcal{N}_{S}$. We find that the evolution can be classified into two cases depending on $\mathcal{U}$. In low $\mathcal{U}$ cases ($\mathcal{U}\approx 10^{-2}$), the evolution is mainly driven by photo-evaporation. A approximately spherically-symmetric evaporation flow with velocity of $100\operatorname{-}150\;\mathrm{km\;s^{-1}}$ is launched from the irradiated face. The cloud is compressed by a D-type shock with losing its mass due to photo-evaporation and is finally turned into a dense filament by $t\lesssim 1.5t_{\mathrm{sc}}$. In high $\mathcal{U}$ cases ($\mathcal{U}\approx 5\times 10^{-2}$), radiation pressure suppresses photo-evaporation from the central part of the irradiated face, reducing photo-evaporation rate. A evaporation flow from the outskirts of the irradiated face is turned into a high velocity ($\lesssim 500\;\mathrm{km\;s^{-1}}$) gas wind because of radiation pressure on dust. The cloud is swept by a radiation pressure-driven shock and becomes a dense gas disk by $t\approx t_{\mathrm{sweep}}$. Star formation is expected in these dense regions for both cases of $\mathcal{U}$. We discuss the influences of the AGN radiation on the clumpy torus. A simple estimate suggests that the clumps are destroyed in timescales shorter than their orbital periods. For the clumpy structure to be maintained over long period, the incident radiation field needs to be sufficiently weaken for most of the clumps, or, some mechanism that creates the clumps continuously is needed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

What is the Physical Origin of Strong Lya Emission? I. Demographics of Lya Emitter Structures

We present the results of structure analyses for a large sample of 426 Lya emitters (LAEs) at z~2.2 that are observed with HST/ACS and WFC3-IR by deep extra-galactic legacy surveys. We confirm that the merger fraction and the average ellipticity of LAE's stellar component are 10-30 % and 0.4-0.6, respectively, that are comparable with previous study results. We successfully identify that some LAEs have a spatial offset between Lya and stellar-continuum emission peaks, d_Lya, by ~2.5-4 kpc beyond our statistical errors. To uncover the physical origin of strong Lya emission found in LAEs, we investigate Lya equivalent width (EW) dependences of these three structural parameters, merger fraction, d_Lya, and ellipticity of stellar distribution in the range of EW(Lya)=20-250A. Contrary to expectations, we find that merger fraction does not significantly increase with Lya EW. We reveal an anti-correlation between d_Lya and EW(Lya) by Kolmogorov-Smirnov (KS) test. There is a trend that the LAEs with a large Lya EW have a small ellipticity. This is consistent with the recent theoretical claims that Lya photons can more easily escape from face-on disks having a small ellipticity, due to less inter-stellar gas along the line of sight, although our KS test indicates that this trend is not statistically significant. Our results of Lya-EW dependence generally support the idea that an HI column density is a key quantity determining Lya emissivity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

What is the physical origin of strong Lya emission? II. Gas Kinematics and Distribution of Lya Emitters

We present a statistical study of velocities of Lya, interstellar (IS) absorption, and nebular lines and gas covering fraction for Lya emitters (LAEs) at z~2. We make a sample of 22 LAEs with a large Lya equivalent width (EW) of > 50A based on our deep Keck/LRIS observations, in conjunction with spectroscopic data from the Subaru/FMOS program and the literature. We estimate the average velocity offset of Lya from a systemic redshift determined with nebular lines to be dv_Lya=234+-9 km s-1. Using a Kolmogorv-Smirnov test, we confirm the previous claim of Hashimoto et al. (2013) that the average dv_Lya of LAEs is smaller than that of LBGs. Our LRIS data successfully identify blue-shifted multiple IS absorption lines in the UV continua of four LAEs on an individual basis. The average velocity offset of IS absorption lines from a systemic redshift is dv_IS=204+-27 km s-1, indicating LAE's gas outflow with a velocity comparable to typical LBGs. Thus, the ratio, R^Lya_ IS = dv_Lya/dv_IS of LAEs, is around unity, suggestive of low impacts on Lya transmission by resonant scattering of neutral hydrogen in the IS medium. We find an anti-correlation between Lya EW and the covering fraction, f_c, estimated from the depth of absorption lines, where f_c is an indicator of average neutral hydrogen column density, N_HI. The results of our study support the idea that N_HI is a key quantity determining Lya emissivity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

First Systematic Search for Oxygen-Line Blobs at High Redshift: Uncovering AGN Feedback and Star-Formation Quenching

We present the first systematic search for extended metal-line [OII]λλ3726,3729 nebulae, or [OII] blobs (OIIBs), at z=1.2 using deep narrowband imaging with a survey volume of 1.9x10^5 Mpc^3 on the 0.62 deg^2 sky of Subaru-XMM Deep Survey (SXDS) field. We discover a giant OIIB, dubbed 'OIIB 1', with a spatial extent over ~75 kpc at a spectroscopic redshift of z=1.18, and also identify a total of twelve OIIBs with a size of >30 kpc. Our optical spectrum of OIIB 1 presents [NeV]λ3426 line at the 6σ level, indicating that this object harbors an obscured type-2 AGN. The presence of gas outflows in this object is suggested by two marginal detections of FeIIλ2587 absorption and FeII*λ2613 emission lines both of which are blueshifted at as large as 500-600 km/s, indicating that the heating source of OIIB 1 is AGN or associated shock excitation rather than supernovae produced by starbursts. The number density of OIIB 1-type giant blobs is estimated to be ~5x10^{-6} Mpc^{-3} at z~1.2, which is comparable with that of AGNs driving outflow at a similar redshift, suggesting that giant OIIBs are produced only by AGN activity. On the other hand, the number density of small OIIBs, 6x10^{-5} Mpc^{-3}, compared to that of z~1 galaxies in the blue cloud in the same M_B range, may imply that 3% of star-forming galaxies at z~1 are quenching star formation through outflows involving extended [OII] emission.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Direct Integration of the Collisionless Boltzmann Equation in Six-dimensional Phase Space: Self-gravitating Systems

We present a scheme for numerical simulations of collisionless self-gravitating systems which directly integrates the Vlasov--Poisson equations in six-dimensional phase space. By the results from a suite of large-scale numerical simulations, we demonstrate that the present scheme can simulate collisionless self-gravitating systems properly. The integration scheme is based on the positive flux conservation method recently developed in plasma physics. We test the accuracy of our code by performing several test calculations including the stability of King spheres, the gravitational instability and the Landau damping. We show that the mass and the energy are accurately conserved for all the test cases we study. The results are in good agreement with linear theory predictions and/or analytic solutions. The distribution function keeps the property of positivity and remains non-oscillatory. The largest simulations are run on 64^6 grids. The computation speed scales well with the number of processors, and thus our code performs efficiently on massively parallel supercomputers.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Formation and Radiative Feedback of First Objects and First Galaxies

First, the formation of first objects driven by dark matter is revisited by high-resolution hydrodynamic simulations. It is revealed that dark matter haloes of ~10^4M_sun can produce first luminous objects with the aid of dark matter cusps. Therefore, the mass of first objects is smaller by roughly two orders of magnitude than in the previous prediction. This implies that the number of Pop III stars formed in the early universe could be significantly larger than hitherto thought. Secondly, the feedback by photo-ionization and photo-dissociation photons in the first objects is explored with radiation hydrodynamic simulations, and it is demonstrated that multiple stars can form in a 10^5M_sun halo. Thirdly, the fragmentation of an accretion disk around a primordial protostar is explored with photo-dissociation feedback. As a result, it is found that the photo-dissociation can reduce the mass accretion rate onto protostars. Also, protostars as small as 0.8M_sun may be ejected and evolve with keeping their mass, which might be detected as "real first stars" in the Galactic halo. Finally, state-of-the-art radiation hydrodynamic simulations are performed to investigate the internal ionization of first galaxies and the escape of ionizing photons. We find that UV feedback by forming massive stars enhances the escape fraction even in a halo as massive as > 6* 10^9M_sun, while it reduces the star formation rate significantly. This may have a momentous impact on the cosmic reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

ARGOT: Accelerated radiative transfer on grids using oct-tree

We present two types of numerical prescriptions that accelerate the radiative transfer calculation around point sources within a three-dimensional Cartesian grid by using the oct-tree structure for the distribution of radiation sources. In one prescription, distant radiation sources are grouped as a bright extended source when the group's angular size, $θ_{\rm s}$, is smaller than a critical value, $θ_{\rm crit}$, and radiative transfer is solved on supermeshes whose angular sizes are similar to that of the group of sources. The supermesh structure is constructed by coarse-graining the mesh structure. With this method, the computational time scales with $N_{\rm m} \log(N_{\rm m}) \log(N_{\rm s})$ where $N_{\rm m}$ and $N_{\rm s}$ are the number of meshes and that of radiation sources, respectively. While this method is very efficient, it inevitably overestimates the optical depth when a group of sources acts as an extended powerful radiation source and affects distant meshes. In the other prescription, a distant group of sources is treated as a bright point source ignoring the spatial extent of the group and the radiative transfer is solved on the meshes rather than the supermeshes. This prescription is simply a grid-based version of {\scriptsize START} by Hasegawa & Umemura and yields better results in general with slightly more computational cost ($\propto N_{\rm m}^{4/3} \log(N_{\rm s})$) than the supermesh prescription. Our methods can easily be implemented to any grid-based hydrodynamic codes and are well-suited to the adaptive mesh refinement methods.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Sub-millimeter brightness of early star-forming galaxies

Based on a three-dimensional model of an early star-forming galaxy, we explore the evolution of the sub-millimeter brightness. The model galaxy is employed from an ultra-high-resolution chemodynamic simulation of a primordial galaxy by Mori & Umemura, where the SFR is ~ 10 Msun/yr at t<0.3 Gyr and several Msun/yr at t>0.3 Gyr. The former phase well reproduces the observed properties of LAEs and the latter does LBGs. We solve the three-dimensional radiative transfer in the clumpy interstellar media in this model galaxy, taking the size distributions of dust grains into account, and calculate the dust temperature as a function of galactic evolutionary time. We find that the clumpiness of interstellar media plays an important role for the sub-millimeter brightness. In the LAE phase, dust grains are concentrated on clumpy star-forming regions that are distributed all over the galaxy, and the grains can effectively absorb UV radiation from stars. As a result, the dust is heated up to T>35 K. In the LBG phase, the continuous supernovae drive dust grains far away from star-forming regions. Then, the grains cannot absorb much radiation from stars, and becomes into a cold state close to the CMB temperature. Consequently, the dust temperature decreases with the evolutionary time, where the mass-weighted mean temperature is T=26 K at t=0.1 Gyr and T=21 K at t=1.0 Gyr. By this analysis, it turns out that the sub-millimeter brightness is higher in the LAE phase than that in the LBG phase, although the dust-to-gas ratio increases monotonically as a function of time. We derive the spectral energy distributions by placing the model galaxy at a given redshift. The peak flux at 850 micron is found to be S_850 ~ 0.2 - 0.9 mJy if the model galaxy is placed at 6>z>2. This means that ALMA can detect an early star-forming galaxy with SFR of ~ 10 Msun/yr by less than one hour integration with 16 antennas.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Successive Merger of Multiple Massive Black Holes in a Primordial Galaxy

Using highly-accurate $N$-body simulations, we explore the evolution of multiple massive black holes (hereafter, MBHs) in a primordial galaxy that is composed of stars and MBHs. The evolution is pursed with a fourth-order Hermite scheme, where not only three-body interaction of MBHs but also dynamical friction by stars are incorporated. Initially, ten MBHs with equal mass of $10^7M_\odot$ are set in a host galaxy with $10^{11}M_\odot$. It is found that 4 - 6 MBHs merge successively within 1 Gyr, emitting gravitational wave radiation. The key process for the successive merger of MBHs is the dynamical friction by field stars, which enhances three-body interactions of MBHs when they enter the central regions of the galaxy. The heaviest MBH always composes a close binary at the galactic center, which shrinks owing to the angular momentum transfer by the third MBH and eventually merges. The angular momentum transfer by the third MBH is due to the sling-shot mechanism. We find that the secular Kozai mechanism does not work for a binary to merge if we include the relativistic pericenter shift. The simulations show that a multiple MBH system can produce a heavier MBH at the galactic center purely through $N$-body process. This merger path can be of great significance for the growth of MBHs in a primordial galaxy. The merger of multiple MBHs may be a potential source of gravitational waves for the Laser Interferometer Space Antenna (LISA) and pulsar timing.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

START: Smoothed particle hydrodynamics with tree-based accelerated radiative transfer

We present a novel radiation hydrodynamics code, START, which is a smoothed particle hydrodynamics (SPH) scheme coupled with accelerated radiative transfer. The basic idea for the acceleration of radiative transfer is parallel to the tree algorithm that is hitherto used to speed up the gravitational force calculation in an N-body system. It is demonstrated that the radiative transfer calculations can be dramatically accelerated, where the computational time is scaled as Np log Ns for Np SPH particles and Ns radiation sources. Such acceleration allows us to readily include not only numerous sources but also scattering photons, even if the total number of radiation sources is comparable to that of SPH particles. Here, a test simulation is presented for a multiple source problem, where the results with START are compared to those with a radiation SPH code without tree-based acceleration. We find that the results agree well with each other if we set the tolerance parameter as < 1.0, and then it demonstrates that START can solve radiative transfer faster without reducing the accuracy. One of important applications with START is to solve the transfer of diffuse ionizing photons, where each SPH particle is regarded as an emitter. To illustrate the competence of START, we simulate the shadowing effect by dense clumps around an ionizing source. As a result, it is found that the erosion of shadows by diffuse recombination photons can be solved. Such an effect is of great significance to reveal the cosmic reionization process.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Two types of Lyman-alpha emitters envisaged from hierarchical galaxy formation

From multi-wavelength observations of LAEs,we know that while many LAEs appear to be young and less massive,a noticeable fraction of LAEs possess much older populations of stars and larger stellar mass.How these two classes of LAEs are concordant with the hierarchical galaxy formation scenario has not been understood clearly so far.In this paper,we model LAEs by three-dimensional cosmological simulations of dark halo merger in a CDM universe.As a result,it is shown that the age of simulated LAEs can spread over a wide range from 2*10^6yr to 9*10^8yr.Also,we find that there are two types of LAEs, in one of which the young half-mass age is comparable to the mean age of stellar component,and in the other of which the young half-mass age is appreciably shorter than the mean age.We define the former as Type 1 LAEs and the latter as Type 2 LAEs.A Type 1 corresponds to early starburst in a young galaxy,whereas a Type 2 does to delayed starburst in an evolved galaxy,as a consequence of delayed accretion of a subhalo onto a larger parent halo.Thus,the same halo can experience a Type 2 LAE-phase as well as a Type 1 LAE-phase in the merger history.Type 1s are expected to be younger than 1.5*10^8yr,less dusty,and less massive with stellar mass M*<10^8 Msun,while Type 2s are older than 1.5*10^8yr,even dustier,and as massive as M*~10^8-10^10Msun.The fraction of Type 2s in all LAEs is a function of redshift.Type 2s discriminated clearly from Type 1s in two color diagram of z'-H vs J-K.We find that the brightness distribution of LyA in Type 2s is more extended than the main stellar component,in contrast to Type 1s.This is not only because delayed starbursts tend to occur in the outskirts of a parent galaxy,but also because LyA photons are effectively absorbed by dust in an evolved galaxy.Hence,the extent of LyA emission may be an additional measure to distinguish Type 2s from Type 1s

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Cosmological Radiative Transfer Comparison Project II: The Radiation-Hydrodynamic Tests

The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer self-consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely allowing even for approximate analytical solutions against which numerical codes can be tested. An alternative validation procedure is to compare different methods against each other on common problems, in order to assess the robustness of the results and establish a range of validity for the methods. Previously, we presented such a comparison for a set of pure radiative transfer tests (i.e. for fixed, non-evolving density fields). This is the second paper of the Cosmological Radiative Transfer (RT) Comparison Project, in which we compare 9 independent RT codes directly coupled to gasdynamics on 3 relatively simple astrophysical hydrodynamics problems: (5) the expansion of an H II region in a uniform medium; (6) an ionization front (I-front) in a 1/r^2 density profile with a flat core, and (7), the photoevaporation of a uniform dense clump. Results show a broad agreement between the different methods and no big failures, indicating that the participating codes have reached a certain level of maturity and reliability. However, many details still do differ, and virtually every code has showed some shortcomings and has disagreed, in one respect or another, with the majority of the results. This underscores the fact that no method is universal and all require careful testing of the particular features which are most relevant to the specific problem at hand.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Formation of globular clusters induced by external ultraviolet radiation

We present a novel scenario for globular cluster (GC) formation, where the ultraviolet (UV) background radiation effectively works so as to produce compact star clusters. Here, we explore the formation of GCs in UV radiation fields. For this purpose, we calculate baryon and dark matter (DM) dynamics in spherical symmetry, incorporating the self-shielding effects by solving the radiative transfer of UV radiation. In addition, we prescribe the star formation in cooled gas components and pursue the dynamics of formed stars. As a result, we find that the evolution of subgalactic objects in UV background radiation are separated into three types, that is, (1) prompt star formation, where less massive clouds ~10^{5-8} M_sun are promptly self-shielded and undergo star formation, (2) delayed star formation, where photoionized massive clouds >10^8 M_sun collapse despite high thermal pressure and are eventually self-shielded to form stars in a delayed fashion, and (3) supersonic infall, where photoionized less massive clouds ~10^{5-8} M_sun contract with supersonic infall velocity and are self-shielded when a compact core forms. In particular, the type (3) is a novel type found in the present simulations, and eventually produces a very compact star cluster. The resultant mass-to-light ratios, half-mass radii, and velocity dispersions for the three types are compared to the observations of GCs, dwarf spheroidals (dSphs), and ultra-compact dwarfs (UCDs). It turns out that the properties of star clusters resulting from supersonic infall match well with those of observed GCs, whereas the other two types are distinct from GCs. Hence, we conclude that supersonic infall in a UV background is a promising mechanism to form GCs.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The escape of ionizing photons from supernova-dominated primordial galaxies

In order to assess the contribution of Lyman break galaxies (LBGs) and Lyman alpha emitters (LAEs) at redshifts 3<z<7 to the ionization of intergalactic medium (IGM), we investigate the escape fractions of ionizing photons from supernova-dominated primordial galaxies by solving the three-dimensional radiative transfer. The model galaxy is employed from an ultra-high-resolution chemodynamic simulation of a primordial galaxy by Mori & Umemura (2006),which well reproduces the observed properties of LAEs and LBGs. The total mass of model galaxy is 10^{11}M_sun. We solve not only photo-ionization but also collisional ionization by shocks. In addition, according to the chemical enrichment, we incorporate the effect of dust extinction, taking the size distributions of dust into account. As a result, we find that dust extinction reduces the escape fractions by a factor 1.5-8.5 in the LAE phase and by a factor 2.5-11 in the LBG phase, while the collisional ionization by shocks increases the escape fractions by a factor ~ 2. The resultant escape fractions are 0.07-0.47 in the LAE phase and 0.06-0.17 in the LBG phase. These results are well concordant with the recent estimations derived from the flux ratio at 1500 angstrom to 900 angstrom of LAEs and LBGs. Combining the resultant escape fractions with the luminosity functions of LAEs and LBGs, we find that high-z LAEs and LBGs can ionize the IGM at z=3-5. However, ionizing radiation from LAEs as well as LBGs falls short to ionize the IGM at z>6. That implies that additional ionization sources may required at z>6.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Quasar Mesolensing - Direct Probe to Substructures around Galaxies -

Recently, ``CDM crisis'' is under discussion. The main point of this crisis is that number of substructures presented by cosmological N-body simulations based on CDM scenario for structure formation is much larger than observed substructures. Therefore, it is crucial for this crisis to discriminate whether expected number of CDM substructures really exist but non-luminous or do not exist. In this paper, we present a new idea to detect such invisible substructures by utilizing a gravitational lensing. Here, we consider quasars that are gravitationally lensed by a foreground galaxy. A substructure around the lensing galaxy may superposed on one of the lensed images of such quasars. In this situation, additional image splitting should occur in the image behind the substructure, and further multiple images are created. This is ``quasar mesolensing''. We estimate separation and time delay between further multiple images due to quasar mesolensing. The expected value is $1 \sim 30$ milli-arcsecond for the separation and future fine resolution imaging enable us to find invisible substructures, and is $1 \sim 10^3$ second for the time delay and high-speed monitoring of such quasar will be able to find ``echo''-like variation due to quasar mesolensing in intrinsic variability of the quasar. Furthermore, we evaluate that the optical depth for the quasar mesolensing is $\sim 0.1$. Consequently, if we monitor a few multiple quasars, we can find ``echo''-like variation in one of the images after intrinsic flux variations of quasars.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

The Stellar IMF in Very Metal-Deficient Gas

In the context of the star formation through the fragmentation of an extremely metal-deficient protogalactic cloud, the gravitational collapse of filamentary gas clouds is explored with H$_2$ and HD chemistry. It is found by 1D hydrodynamical simulations that the cloud evolution is prescribed mainly by the initial density ($n_0$) and H$_2$ abundance ($x_{\rm H_2,0}$). In particular, it turns out that the evolution of low-density filaments ($n_0 \lesssim 10^5$ cm$^{-3}$) bifurcates at a critical H$_2$ abundance of $x_{\rm H_2,cr}\simeq 3\times 10^{-3}$, beyond which HD cooling overwhelms H$_2$ cooling. The numerical results indicate that the stellar IMF is likely to be double-peaked and deficient in sub-solar mass stars, where the high mass peak of the IMF is around $10M_\odot$ or $10^2M_\odot$, dependently on the initial density and H$_2$ abundance. If the gas in protogalactic clouds is photoionized by UV radiation or shock-heated, the H$_2$ abundance could exceed $x_{\rm H_2,cr}\simeq 3\times 10^{-3}$ by H$^-$ reactions. Then, the high mass peak would be $O(10) M_\odot$.

0 citations0 reviewsNo code linkedOpen access