Search

Early galaxy formation, initiated by the dark matter and gas assembly, evolves through frequent mergers and feedback processes into dynamically hot, chaotic structures. In contrast, dynamically cold, smooth rotating disks have been observed in massive evolved galaxies merely 1.4 billion years after the Big Bang, suggesting rapid morphological and dynamical evolution in the early Universe. Probing this evolution mechanism necessitates studies of young galaxies, yet efforts have been hindered by observational limitations in both sensitivity and spatial resolution. Here we report high-resolution observations of a strongly lensed and quintuply imaged, low-luminosity, young galaxy at $z=6.072$ (dubbed the Cosmic Grapes), 930 million years after the Big Bang. Magnified by gravitational lensing, the galaxy is resolved into at least 15 individual star-forming clumps with effective radii of $r_{\rm e}\simeq$ 10--60 parsec (pc), which dominate $\simeq$ 70\% of the galaxy's total flux. The cool gas emission unveils a smooth, underlying rotating disk characterized by a high rotational-to-random motion ratio and a gravitationally unstable state (Toomre $Q \simeq$ 0.2--0.3), with high surface gas densities comparable to local dusty starbursts with $\simeq10^{3-5}$ $M_{\odot}$/pc$^{2}$. These gas properties suggest that the numerous star-forming clumps are formed through disk instabilities with weak feedback effects. The clumpiness of the Cosmic Grapes significantly exceeds that of galaxies at later epochs and the predictions from current simulations for early galaxies. Our findings shed new light on internal galaxy substructures and their relation to the underlying dynamics and feedback mechanisms at play during their early formation phases, potentially explaining the high abundance of bright galaxies observed in the early Universe and the dark matter core-cusp problem., Comment: Submitted. 44 pages, 14 figures, 5 tables. Comments are welcome! See also the companion papers on arXiv. Valentino+2024: arXiv:2402.17845 Gim\'enez-Arteaga+2024: arXiv:2402.17875

The ALMA lensing cluster survey (ALCS) is a 96-hr large program dedicated to uncovering and characterizing intrinsically faint continuum sources and line emitters with the assistance of gravitational lensing. All 33 cluster fields were selected from HST/Spitzer treasury programs including CLASH, Hubble Frontier Fields, and RELICS, which also have Herschel and Chandra coverages. The total sky area surveyed reaches $\sim$133 arcmin$^2$ down to a depth of $\sim$60 $\mu$Jy beam$^{-1}$ (1$\sigma$) at 1.2 mm, yielding 141 secure blind detections of continuum sources and additional 39 sources aided by priors. We present scientific motivation, survey design, the status of spectroscopy follow-up observations, and number counts down to $\sim$7 $\mu$Jy. Synergies with JWST are also discussed., Comment: 6 pages, 4 figures, Proceedings of the 7th Chile-Cologne-Bonn-Symposium: Physics and Chemistry of Star Formation, V. Ossenkopf-Okada, R. Schaaf, I. Breloy (eds.)

We present 850um imaging of the XMM-LSS field observed for 170 hours as part of the James Clerk Maxwell Telescope SCUBA-2 Large eXtragalactic Survey (S2LXS). S2LXS XMM-LSS maps an area of 9 square degrees, reaching a moderate depth of 1-sigma ~ 4 mJy/beam. This is the largest contiguous area of extragalactic sky mapped by JCMT at 850um to date. The wide area of the S2LXS XMM-LSS survey allows us to probe the ultra-bright (S_850um > 15 mJy), yet rare submillimetre population. We present the S2LXS XMM-LSS catalogue, which comprises 40 sources detected at >5-sigma significance, with deboosted flux densities in the range of 7 mJy to 48 mJy. We robustly measure the bright-end of the 850um number counts at flux densities >7 mJy, reducing the Poisson errors compared to existing measurements. The S2LXS XMM-LSS observed number counts show the characteristic upturn at bright fluxes, expected to be motivated by local sources of submillimetre emission and high-redshift strongly lensed galaxies. We find that the observed 850um number counts are best reproduced by model predictions that include either strong lensing or source blending from a 15 arcsec beam, indicating that both may make an important contribution to the observed over-abundance of bright single dish 850um selected sources. We make the S2LXS XMM-LSS 850um map and >5-sigma catalogue presented here publicly available., Comment: Accepted for publication in MNRAS

We present a joint analysis of the galaxy S04590 at $z=8.496$ based on NIRSpec, NIRCam, and NIRISS observations obtained through as part of Early Release Observations programme of the James Webb Space Telescope (JWST) and the far-infrared [CII]-$158\mu$m emission line detected by dedicated Atacama Large Millimeter/submillimeter Array (ALMA) observations. We determine the physical properties of S04590 from modelling of the spectral energy distribution (SED) and through the redshifted optical nebular emission lines detected with JWST/NIRSpec. The best-fit SED model reveals a low-mass ($M_\star = 10^{7.2}-10^{8}\,M_{\odot}$) galaxy with a low oxygen abundance of $12+\log{\rm (O/H)} = 7.16^{+0.10}_{-0.12}$ derived from the strong nebular and auroral emission lines. Assuming that [CII] effectively traces the interstellar medium (ISM), we estimate the total gas mass of the galaxy to be $M_{\rm gas} = (8.0\pm 4.0)\times 10^{8}\,M_\odot$ based on the luminosity and spatial extent of [CII]. This yields an exceptionally high gas fraction, $f_{\rm gas} = M_{\rm gas}/(M_{\rm gas} + M_\star) \gtrsim 90\%$, though still consistent within the range expected for its low metallicity. We further derive the metal mass of the galaxy based on the gas mass and gas-phase metallicity, which we find to be consistent with the expected metal production from Type II supernovae. Finally, we make the first constraints on the dust-to-gas (DTG) and dust-to-metals (DTM) ratios of galaxies in the epoch of reionization at $z\gtrsim 6$, showing overall low mass ratios of logDGT $<-3.8$ and logDTM $<-0.5$, though consistent with local scaling relations and in particular the local metal-poor galaxy I Zwicky 18. Our analysis highlights the synergy between ALMA and JWST in characterizing the gas, metal, and stellar content of the first generation of galaxies., Comment: Accepted for publication in ApJL

