We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

The Large Magellanic Cloud (LMC) is the closest and most studied example of an irregular galaxy. Among its principal defining morphological features, its off-centred bar and single spiral arm stand out, defining a whole family of galaxies known as the Magellanic spirals (Sm). These structures are thought to be triggered by tidal interactions and possibly maintained via gas accretion. However, it is still unknown whether they are long-lived stable structures. In this work, by combining photometry that reaches down to the oldest main sequence turn-off in the colour-magnitude diagrams (CMD, up to a distance of $\sim$4.4 kpc from the LMC centre) from the SMASH survey and CMD fitting techniques, we find compelling evidence supporting the long-term stability of the LMC spiral arm, dating the origin of this structure to more than 2~Gyr ago. The evidence suggests that the close encounter between the LMC and the Small Magellanic Cloud (SMC) that produced the gaseous Magellanic Stream and its Leading Arm (LA) also triggered the formation of the LMC's spiral arm. Given the mass difference between the Clouds and the notable consequences of this interaction, we can speculate that this should have been one of their closest encounters. These results set important constraints on the timing of LMC-SMC collisions, as well as on the physics behind star formation induced by tidal encounters.

We report a simultaneous 44 and 95 GHz class I methanol maser survey toward 144 sources from the 95 GHz class I methanol maser catalog. The observations were made with the three telescopes of the Korean very long baseline interferometry network operating in single-dish mode. The detection rates are 89% at 44 GHz and 77% at 95 GHz. There are 106 new discoveries at 44 GHz. Comparing the previous 95 GHz detections with new observations of the same transitions made using the Purple Mountain Observatory 13.7 m radio telescope shows no clear evidence of variability on a timescale of six years. Emission from the 44 and 95 GHz transitions shows strong correlations in peak velocity, peak flux density, and integrated flux density, indicating that they are likely cospatial. We found that the peak flux density ratio Spk,95/Spk,44 decreases as the 44 GHz peak flux density increases. We found that some class I methanol masers in our sample could be associated with infrared dark clouds, while others are associated with H II regions, indicating that some sources occur at an early stage of high-mass star formation, while others are located toward more evolved sources.

We report the discovery of two ultra-faint stellar systems found in early data from the DECam Local Volume Exploration survey (DELVE). The first system, Centaurus I (DELVE J1238-4054), is identified as a resolved overdensity of old and metal-poor stars with a heliocentric distance of ${\rm D}_{\odot} = 116.3_{-0.6}^{+0.6}$ kpc, a half-light radius of $r_h = 2.3_{-0.3}^{+0.4}$ arcmin, an age of $\tau > 12.85$ Gyr, a metallicity of $Z = 0.0002_{-0.0002}^{+0.0001}$, and an absolute magnitude of $M_V = -5.55_{-0.11}^{+0.11}$ mag. This characterization is consistent with the population of ultra-faint satellites, and confirmation of this system would make Centaurus I one of the brightest recently discovered ultra-faint dwarf galaxies. Centaurus I is detected in Gaia DR2 with a clear and distinct proper motion signal, confirming that it is a real association of stars distinct from the Milky Way foreground; this is further supported by the clustering of blue horizontal branch stars near the centroid of the system. The second system, DELVE 1 (DELVE J1630-0058), is identified as a resolved overdensity of stars with a heliocentric distance of ${\rm D}_{\odot} = 19.0_{-0.6}^{+0.5} kpc$, a half-light radius of $r_h = 0.97_{-0.17}^{+0.24}$ arcmin, an age of $\tau = 12.5_{-0.7}^{+1.0}$ Gyr, a metallicity of $Z = 0.0005_{-0.0001}^{+0.0002}$, and an absolute magnitude of $M_V = -0.2_{-0.6}^{+0.8}$ mag, consistent with the known population of faint halo star clusters. Given the low number of probable member stars at magnitudes accessible with Gaia DR2, a proper motion signal for DELVE 1 is only marginally detected. We compare the spatial position and proper motion of both Centaurus I and DELVE 1 with simulations of the accreted satellite population of the Large Magellanic Cloud (LMC) and find that neither is likely to be associated with the LMC.

Context. The study of non-principal axis (NPA) rotators can provide important clues to the evolution of the spin state of asteroids. However, so far, very few studies have focused on NPA-rotating main-belt asteroids (MBAs). One of MBAs that are known to be in an excited rotation state is asteroid (5247) Krylov. Aims. By using disk-integrated photometric data, we construct a physical model of (5247) Krylov including shape and spin state. Methods. We apply the light curve convex inversion method employing optical light curves obtained by using ground-based telescopes in three apparitions during 2006, 2016, and 2017, along with infrared light curves obtained by the Wide-field Infrared Survey Explorer (WISE) satellite in 2010. Results. Asteroid (5247)~Krylov is spinning in a short axis mode (SAM) characterized by rotation and precession periods of 368.7 hr and 67.27 hr, respectively. The angular momentum vector orientation of Krylov is found to be $\lambda_{L} = 298^\circ$ and $\beta_{L} = -58^\circ$. The ratio of the rotational kinetic energy to the basic spin state energy $E/E_{0} \simeq 1.02$ shows that the (5247) Krylov is about 2% excited state compared to the Principal Axis (PA) rotation state. The shape of (5247) Krylov can be approximated by an elongated prolate ellipsoid with a ratio of moments of inertia of $I_{a}:I_{b}:I_{c}=0.36:0.96:1$. This is the first physical model of NPA rotator among MBAs. The physical processes that led to the current NPA rotation cannot be unambiguously reconstructed.

