Gaia Data Release 2 revealed detailed structures of nearby stars in phase space. These include the Hercules stream, whose origin is still debated. Earlier numerical studies conjectured that the observed structures originate from orbits in resonance with the bar, based on static potential models for the Milky Way. We, in contrast, approach the problem via a self-consistent, dynamic, and morphologically well-resolved model, namely a full $N$-body simulation of the Milky Way. Our simulation comprises about 5.1 billion particles in the galactic stellar bulge, bar, disk, and dark-matter halo and is evolved to 10 Gyr. Our model's disk component is composed of 200 million particles, and its simulation snapshots are stored every 10 Myr, enabling us to resolve and classify resonant orbits of representative samples of stars. After identifying the Sun's position in the simulation, we compare the distribution of stars in its neighborhood with Gaia's astrometric data, thereby establishing the role of identified resonantly trapped stars in the formation of Hercules-like structures. From our orbital spectral-analysis we identify multiple resonances and conclude that the Hercules stream is dominated by the 4:1 and 5:1 outer Lindblad and corotation resonances. In total, this yields a trimodal structure of the Hercules stream. From the relation between resonances and ridges in phase space, we conclude that the pattern speed of the Milky-Way bar is 40$-$45 km s$^{-1}$ kpc$^{-1}$.

Most stars are formed as star clusters in galaxies, which then disperse into galactic disks. Upcoming exascale supercomputational facilities will enable performing simulations of galaxies and their formation by resolving individual stars (star-by-star simulations). This will substantially advance our understanding of star formation in galaxies, star cluster formation, and assembly histories of galaxies. In previous galaxy simulations, a simple stellar population approximation was used. It is, however, difficult to improve the mass resolution with this approximation. Therefore, a model for forming individual stars that can be used in simulations of galaxies must be established. In this first paper of a series of the SIRIUS (SImulations Resolving IndividUal Stars) project, we propose and demonstrate a new stochastic star formation model for star-by-star simulations. In this model, an assumed stellar initial mass function (IMF) is randomly assigned to newly formed stars. We then introduce a maximum search radius to assemble the mass from surrounding gas particles to form star particles. In this study, we perform a series of $N$-body/smoothed particle hydrodynamics simulations of star cluster formations from turbulent molecular clouds as test cases. The IMF can be correctly sampled if a maximum search radius that is larger than the value estimated from the threshold density for star formation is adopted. In small clouds, the formation of massive stars is highly stochastic because of the small number of stars. We confirm that the star formation efficiency and threshold density do not strongly affect the results. Herein, we demonstrate that our models can be applied to simulations varying from star clusters to galaxies for a wide range of resolutions.

Several binary black holes (BBHs) have been observed using gravitational wave detectors. For the formation mechanism of BBHs, two main mechanisms, isolated binary evolution and dynamical formation in dense star clusters, have been suggested. Future observations are expected to provide more information about BBH distributions, and it will help us to distinguish the two formation mechanisms. For the star cluster channel, globular clusters have mainly been investigated. However, recent simulations have suggested that BBH formation in open clusters is not negligible. We estimate a local merger rate density of BBHs originated from open clusters using the results of our $N$-body simulations of open clusters with four different metallicities. We find that the merger rate per cluster is the highest for our 0.1 solar metallicity model. Assuming a cosmic star formation history and a metallicity evolution with dispersion, we estimate the local merger rate density of BBHs originated from open clusters to be $\sim 70~{\rm yr}^{-1} {\rm Gpc}^{-3}$. This value is comparable to the merger rate density expected from the first and second observation runs of LIGO and Virgo. In addition, we find that BBH mergers obtained from our simulations can reproduce the distribution of primary mass and mass ratio of merging BBHs estimated from the LIGO and Virgo observations.

Black hole-main sequence star (BH-MS) binaries are one of the targets of the future data releases of the astrometric satellite {\it Gaia}. They are supposed to be formed in two main sites: a galactic field and star clusters. However, previous work has never predicted the number of BH-MS binaries originating in the latter site. In this paper, we estimate the number of BH-MS binaries formed in open clusters and detectable with {\it Gaia} based on the results of {\it N}-body simulations. By considering interstellar extinction in the Milky Way (MW) and observational constraints, we predict $\sim 10$ BH-MS binaries are observable. We also find that chemical abundance patterns of companion MSs will help us to identify the origin of the binaries as star clusters. Such MSs are not polluted by outflows of the BH progenitors, such as stellar winds and supernova ejecta. Chemical anomalies might be a good test to confirm the origin of binaries with relatively less massive MSs ($\lesssim 5M_{\odot}$), orbital periods ($\sim 1.5\;$year) and higher eccentricities ($e \gtrsim 0.1$).

We estimate formation rates of LB-1-like systems through dynamical interactions in the framework of the theory of stellar evolution before the discovery of the LB-1 system. The LB-1 system contains $\sim 70M_\odot$ black hole (BH), so-called pair instability (PI)-gap BH, and B-type star with solar metallicity, and has nearly zero eccentricity. The most efficient formation mechanism is as follows. In an open cluster, a naked helium (He) star (with $\sim 20M_\odot$) collides with a heavy main-sequence (MS) star (with $\sim 50M_\odot$) which has a B-type companion. The collision results in a binary consisting of the collision product and B-type star with a high eccentricity. The binary can be circularized through the dynamical tide with radiative damping of the collision-product envelope. Finally, the collision product collapses to a PI-gap BH, avoiding pulsational pair instability and pair instability supernovae because its He core is as massive as the pre-colliding naked He star. We find that the number of LB-1-like systems in the Milky Way galaxy is $\sim 0.01 (\rho_{\rm oc} / 10^4 M_\odot \mbox{pc}^{-3})$, where $\rho_{\rm oc}$ is the initial mass densities of open clusters. If we take into account LB-1-like systems with O-type companion stars, the number increases to $\sim 0.03 (\rho_{\rm oc} / 10^4 M_\odot \mbox{pc}^{-3})$. This mechanism can form LB-1-like systems at least 10 times more efficiently than the other mechanisms: captures of B-type stars by PI-gap BHs, stellar collisions between other type stars, and stellar mergers in hierarchical triple systems. We conclude that no dynamical mechanism can explain the presence of the LB-1 system.

We investigate the role of the Keplerian tidal field generated by a supermassive black hole (SMBH) on the three-body dynamics of stellar mass black holes. We consider two scenarios occurring close to the SMBH: the breakup of unstable triples and three-body encounters between a binary and a single. These two cases correspond to the hard and soft binary cases, respectively. The tidal field affects the breakup of triples by tidally limiting the system, so that the triples break earlier with lower breakup velocity, leaving behind slightly larger binaries (relative to the isolated case). The breakup direction becomes anisotropic and tends to follow the shape of the Hill region of the triple, favouring breakups in the radial direction. Furthermore, the tidal field can torque the system, leading to angular momentum exchanges between the triple and its orbit about the SMBH. This process changes the properties of the final binary, depending on the initial angular momentum of the triple. Finally, the tidal field also affects binary-single encounters: binaries tend to become both harder and more eccentric with respect to encounters that occur in isolation. Consequently, single-binary scattering in a deep Keplerian potential produces binaries with shorter gravitational wave merger timescales.