We introduce an ALMA band 3 spectroscopic survey, targeting the brightest submillimeter galaxies (SMGs) in the COSMOS field. Here we present the first results based on the 18 primary SMGs that have 870 $\mu$m flux densities of $S_{870}=12.4-19.3$ mJy and are drawn from a parent sample of 260 ALMA-detected SMGs from the AS2COSMOS survey. We detect emission lines in 17 and determine their redshifts to be in the range of $z=2-5$ with a median of ${3.3\pm0.3}$. We confirm that SMGs with brighter $S_{870}$ are located at higher redshifts. The data additionally cover five fainter companion SMGs, and we obtain line detection in one. Together with previous studies, our results indicate that for SMGs that satisfy our selection, their brightest companion SMGs are physically associated with their corresponding primary SMGs in $\ge40$% of the time, suggesting that mergers play a role in the triggering of star formation. By modeling the foreground gravitational fields, $<10$% of the primary SMGs can be strongly lensed with a magnification $\mu>2$. We determine that about 90\% of the primary SMGs have lines that are better described by double Gaussian profiles, and the median separation of the two Gaussian peaks is 430$\pm$40 km s$^{-1}$. This allows estimates of an average baryon mass, which together with the line dispersion measurements puts our primary SMGs on the similar mass-$\sigma$ correlation found on local early-type galaxies. Finally, the number density of our $z>4$ primary SMGs is found to be $1^{+0.9}_{-0.6}\times10^6$ cMpc$^{-3}$, suggesting that they can be the progenitors of $z\sim3-4$ massive quiescent galaxies., Comment: 30 pages, 16 figures, ApJ in press

We present a detailed investigation of millimeter-wave line emitters ALMA J010748.3-173028 (ALMA-J0107a) and ALMA J010747.0-173010 (ALMA-J0107b), which were serendipitously uncovered in the background of the nearby galaxy VV114 with spectral scan observations at $\lambda$ = 2 - 3 mm. Via Atacama Large Millimeter/submillimeter Array (ALMA) detection of CO(4-3), CO(3-2), and [CI](1-0) lines for both sources, their spectroscopic redshifts are unambiguously determined to be $z= 2.4666\pm0.0002$ and $z=2.3100\pm0.0002$, respectively. We obtain the apparent molecular gas masses $M_{\rm gas}$ of these two line emitters from [CI] line fluxes as $(11.2 \pm 3.1) \times 10^{10} M_\odot$ and $(4.2 \pm 1.2) \times 10^{10} M_\odot$, respectively. The observed CO(4-3) velocity field of ALMA-J0107a exhibits a clear velocity gradient across the CO disk, and we find that ALMA-J0107a is characterized by an inclined rotating disk with a significant turbulence, that is, a deprojected maximum rotation velocity to velocity dispersion ratio $v_{\rm max}/\sigma_{v}$ of $1.3 \pm 0.3$. We find that the dynamical mass of ALMA-J0107a within the CO-emitting disk computed from the derived kinetic parameters, $(1.1 \pm 0.2) \times 10^{10}\ M_\odot$, is an order of magnitude smaller than the molecular gas mass derived from dust continuum emission, $(3.2\pm1.6)\times10^{11}\ M_{\odot}$. We suggest this source is magnified by a gravitational lens with a magnification of $\mu \gtrsim10$, which is consistent with the measured offset from the empirical correlation between CO-line luminosity and width., Comment: 19 pages, 11 figures, accepted for publication in ApJ

We present a conceptual study of a large format imaging spectrograph for the Large Submillimeter Telescope (LST) and the Atacama Large Aperture Submillimeter Telescope (AtLAST). Recent observations of high-redshift galaxies indicate the onset of earliest star formation just a few 100 million years after the Big Bang (i.e., z = 12--15), and LST/AtLAST will provide a unique pathway to uncover spectroscopically-identified first forming galaxies in the pre-reionization era, once it will be equipped with a large format imaging spectrograph. We propose a 3-band (200, 255, and 350 GHz), medium resolution (R = 2,000) imaging spectrograph with 1.5 M detectors in total based on the KATANA concept (Karatsu et al. 2019), which exploits technologies of the integrated superconducting spectrometer (ISS) and a large-format imaging array. A 1-deg2 drilling survey (3,500 hr) will capture a large number of [O III] 88 um (and [C II] 158 um) emitters at z = 8--9, and constrain [O III] luminosity functions at z > 12., Comment: Presented at SPIE Astronomical Telescopes + Instrumentation 2020. Full published paper and video available at http://dx.doi.org/10.1117/12.2561238 See also https://www.lstobservatory.org for the LST project

We report the discovery of an intrinsically faint, quintuply-imaged, dusty galaxy MACS0600-z6 at a redshift $z=$6.07 viewed through the cluster MACSJ0600.1-2008 ($z$=0.46). A $\simeq4\sigma$ dust detection is seen at 1.2mm as part of the ALMA Lensing Cluster Survey (ALCS), an on-going ALMA Large program, and the redshift is secured via [C II] 158 $\mu$m emission described in a companion paper. In addition, spectroscopic follow-up with GMOS/Gemini-North shows a break in the galaxy's spectrum, consistent with the Lyman break at that redshift. We use a detailed mass model of the cluster and infer a magnification $\mu\gtrsim$30 for the most magnified image of this galaxy, which provides an unprecedented opportunity to probe the physical properties of a sub-luminous galaxy at the end of cosmic reionisation. Based on the spectral energy distribution, we infer lensing-corrected stellar and dust masses of $\rm{2.9^{+11.5}_{-2.3}\times10^9}$ and $\rm{4.8^{+4.5}_{-3.4}\times10^6}$ $\rm{M_{\odot}}$ respectively, a star formation rate of $\rm{9.7^{+22.0}_{-6.6} M_{\odot} yr^{-1}}$, an intrinsic size of $\rm{0.54^{+0.26}_{-0.14}}$ kpc, and a luminosity-weighted age of 200$\pm$100 Myr. Strikingly, the dust production rate in this relatively young galaxy appears to be larger than that observed for equivalent, lower redshift sources. We discuss if this implies that early supernovae are more efficient dust producers and the consequences for using dust mass as a probe of earlier star formation., Comment: submitted to MNRAS, revised version after including referee's comments