We present the first astrometry catalog from the Japanese VLBI (very long baseline interferometer) project VERA (VLBI Exploration of Radio Astrometry). We have compiled all the astrometry results from VERA, providing accurate trigonometric annual parallax and proper motion measurements. In total, 99 maser sources are listed in the VERA catalog. Among them, 21 maser sources are newly reported while the rest of 78 sources are referred to previously published results or those in preparation for forthcoming papers. The accuracy in the VERA astrometry are revisited and compared with those from the other VLBI astrometry projects such as BeSSeL (The Bar and Spiral Structure Legacy) Survey and GOBELINS (the Gould's Belt Distances Survey) with the VLBA (Very Long Baseline Array). We have confirmed that most of the astrometry results are consistent with each other, and the largest error sources are due to source structure of the maser features and their rapid variation, along with the systematic calibration errors and different analysis methods. Combined with the BeSSeL results, we estimate the up-to-date fundamental Galactic parameter of $R_{0}=7.92\pm0.16_{\rm{stat.}}\pm0.3_{\rm{sys.}}$~kpc and $\Omega_{\odot}=30.17\pm0.27_{\rm{stat.}}\pm0.3_{\rm{sys.}}$~km~s$^{-1}$~kpc$^{-1}$, where $R_{0}$ and $\Omega_{\odot}$ are the distance from the Sun to the Galactic center and the Sun's angular velocity of the Galactic circular rotation, respectively.

We report the discovery of a young (${\rm age} \sim 130~{\rm Myr}$), low-mass ($M \sim 1200~{\rm M}_\odot$), metal-poor ($[{\rm Fe}/{\rm H}] \sim -1.1$) stellar association at a heliocentric distance $D \approx 29~{\rm kpc}$, placing it far into the Milky Way halo. At its present Galactocentric position $(R, z) \sim (23, 15)~{\rm kpc}$, the association is (on the sky) near the leading arm of the gas stream emanating from the Magellanic cloud system, but is located $\approx 60^\circ$ from the Large Magellanic Cloud (LMC) center on the other side of the Milky Way disk. If we assume that the cluster is co-located with HI gas in the stream, we directly measure the distance to the leading arm of the Magellanic stream. The measured distance is inconsistent with stream predictions from models of the LMC/SMC interaction and infall into the Milky Way that do not account for ram pressure and gas interaction with Milky Way disk. The estimated age of the cluster is consistent with the time of last passage of the leading arm gas through the Galactic midplane, and we therefore speculate that this star-formation event was triggered by its the last disk midplane passage. Most details of this idea remain a puzzle: the Magellanic stream has low column density, the Milky Way disk at this large radius has low gas density, and the relative velocity of the leading arm gas and Milky Way gas is large. However it formed, the discovery of a young stellar cluster in the Milky Way halo presents an interesting opportunity for study. This cluster was discovered with Gaia astrometry and photometry alone, but folow-up DECam photometry was crucial for measuring its properties.

A galaxy's stellar mass-to-light ratio ($M_\star/L$) is a useful tool for converting luminosity to stellar mass ($M_\star$). However, the practical utility of $M_\star/L$ inferred from stellar population synthesis (SPS) models is limited by mismatches between the real and assumed models for star formation history (SFH) and dust geometry, both of which vary within galaxies. Here, we measure spatial variations in $M_\star/L$ and their dependence on color, star formation history, and dust across the disk of M31, using a map of $M^\mathrm{CMD}_\star$ derived from color-magnitude diagrams of resolved stars in the Panchromatic Hubble Andromeda Treasury (PHAT) survey. First, we find comparable scatter in $M_\star/L$ for the optical and mid-IR, contrary to the common idea that $M_\star/L$ is less variable in the IR. Second, we confirm that $M_\star/L$ is correlated with color for both the optical and mid-IR and report color vs. $M_\star/L$ relations (CMLRs) in M31 for filters used in the Sloan Digital Sky Survey (SDSS) and Widefield Infrared Survey Explorer (WISE). Third, we show that the CMLR residuals correlate with recent SFH, such that quiescent regions are offset to higher $M_\star/L$ than star-forming regions at a fixed color. The mid-IR CMLR, however, is not linear due to the high scatter of $M_\star/L$ in star-forming regions. Finally, we find a flatter optical CMLR than any SPS-based CMLRs in the literature. We show this is an effect of dust geometry, which is typically neglected but should be accounted for when using optical data to map $M_\star/L$.

We present results from a 2-epoch HST H$\alpha$ emission line survey of the Andromeda Galaxy that overlaps the footprint of the Panchromatic Hubble Andromeda Treasury (PHAT) survey. We find 552 (542) classical Be stars and 8429 (8556) normal B-type stars in epoch # 1 (epoch # 2), yielding an overall fractional Be content of 6.15% $\pm$0.26% (5.96% $\pm$0.25%). The fractional Be content decreased with spectral sub-type from $\sim$23.6% $\pm$2.0% ($\sim$23.9% $\pm$2.0%) for B0-type stars to $\sim$3.1% $\pm$0.34% ($\sim$3.4% $\pm$0.35%) for B8-type stars in epoch # 1 (epoch # 2). We observe a clear population of cluster Be stars at early fractional main sequence lifetimes, indicating that a subset of Be stars emerge onto the ZAMS as rapid rotators. Be stars are 2.8x rarer in M31 for the earliest sub-types compared to the SMC, confirming that the fractional Be content decreases in significantly more metal rich environments (like the Milky Way and M31). However, M31 does not follow a clear trend of Be fraction decreasing with metallicity compared to the Milky Way, which may reflect that the Be phenomenon is enhanced with evolutionary age. The rate of disk-loss or disk-regeneration episodes we observe, 22% $\pm$ 2% yr$^{-1}$, is similar to that observed for seven other Galactic clusters reported in the literature, assuming these latter transient fractions scale by a linear rate. The similar number of disk-loss events (57) as disk-renewal events (43) was unexpected since disk dissipation time-scales can be $\sim$2x the typical time-scales for disk build-up phases.

We examine whether galaxy environments directly affect triggering nuclear activity in Sloan Digital Sky Survey (SDSS) local spiral galaxies using a volume-limited sample with the $r$-band absolute magnitude $M_{r} < -19.0$ and $0.02 < z < 0.055$ selected from the SDSS Data Release 7. To avoid incompleteness of the central velocity dispersion $\sigma$ of the volume-limited sample and to fix the black hole mass affecting AGN activity, we limit the sample to a narrow $\sigma$ range of $130$ km s$^{-1}<\sigma<200$ km s$^{-1}$. We define a variety of environments as a combination of neighbour interactions and local density on a galaxy. After the central star formation rate (which is closely related to AGN activity level) is additionally restricted, the direct impact of the environment is unveiled. In the outskirts of rich clusters, red spiral galaxies show a significant excess of the AGN fraction despite the lack of central gas. We argue that they have been pre-processed before entering the rich clusters, and due to mergers or strong encounters in the in-fall region, their remaining gases efficiently lose angular momentum. We investigate an environment in which many star-forming galaxies coexist with a few starburst-AGN composite hosts having the highest [OIII] luminosity. We claim that they are a gas-rich merger product in groups or are group galaxies in-falling into clusters, indicating that many AGN signatures may be obscured following the merger events.

We report a search for astronomical neutrinos in the energy region from several GeV to TeV in the direction of the blazar TXS0506+056 using the Super-Kamiokande detector following the detection of a 100 TeV neutrino from the same location by the IceCube collaboration. Using Super-Kamiokande neutrino data across several data samples observed from April 1996 to February 2018 we have searched for both a total excess above known backgrounds across the entire period as well as localized excesses on smaller time scales in that interval. No significant excess nor significant variation in the observed event rate are found in the blazar direction. Upper limits are placed on the electron and muon neutrino fluxes at 90\% confidence level as $6.03 \times 10^{-7}$ and $4.52 \times 10^{-7}$ to $9.26 \times 10^{-10}$ [${\rm erg}/{\rm cm}^2/{\rm s}$], respectively.

We report spectroscopic measurements of stars in the recently discovered young stellar association Price-Whelan 1 (PW 1), which was found in the vicinity of the Leading Arm (LA) of the Magellanic Stream. We obtained Magellan+MIKE high-resolution spectra of the 28 brightest stars in PW 1 and used The Cannon to determine their stellar parameters. We find that the mean metallicity of PW 1 is [Fe/H]=-1.23 with a small scatter of 0.06 dex and the mean radial velocity is Vhelio=276.7 km/s with a dispersion of 11.0 km/s. Our results are consistent in Teff, logg, and [Fe/H] with the young and metal-poor characteristics (116 Myr and [Fe/H]=-1.1) determined for PW 1 from our discovery paper. We find a strong correlation between the spatial pattern of the PW 1 stars and the LA II gas with an offset of -10.15 deg in L_MS and +1.55 deg in B_MS. The similarity in metallicity, velocity, and spatial patterns indicates that PW 1 likely originated in LA II. We find that the spatial and kinematic separation between LA II and PW 1 can be explained by ram pressure from Milky Way gas. Using orbit integrations that account for the LMC and MW halo and outer disk gas, we constrain the halo gas density at the orbital pericenter of PW 1 to be n_halo (17 kpc) = 2.7 (3) x 10^-3 atoms/cm^3 and the disk gas density at the midplane at 20 kpc to be n_disk (20 kpc,0) = 6.0 (1.8) x 10^-2 atoms/cm^3. We, therefore, conclude that PW 1 formed from the LA II of the Magellanic Stream, making it a powerful constraint on the Milky Way-Magellanic interaction.

Supernova detection is a major objective of the Super-Kamiokande (SK) experiment. In the next stage of SK (SK-Gd), gadolinium (Gd) sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons. A core-collapse supernova will be preceded by an increasing flux of neutrinos and anti-neutrinos, from thermal and weak nuclear processes in the star, over a timescale of hours; some of which may be detected at SK-Gd. This could provide an early warning of an imminent core-collapse supernova, hours earlier than the detection of the neutrinos from core collapse. Electron anti-neutrino detection will rely on inverse beta decay events below the usual analysis energy threshold of SK, so Gd loading is vital to reduce backgrounds while maximising detection efficiency. Assuming normal neutrino mass ordering, more than 200 events could be detected in the final 12 hours before core collapse for a 15-25 solar mass star at around 200 pc, which is representative of the nearest red supergiant to Earth, $\mathrm{\alpha}$Ori (Betelgeuse). At a statistical false alarm rate of 1 per century, detection could be up to 10 hours before core collapse, and a pre-supernova star could be detected by SK-Gd up to 600 pc away. A pre-supernova alert could be provided to the astrophysics community following gadolinium loading.

We present a method to map the total intrinsic reddening of a foreground extinguishing medium via the analysis of spectral energy distributions (SEDs) of background galaxies. In this pilot study, we implement this technique in two distinct regions of the Small Magellanic Cloud (SMC) - the bar and the southern outskirts - using a combination of optical and near-infrared $ugrizYJK_{\mathrm{s}}$ broadband imaging. We adopt the LePhare $\chi^{2}$-minimisation SED-fitting routine and various samples of galaxies and/or quasi-stellar objects to investigate the intrinsic reddening. We find that only when we construct reddening maps using objects classified as galaxies with low levels of intrinsic reddening (i.e. ellipticals/lenticulars and early-type spirals), the resultant maps are consistent with previous literature determinations i.e. the intrinsic reddening of the SMC bar is higher than that in the outer environs. We employ two sets of galaxy templates - one theoretical and one empirical - to test for template dependencies in the resulting reddening maps and find that the theoretical templates imply systematically higher reddening values by up to 0.20 mag in $E(B-V)$. A comparison with previous reddening maps, based on the stellar components of the SMC, typically shows reasonable agreement. There is, however, significant variation amongst the literature reddening maps as to the level of intrinsic reddening associated with the bar. Thus, it is difficult to unambiguously state that instances of significant discrepancies are the result of appreciable levels of dust not accounted for in some literature reddening maps or whether they reflect issues with our adopted methodology.

We derive the geometric albedo of a near-Earth asteroid, (4179) Toutatis, to investigate its surface physical conditions. The asteroid has been studied rigorously not only via ground-based photometric, spectrometric, polarimetric, and radar observations but also via \textit{in situ} observation by the Chinese Chang'e-2 space probe; however, its geometric albedo is not well understood. We conducted V-band photometric observations when the asteroid was at opposition in April 2018 using the three telescopes in the southern hemisphere that compose the Korea Microlensing Telescope Network (KMTNet). The observed time-variable cross section was corrected using the radar shape model. We find that Toutatis has a geometric albedo $p_\mathrm{V} = 0.185^{+0.045}_{-0.039} $, which is typical of S-type asteroids. We compare the geometric albedo with archival polarimetric data and further find that the polarimetric slope--albedo law provides a reliable estimate for the albedo of this S-type asteroid. The thermal infrared observation also produced similar results if the size of the asteroid is updated to match the results from Chang'e-2. We conjecture that the surface of Toutatis is covered with grains smaller than that of the near-Sun asteroids including (1566) Icarus and (3200) Phaethon.

Understanding the evolutionary history of the Magellanic Clouds requires an in-depth exploration and characterization of the stellar content in their outer regions, which ultimately are key to tracing the epochs and nature of past interactions. We present new deep images of a shell-like over-density of stars in the outskirts of the Small Magellanic Cloud (SMC). The shell, also detected in photographic plates dating back to the fifties, is located at ~1.9 degr from the center of the SMC in the north-east direction.The structure and stellar content of this feature were studied with multi-band, optical data from the Survey of the MAgellanic Stellar History (SMASH) carried out with the Dark Energy Camera on the Blanco Telescope at Cerro Tololo Inter-American Observatory. We also investigate the kinematic of the stars in the shell using the Gaia Data Release 2. The shell is composed of a young population with an age ~ 150 Myr, with no contribution from an old population. Thus, it is hard to explain its origin as the remnant of a tidally disrupted stellar system. The spatial distribution of the young main-sequence stars shows a rich sub-structure, with a spiral arm-like feature emanating from the main shell and a separated small arc of young stars close to the globular cluster NGC 362. We find that the absolute g-band magnitude of the shell is M_{g,shell} = -10.78+/- 0.02, with a surface brightness of mu_{g,shell} = 25.81+/- 0.01 mag/arcsec^{2}. We have not found any evidence that this feature is of tidal origin or a bright part of a spiral arm-like structure. Instead, we suggest that the shell formed in a recent star formation event, likely triggered by an interaction with the Large Magellanic Cloud and/or the Milky Way, ~150 Myr ago.

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

We explore the stellar structure of the Large Magellanic Cloud (LMC) disk using data from the Survey of the MAgellanic Stellar History (SMASH) and the Dark Energy Survey. We detect a ring-like stellar overdensity in the red clump star count map at a radius of ~6 degrees (~5.2 kpc at the LMC distance) that is continuous over ~270 degrees in position angle and is only limited by the current data coverage. The overdensity shows an amplitude up to 2.5 times higher than that of the underlying smooth disk. This structure might be related to the multiple arms found by de Vaucouleurs. We find that the overdensity shows spatial correlation with intermediate-age star clusters, but not with young (< 1 Gyr) main-sequence stars, indicating the stellar populations associated with the overdensity are intermediate in age or older. Our findings on the LMC overdensity can be explained by either of two distinct formation mechanisms of a ring-like overdensity: (1) the overdensity formed out of an asymmetric one-armed spiral wrapping around the LMC main body, which is induced by repeated encounters with the Small Magellanic Cloud (SMC) over the last Gyr, or (2) the overdensity formed very recently as a tidal response to a direct collision with the SMC. Although the measured properties of the overdensity alone cannot distinguish between the two candidate scenarios, the consistency with both scenarios suggests that the ring-like overdensity is likely a product of tidal interaction with the SMC, but not with the Milky Way halo.

We present a study of the three-dimensional (3D) structure of the Large Magellanic Cloud (LMC) using ~2.2 million red clump (RC) stars selected from the Survey of the MAgellanic Stellar History. To correct for line-of-sight dust extinction, the intrinsic RC color and magnitude and their radial dependence are carefully measured by using internal nearly dust-free regions. These are then used to construct an accurate 2D reddening map (165 square degrees with ~10 arcmin resolution) of the LMC disk and the 3D spatial distribution of RC stars. An inclined disk model is fit to the 2D distance map yielding a best-fit inclination angle i = 25.86(+0.73,-1.39) degrees with random errors of +\-0.19 degrees and line-of-nodes position angle theta = 149.23(+6.43,-8.35) degrees with random errors of +/-0.49 degrees. These angles vary with galactic radius, indicating that the LMC disk is warped and twisted likely due to the repeated tidal interactions with the Small Magellanic Cloud (SMC). For the first time, our data reveal a significant warp in the southwestern part of the outer disk starting at rho ~ 7 degrees that departs from the defined LMC plane up to ~4 kpc toward the SMC, suggesting that it originated from a strong interaction with the SMC. In addition, the inner disk encompassing the off-centered bar appears to be tilted up to 5-15 degrees relative to the rest of the LMC disk. These findings on the outer warp and the tilted bar are consistent with the predictions from the Besla et al. simulation of a recent direct collision with the SMC.

Time-series multi-band photometry and spectrometry were performed in Nov.-Dec. 2017 to investigate the homogeneity of the surface of asteroid (3200) Phaethon. We found that Phaethon is a B-type asteroid, in agreement with previous studies, and that it shows no evidence for rotational color variation. The sub-solar latitude during our observation period was approximately 55 degree S, which corresponded to the southern hemisphere of Phaethon. Thus, we found that the southern hemisphere of Phaethon has a homogeneous surface. We compared our spectra with existing spectral data to examine the latitudinal surface properties of Phaethon. The result showed that it doesn't have a latitudinal color variation. To explain this observation, we investigated the solar-radiation heating effect on Phaethon, and the result suggested that Phaethon underwent a uniform thermal metamorphism regardless of latitude, which was consistent with our observations. Based on this result, we discuss the homogeneity of the surface of Phaethon.

We obtained, for the first time, astrometrically registered maps of the 22.2 GHz H2O and 42.8, 43.1, and 86.2 GHz SiO maser emission toward the semiregular b-type variable (SRb) R Crateris, at three epochs (2015 May 21, and 2016 January 7 and 26) using the Korean Very-long-baseline Interferometry Network. The SiO masers show a ring-like spatial structure, while the H2O maser shows a very asymmetric one-side outflow structure, which is located at the southern part of the ring-like SiO maser feature. We also found that the 86.2 GHz SiO maser spots are distributed in an inner region, compared to those of the 43.1 GHz SiO maser, which is different from all previously known distributions of the 86.2 GHz SiO masers in variable stars. The different distribution of the 86.2 GHz SiO maser seems to be related to the complex dynamics caused by the overtone pulsation mode of the SRb R Crateris. Furthermore, we estimated the position of the central star based on the ring fitting of the SiO masers, which is essential for interpreting the morphology and kinematics of a circumstellar envelope. The estimated stellar coordinate corresponds well to the position measured by Gaia.

Sample return from the near-Earth asteroid known as 25143 Itokawa was conducted as part of the Hayabusa mission, with a large number of scientific findings being derived from the returned samples. Following the Hayabusa mission, Hayabusa2 was planned, targeting sample return from a primitive asteroid. The primary target body of Hayabusa2 was asteroid 162173 Ryugu; however, it was also necessary to gather physical information for backup target selection. Therefore, we examined five asteroids spectroscopically, 43 asteroids spectrophotometrically, and 41 asteroids through periodic analysis. Hence, the physical properties of 74 near-Earth asteroids were obtained, which helped the Hayabusa2 backup target search and, also, furthered understanding of the physical properties of individual asteroids and their origins.

The DEEP-South photometric census of small Solar System bodies produces massive time-series data of variable, transient or moving objects as a by-product. To fully investigate unexplored variable phenomena, we present an application of multi-aperture photometry and FastBit indexing techniques for faster access to a portion of the DEEP-South year-one data. Our new pipeline is designed to perform automated point source detection, robust high-precision photometry and calibration of non-crowded fields which have overlap with previously surveyed areas. In this paper, we show some examples of catalog-based variability searches to find new variable stars and to recover targeted asteroids. We discover 21 new periodic variables with period ranging between 0.1 and 31 days, including four eclipsing binary systems (detached, over-contact, and ellipsoidal variables), one white dwarf/M dwarf pair candidate, and rotating variable stars. We also recover astrometry ($<\pm$1--2 arcsec level accuracy) and photometry of two targeted near-earth asteroids, 2006 DZ169 and 1996 SK, along with the small- ($\sim$0.12 mag) and relatively large-amplitude ($\sim$0.5 mag) variations of their dominant rotational signals in $R$-band.

The near-Earth asteroid 3200 Phaethon (1983 TB) is an attractive object not only from a scientific viewpoint but also because of JAXA's DESTINY+ target. The rotational lightcurve and spin properties were investigated based on the data obtained in the ground-based observation campaign of Phaethon. We aim to refine the lightcurves and shape model of Phaethon using all available lightcurve datasets obtained via optical observation, as well as our time-series observation data from the 2017 apparition. Using eight 1-2-m telescopes and an optical imager, we acquired the optical lightcurves and derived the spin parameters of Phaethon. We applied the lightcurve inversion method and SAGE algorithm to deduce the convex and non-convex shape model and pole orientations. We analysed the optical lightcurve of Phaethon and derived a synodic and a sidereal rotational period of 3.6039 h, with an axis ratio of a/b = 1.07. The ecliptic longitude (lambda) and latitude (beta) of the pole orientation were determined as (308, -52) and (322, -40) via two independent methods. A non-convex model from the SAGE method, which exhibits a concavity feature, is also presented.

The supergiant VX Sagittarii is a strong emitter of both H2O and SiO masers. However, previous VLBI observations have been performed separately, which makes it difficult to spatially trace the outward transfer of the material consecutively. Here we present the astrometrically registered, simultaneous maps of 22.2 GHz H2O and 43.1/42.8/86.2/129.3 GHz SiO masers toward VX Sagittarii. The H2O masers detected above the dust-forming layers have an asymmetric distribution. The multi-transition SiO masers are nearly circular ring, suggesting spherically symmetric wind within a few stellar radii. These results provide the clear evidence that the asymmetry in the outflow is enhanced after the smaller molecular gas clump transform into the inhomogeneous dust layers. The 129.3 GHz maser arises from the outermost region compared to that of 43.1/42.8/86.2 GHz SiO masers. The ring size of the 129.3 GHz maser is maximized around the optical maximum, suggesting that radiative pumping is dominant.

We study the effects of bulge elongation on the star formation activity in the centers of spiral galaxies using the data from the Sloan Digital Sky Survey Data Release 7. We construct a volume-limited sample of face-on spiral galaxies with $M_r < -$19.5 mag at 0.02 $\leq z <$ 0.055 by excluding barred galaxies, where the aperture of the SDSS spectroscopic fibre covers the bulges of the galaxies. We adopt the ellipticity of bulges measured by Simard et al. (2011) who performed two-dimensional bulge+disc decompositions using the SDSS images of galaxies, and identify nuclear starbursts using the fibre specific star formation rates derived from the SDSS spectra. We find a statistically significant correlation between bulge elongation and nuclear starbursts in the sense that the fraction of nuclear starbursts increases with bulge elongation. This correlation is more prominent for fainter and redder galaxies, which exhibit higher ratios of elongated bulges. We find no significant environmental dependence of the correlation between bulge elongation and nuclear starbursts. These results suggest that non-axisymmetric bulges can efficiently feed the gas into the centre of galaxies to trigger nuclear starburst activity.

We present the detection of very extended stellar populations around the Large Magellanic Cloud (LMC) out to R~21 degrees, or ~18.5 kpc at the LMC distance of 50 kpc, as detected in the Survey of the MAgellanic Stellar History (SMASH) performed with the Dark Energy Camera on the NOAO Blanco 4m Telescope. The deep (g~24) SMASH color magnitude diagrams (CMDs) clearly reveal old (~9 Gyr), metal-poor ([Fe/H]=-0.8 dex) main-sequence stars at a distance of 50 kpc. The surface brightness of these detections is extremely low with our most distant detection having 34 mag per arcsec squared in g-band. The SMASH radial density profile breaks from the inner LMC exponential decline at ~13-15 degrees and a second component at larger radii has a shallower slope with power-law index of -2.2 that contributes ~0.4% of the LMC's total stellar mass. In addition, the SMASH densities exhibit large scatter around our best-fit model of ~70% indicating that the envelope of stellar material in the LMC periphery is highly disturbed. We also use data from the NOAO Source catalog to map the LMC main-sequence populations at intermediate radii and detect a steep dropoff in density on the eastern side of the LMC (at R~8 deg) as well as an extended structure to the far northeast. These combined results confirm the existence of a very extended, low-density envelope of stellar material with disturbed shape around the LMC. The exact origin of this structure remains unclear but the leading options include a classical accreted halo or tidally stripped outer disk material.

We report on the astrometric registration of VLBI images of the SiO and water masers in OH231.8+4.2, the iconic Proto-Planetary Nebula also known as the Calabash nebula, using the KVN and Source/Frequency Phase Referencing. This, for the first time, robustly confirms the alignment of the SiO masers, close to the AGB star, driving the bi-lobe structure with the water masers in the out-flow. We are able to trace the bulk motions for the water masers over the last few decades to be 19 km/s and deduce that the age of this expansion stage is 38$\pm$2 years. The combination of this result with the distance allows a full 3D reconstruction, and confirms that the water masers lie on and expand along the known large-scale symmetry axis and that the outflow is only a few decades old, so mass loss is almost certainly on-going. Therefore we conclude that the SiO emission marks the stellar core of the nebular, the water emission traces the expansion, and that there must be multiple epochs of ejection to drive the macro-scale structure.

We report the results of a neutrino search in Super-Kamiokande for coincident signals with the first detected gravitational wave produced by a binary neutron star merger, GW170817, which was followed by a short gamma-ray burst, GRB170817A, and a kilonova/macronova. We searched for coincident neutrino events in the range from 3.5 MeV to $\sim$100 PeV, in a time window $\pm$500 seconds around the gravitational wave detection time, as well as during a 14-day period after the detection. No significant neutrino signal was observed for either time window. We calculated 90% confidence level upper limits on the neutrino fluence for GW170817. From the upward-going-muon events in the energy region above 1.6 GeV, the neutrino fluence limit is $16.0^{+0.7}_{-0.6}$ ($21.3^{+1.1}_{-0.8}$) cm$^{-2}$ for muon neutrinos (muon antineutrinos), with an error range of $\pm5^{\circ}$ around the zenith angle of NGC4993, and the energy spectrum is under the assumption of an index of $-2$. The fluence limit for neutrino energies less than 100 MeV, for which the emission mechanism would be different than for higher-energy neutrinos, is also calculated. It is $6.6 \times 10^7$ cm$^{-2}$ for anti-electron neutrinos under the assumption of a Fermi-Dirac spectrum with average energy of 20 MeV.

We present the results of a search in the Super-Kamiokande (SK) detector for excesses of neutrinos with energies above a few GeV that are in the direction of the track events reported in IceCube. Data from all SK phases (SK-I through SK-IV) were used, spanning a period from April 1996 to April 2016 and corresponding to an exposure of 225 kilotonne-years . We considered the 14 IceCube track events from a data set with 1347 livetime days taken from 2010 to 2014. We use Poisson counting to determine if there is an excess of neutrinos detected in SK in a 10 degree search cone (5 degrees for the highest energy data set) around the reconstructed direction of the IceCube event. No significant excess was found in any of the search directions we examined. We also looked for coincidences with a recently reported IceCube multiplet event. No events were detected within a $\pm$ 500 s time window around the first detected event, and no significant excess was seen from that direction over the lifetime of SK.

The Large and Small Magellanic Clouds (LMC and SMC) are unique local laboratories for studying the formation and evolution of small galaxies in exquisite detail. The Survey of the MAgellanic Stellar History (SMASH) is an NOAO community DECam survey of the Clouds mapping 480 square degrees (distributed over ~2400 square degrees at ~20% filling factor) to ~24th mag in ugriz with the goal of identifying broadly distributed, low surface brightness stellar populations associated with the stellar halos and tidal debris of the Clouds. SMASH will also derive spatially-resolved star formation histories covering all ages out to large radii from the MCs that will further complement our understanding of their formation. Here, we present a summary of the survey, its data reduction, and a description of the first public Data Release (DR1). The SMASH DECam data have been reduced with a combination of the NOAO Community Pipeline, PHOTRED, an automated PSF photometry pipeline based mainly on the DAOPHOT suite, and custom calibration software. The attained astrometric precision is ~15 mas and the accuracy is ~2 mas with respect to the Gaia DR1 astrometric reference frame. The photometric precision is ~0.5-0.7% in griz and ~1% in u with a calibration accuracy of ~1.3% in all bands. The median 5 sigma point source depths in ugriz bands are 23.9, 24.8, 24.5, 24.2, 23.5 mag. The SMASH data already have been used to discover the Hydra II Milky Way satellite, the SMASH 1 old globular cluster likely associated with the LMC, and very extended stellar populations around the LMC out to R~18.4 kpc. SMASH DR1 contains measurements of ~100 million objects distributed in 61 fields. A prototype version of the NOAO Data Lab provides data access, including a data discovery tool, SMASH database access, an image cutout service, and a Jupyter notebook server with example notebooks for exploratory analysis.

We map the star formation history across M31 by fitting stellar evolution models to color-magnitude diagrams of each 83$''{\times}$83$''$ (0.3$\times$1.4 kpc, deprojected) region of the PHAT survey outside of the innermost 6$'{\times}$12$'$ portion. We find that most of the star formation occurred prior to $\sim$8 Gyr ago, followed by a relatively quiescent period until $\sim$4 Gyr ago, a subsequent star formation episode about 2 Gyr ago and a return to relative quiescence. There appears to be little, if any, structure visible for populations with ages older than 2 Gyr, suggesting significant mixing since that epoch. Finally, assuming a Kroupa IMF from 0.1$-$100 M$_{\odot}$, we find that the total amount of star formation over the past 14 Gyr in the area over which we have fit models is 5${\times}$10$^{10}$ M$_{\odot}$. Fitting the radial distribution of this star formation and assuming azimuthal symmetry, (1.5$\pm$0.2)${\times}$10$^{11}$ M$_{\odot}$ of stars have formed in the M31 disk as a whole, (9$\pm$2)${\times}$10$^{10}$ M$_{\odot}$ of which has likely survived to the present after accounting for evolutionary effects. This mass is about one fifth of the total dynamical mass of M31.

Context. The so-called Barbarian asteroids share peculiar, but common polarimetric properties, probably related to both their shape and composition. They are named after (234) Barbara, the first on which such properties were identified. As has been suggested, large scale topographic features could play a role in the polarimetric response, if the shapes of Barbarians are particularly irregular and present a variety of scattering/incidence angles. This idea is supported by the shape of (234) Barbara, that appears to be deeply excavated by wide concave areas revealed by photometry and stellar occultations. Aims. With these motivations, we started an observation campaign to characterise the shape and rotation properties of Small Main- Belt Asteroid Spectroscopic Survey (SMASS) type L and Ld asteroids. As many of them show long rotation periods, we activated a worldwide network of observers to obtain a dense temporal coverage. Methods. We used light-curve inversion technique in order to determine the sidereal rotation periods of 15 asteroids and the con- vergence to a stable shape and pole coordinates for 8 of them. By using available data from occultations, we are able to scale some shapes to an absolute size. We also study the rotation periods of our sample looking for confirmation of the suspected abundance of asteroids with long rotation periods. Results. Our results show that the shape models of our sample do not seem to have peculiar properties with respect to asteroids with similar size, while an excess of slow rotators is most probably confirmed.

The BeSSeL Survey is mapping the spiral structure of the Milky Way by measuring trigonometric parallaxes of hundreds of maser sources associated with high-mass star formation. While parallax techniques for water masers at high frequency (22 GHz) have been well documented, recent observations of methanol masers at lower frequency (6.7 GHz) have revealed astrometric issues associated with signal propagation through the ionosphere that could significantly limit parallax accuracy. These problems displayed as a "parallax gradient" on the sky when measured against different background quasars. We present an analysis method in which we generate position data relative to an "artificial quasar" at the target maser position at each epoch. Fitting parallax to these data can significantly mitigate the problems and improve parallax accuracy.

We present observations and analysis of the low-mass star-forming region, Taurus Molecular Cloud-1 (TMC-1). CS ($J$=2-1)/N$_2$H$^+$ ($J$=1-0) and C$^{17}$O ($J$=2-1)/C$^{18}$O ($J$=2-1) were observed with FCRAO (Five College Radio Astronomy Observatory) and SRAO (Seoul Radio Astronomy Observatory), respectively. In addition, Spitzer infrared data and 1.2 mm continuum data observed with MAMBO (Max-Planck Millimetre Bolometer) are used. We also perform chemical modeling to investigate the relative molecular distributions of the TMC-1 filament. Based on Spitzer observations, there is no young stellar object along the TMC-1 filament, while five Class II and one Class I young stellar objects are identified outside the filament. The comparison between column densities calculated from dust continuum and C$^{17}$O 2-1 line emission shows that CO is depleted much more significantly in the ammonia peak than in the cyanopolyyne peak, while the column densities calculated from the dust continuum are similar at the two peaks. N$_2$H$^+$ is not depleted much in either peak. According to our chemical calculation, the differential chemical distribution in the two peaks can be explained by different timescales required to reach the same density, i.e., by different dynamical processes.

We report the identification of a bright hard X-ray source dominating the M31 bulge above 25 keV from a simultaneous NuSTAR-Swift observation. We find that this source is the counterpart to Swift J0042.6+4112, which was previously detected in the Swift BAT All-sky Hard X-ray Survey. This Swift BAT source had been suggested to be the combined emission from a number of point sources; our new observations have identified a single X-ray source from 0.5 to 50 keV as the counterpart for the first time. In the 0.5-10 keV band, the source had been classified as an X-ray binary candidate in various Chandra and XMM studies; however, since it was not clearly associated with Swift J0042.6+4112, the previous E < 10 keV observations did not generate much attention. This source has a spectrum with a soft X-ray excess (kT~ 0.2 keV) plus a hard spectrum with a power law of Gamma ~ 1 and a cutoff around 15-20 keV, typical of the spectral characteristics of accreting pulsars. Unfortunately, any potential pulsation was undetected in the NuSTAR data, possibly due to insufficient photon statistics. The existing deep HST images exclude high-mass (>3 Msun) donors at the location of this source. The best interpretation for the nature of this source is an X-ray pulsar with an intermediate-mass (<3 Msun) companion or a symbiotic X-ray binary. We discuss other possibilities in more detail.

Photon Dominated Regions (PDRs) are interfaces between the mainly ionized and mainly molecular material around young massive stars. Analysis of the physical and chemical structure of such regions traces the impact of far-ultraviolet radiation of young massive stars on their environment. We present results on the physical and chemical structure of the prototypical high UV-illumination edge-on Orion Bar PDR from an unbiased spectral line survey with a wide spectral coverage. A spectral scan from 480-1250 GHz and 1410-1910 GHz at 1.1 MHz resolution was obtained by the HIFI instrument onboard the Herschel Space Observatory. For molecules with multiple transitions we used rotational diagrams to obtain excitation temperatures and column densities. For species with a single detected transition we used an optically thin LTE approximation. In case of species with available collisional rates, we also performed a non-LTE analysis to obtain kinetic temperatures, H2 volume densities, and column densities. About 120 lines corresponding to 29 molecules (including isotopologues) have been detected in the Herschel/HIFI line survey, including 11 transitions of CO, 7 transitions of 13CO, 6 transitions of C18O, 10 transitions of H2CO, and 6 transitions of H2O. Most species trace kinetic temperatures in the range between 100 and 150 K and H2 volume densities in the range between 10^5 and 10^6 cm^-3. The species with temperatures and / or densities outside of this range include the H2CO transitions tracing a very high temperature (315 K) and density (1.4x10^6 cm^-3) component and SO corresponding to the lowest temperature (56 K) measured as a part of this line survey. The observed lines/species reveal a range of physical conditions (gas density /temperature) involving structures at high density / high pressure, obsoleting the traditional 'clump / interclump' picture of the Orion Bar.

The JAXA Hayabusa-2 mission was approved in 2010 and launched on December 3, 2014. The spacecraft will arrive at the near-Earth asteroid 162173 Ryugu in 2018 where it will perform a survey, land and obtain surface material, then depart in Dec 2019 and return to Earth in Dec 2020. We observed Ryugu with the Herschel Space Observatory in Apr 2012 at far-IR thermal wavelengths, supported by several ground-based observations to obtain optical lightcurves. We reanalysed previously published Subaru-COMICS and AKARI-IRC observations and merged them with a Spitzer-IRS data set. In addition, we used a large set of Spitzer-IRAC observations obtained in the period Jan to May, 2013. The data set includes two complete rotational lightcurves and a series of ten "point-and-shoot" observations. The almost spherical shape of the target together with the insufficient lightcurve quality forced us to combine radiometric and lightcurve inversion techniques in different ways to find the object's key physical and thermal parameters. We find that the solution which best matches our data sets leads to this C class asteroid having a retrograde rotation with a spin-axis orientation of (lambda = 310-340 deg; beta = -40+/-15 deg) in ecliptic coordinates, an effective diameter (of an equal-volume sphere) of 850 to 880 m, a geometric albedo of 0.044 to 0.050 and a thermal inertia in the range 150 to 300 Jm-2s-0.5K-1. Based on estimated thermal conductivities of the top-layer surface in the range 0.1 to 0.6 WK-1m-1, we calculated that the grain sizes are approximately equal to between 1 and 10 mm. The finely constrained values for this asteroid serve as a `design reference model', which is currently used for various planning, operational and modelling purposes by the Hayabusa2 team.

A comprehensive study on the atmospheric neutrino flux in the energy region from sub-GeV up to several TeV using the Super-Kamiokande water Cherenkov detector is presented in this paper. The energy and azimuthal spectra of the atmospheric ${\nu}_e+{\bar{\nu}}_e$ and ${\nu}_{\mu}+{\bar{\nu}}_{\mu}$ fluxes are measured. The energy spectra are obtained using an iterative unfolding method by combining various event topologies with differing energy responses. The azimuthal spectra depending on energy and zenith angle, and their modulation by geomagnetic effects, are also studied. A predicted east-west asymmetry is observed in both the ${\nu}_e$ and ${\nu}_{\mu}$ samples at 8.0 {\sigma} and 6.0 {\sigma} significance, respectively, and an indication that the asymmetry dipole angle changes depending on the zenith angle was seen at the 2.2 {\sigma} level. The measured energy and azimuthal spectra are consistent with the current flux models within the estimated systematic uncertainties. A study of the long-term correlation between the atmospheric neutrino flux and the solar magnetic activity cycle is also performed, and a weak indication of a correlation was seen at the 1.1 {\sigma} level, using SK I-IV data spanning a 20 year period. For particularly strong solar activity periods known as Forbush decreases, no theoretical prediction is available, but a deviation below the typical neutrino event rate is seen at the 2.4 {\sigma} level.

We report the results from a search in Super-Kamiokande for neutrino signals coincident with the first detected gravitational wave events, GW150914 and GW151226, using a neutrino energy range from 3.5 MeV to 100 PeV. We searched for coincident neutrino events within a time window of $\pm$500 seconds around the gravitational wave detection time. Four neutrino candidates are found for GW150914 and no candidates are found for GW151226. The remaining neutrino candidates are consistent with the expected background events. We calculated the 90\% confidence level upper limits on the combined neutrino fluence for both gravitational wave events, which depends on event energy and topologies. Considering the upward going muon data set (1.6 GeV - 100 PeV) the neutrino fluence limit for each gravitational wave event is 14 - 37 (19 - 50) cm$^{-2}$ for muon neutrinos (muon antineutrinos), depending on the zenith angle of the event. In the other data sets, the combined fluence limits for both gravitational wave events range from 2.4$\times 10^{4}$ to 7.0$\times 10^{9}$ cm$^{-2}$.

We present the Bayesian Extinction And Stellar Tool (BEAST), a probabilistic approach to modeling the dust extinguished photometric spectral energy distribution of an individual star while accounting for observational uncertainties common to large resolved star surveys. Given a set of photometric measurements and an observational uncertainty model, the BEAST infers the physical properties of the stellar source using stellar evolution and atmosphere models and constrains the line of sight extinction using a newly developed mixture model that encompasses the full range of dust extinction curves seen in the Local Group. The BEAST is specifically formulated for use with large multi-band surveys of resolved stellar populations. Our approach accounts for measurement uncertainties and any covariance between them due to stellar crowding (both systematic biases and uncertainties in the bias) and absolute flux calibration, thereby incorporating the full information content of the measurement. We illustrate the accuracy and precision possible with the BEAST using data from the Panchromatic Hubble Andromeda Treasury. While the BEAST has been developed for this survey, it can be easily applied to similar existing and planned resolved star surveys.

We present a real-time supernova neutrino burst monitor at Super-Kamiokande (SK). Detecting supernova explosions by neutrinos in real time is crucial for giving a clear picture of the explosion mechanism. Since the neutrinos are expected to come earlier than light, a fast broadcasting of the detection may give astronomers a chance to make electromagnetic radiation observations of the explosions right at the onset. The role of the monitor includes a fast announcement of the neutrino burst detection to the world and a determination of the supernova direction. We present the online neutrino burst detection system and studies of the direction determination accuracy based on simulations at SK.

The gas-solid budget of carbon in protoplanetary disks is related to the composition of the cores and atmospheres of the planets forming in them. The key gas-phase carbon carriers CO, C$^{0}$ and C$^{+}$ can now be observed in disks. The gas-phase carbon abundance in disks has not yet been well characterized, we aim to obtain new constraints on the [C]/[H] ratio in a sample of disks, and to get an overview of the strength of [CI] and warm CO emission. We carried out a survey of the CO$\,6$--$5$ and [CI]$\,1$--$0$ and $2$--$1$ lines towards $37$ disks with APEX, and supplemented it with [CII] data from the literature. The data are interpreted using a grid of models produced with the DALI code. We also investigate how well the gas-phase carbon abundance can be determined in light of parameter uncertainties. The CO$\,6$--$5$ line is detected in $13$ out of $33$ sources, the [CI]$\,1$--$0$ in $6$ out of $12$, and the [CI]$\,2$--$1$ in $1$ out of $33$. With deep integrations, the first unambiguous detections of [CI]~$1$--$0$ in disks are obtained, in TW~Hya and HD~100546. Gas-phase carbon abundance reductions of a factor $5$--$10$ or more can be identified robustly based on CO and [CI] detections. The atomic carbon detection in TW~Hya confirms a factor $100$ reduction of [C]/[H]$_{\rm gas}$ in that disk, while the data are consistent with an ISM-like carbon abundance for HD~100546. In addition, BP~Tau, T~Cha, HD~139614, HD~141569, and HD~100453 are either carbon-depleted or gas-poor disks. The low [CI]~$2$--$1$ detection rates in the survey mostly reflect insufficient sensitivity to detect T~Tauri disks. The Herbig~Ae/Be disks with CO and [CII] upper limits below the models are debris disk like systems. A roughly order of magnitude increase in sensitivity compared to our survey is required to obtain useful constraints on the gas-phase [C]/[H] ratio in most of the targeted systems.

We performed a monitoring observation of a Jupiter-Family comet, 17P/Holmes, during its 2014 perihelion passage to investigate its secular change in activity. The comet has drawn the attention of astronomers since its historic outburst in 2007, and this occasion was its first perihelion passage since then. We analyzed the obtained data using aperture photometry package and derived the Afrho parameter, a proxy for the dust production rate. We found that Afrho showed asymmetric properties with respect to the perihelion passage: it increased moderately from 100 cm at the heliocentric distance r_h=2.6-3.1 AU to a maximal value of 185 cm at r_h = 2.2 AU (near the perihelion) during the inbound orbit, while dropping rapidly to 35 cm at r_h = 3.2 AU during the outbound orbit. We applied a model for characterizing dust production rates as a function of r_h and found that the fractional active area of the cometary nucleus had dropped from 20%-40% in 2008-2011 (around the aphelion) to 0.1%-0.3% in 2014-2015 (around the perihelion). This result suggests that a dust mantle would have developed rapidly in only one orbital revolution around the sun. Although a minor eruption was observed on UT 2015 January 26 at r_h = 3.0 AU, the areas excavated by the 2007 outburst would be covered with a layer of dust (<~ 10 cm depth) which would be enough to insulate the subsurface ice and to keep the nucleus in a state of low activity.

In this paper, we present a catalogue of the spectra of bright stars observed during the sky survey using the Far-Ultraviolet Imaging Spectrograph (FIMS), which was designed primarily to observe diffuse emissions. By carefully eliminating the contamination from the diffuse background, we obtain the spectra of 70 bright stars observed for the first time with a spectral resolution of 2--3 {\AA} over the wavelength of 1370--1710 {\AA}. The far-ultraviolet spectra of an additional 139 stars are also extracted with a better spectral resolution and/or higher reliability than those of the previous observations. The stellar spectral type of the stars presented in the catalogue spans from O9 to A3. The method of spectral extraction of the bright stars is validated by comparing the spectra of 323 stars with those of the International Ultraviolet Explorer (IUE) observations.

We employ a stochastic approach to probing the origin of the log-normal distributions of halo spin in N-body simulations. After analyzing spin evolution in halo merging trees, it was found that a spin change can be characterized by a stochastic random walk of angular momentum. Also, spin distributions generated by random walks are fairly consistent with those directly obtained from N-body simulations. We derived a stochastic differential equation from a widely used spin definition and measured the probability distributions of the derived angular momentum change from a massive set of halo merging trees. The roles of major merging and accretion are also statistically analyzed in evolving spin distributions. Several factors (local environment, halo mass, merging mass ratio, and redshift) are found to influence the angular momentum change. The spin distributions generated in the mean-field or void regions tend to shift slightly to a higher spin value compared with simulated spin distributions, which seems to be caused by the correlated random walks. We verified the assumption of randomness in the angular momentum change observed in the N-body simulation and detected several degrees of correlation between walks, which may provide a clue for the discrepancies between the simulated and generated spin distributions in the voids. However, the generated spin distributions in the group and cluster regions successfully match the simulated spin distribution. We also demonstrated that the log-normality of the spin distribution is a natural consequence of the stochastic differential equation of the halo spin, which is well described by the Geometric Brownian Motion model.

Stationary density waves rotating at a constant pattern speed $\Omega_{\rm P}$ would produce age gradients across spiral arms. We test whether such age gradients are present in M81 by deriving the recent star formation histories (SFHs) of 20 regions around one of M81's grand-design spiral arms. For each region, we use resolved stellar populations to determine the SFH by modeling the observed color-magnitude diagram (CMD) constructed from archival Hubble Space Telescope (HST) F435W and F606W imaging. Although we should be able to detect systematic time delays in our spatially-resolved SFHs, we find no evidence of star formation propagation across the spiral arm. Our data therefore provide no convincing evidence for a stationary density wave with a single pattern speed in M81, and instead favor the scenario of kinematic spiral patterns that are likely driven by tidal interactions with the companion galaxies M82 and NGC 3077.

Intermediate mass protostarsprovide a bridge between theories of low- and high-mass star formation. Emerging molecular outflows can be used to determine the influence of fragmentation and multiplicity on protostellar evolution through the correlation of outflow forces of intermediate mass protostars with the luminosity. The aim of this paper is to derive outflow forces from outflows of six intermediate mass protostellar regions and validate the apparent correlation between total luminosity and outflow force seen in earlier work, as well as remove uncertainties caused by different methodology. By comparing CO 6--5 observations obtained with APEX with non-LTE radiative transfer model predictions, optical depths, temperatures, densities of the gas of the molecular outflows are derived. Outflow forces, dynamical timescales and kinetic luminosities are subsequently calculated. Outflow parameters, including the forces, were derived for all sources. Temperatures in excess of 50 K were found for all flows, in line with recent low-mass results. However, comparison with other studies could not corroborate conclusions from earlier work on intermediate mass protostars which hypothesized that fragmentation enhances outflow forces in clustered intermediate mass star formation. Any enhancement in comparison with the classical relation between outflow force and luminosity can be attributed the use of a higher excitation line and improvement in methods; They are in line with results from low-mass protostars using similar techniques. The role of fragmentation on outflows is an important ingredient to understand clustered star formation and the link between low and high-mass star formation. However, detailed information on spatial scales of a few 100 AU, covering all individual members is needed to make the necessary progress.

We measure the recent star formation history (SFH) across M31 using optical images taken with the \texit{Hubble Space Telescope} as part of the Panchromatic Hubble Andromeda Treasury (PHAT). We fit the color-magnitude diagrams in ~9000 regions that are ~100 pc $\times$ 100 pc in projected size, covering a 0.5 square degree area (~380 kpc$^2$, deprojected) in the NE quadrant of M31. We show that the SFHs vary significantly on these small spatial scales but that there are also coherent galaxy-wide fluctuations in the SFH back to ~500 Myr, most notably in M31's 10-kpc star-forming ring. We find that the 10-kpc ring is at least 400 Myr old, showing ongoing star formation over the past ~500 Myr. This indicates the presence of molecular gas in the ring over at least 2 dynamical times at this radius. We also find that the ring's position is constant throughout this time, and is stationary at the level of 1 km/s, although there is evidence for broadening of the ring due to diffusion of stars into the disk. Based on existing models of M31's ring features, the lack of evolution in the ring's position makes a purely collisional ring origin highly unlikely. We find that the global SFR has been fairly constant over the last ~500 Myr, though it does show a small increase at 50 Myr that is 1.3 times the average SFR over the past 100 Myr. During the last ~500 Myr, ~60% of all SF occurs in the 10-kpc ring. Finally, we find that in the past 100 Myr, the average SFR over the PHAT survey area is $0.28\pm0.03$ M$_\odot$ yr$^{-1}$ with an average deprojected intensity of $7.3 \times 10^{-4}$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$, which yields a total SFR of ~0.7 M$_\odot$ yr$^{-1}$ when extrapolated to the entire area of M31's disk. This SFR is consistent with measurements from broadband estimates. [abridged]

We report parallaxes and proper motions of three water maser sources in high-mass star-forming regions in the Outer Spiral Arm of the Milky Way. The observations were conducted with the Very Long Baseline Array as part of Bar and Spiral Structure Legacy Survey and double the number of such measurements in the literature. The Outer Arm has a pitch angle of 14.9 +/- 2.7 deg and a Galactocentric distance of 14.1 +/- 0.6 kpc toward the Galactic anticenter. The average motion of these sources toward the Galactic center is 10.7 +/- 2.1 km/s and we see no sign of a significant fall in the rotation curve out to 15 kpc from the Galactic center. The three-dimensional locations of these star-forming regions are consistent with a Galactic warp of several hundred parsecs from the plane.

Water is a sensitive tracer of physical conditions in star-forming regions because of its large abundance variations between hot and cold regions. We use spectrally resolved observations of rotational lines of H$_2$O and its isotopologs with Herschel/HIFI to constrain the physical conditions of the water emitting region toward the high-mass protostar AFGL2591. We use analytical estimates and rotation diagrams to estimate T$_{ex}$ and column densities of H$_2$O of the envelope, the outflow, and a foreground cloud. Furthermore, we use the non-LTE radiative transfer code to estimate the temperature and volume density of the H$_2$O emitting gas. Assuming LTE, we estimate an T$_{ex}$ of 42 K and a column density of 2$\times$10$^{14}$ cm$^{-2}$ for the envelope and 45 K and 4$\times$10$^{13}$ cm$^{-2}$ for the outflow, in beams of 4" and 30", respectively. Non-LTE models indicate a kinetic temperature of 60-230 K and a volume density of 7$\times$10$^6$-10$^8$ cm$^{-3}$ for the envelope, and a kinetic temperature of 70-90 K and a gas density of 10$^7$-10$^8$ cm$^{-3}$ for the outflow. The o/p ratio of the foreground absorption is 1.9$\pm$0.4, suggesting a low temperature. In contrast, the o/p ratio seen in absorption by the outflow is 3.5$\pm$1.0, as expected for warm gas. The water abundance in the envelope is 10$^{-9}$, similar to the low values found for other high- and low-mass protostars, suggesting that this abundance is constant during the embedded phase of high-mass star formation. The water abundance in the outflow is 10$^{-10}$, which is 10$\times$ lower than in the envelope and in the outflows of high- and low-mass protostars. Since beam size effects can only increase this estimate by a factor of 2, we suggest that the water in the outflow is affected by dissociating UV radiation due to the low extinction in the outflow lobe.

We report a detection of the effect of the large-scale velocity shear on the spatial distributions of the galactic satellites around the isolated hosts. Identifying the isolated galactic systems each of which consists of a single host galaxy and its satellites from the Seventh Data Release of the Sloan Digital Sky Survey and reconstructing linearly the velocity shear field in the local universe, we measure the alignments between the relative positions of the satellites from their isolated hosts and the principal axes of the local velocity shear tensors projected on to the plane of sky. We find a clear signal that the galactic satellites in isolated systems are located preferentially along the directions of the minor principal axes of the large-scale velocity shear field. Those galactic satellites which are spirals, brighter, located at distances larger than the projected virial radii of the hosts and belonging to the spiral hosts yield stronger alignment signals, which implies that the alignment strength depends on the formation and accretion epochs of the galactic satellites. It is also shown that the alignment strength is quite insensitive to the cosmic web environment as well as the size and luminosity of the isolated hosts. Although this result is consistent with the numerical finding of Libeskind et al. based on a N-body experiment, due to the very low significance of the observed signals it remains inconclusive whether or not the velocity shear effect on the satellite distribution is truly universal.

The ortho-to-para ratio (OPR) of H$_2$O is thought to be sensitive to the temperature of water formation. The OPR of H$_2$O is thus useful to study the formation mechanism of water. We investigate the OPR of water in the Orion PDR (Photon-dominated region), at the Orion Bar and Orion S positions, using data from {\it Herschel}/HIFI. We detect the ground-state lines of ortho- and para-H$_2$$^{18}$O in the Orion Bar and Orion S and we estimate the column densities using LTE and non-LTE methods. Based on our calculations, the ortho-to-para ratio (OPR) in the Orion Bar is 0.1 $-$ 0.5, which is unexpectedly low given the gas temperature of $\sim$ 85 K, and also lower than the values measured for other interstellar clouds and protoplanetary disks. Toward Orion S, our OPR estimate is below 2. This low OPR at 2 positions in the Orion PDR is inconsistent with gas phase formation and with thermal evaporation from dust grains, but it may be explained by photodesorption.

We investigated the magnitude-phase relation of (162173) 1999 JU3, a target asteroid for the JAXA Hayabusa 2 sample return mission. We initially employed the international Astronomical Union's H-G formalism but found that it fits less well using a single set of parameters. To improve the inadequate fit, we employed two photometric functions, the Shevchenko and Hapke functions. With the Shevchenko function, we found that the magnitude-phase relation exhibits linear behavior in a wide phase angle range (alpha = 5-75 deg) and shows weak nonlinear opposition brightening at alpha< 5 deg, providing a more reliable absolute magnitude of Hv = 19.25 +- 0.03. The phase slope (0.039 +- 0.001 mag/deg) and opposition effect amplitude (parameterized by the ratio of intensity at alpha=0.3 deg to that at alpha=5 deg, I(0.3)/I(5)=1.31+-0.05) are consistent with those of typical C-type asteroids. We also attempted to determine the parameters for the Hapke model, which are applicable for constructing the surface reflectance map with the Hayabusa 2 onboard cameras. Although we could not constrain the full set of Hapke parameters, we obtained possible values, w=0.041, g=-0.38, B0=1.43, and h=0.050, assuming a surface roughness parameter theta=20 deg. By combining our photometric study with a thermal model of the asteroid (Mueller et al. in preparation), we obtained a geometric albedo of pv = 0.047 +- 0.003, phase integral q = 0.32 +- 0.03, and Bond albedo AB = 0.014 +- 0.002, which are commensurate with the values for common C-type asteroids.

We report trigonometric parallaxes and proper motions of water masers for 12 massive star forming regions in the Perseus spiral arm of the Milky Way as part of the Bar and Spiral Structure Legacy (BeSSeL) Survey. Combining our results with 14 parallax measurements in the literature, we estimate a pitch angle of 9.9 +/- 1.5 degrees for a section of the Perseus arm. The three-dimensional Galactic peculiar motions of these sources indicate that on average they are moving toward the Galactic center and slower than the Galactic rotation.

The Korean very-long-baseline interferometry (VLBI) network (KVN) and VLBI Exploration of Radio Astrometry (VERA) Array (KaVA) is the first international VLBI array dedicated to high-frequency (23 and 43 GHz bands) observations in East Asia. Here, we report the first imaging observations of three bright active galactic nuclei (AGNs) known for their complex morphologies: 4C 39.25, 3C 273, and M 87. This is one of the initial result of KaVA early science. Our KaVA images reveal extended outflows with complex substructure such as knots and limb brightening, in agreement with previous Very Long Baseline Array (VLBA) observations. Angular resolutions are better than 1.4 and 0.8 milliarcsecond at 23 GHz and 43 GHz, respectively. KaVA achieves a high dynamic range of ~1000, more than three times the value achieved by VERA. We conclude that KaVA is a powerful array with a great potential for the study of AGN outflows, at least comparable to the best existing radio interferometric arrays.

Large-scale structure of the universe is a useful cosmological probe of the primordial non-Gaussianity and the expansion history of the universe because its topology does not change with time in the linear regime in the standard paradigm of structure formation. However, when the topology of iso-density contour surfaces is measured from an observational data, many systematic effects are introduced due to the finite size of pixels used to define the density field, non-linear gravitational evolution, redshift-space distortion, shot noise (discrete sampling), and bias in the distribution of the density field tracers. We study the various systematic effects on the genus curve to a great accuracy by using the Horizon Run 2 simulation of a {\Lambda}CDM cosmology. We numerically measure the genus curve from the gravitationally evolved matter and dark matter halo density fields. It is found that all the non-Gaussian deviations due to the systematic effects can be modeled by using a few low-order Hermite polynomials from H0 to H4. We compare our results with the analytic theories whenever possible, and find many new terms in the Hermite series that are making significant contributions to the non-Gaussian deviations. In particular, it is found that the amplitude drop of the genus curve due to the non-linear gravitational evolution can be accurately modeled by two terms H0 and H2 with coefficients both proportional to \sigma_0^2, the mean-square density fluctuation.

Over 100 trigonometric parallaxes and proper motions for masers associated with young, high-mass stars have been measured with the BeSSeL Survey, a VLBA key science project, the EVN, and the Japanese VERA project. These measurements provide strong evidence for the existence of spiral arms in the Milky Way, accurately locating many arm segments and yielding spiral pitch angles ranging from 7 to 20 degrees. The widths of spiral arms increase with distance from the Galactic center. Fitting axially symmetric models of the Milky Way with the 3-D position and velocity information and conservative priors for the solar and average source peculiar motions, we estimate the distance to the Galactic center, Ro, to be 8.34 +/- 0.16 kpc, a circular rotation speed at the Sun, To, to be 240 +/- 8 km/s, and a rotation curve that is nearly flat (a slope of -0.2 +/- 0.4 km/s/kpc) between Galactocentric radii of 5 and 16 kpc. Assuming a "universal" spiral galaxy form for the rotation curve, we estimate the thin disk scale length to be 2.44 +/- 0.16 kpc. The parameters Ro and To are not highly correlated and are relatively insensitive to different forms of the rotation curve. Adopting a theoretically motivated prior that high-mass star forming regions are in nearly circular Galactic orbits, we estimate a global solar motion component in the direction of Galactic rotation, Vsun = 14.6 +/- 5.0 km/s. While To and Vsun are significantly correlated, the sum of these parameters is well constrained, To + Vsun = 255.2 +/- 5.1 km/s, as is the angular speed of the Sun in its orbit about the Galactic center, (To + Vsun)/Ro = 30.57 +/- 0.43 km/s/kpc. These parameters improve the accuracy of estimates of the accelerations of the Sun and the Hulse-Taylor binary pulsar in their Galactic orbits, significantly reducing the uncertainty in tests of gravitational radiation predicted by general relativity.

Procedures and results on hardware level detector calibration in Super-Kamiokande (SK) are presented in this paper. In particular, we report improvements made in our calibration methods for the experimental phase IV in which new readout electronics have been operating since 2008. The topics are separated into two parts. The first part describes the determination of constants needed to interpret the digitized output of our electronics so that we can obtain physical numbers such as photon counts and their arrival times for each photomultiplier tube (PMT). In this context, we developed an in-situ procedure to determine high-voltage settings for PMTs in large detectors like SK, as well as a new method for measuring PMT quantum efficiency and gain in such a detector. The second part describes the modeling of the detector in our Monte Carlo simulation, including in particular the optical properties of its water target and their variability over time. Detailed studies on the water quality are also presented. As a result of this work, we achieved a precision sufficient for physics analysis over a wide energy range (from a few MeV to above a TeV). For example, the charge determination was understood at the 1% level, and the timing resolution was 2.1 nsec at the one-photoelectron charge level and 0.5 nsec at the 100-photoelectron charge level.

With upcoming all sky surveys such as LSST poised to generate a deep digital movie of the optical sky, variability-based AGN selection will enable the construction of highly-complete catalogs with minimum contamination. In this study, we generate $g$-band difference images and construct light curves for QSO/AGN candidates listed in SDSS Stripe 82 public catalogs compiled from different methods, including spectroscopy, optical colors, variability, and X-ray detection. Image differencing excels at identifying variable sources embedded in complex or blended emission regions such as Type II AGNs and other low-luminosity AGNs that may be omitted from traditional photometric or spectroscopic catalogs. To separate QSOs/AGNs from other sources using our difference image light curves, we explore several light curve statistics and parameterize optical variability by the characteristic damping timescale ($\tau$) and variability amplitude. By virtue of distinguishable variability parameters of AGNs, we are able to select them with high completeness of 93.4% and efficiency (i.e., purity) of 71.3%. Based on optical variability, we also select highly variable blazar candidates, whose infrared colors are consistent with known blazars. One third of them are also radio detected. With the X-ray selected AGN candidates, we probe the optical variability of X-ray detected optically-extended sources using their difference image light curves for the first time. A combination of optical variability and X-ray detection enables us to select various types of host-dominated AGNs. Contrary to the AGN unification model prediction, two Type II AGN candidates (out of 6) show detectable variability on long-term timescales like typical Type I AGNs. This study will provide a baseline for future optical variability studies of extended sources.

Maria family is regarded as an old-type (~3 +/- 1 Gyr) asteroid family which has experienced substantial collisional and dynamical evolution in the Main-belt. It is located nearby the 3:1 Jupter mean motion resonance area that supplies Near-Earth asteroids (NEAs) to the inner Solar System. We carried out observations of Maria family asteroids during 134 nights from 2008 July to 2013 May, and derived synodic rotational periods for 51 objects, including newly obtained periods of 34 asteroids. We found that there is a significant excess of fast and slow rotators in observed rotation rate distribution. The two-sample Kolmogorov-Smirnov test confirms that the spin rate distribution is not consistent with a Maxwellian at a 92% confidence level. From correlations among rotational periods, amplitudes of lightcurves, and sizes, we conclude that the rotational properties of Maria family asteroids have been changed considerably by non-gravitational forces such as the YORP effect. Using a lightcurve inversion method (Kaasalainen & Torppa 2001; Kaasalainen et al. 2001), we successfully determined the pole orientations for 13 Maria members, and found an excess of prograde versus retrograde spins with a ratio (N_p/N_r) of 3. This implies that the retrograde rotators could have been ejected by the 3:1 resonance into the inner Solar System since the formation of Maria family. We estimate that approximately 37 to 75 Maria family asteroids larger than 1 km have entered the near-Earth space every 100 Myr.

The results of a maser survey toward ninety-nine protostars in the Orion molecular cloud complex are presented. The target sources are low-mass protostars identified from infrared observations. Single-dish observations were carried out in the water maser line at 22 GHz and the methanol class I maser lines at 44, 95, and 133 GHz. Most of the detected sources were mapped to determine the source positions. Five water maser sources were detected, and they are excited by HH 1-2 VLA 3, HH 1-2 VLA 1, L1641N MM1/3, NGC 2071 IRS 1/3, and an object in the OMC 3 region. The water masers showed significant variability in intensity and velocity with time scales of a month or shorter. Four methanol emission sources were detected, and those in the OMC 2 FIR 3/4 and L1641N MM1/3 regions are probably masers. The methanol emission from the other two sources in the NGC 2071 IRS 1-3 and V380 Ori NE regions are probably thermal. For the water masers, the number of detections per protostar in the survey region is about 2%, which suggests that the water masers of low-mass protostars are rarely detectable. The methanol class I maser of low-mass protostars is an even rarer phenomenon, with a detection rate much smaller than 1%.

We present measurements of the genus topology of luminous red galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS) Data Release 7 catalog, with unprecedented statistical significance. To estimate the uncertainties in the measured genus, we construct 81 mock SDSS LRG surveys along the past light cone from the Horizon Run 3, one of the largest N-body simulations to date that evolved 7210^3 particles in a 10815 Mpc/h size box. After carefully modeling and removing all known systematic effects due to finite pixel size, survey boundary, radial and angular selection functions, shot noise and galaxy biasing, we find the observed genus amplitude to reach 272 at 22 Mpc/h smoothing scale with an uncertainty of 4.2%; the estimated error fully incorporates cosmic variance. This is the most accurate constraint of the genus amplitude to date, which significantly improves on our previous results. In particular, the shape of the genus curve agrees very well with the mean topology of the SDSS LRG mock surveys in the LCDM universe. However, comparison with simulations also shows small deviations of the observed genus curve from the theoretical expectation for Gaussian initial conditions. While these discrepancies are mainly driven by known systematic effects such as those of shot noise and redshift-space distortions, they do contain important cosmological information on the physical effects connected with galaxy formation, gravitational evolution and primordial non-Gaussianity. We address here the key role played by systematics on the genus curve, and show how to accurately correct for their effects to recover the topology of the underlying matter. In a forthcoming paper, we provide an interpretation of those deviations in the context of the local model of non-Gaussianity.