In clustered environments, stellar encounters can liberate planets from their host stars via close encounters. Although the detection probability of planets suggests that the planet population in open clusters resembles that in the field, only a few dozen planet-hosting stars have been discovered in open clusters. We explore the survival rates of planets against stellar encounters in open clusters similar to the Pleiades, Hyades, and Praesepe and embedded clusters. We performed a series of N-body simulations of high-density and low-density open clusters, open clusters that grow via mergers of subclusters, and embedded clusters. We semi-analytically calculated the survival rate of planets in star clusters up to 1Gyr using relative velocities, masses, and impact parameters of intruding stars. Less than 1.5% of close-in planets within 1 AU and at most 7% of planets with 1-10 AU are ejected by stellar encounters in clustered environments after the dynamical evolution of star clusters. If a planet population from 0.01-100 AU in an open cluster initially follows the probability distribution function of exoplanets with semi-major axis ($a_p$) between 0.03-3 AU in the field discovered by RV surveys, the PDF of surviving planets beyond ~10 AU in open clusters can be slightly modified to $\propto a_p^{-0.76}$. The production rate of free-floating planets (FFPs) per star is 0.0096-0.18, where we have assumed that all the stars initially have one giant planet with a mass of 1--13 MJ in a circular orbit. The expected frequency of FFPs is compatible with the upper limit on that of FFPs indicated by recent microlensing surveys. Our survival rates of planets in open clusters suggest that planets within 10 AU around FGKM-type stars are rich in relatively-young (<~10-100 Myr for open clusters and ~1-10 Myr for embedded clusters), less massive open clusters, which are promising targets for planet searches.

Star cluster complexes such as the Carina Nebula can have formed in turbulent giant molecular clouds. We perform a series of $N$-body simulations starting from subclustering initial conditions based on hydrodynamic simulations of turbulent molecular clouds. These simulations finally result in the formation of star cluster complexes consisting of several subclusters (clumps). We obtain the inter-clump velocity distribution, the size of the region, and the mass of the most massive cluster in our simulated complex and compare the results with observed ones (the Carina Nebula and NGC 2264). The one-dimensional inter-clump velocity dispersion obtained from our simulations is $2.9\pm0.3$ and $1.4\pm0.4$ km s$^{-1}$ for the Carina- and NGC 2264-like models, respectively, which are consistent with those obtained from Gaia Data Release 2: 2.35 and 0.99 km s$^{-1}$ for the Carina Nebula and NGC 2264, respectively. We estimate that the masses of the parental molecular clouds for the Carina Nebula and the NGC 2264 are $4\times 10^5$ and $4\times 10^4M_{\odot}$, respectively.

Recent spectroscopic analysis has set an upper limit to the age of the S-stars, the ~30 B-type stars in highly eccentric orbits around the supermassive black hole (SMBH) in the Galactic center. The inferred age (<15 Myr) is in tension with the binary break-up scenario proposed to explain their origin. However, the new estimate is compatible with the age of the disk of O-type stars that lies at a farther distance from the SMBH. Here we investigate a new formation scenario, assuming that both S-stars and the O-type stars were born in the same disk around SgrA*. We simulate encounters between binaries of the stellar disk and stellar black holes from a dark cusp around SgrA*. We find that B-type binaries can be easily broken up by the encounters and their binary components are kicked into highly eccentric orbits around the SMBH. In contrast, O-type binaries are less frequently disrupted and their members remain in low eccentricity orbits. This mechanism can reproduce 12 S-stars just by assuming that the binaries initially lie within the stellar disk as observed nowadays. To reproduce all the S-stars, the original disk must have been extended down to 0.006 pc. However in this case many B- and O-type stars remain in low eccentricity orbits below 0.03 pc, in contrast with the observations. Therefore, some other mechanism is necessary to disrupt the disk below 0.03 pc. This scenario can also explain the high eccentricity of the G-objects, if they have a stellar origin.

We study the stellar dynamics of the first star clusters after intermediate-mass black holes (IMBHs) are formed via runaway stellar collisions. We use the outputs of cosmological simulations of Sakurai et al. (2017) to follow the star cluster evolution in a live dark matter (DM) halo. Mass segregation within a cluster promotes massive stars to be captured by the central IMBH occasionally, causing tidal disruption events (TDEs). We find that the TDE rate scales with the IMBH mass as $\dot{N}_{\rm TDE}\sim0.3\,{\rm Myr}^{-1}(M_{\rm IMBH}/1000\,{\rm M}_{\odot})^2$. The DM component affects the star cluster evolution by stripping stars from the outer part. When the DM density within the cluster increases, the velocity dispersion of the stars increases, and then the TDE rate decreases. By the TDEs, the central IMBHs grow to as massive as $700-2500\,{\rm M}_{\odot}$ in 15 million years. The IMBHs are possible seeds for the formation of supermassive BHs observed at $z\gtrsim 6-7$, if a large amount of gas is supplied through galaxy mergers and/or large-scale gas accretion, or they might remain as IMBHs from the early epochs to the present-day Universe.

In order to investigate the formation rate of binary black holes (BBHs) in stellar clusters with a mass comparable to open clusters, we performed a series of direct $N$-body simulations of open clusters with a mass of $2.5\times10^3$ (Model A) and $10^4 M_{\odot}$ (Model B). Since such low-mass clusters would have been more populous than globular clusters when they were born, low-mass clusters are also candidates as the origin of BBHs which are the source of the gravitational waves. In model A, most of BBHs merged within 10 Gyr formed via dynamically formed main-sequence binary stars and stable and unstable mass transfer between them since open clusters collapse within the main-sequence life-time of massive stars. These binaries, therefore, have little eccentricities. The fraction of such binaries among all merging BBHs increases as the cluster mass decrease due to the shorter relaxation time. In our simulations, $4.0\times10^{-5}$ and $1.7\times10^{-5}$ BBHs per solar mass merged within 10 Gyr for models A and B, respectively. These values correspond to $\sim$ 20-50% of the number of mergers per solar mass originated from globular clusters with a mass of $10^5$-$10^6M_{\odot}$. Thus, the contribution of BBHs originated from open clusters is not negligible. The estimated mergers rate density in the local universe is about 0.3 yr$^{-1}$ Gpc$^{-3}$ assuming a cluster mass function with a power of $-$2.

We have newly developed a Parallelized Particle-Particle Particle-tree code for Planet formation, PENTACLE, which is a parallelized hybrid $N$-body integrator executed on a CPU-based (super)computer. PENTACLE uses a 4th-order Hermite algorithm to calculate gravitational interactions between particles within a cutoff radius and a Barnes-Hut tree method for gravity from particles beyond. It also implements an open-source library designed for full automatic parallelization of particle simulations, FDPS (Framework for Developing Particle Simulator) to parallelize a Barnes-Hut tree algorithm for a memory-distributed supercomputer. These allow us to handle $1-10$ million particles in a high-resolution $N$-body simulation on CPU clusters for collisional dynamics, including physical collisions in a planetesimal disc. In this paper, we show the performance and the accuracy of PENTACLE in terms of $\tilde{R}_{\rm cut}$ and a time-step $\Delta t$. It turns out that the accuracy of a hybrid $N$-body simulation is controlled through $\Delta t / \tilde{R}_{\rm cut}$ and $\Delta t / \tilde{R}_{\rm cut} \sim 0.1$ is necessary to simulate accurately accretion process of a planet for $\geq 10^6$ years. For all those who interested in large-scale particle simulations, PENTACLE customized for planet formation will be freely available from https://github.com/PENTACLE-Team/PENTACLE under the MIT lisence.

