Search

Studying the early stages of transient events provides crucial information about the fundamental physical processes in cataclysmic variables (CVs). However, determining an appropriate observation mode immediately after the discovery of a new transient presents challenges due to significant uncertainties regarding its nature. We developed a framework designed for autonomous decision making in prompt follow-up observations of CVs using the Kanata 1.5-m telescope. The system, named Smart Kanata, first estimates the class probabilities of variable star types using a generative model. It then selects the optimal observation mode from three possible options based on the mutual information calculated from the class probabilities. We have operated the system for ~300 days and obtained 21 samples, among which automated observations were successfully performed for a nova and a microlensing event. In the time-series spectra of the nova V4370 Oph, we detected a rapid deepening of the absorption component of the H_alpha line. These initial results demonstrate the capability of Smart Kanata in facilitating rapid observations and improving our understanding of outbursts and eruptions of CVs and other galactic transients., Comment: 18 pages, 12 figures, accepted for publication in PASJ

We present a detailed study of the MASTER OT J030227.28+191754.5 outburst in 2021-2022, reaching an amplitude of 10.2 mag and a duration of 60 d. The detections of (1) the double-peaked optical emission lines, and (2) the early and ordinary superhumps, established that MASTER OT J030227.28+191754.5 is an extremely energetic WZ Sge-type dwarf nova (DN). Based on the superhump observations, we obtained its orbital period and mass ratio as 0.05986(1) d and 0.063(1), respectively. These are within a typical range of low-mass-ratio DNe. According to the binary parameters derived based on the thermal-tidal instability model, our analyses showed that (1) the standard disk model requires an accretion rate $\simeq$ 10$^{20}$ g s$^{-1}$ to explain its peak optical luminosity and (2) large mass was stored in the disk at the outburst onset. These cannot be explained solely by the impact of its massive ($\gtrsim$ 1.15 M$_\odot$) primary white dwarf implied by Kimura et al. (2023). Instead, we propose that the probable origin of this enormously energetic DN outburst is the even lower quiescence viscosity than other WZ Sge-type DNe. This discussion is qualitatively valid for most possible binary parameter spaces unless the inclination is low ($\lesssim 40^\circ$) enough for the disk to be bright explaining the outburst amplitude. Such low inclinations, however, would not allow detectable amplitude of early superhumps in the current thermal-tidal instability model. The optical spectra at outburst maximum showed the strong emission lines of Balmer, He I, and He II series whose core is narrower than $\sim 800$ km s$^{-1}$. Considering its binary parameters, a Keplerian disk cannot explain this narrow component, but the presumable origin is disk winds., Comment: 23 pages, 13 figures, 2 tables. Accepted by PASJ. Part of the online supplemental information is included

We present photometric, spectroscopic and polarimetric observations of the intermediate-luminosity Type IIP supernova (SN) 2021gmj from 1 to 386 days after the explosion. The peak absolute V-band magnitude of SN 2021gmj is -15.5 mag, which is fainter than that of normal Type IIP SNe. The spectral evolution of SN 2021gmj resembles that of other sub-luminous supernovae: the optical spectra show narrow P-Cygni profiles, indicating a low expansion velocity. We estimate the progenitor mass to be about 12 Msun from the nebular spectrum and the 56Ni mass to be about 0.02 Msun from the bolometric light curve. We also derive the explosion energy to be about 3 x 10^{50} erg by comparing numerical light curve models with the observed light curves. Polarization in the plateau phase is not very large, suggesting nearly spherical outer envelope. The early photometric observations capture the rapid rise of the light curve, which is likely due to the interaction with a circumstellar material (CSM). The broad emission feature formed by highly-ionized lines on top of a blue continuum in the earliest spectrum gives further indication of the CSM at the vicinity of the progenitor. Our work suggests that a relatively low-mass progenitor of an intermediate-luminosity Type IIP SN can also experience an enhanced mass loss just before the explosion, as suggested for normal Type IIP SNe., Comment: 18 pages, 16 figures, resubmitted to MNRAS after addressing referee comments