We report the analysis results of far ultraviolet (FUV) observations, made for a broad region around $\alpha$ Vir (Spica) including the interaction zone of Loop I and the Local Bubble. The whole region was optically thin and a general correlation was seen between the FUV continuum intensity and the dust extinction, except in the neighborhood of the bright central star, indicating the dust scattering nature of the FUV continuum. We performed Monte-Carlo radiative transfer simulations to obtain the optical parameters related to the dust scattering as well as the geometrical structure of the region. The albedo and asymmetry factor were found to be 0.38$\pm$0.06 and 0.46$\pm$0.06, respectively, in good agreement with the Milky Way dust grain models. The distance to and the thickness of the interaction zone were estimated to be 70$^{+4}_{-8}$ pc and 40$^{+8}_{-10}$ pc, respectively. The diffuse FUV continuum in the northern region above Spica was mostly the result of scattering of the starlight from Spica, while that in the southern region was mainly due to the background stars. The \ion{C}{4} $\lambda\lambda$1548, 1551 emission was found throughout the whole region, in contrast to the \ion{Si}{2}* $\lambda$1532 emission which was bright only within the \ion{H}{2} region. This indicates that the \ion{C}{4} line arises mostly at the shell boundaries of the bubbles, with a larger portion likely from the Loop I than from the Local Bubble side, whereas the \ion{Si}{2}* line is from the photoionized Spica nebula.

