著者
Miki Yohei Umemura Masayuki
出版者
Oxford University Press
雑誌
Monthly notices of the Royal Astronomical Society (ISSN:00358711)
巻号頁・発行日
vol.475, no.2, pp.2269-2281, 2018-04
被引用文献数
20

Providing initial conditions is an essential procedure for numerical simulations of galaxies. The initial conditions for idealized 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 = 2097 152 and 8388 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 https://bitbucket.org/ymiki/magi.
著者
Takahashi Rohta Umemura Masayuki
出版者
Oxford University Press
雑誌
Monthly Notices of the Royal Astronomical Society (ISSN:00358711)
巻号頁・発行日
vol.464, no.4, pp.4567-4585, 2017-02
被引用文献数
12

We present a general relativistic radiative transfer code, ARTIST (Authentic Radiative Transfer In Space–Time), that is a perfectly causal scheme to pursue the propagation of radiation with absorption and scattering around a Kerr black hole. The code explicitly solves the invariant radiation intensity along null geodesics in the Kerr–Schild coordinates, and therefore properly includes light bending, Doppler boosting, frame dragging, and gravitational redshifts. The notable aspect of ARTIST is that it conserves the radiative energy with high accuracy, and is not subject to the numerical diffusion, since the transfer is solved on long characteristics along null geodesics. We first solve the wavefront propagation around a Kerr black hole that was originally explored by Hanni. This demonstrates repeated wavefront collisions, light bending, and causal propagation of radiation with the speed of light. We show that the decay rate of the total energy of wavefronts near a black hole is determined solely by the black hole spin in late phases, in agreement with analytic expectations. As a result, the ARTIST turns out to correctly solve the general relativistic radiation fields until late phases as t ∼ 90 M. We also explore the effects of absorption and scattering, and apply this code for a photon wall problem and an orbiting hotspot problem. All the simulations in this study are performed in the equatorial plane around a Kerr black hole. The ARTIST is the first step to realize the general relativistic radiation hydrodynamics.
著者
Umemura Masayuki Susa Hajime Hasegawa Kenji Suwa Tamon Semelin Benoit
出版者
Oxford University Press
雑誌
Progress of theoretical and experimental physics (ISSN:20503911)
巻号頁・発行日
vol.1, pp.01A306, 2012-10
被引用文献数
8

Firstly, the formation of first objects driven by dark matter is revisited by high-resolution hydrodynamic simulations. It is revealed that dark matter halos of ∼104 M⊙ 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 Population III stars formed in the early universe could be significantly larger than hitherto thought. Secondly, 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 105 M⊙ 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 photo-dissociation can reduce the mass accretion rate onto protostars. Also, protostars as small as 0.8 M⊙ may be ejected and evolve while keeping their mass, which might be detected as the “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 × 109 M⊙, while it reduces the star formation rate significantly. This may have a momentous impact on cosmic reionization.