We report the discovery of a Multi Unit Spectroscopic Explorer (MUSE) galaxy group at z=4.32 lensed by the massive galaxy cluster ACT-CL J0102-4915 (aka El Gordo) at z=0.87, associated with a 1.2 mm source which is at a 2.07+/-0.88 kpc projected distance from one of the group galaxies. Three images of the whole system appear in the image plane. The 1.2 mm source has been detected within the Atacama Large Millimetre/submillimetre Array (ALMA) Lensing Cluster Survey (ALCS). As this ALMA source is undetected at wavelengths lambda < 2 microns, its redshift cannot be independently determined, however, the three lensing components indicate that it belongs to the same galaxy group at z=4.32. The four members of the MUSE galaxy group have low to intermediate stellar masses (~ 10^7-10^{10} Msun) and star formation rates (SFRs) of 0.4-24 Msun/yr, resulting in high specific SFRs (sSFRs) for two of them, which suggest that these galaxies are growing fast (with stellar-mass doubling times of only ~ 2x10^7 years). This high incidence of starburst galaxies is likely a consequence of interactions within the galaxy group, which is compact and has high velocity dispersion. Based on the magnification-corrected sub-/millimetre continuum flux density and estimated stellar mass, we infer that the ALMA source is classified as an ordinary ultra-luminous infrared galaxy (with associated dust-obscured SFR~200-300 Msun/yr) and lies on the star-formation main sequence. This reported case of an ALMA/MUSE group association suggests that some presumably isolated ALMA sources are in fact signposts of richer star-forming environments at high redshifts., Comment: 13 pages including 7 figures and 1 table. Accepted for publication at the ApJ. Minor changes with respect to version 1. Figure 6 has been expanded to broaden comparison with the literature

We make use of the ASAGAO, deep 1.2 mm continuum observations of a 26 arcmin$^2$ region in the GOODS-South field obtained with ALMA, to probe dust-enshrouded star formation in $K$-band selected (i.e., stellar mass selected) galaxies, which are drawn from the ZFOURGE catalog. Based on the ASAGAO combined map, which was created by combining ASAGAO and ALMA archival data in the GOODS-South field, we find that 24 ZFOURGE sources have 1.2 mm counterparts with a signal-to-noise ratio $>$ 4.5 (1$\sigma\simeq$ 30 - 70 $\mu$Jy beam$^{-1}$ at 1.2 mm). Their median redshift is estimated to be $z_\mathrm{median}=$ 2.38 $\pm$ 0.14. They generally follow the tight relationship of the stellar mass versus star formation rate (i.e., the main sequence of star-forming galaxies). ALMA-detected ZFOURGE sources exhibit systematically larger infrared (IR) excess (IRX $\equiv L_\mathrm{IR}/L_\mathrm{UV}$) compared to ZFOURGE galaxies without ALMA detections even though they have similar redshifts, stellar masses, and star formation rates. This implies the consensus stellar-mass versus IRX relation, which is known to be tight among rest-frame-UV-selected galaxies, can not fully predict the ALMA detectability of stellar-mass-selected galaxies. We find that ALMA-detected ZFOURGE sources are the main contributors to the cosmic IR star formation rate density at $z$ = 2 - 3., Comment: Accepted for publication in PASJ. A version with a high resolution figure and ALMA fits files are available from https://sites.google.com/view/asagao26/

Cosmological simulations predict the Universe contains a network of intergalactic gas filaments, within which galaxies form and evolve. However, the faintness of any emission from these filaments has limited tests of this prediction. We report the detection of rest-frame ultraviolet Lyman-alpha radiation from multiple filaments extending more than one megaparsec between galaxies within the SSA 22 proto-cluster at a redshift of 3.1. Intense star formation and supermassive black-hole activity is occurring within the galaxies embedded in these structures, which are the likely sources of the elevated ionizing radiation powering the observed Lyman-alpha emission. Our observations map the gas in filamentary structures of the type thought to fuel the growth of galaxies and black holes in massive proto-clusters., Comment: Published in Science on 2019 October 3 (Main: 3 figures, Supp: 10 figures, 2 tables)

We present observations at 1.1 mm towards 16 powerful radio galaxies and a radio-quiet quasar at 0.5

We present the results of 3 GHz radio continuum observations of the 8 host galaxies of super-luminous supernovae (SLSNe) at $0.1 < z < 0.3$ by using the Karl G. Jansky Very Large Array. Four host galaxies are detected significantly, and two of them are found to have high star-formation rates (SFRs $>$ 20 $M_{\odot}$ yr$^{-1}$) derived from radio emission, making them the most intensely star-forming host galaxies among SLSN host galaxies. We compare radio SFRs and optical SFRs, and find that three host galaxies have an excess in radio SFRs by a factor of $>$2, suggesting the existence of dust-obscured star formation, which cannot be traced by optical studies. Two of the three host galaxies, which are located in the galaxy main sequence based on optical SFRs, are found to be above the main sequence based on their radio SFRs. This suggests a higher fraction of starburst galaxies in SLSN hosts than estimated in previous studies. We calculate extinction from the ratio between radio SFRs and dust-uncorrected optical SFRs and find that the hosts are on the trend of increasing extinction with metallicity, which is consistent with the relation in local star-forming galaxies. We also place a constraint on a pulsar-driven SN model, which predicts quasi-steady synchrotron radio emission., Comment: 7 pages, 8 figures, accepted for publication in ApJ

We investigate the X-ray active galactic nucleus (AGN) properties of millimeter galaxies in the Great Observatories Origins Deep Survey South (GOODS-S) field detected with the Atacama Large Millimeter/submillimeter Array (ALMA), by utilizing the Chandra 7-Ms data, the deepest X-ray survey to date. Our millimeter galaxy sample comes from the ASAGAO survey covering 26 arcmin$^2$ (12 sources at a 1.2-mm flux-density limit of $\approx$0.6 mJy), supplemented by the deeper but narrower 1.3-mm survey of a part of the ASAGAO field by Dunlop et al.\ (2017). Fourteen out of the total 25 millimeter galaxies have Chandra counterparts. The observed AGN fractions at $z=1.5-3$ is found to be 90$^{+8}_{-19}$\% and $57^{+23}_{-25}$\% for the ultra/luminous infrared galaxies with $\log L_{\rm IR}/L_{\odot} = 12-12.8$ and $\log L_{\rm IR}/L_{\odot} = 11.5-12$, respectively. The majority ($\sim$2/3) of the ALMA and/or Herschel detected X-ray AGNs at $z=1.5-3$ appear to be star-formation dominant populations, having $L_{\rm X}$/ $L_{\rm IR}$ ratios smaller than the "simultaneous evolution" value expected from the local black-hole mass to stellar mass ($M_{\rm BH}$-$M_*$) relation. On the basis of the $L_{\rm X}$ and stellar mass relation, we infer that a large fraction of star-forming galaxies at $z=1.5-3$ have black hole masses smaller than those expected from the local $M_{\rm BH}$-$M_*$ relation. This is opposite to previous reports on luminous AGN at same redshifts detected in wider and shallower surveys, which are subject to selection biases against lower luminosity AGN. Our results are consistent with an evolutionary scenario that star formation occurs first, and an AGN-dominant phase follows later, in objects finally evolving into galaxies with classical bulges., Comment: 14 pages including 5 figures and 1 table, accepted for publication in ApJ

