We propose a deep learning analyzing technique with convolutional neural network (CNN) to predict the evolutionary track of the Epoch of Reionization (EoR) from the 21-cm differential brightness temperature tomography images. We use 21cmFAST, a fast semi-numerical cosmological 21-cm signal simulator, to produce mock 21-cm maps between $z=6 \sim 13$. We then apply two observational effects into those 21-cm maps, such as instrumental noise and limit of (spatial and depth) resolution somewhat suitable for realistic choices of the Square Kilometre Array (SKA). We design our deep learning model with CNN to predict the sliced-averaged neutral hydrogen fraction from the given 21-cm map. The estimated neutral fraction from our CNN model has a great agreement with its true value even after coarsely smoothing with broad beamsize and frequency bandwidth, and also heavily covered by noise with narrow. Our results have shown that deep learning analyzing method has a large potential to efficiently reconstruct the EoR history from the 21-cm tomography surveys in future.

We analyze an extremely deep 450-$\mu$m image ($1\sigma=0.56$\,mJy\,beam$^{-1}$) of a $\simeq 300$\,arcmin$^{2}$ area in the CANDELS/COSMOS field as part of the SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). We select a robust (signal-to-noise ratio $\geqslant 4$) and flux-limited ($\geqslant 4$\,mJy) sample of 164 sub-millimeter galaxies (SMGs) at 450-$\mu$m that have $K$-band counterparts in the COSMOS2015 catalog identified from radio or mid-infrared imaging. Utilizing this SMG sample and the 4705 $K$-band-selected non-SMGs that reside within the noise level $\leqslant 1$\,mJy\,beam$^{-1}$ region of the 450-$\mu$m image as a training set, we develop a machine-learning classifier using $K$-band magnitude and color-color pairs based on the thirteen-band photometry available in this field. We apply the trained machine-learning classifier to the wider COSMOS field (1.6\,deg$^{2}$) using the same COSMOS2015 catalog and identify a sample of 6182 450-$\mu$m SMG candidates with similar colors. The number density, radio and/or mid-infrared detection rates, redshift and stellar mass distributions, and the stacked 450-$\mu$m fluxes of these SMG candidates, from the S2COSMOS observations of the wide field, agree with the measurements made in the much smaller CANDELS field, supporting the effectiveness of the classifier. Using this 450-$\mu$m SMG candidate sample, we measure the two-point autocorrelation functions from $z=3$ down to $z=0.5$. We find that the 450-$\mu$m SMG candidates reside in halos with masses of $\simeq (2.0\pm0.5) \times10^{13}\,h^{-1}\,\rm M_{\odot}$ across this redshift range. We do not find evidence of downsizing that has been suggested by other recent observational studies.

Horizon Run 5 (HR5) is a cosmological hydrodynamics simulation which captures the properties of the Universe on a Gpc scale while achieving a resolution of 1 kpc. This enormous dynamic range allows us to simultaneously capture the physics of the cosmic web on very large scales and account for the formation and evolution of dwarf galaxies on much smaller scales. Inside the simulation box we zoom-in on a high-resolution cuboid region with a volume of 1049x114x114 Mpc^3. The sub-grid physics chosen to model galaxy formation includes radiative heating/cooling, reionization, star formation, supernova feedback, chemical evolution tracking the enrichment of oxygen and iron, the growth of supermassive black holes and active galactic nuclei (AGN) feedback in the form of a dual jet-heating mode. For this simulation we implemented a hybrid MPI-OpenMP version of the RAMSES code, specifically targeted for modern many-core many thread parallel architectures. For the post-processing, we extended the Friends-of-Friend (FoF) algorithm and developed a new galaxy finder to analyse the large outputs of HR5. The simulation successfully reproduces many observations, such as the cosmic star formation history, connectivity of galaxy distribution and stellar mass functions. The simulation also indicates that hydrodynamical effects on small scales impact galaxy clustering up to very large scales near and beyond the baryonic acoustic oscillation (BAO) scale. Hence, caution should be taken when using that scale as a cosmic standard ruler: one should carefully understand the corresponding biases. The simulation is expected to be an invaluable asset for the interpretation of upcoming deep surveys of the Universe.

The redshift-space distortion (RSD) in the observed distribution of galaxies is known as a powerful probe of cosmology. Observations of large-scale RSD have given tight constraints on the linear growth rate of the large-scale structures in the universe. On the other hand, the small-scale RSD, caused by galaxy random motions inside clusters, has not been much used in cosmology, but also has cosmological information because universes with different cosmological parameters have different halo mass functions and virialized velocities. We focus on the projected correlation function $w(r_p)$ and the multipole moments $\xi_l$ on small scales ($1.4$ to $30\ h^{-1}\rm{Mpc}$). Using simulated galaxy samples generated from a physically motivated most bound particle (MBP)-galaxy correspondence scheme in the Multiverse Simulation, we examine the dependence of the small-scale RSD on the cosmological matter density parameter $\Omega_m$, the satellite velocity bias with respect to MBPs, $b_v^s$, and the merger-time-scale parameter $\alpha$. We find that $\alpha=1.5$ gives an excellent fit to the $w(r_p)$ and $\xi_l$ measured from the SDSS-KIAS value added galaxy catalog. We also define the ``strength'' of Fingers-of-God as the ratio of the parallel and perpendicular size of the contour in the two-point correlation function set by a specific threshold value and show that the strength parameter helps constraining $(\Omega_m, b_v^s, \alpha)$ by breaking the degeneracy among them. The resulting parameter values from all measurements are $(\Omega_m,b_v^s)=(0.272\pm0.013,0.982\pm0.040)$, indicating a slight reduction of satellite galaxy velocity relative to the MBP. However, considering that the average MBP speed inside haloes is $0.94$ times the dark matter velocity dispersion, the main drivers behind the galaxy velocity bias are gravitational interactions, rather than baryonic effects.

We present measurements of the two-dimensional genus of the SDSS-III BOSS catalogs to constrain cosmological parameters governing the shape of the matter power spectrum. The BOSS data are divided into twelve concentric shells over the redshift range $0.2 < z < 0.6$, and we extract the genus from the projected two-dimensional galaxy density fields. We compare the genus amplitudes to their Gaussian expectation values, exploiting the fact that this quantity is relatively insensitive to non-linear gravitational collapse. The genus amplitude provides a measure of the shape of the linear matter power spectrum, and is principally sensitive to $\Omega_{\rm c}h^{2}$ and scalar spectral index $n_{\rm s}$. A strong negative degeneracy between $\Omega_{\rm c}h^{2}$ and $n_{\rm s}$ is observed, as both can increase small scale power by shifting the peak and tilting the power spectrum respectively. We place a constraint on the particular combination $n_{\rm s}^{3/2} \Omega_{\rm c}h^{2}$ -- we find $n_{\rm s}^{3/2} \Omega_{\rm c}h^{2} = 0.1121 \pm 0.0043$ after combining the LOWZ and CMASS data sets, assuming a flat $\Lambda$CDM cosmology. This result is practically insensitive to reasonable variations of the power spectrum amplitude and linear galaxy bias. Our results are consistent with the Planck best fit $n_{\rm s}^{3/2}\Omega_{\rm c}h^{2} = 0.1139 \pm 0.0009$.

We develop an improved Alcock-Paczynski (AP) test method that uses the redshift-space two-point correlation function (2pCF) of galaxies. Cosmological constraints can be obtained by examining the redshift dependence of the normalized 2pCF, which should not change apart from the expected small non-linear evolution. An incorrect choice of cosmology used to convert redshift to comoving distance will manifest itself as redshift-dependent 2pCF. Our method decomposes the redshift difference of the two-dimensional correlation function into the Legendre polynomials whose amplitudes are modeled by radial fitting functions. Our likelihood analysis with this 2-D fitting scheme tightens the constraints on $\Omega_m$ and ${w}$ by $\sim 40\%$ compared to the method of Li et al. (2016, 2017, 2018) that uses one dimensional angular dependence only. We also find that the correction for the non-linear evolution in the 2pCF has a non-negligible cosmology dependence, which has been neglected in previous similar studies by Li et al.. With an accurate accounting for the non-linear systematics and use of full two-dimensional shape information of the 2pCF down to scales as small as $5~h^{-1}{\rm Mpc}$ it is expected that the AP test with redshift-space galaxy clustering anisotropy can be a powerful method to constrain the expansion history of the universe.

By utilizing large-scale graph analytic tools implemented in the modern Big Data platform, Apache Spark, we investigate the topological structure of gravitational clustering in five different universes produced by cosmological $N$-body simulations with varying parameters: (1) a WMAP 5-year compatible $\Lambda$CDM cosmology, (2) two different dark energy equation of state variants, and (3) two different cosmic matter density variants. For the Big Data calculations, we use a custom build of stand-alone Spark/Hadoop cluster at Korea Institute for Advanced Study (KIAS) and Dataproc Compute Engine in Google Cloud Platform (GCP) with the sample size ranging from 7 millions to 200 millions. We find that among the many possible graph-topological measures, three simple ones: (1) the average of number of neighbors (the so-called average vertex degree) $\alpha$, (2) closed-to-connected triple fraction (the so-called transitivity) $\tau_\Delta$, and (3) the cumulative number density $n_{s\ge5}$ of subcomponents with connected component size $s \ge 5$, can effectively discriminate among the five model universes. Since these graph-topological measures are in direct relation with the usual $n$-points correlation functions of the cosmic density field, graph-topological statistics powered by Big Data computational infrastructure opens a new, intuitive, and computationally efficient window into the dark Universe.

By utilizing large-scale graph analytic tools implemented in the modern Big Data platform, Apache Spark, we investigate the topological structure of gravitational clustering in five different universes produced by cosmological $N$-body simulations with varying parameters: (1) a WMAP 5-year compatible $\Lambda$CDM cosmology, (2) two different dark energy equation of state variants, and (3) two different cosmic matter density variants. For the Big Data calculations, we use a custom build of stand-alone Spark/Hadoop cluster at Korea Institute for Advanced Study (KIAS) and Dataproc Compute Engine in Google Cloud Platform (GCP) with the sample size ranging from 7 millions to 200 millions. We find that among the many possible graph-topological measures, three simple ones: (1) the average of number of neighbors (the so-called average vertex degree) $\alpha$, (2) closed-to-connected triple fraction (the so-called transitivity) $\tau_\Delta$, and (3) the cumulative number density $n_{s\ge5}$ of subcomponents with connected component size $s \ge 5$, can effectively discriminate among the five model universes. Since these graph-topological measures are in direct relation with the usual $n$-points correlation functions of the cosmic density field, graph-topological statistics powered by Big Data computational infrastructure opens a new, intuitive, and computationally efficient window into the dark Universe.

We present a detailed census of galaxies in and around PC217.96+32.3, a spectroscopically confirmed Coma analog at z=3.78. Diverse galaxy types identified in the field include Lya emitters (LAEs), massive star-forming galaxies, and ultra-massive galaxies (log (Mstar/Msun)>= 11) which may have already halted their star formation. The sky distribution of the star-forming galaxies suggests the presence of a significant overdensity (delta_g=8+/-2), which is spatially offset from the previously confirmed members by 3-4 Mpc to the west. Candidate quiescent and post-starburst galaxies are also found in large excess (a factor of ~ 8-15 higher surface density than the field) although their redshifts are less certain. We estimate that the total enclosed mass traced by star-forming galaxy candidates is roughly comparable to that of PC217.96+32.3 traced by the LAEs. We speculate that the true extent of PC217.96+32.3 may be larger than previously known, a half of which is missed by our LAE selection. Alternatively, the newly discovered overdensity may belong to another Coma progenitor not associated with PC217.96+32.3. Expectations from theory suggest that both scenarios are equally unlikely (<1%), particularly in the cosmic volume probed in our survey. If confirmed as a single structure, its total mass will be well in excess of Coma, making this a singularly large cosmic structure rarely seen even in large cosmological simulations. Finally, we find that the protocluster galaxies follow the same SFR-M_star scaling relation as the field galaxies, suggesting that the environmental effect at z~4 is a subtle one at best for normal star-forming galaxies.

