Planetary nebulae (PNe) have diverse morphological shapes, including point-symmetric and multipolar structures. Many PNe also have complicated internal structures such as torus, lobes, knots, and ansae. A complete accounting of all the morphological structures through physical models is difficult. A first step toward such an understanding is to derive the true three-dimensional structure of the nebulae. In this paper, we show that a multipolar nebula with three pairs of lobes can explain many of such features, if orientation and sensitivity effects are taken into account. Using only six parameters - the inclination and position angles of each pair - we are able to simulate the observed images of 20 PNe with complex structures. We suggest that the multipolar structure is an intrinsic structure of PNe and the statistics of multipolar PNe has been severely underestimated in the past.

In this Letter, we report detections of SiO v=3 J=1--0 maser emission in very long baseline interferometric (VLBI) observations towards 4 out of 12 long-period variable stars: WX Psc, R Leo, W Hya, and T Cep. The detections towards WX Psc and T Cep are new ones. We also present successful astrometric observations of SiO v=2 and v=3 J=1--0 maser emissions associated with two stars: WX Psc and W Hya and their position-reference continuum sources: J010746.0+131205 and J135146.8-291218 with the VLBI Exploration of Radio Astrometry (VERA). The relative coordinates of the position-reference continuum source and SiO v=3 maser spots were measured with respect to those of an SiO v=2 maser spot adopted as fringe-phase reference. Thus the faint continuum sources were inversely phase-referenced to the bright maser sources. It implies possible registration of multiple SiO maser line maps onto a common coordinate system with 10 microarcsecond-level accuracy.

We observed four "water fountain" sources in the CO J=3-2 line emission with the Atacama Submillimeter Telescope Experiment (ASTE) 10 m telescope in 2010-2011. The water fountain sources are evolved stars that form high-velocity collimated jets traced by water maser emission. The CO line was detected only from IRAS 16342-3814. The present work confirmed that the ^{12}CO to ^{13}CO line intensity ratio is ~1.5 at the systemic velocity. We discuss the origins of the very low ^{12}CO to ^{13}CO intensity ratio, as possible evidence for the "hot-bottom burning" in an oxygen-rich star, and the CO intensity variation in IRAS 16342-3814.

We present the results of multi-epoch very long baseline interferometry (VLBI) water (H2O) maser observations carried out with the VLBI Exploration of Radio Astrometry (VERA) toward the Cepheus A HW3d object. We measured for the first time relative proper motions of the H2O maser features, whose spatio-kinematics traces a compact bipolar outflow. This outflow looks highly collimated and expanding through ~ 280 AU (400 mas) at a mean velocity of ~ 21 km/s (~ 6 mas/yr) without taking into account the turbulent central maser cluster. The opening angle of the outflow is estimated to be ~ 30{\circ}. The dynamical time-scale of the outflow is estimated to be ~ 100 years. Our results provide strong support that HW3d harbors an internal massive young star, and the observed outflow could be tracing a very early phase of star formation. We also have analyzed Very Large Array (VLA) archive data of 1.3 cm continuum emission obtained in 1995 and 2006 toward Cepheus A. The comparative result of the HW3d continuum emission suggests the possibility of the existence of distinct young stellar objects (YSOs) in HW3d and/or strong variability in one of their radio continuum emission components.

We present the result of radio observations of red supergiants in the star cluster, Stephenson's #2, and candidates for red supergiants in the star clusters, Mercer et al. (2005)'s #4, #8, and #13, in the SiO and H$_2$O maser lines.The Stephenson's #2 cluster and nearby aggregation at the South-West contain more than 15 red supergiants. We detected one at the center of Stephenson's #2 and three in the south-west aggregation in the SiO maser line, and three of these 4 were also detected in the H2O maser line. The average radial velocity of the 4 detected objects is 96 km s^{-1}, giving a kinematic distance of 5.5 kpc, which locates this cluster near the base of the Scutum-Crux spiral arm. We also detected 6 SiO emitting objects associated with the other star clusters. In addition, mapping observations in the CO J=1--0 line toward these clusters revealed that an appreciable amount of molecular gas still remains around Stephenson's #2 cluster in contrast to the prototypical red-supergiant cluster, Bica et al.'s #122. It indicates that a time scale of gas expulsion differs considerably in individual clusters.