We present early-phase panchromatic photometric and spectroscopic coverage spanning far-ultraviolet (FUV) to the near-infrared (NIR) regime of the nearest hydrogen-rich core-collapse supernova in the last 25 years, SN 2023ixf. We observe early 'flash' features in the optical spectra due to a confined dense circumstellar material (CSM). We observe high-ionization absorption lines (FeII, MgII) in the ultraviolet spectra from very early on. We also observe a multi-peaked emission profile of H-alpha in the spectrum beginning ~16 d, which indicates ongoing interaction of the SN ejecta with a pre-existing shell-shaped CSM having an inner radius of ~75 AU and an outer radius of ~140 AU. The shell-shaped CSM is likely a result of enhanced mass loss ~35-65 years before the explosion assuming a standard Red-Supergiant wind. The UV spectra are dominated by multiple highly ionized narrow absorption features and broad emission features from elements such as C, N, O, Si, Fe, and Ni. Based on early light curve models of Type II SNe, we infer that the nearby dense CSM confined to (7+-3)e14cm (~45 AU) is a result of enhanced mass loss (10^{-3.0+-0.5} Msol/yr) two decades before the explosion., Comment: Accepted for publication in ApJL (Revised, 3 figures and 1 table)

We present an extensive, panchromatic photometric (UV, Optical, and NIR) and low-resolution optical spectroscopic coverage of a Type IIP supernova SN 2018gj that occurred on the outskirts of the host galaxy NGC 6217. From the V-band light curve, we estimate the plateau length to be ~ 70 +- 2 d, placing it among the very few well-sampled short plateau supernovae (SNe). With V-band peak absolute magnitude Mv < -17.0 +- 0.1 mag, it falls in the middle of the luminosity distribution of the Type II SNe. The colour evolution is typical to other Type II SNe except for an early elbow-like feature in the evolution of V-R colour owing to its early transition from the plateau to the nebular phase. Using the expanding photospheric method, we present an independent estimate of the distance to SN 2018gj. We report the spectral evolution to be typical of a Type II SNe. However, we see a persistent blue shift in emission lines until the late nebular phase, not ordinarily observed in Type II SNe. The amount of radioactive nickel (56Ni) yield in the explosion was estimated to be 0.026 +- 0.007 Msol. We infer from semi-analytical modelling, nebular spectrum, and 1-D hydrodynamical modelling that the probable progenitor was a red supergiant with a zero-age-main-sequence mass < 13 Msol. In the simulated hydrodynamical model light curves, reproducing the early optical bolometric light curve required an additional radiation source, which could be the interaction with the proximal circumstellar matter (CSM)., Comment: Accepted for publication in ApJ (31 pages, 23 figures and 7 tables)

Type IIn/Ia-CSM supernovae (SNe IIn/Ia-CSM) are classified by their characteristic spectra, which exhibit narrow hydrogen emission lines originating from a strong interaction with a circumstellar medium (CSM) together with broad lines of intermediate-mass elements. We performed intensive follow-up observations of SN IIn/Ia-CSM 2020uem, including photometry, spectroscopy, and polarimetry. In this paper, we focus on the results of polarimetry. We performed imaging polarimetry at $66$ days and spectropolarimetry at $103$ days after the discovery. SN 2020uem shows a high continuum polarization of $1.0-1.5\%$ without wavelength dependence. Besides, the polarization degree and position angle keep roughly constant. These results suggest that SN 2020uem is powered by a strong interaction with a confined and aspherical CSM. We performed a simple polarization modeling, based on which we suggest that SN 2020uem has an equatorial-disk/torus CSM. Besides, we performed semi-analytic light-curve modeling and estimated the CSM mass. We revealed that the mass-loss rate in the final few hundred years immediately before the explosion of SN 2020uem is in the range of $0.01 - 0.05 {\rm ~M_{\odot}~yr^{-1}}$, and that the total CSM mass is $0.5-4 {\rm ~M_{\odot}}$. The CSM mass can be accommodated by not only a red supergiant (RSG) but a red giant (RG) or an asymptotic-giant-branch (AGB) star. As a possible progenitor scenario of SN 2020uem, we propose a white-dwarf binary system including an RG, RSG or AGB star, especially a merger scenario via common envelope evolution, i.e., the core-degenerate scenario or its variant., Comment: 16pages, 9 figures, 1 table. Accepted for publication in ApJ