Near-Earth asteroid 162173 (1999 JU3) is a potential target of two asteroid sample return missions, not only because of its accessibility but also because of the first C-type asteroid for exploration missions. The lightcurve-related physical properties of this object were investigated during the 2011-2012 apparition. We aim to confirm the physical parameters useful for JAXA's Hayabusa 2 mission, such as rotational period, absolute magnitude, and phase function. Our data complement previous studies that did not cover low phase angles. With optical imagers and 1-2 m class telescopes, we acquired the photometric data at different phase angles. We independently derived the rotational lightcurve and the phase curve of the asteroid. We have analyzed the lightcurve of 162173 (1999 JU3), and derived a synodic rotational period of 7.625 +/- 0.003 h, the axis ratio a/b = 1.12. The absolute magnitude H_R = 18.69 +/- 0.07 mag and the phase slope of G = -0.09 +/- 0.03 were also obtained based on the observations made during the 2011-2012 apparition.

The Swift Gamma-ray Burst (GRB) observatory responds to GRB triggers with optical observations in ~ 100 s, but cannot respond faster than ~ 60 s. While some ground-based telescopes respond quickly, the number of sub-60 s detections remains small. In mid- to late-2013, the Ultra-Fast Flash Observatory-Pathfinder is to be launched on the Lomonosov spacecraft to investigate early optical GRB emission. This pathfinder mission is necessarily limited in sensitivity and event rate; here we discuss a next generation rapid-response space observatory. We list science topics motivating our instruments, those that require rapid optical-IR GRB response, including: A survey of GRB rise shapes/times, measurements of optical bulk Lorentz factors, investigation of magnetic dominated (vs. non-magnetic) jet models, internal vs. external shock origin of prompt optical emission, the use of GRBs for cosmology, and dust evaporation in the GRB environment. We also address the impacts of the characteristics of GRB observing on our instrument and observatory design. We describe our instrument designs and choices for a next generation observatory as a second instrument on a low-earth orbit spacecraft, with a 120 kg instrument mass budget. Restricted to relatively modest mass and power, we find that a coded mask X-ray camera with 1024 cm2 of detector area could rapidly locate about 64 GRB triggers/year. Responding to the locations from the X-ray camera, a 30 cm aperture telescope with a beam-steering system for rapid (~ 1 s) response and a near-IR camera should detect ~ 29 GRB, given Swift GRB properties. Am additional optical camera would give a broadband optical-IR slope, allowing dynamic measurement of dust extinction at the source, for the first time.