We present a catalog of 100,563 unresolved, UV-excess (UVX) quasar candidates to g=21 from 2099 deg^2 of the Sloan Digital Sky Survey (SDSS) Data Release One (DR1) imaging data. Existing spectra of 22,737 sources reveals that 22,191 (97.6%) are quasars; accounting for the magnitude dependence of this efficiency, we estimate that 95,502 (95.0%) of the objects in the catalog are quasars. Such a high efficiency is unprecedented in broad-band surveys of quasars. This ``proof-of-concept'' sample is designed to be maximally efficient, but still has 94.7% completeness to unresolved, g<~19.5, UVX quasars from the DR1 quasar catalog. This efficient and complete selection is the result of our application of a probability density type analysis to training sets that describe the 4-D color distribution of stars and spectroscopically confirmed quasars in the SDSS. Specifically, we use a non-parametric Bayesian classification, based on kernel density estimation, to parameterize the color distribution of astronomical sources -- allowing for fast and robust classification. We further supplement the catalog by providing photometric redshifts and matches to FIRST/VLA, ROSAT, and USNO-B sources. Future work needed to extend the this selection algorithm to larger redshifts, fainter magnitudes, and resolved sources is discussed. Finally, we examine some science applications of the catalog, particularly a tentative quasar number counts distribution covering the largest range in magnitude (14.2<g<21.0) ever made within the framework of a single quasar survey.

We investigate the evolution of the fractions of late type cluster galaxies as a function of redshift, using one of the largest, most uniform cluster samples available. The sample consists of 514 clusters of galaxies in the range 0.02<z<0.3 from the Sloan Digital Sky Survey Cut & Enhance galaxy cluster catalog. This catalog was created using a single automated cluster finding algorithm on uniform data from a single telescope, with accurate CCD photometry, thus, minimizing selection biases. We use four independent methods to analyze the evolution of the late type galaxy fraction. Specifically, we select late type galaxies based on: restframe g-r color, u-r color, galaxy profile fitting and concentration index. The first criterion corresponds to the one used in the classical Butcher-Oemler analyses. The last three criteria are more sensitive to the morphological type of the galaxies. In all four cases, we find an increase in the fraction of late type galaxies with increasing redshift, significant at the 99.9% level. The results confirm that cluster galaxies do change colors with redshift (the Butcher-Oemler effect) and, in addition, they change their morphology to later-type toward higher redshift -- indicating a morphological equivalent of the Butcher-Oemler effect. We also find a tendency of richer clusters to have lower fractions of late type galaxies. The trend is consistent with a ram pressure stripping model, where richer clusters have more effective ram pressure due to their higher temperature.

I present here a review of past and present multi-disciplinary research of the Pittsburgh Computational AstroStatistics (PiCA) group. This group is dedicated to developing fast and efficient statistical algorithms for analysing huge astronomical data sources. I begin with a short review of multi-resolutional kd-trees which are the building blocks for many of our algorithms. For example, quick range queries and fast n-point correlation functions. I will present new results from the use of Mixture Models (Connolly et al. 2000) in density estimation of multi-color data from the Sloan Digital Sky Survey (SDSS). Specifically, the selection of quasars and the automated identification of X-ray sources. I will also present a brief overview of the False Discovery Rate (FDR) procedure (Miller et al. 2001a) and show how it has been used in the detection of ``Baryon Wiggles'' in the local galaxy power spectrum and source identification in radio data. Finally, I will look forward to new research on an automated Bayes Network anomaly detector and the possible use of the Locally Linear Embedding algorithm (LLE; Roweis & Saul 2000) for spectral classification of SDSS spectra.

We outline here the next generation of cluster-finding algorithms. We show how advances in Computer Science and Statistics have helped develop robust, fast algorithms for finding clusters of galaxies in large multi-dimensional astronomical databases like the Sloan Digital Sky Survey (SDSS). Specifically, this paper presents four new advances: (1) A new semi-parametric algorithm - nicknamed ``C4'' - for jointly finding clusters of galaxies in the SDSS and ROSAT All-Sky Survey databases; (2) The introduction of the False Discovery Rate into Astronomy; (3) The role of kernel shape in optimizing cluster detection; (4) A new determination of the X-ray Cluster Luminosity Function which has bearing on the existence of a ``deficit'' of high redshift, high luminosity clusters. This research is part of our ``Computational AstroStatistics'' collaboration (see Nichol et al. 2000) and the algorithms and techniques discussed herein will form part of the ``Virtual Observatory'' analysis toolkit.