We have performed intensive follow-up observations of a Type IIn/Ia-CSM SN (SN IIn/Ia-CSM), 2020uem, with photometry, spectroscopy, and polarimetry. In this paper, we report on the results of our observations focusing on optical/near-infrared (NIR) photometry and spectroscopy. The maximum V-band magnitude of SN 2020uem is over $-19.5$ mag. The light curves decline slowly with a rate of $\sim 0.75 {\rm ~mag}/100 {\rm ~days}$. In the late phase ($\gtrsim 300$ days), the light curves show accelerated decay ($\sim 1.2 {\rm ~mag}/100 {\rm ~days}$). The optical spectra show prominent hydrogen emission lines and broad features possibly associated with Fe-peak elements. In addition, the $\rm H\alpha$ profile exhibits a narrow P-Cygni profile with the absorption minimum of $\sim 100 {\rm ~km~s^{-1}}$. SN 2020uem shows a higher $\rm H\alpha/H\beta$ ratio ($\sim 7$) than those of SNe IIn, which suggests a denser CSM. The NIR spectrum shows the Paschen and Brackett series with continuum excess in the H and Ks bands. We conclude that the NIR excess emission originates from newly-formed carbon dust. The dust mass ($M_{\rm d}$) and temperature ($T_{\rm d}$) are derived to be $(M_{\rm d}, T_{\rm d}) \sim (4-7 \times 10^{-5} {\rm ~M_{\odot}}, 1500-1600 {\rm ~K})$. We discuss the differences and similarities between the observational properties of SNe IIn/Ia-CSM and those of other SNe Ia and interacting SNe. In particular, spectral features around $\sim 4650$ {\text \AA} and $\sim 5900$ {\text \AA} of SNe IIn/Ia-CSM are more suppressed than those of SNe Ia; these lines are possibly contributed, at least partly, by \ion{Mg}{1}] and \ion{Na}{1}, and may be suppressed by high ionization behind the reverse shock caused by the massive CSM., Comment: 22 pages, 15 figures, 3 tables. Accepted for publication in ApJ

In this Letter we report a discovery of a prominent flash of a peculiar overluminous Type Ia supernova, SN 2020hvf, in about 5 hours of the supernova explosion by the first wide-field mosaic CMOS sensor imager, the Tomo-e Gozen Camera. The fast evolution of the early flash was captured by intensive intranight observations via the Tomo-e Gozen high-cadence survey. Numerical simulations show that such a prominent and fast early emission is most likely generated from an interaction between $0.01~M_{\odot}$ circumstellar material (CSM) extending to a distance of $\sim$$10^{13}~\text{cm}$ and supernova ejecta soon after the explosion, indicating a confined dense CSM formation at the final evolution stage of the progenitor of SN 2020hvf. Based on the CSM-ejecta interaction-induced early flash, the overluminous light curve, and the high ejecta velocity of SN 2020hvf, we suggest that the SN 2020hvf may originate from a thermonuclear explosion of a super-Chandrasekhar-mass white dwarf ("super-$M\rm_{Ch}$ WD"). Systematical investigations on explosion mechanisms and hydrodynamic simulations of the super-$M\rm_{Ch}$ WD explosion are required to further test the suggested scenario and understand the progenitor of this peculiar supernova., Comment: 20pages, 6 figures, accepted for publication in the Astrophysical Journal Letters

We report photometric and spectroscopic observations of the eclipsing SU UMa-type dwarf nova ASASSN-18aan. We observed the 2018 superoutburst with 2.3 mag brightening and found the orbital period ($P_{\rm orb}$) to be 0.149454(3) d, or 3.59 hr. This is longward of the period gap, establishing ASASSN-18aan as one of a small number of long-$P_{\rm orb}$ SU UMa-type dwarf novae. The estimated mass ratio, ($q=M_2/M_1 = 0.278(1)$), is almost identical to the upper limit of tidal instability by the 3:1 resonance. From eclipses, we found that the accretion disk at the onset of the superoutburst may reach the 3:1 resonance radius, suggesting that the superoutburst of ASASSN-18aan results from the tidal instability. Considering the case of long-$P_{\rm orb}$ WZ Sge-type dwarf novae, we suggest that the tidal dissipation at the tidal truncation radius is enough to induce SU UMa-like behavior in relatively high-$q$ systems such as SU UMa-type dwarf novae, but that this is no longer effective in low-$q$ systems such as WZ Sge-type dwarf novae. The unusual nature of the system extends to the secondary star, for which we find a spectral type of G9, much earlier than typical for the orbital period, and a secondary mass $M_2$ of around 0.18 M$_{\odot}$, smaller than expected for the orbital period and the secondary's spectral type. We also see indications of enhanced sodium abundance in the secondary's spectrum. Anomalously hot secondaries are seen in a modest number of other CVs and related objects. These systems evidently underwent significant nuclear evolution before the onset of mass transfer. In the case of ASASSN-18aan, this apparently resulted in a mass ratio lower than typically found at the system's $P_{\rm orb}$, which may account for the occurrence of a superoutburst at this relatively long period., Comment: 17 pages, 14 figures, accepted for publication in PASJ