We present analyses to determine the fundamental parameters of the Galaxy based on VLBI astrometry of 52 Galactic maser sources obtained with VERA, VLBA and EVN. We model the Galaxy's structure with a set of parameters including the Galaxy center distance R_0, the angular rotation velocity at the LSR Omega_0, mean peculiar motion of the sources with respect to Galactic rotation (U_src, V_src, W_src), rotation-curve shape index, and the V component of the Solar peculiar motions V_sun. Based on a Markov chain Monte Carlo method, we find that the Galaxy center distance is constrained at a 5% level to be R_0 = 8.05 +/- 0.45 kpc, where the error bar includes both statistical and systematic errors. We also find that the two components of the source peculiar motion U_src and W_src are fairly small compared to the Galactic rotation velocity, being U_src = 1.0 +/- 1.5 km/s and W_src = -1.4 +/- 1.2 km/s. Also, the rotation curve shape is found to be basically flat between Galacto-centric radii of 4 and 13 kpc. On the other hand, we find a linear relation between V_src and V_sun as V_src = V_sun -19 (+/- 2) km/s, suggesting that the value of V_src is fully dependent on the adopted value of V_sun. Regarding the rotation speed in the vicinity of the Sun, we also find a strong correlation between Omega_0 and V_sun. We find that the angular velocity of the Sun, Omega_sun, which is defined as Omega_sun = Omega_0 + V_sun/R_0, can be well constrained with the best estimate of Omega_sun = 31.09 +/- 0.78 km/s/kpc. This corresponds to Theta_0 = 238 +/- 14 km/s if one adopts the above value of R_0 and recent determination of V_sun ~ 12 km/s.

Large galaxy redshift surveys have long been used to constrain cosmological models and structure formation scenarios. In particular, the largest structures discovered observationally are thought to carry critical information on the amplitude of large-scale density fluctuations or homogeneity of the universe, and have often challenged the standard cosmological framework. The Sloan Great Wall (SGW) recently found in the Sloan Digital Sky Survey (SDSS) region casts doubt on the concordance cosmological model with a cosmological constant (i.e. the flat LCDM model). Here we show that the existence of the SGW is perfectly consistent with the LCDM model, a result that only our very large cosmological N-body simulation (the Horizon Run 2, HR2) could supply. In addition, we report on the discovery of a void complex in the SDSS much larger than the SGW, and show that such size of the largest void is also predicted in the LCDM paradigm. Our results demonstrate that an initially homogeneous isotropic universe with primordial Gaussian random phase density fluctuations growing in accordance with the General Relativity, can explain the richness and size of the observed large-scale structures in the SDSS. Using the HR2 simulation we predict that a future galaxy redshift survey about four times deeper or with 3 magnitude fainter limit than the SDSS should reveal a largest structure of bright galaxies about twice as big as the SGW.

We present the result of an indirect search for high energy neutrinos from WIMP annihilation in the Sun using upward-going muon (upmu) events at Super-Kamiokande. Datasets from SKI-SKIII (3109.6 days) were used for the analysis. We looked for an excess of neutrino signal from the Sun as compared with the expected atmospheric neutrino background in three upmu categories: stopping, non-showering, and showering. No significant excess was observed. The 90% C.L. upper limits of upward-going muon flux induced by WIMPs of 100 GeV/c$^2$ were 6.4$\times10^{-15}$ cm$^{-2}$ sec$^{-1}$ and 4.0$\times10^{-15}$ cm$^{-2}$ sec$^{-1}$ for the soft and hard annihilation channels, respectively. These limits correspond to upper limits of 4.5$\times10^{-39}$ cm$^{-2}$ and 2.7$\times10^{-40}$ cm$^{-2}$ for spin-dependent WIMP-nucleon scattering cross sections in the soft and hard annihilation channels, respectively.

The influence of a galaxy's environment on its evolution has been studied and compared extensively in the literature, although differing techniques are often used to define environment. Most methods fall into two broad groups: those that use nearest neighbours to probe the underlying density field and those that use fixed apertures. The differences between the two inhibit a clean comparison between analyses and leave open the possibility that, even with the same data, different properties are actually being measured. In this work we apply twenty published environment definitions to a common mock galaxy catalogue constrained to look like the local Universe. We find that nearest neighbour-based measures best probe the internal densities of high-mass haloes, while at low masses the inter-halo separation dominates and acts to smooth out local density variations. The resulting correlation also shows that nearest neighbour galaxy environment is largely independent of dark matter halo mass. Conversely, aperture-based methods that probe super-halo scales accurately identify high-density regions corresponding to high mass haloes. Both methods show how galaxies in dense environments tend to be redder, with the exception of the largest apertures, but these are the strongest at recovering the background dark matter environment. We also warn against using photometric redshifts to define environment in all but the densest regions. When considering environment there are two regimes: the 'local environment' internal to a halo best measured with nearest neighbour and 'large-scale environment' external to a halo best measured with apertures. This leads to the conclusion that there is no universal environment measure and the most suitable method depends on the scale being probed.

We present an analytic model for the local bias of dark matter halos in a LCDM universe. The model uses the halo mass density instead of the halo number density and is searched for various halo mass cuts, smoothing lengths, and redshift epoches. We find that, when the logarithmic density is used, the second-order polynomial can fit the numerical relation between the halo mass distribution and the underlying matter distribution extremely well. In this model the logarithm of the dark matter density is expanded in terms of log halo mass density to the second order. The model remains excellent for all halo mass cuts (from M_{cut}=3\times10^{11}$ to $3\times10^{12}h^{-1}M_{\odot}$), smoothing scales (from $R=5h^{-1}$Mpc to $50h^{-1}$Mpc), and redshift ranges (from z=0 to 1.0) considered in this study. The stochastic term in the relation is found not entirely random, but a part of the term can be determined by the magnitude of the shear tensor.

GUT monopoles captured by the Sun's gravitation are expected to catalyze proton decays via the Callan-Rubakov process. In this scenario, protons, which initially decay into pions, will ultimately produce \nu_{e}, \nu_{\mu} and \bar{\nu}_{\mu}. After undergoing neutrino oscillation, all neutrino species appear when they arrive at the Earth, and can be detected by a 50,000 metric ton water Cherenkov detector, Super-Kamiokande (SK). A search for low energy neutrinos in the electron total energy range from 19 to 55 MeV was carried out with SK and gives a monopole flux limit of F_M(\sigma_0/1 mb) < 6.3 \times 10^{-24} (\beta_M/10^{-3})^2 cm^{-2} s^{-1} sr^{-1} at 90% C.L., where \beta_M is the monopole velocity in units of the speed of light and \sigma_0 is the catalysis cross section at \beta_M=1. The obtained limit is more than eight orders of magnitude more stringent than the current best cosmic-ray supermassive monopole flux limit, F_M < 1 \times 10^{-15} cm^{-2} s^{-1} sr^{-1} for \beta_M < 10^{-3} and also two orders of magnitude lower than the result of the Kamiokande experiment, which used a similar detection method.

We study the environmental effects on the activity in galactic nuclei by comparing galaxies in clusters and in the field. Using a spectroscopic sample of galaxies in Abell clusters from the SDSS DR7, we investigate the dependence of nuclear activity on the physical parameters of clusters as well as the nearest neighbor galaxy. We also compare galaxy properties between AGN hosts and non-AGN galaxies. We find that the AGN fraction of early-type galaxies starts to decrease around one virial radius of clusters (r_{200,cl}) as decreasing clustercentric radius, while that of late types starts to decrease close to the cluster center (R~0.1-0.5r_{200,cl}). The AGN fractions of early-type cluster galaxies on average are found to be lower than those of early-type field galaxies by a factor ~3. However, the mean AGN fractions of late-type cluster galaxies are similar to those of late-type field galaxies. The AGN fraction of early-type BCGs lies between those of other early-type, cluster and field galaxies with similar luminosities. In the field, the AGN fraction strongly depends on the morphology of and the distance to the nearest neighbor galaxy. We find an anti-correlation between the AGN fraction and the velocity dispersion of clusters for all subsamples divided by morphology and luminosity of host galaxies. The AGN power indicated by L_[OIII]/M_BH is found to depend strongly on the mass of host galaxies rather than the clustercentric radius. The difference in physical parameters such as luminosity, (u-r) colors, SFRs, and (g-i) color gradients between AGN hosts and non-AGN galaxies is seen for both early and late types at all clustercentric radii, while the difference in structure parameters between the two is significant only for late types. These results support the idea that the nuclear activity is triggered through galaxy-galaxy interactions and mergers when gas supply for AGN is available.

We investigate the dependence of occurrence of bars in galaxies on galaxy properties and environment. We use a volume-limited sample of 33,391 galaxies brighter than $M_{r}=-19.5+5$log$h$ at $0.02\le z\le0.05489$, drawn from the SDSS DR 7. We classify the galaxies into early and late types, and identify bars by visual inspection. Among 10,674 late-type galaxies with axis ratio $b/a>0.60$, we find 3,240 barred galaxies ($f_{bar}=30.4%$) which divide into 2,542 strong bars ($f_{SB1}=23.8%$) and 698 weak bars ($f_{SB2}=6.5%$). We find that $f_{SB1}$ increases as $u-r$ color becomes redder, and that it has a maximum value at intermediate velocity dispersion ($\sigma\simeq$150 km s$^{-1}$). This trend suggests that strong bars are dominantly hosted by intermediate-mass systems. Weak bars prefer bluer galaxies with lower mass and lower concentration. In the case of strong bars, their dependence on the concentration index appears only for massive galaxies with $\sigma>150$ km s$^{-1}$. We also find that $f_{bar}$ does not directly depend on the large-scale background density when other physical parameters ($u-r$ color or $\sigma$) are fixed. We discover that $f_{SB1}$ decreases as the separation to the nearest neighbor galaxy becomes smaller than 0.1 times the virial radius of the neighbor regardless of neighbor's morphology. These results imply that strong bars are likely to be destroyed during strong tidal interactions, and that the mechanism for this phenomenon is gravitational and not hydrodynamical. The fraction of weak bars has no correlation with environmental parameters. We do not find any direct evidence for environmental stimulation of bar formation.

The Ultra-Fast Flash Observatory (UFFO) Pathfinder is a payload on the Russian Lomonosov satellite, scheduled to be launched in November 2011. The Observatory is designed to detect early UV/Optical photons from Gamma-Ray Bursts (GRBs). There are two telescopes and one main data acquisition system: the UFFO Burst Alert & Trigger Telescope (UBAT), the Slewing Mirror Telescope (SMT), and the UFFO Data Acquisition (UDAQ) system. The UDAQ controls and manages the operation and communication of each telescope, and is also in charge of the interface with the satellite. It will write the data taken by each telescope to the NOR flash memory and sends them to the satellite via the Bus-Interface system (BI). It also receives data from the satellite including the coordinates and time of an external trigger from another payload, and distributes them to two telescopes. These functions are implemented in field programmable gates arrays (FPGA) for low power consumption and fast processing without a microprocessor. The UDAQ architecture, control of the system, and data flow will be presented.

The Ultra-Fast Flash Observatory (UFFO) pathfinder is a space mission devoted to the measurement of Gamma-Ray Bursts (GRBs), especially their early light curves which will give crucial information on the progenitor stars and central engines of the GRBs. It consists of two instruments: the UFFO Burst Alert & Trigger telescope (UBAT) for the detection of GRB locations and the Slewing Mirror Telescope (SMT) for the UV/optical afterglow observations, upon triggering by UBAT. The UBAT employs a coded-mask {\gamma}/X-ray camera with a wide field of view (FOV), and is comprised of three parts: a coded mask, a hopper, and a detector module (DM). The UBAT DM consists of a LYSO scintillator crystal array, multi-anode photo multipliers, and analog and digital readout electronics. We present here the design and fabrication of the UBAT DM, as well as its preliminary test results.