We investigate the spatial distribution of Lyman alpha emitting galaxies (LAEs) at $z \approx 2.67$, selected from the NOAO Deep Wide-Field Survey (NDWFS), using two-point statistics and topological diagnostics adopted from network science. We measure the clustering length, $r_0 \approx 4 h^{-1}$ Mpc, and the bias, $b_{LAE} = 2.2^{+0.2}_{-0.1}$. Fitting the clustering with halo occupation distribution (HOD) models results in two disparate possibilities: (1) where the fraction of central galaxies is $<$1% in halos of mass $>10^{12}$$M_\odot$; and (2) where the fraction is $\approx$20%. We refer to these two scenarios as the `Dusty Core Scenario' for Model#1 since most of central galaxies in massive halos are dead in Ly$\alpha$ emission, and the `Pristine Core Scenario' for Model#2 since the central galaxies are bright in Ly$\alpha$ emission. Traditional two-point statistics cannot distinguish between these disparate models given the current data sets. To overcome this degeneracy, we generate mock catalogs for each HOD model using a high resolution $N$-body simulation and adopt a network statistics approach, which provides excellent topological diagnostics for galaxy point distributions. We find three topological anomalies from the spatial distribution of observed LAEs, which are not reproduced by the HOD mocks. We find that Model#2 matches better all network statistics than Model#1, suggesting that the central galaxies in $> 10^{12} h^{-1} M_\odot$ halos at $z \approx 2.67$ need to be less dusty to be bright as LAEs, potentially implying some replenishing channels of pristine gas such as the cold mode accretion.

Devasthal Optical Telescope Integral Field Spectrograph (DOTIFS) is a new multi-Integral Field Unit (IFU) instrument, planned to be mounted on the 3.6m Devasthal optical telescope in Nainital, India. It has eight identical, fiber-fed spectrographs to disperse light coming from 16 IFUs. The spectrographs produce 2,304 spectra over a 370-740nm wavelength range simultaneously with a spectral resolution of R=1200-2400. It is composed of all-refractive, all spherical optics designed to achieve on average 26.0% throughput from the telescope to the CCD with the help of high transmission spectrograph optics, volume phase holographic grating, and graded coated e2v 2K by 4K CCD. We present the optical and opto-mechanical design of the spectrograph as well as current development status. Optics and optomechanical components for the spectrographs are being fabricated.

IUCAA Digital Sampling Array Controller (IDSAC) is a flexible and generic yet powerful CCD controller which can handle a wide range of scientific detectors. Based on an easily scalable modular backplane architecture consisting of Single Board Controllers (SBC), IDSAC can control large detector arrays and mosaics. Each of the SBCs offers the full functionality required to control a CCD independently. The SBCs can be cold swapped without the need to reconfigure them. Each SBC can handle data from up to four video channels with or without dummy outputs at speeds up to 500 kilo Pixels Per Second (kPPS) Per Channel with a resolution of 16 bits. Communication with Linux based host computer is through a USB3.0 interface, with the option of using copper or optical fibers. A Field Programmable Gate Array (FPGA) is used as the master controller in each SBC which allows great flexibility in optimizing performance by adjusting gain, timing signals, bias levels, etc. using user-editable configuration files without altering the circuit topology. Elimination of thermal kTC noise is achieved via Digital Correlated Double Sampling (DCDS). We present the results of noise performance characterization of IDSAC through simulation, theoretical modeling, and actual measurements. The contribution of different types of noise sources is modeled using a tool to predict noise of a generic DCDS signal chain analytically. The analytical model predicts the net input referenced noise of the signal chain to be 5 electrons for 200k pixels per second per channel readout rate with 3 samples per pixel. Using a cryogenic test set up in the lab, the noise is measured to be 5.4 e (24.3 \muV), for the same readout configuration.

A new field re-configuration technique, Multiple Rooks of Chess (MRC), for multiple deployable Integral Field Spectrographs has been developed. The method involves mechanical geometry as well as an optimized deployment algorithm. The geometry is found to be simple for mechanical implementation. The algorithm initially assigns the IFUs to the target objects and then devises the movement sequence based on the current and the desired IFU positions. The reconfiguration time using the suitable actuators which runs at 20 cm/s is found to be a maximum of 25 seconds for the circular DOTIFS focal plane (180 mm diameter). The Geometry Algorithm Combination (GAC) has been tested on several million mock target configurations with object-to-IFU ({\tau} ) ratio varying from 0.25 to 16. The MRC method is found to-be efficient in target acquisition in terms of field revisit and deployment time without any collision or entanglement of the fiber bundles. The efficiency of the technique does not get affected by the increase in number density of target objects. The technique is compared with other available methods based on sky coverage, flexibility and overhead time. The proposed geometry and algorithm combination is found to have an advantage in all of the aspects.

We report the discovery and spectroscopic confirmation of a very large star-forming Lyman Break galaxy, G6025, at z_spec=3.721+/-0.003. In the rest-frame ~2100A, G6025 subtends ~24 kpc in physical extent when measured from the 1.5-sigma isophote, in agreement with the parametric size measurements which yield the half-light radius of 4.9+/-0.5 kpc and the semi-major axis of 12.5+/-0.1 kpc. G6025 is also very UV-luminous (~5L*(z~4}) and young (~140+/-60 Myr). Despite its unusual size and luminosity, the stellar population parameters and dust reddening (M_star~M*(z~4)$, and E(B-V)=0.18+/-0.05) estimated from the integrated light, are similar to those of smaller galaxies at comparable redshifts. The ground-based morphology and spectroscopy show two dominant components, both located off-center, embedded in more diffuse emission. We speculate that G6025 may be a scaled-up version of chain galaxies seen in deep HST imaging, or alternatively, a nearly equal-mass merger involving two super-L* galaxies in its early stage. G6025 lies close to but not within a known massive protocluster at z=3.78. We find four companions within 6 Mpc from G6025, two of which lie within 1.6 Mpc. While the limited sensitivity of the existing spectroscopy does not allow us to robustly characterize the local environment of G6025, it likely resides in a locally overdense environment. The luminosity, size, and youth of G6025 make it uniquely suited to study the early formation of massive galaxies in the universe.

The two red clumps (RCs) observed in the color-magnitude diagram of the Milky Way bulge is widely accepted as evidence for an X-shaped structure originated from the bar instability. A drastically different interpretation has been suggested, however, based on the He-enhanced multiple stellar population phenomenon as is observed in globular clusters (GCs). Because these two scenarios imply very different pictures on the formation of the bulge and elliptical galaxies, understanding the origin of the double RC is of crucial importance. Here we report our discovery that the stars in the two RCs show a significant (> 5.3 {\sigma}) difference in CN-band strength, in stark contrast to that expected in the X-shaped bulge scenario. The difference in CN abundance and the population ratio between the two RCs are comparable to those observed in GCs between the first- and later generation stars. Since CN-strong stars trace a population with enhanced N, Na, and He abundances originated in GCs, this is direct evidence that the double RC is due to the multiple population phenomenon, and that a significant population of stars in the Milky Way bulge were assembled from disrupted proto-GCs. Our result also calls for the major revision of the 3D structure of the Milky Way bulge given that the current view is based on the previous interpretation of the double RC phenomenon.

We apply the Minkowski Tensor statistics to two dimensional slices of the three dimensional density field. The Minkowski Tensors are a set of functions that are sensitive to directionally dependent signals in the data, and furthermore can be used to quantify the mean shape of density peaks. We begin by introducing our algorithm for constructing bounding perimeters around subsets of a two dimensional field, and reviewing the definition of Minkowski Tensors. Focusing on the translational invariant statistic $W^{1,1}_{2}$ - a $2 \times 2$ matrix - we calculate its eigenvalues for both the entire excursion set ($\Lambda_{1},\Lambda_{2}$) and for individual connected regions and holes within the set ($\lambda_{1},\lambda_{2}$). The ratio of eigenvalues $\Lambda_{2}/\Lambda_{1}$ informs us of the presence of global anisotropies in the data, and $\langle \lambda_{2}/\lambda_{1} \rangle$ is a measure of the mean shape of peaks and troughs in the density field. We study these quantities for a Gaussian field, then consider how they are modified by the effect of gravitational collapse using the latest Horizon Run 4 cosmological simulation. We find $\Lambda_{1,2}$ are essentially independent of gravitational collapse, as the process maintains statistical isotropy. However, the mean shape of peaks is modified significantly - overdensities become relatively more circular compared to underdensities of the same area. When applying the statistic to a redshift space distorted density field, we find a significant signal in the eigenvalues $\Lambda_{1,2}$, suggesting that they can be used to probe the large-scale velocity field.

We perform an anisotropic clustering analysis of 1,133,326 galaxies from the Sloan Digital Sky Survey (SDSS-III) Baryon Oscillation Spectroscopic Survey (BOSS) Data Release (DR) 12 covering the redshift range $0.15<z<0.69$. The geometrical distortions of the galaxy positions, caused by incorrect cosmological model assumptions, are captured in the anisotropic two-point correlation function on scales 6 -- 40 $h^{-1}\rm Mpc$. The redshift evolution of this anisotropic clustering is used to place constraints on the cosmological parameters. We improve the methodology of Li et al. 2016, to enable efficient exploration of high dimensional cosmological parameter spaces, and apply it to the Chevallier-Polarski-Linder parametrization of dark energy, $w=w_0+w_a{z}/({1+z})$. In combination with the CMB, BAO, SNIa and $H_0$ from Cepheid data, we obtain $\Omega_m = 0.301 \pm 0.008,\ w_0 = -1.042 \pm 0.067,\ $ and $w_a = -0.07 \pm 0.29$ (68.3\% CL). Adding our new AP measurements to the aforementioned results reduces the error bars by $\sim$30 -- 40\% and improves the dark energy figure of merit by a factor of $\sim$2. We check the robustness of the results using realistic mock galaxy catalogues.

We study the topology of the matter density field in two dimensional slices, and consider how we can use the amplitude $A$ of the genus for cosmological parameter estimation. Using the latest Horizon Run 4 simulation data, we calculate the genus of the smoothed density field constructed from lightcone mock galaxy catalogs. Information can be extracted from the amplitude of the genus by considering both its redshift evolution and magnitude. The constancy of the genus amplitude with redshift can be used as a standard population, from which we derive constraints on the equation of state of dark energy $w_{\rm de}$ - by measuring $A$ at $z \sim 0.1$ and $z \sim 1$, we can place an order $\Delta w_{\rm de} \sim {\cal O}(15\%)$ constraint on $w_{\rm de}$. By comparing $A$ to its Gaussian expectation value we can potentially derive an additional stringent constraint on the matter density $\Delta \Omega_{\rm mat} \sim 0.01$. We discuss the primary sources of contamination associated with the two measurements - redshift space distortion and shot noise. With accurate knowledge of galaxy bias, we can successfully remove the effect of redshift space distortion, and the combined effect of shot noise and non-linear gravitational evolution is suppressed by smoothing over suitably large scales $R_{\rm G} \ge 15 {\rm Mpc}/h$. Without knowledge of the bias, we discuss how joint measurements of the two and three dimensional genus can be used to constrain the growth factor $\beta = f/b$. The method can be applied optimally to redshift slices of a galaxy distribution generated using the drop-off technique.