We study the heating mechanisms and Ly{\alpha} escape fractions of 35 Ly{\alpha} blobs (LABs) at z = 3.1 in the SSA22 field. Dust continuum sources have been identified in 11 of the 35 LABs, all with star formation rates (SFRs) above 100 Msun/yr. Likely radio counterparts are detected in 9 out of 29 investigated LABs. The detection of submm dust emission is more linked to the physical size of the Ly{\alpha} emission than to the Ly{\alpha} luminosities of the LABs. A radio excess in the submm/radio detected LABs is common, hinting at the presence of active galactic nuclei. Most radio sources without X-ray counterparts are located at the centers of the LABs. However, all X-ray counterparts avoid the central regions. This may be explained by absorption due to exceptionally large column densities along the line-of-sight or by LAB morphologies, which are highly orientation dependent. The median Ly{\alpha} escape fraction is about 3\% among the submm-detected LABs, which is lower than a lower limit of 11\% for the submm-undetected LABs. We suspect that the large difference is due to the high dust attenuation supported by the large SFRs, the dense large-scale environment as well as large uncertainties in the extinction corrections required to apply when interpreting optical data., Comment: 22 pages, 7 figures, ready for publication on ApJ

We report the identification of four millimeter line emitting galaxies with the Atacama Large Milli/submillimeter Array (ALMA) in SSA22 Field (ADF22). We analyze the ALMA 1.1 mm survey data, with an effective survey area of 5 arcmin$^2$, a frequency range of 253.1--256.8 and 269.1--272.8 GHz, angular resolution of 0".7 and RMS noise of 0.8 mJy beam$^{-1}$ at 36 km s$^{-1}$ velocity resolution. We detect four line emitter candidates with significance levels above $6 \sigma$. We identify one of the four sources as a CO(9-8) emitter at $z = 3.1$ in a member of the proto-cluster known in this field. Another line emitter with an optical counterpart is likely a CO(4-3) emitter at $z = 0.7$. The other two sources without any millimeter continuum or optical/near-infrared counterpart are likely to be [CII] emitter candidates at $z = 6.0$ and $6.5$. The equivalent widths of the [CII] candidates are consistent with those of confirmed high-redshift [CII] emitters and candidates, and are a factor of 10 times larger than that of the CO(9-8) emitter detected in this search. The [CII] luminosity of the candidates are $4-7 \times 10^8~\rm L_\odot$. The star formation rates (SFRs) of these sources are estimated to be $10-20~\rm M_{\odot}~yr^{-1}$ if we adopt an empirical [CII] luminosity - SFR relation. One of them has a relatively low-S/N ratio, but shows features characteristic of emission lines. Assuming that at least one of the two candidates is a [CII] emitter, we derive a lower limit of [CII]-based star formation rate density (SFRD) at $z~\sim~6$. The resulting value of $> 10^{-2}$ $\rm M_\odot yr^{-1} Mpc^{-3}$ is consistent with the dust-uncorrected UV-based SFRD. Future millimeter/submillimeter surveys can be used to detect a number of high redshift line emitters, with which to study the star formation history in the early Universe., Comment: 14 pages, 10 figures. Accepted to PASJ

We present a Herschel/SPIRE survey of three protoclusters at z=2-3 (2QZCluster, HS1700, SSA22). Based on the SPIRE colours (S350/S250 and S500/S350) of 250 $\mu$m sources, we selected high redshift dusty star-forming galaxies potentially associated with the protoclusters. In the 2QZCluster field, we found a 4-sigma overdensity of six SPIRE sources around 4.5' (~2.2 Mpc) from a density peak of H$\alpha$ emitters at z=2.2. In the HS1700 field, we found a 5-sigma overdensity of eight SPIRE sources around 2.1' (~1.0 Mpc) from a density peak of LBGs at z=2.3. We did not find any significant overdensities in SSA22 field, but we found three 500 $\mu$m sources are concentrated 3' (~1.4 Mpc) east to the LAEs overdensity. If all the SPIRE sources in these three overdensities are associated with protoclusters, the inferred star-formation rate densities are 10$^3$-10$^4$ times higher than the average value at the same redshifts. This suggests that dusty star-formation activity could be very strongly enhanced in z~2-3 protoclusters. Further observations are needed to confirm the redshifts of the SPIRE sources and to investigate what processes enhance the dusty star-formation activity in z~2-3 protoclusters., Comment: 15 pages, 5 figures, 12 tables, Accepted for publication in MNRAS

We have conducted 1.1 mm ALMA observations of a contiguous $105'' \times 50''$ or 1.5 arcmin$^2$ window in the SXDF-UDS-CANDELS. We achieved a 5$\sigma$ sensitivity of 0.28 mJy, providing a flat sensus of dusty star-forming galaxies with $L_{\rm IR} \sim6\times10^{11}$ $L_\odot$ (for $T_{\rm dust}$ =40K) up to $z\sim10$ thanks to the negative K-correction at this wavelength. We detected 5 brightest sources (S/N$>$6) and 18 low-significance sources (5$>$S/N$>$4; these may contain spurious detections, though). One of the 5 brightest ALMA sources ($S_{\rm 1.1mm} = 0.84 \pm 0.09$ mJy) is extremely faint in the WFC3 and VLT/HAWK-I images, demonstrating that a contiguous ALMA imaging survey is able to uncover a faint dust-obscured population that is invisible in deep optical/near-infrared surveys. We found a possible [CII]-line emitter at $z=5.955$ or a low-$z$ CO emitting galaxy within the field, which may allow us to constrain the [CII] and/or the CO luminosity functions across the history of the universe., Comment: 4 pages, 2 figures, 1 table, to appear in the proceedings of IAU Symposium 319 "Galaxies at High Redshift and Their Evolution over Cosmic Time", eds. S. Kaviraj & H. Ferguson