We report on the multi-wavelength photometry of the 2018 superoutburst in EG Cnc. We have detected stage A superhumps and long-lasting late-stage superhumps via the optical photometry and have constrained the binary mass ratio and its possible range. The median value of the mass ratio is 0.048 and the upper limit is 0.057, which still implies that EG Cnc is one of the possible candidates for the period bouncer. This object also showed multiple rebrightenings in this superoutburst, which are the same as those in its previous superoutburst in 1996--1997 despite the difference in the main superoutburst. This would represent that the rebrightening type is inherent to each object and is independent of the initial disk mass at the beginning of superoutbursts. We also found that $B-I$ and $J-K_{\rm S}$ colors were unusually red just before the rebrightening phase and became bluer during the quiescence between rebrightenings, which would mean that the low-temperature mass reservoir at the outermost disk accreted with time after the main superoutburst. Also, the ultraviolet flux was sensitive to rebrightenings as well as the optical flux, and the $U-B$ color became redder during the rebrightening phase, which would indicate that the inner disk became cooler when this object repeated rebrightenings. Our results thus basically support the idea that the cool mass reservoir in the outermost disk is responsible for rebrightenings., Comment: Accepted for publication in PASJ. 13 pages and 9 figures

PTFO 8-8695b (CVSO 30b) is a young planet candidate whose host star is a $\sim$ 2.6 Myr-old T-Tauri star, and there have been continuous discussions about the nature of this system. To unveil the mystery of this system, we observed PTFO8-8695 for around five years at optical and infrared bands simultaneously using Kanata telescope at the Higashi-Hiroshima Observatory. Through our observations, we found that the reported fading event split into two: deeper but phase-shifted "dip-A" and shallower but equiphase "dip-B". These dips disappeared at different epochs, and then, dip-B reappeared. Based on the observed wavelength dependence of dip depths, a dust clump and a precessing planet are likely origins of dip-A and B, respectively. Here we propose "a precessing planet associated with a dust cloud" scenario for this system. This scenario is consistent with the reported change in the depth of fading events, and even with the reported results, which were thought to be negative evidence to the planetary hypothesis, such as the past non-detection of the Rossiter-McLaughlin effect. If this scenario is correct, this is the third case of a young (<3 Myr) planet around a pre-main sequence star. This finding implies that a planet can be formed within a few Myr., Comment: 17 pages, 19 figures, 6 tables, accepted for publication in PASJ

We followed-up the massive young stellar object (MYSO) S255-NIRS3 (=S255-IRS1b) during its recent accretion outburst event in the Ks band with Kanata/HONIR for four years after its burst and obtained a long-term light curve. This is the most complete NIR light-curve of the S255-NIRS3 burst event that has ever been presented. The light curve showed a steep increase reaching a peak flux that was 3.4 mag brighter than the quiescent phase and then a relatively moderate year-scale fading until the last observation, similar to that of the accretion burst events such as EXors found in lower-mass young stellar objects. The behavior of the Ks band light curve is similar to that observed in 6.7 GHz class II methanol maser emission, with a sudden increase followed by moderate year-scale fading. However, the maser emission peaks appear 30-50 days earlier than that of the Ks band emission. The similarities confirmed that the origins of the maser emission and the Ks band continuum emission is common as previously shown from another infrared and radio observations by Stecklum et al. (2016); Caratti o Garatti et al. (2017a); Moscadelli et al. (2017). However, the differences in energy transfer paths, such as the exciting/emitting/scattering structures, may cause the delay in the flux-peak dates., Comment: 8 pages, 4 figures, 2 tables, accepted to PASJ