We construct a model for the Milky Way Galaxy composed of a stellar disc and bulge embedded in a dark-matter halo. All components are modelled as $N$-body systems with up to 8 billion equal-mass particles and integrated up to an age of 10\,Gyr. We find that net angular-momentum of the dark-matter halo with a spin parameter of $\lambda=0.06$ is required to form a relatively short bar ($\sim 4$\,kpc) with a high pattern speed (40--50\,km\,s$^{-1}$). By comparing our model with observations of the Milky Way Galaxy, we conclude that a disc mass of $\sim 3.7\times10^{10}M_{\odot}$ and an initial bulge scale length and velocity of $\sim 1$\,kpc and $\sim 300$\,km\,s$^{-1}$, respectively, fit best to the observations. The disc-to-total mass fraction ($f_{\rm d}$) appears to be an important parameter for the evolution of the Galaxy and models with $f_{\rm d}\sim 0.45$ are most similar to the Milky Way Galaxy. In addition, we compare the velocity distribution in the solar neighbourhood in our simulations with observations in the Milky Way Galaxy. In our simulations the observed gap in the velocity distribution, which is expected to be caused by the outer Lindblad resonance (the so-called Hercules stream), appears to be a time-dependent structure. The velocity distribution changes on a time scale of 20--30\,Myr and therefore it is difficult to estimate the pattern speed of the bar from the shape of the local velocity distribution alone.

Recent developments in computer hardware and software enable researchers to simulate the self-gravitating evolution of galaxies at a resolution comparable to the actual number of stars. Here we present the results of a series of such simulations. We performed $N$-body simulations of disk galaxies with between 100 and 500 million particles over a wide range of initial conditions. Our calculations include a live bulge, disk, and dark matter halo, each of which is represented by self-gravitating particles in the $N$-body code. The simulations are performed using the gravitational $N$-body tree-code Bonsai running on the Piz Daint supercomputer. We find that the time scale over which the bar forms increases exponentially with decreasing disk-mass fraction and that the bar formation epoch exceeds a Hubble time when the disk-mass fraction is $\sim0.35$. These results can be explained with the swing-amplification theory. The condition for the formation of $m=2$ spirals is consistent with that for the formation of the bar, which is also an $m=2$ phenomenon. We further argue that the non-barred grand-design spiral galaxies are transitional, and that they evolve to barred galaxies on a dynamical timescale. We also confirm that the disk-mass fraction and shear rate are important parameters for the morphology of disk galaxies. The former affects the number of spiral arms and the bar formation epoch, and the latter determines the pitch angle of the spiral arms.

We perform a set of cosmological simulations of early structure formation with incorporating baryonic streaming motions. We present a case where a significantly elongated gas cloud with $\sim\!10^4\,$solar masses (${\rm M_\odot}$) is formed in a pre-galactic ($\sim\!10^7\,{\rm M_\odot}$) dark halo. The gas streaming into the halo compresses and heats the massive filamentary cloud to a temperature of $\sim\!10,000\,$Kelvin. The gas cloud cools rapidly by atomic hydrogen cooling, and then by molecular hydrogen cooling down to $\sim\!400\,$Kelvin. The rapid decrease of the temperature and hence of the Jeans mass triggers fragmentation of the filament to yield multiple gas clumps with a few hundred solar masses. We estimate the mass of the primordial star formed in each fragment by adopting an analytic model based on a large set of radiation hydrodynamics simulations of protostellar evolution. The resulting stellar masses are in the range of $\sim\!50$-$120\,{\rm M_\odot}$. The massive stars gravitationally attract each other and form a compact star cluster. We follow the dynamics of the star cluster using a hybrid $N$-body simulation. We show that massive star binaries are formed in a few million years through multi-body interactions at the cluster center. The eventual formation of the remnant black holes will leave a massive black hole binary, which can be a progenitor of strong gravitational wave sources similar to those recently detected by the Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO).

Advanced LIGO achieved the first detection of the gravitational wave, which was from a merging binary black hole (BBH). In the near future, more merger events will be observed, and the mass distribution of them will become available. The mass distribution of merger events reflects the evolutionary path of BBHs: dynamical formation in dense star clusters or common envelope evolution from primordial binaries. In this paper, we estimate the detection rate of merging BBHs which dynamically formed in dense star clusters by combining the results of $N$-body simulations, modeling of globular clusters, and cosmic star-cluster formation history. We estimate that the merger rate density in the local universe within the redshift of 0.1 is 13--57 Gpc$^{-3}$ yr$^{-1}$. We find that the detection rate is 0.23--4.6 per year for the current sensitivity limit and that it would increase to 5.1--99 per year for the designed sensitivity which will be achieved in 2019. The distribution of merger rate density in the local universe as a function of redshifted chirp mass has a peak close to the low-mass end.The chirp mass function of the detected mergers, on the other hand, has a peak at the high-mass end, but is almost flat. This difference is simply because the detection range is larger for more massive BBHs.

We study the formation of massive black holes in the first star clusters. We first locate star-forming gas clouds in proto-galactic haloes of $\gtrsim \!10^7\,{\rm M}_{\odot}$ in cosmological hydrodynamics simulations and use them to generate the initial conditions for star clusters with masses of $\sim \!10^5\,{\rm M}_{\odot}$. We then perform a series of direct-tree hybrid $N$-body simulations to follow runaway stellar collisions in the dense star clusters. In all the cluster models except one, runaway collisions occur within a few million years, and the mass of the central, most massive star reaches $\sim \!400-1900\,{\rm M}_{\odot}$. Such very massive stars collapse to leave intermediate-mass black holes (IMBHs). The diversity of the final masses may be attributed to the differences in a few basic properties of the host haloes such as mass, central gas velocity dispersion, and mean gas density of the central core. Finally, we derive the IMBH mass to cluster mass ratios, and compare them with the observed black hole to bulge mass ratios in the present-day Universe.

AG Peg is known as the slowest symbiotic nova, which experienced its nova-like outburst around 1850. After 165 years, during June of 2015, it erupted again showing characteristics of the Z And-type outburst. The primary objective is to determine basic characteristics, the nature and type of the 2015 outburst of AG Peg. We achieved this aim by modelling the spectral energy distribution using low-resolution spectroscopy (330-750 nm), medium-resolution spectroscopy (420-720 nm; R=11000), and $UBVR_{\rm C}I_{\rm C}$ photometry covering the 2015 outburst with a high cadence. Optical observations were complemented with the archival HST and FUSE spectra from the preceding quiescence. During the outburst, the luminosity of the hot component was in the range of 2-11$\times 10^{37}(d/0.8{\rm kpc})^2$ erg/s. To generate the maximum luminosity the white dwarf (WD) had to accrete at $\sim 3\times 10^{-7}$ M$_{\odot}yr^{-1}$, which exceeds the stable-burning limit and thus led to blowing optically thick wind from the WD. We determined its mass-loss rate to a few $\times 10^{-6}$ M$_{\odot}yr^{-1}$. At the high temperature of the ionising source, $1.5-2.3\times 10^5$ K, the wind converted a fraction of the WD's photospheric radiation into the nebular emission that dominated the optical. A one order of magnitude increase of the emission measure, from a few $\times 10^{59}(d/0.8 {\rm kpc})^2$ cm$^{-3}$ during quiescence, to a few $\times 10^{60}(d/0.8\,{\rm kpc})^2$ cm$^{-3}$ during the outburst, caused a 2 mag brightening in the LC, which is classified as the Z And-type of the outburst. The very high nebular emission and the presence of a disk-like HI region encompassing the WD, as indicated by a significant broadening and high flux of the Raman-scattered OVI 6825 \AA\ line during the outburst, is consistent with the ionisation structure of hot components in symbiotic stars during active phases.