The Ultra-Fast Flash Observatory (UFFO) is a new space-based experiment to observe Gamma-Ray Bursts (GRBs). GRBs are the most luminous electromagnetic events in the universe and occur randomly in any direction. Therefore the UFFO consists of two telescopes; UFFO Burst Alert & Trigger Telescope (UBAT) to detect GRBs using a wide field-of-view (FOV), and a Slewing Mirror Telescope (SMT) to observe UV/optical events rapidly within the narrow, targeted FOV. The SMT is a Ritchey-Chretien telescope that uses a motorized mirror system and an Intensified Charge-Coupled Device (ICCD). When the GRB is triggered by the UBAT, the SMT receives the position information and rapidly tilts the mirror to the target. The ICCD start to take the data within a second after GRB is triggered. Here we give the details about the SMT readout electronics that deliver the data.

Hundreds of gamma-ray burst (GRB) optical light curves have been measured since the discovery of optical afterglows. However, even after nearly 7 years of operation of the Swift Observatory, only a handful of measurements have been made soon (within a minute) after the gamma ray signal. This lack of early observations fails to address burst physics at short time scales associated with prompt emissions and progenitors. Because of this lack of sub-minute data, the characteristics of the rise phase of optical light curve of short-hard type GRB and rapid-rising GRB, which may account for ~30% of all GRB, remain practically unknown. We have developed methods for reaching sub-minute and sub-second timescales in a small spacecraft observatory. Rather than slewing the entire spacecraft to aim the optical instrument at the GRB position, we use rapidly moving mirror to redirect our optical beam. As a first step, we employ motorized slewing mirror telescope (SMT), which can point to the event within 1s, in the UFFO Pathfinder GRB Telescope onboard the Lomonosov satellite to be launched in Nov. 2011. UFFO's sub-minute measurements of the optical emission of dozens of GRB each year will result in a more rigorous test of current internal shock models, probe the extremes of bulk Lorentz factors, provide the first early and detailed measurements of fast-rise GRB optical light curves, and help verify the prospect of GRB as a new standard candle. We will describe the science and the mission of the current UFFO Pathfinder project, and our plan of a full-scale UFFO-100 as the next step.

Using an extensive sample of galaxies selected from the Sloan Digital Sky Survey Data Release 5, we compare the angular momentum distribution of active galactic nucleus (AGN) with non-AGN hosting late-type galaxies. To this end we characterize galactic spin through the dimensionless angular momentum parameter $\lambda$, which we estimate approximately through simple dynamical considerations. Using a volume limited sample, we find a considerable difference when comparing the empirical distributions of $\lambda$ for AGNs and non-AGN galaxies, the AGNs showing typically low $\lambda$ values and associated dispersions, while non-AGNs present higher $\lambda$ values and a broader distribution. A more striking difference is found when looking at $\lambda$ distributions in thin $M_{r}$ cuts, while the spin of non-AGN galaxies presents an anti-correlation with $M_{r}$, with bright (massive) galaxies having low spins, AGN host galaxies present uniform values of $\lambda$ at all magnitudes, a behavior probably imposed by the fact that most late-type AGN galaxies present a narrow range in color, with a typical constant $\lambda$ value. We also find that the fraction of AGN hosting galaxies in our sample strongly depends on galactic spin, increasing dramatically for decreasing $\lambda$. For AGN host galaxies we compute the mass of their supermassive black holes and find that this value tends to be higher for low spin galaxies, even at fixed luminosity, a result that could account, to a certain extent, for the spread on the luminosity black-hole mass relation.

We present the results of a survey of radio continuum sources near the Galactic plane using the Very Long Baseline Array (VLBA). Our observations are designed to identify compact extragalactic sources of milliarcsecond size that can be used for parallax measurements in the Bar and Spiral Structure Legacy Survey. We selected point sources from the NVSS and CORNISH catalogs with flux densities above 30 mJy and within $1.5\degr$ of known maser targets. Of the 1529 sources observed, 199 were detected. For sources detected on 3 or more baselines, we determined accurate positions and evaluated their quality as potential calibrators. Most of the 1330 sources that were not detected with the VLBA are probably of extragalactic origin.

We study the distribution of cosmic voids and void galaxies using Sloan Digital Sky Survey Data Release 7 (SDSS DR7). Using the VoidFinder algorithm as described by Hoyle 2002, we identify 1054 statistically significant voids in the northern galactic hemisphere with radii > 10 h^{-1} Mpc. The filling factor of voids in the sample volume is 62%. The largest void is just over 30 h^{-1} Mpc in effective radius. The median effective radius is 17 h^{-1} Mpc. The voids are found to be significantly underdense, with density contrast \delta < -0.85 at the edges of the voids. The radial density profiles of these voids are similar to predictions of dynamically distinct underdensities in gravitational theory. We find 8,046 galaxies brighter than M_r = -20.09 within the voids, accounting for 7% of the galaxies. We compare the results of VoidFinder on SDSS DR7 to mock catalogs generated from a SPH halo model simulation as well as other \Lambda -CDM simulations and find similar void fractions and void sizes in the data and simulations. This catalog is made publicly available at http://www.physics.drexel.edu/~pan/voidcatalog.html for download.

Astrometric Very Long Baseline Interferometry (VLBI) observations of maser sources in the Milky Way are used to map the spiral structure of our Galaxy and to determine fundamental parameters such as the rotation velocity ($\Theta_0$) and curve and the distance to the Galactic center (R$_0$). Here, we present an update on our first results, implementing a recent change in the knowledge about the Solar motion. It seems unavoidable that the IAU recommended values for R$_0$ and $\Theta_0$ need a substantial revision. In particular the combination of 8.5 kpc and 220 \kms\, can be ruled out with high confidence. Combining the maser data with the distance to the Galactic center from stellar orbits and the proper motion of Sgr\,A* gives best values of R$_0$ = 8.3 $\pm$ 0.23 kpc and $\Theta_0$ = 239 or 246 $\pm$ 7 \kms, for Solar motions of V$_ \odot$ = 12.23 and 5.25 \kms, respectively. Finally, we give an outlook to future observations in the Bar and Spiral Structure Legacy (BeSSeL) Survey.

Near-Earth asteroid 162173 (1999 JU3) is a potential flyby and rendezvous target for interplanetary missions because of its easy to reach orbit. The physical and thermal properties of the asteroid are relevant for establishing the scientific mission goals and also important in the context of near-Earth object studies in general. Our goal was to derive key physical parameters such as shape, spin-vector, size, geometric albedo, and surface properties of 162173 (1999 JU3). With three sets of published thermal observations (ground-based N-band, Akari IRC, Spitzer IRS), we applied a thermophysical model to derive the radiometric properties of the asteroid. The calculations were performed for the full range of possible shape and spin-vector solutions derived from the available sample of visual lightcurve observations. The near-Earth asteroid 162173 (1999 JU3) has an effective diameter of 0.87 +/- 0.03 km and a geometric albedo of 0.070 +/- 0.006. The chi2-test reveals a strong preference for a retrograde sense of rotation with a spin-axis orientation of lambda_ecl = 73 deg, beta_ecl = -62 deg and P_sid = 7.63 +/- 0.01 h. The most likely thermal inertia ranges between 200 and 600 Jm-2s-0.5K-1, about a factor of 2 lower than the value for 25143 Itokawa. This indicates that the surface lies somewhere between a thick-dust regolith and a rock/boulder/cm-sized, gravel-dominated surface like that of 25143 Itokawa. Our analysis represents the first time that shape and spin-vector information has been derived from a combined data set of visual lightcurves (reflected light) and mid-infrared photometry and spectroscopy (thermal emission).

We investigate the properties and environments of Type Ia Supernova (SN Ia) host galaxies in the Stripe 82 of the Sloan Digital Sky Survey-II Supernova Survey centered on the celestial equator. Host galaxies are defined as the galaxy nearest to the supernova (SN) in terms of angular distance whose velocity difference from the SN is less than 1000 km s^{-1}. Eighty seven SN Ia host galaxies are selected from the SDSS Main galaxy sample with the apparent r-band magnitude m_r < 17.77, and compared with the SDSS Main galaxies. The SN Ia rates for early and late-type galaxies are 0.81 +- 0.19 SN (100yr)^{-1} and 0.99 +- 0.21 SN (100yr)^{-1}, respectively. We find that the host galaxies have a color distribution consistent with that of the Main galaxies, regardless of their morphology. However, host galaxies are on average brighter than the Main galaxies by ~ 0.3 mag over the range of -18.3 > M_r > -21.3. But the brighter ends of their luminosity distributions are similar. The distribution of the distance to the nearest neighbor galaxy shows that SNe Ia are more likely to occur in isolated galaxies without close neighbors. We also find that the SN Ia host galaxies are preferentially located in a region close to massive galaxy clusters compared to the Main galaxies.

We measure the topology of the main galaxy distribution using the Seventh Data Release of the Sloan Digital Sky Survey, examining the dependence of galaxy clustering topology on galaxy properties. The observational results are used to test galaxy formation models. A volume-limited sample defined by $M_r<-20.19$ enables us to measure the genus curve with amplitude of $G=378$ at $6h^{-1}$Mpc smoothing scale, with 4.8\% uncertainty including all systematics and cosmic variance. The clustering topology over the smoothing length interval from 6 to $10 h^{-1}$Mpc reveals a mild scale-dependence for the shift ($\Delta\nu$) and void abundance ($A_V$) parameters of the genus curve. We find substantial bias in the topology of galaxy clustering with respect to the predicted topology of the matter distribution, which varies with luminosity, morphology, color, and the smoothing scale of the density field. The distribution of relatively brighter galaxies shows a greater prevalence of isolated clusters and more percolated voids. Even though early (late)-type galaxies show topology similar to that of red (blue) galaxies, the morphology dependence of topology is not identical to the color dependence. In particular, the void abundance parameter $A_V$ depends on morphology more strongly than on color. We test five galaxy assignment schemes applied to cosmological N-body simulations of a $\Lambda$CDM universe to generate mock galaxies: the Halo-Galaxy one-to-one Correspondence model, the Halo Occupation Distribution model, and three implementations of Semi-Analytic Models (SAMs). None of the models reproduces all aspects of the observed clustering topology; the deviations vary from one model to another but include statistically significant discrepancies in the abundance of isolated voids or isolated clusters and the amplitude and overall shift of the genus curve. (Abridged)

A new catalog of isolated galaxies from The Sloan Digital Sky Survey (DR5) is presented. 1520 isolated galaxies were found in 1.4 steradians of sky. The selection criteria in this so called UNAM-KIAS catalog was implemented from a variation on the criteria developed by Karachentseva 1973 including full redshift information. Through an image processing pipeline that takes advantage from the high resolution (~ 0.4 ''/pix) and high dynamic range of the SDSS images, a uniform g band morphological classification for all these galaxies is presented. We identify 80% (SaSm) spirals (50% later than Sbc types) on one hand, and a scarce population of early-type E(6.5%) and S0(8%) galaxies amounting to 14.5% on the other hand. This magnitude-limited catalog is ~ 80% complete at 16.5, 15.6, 15.0, 14.6 and 14.4 magnitudes in the ugriz bands respectively. Some representative physical properties including SDSS magnitudes and color distributions, color-color diagrams, absolute magnitude-color, and concentration-color diagrams as a function of morphological type are presented. The UNAM-KIAS Morphological Atlas is also released along with this paper. For each galaxy of type later than Sa, a mosaic is presented that includes: (1) a g-band logarithmic image, (2) a g band filtered-enhanced image where a Gaussian kernel of various sizes was applied and (3) an RGB color image from the SDSS database. For E/S0/Sa galaxies, in addition to the images in (1), (2) and (3), plots of r band surface brightness and geometric profiles (ellipticity, Position Angle PA and A4/B4 coefficients of the Fourier series expansions of deviations of a pure ellipse) are provided...

A long infrared jet has been discovered by the Spitzer c2d legacy program in core A of L1251. It is associated with a very embedded Class 0 object with an accretion luminosity of about 0.9 Lsun derived by radiative transfer model fitting to the observed SED. Comparing the observed IRAC colors along the infrared jet with those calculated from a model of an admixture of gas with a power-law temperature distribution indicates that the jet is possibly created by a paraboloidal bow shock propagating into the ambient medium of n(H_2)=10^5 cm^{-3}. In addition, the variation of the power-law index along the jet suggests that the portion of hot gas decreases with distance from the jet engine. The molecular outflow in this region has been mapped for the first time using CO data. From the calculated outflow momentum flux, a very strong lower limit to the average accretion luminosity is 3.6 sin i/cos^3 i Lsun, indicative of a decrease in the accretion rate with time.

We present a study on the environments of the SDSS galaxies divided into fine classes based on their morphology, colour and spectral features. The SDSS galaxies are classified into early-type and late-type; red and blue; passive, HII, Seyfert and LINER, which returns a total of 16 fine classes of galaxies. We estimate the local number density, target-excluded local luminosity density, local colour, close pair fraction and the luminosity and colour of the brightest neighbour, which are compared between the fine classes comprehensively. The morphology-colour class of galaxies strongly depends on the local density, with the approximate order of high-density preference: red early-type galaxies (REGs) -- red late-type galaxies (RLGs) -- blue early-type galaxies (BEGs) -- blue late-type galaxies (BLGs). We find that high-density environments (like cluster environments) seem to suppress AGN activity. The pair fraction of HII REGs does not show statistically significant difference from that of passive REGs, while the pair fraction of HII BLGs is smaller than that of non-HII BLGs. HII BLGs show obvious double (red + blue) peaks in the distribution of the brightest neighbour colour, while red galaxies show a single red peak. The brightest neighbours of Seyfert BLGs tend to be blue, while those of LINER BLGs tend to be red, which implies that the difference between Seyfert and LINER may be related to the pair interaction. Other various environments of the fine classes are investigated, and their implication on galaxy evolution is discussed.

We present a multi-wavelength study of the nature of the SDSS galaxies divided into fine classes based on their morphology, colour and spectral features. The SDSS galaxies are classified into early-type and late-type; red and blue; passive, HII, Seyfert and LINER, which returns a total of 16 fine classes of galaxies. The properties of galaxies in each fine class are investigated from radio to X-ray, using 2MASS, IRAS, FIRST, NVSS, GALEX and ROSAT data. The UV - optical - NIR colours of blue early-type galaxies (BEGs) seem to result from the combination of old stellar population and recent star formation (SF). Non-passive red early-type galaxies (REGs) have larger metallicity and younger age than passive REGs, which implies that non-passive REGs have suffered recent SF adding young and metal-rich stars to them. The radio detection fraction of REGs strongly depends on their optical absolute magnitudes, while that of most late-type galaxies does not, implying the difference in their radio sources: AGN and SF. The UV - optical colours and the radio detection fraction of passive RLGs show that they have properties similar to REGs rather than non-passive RLGs. Dust extinction may not be a dominant factor making RLGs red, because RLGs are detected in the mid- and far-infrared bands less efficiently than blue late-type galaxies (BLGs). The passive BLGs have very blue UV - optical - NIR colours, implying either recent SF quenching or current SF in their outskirts. Including star formation rate, other multi-wavelength properties in each fine class are investigated, and their implication on the identity of each fine class is discussed (abridged).

We search for isolated galaxies using a volume-limited sample of galaxies with 0.02< z < 0.04742 from SDSS DR7 supplemented by bright galaxies. We devise a diagnostic tool to select isolated galaxies in different environments using the projected separation (r_p) normalized by the virial radius of the nearest neighbor (r_{vir,nei}) and the local background density. We find that the isolation condition of r_p>r_{vir,nei} and \rho <\bar{\rho} well segregates the CIG galaxies. We confirm the morphology conformity between the host and their satellites, which suggests importance of hydrodynamic interaction among galaxies within their virial radii in galaxy evolution.

Using a sample of galaxy groups found in the Sloan Digital Sky Survey Data Release 4, we measure the following four types of alignment signals: (1) the alignment between the distributions of the satellites of each group relative to the direction of the nearest neighbor group (NNG); (2) the alignment between the major axis direction of the central galaxy of the host group (HG) and the direction of the NNG; (3) the alignment between the major axes of the central galaxies of the HG and the NNG; and (4) the alignment between the major axes of the satellites of the HG and the direction of the NNG. We find strong signal of alignment between the satellite distribution and the orientation of central galaxy relative to the direction of the NNG, even when the NNG is located beyond $3r_{\rm vir}$ of the host group. The major axis of the central galaxy of the HG is aligned with the direction of the NNG. The alignment signals are more prominent for groups that are more massive and with early type central galaxies. We also find that there is a preference for the two major axis of the central galaxies of the HG and NNG to be parallel for the system with both early central galaxies, however not for the systems with both late type central galaxies. For the orientation of satellite galaxies, we do not find any significant alignment signals relative to the direction of the NNG. From these four types of alignment measurements, we conclude that the large scale environment traced by the nearby group affects primarily the shape of the host dark matter halo, and hence also affects the distribution of satellite galaxies and the orientation of central galaxies. (Shortened version)