High cadence ultraviolet, optical and near-infrared photometric and low-resolution spectroscopic observations of the peculiar Type II supernova (SN) 2018hna are presented. The early phase multiband light curves exhibit the adiabatic cooling envelope emission following the shock breakout up to ~14 days from the explosion. SN~2018hna has a rise time of $\sim$\,88 days in the V-band, similar to SN 1987A. A $\rm^{56}Ni$ mass of ~0.087$\pm$0.004 $\rm M_{\odot}$ is inferred for SN 2018hna from its bolometric light curve. Hydrodynamical modelling of the cooling phase suggests a progenitor with a radius ~50 $\rm R_{\odot}$, a mass of ~14-20 $\rm M_{\odot}$ and explosion energy of ~1.7-2.9$\rm \times$ $\rm 10^{51}\ erg$. The smaller inferred radius of the progenitor than a standard red supergiant is indicative of a blue supergiant progenitor of SN 2018hna. A sub-solar metallicity (~0.3 $\rm Z_{\odot}$) is inferred for the host galaxy UGC 07534, concurrent with the low-metallicity environments of 1987A-like events., Comment: 8 pages, 5 figures, Accepted for Publication in ApJ Letters

We present the results based on photometric ($Swift$ UVOT), broad-band polarimetric ($V$ and $R$-band) and optical spectroscopic observations of the Type IIn supernova (SN) 2017hcc. Our study is supplemented with spectropolarimetric data available in literature for this event. The post-peak light curve evolution is slow ($\sim$0.2 mag 100 d$^{-1}$ in $b$-band). The spectrum of $\sim$+27 d shows a blue continuum with narrow emission lines, typical of a Type IIn SN. Archival polarization data along with the $Gaia$ DR2 distances have been utilized to evaluate the interstellar polarization (ISP) towards the SN direction which is found to be $P_{ISP}$ = 0.17 $\pm$ 0.02 per cent and $\theta_{ISP}$ = 140$^{\circ}$ $\pm$ 3$^{\circ}$. To extract the intrinsic polarization of SN 2017hcc, both the observed and the literature polarization measurements were corrected for ISP. We noticed a significant decline of $\sim$3.5 per cent ($V$-band) in the intrinsic level of polarization spanning a period of $\sim$2 months. In contrast, the intrinsic polarization angles remain nearly constant at all epochs. Our study indicates a substantial variation in the degree of asymmetry in either the ejecta and/or the surrounding medium of SN 2017hcc. We also estimate a mass-loss rate of $\dot M$ = 0.12 M$_{\odot}$ yr$^{-1}$ (for $v_w$ = 20 km s$^{-1}$) which suggests that the progenitor of SN 2017hcc is most likely a Luminous Blue Variable., Comment: 11 pages, 7 figures, 4 tables, Accepted for publication in MNRAS

We present optical and near-infrared observations of a low-luminosity Type IIP supernova (SN) 2016bkv from the initial rising phase to the plateau phase. Our observations show that the end of the plateau is extended to $\gtrsim 140$ days since the explosion, indicating that this SN takes one of the longest time to finish the plateau phase. among Type IIP SNe (SNe IIP), including low-luminosity (LL) SNe IIP. The line velocities of various ions at the middle of the plateau phase are as low as 1,000--1,500 km s$^{-1}$, which is the lowest even among LL SNe IIP. These measurements imply that the ejecta mass in SN 2016bkv is larger than that of the well-studied LL IIP SN 2003Z. In the early phase, SN 2016bkv shows a strong bump in the light curve. In addition, the optical spectra in this bump phase exhibit a blue continuum accompanied with a narrow H$\alpha$ emission line. These features indicate an interaction between the SN ejecta and the circumstellar matter (CSM) as in SNe IIn. Assuming the ejecta-CSM interaction scenario, the mass loss rate is estimated to be $\sim 1.7 \times 10^{-2} M_{\odot}$ yr$^{-1}$ within a few years before the SN explosion. This is comparable to or even larger than the largest mass loss rate observed for the Galactic red supergiants ($\sim 10^{-3} M_{\odot}$ yr$^{-1}$ for VY CMa). We suggest that the progenitor star of SN 2016bkv experienced a violent mass loss just before the SN explosion., Comment: 9 pages, 12 figures. Accepted for publication in ApJ