The abundance of elements synthesized by the rapid neutron-capture process (r-process elements) of extremely metal-poor (EMP) stars in the Local Group galaxies gives us clues to clarify the early evolutionary history of the Milky Way halo. The Local Group dwarf galaxies would have similarly evolved with building blocks of the Milky Way halo. However, how the chemo-dynamical evolution of the building blocks affects the abundance of r-process elements is not yet clear. In this paper, we perform a series of simulations using dwarf galaxy models with various dynamical times and total mass, which determine star-formation histories. We find that galaxies with dynamical times longer than 100 Myr have star formation rates less than $10^{-3} M_{\odot}$ yr$^{-1}$ and slowly enrich metals in their early phase. These galaxies can explain the observed large scatters of r-process abundance in EMP stars in the Milky Way halo regardless of their total mass. On the other hand, the first neutron star merger appears at a higher metallicity in galaxies with a dynamical time shorter than typical neutron star merger times. The scatters of r-process elements mainly come from inhomogeneity of the metals in the interstellar medium whereas the scatters of $\alpha$-elements are mostly due to the difference in the yield of each supernova. Our results demonstrate that the future observations of r-process elements in EMP stars will be able to constrain the early chemo-dynamical evolution of the Local Group galaxies.

The scaling relations and the star formation laws for molecular cloud complexes in the Milky Way is investigated. We compare their masses $M_{\rm gas}$, mass surface densities $\Sigma_{M_{\rm gas}}$, radii $R$, velocity dispersions $\sigma$, star formation rates $SFR$, and SFR densities $\Sigma_{\rm SFR}$ with those of structures ranging from cores, clumps, Giant Molecular Clouds (GMCs), to Molecular Cloud Complexes (MCCs), and to Galaxies, spanning 8 orders of magnitudes in size and 13 orders of magnitudes in mass. MCC are mostly large ($R>50$ pc), massive ($\sim 10^{6}$\,\msun) gravitationally unbound cloud structures. This results in the following universal relations: $\sigma\sim R^{0.5}$, $M_{\rm gas}\sim R^{2}$, $\Sigma_{\rm SFR}\sim \Sigma_{M_{\rm gas}}^{1.5}$, ${SFR}\sim {M_{\rm gas}}^{0.9}$, and ${SFR}\sim {\sigma}^{2.7}$. Variations in the slopes and the coefficients of these relations are found at individual scales signifying different physics acting at different scales. Additionally, there are breaks at the MCC scale in the $\sigma-R$ relation and between the starburst and the normal star-forming objects in the $SFR-M_{\rm gas}$ and $\Sigma_{\rm SFR}$-$\Sigma_{\rm M_{\rm gas}}$ relations. We propose to use the Schmidt-Kennicutt diagram to distinguish the starburst from the normal star-forming structures by applying a $\Sigma_{M_{\rm gas}}$ threshold of $\sim100$\,\msun pc$^{-2}$ and a $\Sigma_{\rm SFR}$ threshold of 1\,\msun yr$^{-1}$ kpc$^{-2}$. Mini-starburst complexes have enhanced $\Sigma_{\rm SFR}$ ($>$1\,\msun yr$^{-1}$ kpc$^{-2}$), probably caused by dynamic events such as radiation pressure, colliding flows, or spiral arm gravitational instability. Because of the dynamical evolution, gravitational boundedness does not play a significant role in characterizing the star formation activity of MCCs, especially the mini-starburst complexes.

Recent observations have revealed a variety of young star clusters, including embedded systems, young massive clusters, and associations. We study the formation and dynamical evolution of these clusters using a combination of simulations and theoretical models. Our simulations start with a turbulent molecular cloud that collapses under its own gravity. The stars are assumed to form in the densest regions in the collapsing cloud after an initial free-fall times of the molecular cloud. The dynamical evolution of these stellar distributions are continued by means of direct $N$-body simulations. The molecular clouds typical for the Milky Way Galaxy tend to form embedded clusters which evolve to resemble open clusters. The associations were initially considerably more clumpy, but lost their irregularity in about a dynamical time scale due to the relaxation process. The densest molecular clouds, which are absent in the Milky Way but are typical in starburst galaxies, form massive young star clusters. They indeed are rare in the Milky Way. Our models indicate a distinct evolutionary path from molecular clouds to open clusters and associations or to massive star clusters. The mass-radius relation for both types of evolutionary tracks excellently matches the observations. According to our calculations the time evolution of the half-mass radius for open clusters and associations follows $r_{\rm h}/{\rm pc}=2.7(t_{\rm age}/{\rm pc})^{2/3}$, whereas for massive star clusters $r_{\rm h}/{\rm pc}=0.34(t_{\rm age}/{\rm Myr})^{2/3}$. Both trends are consistent with the observed age-mass-radius relation for clusters in the Milky Way.

The rapid neutron-capture process (r-process) is a major process to synthesize elements heavier than iron, but the astrophysical site(s) of r-process is not identified yet. Neutron star mergers (NSMs) are suggested to be a major r-process site from nucleosynthesis studies. Previous chemical evolution studies however require unlikely short merger time of NSMs to reproduce the observed large star-to-star scatters in the abundance ratios of r-process elements relative to iron, [Eu/Fe], of extremely metal-poor stars in the Milky Way (MW) halo. This problem can be solved by considering chemical evolution in dwarf spheroidal galaxies (dSphs) which would be building blocks of the MW and have lower star formation efficiencies than the MW halo. We demonstrate that enrichment of r-process elements in dSphs by NSMs using an N-body/smoothed particle hydrodynamics code. Our high-resolution model reproduces the observed [Eu/Fe] by NSMs with a merger time of 100 Myr when the effect of metal mixing is taken into account. This is because metallicity is not correlated with time up to ~ 300 Myr from the start of the simulation due to low star formation efficiency in dSphs. We also confirm that this model is consistent with observed properties of dSphs such as radial profiles and metallicity distribution. The merger time and the Galactic rate of NSMs are suggested to be <~ 300 Myr and ~ $10^{-4}$ yr$^{-1}$, which are consistent with the values suggested by population synthesis and nucleosynthesis studies. This study supports that NSMs are the major astrophysical site of r-process.

We simulate the formation and evolution of young star clusters using the combination of SPH simulations and direct N-body simulations. We start by performing SPH simulations of the giant molecular cloud with a turbulent velocity field, a mass of $4\times10^4$ to $5\times10^6M_{\odot}$, and a density between $1.7\times10^3$ and $170cm^{-3}$. We continue the SPH simulations for a free-fall time scale, and analyze the resulting structure of the collapsed cloud. We subsequently replace a density-selected subset of SPH particles with stars by adopting a local star-formation efficiency proportional to $\rho^{1/2}$. As a consequence, the local star formation efficiency exceeds 30 %, whereas globally only a few % of the gas is converted to stars. The stellar distribution by the time gas is converted to stars is very clumpy, with typically a dozen bound conglomerates that consist of 100 to $10^4$ stars. We continue to evolve the stars dynamically using the collisional N-body method, which accurately treats all pairwise interactions, stellar collisions and stellar evolution. We analyze the results of the N-body simulations when the stars have an age of 2 Myr and 10 Myr. During the dynamical simulations, massive clusters grow via hierarchical merging of smaller clusters. The shape of the cluster mass function that originates from an individual molecular cloud is consistent with a Schechter function with a power-law slope of -1.73 at 2 Myr and -1.67 at 10 Myr, which fits to observed cluster mass function of the Carina region. The superposition of mass functions have a power-law slope of < -2, which fits the observed mass function of star clusters in the Milky Way, M31 and M83. We further find that the mass of the most massive cluster formed in a single molecular cloud with a mass of $M_g$ scales with $6.1M_g^{0.51}$ which also agrees with recent observation of the GMC and young clusters in M51.