We report the results of $1^{\prime}.5 \times3^{\prime}$ mapping at 1.1~mm with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the central region of the $z=3.09$ SSA22 protocluster. By combining our source catalog with archival spectroscopic redshifts, we find that eight submillimeter galaxies (SMGs) with flux densities, $S_{\rm 1.1~mm}=0.7-6.4$~mJy ($L_{\rm IR}\sim10^{12.1}-10^{13.1}L_\odot$) are at $z=3.08-3.10$. Not only are these SMGs members of the protocluster but they in fact reside within the node at the junction of the 50 Mpc-scale filamentary three-dimensional structure traced by Lyman-$\alpha$ emitters (LAEs) in this field. The eight SMGs account for a star formation rate density (SFRD) $\sim$10 $M_\odot$ yr$^{-1}$ Mpc$^{-3}$ in the node, which is two orders of magnitudes higher than the global SFRD at this redshift. We find that four of the eight SMGs host a X-ray luminous active galactic nuclei (AGN). Our results suggest that the vigorous star formation activity and the growth of super massive black holes (SMBHs) occurred simultaneously in the densest regions at $z\sim3$, which may correspond to the most active historical phase of the massive galaxy population found in the core of the clusters in the present universe. Two SMGs are associated with Lyman-$\alpha$ blobs (LABs), implying that the two populations coexist in high density environments for a few cases., Comment: Accepted for Publication in ApJL. 6 pages, 3 figures

We report detections of two candidate distant submillimeter galaxies (SMGs), MM J154506.4$-$344318 and MM J154132.7$-$350320, which are discovered in the AzTEC/ASTE 1.1 mm survey toward the Lupus-I star-forming region. The two objects have 1.1 mm flux densities of 43.9 and 27.1 mJy, and have Herschel/SPIRE counterparts as well. The Submillimeter Array counterpart to the former SMG is identified at 890 $\mu$m and 1.3 mm. Photometric redshift estimates using all available data from the mid-infrared to the radio suggest that the redshifts of the two SMGs are $z_{\rm photo} \simeq$ 4-5 and 3, respectively. Near-infrared objects are found very close to the SMGs and they are consistent with low-$z$ ellipticals, suggesting that the high apparent luminosities can be attributed to gravitational magnification. The cumulative number counts at $S_{\rm 1.1mm} \ge 25$ mJy, combined with other two 1.1-mm brightest sources, are $0.70 ^{+0.56}_{-0.34}$ deg$^{-2}$, which is consistent with a model prediction that accounts for flux magnification due to strong gravitational lensing. Unexpectedly, a $z > 3$ SMG and a Galactic dense starless core (e.g., a first hydrostatic core) could be similar in the mid-infrared to millimeter spectral energy distributions and spatial structures at least at $\gtrsim 1"$. This indicates that it is necessary to distinguish the two possibilities by means of broad band photometry from the optical to centimeter and spectroscopy to determine the redshift, when a compact object is identified toward Galactic star-forming regions., Comment: Accepted for publication in ApJ. 16 pages, 9 figures, 2 tables

We present a result of a blind search for [CII] 158 $\mu$m emitters at $z\sim 4.5$ using ALMA Cycle~0 archival data. We collected extra-galactic data covering at 330-360 GHz (band~7) from 8 Cycle~0 projects from which initial results have been already published. The total number of fields is 243 and the total on-source exposure time is 19.2 hours. We searched for line emitters in continuum-subtracted data cubes with spectral resolutions of 50, 100, 300 and 500 km/s. We could not detect any new line emitters above a 6-$\sigma$ significance level. This result provides upper limits to the [CII] luminosity function at $z\sim 4.5$ over $L_{\rm [CII]} \sim 10^8 - 10^{10} L_{\odot}$ or star formation rate, SFR $\sim$ 10-1000 M$_{^\odot}$/yr. These limits are at least 2 orders of magnitude larger than the [CII] luminosity functions expected from the $z \sim 4$ UV luminosity function or from numerical simulation. However, this study demonstrates that we would be able to better constrain the [CII] luminosity function and to investigate possible contributions from dusty galaxies to the cosmic star-formation rate density by collecting Cycle~1+2 archival data as the ALMA Patchy Deep Survey., Comment: 6 pages, 3 figures, accepted for publication in MNRAS

We present initial results of very high resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of the $z$=3.042 gravitationally lensed galaxy HATLAS J090311.6+003906 (SDP.81). These observations were carried out using a very extended configuration as part of Science Verification for the 2014 ALMA Long Baseline Campaign, with baselines of up to 15 km. We present continuum imaging at 151, 236 and 290 GHz, at unprecedented angular resolutions as fine as 23 milliarcseconds (mas), corresponding to an un-magnified spatial scale of ~180 pc at z=3.042. The ALMA images clearly show two main gravitational arc components of an Einstein ring, with emission tracing a radius of ~1.5". We also present imaging of CO(10-9), CO(8-7), CO(5-4) and H2O line emission. The CO emission, at an angular resolution of ~170 mas, is found to broadly trace the gravitational arc structures but with differing morphologies between the CO transitions and compared to the dust continuum. Our detection of H2O line emission, using only the shortest baselines, provides the most resolved detection to date of thermal H2O emission in an extragalactic source. The ALMA continuum and spectral line fluxes are consistent with previous Plateau de Bure Interferometer and Submillimeter Array observations despite the impressive increase in angular resolution. Finally, we detect weak unresolved continuum emission from a position that is spatially coincident with the center of the lens, with a spectral index that is consistent with emission from the core of the foreground lensing galaxy., Comment: 9 pages, 5 figures and 3 tables, accepted for publication in the Astrophysical Journal Letters

Long-duration gamma-ray bursts are associated with the explosions of massive stars and are accordingly expected to reside in star-forming regions with molecular gas (the fuel for star formation). Previous searches for carbon monoxide (CO), a tracer of molecular gas, in burst host galaxies did not detect any emission. Molecules have been detected as absorption in the spectra of gamma-ray burst afterglows, and the molecular gas is similar to the translucent or diffuse molecular clouds of the Milky Way. Absorption lines probe the interstellar medium only along the line of sight, so it is not clear whether the molecular gas represents the general properties of the regions where the bursts occur. Here we report spatially resolved observations of CO line emission and millimetre-wavelength continuum emission in two galaxies hosting gamma-ray bursts. The bursts happened in regions rich in dust, but not particularly rich in molecular gas. The ratio of molecular gas to dust (<9-14) is significantly lower than in star-forming regions of the Milky Way and nearby star-forming galaxies, suggesting that much of the dense gas where stars form has been dissipated by other massive stars., Comment: To appear in Nature June 12, 2014. 17 pages, 4 figures and 1 table, including Methods section, Extended Data Figures, and Supplementary Information