In order to search for the evidence of quasi-periodic oscillation (QPO) in blazar, multicolor optical observations of the BL Lacertae object S5 0716+714 were performed from 2005 to 2012. For $I$ band on March 19 2010 with low variability amplitude and low flux level, the same quasi-periodic oscillation $\simeq50$ minutes with 99% significance levels is confirmed by ZDCF method, Lomb-Scargle method, REDFIT and fitting sinusoidal curves. The observed QPO is likely to be explained by accretion disk variability. If the observed QPO indicates an innermost stable orbital period from the accretion disk, the QPO $\simeq50$ min corresponds to a black hole mass of $5.03\times10^6 M_\odot$ for a non-rotating Schwarzschild black hole and $3.2\times10^7 M_\odot$ for a maximally rotating Kerr black hole.

We study the impact of thermal inflation on the formation of cosmological structures and present astrophysical observables which can be used to constrain and possibly probe the thermal inflation scenario. These are dark matter halo abundance at high redshifts, satellite galaxy abundance in the Milky Way, and fluctuation in the 21-cm radiation background before the epoch of reionization. The thermal inflation scenario leaves a characteristic signature on the matter power spectrum by boosting the amplitude at a specific wavenumber determined by the number of e-foldings during thermal inflation ($N_{\rm bc}$), and strongly suppressing the amplitude for modes at smaller scales. For a reasonable range of parameter space, one of the consequences is the suppression of minihalo formation at high redshifts and that of satellite galaxies in the Milky Way. While this effect is substantial, it is degenerate with other cosmological or astrophysical effects. The power spectrum of the 21-cm background probes this impact more directly, and its observation may be the best way to constrain the thermal inflation scenario due to the characteristic signature in the power spectrum. The Square Kilometre Array (SKA) in phase 1 (SKA1) has sensitivity large enough to achieve this goal for models with $N_{\rm bc}\gtrsim 26$ if a 10000-hr observation is performed. The final phase SKA, with anticipated sensitivity about an order of magnitude higher, seems more promising and will cover a wider parameter space.

The standard paradigm of collisionless cold dark matter is in tension with measurements on large scales. In particular, the best fit values of the Hubble rate $H_0$ and the matter density perturbation $\sigma_8$ inferred from the cosmic microwave background seem inconsistent with the results from direct measurements. We show that both problems can be solved in a framework in which dark matter consists of two distinct components, a dominant component and a subdominant component. The primary component is cold and collisionless. The secondary component is also cold, but interacts strongly with dark radiation, which itself forms a tightly coupled fluid. The growth of density perturbations in the subdominant component is inhibited by dark acoustic oscillations due to its coupling to the dark radiation, solving the $\sigma_8$ problem, while the presence of tightly coupled dark radiation ameliorates the $H_0$ problem. The subdominant component of dark matter and dark radiation continue to remain in thermal equilibrium until late times, inhibiting the formation of a dark disk. We present an example of a simple model that naturally realizes this scenario in which both constituents of dark matter are thermal WIMPs. Our scenario can be tested by future stage-IV experiments designed to probe the CMB and large scale structure.

Thermal inflation is a second epoch of exponential expansion at typical energy scales $V^{1/4} \sim 10^{6 \sim 8} \mathrm{GeV}$. If the usual primordial inflation is followed by thermal inflation, the primordial power spectrum is only modestly redshifted on large scales, but strongly suppressed on scales smaller than the horizon size at the beginning of thermal inflation, $k > k_{\rm b} = a_{\rm b} H_{\rm b}$. We calculate the spectral distortion of the cosmic microwave background generated by the dissipation of acoustic waves in this context. For $k_{\rm b} \ll 10^3 \mathrm{Mpc}^{-1}$, thermal inflation results in a large suppression of the $\mu$-distortion amplitude, predicting that it falls well below the standard value of $\mu \simeq 2\times 10^{-8}$. Thus, future spectral distortion experiments, similar to PIXIE, can place new limits on the thermal inflation scenario, constraining $k_{\rm b} \gtrsim 10^3 \mathrm{Mpc}^{-1}$ if $\mu \simeq 2\times 10^{-8}$ were found.

We develop a methodology to use the redshift dependence of the galaxy 2-point correlation function (2pCF) across the line-of-sight, $\xi(r_{\bot})$, as a probe of cosmological parameters. The positions of galaxies in comoving Cartesian space varies under different cosmological parameter choices, inducing a {\it redshift-dependent scaling} in the galaxy distribution. This geometrical distortion can be observed as a redshift-dependent rescaling in the measured $\xi(r_{\bot})$. We test this methodology using a sample of 1.75 billion mock galaxies at redshifts 0, 0.5, 1, 1.5, 2, drawn from the Horizon Run 4 N-body simulation. The shape of $\xi(r_{\bot})$ can exhibit a significant redshift evolution when the galaxy sample is analyzed under a cosmology differing from the true, simulated one. Other contributions, including the gravitational growth of structure, galaxy bias, and the redshift space distortions, do not produce large redshift evolution in the shape. We show that one can make use of this geometrical distortion to constrain the values of cosmological parameters governing the expansion history of the universe. This method could be applicable to future large scale structure surveys, especially photometric surveys such as DES, LSST, to derive tight cosmological constraints. This work is a continuation of our previous works as a strategy to constrain cosmological parameters using redshift-invariant physical quantities.

We performed low-resolution spectroscopy for the red giant stars in the Galactic globular cluster (GC) NGC 5286, which is known to show intrinsic heavy element abundance variations. We found that the observed stars in this GC are clearly divided into three subpopulations by CN index (CN-weak, CN-intermediate, and CN-strong). The CN-strong stars are also enhanced in the calcium HK' (7.4{\sigma}) and CH (5.1{\sigma}) indices, while the CN-intermediate stars show no significant difference in the strength of HK' index with CN-weak stars. From the comparison with high-resolution spectroscopic data, we found that the CN- and HK'-strong stars are also enhanced in the abundances of Fe and s-process elements. It appears, therefore, that these stars are later generation stars affected by some supernovae enrichment in addition to the asymptotic giant branch ejecta. In addition, unlike normal GCs, sample stars in NGC 5286 show the CN-CH positive correlation, strengthening our previous suggestion that this positive correlation is only discovered in GCs with heavy element abundance variations such as M22 and NGC 6273.

We have monitored the BL Lacertae object S5 0716+714 in the optical bands during 2012 January and February with long time spans on intraday timescales ($>$5 hr) and high time resolutions. During this monitoring period, the object shows violent flaring activity both in short and intraday timescales. The object has high value of duty cycle. The light curves detected as intraday variability (IDV) show variability of various shapes. The variability amplitude is from 12.81 per cent to 33.22 per cent, and the average value $19.92\pm5.87$ per cent. The overall magnitude variabilities are $\bigtriangleup B=1^{\rm m}.24$, $\bigtriangleup V=1^{\rm m}.42$, $\bigtriangleup R=1^{\rm m}.3$, $\bigtriangleup I=1^{\rm m}.23$ respectively. During the observations, the average change rate is $<CR>=0.035\pm0.009$ Mag/h during the ascent and $<CR>=0.035\pm0.014$ Mag/h during the descent. However, different cases are found on certain nights. There are good inter-bands correlations but not significant time lags for intraday and short timescales. The results of the autocorrelation function show that the variability timescales range from 0.054 day to 0.134 day. Most of nights show bluer when brighter (BWB) chromatic trend; a weak redder with brighter (RWB) trend is found; a few nights show no correlations between magnitude and color index. The BWB trend appears in the short timescales. During the flare the spectral index exhibits a clockwise loop for inter-nights. A shock-in-jet model and the shock wave propagating along a helical path are likely to explain the variability and color index variability.

Starting from a very accurate model for density-in-cells statistics of dark matter based on large deviation theory, a bias model for the tracer density in spheres is formulated. It adopts a mean bias relation based on a quadratic bias model to relate the log-densities of dark matter to those of mass-weighted dark haloes in real and redshift space. The validity of the parametrised bias model is established using a parametrisation-independent extraction of the bias function. This average bias model is then combined with the dark matter PDF, neglecting any scatter around it: it nevertheless yields an excellent model for densities-in-cells statistics of mass tracers that is parametrised in terms of the underlying dark matter variance and three bias parameters. The procedure is validated on measurements of both the one and two point statistics of subhalo densities in the state-of-the-art Horizon Run 4 simulation showing excellent agreement for measured dark matter variance and bias parameters. Finally, it is demonstrated that this formalism allows for a joint estimation of the nonlinear dark matter variance and the bias parameters using solely the statistics of subhaloes. Having verified that galaxy counts in hydrodynamical simulations sampled on a scale of 10 Mpc/h closely resemble those of subhaloes, this work provides important steps towards making theoretical predictions for density-in-cells statistics applicable to upcoming galaxy surveys like Euclid or WFIRST.

We report the detection of diffuse Ly$\alpha$ emission, or Ly$\alpha$ halos (LAHs), around star-forming galaxies at $z\approx3.78$ and $2.66$ in the NOAO Deep Wide-Field Survey Bo\"otes field. Our samples consist of a total of $\sim$1400 galaxies, within two separate regions containing spectroscopically confirmed galaxy overdensities. They provide a unique opportunity to investigate how the LAH characteristics vary with host galaxy large-scale environment and physical properties. We stack Ly$\alpha$ images of different samples defined by these properties and measure their median LAH sizes by decomposing the stacked Ly$\alpha$ radial profile into a compact galaxy-like and an extended halo-like component. We find that the exponential scale-length of LAHs depends on UV continuum and Ly$\alpha$ luminosities, but not on Ly$\alpha$ equivalent widths or galaxy overdensity parameters. The full samples, which are dominated by low UV-continuum luminosity Ly$\alpha$ emitters ($M_{\rm UV} \gtrsim -21$), exhibit LAH sizes of 5$\,-\,6\,$kpc. However, the most UV- or Ly$\alpha$-luminous galaxies have more extended halos with scale-lengths of 7$\,-\,9\,$kpc. The stacked Ly$\alpha$ radial profiles decline more steeply than recent theoretical predictions that include the contributions from gravitational cooling of infalling gas and from low-level star formation in satellites. On the other hand, the LAH extent matches what one would expect for photons produced in the galaxy and then resonantly scattered by gas in an outflowing envelope. The observed trends of LAH sizes with host galaxy properties suggest that the physical conditions of the circumgalactic medium (covering fraction, HI column density, and outflow velocity) change with halo mass and/or star-formation rates.

We study the two-dimensional topology of the galactic distribution when projected onto two-dimensional spherical shells. Using the latest Horizon Run 4 simulation data, we construct the genus of the two-dimensional field and consider how this statistic is affected by late-time nonlinear effects -- principally gravitational collapse and redshift space distortion (RSD). We also consider systematic and numerical artifacts such as shot noise, galaxy bias, and finite pixel effects. We model the systematics using a Hermite polynomial expansion and perform a comprehensive analysis of known effects on the two-dimensional genus, with a view toward using the statistic for cosmological parameter estimation. We find that the finite pixel effect is dominated by an amplitude drop and can be made less than $1\%$ by adopting pixels smaller than $1/3$ of the angular smoothing length. Nonlinear gravitational evolution introduces time-dependent coefficients of the zeroth, first, and second Hermite polynomials, but the genus amplitude changes by less than $1\%$ between $z=1$ and $z=0$ for smoothing scales $R_{\rm G} > 9 {\rm Mpc/h}$. Non-zero terms are measured up to third order in the Hermite polynomial expansion when studying RSD. Differences in shapes of the genus curves in real and redshift space are small when we adopt thick redshift shells, but the amplitude change remains a significant $\sim {\cal O}(10\%)$ effect. The combined effects of galaxy biasing and shot noise produce systematic effects up to the second Hermite polynomial. It is shown that, when sampling, the use of galaxy mass cuts significantly reduces the effect of shot noise relative to random sampling.