We have simulated, for the first time, the long term evolution of the Milky Way Galaxy using 51 billion particles on the Swiss Piz Daint supercomputer with our $N$-body gravitational tree-code Bonsai. Herein, we describe the scientific motivation and numerical algorithms. The Milky Way model was simulated for 6 billion years, during which the bar structure and spiral arms were fully formed. This improves upon previous simulations by using 1000 times more particles, and provides a wealth of new data that can be directly compared with observations. We also report the scalability on both the Swiss Piz Daint and the US ORNL Titan. On Piz Daint the parallel efficiency of Bonsai was above 95%. The highest performance was achieved with a 242 billion particle Milky Way model using 18600 GPUs on Titan, thereby reaching a sustained GPU and application performance of 33.49 Pflops and 24.77 Pflops respectively.

Classical Cepheids are useful tracers of the Galactic young stellar population because their distances and ages can be determined from their period-luminosity and period-age relations. In addition, the radial velocities and chemical abundance of the Cepheids can be derived from spectroscopic observations, providing further insights into the structure and evolution of the Galaxy. Here, we report the radial velocities of classical Cepheids near the Galactic Center, three of which were reported in 2011, the other reported for the first time. The velocities of these Cepheids suggest that the stars orbit within the Nuclear Stellar Disk, a group of stars and interstellar matter occupying a region of 200 pc around the Center, although the three-dimensional velocities cannot be determined until the proper motions are known. According to our simulation, these four Cepheids formed within the Nuclear Stellar Disk like younger stars and stellar clusters therein.

Young massive clusters are as young as open clusters but more massive and compact compared with typical open clusters. The formation process of young massive clusters is still unclear, and it is an open question whether the formation process is the same as typical open clusters or not. We perform a series of $N$-body simulations starting from initial conditions constructed from the results of hydrodynamical simulations of turbulent molecular clouds. In our simulations, both open clusters and young massive clusters form when we assume a density-dependent star formation efficiency. We find that a local star formation efficiency higher than 50 % is necessary for the formation of young massive clusters, but open clusters forms from less dense regions with a local star formation efficiency of $<50$ %. We confirm that the young massive clusters formed in our simulations have mass, size, and density profile similar to those of observed young massive clusters such as NGC 3603 and Trumpler 14. We also find that these simulated clusters evolve via hierarchical mergers of sub-clusters within a few Myr, as is suggested by recent simulations and observations. Although we do not assume initial mass segregation, we observe that the simulated massive clusters show a shallower slope of the mass function ($\Gamma\sim-1$) in the cluster center compared to that of the entire cluster ($\Gamma\sim-1.3$). These values are consistent with those of some young massive clusters in the Milky Way such as Westerlund 1 and Arches.

We determine the temporal evolution of the luminosity L(WD), radius R(WD) and effective temperature Teff of the white dwarf (WD) pseudophotosphere of V339 Del from its discovery to around day 40. Another main objective was studying the ionization structure of the ejecta. These aims were achieved by modelling the optical/near-IR spectral energy distribution (SED) using low-resolution spectroscopy (3500 - 9200 A), UBVRcIc and JHKLM photometry. During the fireball stage (Aug. 14.8 - 19.9, 2013), Teff was in the range of 6000 - 12000 K, R(WD) was expanding non-uniformly in time from around 66 to around 300 (d/3 kpc) R(Sun), and L(WD) was super-Eddington, but not constant. After the fireball stage, a large emission measure of 1.0-2.0E+62 (d/3 kpc)**2 cm**(-3) constrained the lower limit of L(WD) to be well above the super-Eddington value. The evolution of the H-alpha line and mainly the transient emergence of the Raman-scattered O VI 1032 A line suggested a biconical ionization structure of the ejecta with a disk-like H I region persisting around the WD until its total ionization, around day 40. It is evident that the nova was not evolving according to the current theoretical prediction. The unusual non-spherically symmetric ejecta of nova V339 Del and its extreme physical conditions and evolution during and after the fireball stage represent interesting new challenges for the theoretical modelling of the nova phenomenon.

Star clusters with multi-mass components dynamically evolve faster than those modeled with equal-mass components. Using a series of direct $N$-body simulations, we investigate the dynamical evolution of star clusters with mass functions, especially their core collapse time. Multi-mass clusters tend to behave like systems with a smaller number of particles, which we call the effective number of particles ($N_{eff}$) and for which $N_{eff} = M/m_{max}$ (here $M$ and $m_{max}$ are the total cluster mass and the mass of the most massive star in the cluster, respectively). We find that the time of core collapse is inversely proportional to the mass of the most massive star in the cluster and analytically confirm that this is because the core collapse of clusters with a mass function proceeds on the dynamical friction timescale of the most massive stars. As the mass of the most massive star increases, however, the core-collapse time, which is observed as a core bounce of the cluster core from the evolution of the core density or core radius, becomes ambiguous. We find that in that case the total binding energy of the hard binaries gives a good diagnosis for determining the moment of core collapses. Based on the results of our simulations, we argue that the core bounce becomes ambiguous when the mass of the most massive star exceeds 0.1\% of the total mass of the cluster.

Recent simulations and observations suggest that star clusters form via the assembling of smaller sub-clusters. Because of their short relaxation time, sub-clusters experience core collapse much earlier than virialized solo-clusters, which have similar properties of the merger remnant of the assembling clusters. As a consequence it seems that the assembling clusters result in efficient multiple collisions of stars in the cluster core. We performed a series of $N$-body simulations of ensemble and solitary clusters including stellar collisions and found that the efficiency of multiple collisions between stars are suppressed if sub-clusters assemble after they experience core collapse individually. In this case, sub-clusters form their own multiple collision stars which experienced a few collisions, but they fail to collide with each other after their host sub-clusters assemble. The multiple collision stars scatter each other and escape, and furthermore the central density of the remnant clusters had already been depleted for the stars to experience more collisions. On the other hand, if sub-clusters assemble before they experience core collapse, the multiple collisions of stars proceed efficiently in the remnant cluster, and the collision products are more massive than virialized solo-clusters and comparable in mass to cold solo-clusters.

We have found a radio-wave-reflection effect in rock salt for the detection of ultra-high energy neutrinos which are expected to be generated in Greisen, Zatsepin, and Kuzmin (GZK) processes in the universe. When an UHE neutrino interacts with rock salt or ice as a detection medium, a shower is generated. That shower is formed by hadronic and electromagnetic avalanche processes. The energy of the UHE neutrino shower converts to thermal energy through ionization processes. Consequently, the temperature rises along the shower produced by the UHE neutrino. The refractive index of the medium rises with temperature. The irregularity of the refractive index in the medium leads to a reflection of radio waves. This reflection effect combined with the long attenuation length of radio waves in rock salt and ice would yield a new method to detect UHE neutrinos. We measured the phase of the reflected radio wave under irradiation with an electron beam on ice and rock salt powder. The measured phase showed excellent consistence with the power reflection fraction which was measured directly. A model taking into account the temperature change explained the phase and the amplitude of the reflected wave. Therefore the reflection mechanism was confirmed. The power reflection fraction was compared with that calculated with the Fresnel equations, the ratio between the measured result and that obtained with the Fresnel equations in ice was larger than that of rock salt.

We present multi-color light curves and optical spectra of V1280 Scorpii obtained from 2007 to 2012. It is shown that V1280 Sco is the extremely slow nova and the mass of white dwarf appears to be $\sim$ 0.6 M$\odot$ or lower. Blue-shifted multiple absorption lines of Na {\sc i} D, Ca {\sc ii} HK, and He {\sc i*} are detected on high-resolution spectra. We also discuss that an approach using metastable He absorption lines is useful to investigate structures of nova shells.