We present the results from a 1.1 mm imaging survey of the SSA22 field, known for having an overdensity of z=3.1 Lyman-alpha emitting galaxies (LAEs), taken with the AzTEC camera on the Atacama Submillimeter Telescope Experiment (ASTE). We imaged a 950 arcmin$^2$ field down to a 1 sigma sensitivity of 0.7-1.3 mJy/beam to find 125 submillimeter galaxies (SMGs) with a signal to noise ratio >= 3.5. Counterpart identification using radio and near/mid-infrared data was performed and one or more counterpart candidates were found for 59 SMGs. Photometric redshifts based on optical to near-infrared images were evaluated for 45 SMGs of these SMGs with Spitzer/IRAC data, and the median value is found to be z=2.4. By combining these estimation with estimates from the literature we determined that 10 SMGs might lie within the large-scale structure at z=3.1. The two-point angular cross-correlation function between LAEs and SMGs indicates that the positions of the SMGs are correlated with the z=3.1 protocluster. These results suggest that the SMGs were formed and evolved selectively in the high dense environment of the high redshift universe. This picture is consistent with the predictions of the standard model of hierarchical structure formation., Comment: 17 pages, 9 figures, 7 tables, Accepted for publication in MNRAS

We present initial results from the CO survey toward high redshift galaxies using the Nobeyama 45m telescope. Using the new wide bandwidth spectrometer equipped with a two-beam SIS receiver, we have robust new detections of three high redshift (z=1.6-3.4) submillimeter galaxies (SXDF 1100.001, SDP9, and SDP17), one tentative detection (SDSS J160705+533558), and one non-detection (COSMOS-AzTEC1). The galaxies observed during the commissioning phase are sources with known spectroscopic redshifts from previous optical or from wide-band submm spectroscopy. The derived molecular gas mass and line widths from Gaussian fits are ~10^11 Msun and 430-530 km/s, which are consistent with previous CO observations of distant submm galaxies and quasars. The spectrometer that allows a maximum of 32 GHz instantaneous bandwidth will provide new science capabilities at the Nobeyama 45m telescope, allowing us to determine redshifts of bright submm selected galaxies without any prior redshift information., Comment: 4 pages, 1 figure, PASJ Letter Accepted

The source counts of galaxies discovered at sub-millimetre and millimetre wavelengths provide important information on the evolution of infrared-bright galaxies. We combine the data from six blank-field surveys carried out at 1.1 mm with AzTEC, totalling 1.6 square degrees in area with root-mean-square depths ranging from 0.4 to 1.7 mJy, and derive the strongest constraints to date on the 1.1 mm source counts at flux densities S(1100) = 1-12 mJy. Using additional data from the AzTEC Cluster Environment Survey to extend the counts to S(1100) ~ 20 mJy, we see tentative evidence for an enhancement relative to the exponential drop in the counts at S(1100) ~ 13 mJy and a smooth connection to the bright source counts at >20 mJy measured by the South Pole Telescope; this excess may be due to strong lensing effects. We compare these counts to predictions from several semi-analytical and phenomenological models and find that for most the agreement is quite good at flux densities > 4 mJy; however, we find significant discrepancies (>3sigma) between the models and the observed 1.1 mm counts at lower flux densities, and none of them are consistent with the observed turnover in the Euclidean-normalised counts at S(1100) < 2 mJy. Our new results therefore may require modifications to existing evolutionary models for low luminosity galaxies. Alternatively, the discrepancy between the measured counts at the faint end and predictions from phenomenological models could arise from limited knowledge of the spectral energy distributions of faint galaxies in the local Universe., Comment: 16 pages, 3 figures, 4 tables; accepted for publication in MNRAS

We report the results of the counterpart identification and a detailed analysis of the physical properties of the 48 sources discovered in our deep 1.1mm wavelength imaging survey of the GOODS-South field using the AzTEC instrument on the Atacama Submillimeter Telescope Experiment (ASTE). One or more robust or tentative counterpart candidate is found for 27 and 14 AzTEC sources, respectively, by employing deep radio continuum, Spitzer MIPS & IRAC, and LABOCA 870 micron data. Five of the sources (10%) have two robust counterparts each, supporting the idea that these galaxies are strongly clustered and/or heavily confused. Photometric redshifts and star formation rates (SFRs) are derived by analyzing UV-to-optical and IR-to-radio SEDs. The median redshift of z~2.6 is similar to other earlier estimates, but we show that 80% of the AzTEC-GOODS sources are at z>2, with a significant high redshift tail (20% at z>3.3). Rest-frame UV and optical properties of AzTEC sources are extremely diverse, spanning 10 magnitude in the i- and K-band photometry with median values of i=25.3 and K=22.6 and a broad range of red colour (i-K=0-6). These AzTEC sources are some of the most luminous galaxies in the rest-frame optical bands at z>2, with inferred stellar masses of (1-30) x 10^{10} solar masses and UV-derived star formation rates of SFR(UV) > 10-1000 solar masses per year. The IR-derived SFR, 200-2000 solar masses per year, is independent of redshift or stellar mass. The resulting specific star formation rates, SSFR = 1-100 per Gyr, are 10-100 times higher than similar mass galaxies at z=0, and they extend the previously observed rapid rise in the SSFR with redshift to z=2-5. These galaxies have a SFR high enough to have built up their entire stellar mass within their Hubble time. We find only marginal evidence for an AGN contribution to the near-IR and mid-IR SEDs. (abridged), Comment: 31 pages including 14 figures, accepted for publication in the MNRAS. A higher quality Figure 1 is also included as Figure1.jpg