In cosmological scenarios with thermal inflation, extra eras of moduli matter domination, thermal inflation and flaton matter domination exist between primordial inflation and the radiation domination of Big Bang nucleosynthesis. During these eras, cosmological perturbations on small scales can enter and re-exit the horizon, modifying the power spectrum on those scales. The largest modified scale, $k_\mathrm{b}$, touches the horizon size when the expansion changes from deflation to inflation at the transition from moduli domination to thermal inflation. We analytically calculate the evolution of perturbations from moduli domination through thermal inflation and evaluate the curvature perturbation on the constant radiation density hypersurface at the end of thermal inflation to determine the late time curvature perturbation. Our resulting transfer function suppresses the power spectrum by a factor $\sim 50$ at $k \gg k_\mathrm{b}$, with $k_\mathrm{b}$ corresponding to anywhere from megaparsec to subparsec scales depending on the parameters of thermal inflation. Thus, thermal inflation might be constrained or detected by small scale observations such as CMB distortions or 21cm hydrogen line observations.

It has been known that some solar activity indicators show a double-peak feature in their evolution through a solar cycle, which is not conspicuous in sunspot number. In this letter, we investigate the high solar dynamic activity in the declining phase of the sunspot cycle by examining the evolution of polar and low latitude coronal hole areas and the statistics of splitting and merging events of coronal holes and coronal mass ejections detected by SOHO/LASCO C3 in solar cycle 23. Although the total coronal hole area is at its maximum near the sunspot minimum, in which polar coronal holes prevail, it shows a comparable second maximum in the declining phase of the cycle, in which low latitude coronal holes are dominant. The events of coronal hole splitting or merging, which are attributed to surface motions of magnetic fluxes, are also mostly populated in the declining phase of the cycle. The far-reaching C3 coronal mass ejections are also over-populated in the declining phase of the cycle. From these results we suggest that solar dynamic activities due to the horizontal motions of magnetic fluxes extend far in the declining phase of the sunspot cycle.

We apply the methodology developed in \cite{Li2014,Li2015} to BOSS DR12 galaxies and derive cosmological constraints from the redshift dependence of the Alcock-Paczynski (AP) effect. The apparent anisotropy in the distribution of observed galaxies arise from two main sources, the redshift-space distortion (RSD) effect due to the galaxy peculiar velocities, and the geometric distortion when incorrect cosmological models are assumed for transforming redshift to comoving distance, known as the AP effect. Anisotropies produced by the RSD effect are, although large, maintaining a nearly uniform magnitude over a large range of redshift, while the degree of anisotropies from the AP effect varies with redshift by much larger magnitude. We split the DR12 galaxies into six redshift bins, measure the 2-point correlation function in each bin, and assess the redshift evolution of anisotropies. We obtain constraints of $\Omega_m=0.290 \pm 0.053,\ \ w = -1.07 \pm 0.15$, which are comparable with the current constraints from other cosmological probes such as type Ia supernovae, cosmic microwave background, and baryon acoustic oscillation (BAO). Combining these cosmological probes with our method yield tight constraints of $ \Omega_m = 0.301 \pm 0.006,\ w=-1.054 \pm 0.025$. Our method is complementary to the other large scale structure probes like BAO and topology. We expect this technique will play an important role in deriving cosmological constraints from large scale structure surveys.

The Australian Astronomical Observatory's TAIPAN instrument deploys 150 Starbug robots to position optical fibres to accuracies of 0.3 arcsec, on a 32 cm glass field plate on the focal plane of the 1.2 m UK-Schmidt telescope. This paper describes the software system developed to control and monitor the Starbugs, with particular emphasis on the automated path-finding algorithms, and the metrology software which keeps track of the position and motion of individual Starbugs as they independently move in a crowded field. The software employs a tiered approach to find a collision-free path for every Starbug, from its current position to its target location. This consists of three path-finding stages of increasing complexity and computational cost. For each Starbug a path is attempted using a simple method. If unsuccessful, subsequently more complex (and expensive) methods are tried until a valid path is found or the target is flagged as unreachable.

We combine a semi-analytic model of galaxy evolution with constraints on circumstellar habitable zones and the distribution of terrestrial planets to probe the suitability of galaxies of different mass and type to host habitable planets, and how it evolves with time. We find that the fraction of stars with terrestrial planets in their habitable zone (known as habitability) depends only weakly on galaxy mass, with a maximum around 4e10 Msun. We estimate that 0.7% of all stars in Milky Way type galaxies to host a terrestrial planet within their habitable zone, consistent with the value derived from Kepler observations. On the other hand, the habitability of passive galaxies is slightly but systematically higher, unless we assume an unrealistically high sensitivity of planets to supernovae. We find that the overall habitability of galaxies has not changed significantly in the last ~8 Gyr, with most of the habitable planets in local disk galaxies having formed ~1.5 Gyr before our own solar system. Finally, we expect that ~1.4e9 planets similar to present-day Earth have existed so far in our galaxy.

We develop a galaxy assignment scheme that populates dark matter halos with galaxies by tracing the most bound member particles (MBPs) of simulated halos. Several merger-timescale models based on analytic calculations and numerical simulations are adopted as the survival time of mock satellite galaxies. We build mock galaxy samples from halo merger data of the Horizon Run 4 $N$-body simulation from $z = 12-0$. We compare group properties and two-point correlation functions (2pCFs) of mock galaxies with those of volume-limited SDSS galaxies, with $r$-band absolute magnitudes of $\mathcal{M}_r - 5 \log h < -21$ and $-20$ at $z=0$. It is found that the MBP-galaxy correspondence scheme reproduces the observed population of SDSS galaxies in massive galaxy groups ($M > 10^{14} h^{-1} M_{\odot}$) and the small-scale 2pCF ($r_{\rm p} < 10 h^{-1} {\rm Mpc}$) quite well for the majority of the merger timescale models adopted. The new scheme outperforms the previous subhalo-galaxy correspondence scheme by more than $2\sigma$.

We present new observations of the field containing the z=3.786 protocluster, PC217.96+32.3. We confirm that it is one of the largest and most overdense high-redshift structures known. Such structures are rare even in the largest cosmological simulations. We used the Mayall/MOSAIC1.1 imaging camera to image a 1.2x0.6 deg area (~150x75 comoving Mpc) surrounding the protocluster's core and discovered 165 candidate Lyman Alpha emitting galaxies (LAEs) and 788 candidate Lyman Break galaxies (LBGs). There are at least 2 overdense regions traced by the LAEs, the largest of which shows an areal overdensity in its core (i.e., within a radius of 2.5 comoving Mpc) of 14+/-7 relative to the average LAE spatial density in the imaged field. Further, the average LAE spatial density in the imaged field is twice that derived by other field LAE surveys. Spectroscopy with Keck/DEIMOS yielded redshifts for 164 galaxies (79 LAEs and 85 LBGs); 65 lie at a redshift of 3.785+/-0.010. The velocity dispersion of galaxies near the core is 350+/-40 km/s, a value robust to selection effects. The overdensities are likely to collapse into systems with present-day masses of >10^{15} solar masses and >6x10^{14} solar masses. The low velocity dispersion may suggest a dynamically young protocluster. We find a weak trend between narrow-band (Lyman Alpha) luminosity and environmental density: the Lyman Alpha luminosity is enhanced on average by 1.35X within the protocluster core. There is no evidence that the Lyman Alpha equivalent width depends on environment. These suggest that star-formation and/or AGN activity is enhanced in the higher density regions of the structure. PC217.96+32.3 is a Coma cluster analog, witnessed in the process of formation.

The concept of the cosmic web, viewing the Universe as a set of discrete galaxies held together by gravity, is deeply engrained in cosmology. Yet, little is known about the most effective construction and the characteristics of the underlying network. Here we explore seven network construction algorithms that use various galaxy properties, from their location, to their size and relative velocity, to assign a network to galaxy distributions provided by both simulations and observations. We find that a model relying only on spatial proximity offers the best correlations between the physical characteristics of the connected galaxies. We show that the properties of the networks generated from simulations and observations are identical, unveiling a deep universality of the cosmic web.

(abridged) The large-scale distribution of galaxies is generally analyzed using the two-point correlation function. However, this statistic does not capture the topology of the distribution, and it is necessary to resort to higher order correlations to break degeneracies. We demonstrate that an alternate approach using network analysis can discriminate between topologically different distributions that have similar two-point correlations. We investigate two galaxy point distributions, one produced by a cosmological simulation and the other by a L\'evy walk. For the cosmological simulation, we adopt the redshift $z = 0.58$ slice from Illustris (Vogelsberger et al. 2014A) and select galaxies with stellar masses greater than $10^8$$M_\odot$. The two point correlation function of these simulated galaxies follows a single power-law, $\xi(r) \sim r^{-1.5}$. Then, we generate L\'evy walks matching the correlation function and abundance with the simulated galaxies. We find that, while the two simulated galaxy point distributions have the same abundance and two point correlation function, their spatial distributions are very different; most prominently, \emph{filamentary structures}, absent in L\'evy fractals. To quantify these missing topologies, we adopt network analysis tools and measure diameter, giant component, and transitivity from networks built by a conventional friends-of-friends recipe with various linking lengths. Unlike the abundance and two point correlation function, these network quantities reveal a clear separation between the two simulated distributions; therefore, the galaxy distribution simulated by Illustris is not a L\'evy fractal quantitatively. We find that the described network quantities offer an efficient tool for discriminating topologies and for comparing observed and theoretical distributions.

We consider theories where dark matter is composed of a thermal relic of weak scale mass, whose couplings to the Standard Model (SM) are however too small to give rise to the observed abundance. Instead, the abundance is set by annihilation to light hidden sector states that carry no charges under the SM gauge interactions. In such a scenario the constraints from direct and indirect detection, and from collider searches for dark matter, can easily be satisfied. The masses of such light hidden states can be protected by symmetry if they are Nambu-Goldstone bosons, fermions, or gauge bosons. These states can then contribute to the cosmic energy density as dark radiation, leading to observable signals in the cosmic microwave background (CMB). Furthermore, depending on whether or not the light hidden sector states self-interact, the fraction of the total energy density that free-streams is either decreased or increased, leading to characteristic effects on both the scalar and tensor components of the CMB anisotropy that allows these two cases to be distinguished. The magnitude of these signals depends on the number of light degrees of freedom in the hidden sector, and on the temperature at which it kinetically decouples from the SM. We consider a simple model that realizes this scenario, based on a framework in which the SM and hidden sector are initially in thermal equilibrium through the Higgs portal, and show that the resulting signals are compatible with recent Planck results, while large enough to be detected in upcoming experiments such as CMBPol and CMB Stage-IV. Invisible decays of the Higgs into hidden sector states at colliders can offer a complementary probe of this model.

Radio relics detected in the outskirts of galaxy clusters are thought to trace radio-emitting relativistic electrons accelerated at cosmological shocks. In this study, using the cosmological hydrodynamic simulation data for the large-scale structure formation and adopting a diffusive shock acceleration (DSA) model for the production of cosmic-ray (CR) electrons, we construct mock radio and X-ray maps of simulated galaxy clusters that are projected in the sky plane. Various properties of shocks and radio relics, including the shock Mach number, radio spectral index and luminosity are extracted from the synthetic maps and compared with observations. A substantial fraction of radio and X-ray shocks identified in these maps involve multiple shock surfaces along line of sights (LoSs), and the morphology of shock distributions in the maps depends on the projection direction. Among multiple shocks in a given LoS, radio observations tend to pick up stronger shocks with flatter radio spectra, while X-ray observations preferentially select weaker shocks with larger kinetic energy flux. As a result, the shock Mach numbers and locations derived from radio and X-ray observations could differ from each other in some cases. We also find that the distributions of the spectral index and radio power of the synthetic radio relics are somewhat inconsistent with those of observed real relics; a bit more radio relics have been observed closer to the cluster core and with steeper spectral indices. We suggest the inconsistency could be explained, if very weak shocks with $M_s \la 2$ accelerate CR electrons more efficiently, compared to the DSA model adopted here.