We present visible spectroscopic observations of a near-earth object (NEO) 2012 DA14. The asteroid 2012 DA14 came close to the surface of the Earth on February 15, 2013 at a distance of 27,700 km. Its estimated diameter is around 45 m. The physical properties of such a small asteroid have not yet been well determined. The close encounter is a good opportunity to conduct a variety of observations. The purpose of this paper is to deduce the taxonomy of 2012 DA14 by visible spectroscopic observations using the 0.4 m f/10 telescope at the Fujii Kurosaki Observatory. We conclude that the taxonomy of 2012 DA14 is an L-type in the visible wavelength region. In addition, we refer to the availability of a small, accessible telescope for NEOs smaller than 100 m.

Most stars in galactic disks are believed to be born as a member of star clusters or associations. Star clusters formed in disks are disrupted due to the tidal stripping and the evolution of star clusters themselves, and as a results new stars are supplied to the galactic disks. We performed $N$-body simulations of star clusters in galactic disks, in which both star clusters and galactic disks are modeled as $N$-body ("live") systems, and as a consequence the disks form transient and recurrent spiral arms. In such non-steady spiral arms, star clusters migrate radially due to the interaction with spiral arms. We found that the migration timescale is a few hundreds Myr and that the angular momentum changes of star clusters are at most $\sim 50$% in 1 Gyr. Radial migration of star clusters to the inner region of galaxies results in a fast disruption of the star clusters because of a stronger tidal field in the inner region of the galaxy. This effect is not negligible for the disruption timescale of star clusters in galactic disks. Stars stripped from clusters form tidal tails which spread over 1--2 kpc. While the spatial distribution of tidal tails change in a complicated way due to the non-steady spiral arms, the velocity distribution conserve well even if the tidal tails are located at a few kpc from their parent clusters. Tidal tails of clusters in galactic disks might be detected using velocity plots.

Young star clusters like R136 in the Large Magellanic Cloud and NGC 3603, Westerlund 1, and 2 in the Milky Way are dynamically more evolved than expected based on their current relaxation times. In particular, the combination of a high degree of mass segregation, a relatively low central density, and the large number of massive runaway stars in their vicinity are hard to explain with the monolithic formation of these clusters. Young star clusters can achieve such a mature dynamical state if they formed through the mergers of a number of less massive clusters. The shorter relaxation times of less massive clusters cause them to dynamically evolve further by the time they merge, and the merger product preserves the memory of the dynamical evolution of its constituent clusters. With a series of $N$-body simulations, we study the dynamical evolution of single massive clusters and those that are assembled through merging smaller clusters together. We find that the formation of massive star clusters through the mergers of smaller clusters can reproduce the currently observed spatial distribution of massive stars, the density, and the characteristics (number and mass distribution) of the stars ejected as runaways from young dense clusters. We therefore conclude that these clusters and possibly other young massive star clusters formed through the mergers of smaller clusters.

We present optical ($B$, $V$, $R_{\rm c}$, $I_{\rm c}$ and $y$) and near infrared ($J$, $H$ and $K_{\rm s}$) photometric and spectroscopic observations of a classical nova V1280 Scorpii for five years from 2007 to 2011. Our photometric observations show a declining event in optical bands shortly after the maximum light which continues $\sim$ 250 days. The event is most probably caused by a dust formation. The event is accompanied by a short ($\sim$ 30 days) re-brightening episode ($\sim$ 2.5 mag in $V$), which suggests a re-ignition of the surface nuclear burning. After 2008, the $y$ band observations show a very long plateau at around $y$ = 10.5 for more than 1000 days until April 2011 ($\sim$ 1500 days after the maximum light). The nova had taken a very long time ($\sim$ 50 months) before entering the nebular phase (clear detection of both [\ion{O}{iii}] 4959 and 5007) and is still continuing to generate the wind caused by H-burning. The finding suggests that V1280 Sco is going through the historically slowest evolution. The interval from the maximum light (2007 February 16) to the beginning of the nebular phase is longer than any previously known slow novae: V723 Cas (18 months), RR Pic (10 months), or HR Del (8 months). It suggests that the mass of a white dwarf in the V1280 Sco system might be 0.6 $M_\mathrm{\sun}$ or smaller. The distance, based on our measurements of the expansion velocity combined with the directly measured size of the dust shell, is estimated to be 1.1 $\pm$ 0.5 kpc.

About 20% of all massive stars in the Milky Way have unusually high velocities, the origin of which has puzzled astronomers for half a century. We argue that these velocities originate from strong gravitational interactions between single stars and binaries in the centers of star clusters. The ejecting binary forms naturally during the collapse of a young ($\aplt 1$\,Myr) star cluster. This model replicates the key characteristics of OB runaways in our galaxy and it explains the $\apgt 100$\,\Msun\, runaway stars around young star clusters, e.g. R136 and Westerlund~2. The high proportion and the distributions in mass and velocity of runaways in the Milky Way is reproduced if the majority of massive stars are born in dense and relatively low-mass (5000-10000 \Msun) clusters.

We report on the spectral evolution of the enigmatic, very slow nova V5558 Sgr, based on the low-resolution spectra obtained at the Fujii-Bisei Observatory and the Bisei Astronomical Observatory, Japan during a period of 2007 April 6 to 2008 May 3. V5558 Sgr shows a pre-maximum halt and then several flare-like rebrightenings, which is similar to another very slow nova V723 Cas. In our observations, the spectral type of V5558 Sgr evolved from the He/N type toward the Fe II type during the pre-maximum halt, and then toward the He/N type again. This course of spectral transition was observed for the first time in the long history of the nova research. In the rebrightening stage after the initial brightness maximum, we could identify many emission lines accompanied by a stronger absorption component of the P-Cygni profile at the brightness maxima. We found that the velocity of the P-Cygni absorption component measured from the emission peak decreased at the brightness maxima. Furthermore, we compared the spectra of V5558 Sgr with V723 Cas, and other novae which exhibited several rebrightenings during the early phase.

It has been believed that spirals in pure stellar disks, especially the ones spontaneously formed, decay in several galactic rotations due to the increase of stellar velocity dispersions. Therefore, some cooling mechanism, for example dissipational effects of the interstellar medium, was assumed to be necessary to keep the spiral arms. Here we show that stellar disks can maintain spiral features for several tens of rotations without the help of cooling, using a series of high-resolution three-dimensional $N$-body simulations of pure stellar disks. We found that if the number of particles is sufficiently large, e.g., $3\times 10^6$, multi-arm spirals developed in an isolated disk can survive for more than 10 Gyrs. We confirmed that there is a self-regulating mechanism that maintains the amplitude of the spiral arms. Spiral arms increase Toomre's $Q$ of the disk, and the heating rate correlates with the squared amplitude of the spirals. Since the amplitude itself is limited by the value of $Q$, this makes the dynamical heating less effective in the later phase of evolution. A simple analytical argument suggests that the heating is caused by gravitational scattering of stars by spiral arms, and that the self-regulating mechanism in pure-stellar disks can effectively maintain spiral arms on a cosmological timescale. In the case of a smaller number of particles, e.g., $3\times 10^5$, spiral arms grow faster in the beginning of the simulation (while $Q$ is small) and they cause a rapid increase of $Q$. As a result, the spiral arms become faint in several Gyrs.