We present a 0.72 sq. deg. contiguous 1.1mm survey in the central area of the COSMOS field carried out to a 1sigma ~ 1.26 mJy/beam depth with the AzTEC camera mounted on the 10m Atacama Submillimeter Telescope Experiment (ASTE). We have uncovered 189 candidate sources at a signal-to-noise ratio S/N >= 3.5, out of which 129, with S/N >= 4, can be considered to have little chance of being spurious (< 2 per cent). We present the number counts derived with this survey, which show a significant excess of sources when compared to the number counts derived from the ~0.5 sq. deg. area sampled at similar depths in the Scuba HAlf Degree Extragalactic Survey (SHADES, Austermann et al. 2010). They are, however, consistent with those derived from fields that were considered too small to characterize the overall blank-field population. We identify differences to be more significant in the S > 5 mJy regime, and demonstrate that these excesses in number counts are related to the areas where galaxies at redshifts z < 1.1 are more densely clustered. The positions of optical-IR galaxies in the redshift interval 0.6 < z < 0.75 are the ones that show the strongest correlation with the positions of the 1.1mm bright population (S > 5 mJy), a result which does not depend exclusively on the presence of rich clusters within the survey sampled area. The most likely explanation for the observed excess in number counts at 1.1mm is galaxy-galaxy and galaxy-group lensing at moderate amplification levels, that increases in amplitude as one samples larger and larger flux densities. This effect should also be detectable in other high redshift populations., Comment: 21 pages, 17 figures, accepted for publication in MNRAS

We present results from Submillimeter Array (SMA) 860-micron sub-arcsec astrometry and multiwavelength observations of the brightest millimeter (S_1.1mm = 8.4 mJy) source, SSA22-AzTEC1, found near the core of the SSA22 protocluster that is traced by Ly\alpha emitting galaxies at z = 3.09. We identify a 860-micron counterpart with a flux density of S_860um = 12.2 +/- 2.3 mJy and absolute positional accuracy that is better than 0.3". At the SMA position, we find radio to mid-infrared counterparts, whilst no object is found in Subaru optical and near-infrared deep images at wavelengths \le 1 micron (J > 25.4 in AB, 2\sigma). The photometric redshift estimate, using flux densities at \ge 24 microns, indicates z_phot = 3.19^{+0.26}_{-0.35}, consistent with the protocluster redshift. We then model the near-to-mid-infrared spectral energy distribution (SED) of SSA22-AzTEC1, and find that the SED modeling requires a large extinction (A_V \approx 3.4 mag) of starlight from a stellar component with M_star ~ 10^{10.9} M_sun, assuming z = 3.1. Additionally, we find a significant X-ray counterpart with a very hard spectrum (Gamma_eff = -0.34 ^{+0.57}_{-0.61}), strongly suggesting that SSA22-AzTEC1 harbors a luminous AGN (L_X ~ 3*10^{44} ergs s^{-1}) behind a large hydrogen column (N_H ~ 10^{24} cm^{-2}). The AGN, however, is responsible for only ~10% of the bolometric luminosity of the host galaxy, and therefore the star-formation activity likely dominates the submillimeter emission. It is possible that SSA22-AzTEC1 is the first example of a protoquasar growing at the bottom of the gravitational potential underlying the SSA22 protocluster., Comment: 16 pages, 2 tables, 6 figures, accepted for publication in the Astrophysical Journal

We present results of a 1.1 mm deep survey of the AKARI Deep Field South (ADF-S) with AzTEC mounted on the Atacama Submillimetre Telescope Experiment (ASTE). We obtained a map of 0.25 sq. deg area with an rms noise level of 0.32-0.71 mJy. This is one of the deepest and widest maps thus far at millimetre and submillimetre wavelengths. We uncovered 198 sources with a significance of 3.5-15.6 sigma, providing the largest catalog of 1.1 mm sources in a contiguous region. Most of the sources are not detected in the far-infrared bands of the AKARI satellite, suggesting that they are mostly at z ~ 1.5 given the detection limits. We constructed differential and cumulative number counts in the ADF-S, the Subaru/XMM Newton Deep Field (SXDF), and the SSA 22 field surveyed by AzTEC/ASTE, which provide currently the tightest constraints on the faint end. The integration of the best-fit number counts in the ADF-S find that the contribution of 1.1 mm sources with fluxes >=1 mJy to the cosmic infrared background (CIB) at 1.1 mm is 12-16%, suggesting that the large fraction of the CIB originates from faint sources of which the number counts are not yet constrained. We estimate the cosmic star-formation rate density contributed by 1.1 mm sources with >=1 mJy using the best-fit number counts in the ADF-S and find that it is lower by about a factor of 5-10 compared to those derived from UV/optically-selected galaxies at z ~ 2-3. The fraction of stellar mass of the present-day universe produced by 1.1 mm sources with >=1 mJy at z >= 1 is ~20%, calculated by the time integration of the star-formation rate density. If we consider the recycled fraction of >0.4, which is the fraction of materials forming stars returned to the interstellar medium, the fraction of stellar mass produced by 1.1 mm sources decrease to <~10%., Comment: 15 pages, 12 figure, accepted for publication in MNRAS

We report the detection of an extremely bright ($\sim$37 mJy at 1100 $\mu$m and $\sim$91 mJy at 880 $\mu$m) submillimeter galaxy (SMG), AzTEC-ASTE-SXDF1100.001 (hereafter referred to as SXDF1100.001 or Orochi), discovered in 1100 $\mu$m observations of the Subaru/XMM-Newton Deep Field using AzTEC on ASTE. Subsequent CARMA 1300 $\mu$m and SMA 880 $\mu$m observations successfully pinpoint the location of Orochi and suggest that it has two components, one extended (FWHM of $\sim$ 4$^{\prime\prime}$) and one compact (unresolved). Z-Spec on CSO has also been used to obtain a wide band spectrum from 190 to 308 GHz, although no significant emission/absorption lines are found. The derived upper limit to the line-to-continuum flux ratio is 0.1--0.3 (2 $\sigma$) across the Z-Spec band. Based on the analysis of the derived spectral energy distribution from optical to radio wavelengths of possible counterparts near the SMA/CARMA peak position, we suggest that Orochi is a lensed, optically dark SMG lying at $z \sim 3.4$ behind a foreground, optically visible (but red) galaxy at $z \sim 1.4$. The deduced apparent (i.e., no correction for magnification) infrared luminosity ($L_{\rm IR}$) and star formation rate (SFR) are $6 \times 10^{13}$ $L_{\odot}$ and 11000 $M_{\odot}$ yr$^{-1}$, respectively, assuming that the $L_{\rm IR}$ is dominated by star formation. These values suggest that Orochi will consume its gas reservoir within a short time scale ($3 \times 10^{7}$ yr), which is indeed comparable to those in extreme starbursts like the centres of local ULIRGs., Comment: 18 pages, 13 figures