We present the early-phase spectra and the light curves of the broad-lined supernova (SN) 2016coi from $t=7$ to $67$ days after the estimated explosion date. This SN was initially reported as a broad-lined Type SN Ic (SN Ic-BL). However, we found that spectra up to $t=12$ days exhibited the He~{\sc i} $\lambda$5876, $\lambda$6678, and $\lambda$7065 absorption lines. We show that the smoothed and blueshifted spectra of normal SNe Ib are remarkably similar to the observed spectrum of SN 2016coi. The line velocities of SN 2016coi were similar to those of SNe Ic-BL and significantly faster than those of SNe Ib. Analyses of the line velocity and light curve suggest that the kinetic energy and the total ejecta mass of SN 2016coi are similar to those of SNe Ic-BL. Together with broad-lined SNe 2009bb and 2012ap for which the detection of He~{\sc i} were also reported, these SNe could be transitional objects between SNe Ic-BL and SNe Ib, and be classified as broad-lined Type `Ib' SNe (SNe `Ib'-BL). Our work demonstrates the diversity of the outermost layer in broad-lined SNe, which should be related to the variety of the evolutionary paths., Comment: 9 pages, 10 figures. Accepted for publication in ApJ

We report the results of optical--infrared follow-up observations of the gravitational wave (GW) event GW151226 detected by the Advanced LIGO in the framework of J-GEM (Japanese collaboration for Gravitational wave ElectroMagnetic follow-up). We performed wide-field optical imaging surveys with Kiso Wide Field Camera (KWFC), Hyper Suprime-Cam (HSC), and MOA-cam3. The KWFC survey started at 2.26 days after the GW event and covered 778 deg$^2$ centered at the high Galactic region of the skymap of GW151226. We started the HSC follow-up observations from 12 days after the event and covered an area of 63.5 deg$^2$ of the highest probability region of the northern sky with the limiting magnitudes of 24.6 and 23.8 for i band and z band, respectively. MOA-cam3 covered 145 deg$^2$ of the skymap with MOA-red filter 2.5 months after the GW alert. Total area covered by the wide-field surveys was 986.5 deg$^2$. The integrated detection probability of all the observed area was $\sim$29%. We also performed galaxy-targeted observations with six optical and near-infrared telescopes from 1.61 days after the event. Total of 238 nearby (<100 Mpc) galaxies were observed with the typical I band limiting magnitude of $\sim$19.5. We detected 13 supernova candidates with the KWFC survey, and 60 extragalactic transients with the HSC survey. Two third of the HSC transients were likely supernovae and the remaining one third were possible active galactic nuclei. With our observational campaign, we found no transients that are likely to be associated with GW151226., Comment: 12 pages, 6 figures, 4 tables, accepted for publication in PASJ

We present extensively dense observations of the super-Chandrasekhar supernova (SC SN) candidate SN 2012dn from $-11$ to $+140$ days after the date of its $B$-band maximum in the optical and near-infrared (NIR) wavelengths conducted through the OISTER ToO program. The NIR light curves and color evolutions up to 35 days after the $B$-band maximum provided an excellent match with those of another SC SN 2009dc, providing a further support to the nature of SN 2012dn as a SC SN. We found that SN 2012dn exhibited strong excesses in the NIR wavelengths from $30$ days after the $B$-band maximum. The $H$ and $K_{s}$-band light curves exhibited much later maximum dates at $40$ and $70$ days after the $B$-band maximum, respectively, compared with those of normal SNe Ia. The $H$ and $K_{s}$-band light curves subtracted by those of SN 2009dc displayed plateaued evolutions, indicating a NIR echo from the surrounding dust. The distance to the inner boundary of the dust shell is limited to be $4.8 - 6.4\times10^{-2}$ pc. No emission lines were found in its early phase spectrum, suggesting that the ejecta-CSM interaction could not occur. On the other hand, we found no signature that strongly supports the scenario of dust formation. The mass loss rate of the pre-explosion system is estimated to be $10^{-6}-10^{-5}$ M$_{\odot}$ yr$^{-1}$, assuming that the wind velocity of the system is $10-100$ km~s$^{-1}$, which suggests that the progenitor of SN 2012dn could be a recurrent nova system. We conclude that the progenitor of this SC SN could be explained by the single-degenerate scenario., Comment: 22 pages, 10 pages, 8 tables, accepted for publication in PASJ. Reference corrected