The Horizon Run 4 is a cosmological $N$-body simulation designed for the study of coupled evolution between galaxies and large-scale structures of the Universe, and for the test of galaxy formation models. Using $6300^3$ gravitating particles in a cubic box of $L_{\rm box} = 3150 ~h^{-1}{\rm Mpc}$, we build a dense forest of halo merger trees to trace the halo merger history with a halo mass resolution scale down to $M_s = 2.7 \times 10^{11} h^{-1}{\rm M_\odot}$. We build a set of particle and halo data, which can serve as testbeds for comparison of cosmological models and gravitational theories with observations. We find that the FoF halo mass function shows a substantial deviation from the universal form with tangible redshift evolution of amplitude and shape. At higher redshifts, the amplitude of the mass function is lower, and the functional form is shifted toward larger values of $\ln (1/\sigma)$. We also find that the baryonic acoustic oscillation feature in the two-point correlation function of mock galaxies becomes broader with a peak position moving to smaller scales and the peak amplitude decreasing for increasing directional cosine $\mu$ compared to the linear predictions. From the halo merger trees built from halo data at 75 redshifts, we measure the half-mass epoch of halos and find that less massive halos tend to reach half of their current mass at higher redshifts. Simulation outputs including snapshot data, past lightcone space data, and halo merger data are available at http://sdss.kias.re.kr/astro/Horizon-Run4/.

We apply simple analyses techniques developed for the study of complex networks to the study of the cosmic web, the large scale galaxy distribution. In this paper, we measure three network centralities (ranks of topological importance), Degree Centrality (DC), Closeness Centrality (CL), and Betweenness Centrality (BC) from a network built from the Cosmological Evolution Survey (COSMOS) catalog. We define 8 galaxy populations according to the centrality measures; Void, Wall, and Cluster by DC, Main Branch and Dangling Leaf by BC, and Kernel, Backbone, and Fracture by CL. We also define three populations by voronoi tessellation density to compare these with the DC selection. We apply the topological selections to galaxies in the (photometric) redshift range $0.91<z<0.94$ from the COSMOS survey, and explore whether the red and blue galaxy populations show differences in color, star-formation rate (SFR) and stellar mass in the different topological regions. Despite the limitations and uncertainties associated with using photometric redshift and indirect measurements of galactic parameters, the preliminary results illustrate the potential of network analysis. The coming future surveys will provide better statistical samples to test and improve this "network cosmology".

The TAIPAN instrument, currently being developed for the Australian Astronomical Observatory's UK Schmidt telescope at Siding Spring Observatory, makes use of the AAO's Starbug technology to deploy 150 science fibres to target positions on the optical plane. This paper describes the software system for controlling and deploying the fibre-bearing Starbug robots. The TAIPAN software is responsible for allocating each Starbug to its next target position based on its current position and the distribution of targets, finding a collision-free path for each Starbug, and then simultaneously controlling the Starbug hardware in a closed loop, with a metrology camera used to determine the position of each Starbug in the field during reconfiguration. The software is written in C++ and Java and employs a DRAMA middleware layer (Farrell et al. 1995).

Starbugs are miniaturised robotic devices that position optical fibres over a telescope's focal plane in parallel operation for high multiplex spectroscopic surveys. The key advantage of the Starbug positioning system is its potential to configure fields of hundreds of targets in a few minutes, consistent with typical detector readout times. Starbugs have been selected as the positioning technology for the TAIPAN (Transforming Astronomical Imaging surveys through Polychromatic Analysis of Nebulae) instrument, a prototype for MANIFEST (Many Instrument Fiber System) on the GMT (Giant Magellan Telescope). TAIPAN consists of a 150-fibre Starbug positioner accessing the 6 degree field-of-view of the AAO's UK Schmidt Telescope at Siding Spring Observatory. For TAIPAN, it is important to optimise the target allocation and routing algorithms to provide the fastest configurations times. We present details of the algorithms and results of the simulated performance.

A novel method to characterize the topology of the early-universe intergalactic medium during the epoch of cosmic reionization is presented. The 21-cm radiation background from high redshift is analyzed through the calculation of the 2-dimensional (2D) genus. The radiative transfer of hydrogen-ionizing photons and ionization-rate equations are calculated in a suite of numerical simulations under various input parameters. The 2D genus is calculated from the mock 21-cm images of the high-redshift Universe. We construct the 2D genus curve by varying the threshold differential brightness temperature, and compare this to the 2D genus curve of the underlying density field. We find that (1) the 2D genus curve reflects the evolutionary track of cosmic reionization and (2) the 2D genus curve can discriminate between certain reionization scenarios and thus indirectly probe the properties of radiation-sources. Choosing the right beam shape of a radio antenna is crucial for this analysis. To this end, the Square Kilometre Array (SKA) is found to be a suitable apparatus for this analysis in terms of sensitivity, even though some deterioration of the data for this purpose is unavoidable under the planned size of the antenna core.

We report the discovery of a large-scale structure containing multiple protoclusters at z=3.78 in the Bo\"otes field. The spectroscopic discovery of five galaxies at z=3.783+/-0.002 lying within 1 Mpc of one another led us to undertake a deep narrow- and broad-band imaging survey of the surrounding field. Within a comoving volume of 72x72x25 Mpc^3, we have identified 65 Lyman alpha emitter (LAE) candidates at z=3.795+/-0.015, and four additional galaxies at z_spec=3.730,3.753,3.780,3.835. The galaxy distribution within the field is highly non-uniform, exhibiting three large (~3-5x) overdensities separated by 8-14 Mpc (physical) and possibly connected by filamentary structures traced by LAEs. The observed number of LAEs in the entire field is nearly twice the average expected in field environments, based on estimates of the Lya luminosity function at these redshifts. We estimate that by z=0 the largest overdensity will grow into a cluster of mass 10^15 Msun; the two smaller overdensities will grow into clusters of mass (2-6)x10^14 Msun. The highest concentration of galaxies is located at the southern end of the image, suggesting that the current imaging may not map the true extent of the large scale structure. Finding three large protocluster candidates within a single 0.3 deg^2 field is highly unusual; expectations from theory suggest that such alignments should occur less than 2% of the time. Searching for and characterizing such structures and accurately measuring their volume space density can therefore place constraints on the theory of structure formation. Such regions can also serve as laboratories for the study of galaxy formation in dense environments.

The outskirts of galaxy clusters are continuously disturbed by mergers and gas infall along filaments, which in turn induce turbulent flow motions and shock waves. We examine the properties of shocks that form within $r_{200}$ in sample galaxy clusters from structure formation simulations. While most of these shocks are weak and inefficient accelerators of cosmic rays (CRs), there are a number of strong, energetic shocks which can produce large amounts of CR protons via diffusive shock acceleration. We show that the energetic shocks reside mostly in the outskirts and a substantial fraction of them are induced by infall of the warm-hot intergalactic medium from filaments. As a result, the radial profile of the CR pressure in the intracluster medium is expected to be broad, dropping off more slowly than that of the gas pressure, and might be even temporarily inverted, peaking in the outskirts. The volume-integrated momentum spectrum of CR protons inside $r_{200}$ has the power-law slope of $4.25 - 4.5$, indicating that the average Mach number of the shocks of main CR production is in the range of $\left< M_s \right>_{\rm CR} \approx 3 - 4$. We suggest that some radio relics with relatively flat radio spectrum could be explained by primary electrons accelerated by energetic infall shocks with $M_s >~ 3$ induced in the cluster outskirts.

We present HST/WFC3 narrow-band imaging of the starburst galaxy M83 targeting the hydrogen recombination lines (H$\beta$, H$\alpha$ and Pa$\beta$), which we use to investigate the dust extinction in the HII regions. We derive extinction maps with 6 parsec spatial resolution from two combinations of hydrogen lines (H$\alpha$/H$\beta$ and H$\alpha$/Pa$\beta$), and show that the longer wavelengths probe larger optical depths, with $A_V$ values larger by $\gtrsim$1 mag than those derived from the shorter wavelengths. This difference leads to a factor $\gtrsim$2 discrepancy in the extinction-corrected H$\alpha$ luminosity, a significant effect when studying extragalactic HII regions. By comparing these observations to a series of simple models, we conclude that a large diversity of absorber/emitter geometric configurations can account for the data, implying a more complex physical structure than the classical foreground "dust screen" assumption. However, most data points are bracketed by the foreground screen and a model where dust and emitters are uniformly mixed. When averaged over large ($\gtrsim$100--200 pc) scales, the extinction becomes consistent with a "dust screen", suggesting that other geometries tend to be restricted to more local scales. Moreover, the extinction in any region can be described by a combination of the foreground screen and the uniform mixture model with weights of 1/3 and 2/3 in the center ($\lesssim$2 kpc), respectively, and 2/3 and 1/3 for the rest of the disk. This simple prescription significantly improves the accuracy of the dust extinction corrections and can be especially useful for pixel-based analyses of galaxies similar to M83.