We report on the spectral evolution of 6 classical novae, V1186 Sco, V2540 Oph, V4745 Sgr, V5113 Sgr, V458 Vul, and V378 Ser, based on the low-resolution spectra obtained at the Fujii-Bisei Observatory and the Bisei Astronomical Observatory, Japan. In the light curves, these 6 novae show several rebrightenings during the early phase lasting ~10 days after the first maximum in fast novae, and ~100 days in slow novae. The early spectra of all of these novae had emission lines with a P-Cygni profile at the maximum brightness. The absorption component of the P-Cygni profiles then disappeared after the maximum, and reappeared when the novae brightened to the next maximum. We suggest that the re-appearance of the absorption component at the rebrightening is attributable to re-expansion of the photosphere after it once shifts sufficiently inside. From the light curves, we found that the time intervals of the rebrightenings of these 6 novae show a similar systematic trend, which is applicable to all types of novae: fast and slow, and Fe II type and hybrid type. Moreover, we note the difference between the spectra at the rebrightenings during the early phase and at the rebrightening in V2362 Cyg, and at the oscillation during the transition phase in V1494 Aql, which means difference of the physical mechanism of the rebrightening during the early phase and the later oscillations.

We present early spectroscopy of the recurrent nova U~Sco during the outburst in 2010. We successfully obtained time-series spectra at $t_{\rm d}=$0.37--0.44~d, where $t_{\rm d}$ denotes the time from the discovery of the present outburst. This is the first time-resolved spectroscopy on the first night of U Sco outbursts. At $t_{\rm d}\sim 0.4$~d the H$\alpha$ line consists of a blue-shifted ($-5000$ km s$^{-1}$) narrow absorption component and a wide emission component having triple peaks, a blue ($\sim -3000$ km s$^{-1}$), a central ($\sim 0$ km s$^{-1}$) and a red ($\sim +3000$ km s$^{-1}$) ones. The blue and red peaks developed more rapidly than the central one during the first night. This rapid variation would be caused by the growth of aspherical wind produced during the earliest stage of the outburst. At $t_{\rm d}=1.4$~d the H$\alpha$ line has a nearly flat-topped profile with weak blue and red peaks at $\sim \pm 3000$ km s$^{-1}$. This profile can be attributed to a nearly spherical shell, while the asphericity growing on the first night still remains. The wind asphericity is less significant after $t_{\rm d}=9$ d.

Within the distance of 1 pc from the Galactic center (GC), more than 100 young massive stars have been found. The massive stars at 0.1-1 pc from the GC are located in one or two disks, while those within 0.1 pc from the GC, S-stars, have an isotropic distribution. How these stars are formed is not well understood, especially for S-stars. Here we propose that a young star cluster with an intermediate-mass black hole (IMBH) can form both the disks and S-stars. We performed a fully self-consistent $N$-body simulation of a star cluster near the GC. Stars escaped from the tidally disrupted star cluster were carried to the GC due to an 1:1 mean motion resonance with the IMBH formed in the cluster. In the final phase of the evolution, the eccentricity of the IMBH becomes very high. In this phase, stars carried by the 1:1 resonance with the IMBH were dropped from the resonance and their orbits are randomized by a chaotic Kozai mechanism. The mass function of these carried stars is extremely top-heavy within 10". The surface density distributions of young massive stars has a slope of -1.5 within 10" from the GC. The distribution of stars in the most central region is isotropic. These characteristics agree well with those of stars observed within 10" from the GC.

We carried out an international spectroscopic observation campaign of the dwarf nova GW Librae (GW Lib) during the 2007 superoutburst. Our observation period covered the rising phase of the superoutburst, maximum, slowly decaying phase (plateau), and long fading tail after the rapid decline from the plateau. The spectral features dramatically changed during the observations. In the rising phase, only absorption lines of H$\alpha$, H$\beta$, and H$\gamma$ were present. Around the maximum, the spectrum showed singly-peaked emission lines of H$\alpha$, He I 5876, He I 6678, He II 4686, and C III/N III as well as absorption lines of Balmer components and He I. These emission lines significantly weakened in the latter part of the plateau phase. In the fading tail, all the Balmer lines and He I 6678 were in emission, as observed in quiescence. We find that the center of the H$\alpha$ emission component was mostly stable over the whole orbital phase, being consistent with the low inclination of the system. Comparing with the observational results of WZ Sge during the 2001 superoutburst, the same type of stars as GW Lib seen with a high inclination angle, we interpret that the change of the H$\alpha$ profile before the fading tail phase is attributed to a photoionized region formed at the outer edge of the accretion disk, irradiated from the white dwarf and inner disk.

We performed, for the first time, the simulation of spiral-in of a star cluster formed close to the Galactic center (GC) using a fully self-consistent $N$-body model. In our model, the central super-massive black hole (SMBH) is surrounded by stars and the star cluster. Not only are the orbits of stars and the cluster stars integrated self-consistently, but the stellar evolution, collisions and merging of the cluster stars are also included. We found that an intermediate-mass black hole (IMBH) is formed in the star cluster and stars escaped from the cluster are captured into a 1:1 mean motion resonance with the IMBH. These "Trojan" stars are brought close to the SMBH by the IMBH, which spirals into the GC due to the dynamical friction. Our results show that, once the IMBH is formed, it brings the massive stars to the vicinity of the central SMBH even after the star cluster itself is disrupted. Stars carried by the IMBH form a disk similar to the observed disks and the core of the cluster including the IMBH has properties similar to those of IRS13E, which is a compact assembly of several young stars.

We present MUSE, a software framework for combining existing computational tools for different astrophysical domains into a single multiphysics, multiscale application. MUSE facilitates the coupling of existing codes written in different languages by providing inter-language tools and by specifying an interface between each module and the framework that represents a balance between generality and computational efficiency. This approach allows scientists to use combinations of codes to solve highly-coupled problems without the need to write new codes for other domains or significantly alter their existing codes. MUSE currently incorporates the domains of stellar dynamics, stellar evolution and stellar hydrodynamics for studying generalized stellar systems. We have now reached a "Noah's Ark" milestone, with (at least) two available numerical solvers for each domain. MUSE can treat multi-scale and multi-physics systems in which the time- and size-scales are well separated, like simulating the evolution of planetary systems, small stellar associations, dense stellar clusters, galaxies and galactic nuclei. In this paper we describe three examples calculated using MUSE: the merger of two galaxies, the merger of two evolving stars, and a hybrid N-body simulation. In addition, we demonstrate an implementation of MUSE on a distributed computer which may also include special-purpose hardware, such as GRAPEs or GPUs, to accelerate computations. The current MUSE code base is publicly available as open source at http://muse.li

We carried out intensive spectroscopic observations of two WZ Sge-type dwarf novae, GW Lib, and V455 And during their superoutbursts in 2007, at 6 observatories. The observations covered the whole of both superoutbursts from the very maximum to the fading tail. We found evidence of the winds having a speed of $\sim$1000 km s$^{-1}$ which blew in GW Lib during the rising phase. The evolution of the hydrogen, helium, and carbon lines suggests flaring of the accretion disk and emergence of the temperature inversion layer on the disk.

We have performed fully self-consistent $N$-body simulations of star clusters near the Galactic center (GC). Such simulations have not been performed because it is difficult to perform fast and accurate simulations of such systems using conventional methods. We used the Bridge code, which integrates the parent galaxy using the tree algorithm and the star cluster using the fourth-order Hermite scheme with individual timestep. The interaction between the parent galaxy and the star cluster is calculate with the tree algorithm. Therefore, the Bridge code can handle both the orbital and internal evolutions of star clusters correctly at the same time. We investigated the evolution of star clusters using the Bridge code and compared the results with previous studies. We found that 1) the inspiral timescale of the star clusters is shorter than that obtained with "traditional" simulations, in which the orbital evolution of star clusters is calculated analytically using the dynamical friction formula and 2) the core collapse of the star cluster increases the core density and help the cluster survive. The initial conditions of star clusters is not so severe as previously suggested.