We are using the VLBA and the Japanese VERA project to measure trigonometric parallaxes and proper motions of masers found in high-mass star-forming regions across the Milky Way. Early results from 18 sources locate several spiral arms. The Perseus spiral arm has a pitch angle of 16 +/- 3 degrees, which favors four rather than two spiral arms for the Galaxy. Combining positions, distances, proper motions, and radial velocities yields complete 3-dimensional kinematic information. We find that star forming regions on average are orbiting the Galaxy ~15 km/s slower than expected for circular orbits. By fitting the measurements to a model of the Galaxy, we estimate the distance to the Galactic center R_o = 8.4 +/- 0.6 kpc and a circular rotation speed Theta_o = 254 +/- 16 km/s. The ratio Theta_o/R_o can be determined to higher accuracy than either parameter individually, and we find it to be 30.3 +/- 0.9 km/s/kpc, in good agreement with the angular rotation rate determined from the proper motion of Sgr A*. The data favor a rotation curve for the Galaxy that is nearly flat or slightly rising with Galactocentric distance. Kinematic distances are generally too large, sometimes by factors greater than two; they can be brought into better agreement with the trigonometric parallaxes by increasing Theta_o/R_o from the IAU recommended value of 25.9 km/s/kpc to a value near 30 km/s/kpc. We offer a "revised" prescription for calculating kinematic distances and their uncertainties, as well as a new approach for defining Galactic coordinates. Finally, our estimates of Theta_o and To/R_o, when coupled with direct estimates of R_o, provide evidence that the rotation curve of the Milky Way is similar to that of the Andromeda galaxy, suggesting that the dark matter halos of these two dominant Local Group galaxy are comparably massive.

Using a large sample of local galaxies (144,940) with -17.5< Mr <-22 and 0.025< z <0.107, selected from SDSS DR5, we compare AGN host galaxies with non-AGN galaxies at matched luminosity, velocity dispersion, color, color gradient, or concentration index, to investigate how AGN activity is related with galaxy properties. The AGN sample is composed of Type II AGNs identified with flux ratios of narrow-emission lines with S/N > 6. We find that the fraction of galaxies hosting an AGN (f_AGN) depends strongly on morphology together with color, and very weakly on luminosity or velocity dispersion of host galaxies. In particular, f_AGN of early-type galaxies is almost independent of luminosity nor velocity dispersion when color is fixed. The host galaxy color preferred by AGNs is u-r ~2.0 for early-type hosts and u-r=2.0-2.4 for late-type hosts. This trend suggests that AGNs are dominantly hosted by intermediate-mass late-type galaxies. We also investigate how the accretion power varies with galaxy properties. While the Eddington ratio ([OIII] line luminosity normalized by black hole mass) ranges over three orders of magnitude for both morphological types, late-type galaxies are the dominant hosts over all AGN power. Among late-type galaxies, bluer color galaxies host higher power AGNs. These results are consistent with a scenario that more massive and redder galaxies are harder to host AGNs since these galaxiesalready consumed gas at the center or do not have sufficient gas supply to feed the black hole. In contrast, intermediate-mass, intermediate-color, and more concentrated late-type galaxies are more likely to host AGNs, implying that perhaps some fraction of low-mass, blue, and less concentrated lat e-type galaxiesmay not host massive black holes or may host very low-power AGNs.

Strong Balmer absorption lines and the lack of Ha and [OII] emission lines signify that E+As are post-starburst systems. Recent studies suggest that E+As may undergo the transition from the `blue cloud' to the `red sequence' and eventually migrate to red sequence ETGs. An observational validation of this scenario is to identify the intervening galaxy population between E+As and the red-sequence. Motivated by recent findings with GALEX that a large fraction of ETGs exhibit UV-excess as a sign of RSF, we investigate the possible connection of the UV-excess galaxies to E+As. In particular, we examine the FP scaling relations of the largest sample of ~1,000 E+As selected from the SDSS and ~20,000 morphologically-selected SDSS ETGs with GALEX UV data. The FP parameters, combined with stellar population indicators, reveal a certain group of UV-excess ETGs that bridges between E+As and quiescent red galaxies. The newly identified galaxies are the post-starburst systems characterized by UV-excess but no Ha emission. This is a conceptual generalisation of "E+A", in that the Balmer absorption line in the "E+A" definition is replaced with UV-optical colours that are far more sensitive to RSF than the Balmer lines. We refer to these UV-excess galaxies as "E+a" galaxies, which stands for elliptical ("E") galaxies with a minority of A-type ("a") young stars. The species are either (1) galaxies that experienced starbursts weaker than those observed in E+As (1~10% of E+As, "mild E+As") or (2) the products of passively evolved E+As after quenching star formation quite a while ago (~1 Gyr, "old E+As"). We suggest that the latter type of E+a galaxies represents the most recent arrival to the red sequence in the final phase of the "E+A" to "red early-type" transition. (Abridged)

We measure the 3D genus topology of large scale structure using Luminous Red Galaxies (LRGs) in the Sloan Digital Sky Survey and find it consistent with the Gaussian random phase initial conditions expected from the simplest scenarios of inflation. This studies 3D topology on the largest scales ever obtained. The topology is sponge-like. We measure topology in two volume-limited samples: a dense shallow sample studied with smoothing length of 21h^{-1}Mpc, and a sparse deep sample studied with a smoothing length of 34h^{-1}Mpc. The amplitude of the genus curve is measured with 4% uncertainty. Small distortions in the genus curve expected from non-linear biasing and gravitational effects are well explained (to about 1-\sigma accuracy) by N-body simulations using a subhalo-finding technique to locate LRGs. This suggests the formation of LRGs is a clean problem that can be modeled well without any free fitting parameters. This bodes well for using LRGs to measure the characteristic scales such as the baryon oscillation scale in future deep redshift surveys.

We inspect the coupled dependence of physical parameters of the SDSS galaxies on the small-scale (distance to and morphology of the nearest neighbor galaxy) and the large-scale (background density smoothed over 20 nearby galaxies) environments. The impacts of interaction on galaxy properties are detected at least out to the neighbor separation corresponding to the virial radius of galaxies, which is typically between 200 and 400 h^{-1}kpc for the galaxies in our sample. To detect these long-range interaction effects it is crucial to divide galaxy interactions into four cases dividing the morphology of target and neighbor galaxies into early and late types. We show that there are two characteristic neighbor-separation scales where the galaxy interactions cause abrupt changes in the properties of galaxies. The first scale is the virial radius of the nearest neighbor galaxy r_{vir,nei}. Many physical parameters start to deviate from those of extremely isolated galaxies at the projected neighbor separation r_p of about r_{vir,nei}. The second scale is at r_p \approx 0.05 r_{vir,nei} = 10--20 h^{-1} kpc, and is the scale at which the galaxies in pairs start to merge. We find that late-type neighbors enhance the star formation activity of galaxies while early-type neighbors reduce it, and that these effects occur within r_{vir,nei}. The hot halo gas and cold disk gas must be participating in the interactions at separations less than the virial radius of the galaxy plus dark halo system. Our results also show that the role of the large-scale density in determining galaxy properties is minimal once luminosity and morphology are fixed. We propose that the weak residual dependence of galaxy properties on the large-scale density is due to the dependence of the halo gas property on the large-scale density.

We present a distance measurement for the semiregular variable S Crateris (S Crt) based on its annual parallax. With the unique dual beam system of the VLBI Exploration for Radio Astrometry (VERA) telescopes, we measured the absolute proper motion of a water maser spot associated with S Crt, referred to the quasar J1147-0724 located at an angular separation of 1.23$^{\circ}$. In observations spanning nearly two years, we have detected the maser spot at the LSR velocity of 34.7 km s$^{-1}$, for which we measured the annual parallax of 2.33$\pm$0.13 mas corresponding to a distance of 430$^{+25}_{-23}$ pc. This measurement has an accuracy one order of magnitude better than the parallax measurements of HIPPARCOS. The angular distribution and three-dimensional velocity field of maser spots indicate a bipolar outflow with the flow axis along northeast-southwest direction. Using the distance and photospheric temperature, we estimate the stellar radius of S Crt and compare it with those of Mira variables.

We present results of multi-epoch VLBI observations with VERA (VLBI Exploration of Radio Astrometry) of the 22 GHz H$_{2}$O masers associated with a young stellar object (YSO) IRAS 22198+6336 in a dark cloud L1204G. Based on the phase-referencing VLBI astrometry, we derive an annual parallax of IRAS 22198+6336 to be 1.309$\pm$0.047 mas, corresponding to the distance of 764$\pm$27 pc from the Sun. Although the most principal error source of our astrometry is attributed to the internal structure of the maser spots, we successfully reduce the errors in the derived annual parallax by employing the position measurements for all of the 26 detected maser spots. Based on this result, we reanalyze the spectral energy distribution (SED) of IRAS 22198+6336 and find that the bolometric luminosity and total mass of IRAS 22198+6336 are 450$L_{\odot}$ and 7$M_{\odot}$, respectively. These values are consistent with an intermediate-mass YSO deeply embedded in the dense dust core, which has been proposed to be an intermediate-mass counterpart of a low-mass Class 0 source. In addition, we obtain absolute proper motions of the H$_{2}$O masers for the most blue-shifted components. We propose that the collimated jets aligned along the east-west direction are the most plausible explanation for the origin of the detected maser features.

We report astrometric observations of H2O masers around the red supergiant VY Canis Majoris (VY CMa) carried out with VLBI Exploration of Radio Astrometry (VERA). Based on astrometric monitoring for 13 months, we successfully measured a trigonometric parallax of 0.88 +/- 0.08 mas, corresponding to a distance of 1.14 +0.11/-0.09 kpc. This is the most accurate distance to VY CMa and the first one based on an annual parallax measurement. The luminosity of VY CMa has been overestimated due to a previously accepted distance. With our result, we re-estimate the luminosity of VY CMa to be (3 +/- 0.5) x 10^5 L_sun using the bolometric flux integrated over optical and IR wavelengths. This improved luminosity value makes location of VY CMa on the Hertzsprung-Russel (HR) diagram much closer to the theoretically allowable zone (i.e. the left side of the Hayashi track) than previous ones, though uncertainty in the effective temperature of the stellar surface still does not permit us to make a final conclusion.

We present a comprehensive study of the nature of the SDSS galaxies divided into various classes based on their morphology, colour, and spectral features. The SDSS galaxies are classified into early-type and late-type; red and blue; passive, HII, Seyfert, and LINER, returning a total of 16 fine classes of galaxies. We examine the luminosity dependence of seven physical parameters of galaxies in each class. We find that more than half of red early-type galaxies (REGs) have star formation or AGN activity, and that these active REGs have smaller axis ratio and bluer outside compared to the passive REGs. Blue early-type galaxies (BEGs) show structural features similar to those of REGs, but their centres are bluer than REGs. HII BEGs are found to have bluer centres than passive BEGs, but HII REGs have bluer outside than passive REGs. Bulge-dominated late-type galaxies have red colours. Passive red late-types are similar to REGs in several aspects. Most blue late-type galaxies (BLGs) have forming stars, but a small fraction of BLGs do not show evidence for current star formation activity. Differences of other physical parameters among different classes are inspected, and their implication on galaxy evolution is discussed.

We have used the Japanese VLBI array VERA to perform high-precision astrometry of an H2O maser source in the Galactic star-forming region NGC 281 West, which has been considered to be part of a 300-pc superbubble. We successfully detected a trigonometric parallax of 0.355+/-0.030 mas, corresponding to a source distance of 2.82+/-0.24 kpc. Our direct distance determination of NGC 281 has resolved the large distance discrepancy between previous photometric and kinematic studies; likely NGC 281 is in the far side of the Perseus spiral arm. The source distance as well as the absolute proper motions were used to demonstrate the 3D structure and expansion of the NGC 281 superbubble, ~650 pc in size parallel to the Galactic disk and with a shape slightly elongated along the disk or spherical, but not vertically elongated, indicating the superbubble expansion may be confined to the disk. We estimate the expansion velocity of the superbubble as ~20 km/s both perpendicular to and parallel to the Galactic disk with a consistent timescale of ~20 Myr.

We have studied the radial distribution of the early (E/S0) and late (S/Irr) types of satellites around bright host galaxies. We made a volume-limited sample of 4,986 satellites brighter than M_r = -18.0 associated with 2,254 hosts brighter than M_r =-19.0 from the SDSS DR5 sample. The morphology of satellites is determined by an automated morphology classifier, but the host galaxies are visually classified. We found segregation of satellite morphology as a function of the projected distance from the host galaxy. The amplitude and shape of the early-type satellite fraction profile are found to depend on the host luminosity. This is the morphology-radius/density relation at the galactic scale. There is a strong tendency for morphology conformity between the host galaxy and its satellites. The early-type fraction of satellites hosted by early-type galaxies is systematically larger than that of late-type hosts, and is a strong function of the distance from the host galaxies. Fainter satellites are more vulnerable to the morphology transformation effects of hosts. Dependence of satellite morphology on the large-scale background density was detected. The fraction of early-type satellites increases in high density regions for both early and late-type hosts. It is argued that the conformity in morphology of galactic satellite system is mainly originated by the hydrodynamical and radiative effects of hosts on satellites.

We propose a model of allocating galaxies in cosmological N-body simulations. we identify each subhalo with a galaxy, and assign luminosity and morphological type assuming that the galaxy luminosity is a monotonic function of its host subhalo mass. The morphology assignment is made by using two simple relations between subhalo mass and galaxy luminosity of different types. One is using a constant ratio in luminosity of early (E/SO) and late (S/Irr) type galaxies at a fixed subhalo mass. And the other assumes that galaxies of different morphological types but having an equal luminosity have a constant ratio in their subhalo masses. We made a series of comparisons of the properties of these simulated galaxies with those of the SDSS galaxies. The resulting simulated galaxy sample is found to successfully reproduce the observed local number density distribution except for in high density regions. The luminosity function is studied as a function of local density. It was found that the observed luminosity functions in different local density environments are overall well-reproduced by the simulated galaxies. Discrepancy is found at the bright end of the luminosity function of early types in the underdense regions and at the faint end of both morphological types in very high density regions. A significant fraction of the observed early type galaxies in voids seems to have undergone a relatively recent star formation and became brighter. The lack of faint simulated galaxies in dense regions may be due to the strong tidal force of the central halo which destroys less massive satellite subhalos around in the simulation. The mass-to-light ratio is found to depend on the local density in the way similar to that observed in the SDSS sample. (abridged)

We present a unified picture on the evolution of galaxy luminosity and morphology. Galaxy morphology is found to depend critically on the local environment set up by the nearest neighbor galaxy in addition to luminosity and the large scale density. When a galaxy is located farther than the virial radius from its closest neighbor, the probability for the galaxy to have an early morphological type is an increasing function only of luminosity and the local density due to the nearest neighbor ($\rho_n$). The tide produced by the nearest neighbor is thought to be responsible for the morphology transformation toward the early type at these separations. When the separation is less than the virial radius, i.e. when $\rho_n > \rho_{\rm virial}$, its morphology depends also on the neighbor's morphology and the large-scale background density over a few Mpc scales ($\rho_{20}$) in addition to luminosity and $\rho_n$. The early type probability keeps increasing as $\rho_n$ increases if its neighbor is an early type. But the probability decreases as $\rho_n$ increases when the neighbor is a late type. The cold gas streaming from the late type neighbor can be the reason for the morphology transformation toward late type. The overall early-type fraction increases as $\rho_{20}$ increases when $\rho_n > \rho_{\rm virial}$. This can be attributed to the hot halo gas of the neighbor which is confined by the pressure of the ambient medium held by the background mass. We have also found that galaxy luminosity depends on $\rho_n$, and that the isolated bright galaxies are more likely to be recent merger products. We propose a scenario that a series of morphology and luminosity transformation occur through distant interactions and mergers, which results in the morphology--luminosity--local density relation.

We report on the results of multi-epoch VLBI observations with VERA (VLBI Exploration of Radio Astrometry) of the 22 GHz H2O masers associated with the young stellar object SVS 13 in the NGC 1333 region. We have carried out phase-referencing VLBI astrometry and measured an annual parallax of the maser features in SVS 13 of 4.25+/-0.32 mas, corresponding to the distance of 235+/-18 pc from the Sun. Our result is consistent with a photometric distance of 220 pc previously reported. Even though the maser features were detectable only for 6 months, the present results provide the distance to NGC 1333 with much higher accuracy than photometric methods. The absolute positions and proper motions have been derived, revealing that the H2O masers with the LSR (local standard of rest) velocities of 7-8 km s-1 are most likely associated with VLA4A, which is a radio counterpart of SVS 13. The origin of the observed proper motions of the maser features are currently difficult to attribute to either the jet or the rotating circumstellar disk associated with VLA4A, which should be investigated through future high-resolution astrometric observations of VLA4A and other radio sources in NGC 1333.