We present an objective method to remove the stellar continuum emission from narrow-band images to derive emission-line images. The method is based on the skewness of the pixel histogram of the residual images. Specifically, we exploit a transition in the skewness of the signal in the continuum-subtracted image, which appears when the image changes from being under-subtracted to over-subtracted. Tests on one-dimensional artificial images demonstrate that the transition identifies the optimal scaling factor {\mu} to be used on the broad-band image IB in order to produce the optimal line-emission image IE, i.e., IE =IN - {\mu} IB, with IN the original (un-subtracted) narrow-band image. The advantage of this method is that it uses all information-bearing pixels in the final image, and not just a sub-set of those pixels (the latter being common in many traditional approaches to stellar continuum removal from narrow-band images). We apply our method to actual images, both from ground-based and space facilities, in particular to WFPC2 and ACS images from the Hubble Space Telescope, and we show that it is successful irrespective of the nature of the sources (point-like or extended). We also discuss the impact on the accuracy of the method of non-optimal images, such as those containing saturated sources or non-uniform background, and present `workarounds' for those problems.

Modern spectroscopic surveys produce large spectroscopic databases, generally with sizes well beyond the scope of manual investigation. The need arises, therefore, for an automated line detection method with objective indicators for detection significance. In this paper, we present an automated and objective method for emission line detection in spectroscopic surveys and apply this technique to 1574 spectra, obtained with the Hectospec spectrograph on the MMT Observatory (MMTO), to detect Lyman alpha emitters near z ~ 2.7. The basic idea is to generate on-source (signal plus noise) and off-source (noise only) mock observations using Monte Carlo simulations, and calculate completeness and reliability values, (C, R), for each simulated signal. By comparing the detections from real data with the Monte Carlo results, we assign the completeness and reliability values to each real detection. From 1574 spectra, we obtain 881 raw detections and, by removing low reliability detections, we finalize 649 detections from an automated pipeline. Most of high completeness and reliability detections, (C, R) ~ (1.0, 1.0), are robust detections when visually inspected; the low C and R detections are also marginal on visual inspection. This method at detecting faint sources is dependent on the accuracy of the sky subtraction.

(abridged) We investigate the properties of feedback-driven shocks in 8 nearby starburst galaxies using narrow-band imaging data from the Hubble Space Telescope (HST). We identify the shock--ionized component via the line diagnostic diagram \oiii/\hb vs. \sii (or \nii)/\ha, applied to resolved regions 3--15 pc in size. We divide our sample into three sub-samples: sub-solar (Holmberg II, NGC 1569, NGC 4214, NGC 4449, and NGC 5253), solar (He 2-10, NGC 3077) and super-solar (NGC 5236) for consistent shock measurements. For the sub-solar sub-sample, we derive three scaling relations: (1) $L_{shock} \propto {SFR}^{~0.62}$, (2) $L_{shock} \propto {\Sigma_{SFR,HL}}^{~0.92}$, and (3) $L_{shock}/L_{tot} \propto {(L_H/L_{\odot,H})}^{-0.65}$, where $L_{shock}$ is the \ha luminosity from shock--ionized gas, ${\Sigma_{SFR,HL}}$ the SFR per unit half-light area, $L_{tot}$ the total \ha luminosity, and $L_H/L_{\odot,H}$ the absolute H-band luminosity from 2MASS normalized to solar luminosity. The other two sub--samples do not have enough number statistics, but appear to follow the first scaling relation. The energy recovered indicates that the shocks from stellar feedback in our sample galaxies are fully radiative. If the scaling relations are applicable in general to stellar feedback, our results are similar to those by Hopkins et al. (2012) for galactic super winds. This similarity should, however, be taken with caution at this point, as the underlying physics that enables the transition from radiative shocks to gas outflows in galaxies is still poorly understood.

In order to understand the collapse dynamics of observed low-mass starless cores, we revise the conventional stability condition of hydrostatic Bonnor-Ebert spheres to take internal motions into account. Because observed starless cores resemble Bonnor-Ebert density structures, the stability and dynamics of the starless cores are frequently analyzed by comparing to the conventional stability condition of a hydrostatic Bonnor-Ebert sphere. However, starless cores are not hydrostatic but have observed internal motions. In this study, we take gaseous spheres with a homologous internal velocity field and derive stability conditions of the spheres utilizing a virial analysis. We propose two limiting models of spontaneous gravitational collapse: the collapse of critical Bonnor-Ebert spheres and uniform density spheres. The collapse of these two limiting models are intended to provide the lower and the upper limits, respectively, of the infall speeds for a given density structure. The results of our study suggest that the stability condition sensitively depends on internal motions. A homologous inward motion with a transonic speed can reduce the critical size compared to the static Bonnor-Ebert sphere by more than a factor of two. As an application of the two limiting models of spontaneous gravitational collapse, we compare the density structures and infall speeds of the observed starless cores L63, L1544, L1689B, and L694-2 to the two limiting models. L1689B and L694-2 seem to have been perturbed to result in faster infall motions than for spontaneous gravitational collapse.

We investigate gravitational instability (GI) of rotating, vertically-stratified, pressure-confined, polytropic gas disks using linear stability analysis as well as analytic approximations. The disks are initially in vertical hydrostatic equilibrium and bounded by a constant external pressure. We find that GI of a pressure-confined disk is in general a mixed mode of the conventional Jeans and distortional instabilities, and is thus an unstable version of acoustic-surface-gravity waves. The Jeans mode dominates in weakly confined disks or disks with rigid boundaries. When the disk has free boundaries and is strongly pressure-confined, on the other hand, the mixed GI is dominated by the distortional mode that is surface-gravity waves driven unstable under own gravity and thus incompressible. We demonstrate that the Jeans mode is gravity-modified acoustic waves rather than inertial waves and that inertial waves are almost unaffected by self-gravity. We derive an analytic expression for the effective sound speed c_eff of acoustic-surface-gravity waves. We also find expressions for the gravity reduction factors relative to a razor-thin counterpart, appropriate for the Jeans and distortional modes. The usual razor-thin dispersion relation after correcting for c_eff and the reduction factors closely matches the numerical results obtained by solving a full set of linearized equations. The effective sound speed generalizes the Toomre stability parameter of the Jeans mode to allow for the mixed GI of vertically-stratified, pressure-confined disks.

Weinberg et al. calculated the anthropic likelihood of the cosmological constant using a model assuming that the number of observers is proportional to the total mass of gravitationally collapsed objects, with mass greater than a certain threshold, at t \rightarrow \infty. We argue that Weinberg's model is biased toward small \Lambda, and to try to avoid this bias we modify his model in a way that the number of observers is proportional to the number of collapsed objects, with mass and time equal to certain preferred mass and time scales. Compared to Weinberg's model, this model gives a lower anthropic likelihood of \Lambda_0 (T_+(\Lambda_0) ~ 5%). On the other hand, the anthropic likelihood of the primordial density perturbation amplitude from this model is high, while the likelihood from Weinberg's model is low. Furthermore, observers will be affected by the history of the collapsed object, and we introduce a method to calculate the anthropic likelihoods of \Lambda and Q from the mass history using the extended Press-Schechter formalism. The anthropic likelihoods for $\Lambda$ and Q from this method are similar to those from our single mass constraint model, but, unlike models using the single mass constraint which always have degeneracies between \Lambda and Q, the results from models using the mass history are robust even if we allow both \Lambda and Q to vary. In the case of Weinberg's flat prior distribution of \Lambda (pocket based multiverse measure), our mass history model gives T_+(\Lambda_0) ~ 10%, while the scale factor cutoff measure and the causal patch measure give T_+(\Lambda_0) \geq 30%.

We investigate the ionization structure of the nebular gas in M83 using the line diagnostic diagram, [O III](5007 \degA)/H{\beta} vs. [S II](6716 \deg A+6731 \deg A)/H{\alpha} with the newly available narrowband images from the Wide Field Camera 3 (WFC3) of the Hubble Space Telescope (HST). We produce the diagnostic diagram on a pixel-by-pixel (0.2" x 0.2") basis and compare it with several photo- and shock-ionization models. For the photo-ionized gas, we observe a gradual increase of the log([O III]/H{\beta}) ratios from the center to the spiral arm, consistent with the metallicity gradient, as the H II regions go from super solar abundance to roughly solar abundance from the center out. Using the diagnostic diagram, we separate the photo-ionized from the shock-ionized component of the gas. We find that the shock-ionized H{\alpha} emission ranges from ~2% to about 15-33% of the total, depending on the separation criteria used. An interesting feature in the diagnostic diagram is an horizontal distribution around log([O III]/H{\beta}) ~ 0. This feature is well fit by a shock-ionization model with 2.0 Z\odot metallicity and shock velocities in the range of 250 km/s to 350 km/s. A low velocity shock component, < 200 km/s, is also detected, and is spatially located at the boundary between the outer ring and the spiral arm. The low velocity shock component can be due to : 1) supernova remnants located nearby, 2) dynamical interaction between the outer ring and the spiral arm, 3) abnormal line ratios from extreme local dust extinction. The current data do not enable us to distinguish among those three possible interpretations. Our main conclusion is that, even at the HST resolution, the shocked gas represents a small fraction of the total ionized gas emission at less than 33% of the total. However, it accounts for virtually all of the mechanical energy produced by the central starburst in M83.

We use Spitzer 24 $\mu$m, 70 $\mu$m and ground based H$\alpha$ data for a sample of 40 SINGS galaxies to establish a star formation rate (SFR) indicator using 70 $\mu$m emission for sub--galactic ($\sim0.05-2\ \rm{kpc}$) line-emitting regions and to investigate limits in application. A linear correlation between 70 $\mu$m and SFR is found and a star formation indicator SFR(70) is proposed for line-emitting sub-galactic regions as $\rm \Sigma(SFR)\ ({M_{\odot}\cdot yr^{-1}\cdot kpc^{-2}})=9.4\times10^{-44}\ \Sigma(70)\ \rm{(ergs\cdot s^{-1}\cdot kpc^{-2})}$, for regions with $12+\rm{log(O/H)}\gtrsim8.4$ and $\rm \Sigma(SFR)\gtrsim10^{-3}\ (M_{\odot}\cdot yr^{-1}\cdot kpc^{-2})$, with a 1-$\sigma$ dispersion around the calibration of $\sim0.16$ dex. We also discuss the influence of metallicity on the scatter of the data. Comparing with the SFR indicator at 70 $\mu$m for integrated light from galaxies, we find that there is $\sim40%$ excess 70 $\mu$m emission in galaxies, which can be attributed to stellar populations not involved in the current star formation activity.

The kinematics of damped Lyman alpha absorbers (DLAs) are difficult to reproduce in hierarchical galaxy formation models, particularly the preponderance of wide systems. We investigate DLA kinematics at z=3 using high-resolution cosmological hydrodynamical simulations that include a heuristic model for galactic outflows. Without outflows, our simulations fail to yield enough wide DLAs, as in previous studies. With outflows, predicted DLA kinematics are in much better agreement with observations. Comparing two outflow models, we find that a model based on momentum-driven wind scalings provides the best match to the observed DLA kinematic statistics of Prochaska & Wolfe. In this model, DLAs typically arise a few kpc away from galaxies that would be identified in emission. Narrow DLAs can arise from any halo and galaxy mass, but wide ones only arise in halos with mass >10^11 Mo, from either large central or small satellite galaxies. This implies that the success of this outflow model originates from being most efficient at pushing gas out from small satellite galaxies living in larger halos. This increases the cross-section for large halos relative to smaller ones, thereby yielding wider kinematics. Our simulations do not include radiative transfer effects or detailed metal tracking, and outflows are modeled heuristically, but they strongly suggest that galactic outflows are central to understanding DLA kinematics. An interesting consequence is that DLA kinematics may place constraints on the nature and efficiency of gas ejection from high-z galaxies.

We present a study of a large, statistically complete sample of star-forming dwarf galaxies using mid-infrared observations from the {\it Spitzer Space Telescope}. The relationships between metallicity, star formation rate (SFR) and mid-infrared color in these systems show that the galaxies span a wide range of properties. However, the galaxies do show a deficit of 8.0 \um\ polycyclic aromatic hydrocarbon emission as is apparent from the median 8.0 \um\ luminosity which is only 0.004 \lstarf\ while the median $B$-band luminosity is 0.05 \lstarb. Despite many of the galaxies being 8.0 \um\ deficient, there is about a factor of 4 more extremely red galaxies in the [3.6] $-$ [8.0] color than for a sample of normal galaxies with similar optical colors. We show correlations between the [3.6] $-$ [8.0] color and luminosity, metallicity, and to a lesser extent SFRs that were not evident in the original, smaller sample studied previously. The luminosity--metallicity relation has a flatter slope for dwarf galaxies as has been indicated by previous work. We also show a relationship between the 8.0 \um\ luminosity and the metallicity of the galaxy which is not expected given the competing effects (stellar mass, stellar population age, and the hardness of the radiation field) that influence the 8.0 \um\ emission. This larger sample plus a well-defined selection function also allows us to compute the 8.0 \um\ luminosity function and compare it with the one for the local galaxy population. Our results show that below 10$^{9}$ $L$\solar, nearly all the 8.0 \um\ luminosity density of the local universe arises from dwarf galaxies that exhibit strong \ha\ emission -- i.e., 8.0 \um\ and \ha\ selection identify similar galaxy populations despite the deficit of 8.0 \um\ emission observed in these dwarfs.

(Abridged) Spitzer data at 24, 70, and 160 micron and ground-based H-alpha images are analyzed for a sample of 189 nearby star-forming and starburst galaxies to investigate whether reliable star formation rate (SFR) indicators can be defined using the monochromatic infrared dust emission centered at 70 and 160 micron. We compare recently published recipes for SFR measures using combinations of the 24 micron and observed H-alpha luminosities with those using 24 micron luminosity alone. From these comparisons, we derive a reference SFR indicator for use in our analysis. Linear correlations between SFR and the 70 and 160 micron luminosity are found for L(70)>=1.4x10^{42} erg/s and L(160)>=2x10^{42} erg/s, corresponding to SFR>=0.1-0.3 M_sun/yr. Below those two luminosity limits, the relation between SFR and 70 micron (160 micron) luminosity is non-linear and SFR calibrations become problematic. The dispersion of the data around the mean trend increases for increasing wavelength, becoming about 25% (factor ~2) larger at 70 (160) micron than at 24 micron. The increasing dispersion is likely an effect of the increasing contribution to the infrared emission of dust heated by stellar populations not associated with the current star formation. The non-linear relation between SFR and the 70 and 160 micron emission at faint galaxy luminosities suggests that the increasing transparency of the interstellar medium, decreasing effective dust temperature, and decreasing filling factor of star forming regions across the galaxy become important factors for decreasing luminosity. The SFR calibrations are provided for galaxies with oxygen abundance 12+Log(O/H)>8.1. At lower metallicity the infrared luminosity no longer reliably traces the SFR because galaxies are less dusty and more transparent.

Compared to starburst galaxies, normal star forming galaxies have been shown to display a much larger dispersion of the dust attenuation at fixed reddening through studies of the IRX-beta diagram (the IR/UV ratio "IRX" versus the UV color "beta"). To investigate the causes of this larger dispersion and attempt to isolate second parameters, we have used GALEX UV, ground-based optical, and Spitzer infrared imaging of 8 nearby galaxies, and examined the properties of individual UV and 24 micron selected star forming regions. We concentrated on star-forming regions, in order to isolate simpler star formation histories than those that characterize whole galaxies. We find that 1) the dispersion is not correlated with the mean age of the stellar populations, 2) a range of dust geometries and dust extinction curves are the most likely causes for the observed dispersion in the IRX-beta diagram 3) together with some potential dilution of the most recent star-forming population by older unrelated bursts, at least in the case of star-forming regions within galaxies, 4) we also recover some general characteristics of the regions, including a tight positive correlation between the amount of dust attenuation and the metal content. Although generalizing our results to whole galaxies may not be immediate, the possibility of a range of dust extinction laws and geometries should be accounted for in the latter systems as well.

The long, bright gamma-ray burst GRB 070125 was localized by the Interplanetary Network. We present light curves of the prompt gamma-ray emission as observed by Konus-WIND, RHESSI, Suzaku-WAM, and \textit{Swift}-BAT. We detail the results of joint spectral fits with Konus and RHESSI data. The burst shows moderate hard-to-soft evolution in its multi-peaked emission over a period of about one minute. The total burst fluence as observed by Konus is $1.79 \times 10^{-4}$ erg/cm$^2$ (20 keV--10 MeV). Using the spectroscopic redshift $z=1.548$, we find that the burst is consistent with the ``Amati'' $E_{peak,i}-E_{iso}$ correlation. Assuming a jet opening angle derived from broadband modeling of the burst afterglow, GRB 070125 is a significant outlier to the ``Ghirlanda'' $E_{peak,i}-E_\gamma$ correlation. Its collimation-corrected energy release $E_\gamma = 2.5 \times 10^{52}$ ergs is the largest yet observed.

The Wide-band All-sky Monitor (WAM) is a function of the large lateral BGO shield of the Hard X-ray Detector (HXD) onboard Suzaku. Its large geometrical area of 800 cm^2 per side, the large stopping power for the hard X-rays and the wide-field of view make the WAM an ideal detector for gamma-ray bursts (GRBs) observations in the energy range of 50-5000 keV. In fact, the WAM has observed 288 GRBs confirmed by other satellites, till the end of May 2007.

The proton-proton fusion reaction, $pp\to de^+\nu$, is studied in pionless effective field theory (EFT) with di-baryon fields up to next-to leading order. With the aid of the di-baryon fields, the effective range corrections are naturally resummed up to the infinite order and thus the calculation is greatly simplified. Furthermore, the low-energy constant which appears in the axial-current-di-baryon-di-baryon contact vertex is fixed through the ratio of two- and one-body matrix elements which reproduces the tritium lifetime very precisely. As a result we can perform a parameter free calculation for the process. We compare our numerical result with those from the accurate potential model and previous pionless EFT calculations, and find a good agreement within the accuracy better than 1%.

We report the discovery of a star-forming loop around the young, Crab-like supernova remnant (SNR) G54.1+0.3 using the AKARI infrared satellite. The loop consists of at least eleven young stellar objects (YSOs) embedded in a ring-like diffuse emission of radius ~1'. The YSOs are bright in the mid-infrared and are also visible in the Spitzer Space Telescope Galactic plane survey images. Their Spitzer colors are similar to those of class II YSOs in [3.6]-[5.8] but significantly redder in [8]-[24], i.e., 0<[3.6]-[5.8]<1.2 and 5<[8]-[24]<9. Most of them have near-infrared counterparts in the 2MASS JHKs images, and some of them have an optical counterpart too. Their JHKs colors and magnitudes indicate that the YSOs are massive (<= 10 Msun) pre-main-sequence stars at the same distance to the SNR, i.e., 8 kpc, which supports the association of the star-forming loop with the SNR. The dereddened spectral energy distributions are similar to eraly Herbig Be stars, which are early B-type pre-main-sequence stars with inner disks that have been destroyed. The confinement to a loop structure indicates that the YSOs are young, i.e., <= 2 Myr. We propose that their formation is triggered by the progenitor star of G54.1+0.3, which has a mass of <= 15 Msun. The triggering must have occurred near the end of the progenitor's life, possibly after it had evolved off the main sequence.

AKARI, the first Japanese satellite dedicated to infrared astronomy, was launched on 2006 February 21, and started observations in May of the same year. AKARI has a 68.5 cm cooled telescope, together with two focal-plane instruments, which survey the sky in six wavelength bands from the mid- to far-infrared. The instruments also have the capability for imaging and spectroscopy in the wavelength range 2 - 180 micron in the pointed observation mode, occasionally inserted into the continuous survey operation. The in-orbit cryogen lifetime is expected to be one and a half years. The All-Sky Survey will cover more than 90 percent of the whole sky with higher spatial resolution and wider wavelength coverage than that of the previous IRAS all-sky survey. Point source catalogues of the All-Sky Survey will be released to the astronomical community. The pointed observations will be used for deep surveys of selected sky areas and systematic observations of important astronomical targets. These will become an additional future heritage of this mission.

We present a serendipitous detection of the infrared-bright supernova remnant (SNR) B0104-72.3 in the Small Magellanic Cloud by the Infrared Camera (IRC) onboard AKARI. An elongated, partially complete shell is detected in all four observed IRC bands covering 2.6-15 um. The infrared shell surrounds radio, optical, and X-ray emission associated with the SNR and is probably a radiative SNR shell. This is the first detection of a SNR shell in this near/mid-infrared waveband in the Small Magellanic Cloud. The IRC color indicates that the infrared emission might be from shocked H2 molecules with some possible contributions from ionic lines. We conclude that B0104-72.3 is a middle-aged SNR interacting with molecular clouds, similar to the Galactic SNR IC 443. Our results highlight the potential of AKARI IRC observations in studying SNRs, especially for diagnosing SNR shocks.

The in-orbit performance and calibration of the Hard X-ray Detector (HXD) on board the X-ray astronomy satellite Suzaku are described. Its basic performances, including a wide energy bandpass of 10-600 keV, energy resolutions of ~4 keV (FWHM) at 40 keV and ~11% at 511 keV, and a high background rejection efficiency, have been confirmed by extensive in-orbit calibrations. The long-term gains of PIN-Si diodes have been stable within 1% for half a year, and those of scintillators have decreased by 5-20%. The residual non-X-ray background of the HXD is the lowest among past non-imaging hard X-ray instruments in energy ranges of 15-70 and 150-500 keV. We provide accurate calibrations of energy responses, angular responses, timing accuracy of the HXD, and relative normalizations to the X-ray CCD cameras using multiple observations of the Crab Nebula.

The Hard X-ray Detector (HXD) on board Suzaku covers a wide energy range from 10 keV to 600 keV by combination of silicon PIN diodes and GSO scintillators. The HXD is designed to achieve an extremely low in-orbit back ground based on a combination of new techniques, including the concept of well-type active shield counter. With an effective area of 142 cm^2 at 20 keV and 273 cm2 at 150 keV, the background level at the sea level reached ~1x10^{-5} cts s^{-1} cm^{-2} keV^{-1} at 30 keV for the PI N diodes, and ~2x10^{-5} cts s^{-1} cm^{-2} keV^{-1} at 100 keV, and ~7x10^{-6} cts s^{-1} cm^{-2} keV^{-1} at 200 keV for the phoswich counter. Tight active shielding of the HXD results in a large array of guard counters surrounding the main detector parts. These anti-coincidence counters, made of ~4 cm thick BGO crystals, have a large effective area for sub-MeV to MeV gamma-rays. They work as an excellent gamma-ray burst monitor with limited angular resolution (~5 degree). The on-board signal-processing system and the data transmitted to the ground are also described.

The total cross section for radiative neutron capture on a proton, $np \to d \gamma$, is evaluated at big bang nucleosynthesis (BBN) energies. The electromagnetic transition amplitudes are calculated up to next-to leading order within the framework of pionless effective field theory with dibaryon fields. We also calculate the $d\gamma\to np$ cross section and the photon analyzing power for the $d\vec{\gamma}\to np$ process from the amplitudes. The values of low energy constants that appear in the amplitudes are estimated by a Markov Chain Monte Carlo analysis using the relevant low energy experimental data. Our result agrees well with those of other theoretical calculations except for the $np\to d\gamma$ cross section at some energies estimated by an R-matrix analysis. We also study the uncertainties in our estimation of the $np\to d\gamma$ cross section at relevant BBN energies and find that the estimated cross section is reliable to within $\sim$1% error.

The XMM-Newton observation of the nearby FR-II radio galaxy 3C 98 is reported. In two exposures on the target, faint diffuse X-ray emission associated with the radio lobes was significantly detected, together with a bright X-ray active nucleus, of which the 2 -- 10 keV intrinsic luminosity is (4 -- 8) \times 10^{42} erg s-1. The EPIC spectra of the northern and southern lobes are reproduced by a single power law model modified by the Galactic absorption, with a photon index of 2.2-0.5+0.6 and 1.7-0.6+0.7 respectively. These indices are consistent with that of the radio synchrotron spectrum, 1.73 +- 0.01 The luminosity of the northern and southern lobes are measured to be 8.3-2.6+3.3 \times 10^{40} erg s-1 and 9.2-4.3+5.7 \times 10^{40} erg s-1, respectively, in the 0.7 -- 7 keV range. The diffuse X-ray emission is interpreted as an inverse-Compton emission, produced when the synchrotron-emitting energetic electrons in the lobes scatter off the cosmic microwave background photons. The magnetic field in the lobes is calculated to be about 1.7 \mu G, which is about 2.5 times lower than the value estimated under the minimum energy condition. The energy density of the electrons is inferred to exceed that in the magnetic fields by a factor of 40 -- 50.

An analysis was made of 0.3-15 keV X-ray spectra of a Narrow-Line Seyfert 1 Galaxy, Ton S180, using archival data from ASCA, RXTE, and XMM-Newton. At energies above 2.5 keV, a power-law with a photon index of ~ 2.3 successfully and consistently reproduced the spectra from all of these observatories. Assuming this power-law component to extend toward lower energies, a soft excess, which is one of the most remarkable features of Narrow-Line Seyfert 1 Galaxies, is explained by another power-law multiplied by a thermal cutoff at ~ 0.4 keV. Some similarities have been observed between this object and Galactic black hole binaries in very high state, the latter being realized under high accretion rates. Attempts have been made to interpret the soft excess in terms of Comptonization of the disk photons by an electron cloud surrounding the accretion disk, like BHBs in a very high state.

To be a formation mechanism of such large-scale structures as giant molecular clouds (GMCs) and HI superclouds, the classical Parker instability driven by external gravity has to overcome three major obstacles: The convective motion accompanying the instability generates thin sheets than large condensations. The degree of density enhancement achieved by the instability is too low to make dense interstellar clouds. The time and the length scales of the instability are significantly longer and larger than the estimated formation time and the observed mean separation of the GMCs, respectively. This paper examines whether a replacement of the driving agent from the external to the self gravity might remove these obstacles by activating the gravitational instability in the Galactic ISM disk. The self gravity can suppress the convective motions, and a cooperative action of the Jeans and the Parker instabilities can remove all the obstacles confronting the classical version of the Parker instability. The mass and mean separation of the structures resulting from the odd-parity undular mode solution are shown to agree better with the HI superclouds than with the GMCs. We briefly discuss how inclusions of the external gravity and cosmic rays would modify behaviors of the odd-parity undular mode solution.

The Parker instability, which has been considered as a process governing the structure of the interstellar medium, is induced by the buoyancy of magnetic field and cosmic rays. In previous studies, while the magnetic field has been fully incorporated in the context of isothermal magnetohydrodynamics, cosmic rays have been normally treated with the simplifying assumption of infinite diffusion along magnetic field lines but no diffusion across them. The cosmic ray diffusion is, however, finite. In this work, we take into account fully the diffusion process of cosmic rays in a linear stability analysis of the Parker instability. Cosmic rays are described with the diffusion-convection equation. With realistic values of cosmic ray diffusion coefficients expected in the interstellar medium, we show that the result of previous studies with the simplifying assumption on cosmic ray diffusion applies well. Finiteness of parallel diffusion decreases the growth rate of the Parker instability, while the relatively smaller perpendicular diffusion has no significant effect. We discuss the implication of our result on the role of the Parker instability in the interstellar medium.

Here we present three-dimensional MHD models for the Parker instability in a thick magnetized disk, including the presence of a spiral arm. The $B$-field is assumed parallel to the arm, and the model results are applied to the optical segment of the Carina-Sagittarius arm. The characteristic features of the undular and interchange modes are clearly apparent in the simulations. The undular mode creates large gas concentrations distributed along the arm. This results in a clear arm/inter-arm difference: the instability triggers the formation of large interstellar clouds inside the arms, but generates only small structures with slight density enhancements in the inter-arm regions. The resulting clouds are distributed in an antisymmetric way with respect to the midplane, creating an azimuthal corrugation along the arm. For conditions similar to those of the optical segment of the Carina-Sagittarius arm, it has a wavelength of about 2.4 kpc. This structuring can explain the origin of both HI superclouds and the azimuthal corrugations in spiral arms. The wavelength matches the corrugation length derived with the young stellar groups located in the optical segment of the Carina-Sagittarius arm. Keywords: Galaxy: kinematics and dynamics -- Galaxy: structure -- Instabilities -- ISM: clouds -- ISM: magnetic fields -- ISM: structure -- MHD

We present 2D, ideal-MHD numerical simulations of the Parker instability in a multi-component warm disk model. The calculations were done using two numerical codes with different algorithms, TVD and ZEUS-3D. The outcome of the numerical experiments performed with both codes is very similar, and confirms the results of the linear analysis for the undular mode derived by Kim et al. (2000): the most unstable wavelength is about 3 kpc and its growth timescale is between 30-50 Myr (the growth rate is sensitive to the position of the upper boundary of the numerical grid). Thus, the time and length scales of this multicomponent disk model are substantially larger than those derived for thin disk models. We use three different types of perturbations, random, symmetric, and antisymmetric, to trigger the instability. The antisymmetric mode is dominant, and determines the minimum time for the onset of the nonlinear regime. The instability generates dense condensations and the final peak column density value in the antisymmetric case, as also derived by Kim et al. (2000), is about a factor of 3 larger than its initial value. These wavelengths and density enhancement factors indicate that the instability alone cannot be the main formation mechanism of giant molecular clouds in the general interstellar medium. The role of the instability in the formation of large-scale corrugations along spiral arms is briefly discussed.

A linear stability analysis of a multi-component and magnetized Galactic disk model is presented. The disk model uses the observed stratifications for the gas density and gravitational acceleration at the solar neighborhood and, in this sense, it can be called a realistic model. The distribution of the total gas pressure is defined by these observed stratifications, and the gaseous disk is assumed isothermal. The initial magnetic field is taken parallel to the disk, with a midplane value of 5 $\mu$G, and its stratification along the z-axis is derived from the condition of magnetohydrostatic equilibrium in an isothermal atmosphere. The resulting isothermal sound speed is $\sim 8.4$ km s$^{-1}$, similar to the velocity dispersion of the main gas components within 1.5 kpc from midplane. The thermal-to-magnetic pressure ratio decreases with $[z]$ and the warm model is Parker unstable. The dispersion relations show that the fastest growing mode has a wavelength of about 3 kpc, for both symmetric and antisymmetric perturbations, and the corresponding growth time scales are of about $3\times 10^7$ years. The structure of the final equilibrium stage is also derived, and we find that the midplane antisymmetric (MA) mode gathers more gas in the magnetic valleys. The resulting MA gas condensations have larger densities, and the column density enhancement is a factor of about 3 larger than the value of the initial stage. The unstable wavelengths and growth times for the multi-component disk model are substantially larger than those of a thin disk model, and some of the implications of these results are discussed.

A linear stability analysis has been performed onto a self-gravitating magnetized gas disk bounded by external pressure. The resulting dispersion relation is fully explained by three kinds of instability: a Parker-type instability driven by self-gravity, usual Jeans gravitational instability and convection. In the direction parallel to the magnetic fields, the magnetic tension completely suppresses the convection. If the adiabatic index $\gamma$ is less than a certain critical value, the perturbations trigger the Parker as well as the Jeans instability in the disk. Consequently, the growth rate curve has two maxima: one at small wavenumber due to a combination of the Parker and Jeans instabilities, and the other at somewhat larger wavenumber mostly due to the Parker instability. In the horizontal direction perpendicular to the fields, the convection makes the growth rate increase monotonically upto a limiting value as the perturbation wavenumber gets large. However, at small wavenumbers, the Jeans instability becomes effective and develops a peak in the growth rate curve. Depending on the system parameters, the maximum growth rate of the convection may or may not be higher than the peak due to the Jeans-Parker instability. Therefore, a cooperative action of the Jeans and Parker instabilities can have chances to over-ride the convection and may develop large scale structures of cylindrical shape in non-linear stage. In thick disks the cylinder is expected to align its axis perpendicular to the field, while in thin ones parallel to it.

N-\alpha structures of light nuclei with axial symmetry are studied using relativistic Hartree approximation. Metastable excited states are searched in a configuration space which allows linear alpha chain structures. As a result, it is shown that ^{12}C has ^8Be + \alpha resonance state at about 1 MeV above ^8Be-\alpha threshold as an asymmetric 3-\alpha linear-chain structure, which plays an important role in stellar nucleosynthesis.

We present a multi-dimensional numerical code to solve isothermal magnetohydrodynamic (IMHD) equations for use in modeling astrophysical flows. First, we have built a one-dimensional code which is based on an explicit finite-difference method on an Eulerian grid, called the total variation diminishing (TVD) scheme. Recipes for building the one-dimensional IMHD code, including the normalized right and left eigenvectors of the IMHD Jacobian matrix, are presented. Then, we have extended the one-dimensional code to a multi-dimensional IMHD code through a Strang-type dimensional splitting. In the multi-dimensional code, an explicit cleaning step has been included to eliminate non-zero $\nabla\cdot B$ at every time step. To estimate the proformance of the code, one- and two-dimensional IMHD shock tube tests, and the decay test of a two-dimensional Alfv\'{e}n wave have been done. As an example of astrophysical applications, we have simulated the nonlinear evolution of the two-dimensional Parker instability under a uniform gravity.

Using an isothermal MHD code, we have performed three-dimensional, high-resolution simulations of the Parker instability. The initial equilibrium system is composed of exponentially-decreasing isothermal gas and magnetic field (along the azimuthal direction) under a uniform gravity. The evolution of the instability can be divided into three phases: linear, nonlinear, and relaxed. During the linear phase, the perturbations grow exponentially with a preferred scale along the azimuthal direction but with smallest possible scale along the radial direction, as predicted from linear analyses. During the nonlinear phase, the growth of the instability is saturated and flow motion becomes chaotic. Magnetic reconnection occurs, which allows gas to cross field lines. This, in turn, results in the redistribution of gas and magnetic field. The system approaches a new equilibrium in the relaxed phase, which is different from the one seen in two-dimensional works. The structures formed during the evolution are sheet-like or filamentary, whose shortest dimension is radial. Their maximum density enhancement factor relative to the initial value is less than 2. Since the radial dimension is too small and the density enhancement is too low, it is difficult to regard the Parker instability alone as a viable mechanism for the formation of giant molecular clouds.

To examine how non-uniform nature of the Galactic gravity might affect length and time scales of the Parker instability, we took three models of gravity, uniform, linear and realistic ones. To make comparisons of the three gravity models on a common basis, we first fixed the ratio of magnetic pressure to gas pressure at $\alpha$ = 0.25, that of cosmic-ray pressure at $\beta$ = 0.4, and the rms velocity of interstellar clouds at $a_s$ = 6.4 km s$^{-1}$, and then adjusted parameters of the gravity models in such a way that the resulting density scale heights for the three models may all have the same value of 160 pc. Performing linear stability analyses onto equilibrium states under the three models with the typical ISM conditions, we calculate the maximum growth rate and corresponding length scale for each of the gravity models. Under the uniform gravity the Parker instability has the growth time of 1.2$\times10^{8}$ years and the length scale of 1.6 kpc for symmetric mode. Under the realistic gravity it grows in 1.8$\times10^{7}$ years for both symmetric and antisymmetric modes, and develops density condensations at intervals of 400 pc for the symmetric mode and 200 pc for the antisymmetric one. A simple change of the gravity model has thus reduced the growth time by almost an order of magnitude and its length scale by factors of four to eight. These results suggest that an onset of the Parker instability in the ISM may not necessarily be confined to the regions of high $\alpha$ and $\beta$.

A linear stability analysis has been done to a magnetized disk under a linear gravity. We have reduced the linearized perturbation equations to a second-order differential equation which resembles the Schr\"{o}dinger equation with the potential of a harmonic oscillator. Depending on the signs of energy and potential terms, eigensolutions can be classified into ``continuum'' and ``discrete'' families. When magnetic field is ignored, the continuum family is identified as the convective mode, while the discrete family as acoustic-gravity waves. If the effective adiabatic index $\gamma$ is less than unity, the former develops into the convective instability. When a magnetic field is included, the continuum and discrete families further branch into several solutions with different characters. The continuum family is divided into two modes: one is the original Parker mode, which is a slow MHD mode modulated by the gravity, and the other is a stable Alfv\'en mode. The Parker modes can be either stable or unstable depending on $\gamma$. When $\gamma$ is smaller than a critical value $\gamma_{cr}$, the Parker mode becomes unstable. The discrete family is divided into three modes: a stable fast MHD mode modulated by the gravity, a stable slow MHD mode modulated by the gravity, and an unstable mode which is also attributed to a slow MHD mode. The unstable discrete mode does not always exist. Even though the unstable discrete mode exists, the Parker mode dominates it if the Parker mode is unstable. However, if $\gamma \ge \gamma_{cr}$, the discrete mode could be the only unstable one. When $\gamma$ is equal $\gamma_{cr}$, the minimum growth time of the unstable discrete mode is $1.3 \times 10^8$ years with a corresponding length scale of 2.4 kpc. It is suggestive that the corrugatory features seen in the Galaxy and external galaxies are related to the unstable discrete mode.