[Abridged] We present the first results from a 1.1 mm confusion-limited map of the GOODS-S field taken with AzTEC on the ASTE telescope. We imaged a 270 sq. arcmin field to a 1\sigma depth of 0.48 - 0.73 mJy/beam, making this one of the deepest blank-field surveys at mm-wavelengths ever achieved. Although our GOODS-S map is extremely confused, we demonstrate that our source identification and number counts analyses are robust, and the techniques discussed in this paper are relevant for other deeply confused surveys. We find a total of 41 dusty starburst galaxies with S/N >= 3.5 within this uniformly covered region, where only two are expected to be false detections. We derive the 1.1mm number counts from this field using both a "P(d)" analysis and a semi-Bayesian technique, and find that both methods give consistent results. Our data are well-fit by a Schechter function model with (S', N(3mJy), \alpha) = (1.30+0.19 mJy, 160+27 (mJy/deg^2)^(-1), -2.0). Given the depth of this survey, we put the first tight constraints on the 1.1 mm number counts at S(1.1mm) = 0.5 mJy, and we find evidence that the faint-end of the number counts at S(850\mu m) < 2.0 mJy from various SCUBA surveys towards lensing clusters are biased high. In contrast to the 870 \mu m survey of this field with the LABOCA camera, we find no apparent under-density of sources compared to previous surveys at 1.1 mm. Additionally, we find a significant number of SMGs not identified in the LABOCA catalogue. We find that in contrast to observations at wavelengths < 500 \mu m, MIPS 24 \mu m sources do not resolve the total energy density in the cosmic infrared background at 1.1 mm, demonstrating that a population of z > 3 dust-obscured galaxies that are unaccounted for at these shorter wavelengths potentially contribute to a large fraction (~2/3) of the infrared background at 1.1 mm., Comment: 21 pages, 9 figures. Accepted to MNRAS.

We present an 880 micron Submillimeter Array (SMA) detection of the submillimeter galaxy SXDF850.6. SXDF850.6 is a bright source (S(850 micron) = 8 mJy) detected in the SCUBA Half Degree Extragalactic Survey (SHADES), and has multiple possible radio counterparts in its deep radio image obtained at the VLA. Our new SMA detection finds that the submm emission coincides with the brightest radio emission that is found ~8" north of the coordinates determined from SCUBA. Despite the lack of detectable counterparts in deep UV/optical images, we find a source at the SMA position in near-infrared and longer wavelength images. We perform SED model fits to UV-optical-IR photometry (u, B, V, R, i', z', J, H, K, 3.6 micron, 4.5 micron, 5.8 micron, and 8.0 micron) and to submm-radio photometry (850 micron, 880 micron, 1100 micron, and 21 cm) independently, and we find both are well described by starburst templates at a redshift of z ~= 2.2 (+/- 0.3). The best-fit parameters from the UV-optical-IR SED fit are a redshift of z = 1.87 (+0.15/-0.07), a stellar mass of M_star = 2.5 +2.2/-0.3 x 10^11 M_sun, an extinction of A_V = 3.0 (+0.3/-1.0) mag, and an age of 720 (+1880/-210) Myr. The submm-radio SED fit provides a consistent redshift of z ~ 1.8-2.5, an IR luminosity of L_IR = (7-26) x 10^12 L_sun, and a star formation rate of 1300-4500 M_sun/yr. These results suggest that SXDF850.6 is a mature system already having a massive amount of old stellar population constructed before its submm bright phase and is experiencing a dusty starburst, possibly induced by major mergers., Comment: 7 pages, 5 figures, Accepted for publication in Astrophysical Journal

We report our recent progress on extragalactic spectroscopic and continuum observations, including HCN(J=1-0), HCO$^+$(J=1-0), and CN(N=1-0) imaging surveys of local Seyfert and starburst galaxies using the Nobeyama Millimeter Array, high-J CO observations (J=3-2 observations using the Atacama Submillimeter Telescope Experiment (ASTE) and J=2-1 observations with the Submillimeter Array) of galaxies, and $\lambda$ 1.1 mm continuum observations of high-z violent starburst galaxies using the bolometer camera AzTEC mounted on ASTE., Comment: 6 pages, 5 figures, To appear in proceedings of "Far-Infrared and Submillimeter Emission of the Interstellar Medium", EAS Publication Series, Bad Honnef, November 2007, Eds. C. Kramer, S. Aalto, R. Simon. See http://www.nro.nao.ac.jp/~f0212kk/FIR07/kk-ver20.pdf for a version with high resolution figures

We present CO(3-2) emission observations toward the 3'x3' (or 20x20kpc at a distance of 23Mpc) region of the southern barred spiral galaxy NGC 986 using the Atacama Submillimeter Telescope Experiment (ASTE). This effort is a part of our on-going extragalactic CO(3-2) imaging project ADIoS (ASTE Dense gas Imaging of Spiral galaxies). Our CO(3-2) image revealed the presence of a large (the major axis is 14 kpc in total length) gaseous bar filled with dense molecular medium along the dark lanes observed in optical images. This is the largest ``dense-gas rich bar'' known to date. The dense gas bar discovered in NGC 986 could be a huge reservoir of possible ``fuel'' for future starbursts in the central region, and we suggest that the star formation in the central region of NGC 986 could still be in a growing phase. We found a good spatial coincidence between the overall distributions of dense molecular gas traced by CO(3-2) and the massive star formation depicted by H$\alpha$. The global CO(3-2) luminosity $L'_{\rm CO(3-2)}$ of NGC 986 was determined to be $(5.4 \pm 1.1) \times 10^8$ K km s$^{-1}$ pc$^2$. The CO(3-2)/CO(1-0) integrated intensity ratio was found to be 0.60 +/- 0.13 at a spatial resolution of 44'' or 5 kpc, and a CO(3-2)/CO(2-1) ratio was 0.67 +/- 0.14 at a beam size of ~25'' or ~2.8 kpc. These line ratios suggest moderate excitation conditions of CO lines ($n_{\rm H_2} \sim 10^{3-4}$ cm$^{-3}$) in the central a few kpc region of NGC 986., Comment: PASJ, in press. See http://www.nro.nao.ac.jp/~f0212kk/00ASTE986CO32/n986-kk45.pdf for a version with full resolution figures