We have performed high-precision astrometry of H2O maser sources in Galactic star forming region Sharpless 269 (S269) with VERA. We have successfully detected a trigonometric parallax of 189+/-8 micro-arcsec, corresponding to the source distance of 5.28 +0.24/-0.22 kpc. This is the smallest parallax ever measured, and the first one detected beyond 5 kpc. The source distance as well as proper motions are used to constrain the outer rotation curve of the Galaxy, demonstrating that the difference of rotation velocities at the Sun and at S269 (which is 13.1 kpc away from the Galaxy's center) is less than 3%. This gives the strongest constraint on the flatness of the outer rotation curve and provides a direct confirmation on the existence of large amount of dark matter in the Galaxy's outer disk.

We present a study of the connection between star-forming galaxies, AGN host galaxies, and normal early-type galaxies in the Sloan Digital Sky Survey (SDSS). Using the SDSS DR5 and DR4plus data, we select our early-type galaxy sample in the color versus color-gradient space, and we classify the spectral types of the selected early-type galaxies into normal, star-forming, Seyfert, and LINER classes, using several spectral line flux ratios. We investigate the slope in the fundamental space for each class of early-type galaxies and find that there are obvious differences in the slopes of the fundamental planes (FPs) among the different classes of early-type galaxies, in the sense that the slopes for Seyferts and star-forming galaxies are flatter than those for normal galaxies and LINERs. This may be the first identification of the systematic variation of the FP slope among the subclasses of early-type galaxies. The difference in the FP slope might be caused by the difference in the degree of nonhomology among different classes or by the difference of gas contents in their merging progenitors. One possible scenario is that the AGN host galaxies and star-forming galaxies are formed by gas-rich merging and that they may evolve into normal early-type galaxies after finishing their star formation or AGN activities.

We present the initial results of multi-epoch VLBI observations of the 22 GHz H2O masers in the Orion KL region with VERA (VLBI Exploration of Radio Astrometry). With the VERA dual-beam receiving system, we have carried out phase-referencing VLBI astrometry and successfully detected an annual parallax of Orion KL to be 2.29+/-0.10 mas, corresponding to the distance of 437+/-19 pc from the Sun. The distance to Orion KL is determined for the first time with the annual parallax method in these observations. Although this value is consistent with that of the previously reported, 480+/-80 pc, which is estimated from the statistical parallax method using proper motions and radial velocities of the H2O maser features, our new results provide the much more accurate value with an uncertainty of only 4%. In addition to the annual parallax, we have detected an absolute proper motion of the maser feature, suggesting an outflow motion powered by the radio source I along with the systematic motion of source I itself.

We report on absolute proper-motion measurements of an H2O maser source in the NGC 281 West molecular cloud, which is located ~320 pc above the Galactic plane and is associated with an HI loop extending from the Galactic plane. We have conducted multi-epoch phase-referencing observations of the maser source with VERA (VLBI Exploration of Radio Astrometry) over a monitoring period of 6 months since May 2006. We find that the H2O maser features in NGC 281 West are systematically moving toward the southwest and further away from the Galactic plane with a vertical velocity of ~20-30 km/s at its estimated distance of 2.2-3.5 kpc. Our new results provide the most direct evidence that the gas in the NGC 281 region on the HI loop was blown out from the Galactic plane, most likely in a superbubble driven by multiple or sequential supernova explosions in the Galactic plane.

We investigate the dependence of physical properties of galaxies brighter than M_r=-18.0 in the SDSS on environment, as measured by local density using an adaptive smoothing kernel. We find that variations of galaxy properties with environment are almost entirely due to the dependence of morphology and luminosity on environment. Because galaxy properties depend not only on luminosity but also on morphology, it is clear that galaxy properties cannot be determined solely by dark halo mass. When morphology and luminosity are fixed, other physical properties, such as color, color-gradient, concentration, size, velocity dispersion, and star formation rate, are nearly independent of local density. The only feature is the sharp decrease of late type fraction above the critical luminosity of about M_r=-21.3 in the morphology versus luminosity relation. Weak residual dependences on environment include that of the color of late-types (bluer at lower density) and of the L-sigma relation of early-types (larger dispersion at higher density for bright galaxies). The early-type fraction is a monotonically increasing function of local density and luminosity. The morphology-density- luminosity relation should be a key constraint on galaxy formation models. We demonstrate that the dependence on environment of the morphology of galaxie originates from variations in density on effective Gaussian smoothing scales much smaller than 12 h^{-1}Mpc. We find that galaxy morphology varies both with density measured on an effective Gaussian smoothing scale of 4.7 h^{-1} Mpc and on distance to the nearest bright galaxy (particularly, for about 0.2 h^{-1}Mpc). We propose a mechanism that the morphology of galaxies in galaxy systems is transformed by the tidal force.

We examine volume-limited samples from the SDSS galaxies to look for relations among internal and collective physical parameters of galaxies as faint as M_r=- 17.5, which include morphology, luminosity, color, color gradient, concentration, size, velocity dispersion, equivalent width (EW) of H_alpha line,axis ratio, the luminosity and velocity dispersion functions. At fixed morphology and luminosity, we find that bright (M_r<-20) early-types show very small dispersions in color, color gradient, concentration, size, and velocity dispersion. These dispersions increase at fainter magnitudes, where the fraction of blue early-types increases. Late-types show wider dispersions in all physical parameters compared to early types at the same luminosity. Concentration indices of early-types are well-correlated with velocity dispersion, but are insensitive to luminosity and color for bright galaxies. The slope of the Faber-Jackson relation continuously changes from 4.6 +- 0.4 to 2.7+- 0.2 when luminosity changes from M_r = -22 to -20. The size of early- types is well-correlated with stellar velocity dispersion (for >100 km/s). We find that passive spirals are well-separated from star-forming late-types at EW (H_alpha) of about 4. An interesting finding is that many physical parameters of galaxies manifest different behaviors across the absolute magnitude of about M_* +- 1. The morphology fraction as a function of luminosity depends less sensitively on large scale structure than the luminosity function (LF) does, and thus seems to be more universal. The effects of internal extinction in late-types on the completeness of volume limited samples and on the LF and morphology fraction are found to be very important. (abridged)

Asteroid 2867 Steins is one of two flyby targets of ESAs international Rosetta mission, launched in March, 2004. We obtained photometric observations of Steins on April 14-16, 2004 at Table Mountain Observatory, California, in order to characterize the asteroid physically, information that is crucial for planning the Steins flyby. This study includes the first detailed analysis of the physical properties of Steins from time-series R-filter data along with V- and I-filter photometric measurements.

Using simple dimensional arguments for both spiral and elliptical galaxies, we present formulas to derive an estimate of the halo spin parameter $\lambda$ for any real galaxy, in terms of common observational parameters. This allows a rough estimate of $\lambda$, which we apply to a large volume limited sample of galaxies taken from the SDSS data base. The large numbers involved (11,597) allow the derivation of reliable $\lambda$ distributions, as signal adds up significantly in spite of the errors in the inferences for particular galaxies. We find that if the observed distribution of $\lambda$ is modeled with a log-normal function, as often done for this distribution in dark matter halos that appear in cosmological simulations, we obtain parameters $\lambda_{0}=0.04 \pm 0.005$ and $\sigma_{\lambda}=0.51 \pm 0.05$, interestingly consistent with values derived from simulations. For spirals, we find a good correlation between empirical values of $\lambda$ and visually assigned Hubble types, highlighting the potential of this physical parameter as an objective classification tool.

We measure the three-dimensional topology of large-scale structure in the Sloan Digital Sky Survey (SDSS). This allows the genus statistic to be measured with unprecedented statistical accuracy. The sample size is now sufficiently large to allow the topology to be an important tool for testing galaxy formation models. For comparison, we make mock SDSS samples using several state-of-the-art N-body simulations: the Millennium run of Springel et al. (2005)(10 billion particles), Kim & Park (2006) CDM models (1.1 billion particles), and Cen & Ostriker (2006) hydrodynamic code models (8.6 billion cell hydro mesh). Each of these simulations uses a different method for modeling galaxy formation. The SDSS data show a genus curve that is broadly characteristic of that produced by Gaussian random phase initial conditions. Thus the data strongly support the standard model of inflation where Gaussian random phase initial conditions are produced by random quantum fluctuations in the early universe. But on top of this general shape there are measurable differences produced by non-linear gravitational effects (cf. Matsubara 1994), and biasing connected with galaxy formation. The N-body simulations have been tuned to reproduce the power spectrum and multiplicity function but not topology, so topology is an acid test for these models. The data show a ``meatball'' shift (only partly due to the Sloan Great Wall of Galaxies; this shift also appears in a sub-sample not containing the Wall) which differs at the 2.5\sigma level from the results of the Millennium run and the Kim & Park dark halo models, even including the effects of cosmic variance.

We have found the u -r color versus g -i color gradient space can be used for highly successful morphology classification of galaxies in the Sloan Digital Sky Survey. In this space galaxies form early and late type branches well-separated from each other. The location of galaxies along the branches reflects the degree and locality of star formation activity, and monotonically corresponds to the sequence of morphological subclasses. When the concentration index is used together, the completeness and reliability of classification reaches about 91% for a training set of SDSS galaxies brighter than r~15.9. At faintest magnitudes (r~17.5) of the SDSS spectroscopic sample, the performance still remains at about 88%. The new classification scheme will help us find accurate relations of galaxy morphology with spatial and temporal environments, and help one understand the origin of morphology of galaxies.

We measure the topology of volume-limited galaxy samples selected from a parent sample of 314,050 galaxies in the Sloan Digital Sky Survey(SDSS), which is now complete enough to describe the fully three-dimensional topology and its dependence on galaxy properties. We compare the observed genus statistic G\nu_f) to predictions for a Gaussian random field and to the genus measured for mock surveys constructed from new large-volume simulations of the LambdaCDM cosmology. In this analysis we carefully examine the dependence of the observed genus statistic on the Gaussian smoothing scale R_G from 3.5 to 11 h^-1Mpc and on the luminosity of galaxies over the range -22.50<M_r<-18.5. The void multiplicity A_V is less than unity at all smoothing scales. Because A_V cannot become less than 1 through gravitational evolution, this result provides strong evidence for biased galaxy formation in low density environments. We also find clear evidence of luminosity bias of topology within the volume-limited sub-samples. The shift parameter Delta nu indicates that the genus of brighter galaxies shows a negative shift toward a ``meatball'' (i.e. cluster-dominated) topology, while faint galaxies show a positive shift toward a ``bubble'' (i.e. void-dominated) topology. The transition from negative to positive shift occurs approximately at the characteristic absolute magnitude M_{r*}=-20.4. Even in this analysis of the largest galaxy sample to date, we detect the influence of individual large-scale structures, as the shift parameter Delta nu and cluster multiplicity A_C reflect (at ~3 sigma) the presence of the Sloan Great Wall and a X-shaped structure which runs for several hundred Mpc across the survey volume.

We report on the multi-epoch observations of H2O maser emission in star forming region OH 43.8-0.1 carried out with VLBI Exploration of Radio Astrometry (VERA). The large-scale maser distributions obtained by single-beam VLBI mapping reveal new maser spots scattered in area of 0.7 x 1.0 arcsec, in addition to a `shell-like' structure with a scale of 0.3 x 0.5 arcsec which was previously mapped by Downes et al.(1979). Proper motions are also obtained for 43 spots based on 5-epoch monitoring with a time span of 281 days. The distributions of proper motions show a systematic outflow in the north-south direction with an expansion velocity of ~8 km/s, and overall distributions of maser spots as well as proper motions are better represented by a bipolar flow plus a central maser cluster with a complex structure, rather than a shell with a uniform expansion such as those found in Cep A R5 and W75N VLA2. The distance to OH 43.8-0.1 is also estimated based on the statistical parallax, yielding D = 2.8 +/- 0.5 kpc. This distance is consistent with a near kinematic distance and rules out a far kinematics distance (~9 kpc), and the LSR velocity of OH 43.8-0.1 combined with the distance provides a constraint on the flatness of the galactic rotation curve, that there is no systematic difference in rotation speeds at the Sun and at the position of OH 43.8-0.1, which is located at the galacto-centric radius of ~6.3 kpc.

We present the results of indirect searches for Weakly Interacting Massive Particles (WIMPs) with 1679.6 live days of data from the Super-Kamiokande detector using neutrino-induced upward through-going muons. The search is performed by looking for an excess of high energy muon neutrinos from WIMP annihilations in the Sun, the core of the Earth, and the Galactic Center, as compared to the number expected from the atmospheric neutrino background. No statistically significant excess was seen. We calculate flux limits in various angular cones around each of the above celestial objects. We obtain conservative model-independent upper limits on WIMP-nucleon cross-section as a function of WIMP mass and compare these results with the corresponding results from direct dark matter detection experiments.

We report on a strong outburst of the W49N H2O maser observed with VERA. Single-dish monitoring with VERA 20 m telescopes detected a strong outburst of the maser feature at V_LSR = -30.7 km/s in 2003 October. The outburst had a duration of ~100 days and a peak intensity of 7.9 x 10^4 Jy, being one of the strongest outbursts in W49N observed so far. VLBI observations with the VERA array were also carried out near to the maximum phase of the outburst, and the outburst spot was identified in the VLBI map. While the map was in good agreement with previous studies, showing three major concentrations of maser spots, we found a newly formed arc-like structure in the central maser concentration, which may be a shock front powered by a forming star or a star cluster. The outburst spot was found to be located on the arc-like structure, indicating a possible connection of the present outburst to a shock phenomenon.

(abridged) We present a detailed reanalysis of the Chandra data for the galaxy cluster Abell 4059 and its central radio galaxy, PKS2354-35. We also present new 1.4GHz and 4.7GHz CnB-array radio data from the Very Large Array (VLA), as well as a short archival WFPC2 image from the Hubble Space Telescope. The presence of a strong interaction between this radio galaxy and the intracluster medium (ICM) was suggested by Huang & Sarazin (1998) on the basis of a short observation by the High Resolution Imager on ROSAT, and confirmed in our preliminary analysis of the Chandra/ACIS-S data. In particular, X-ray imaging clearly shows two cavities within the ICM that are approximately aligned with the radio-galaxy axis. However, using our new radio maps we fail to find a detailed correspondence between the 1 arcmin scale radio-lobes and the ICM cavities. This suggests that the cavities are ``ghosts'' of a previous burst of powerful activity by PKS 2354-35. We also examine the nature of the central asymmetric ridge (or bar) of X-ray emission extending for 30kpc south-west of the cluster center that has been noted in these previous analyzes. We find the ridge to be denser and cooler than, but probably in pressure balance with, its surroundings. The thermal evolution of this structure seems to be dominated by radiative cooling, possibly enhanced by the radio-galaxy ICM interaction. We discuss several possible models for the formation of this SW ridge and find none of them to be entirely satisfactory.

We present Chandra observations of the galaxy cluster A4059. We find strong evidence that the FR-I radio galaxy PKS 2354-35 at the center of A4059 is inflating cavities with radii ~20 kpc in the intracluster medium, similar to the situation seen in Perseus A and Hydra A. We also find evidence for interaction between the intracluster medium and PKS 2354-35 on small scales in the very center of the cluster. Arguments are presented suggesting that this radio galaxy has faded significantly in radio power (possibly from an FR-II state) over the past 10^8 yrs.

We have studied topology of the distribution of the high redshift galaxies identified in the Hubble Deep Field (HDF) North and South. The two-dimensional genus is measured from the projected distributions of the HDF galaxies at angular scales from $3.8''$ to $ 6.1''$. We have also divided the samples into three redshift slices with roughly equal number of galaxies using photometric redshifts to see possible evolutionary effects on the topology. The genus curve of the HDF North clearly indicates clustering of galaxies over the Poisson distribution while the clustering is somewhat weaker in the HDF South. This clustering is mainly due to the nearer galaxies in the samples. We have also found that the genus curve of galaxies in the HDF is consistent with the Gaussian random phase distribution with no significant redshift dependence.

We investigate a large number of cosmological evolution models for QSOs and Active Galactic Nuclei (AGN). We introduce a spectrum-luminosity correlation as a new input parameter and adopt the estimated mass function (MF) of massive black holes in centers of nearby galactic nuclei as a constraint to distinguish among different QSO/AGN models. We explore three basic types of phenomenological scenarios; (i) Models with multiple short-lived ($\sim$ a few $10^{6-8}$ yrs) populations, (ii) Models with a single long-lived ($\sim 10^9$ yr) QSO population, and (iii) Models with recurrent QSO/AGN activities which are driven by long-term variabilities of the disk instability type. In each model, we derive the expected theoretical luminosity function (LF) and the MF of black holes which grow through mass accretion. We assess the plausibility of each model based on whether each model's LF and MF are compatible with the observed data. We find that the best fits to the observed LFs are obtained in the model with multiple short-lived populations and without any significant spectral evolution. This finding suggests that the QSO populations may be composed of many short-lived generations ($\sim$ a few 10$^8$ yrs) and that there is no significant spectral evolution within each generation. On the other hand, we also show that there is no satisfactory model which can simultaneously account for the observed LF and the estimated MF. We speculate that some of the present-day black holes (BHs) found in galactic nuclei may have formed without undergoing the QSO/AGN phase.

We apply the observed optical/X-ray spectral states of the Galactic black hole candidates (GBHCs) to the cosmological QSO luminosity evolution under the assumptions that QSOs and GBHCs are powered by similar accretion processes and that their emission mechanisms are also similar. The QSO luminosity function (LF) evolution in various energy bands is strongly affected by the spectral evolution which is tightly correlated with the luminosity evolution. We generate a random sample of QSOs born nearly synchronously by allowing the QSOs to have redshifts in a narrow range around an initial high redshift, black hole masses according to a power-law, and mass accretion rates near Eddington rates. The QSOs evolve as a single long-lived population on the cosmological time scale. The pure luminosity evolution results in distinct luminosity evolution features due to the strong spectral evolution. Most notably, different energy bands (optical/UV, soft X-ray, and hard X-ray) show different evolutionary trends and the hard X-ray LF in particular shows an apparent reversal of the luminosity evolution (from decreasing to increasing luminosity) at low redshifts, which is not seen in the conventional pure luminosity evolution scenario without spectral evolution. The resulting mass function of black holes (BHs), which is qualitatively consistent with the observed QSO LF evolution, shows that QSO remnants are likely to be found as BHs with masses in the range 10**8-5x10**10 solar masses. The long-lived single population of QSOs are expected to leave their remnants as supermassive BHs residing in rare, giant elliptical galaxies.