We report the observation of the first gravitational microlensing event in a sparse stellar field, involving the brightest (V=11.4 mag) andclosest (~ 1 kpc) source star to date. This event was discovered by an amateurastronomer, A. Tago, on 2006 October 31 as a transient brightening, by ~4.5 mag during a ~15 day period, of a normal A-type star (GSC 3656-1328) in the Cassiopeia constellation. Analysis of both spectroscopic observations and the light curve indicates that this event was caused by gravitational microlensing rather than an intrinsically variable star. Discovery of this single event over a 30 year period is roughly consistent with the expected microlensing rate for the whole sky down to V = 12 mag stars. However, the probability for finding events with such a high magnification (~ 50) is much smaller, by a factor ~1/50, which implies that the true event rate may be higher than expected. This discovery indicates the potential of all sky variability surveys, employing frequent sampling by telescopes with small apertures and wide fields of view, for finding such rare transient events, and using the observations to explore galactic disk structure and search for exo-planets.

Existence of GZK neutrinos (ultra high energy neutrinos) have been justified although the flux is very low. A new method is desired to use a huge mass of a detector medium to detect them. A fundamental study of radar method was carried out to measure microwave reflection from electromagnetic energy deposit by X-ray irradiation in a small rock salt sample. The reflection rate of 1x10^-6 was found at the energy deposit of 1x10^19 eV which was proportional to square of the X-ray intensity suggesting the effect to be coherent scattering. The decay time of the reflection was several seconds. This effect implies a large scale natural rock salt formation could be utilized like a bubble chamber irradiated by radio wave instead of visible light to detect GZK neutrinos.

We developed a new direct-tree hybrid N-body algorithm for fully self-consistent N-body simulations of star clusters in their parent galaxies. In such simulations, star clusters need high accuracy, while galaxies need a fast scheme because of the large number of the particles required to model it. In our new algorithm, the internal motion of the star cluster is calculated accurately using the direct Hermite scheme with individual timesteps and all other motions are calculated using the tree code with second-order leapfrog integrator. The direct and tree schemes are combined using an extension of the mixed variable symplectic (MVS) scheme. Thus, the Hamiltonian corresponding to everything other than the internal motion of the star cluster is integrated with the leapfrog, which is symplectic. Using this algorithm, we performed fully self-consistent N-body simulations of star clusters in their parent galaxy. The internal and orbital evolutions of the star cluster agreed well with those obtained using the direct scheme. We also performed fully self-consistent N-body simulation for large-N models ($N=2\times 10^6$). In this case, the calculation speed was seven times faster than what would be if the direct scheme was used.

For a rigid model satellite, Chandrasekhar's dynamical friction formula describes the orbital evolution quite accurately, when the Coulomb logarithm is chosen appropriately. However, it is not known if the orbital evolution of a real satellite with the internal degree of freedom can be described by the dynamical friction formula. We performed N-body simulation of the orbital evolution of a self-consistent satellite galaxy within a self-consistent parent galaxy. We found that the orbital decay of the simulated satellite is significantly faster than the estimate from the dynamical friction formula. The main cause of this discrepancy is that the stars stripped out of the satellite are still close to the satellite, and increase the drag force on the satellite through two mechanisms. One is the direct drag force from particles in the trailing tidal arm, a non-axisymmetric force that slows the satellite down. The other is the indirect effect that is caused by the particles remaining close to the satellite after escape. The force from them enhances the wake caused in the parent galaxy by dynamical friction, and this larger wake in turn slows the satellite down more than expected from the contribution of its bound mass. We found these two have comparable effects, and the combined effect can be as large as 20% of the total drag force on the satellite.

We analyze and fit visible spectra of a red horizon at sunrise. The shape of the spectra consist of a blue continuum followed by a red bump. The reddest spectra are well fitted by the product of a spectrum of extinguished sunlight (Rayleigh extinction + ozone absorption) and 1/lambda^4. The former is essentially the radiation field in the outer atmosphere, at the scattering volume location; the latter corresponds to Rayleigh scattering by the gas. Moving to higher altitudes, a second component, corresponding to the spectrum of a blue sky, must be added. The spectra we have obtained are similar to spectra of red nebulae, suggesting there may be other explanations than an emission process to the red color of some nebulae.

The minimal supersymmetric standard model (MSSM) has a truly supersymmetric way to explain both the baryon asymmetry and cold dark matter in the present Universe, that is, ``Affleck-Dine baryo/DM-genesis.'' The associated late-time decay of Q-balls directly connects the origins of the baryon asymmetry and dark matter, and also predicts a specific nature of the LSP. In this paper, we investigate the prospects for indirect detection of these dark matter candidates observing high energy neutrino flux from the Sun, and hard positron flux from the halo. We also update the previous analysis of the direct detection in hep-ph/0205044 by implementing the recent result from WMAP satellite.

We point out that there is no cosmological gravitino problem in a certain class of gauge-mediated supersymmetry-breaking (GMSB) models. The constant term in the superpotential naturally causes small mixings between the standard-model and messenger fields, which give rise to late-time decays of the lightest messenger fields. This decay provides an exquisite amount of entropy, which dilutes the thermal relics of the gravitinos down to just the observed mass density of the dark matter. This remarkable phenomenon takes place naturally, irrespective of the gravitino mass and the reheating temperature of inflation, once the gravitinos and messenger fields are thermalized in the early Universe. In this class of GMSB models, there is no strict upper bound on the reheating temperature of inflation, which makes the standard thermal leptogenesis the most attractive candidate for the origin of the observed baryon asymmetry in the present Universe.

We summarize photometric and spectroscopic observations of V463 Sct (Nova Sct 2000), which was originally thought to be a red variable. The spectrum taken on 2000 March 16.81 UT showed prominent emission lines with a FHWM of 990 km/s (H alpha). The light curve shows a conspicuous premaximum halt lasting at least for 24 d, and a late-phase flare-like maximum. The nova then started rapidly fading at a rate corresponding to t2 = 15 +/- 3 d. Long premaximum halts have been considered as a unique character of the "slowest" novae. The present observation, however, suggests that the long premaximum halts are not a unique character of the slowest novae, but a more general phenomenon spreading over a wider range of nova speed classes than has been previously believed. A recent interpretation of premaximum halts requires that the conditions of thermonuclear runaway was only marginally satisfied. Since such conditions are more difficult to meet in rapidly evolving novae, V463 Sct would provide an unique opportunity in testing this interpretation. The early post-outburst spectrum showed co-existence of Fe II lines and some forbidden lines, which suggests that substantial amount of material may have been ejected before the observed optical maximum. The impact of the modern global network (VSNET) on confirmatory processes of transient objects is briefly discussed.

The formation and late time decays of Q-balls are generic consequences of the Affleck-Dine (AD) baryogenesis. A substantial amount of the lightest supersymmetry (SUSY) particles (LSPs) are produced non-thermally as the decay products of these Q-balls. This requires a significantly large annihilation cross section of the LSP so as not to overclose the universe, which predicts a higgsino- or wino-like LSP instead of the standard bino LSP. We have reexamined the AD baryogenesis with special attention to the late-time decays of the Q-balls, and then specified the parameter regions where the LSPs produced by the Q-ball decays result in a cosmologically interesting mass density of dark matter by adopting several SUSY breaking models. This reveals new cosmologically interesting parameter regions, which have not attracted much attention so far. We have also investigated the prospects of direct and indirect detection of these dark matter candidates, and found that there is an intriguing possibility to detect them in various next generation dark matter searches.