We present a comprehensive analysis (photometric and kinematical) of poorly studied open cluster NGC 4337 using 2MASS, WISE, APASS, and Gaia~DR2 database. By determining the membership probabilities of stars, we identified 624 most probable members with membership probability higher than $50\%$ by using proper motion and parallax data taken from Gaia~DR2. The mean proper motion of the cluster is obtained as $\mu_{x}=-8.83\pm0.01$ and $\mu_{y}=1.49\pm0.006$ mas yr$^{-1}$. We find the normal interstellar extinction towards the cluster region. The radial distribution of members provides a cluster radius of 7.75 arcmin (5.63 pc). The estimated age of $1600\pm180$ Myr indicates that NGC 4337 is an old open cluster with a bunch of red giant stars. The overall mass function slope for main-sequence stars is found as $1.46\pm0.18$ within the mass range 0.75$-$2.0 $M_\odot$, which is in fair agreement with Salpeter's value (x=1.35) within uncertainty. The present study demonstrates that NGC 4337 is a dynamically relaxed open cluster. Using the Galactic potential model, Galactic orbits are obtained for NGC 4337. We found that this object follows a circular path around the Galactic center. Under the kinematical analysis, we compute the apex coordinates $(A, D)$ by using two methods: (i) the classical convergent point method and (ii) the AD-diagram method. The obtained coordinates are: $(A_{conv}, D_{conv})$ = (96$^{\textrm{o}}$.27 $\pm$ 0$^{\textrm{o}}$.10, 13$^{\textrm{o}}$.14 $\pm$ 0$^{\textrm{o}}$.27) $\&$ $(A_\circ, D_\circ)$ = (100$^{\textrm{o}}$.282 $\pm$ 0$^{\textrm{o}}$.10, 9$^{\textrm{o}}$.577 $\pm$ 0$^{\textrm{o}}$.323) respectively. We also computed the Velocity Ellipsoid Parameters (VEPs), matrix elements ($\mu_{ij}$), direction cosines ($l_j$, $m_j$, $n_j$) and the Galactic longitude of the vertex ($l_2$).

Observations recently suggest $4.4\sigma$ discrepancy between Hubble constant measured locally and that inferred from cosmic microwave background in standard cosmology. Either standard cosmology or local measurement might expect something new. We investigate the possibility that the value of Hubble constants might be affected by observers' motional status in the local measurement. Using adapted coordinate for geodesic observers in FLRW space-time and constraints inferred from observation of cosmic shear, we find that the motional status of reference frame could contribute to about $1.1 \pm 0.3\%$ discrepancy of these Hubble constants. Namely, the difference is $\Delta H_0= 0.81\pm0.22$ km s$^{-1}$ Mpc$^{-1}$, if the locally measured Hubble constant is $ 74$ km s$^{-1}$ Mpc$^{-1}$. This effect seems not negligible as an uncertainty for local measurement of Hubble constant

ART$^2$ is a 3D multi-wavelength Monte Carlo radiative transfer (RT) code that couples continuum and emission lines to track the propagation of photons and their interactions with the interstellar medium (ISM). The previous version of ART$^2$, which included continuum and Ly$\alpha$ line, has been extensively applied to hydrodynamics simulations to study multi-band properties of galaxies and ISM. Here, we describe new implementations of non-local thermodynamic equilibrium RT of molecular and atomic fine structure emission lines, and the parallelization of the code using a number of novel methods. The new ART$^2$ can efficiently and self-consistently produce a full spectrum that includes both continuum and lines such as [CII], [NII], [OIII], Ly$\alpha$, and CO. These essential features, together with the multi-phase ISM model and the adaptive grid, make ART$^2$ a multi-purpose code to study multi-wavelength properties of a wide range of astrophysical systems from planetary disks to large-scale structures. To demonstrate the capability of the new ART$^2$, we applied it to two hydrodynamics simulations: the zoom-in Milky Way Simulation to obtain panchromatic properties of individual galaxies, and the large-scale IllustrisTNG100 Simulation to obtain global properties such as the line intensity mappings. These products are essential for a broad array of studies such as the correlations between physical and panchromatic properties and their evolution. By enabling direct comparison between numerical simulations and multi-band observations, ART$^2$ provides a crucial theoretical framework for the interpretations of existing observations, the plan for future surveys, and the synergy between multi-band galaxy surveys and line intensity mappings. Therefore, ART$^2$ is a powerful and versatile tool to bridge the gap between theories and observations of cosmic structures.

This paper presents an investigation on the four open clusters Czernik 14, Haffner 14, Haffner 17 and King 10 located near the Perseus arm of Milky Way Galaxy using Gaia DR2, 2MASS, WISE, APASS and Pan-STARRS1 data sets. We find normal interstellar extinction in twelve photometric bands for these clusters. Likely cluster members are identified as 225, 353, 350 and 395 for Czernik 14, Haffner 14, Haffner 17 and King 10, respectively by using Gaia DR2 proper motion data. Radii are determined as 3.5, 3.7, 6.2 and 5.7 arcmin for Czernik 14, Haffner 14, Haffner 17 and King 10 respectively. Mean proper motions in RA and DEC are estimated as (-0.42 \pm 0.02, -0.38 \pm 0.01), (-1.82 \pm 0.009, 1.73 \pm 0.008), (-1.17 \pm 0.007, 1.88 \pm 0.006) and (-2.75 \pm 0.008, -2.04 \pm 0.006) mas/yr for Czernik 14, Haffner 14, Haffner 17 and King 10, respectively. The comparison of observed CMDs with solar metallicity isochrones leads to an age of 570 \pm 60, 320 \pm 35, 90 \pm 10 and 45 \pm 5 Myr for these clusters. The distances 2.9 \pm 0.1, 4.8 \pm 0.4, 3.6 \pm 0.1 and 3.8 \pm 0.1 kpc determined using parallax are comparable with the values derived by the isochrone fitting method. Mass function slopes are found to be in good agreement with the Salpeter value. The total masses are derived as 348, 595, 763 and 1088 solar mass for the clusters Czernik 14, Haffner 14, Haffner 17 and King 10, respectively. Evidence for the existence of mass-segregation effect is observed in each cluster. Using the Galactic potential model, Galactic orbits are derived for the clusters. The present study indicates that all clusters under study fallow a circular path around the Galactic center.

We have introduced a new method of estimating the electron temperature and density of H II regions by using single dish observations. In this method, multiple hydrogen radio recombination lines of different bands are computed under the assumption of low optical depth. We use evolutionary hydrodynamical models of H II regions to model hydrogen recombination line emission from a variety of H II regions and assess the reliability of the method. According to the simulated results, the error of the estimated temperature is commonly < 13%, and that of the estimated density is < 25% for a < 1% uncertainty of the observed line fluxes. A reasonable estimated value of electron density can be achieved if the uncertainty of the line fluxes are lower than 3%. In addition, the estimated values are more representative of the properties in the relatively high-density region if the gas density gradient is present in the H II region. Our method can be independent of the radio continuum observations. But the accuracy will be improved if a line-to-continuum ratio at millimeter wavelengths is added to the estimation. Our method provides a way to measure the temperature and density in ionized regions without interferometers.

Besides buckminsterfullerene (C60), other fullerenes and their derivatives may also reside in space. In this work, we study the formation and photo-dissociation processes of astronomically relevant fullerene/anthracene (C14H10) cluster cations in the gas phase. Experiments are carried out using a quadrupole ion trap (QIT) in combination with time-of-flight (TOF) mass spectrometry. The results show that fullerene (C60, and C70)/anthracene (i.e., [(C14H10)nC60]+ and [(C14H10)nC70]+), fullerene (C56 and C58)/anthracene (i.e., [(C14H10)nC56]+ and [(C14H10)nC58]+) and fullerene (C66 and C68)/anthracene (i.e., [(C14H10)nC66]+ and [(C14H10)nC68]+) cluster cations, are formed in the gas phase through an ion-molecule reaction pathway. With irradiation, all the fullerene/anthracene cluster cations dissociate into mono$-$anthracene and fullerene species without dehydrogenation. The structure of newly formed fullerene/anthracene cluster cations and the bonding energy for these reaction pathways are investigated with quantum chemistry calculations. Our results provide a growth route towards large fullerene derivatives in a bottom-up process and insight in their photo-evolution behavior in the ISM, and clearly, when conditions are favorable, fullerene/PAH clusters can form efficiently. In addition, these clusters (from 80 to 154 atoms or ~ 2 nm in size) offer a good model for understanding the physical-chemical processes involved in the formation and evolution of carbon dust grains in space, and provide candidates of interest for the DIBs that could motivate spectroscopic studies.

The statistical properties of X-ray flares from two separate locations (nucleus and HST-1) in the M87 jet are investigated to reveal the physical origin of the flares. We analyse the archival \textit{Chandra} data for M87, and identify 14 flares in the nucleus and 9 flares in HST-1. The peak intensity ($I_{\rm{P}}$) and the flaring duration time ($T_{\rm{fl}}$) for each flare are obtained. It is found that the distributions of both $I_{\rm{P}}$ and $T_{\rm{fl}}$ for the nucleus obey a power-law form with a similar index. A similar result is also obtained for HST-1, and no significant inconsistency between the nucleus and HST-1 is found for the indices. Similar to solar X-ray flares, the power-law distributions of the flare event parameters can be well explained by a self-organized criticality (SOC) system, which are triggered by magnetic reconnection. Our results suggest that the flares from nucleus and HST-1 are possibly triggered by magnetic reconnection process. The consistent indices for the distributions of $I_{\rm{P}}$ and $T_{\rm{fl}}$ in the CORE and HST-1 indicate that the dimensions of the energy dissipation of the magnetic reconnection is identical in the two regions. A strong correlation between the flares in the two regions also suggests a similar physical origin for the flares.

We performed a survey in the SiO $J=5\rightarrow4$ line toward a sample of 199 Galactic massive star-forming regions at different evolutionary stages with the SMT 10 m and CSO 10.4 m telescopes. The sample consists of 44 infrared dark clouds (IRDCs), 86 protostellar candidates, and 69 young \HII\ regions. We detected SiO $J=5\rightarrow4$ line emission in 102 sources, with a detection rate of 57\%, 37\%, and 65\% for IRDCs, protostellar candidates, and young \HII\ regions, respectively. We find both broad line with Full Widths at Zero Power (FWZP) $>$ 20 \kms and narrow line emissons of SiO in objects at various evolutionary stages, likely associated with high-velocity shocks and low-velocity shocks, respectively. The SiO luminosities do not show apparent differences among various evolutionary stages in our sample. We find no correlation between the SiO abundance and the luminosity-to-mass ratio, indicating that the SiO abundance does not vary significantly in regions at different evolutionary stages of star formation.

The Five Hundred Meter Aperture Spherical Radio Telescope (FAST) will make contributions to the study of Galactic and extragalactic masers. The telescope, now finished construction and commissioning in China, has an innovative design that leads to the highest sensitivity of any single dish telescope in the world. FAST's potential for OH megamaser research is discussed, including the sky density of masers detectable in surveys. The scientific impact expected from FAST maser studies is also discussed.

The hadronic interaction of cosmic rays with solar atmosphere can produce high energy gamma rays. The gamma-ray luminosity is correlated both with the flux of primary cosmic rays and the intensity of the solar magnetic field. The gamma rays below 200 GeV have been observed by $Fermi$ without any evident energy cutoff. The bright gamma-ray flux above 100 GeV has been detected only during solar minimum. The only available data in TeV range come from the HAWC observations, however outside the solar minimum. The ARGO-YBJ dataset has been used to search for sub-TeV/TeV gamma rays from the Sun during the solar minimum from 2008 to 2010, the same time period covered by the Fermi data. A suitable model containing the Sun shadow, solar disk emission and inverse-Compton emission has been developed, and the chi-square minimization method was used to quantitatively estimate the disk gamma-ray signal. The result shows that no significant gamma-ray signal is detected and upper limits to the gamma-ray flux at 0.3$-$7 TeV are set at 95\% confidence level. In the low energy range these limits are consistent with the extrapolation of the Fermi-LAT measurements taken during solar minimum and are compatible with a softening of the gamma-ray spectrum below 1 TeV. They provide also an experimental upper bound to any solar disk emission at TeV energies. Models of dark matter annihilation via long-lived mediators predicting gamma-ray fluxes > $10^{-7}$ GeV $cm^{-2}$ $s^{-1}$ below 1 TeV are ruled out by the ARGO-YBJ limits.

A relativistic electron-positron pair beam can be produced in the interaction of TeV photons from a blazar with the extragalactic background light (EBL). The relativistic $e^{\pm}$ pairs would loss energy through inverse-Compton scattering (ICS) photons of cosmic microwave background (CMB) or plasma instabilities. The dominant energy-loss process is under debate. Based on the assumption that the dominant energy-loss process is ICS, the resulted cascade GeV radiation is usually used to constrain the intergalactic magnetic field (IGMF). Here, we include the energy-loss due to plasma oblique instability in the calculation of cascade gamma-ray flux, and investigate the impact of the plasma instability on the constraint of IGMF. The up-to-date GeV data and archival TeV data of the blazar 1ES 0229+200 are used. The results indicate that even if the oblique instability cooling is dominating over ICS cooling, the cascade flux could be still used to constrain the IGMF. It is found that with the ratio between the cooling rates of the oblique instability and the ICS varying from 0.1, 1 to 10, the lower limit of the IGMF putted by the cascade flux and the gamma-ray data changes from $8\times10^{-18}\ $G, $5\times10^{-18}\ $G to $10^{-18}\ $G. If the ratio between the two cooling rates is 30, the estimate of IGMF based on the cascade flux is invalid.

The mechanism for producing gamma-ray quasi-periodic oscillation (QPO) in blazar is unknown. One possibility is the geometric model, in which without the need for intrinsic quasi-periodic variation, the relativistic Doppler factor changes periodically, resulting in observed gamma-ray QPO. We propose a method to test this geometric model. We analyze the {\it Fermi}-LAT data of PG 1553+113 spanning from 2008 August until 2018 February. According to 29 four-month average spectral energy distributions (SEDs) in the energy range of 0.1-300 GeV, we split the {\it Fermi}-LAT energy range into three bands: 0.1-1 GeV, 1-10 GeV, and 10-300 GeV. The spectrum in each energy range can be successfully fitted by a power-law. The light curves and photon indices in the three energy ranges are obtained. Then, light curves in three narrow energy ranges, i.e., 0.2-0.5 GeV, 2-5 GeV and 20- 40 GeV, are constructed, and the relative variability amplitudes in the three narrow energy ranges are calculated. A discrete-correlation analysis is performed for the light curves. Our results indicate that (i) the light curves in the different energy ranges follow the same pattern showed in the light curve above 0.1 GeV; (ii) the three groups of photon indices in the energy ranges of 0.1-1 GeV, 1-10 GeV, and 10-300 GeV keep nearly constant; (iii) the ratio between relative variability amplitudes in different narrow energy ranges are equal (within their errors) to the prediction by the Doppler effect. Our results support the scenario of the relativistic boost producing the gamma-ray QPO for PG 1553+113.

Tibet is known as the third pole of the earth, as high as the South Pole and North Pole. The Ngari (Ali) observatory in Tibet has the advantage of plenty of photometric night, low precipitable water vapor, high transmittance, good seeing. It is a good site, and promising to be one of the best place for infrared and submillimeter observations in the world. However, there is no data available for sky background brightness in such place, which restrict the astronomical development of the sites. In the near infrared band of J, H, Ks, a NIR sky brightness monitor (NISBM) is designed based on InGaAs photoelectric diode. By using the method of chopper modulation and digital lock-in amplifier processing, the SNR (Signal Noise Ratio), detectivity and the data acquisition speed of the device is greatly improved. For each band of J, H, Ks, an independent instrument is designed and calibrated in laboratory. The NISBM has been installed in Ngari observatory in July of 2017 and obtained the first data of NIR sky brightness at Ngari observatory.

This paper reports on the observation of the sidereal large-scale anisotropy of cosmic rays using data collected by the ARGO-YBJ experiment over 5 years (2008$-$2012). This analysis extends previous work limited to the period from 2008 January to 2009 December,near the minimum of solar activity between cycles 23 and 24.With the new data sample the period of solar cycle 24 from near minimum to maximum is investigated. A new method is used to improve the energy reconstruction, allowing us to cover a much wider energy range, from 4 to 520 TeV. Below 100 TeV, the anisotropy is dominated by two wide regions, the so-called "tail-in" and "loss-cone" features. At higher energies, a dramatic change of the morphology is confirmed. The yearly time dependence of the anisotropy is investigated. Finally, no noticeable variation of cosmic-ray anisotropy with solar activity is observed for a median energy of 7 TeV.

The galaxy Malin 1 contains the largest stellar disk known but the formation mechanism of this structure has been elusive. In this paper, we report a Malin 1 analogue in the 100 Mpc IllustrisTNG simulation and describe its formation history. At redshift zero, this massive galaxy, having a maximum circular velocity $V_{\rm max}$ of 430 ${\rm km\ s^{-1}}$, contains a 100 kpc gas/stellar disk with morphology similar to Malin 1. The simulated galaxy reproduces well many observed features of Malin 1's vast disk, including its stellar ages, metallicities, and gas rotation curve. We trace the extended disk back in time and find that a large fraction of the cold gas at redshift zero originated from the cooling of hot halo gas, triggered by the merger of a pair of intruding galaxies. Our finding provides a novel way to form large galaxy disks as extreme as Malin 1 within the current galaxy formation framework.

A correlation between the secondary cosmic ray flux and the near-earth electric field intensity, measured during thunderstorms, has been found by analyzing the data of the ARGO-YBJ experiment, a full coverage air shower array located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China). The counting rates of showers with different particle multiplicities, have been found to be strongly dependent upon the intensity and polarity of the electric field measured during the course of 15 thunderstorms. In negative electric fields (i.e. accelerating negative charges downwards), the counting rates increase with increasing electric field strength. In positive fields, the rates decrease with field intensity until a certain value of the field EFmin (whose value depends on the event multiplicity), above which the rates begin increasing. By using Monte Carlo simulations, we found that this peculiar behavior can be well described by the presence of an electric field in a layer of thickness of a few hundred meters in the atmosphere above the detector, which accelerates/decelerates the secondary shower particles of opposite charge, modifying the number of particles with energy exceeding the detector threshold. These results, for the first time, give a consistent explanation for the origin of the variation of the electron/positron flux observed for decades by high altitude cosmic ray detectors during thunderstorms.

$N$-body simulations study the dynamics of $N$ particles under the influence of mutual long-distant forces such as gravity. In practice, $N$-body codes will violate Newton's third law if they use either an approximate Poisson solver or individual timesteps. In this study, we construct a novel $N$-body scheme by combining a fast multipole method (FMM) based Poisson solver and a time integrator using a hierarchical Hamiltonian splitting (HHS) technique. We test our implementation for collision-less systems using several problems in galactic dynamics. As a result of the momentum conserving nature of these two key components, the new $N$-body scheme is also momentum conserving. Moreover, we can fully utilize the $\mathcal O(\textit N)$ complexity of FMM with the integrator. With the restored force symmetry, we can improve both angular momentum conservation and energy conservation substantially. The new scheme will be suitable for many applications in galactic dynamics and structure formation. Our implementation, in the code Taichi, is publicly available at https://bitbucket.org/qirong_zhu/taichi_public/.

A spectral line survey of IRC +10216 between 13.3 and 18.5 GHz is carried out using the Shanghai Tian Ma 65 m Radio Telescope (TMRT-65m) with a sensitivity of < 7 mK. Thirty-five spectral lines of 12 different molecules and radicals are detected in total. Except for SiS, the detected molecules are all carbon-chain molecules, including HC3N, HC5N, HC7N, HC9N, C6H, C6H-, C8H, SiC2, SiC4, c-C3H2 and l-C5H. The presence of rich carbon-bearing molecules is consistent with the identity of IRC +10216 as a carbon-rich AGB star. The excitation temperatures and column densities of the observed species are derived by assuming a local thermodynamic equilibrium and homogeneous conditions.

The ground-breaking detections of gravitational waves from black hole mergers by LIGO have rekindled interest in primordial black holes (PBHs) and the possibility of dark matter being composed of PBHs. It has been suggested that PBHs of tens of solar masses could serve as dark matter candidates. Recent analytical studies demonstrated that compact ultra-faint dwarf galaxies can serve as a sensitive test for the PBH dark matter hypothesis, since stars in such a halo-dominated system would be heated by the more massive PBHs, their present-day distribution can provide strong constraints on PBH mass. In this study, we further explore this scenario with more detailed calculations, using a combination of dynamical simulations and Bayesian inference methods. The joint evolution of stars and PBH dark matter is followed with a Fokker-Planck code PhaseFlow. We run a large suite of such simulations for different dark matter parameters, then use a Markov Chain Monte Carlo approach to constrain the PBH properties with observations of ultra-faint galaxies. We find that two-body relaxation between the stars and PBH drives up the stellar core size, and increases the central stellar velocity dispersion. Using the observed half-light radius and velocity dispersion of stars in the compact ultra-faint dwarf galaxies as joint constraints, we infer that these dwarfs may have a cored dark matter halo with the central density in the range of 1-2 $\rm{M_{\odot}/pc^3}$, and that the PBHs may have a mass range of 2-14 $\rm{M_{\odot}}$ if they constitute all or a substantial fraction of the dark matter.

Stellar feedback from high-mass stars (e.g., H{\sc ii} regions) can strongly influence the surrounding interstellar medium and regulate star formation. Our new ALMA observations reveal sequential high-mass star formation taking place within one sub-virial filamentary clump (the G9.62 clump) in the G9.62+0.19 complex. The 12 dense cores (MM 1-12) detected by ALMA are at very different evolutionary stages, from starless core phase to UC H{\sc ii} region phase. Three dense cores (MM6, MM7/G, MM8/F) are associated with outflows. The mass-velocity diagrams of outflows associated with MM7/G and MM8/F can be well fitted with broken power laws. The mass-velocity diagram of SiO outflow associated with MM8/F breaks much earlier than other outflow tracers (e.g., CO, SO, CS, HCN), suggesting that SiO traces newly shocked gas, while the other molecular lines (e.g., CO, SO, CS, HCN) mainly trace the ambient gas continuously entrained by outflow jets. Five cores (MM1, MM3, MM5, MM9, MM10) are massive starless core candidates whose masses are estimated to be larger than 25 M$_{\sun}$, assuming a dust temperature of $\leq$ 20 K. The shocks from the expanding H{\sc ii} regions ("B" \& "C") to the west may have great impact on the G9.62 clump through compressing it into a filament and inducing core collapse successively, leading to sequential star formation. Our findings suggest that stellar feedback from H{\sc ii} regions may enhance the star formation efficiency and suppress the low-mass star formation in adjacent pre-existing massive clumps.

It has been reported that a 5.7sigma directional muon excess coincident with the 2000 July 14 solar flare was registered by the L3 precision muon spectrometer [Ruiguang Wang, Astroparticle Phys., 31(2009) 149]. Using a same analysis method and similar criteria of event selection, we have analyzed the L3 precision muon spectrometer data during November 2000. The result shows that a 4.7sigma muon excess appeared at a time coincident with the solar flare of 8 of November 2000. This muon excess corresponds to above 40 GeV primary protons which came from a sky cell of solid angle 0.048 sr. The probability of being a background fluctuation is estimated to be about 0.1%. It has been convinced that solar protons could be accelerated to tens of GeV in those Class X solar flares which usually arose solar cosmic ray ground level enhancement (GLE) events. However, whether a Class M solar flare like the non-GLE event of 8 November 2000 may also accelerate solar protons to such high energies? It is interesting and noteworthy.

The "Disk of satellites" (DoS) around Milky Way is a highly debated topic with conflicting interpretations of observations and their theoretical models. We perform a comprehensive analysis of all dwarfs detected in the Milky Way and find that the DoS structure depends strongly on the plane identification method and the sample size. In particular, we demonstrate that a smaller sample size produces a higher anisotropy of the spatial distribution and a stronger clustering of the angular momentum of the satellites. Moreover, we calculate the evolution of the 11 classical satellites with proper motion measurements and find that the thin DoS they currently reside in is transient. Furthermore, we analyze two cosmological simulations using the same initial conditions of a Milky Way-sized galaxy, an N-body run with dark matter only and a hydrodynamic one with both baryonic and dark matter, and find that the hydrodynamic simulation produces more anisotropic distribution of satellites than the N-body one. Our results suggest that an anisotropic distribution of satellites in galaxies can originate from baryonic processes in the hierarchical structure formation model, but the claimed highly-flattened, coherently-rotating DoS of the Milky Way may be a small-number selection effect. These findings may help resolve the contradictory claims of DoS in other galaxies and the discrepancy among numerical simulations.

The ARGO-YBJ experiment, a full coverage extensive air shower (EAS) detector located at high altitude (4300 m a.s.l.) in Tibet, China, has smoothly taken data, with very high stability, since November 2007 to the beginning of 2013. The array consisted of a carpet of about 7000 m$^2$ Resistive Plate Chambers (RPCs) operated in streamer mode and equipped with both digital and analog readout, providing the measurement of particle densities up to few particles per cm$^2$. The unique detector features (full coverage, readout granularity, wide dynamic range, etc) and location (very high altitude) allowed a detailed study of the lateral density profile of charged particles at ground very close to the shower axis and its description by a proper lateral distribution function (LDF). In particular, the information collected in the first 10 m from the shower axis have been shown to provide a very effective tool for the determination of the shower development stage ("age") in the energy range 50 TeV - 10 PeV. The sensitivity of the age parameter to the mass composition of primary Cosmic Rays is also discussed.

Observations of globular clusters show that they have universal lognormal mass functions with a characteristic peak at $\sim 2\times 10^{5}\, {\rm{M_{\odot}}}$, but the origin of this peaked distribution is highly debated. Here we investigate the formation and evolution of star clusters in interacting galaxies using high-resolution hydrodynamical simulations performed with two different codes in order to mitigate numerical artifacts. We find that massive star clusters in the range of $\sim 10^{5.5} - 10^{7.5}\, {\rm{M_{\odot}}}$ form preferentially in the highly-shocked regions produced by galaxy interactions. The nascent cluster-forming clouds have high gas pressures in the range of $P/k \sim 10^8 - 10^{12}\, \rm{K}\,\rm{cm^{-3}}$, which is $\sim 10^4 - 10^8$ times higher than the typical pressure of the interstellar medium but consistent with recent observations of a pre-super star cluster cloud in the Antennae Galaxies. Furthermore, these massive star clusters have quasi-lognormal initial mass functions with a peak around $\sim 10^{6}\, {\rm{M_{\odot}}}$. The number of clusters declines with time due to destructive processes, but the shape and the peak of the mass functions do not change significantly during the course of galaxy collisions. Our results suggest that gas-rich galaxy mergers may provide a favorable environment for the formation of massive star clusters such as globular clusters, and that the lognormal mass functions and the unique peak may originate from the extreme high-pressure conditions of the birth clouds and may survive the dynamical evolution.

We implement and test the exact time integration method proposed by Townsend 2009 for gas cooling in cosmological hydrodynamic simulations. The errors using this time integrator for the internal energy are limited by the resolution of the cooling tables and are insensitive to the size of the timestep, improving accuracy relative to explicit or implicit schemes when the cooling time is short. We compare results with different time integrators for gas cooling in cosmological hydrodynamic simulations. We find that the temperature of the gas in filaments before accreting into dark matter halos to form stars, obtained with the exact cooling integration, lies close to the equilibrium where radiative cooling balances heating from the UV background. For comparison, the gas temperature without the exact integrator shows substantial deviations from the equilibrium relation. Galaxy stellar masses with the exact cooling technique agree reasonably well, but are systematically lower than the results obtained by the other integration schemes, reducing the need for feedback to suppress star formation. Our implementation of the exact cooling technique is provided and can be easily incorporated into any hydrodynamic code.

We present the results of cosmological hydrodynamic simulations with zoom-in initial conditions, and investigate the formation of the first galaxies and their evolution towards observable galaxies at $z \sim 6$. We focus on three different galaxies which end up in halos with masses $M_{h} = 2.4 \times10^{10}~h^{-1}\; M_{\odot}$ (Halo-10), $1.6 \times10^{11}~h^{-1}\; M_{\odot}$ (Halo-11) and $0.7 \times10^{12}~h^{-1} M_{\odot}$ (Halo-12) at z=6. Our simulations also probe impacts of different sub-grid assumptions, i.e., SF efficiency and cosmic reionization, on SF histories in the first galaxies. We find that star formation occurs intermittently due to supernova (SN) feedback at z > 10, and then it proceeds more smoothly as the halo mass grows at lower redshifts. Galactic disks are destroyed due to SN feedback, while galaxies in simulations with no-feedback or lower SF efficiency models can sustain galactic disk for long periods > 10 Myr. The expulsion of gas at the galactic center also affects the inner dark matter density profile. However, SN feedback does not seem to keep the shallow profile of dark matter for a long period. Our simulated galaxies in Halo-11 and Halo-12 reproduce the star formation rates (SFR) and stellar masses of observed Lyman-$\alpha$ emitters (LAEs) at z = 7-8 fairly well given observational uncertainties. In addition, we investigate the effect of UV background radiation on star formation as an external feedback source, and find that earlier reionization extends the quenching time of star formation due to photo-ionization heating, but does not affect the stellar mass at z=6.

The ARGO-YBJ detector, located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China), was a full coverage air shower array dedicated to gamma ray astronomy and cosmic ray studies. The wide field of view (~ 2 sr) and high duty cycle (> 86%), made ARGO-YBJ suitable to search for short and unexpected gamma ray emissions like gamma ray bursts (GRBs). Between 2007 November 6 and 2013 February 7, 156 satellite-triggered GRBs (24 of them with known redshift) occurred within the ARGO-YBJ field of view. A search for possible emission associated to these GRBs has been made in the two energy ranges 10-100 GeV and 10-1000 GeV. No significant excess has been found in time coincidence with the satellite detections nor in a time window of one hour after the bursts. Taking into account the EBL absorption, upper limits to the energy fluence at 99% of confidence level have been evaluated,with values ranging from ~ 10-5 erg cm-2 to ~10-1 erg cm-2.

It has been suggested that the satellite galaxies of the Milky Way reside in a highly-flattened, kinematically-coherent plane called Disk of Satellites (DoS). The origin of the DoS, however, has been hotly debated, and a number of conflicting claims have been reported in the literature on whether or not the DoS is consistent with predictions from the standard Lambda Cold Dark Matter ($\Lambda$CDM) cosmological model. Here we investigate this issue by comparing a high-resolution, hydrodynamic $\Lambda$CDM simulation of a Milky Way sized galaxy with its dark matter only counterpart. We find the following results: (1) The abundance and distribution of satellite galaxies around a host galaxy is significantly different in the hydro simulation compared to its N-body counterpart; (2) No clear coherent rotation is found in the satellite system, as the fractions of corotating and counter-corotating satellites remain comparable across cosmic time; (3) The satellite distribution evolves significantly with time, from nearly isotropic at high redshift to anisotropic at the present day; (4) The DoS properties strongly depend on sample selection and plane identification methods. Our results imply that the spatially-thin and coherently-rotating DoS reported in Milky Way and other galaxies may be a selection effect of small sample size.

We study the orbital properties of dark matter haloes by combining a spectral method and cosmological simulations of Milky Way-sized galaxies. We compare the dynamics and orbits of individual dark matter particles from both hydrodynamic and $N$-body simulations, and find that the fraction of box, tube and resonant orbits of the dark matter halo decreases significantly due to the effects of baryons. In particular, the central region of the dark matter halo in the hydrodynamic simulation is dominated by regular, short-axis tube orbits, in contrast to the chaotic, box and thin orbits dominant in the $N$-body run. This leads to a more spherical dark matter halo in the hydrodynamic run compared to a prolate one as commonly seen in the $N$-body simulations. Furthermore, by using a kernel based density estimator, we compare the coarse-grained phase-space densities of dark matter haloes in both simulations and find that it is lower by $\sim0.5$ dex in the hydrodynamic run due to changes in the angular momentum distribution, which indicates that the baryonic process that affects the dark matter is irreversible. Our results imply that baryons play an important role in determining the shape, kinematics and phase-space density of dark matter haloes in galaxies.

The long-standing challenge of creating a Milky Way-like disk galaxy from cosmological simulations has motivated significant developments in both numerical methods and physical models in recent years. We investigate these two fundamental aspects in a new comparison project using a set of cosmological hydrodynamic simulations of the formation and evolution of a Milky Way-size galaxy. In this study, we focus on the comparison of two particle-based hydrodynamics methods: the improved smoothed particle hydrodynamics (SPH) code Gadget, and the Lagrangian Meshless Finite-Mass (MFM) code GIZMO. All the simulations in this paper use the same initial conditions and physical models, which include physics of both dark matter and baryons, star formation, "energy-driven" outflow, metal-dependent cooling, stellar evolution and metal enrichment from supernovae. We find that both numerical schemes produce a late-type galaxy with extended gaseous and stellar disks. However, notable differences are present in a wide range of galaxy properties and their evolution, including star formation history, gas content, disk structure and kinematics. In particular, there is significant difference in gas properties and their evolution between the two simulations. Compared to GIZMO, Gadget simulation produces a larger fraction of cold, dense gas at high redshift which fuels rapid star formation and results in a higher stellar mass by $20\%$ and a lower gas fraction by $10\%$ at $z = 0$, and the resulting gas disk is smoother and more coherent in rotation due to damping of turbulent motion by the numerical viscosity in SPH, in contrast to the GIZMO simulation which shows more prominent spiral structure. Given its better convergence properties and lower computational cost, we argue that MFM method is a promising alternative to the widely used SPH in cosmological hydrodynamic simulations.

We study the impact of baryons on the distribution of dark matter in a Milky Way-size halo by comparing a high-resolution, moving-mesh cosmological simulation with its dark matter-only counterpart. We identify three main processes related to baryons -- adiabatic contraction, tidal disruption and reionization -- which jointly shape the dark matter distribution in both the main halo and its subhalos. The relative effect of each baryonic process depends strongly on the subhalo mass. For massive subhalos with maximum circular velocity $v_{\rm max} > 35 km/s$, adiabatic contraction increases the dark matter concentration, making these halos less susceptible to tidal disruption. For low-mass subhalos with $v_{\rm max} < 20 km/s$, reionization effectively reduces their mass on average by $\approx$ 30% and $v_{\rm max}$ by $\approx$ 20%. For intermediate subhalos with $20 km/s < v_{\rm max} < 35 km/s$, which share a similar mass range as the classical dwarf spheroidals, strong tidal truncation induced by the main galaxy reduces their $v_{\rm max}$. Moreover, the stellar disk of the main galaxy effectively depletes subhalos near the central region. As a combined result of reionization and increased tidal disruption, the total number of low-mass subhalos in the hydrodynamic simulation is nearly halved compared to that of the $\textit{N-}$body simulation. We do not find dark matter cores in dwarf galaxies, unlike previous studies that employed bursty feedback-driven outflows. The substantial impact of baryons on the abundance and internal structure of subhalos suggests that galaxy formation and evolution models based on $\textit{N}$-body simulations should include these physical processes as major components.

We use the Illustris simulation to study the relative contributions of in situ star formation and stellar accretion to the build-up of galaxies over an unprecedentedly wide range of masses ($M_{\ast} = 10^9-10^{12} \, {\rm M_{\odot}}$), galaxy types, environments, and assembly histories. We find that the `two-phase' picture of galaxy formation predicted by some models is a good approximation only for the most massive galaxies in our simulation -- namely, the stellar mass growth of galaxies below a few times $10^{11} \, {\rm M_{\odot}}$ is dominated by in situ star formation at all redshifts. The fraction of the total stellar mass of galaxies at $z=0$ contributed by accreted stars shows a strong dependence on galaxy stellar mass, ranging from about 10 per cent for Milky Way-sized galaxies to over 80 per cent for $M_{\ast} \approx 10^{12} \, {\rm M_{\odot}}$ objects, yet with a large galaxy-to-galaxy variation. At a fixed stellar mass, elliptical galaxies and those formed at the centres of younger haloes exhibit larger fractions of ex situ stars than disc-like galaxies and those formed in older haloes. On average, $\sim$50 per cent of the ex situ stellar mass comes from major mergers (stellar mass ratio $\mu > 1/4$), $\sim$20 per cent from minor mergers ($1/10 < \mu < 1/4$), $\sim$20 per cent from very minor mergers ($\mu < 1/10$), and $\sim$10 per cent from stars that were stripped from surviving galaxies (e.g. flybys or ongoing mergers). These components are spatially segregated, with in situ stars dominating the innermost regions of galaxies, and ex situ stars being deposited at larger galactocentric distances in order of decreasing merger mass ratio.

We present spectral line mapping observations toward four massive star-forming regions (Cepheus A, DR21S, S76E and G34.26+0.15), with the IRAM 30 meter telescope at 2 mm and 3 mm bands. Totally 396 spectral lines from 51 molecules, one helium recombination line, ten hydrogen recombination lines, and 16 unidentified lines were detected in these four sources. An emission line of nitrosyl cyanide (ONCN, 14$_{0,14}$-13$_{0,13}$) was detected in G34.26+0.15, as first detection in massive star-forming regions. We found that the $c$-C$_{3}$H$_{2}$ and NH$_{2}$D show enhancement in shocked regions as suggested by evidences of SiO and/or SO emission. Column density and rotational temperature of CH$_{3}$CN were estimated with the rotational diagram method for all four sources. Isotope abundance ratios of $^{12}$C/$^{13}$C were derived using HC$_{3}$N and its $^{13}$C isotopologue, which were around 40 in all four massive star-forming regions and slightly lower than the local interstellar value ($\sim$65). $^{14}$N/$^{15}$N and $^{16}$O/$^{18}$O abundance ratios in these sources were also derived using double isotopic method, which were slightly lower than that in local interstellar medium. Except for Cep A, $^{33}$S/$^{34}$S ratio in the other three targets were derived, which were similar to that in the local interstellar medium. The column density ratios of N(DCN)/N(HCN) and N(DCO$^{+}$)/N(HCO$^{+}$) in these sources were more than two orders of magnitude higher than the elemental [D]/[H] ratio, which is 1.5$\times$10$ ^{-5}$. Our results show the later stage sources, G34.26+0.15 in particular, present more molecular species than earlier stage ones. Evidence of shock activity is seen in all stages studied.

We test one-zone synchrotron self-Compton (SSC) models with high-quality multiwavelength spectral energy distribution (SED) data of Mrk 421. We use Markov chain Monte Carlo (MCMC) technique to fit twelve day-scale SEDs of Mrk 421 with one-zone SSC models. Three types of electron energy distribution (EED), a log-parabola (LP) EED, a power-law log-parabola (PLLP) EED and a broken power-law (BPL) EED, are assumed in fits. We find that the one-zone SSC model with the PLLP EED provides successful fits to all the twelve SEDs. However, the one-zone SSC model with the LP and BPL EEDs fail to provide acceptable fits to the highest energy X-ray data or GeV data in several states. We therefore conclude that the one-zone SSC model works well in explaining the SEDs of Mrk 421, and the PLLP EED is preferred over the LP and BPL EEDs for Mrk 421 during the flare in March 2010. We derive magnetic field $B'\sim0.01$\ G, Doppler factor $\delta_{\rm D}\sim$30--50, and the curvature parameter of EED $r\sim1$--$10$ in the model with the PLLP EED. The evolutions of model parameters are explored. The physical implications of our results are discussed.

Recently, Advanced Laser Interferometer Gravitational-wave Observatory (aLIGO) detected gravitational wave (GW) transients from mergers of binary black holes (BHs). The system may also produce a wide-angle, relativistic outflow if the claimed short GRB detected by GBM is in real association with GW 150914. It was suggested that mergers of double neutron stars (or neutron star-black hole binaries), another promising source of GW transients, also produce fast, wide-angle outflows. In this paper, we calculate the high-energy gamma-ray emission arising from the blast waves driven by these wide-angle outflows. We find that TeV emission arising from the inverse-Compton process in the relativistic outflow resulted from mergers of binary BHs similar to those in GW 150914 could be detectable by ground-based IACT telescopes such as Cherenkov Telescope Array (CTA) if the sources occur in { a dense medium with density $n > 0.3 cm^{-3}$}. For neutron star-neutron star (NS-NS) and NS-BH mergers, TeV emission from the wide-angle, mildly-relativistic outflow could be detected as well if they occur in a dense medium {with $n > 10-100 cm^{-3}$.} Thus TeV afterglow emission would be a useful probe of the environment of the GW transients, which could shed light on the evolution channels of the progenitors of GW transients.

We report on a measurement of thermal neutrons, generated by the hadronic component of extensive air showers (EAS), by means of a small array of EN-detectors developed for the PRISMA project (PRImary Spectrum Measurement Array), novel devices based on a compound alloy of ZnS(Ag) and $^{6}$LiF. This array has been operated within the ARGO-YBJ experiment at the high altitude Cosmic Ray Observatory in Yangbajing (Tibet, 4300 m a.s.l.). Due to the tight correlation between the air shower hadrons and thermal neutrons, this technique can be envisaged as a simple way to estimate the number of high energy hadrons in EAS. Coincident events generated by primary cosmic rays of energies greater than 100 TeV have been selected and analyzed. The EN-detectors have been used to record simultaneously thermal neutrons and the air shower electromagnetic component. The density distributions of both components and the total number of thermal neutrons have been measured. The correlation of these data with the measurements carried out by ARGO-YBJ confirms the excellent performance of the EN-detector.

ARGO-YBJ, located at the YangBaJing Cosmic Ray Observatory (4300 m a.s.l., Tibet, China), is a full coverage air shower array, with an energy threshold of 300 GeV for gamma-ray astronomy. Most of the recorded events are single front showers, satisfying the trigger requirement of at least 20 particles detected in a given time window. However, in 13% of the events, two randomly arriving showers may be recorded in the same time window, and the second one, in generally smaller, does not need to satisfy the trigger condition. These events are called double front shower events. By using these small showers, well under the trigger threshold, the detector primary energy threshold can be lowered to a few tens of GeV. In this paper, the angular resolution that can be achieved with these events is evaluated by a full Monte Carlo simulation. The ARGO-YBJ sensitivity in detecting gamma-ray bursts (GRBs) by using double front shower events is also studied for various cutoff energies, time durations, and zenith angles of GRBs in the field view of ARGO.

We report on the extensive multi-wavelength observations of the blazar Markarian 421 (Mrk 421) covering radio to gamma-rays, during the 4.5 year period of ARGO-YBJ and Fermi common operation time, from August 2008 to February 2013. In particular, thanks to the ARGO-YBJ and Fermi data, the whole energy range from 100 MeV to 10 TeV is covered without any gap. In the observation period, Mrk 421 showed both low and high activity states at all wavebands. The correlations among flux variations in different wavebands were analyzed. Seven large flares, including five X-ray flares and two GeV gamma-ray flares with variable durations (3-58 days), and one X-ray outburst phase were identified and used to investigate the variation of the spectral energy distribution with respect to a relative quiescent phase. During the outburst phase and the seven flaring episodes, the peak energy in X-rays is observed to increase from sub-keV to few keV. The TeV gamma-ray flux increases up to 0.9-7.2 times the flux of the Crab Nebula. The behavior of GeV gamma-rays is found to vary depending on the flare, a feature that leads us to classify flares into three groups according to the GeV flux variation. Finally, the one-zone synchrotron self-Compton model was adopted to describe the emission spectra. Two out of three groups can be satisfactorily described using injected electrons with a power-law spectral index around 2.2, as expected from relativistic diffuse shock acceleration, whereas the remaining group requires a harder injected spectrum. The underlying physical mechanisms responsible for different groups may be related to the acceleration process or to the environment properties.

We have developed a full numerical method to study the gas dynamics of cometary ultra-compact (UC) H II regions, and associated photodissociation regions (PDRs). The bow-shock and champagne-flow models with a $40.9/21.9 M_\odot$ star are simulated. In the bow-shock models, the massive star is assumed to move through dense ($n=8000~cm^{-3}$) molecular material with a stellar velocity of $15~km~s^{-1}$. In the champagne-flow models, an exponential distribution of density with a scale height of 0.2 pc is assumed. The profiles of the [Ne II] 12.81\mum and $H_2~S(2)$ lines from the ionized regions and PDRs are compared for two sets of models. In champagne-flow models, emission lines from the ionized gas clearly show the effect of acceleration along the direction toward the tail due to the density gradient. The kinematics of the molecular gas inside the dense shell is mainly due to the expansion of the H II region. However, in bow-shock models the ionized gas mainly moves in the same direction as the stellar motion. The kinematics of the molecular gas inside the dense shell simply reflects the motion of the dense shell with respect to the star. These differences can be used to distinguish two sets of models.

The events recorded by ARGO-YBJ in more than five years of data collection have been analyzed to determine the diffuse gamma-ray emission in the Galactic plane at Galactic longitudes 25{\deg} < l < 100{\deg} and Galactic latitudes . The energy range covered by this analysis, from ~350 GeV to ~2 TeV, allows the connection of the region explored by Fermi with the multi-TeV measurements carried out by Milagro. Our analysis has been focused on two selected regions of the Galactic plane, i.e., 40{\deg} < l < 100{\deg} and 65{\deg} < l < 85{\deg} (the Cygnus region), where Milagro observed an excess with respect to the predictions of current models. Great care has been taken in order to mask the most intense gamma-ray sources, including the TeV counterpart of the Cygnus cocoon recently identified by ARGO-YBJ, and to remove residual contributions. The ARGO-YBJ results do not show any excess at sub-TeV energies corresponding to the excess found by Milagro, and are consistent with the predictions of the Fermi model for the diffuse Galactic emission. From the measured energy distribution we derive spectral indices and the differential flux at 1 TeV of the diffuse gamma-ray emission in the sky regions investigated.

The ARGO-YBJ experiment has been in stable data taking from November 2007 till February 2013 at the YangBaJing Cosmic Ray Observatory (4300 m a.s.l.). The detector consists of a single layer of Resistive Plate Chambers (RPCs) ( about 6700 m^2}) operated in streamer mode. The signal pick-up is obtained by means of strips facing one side of the gas volume. The digital readout of the signals, while allows a high space-time resolution in the shower front reconstruction, limits the measurable energy to a few hundred TeV. In order to fully investigate the 1-10 PeV region, an analog readout has been implemented by instrumenting each RPC with two large size electrodes facing the other side of the gas volume. Since December 2009 the RPC charge readout has been in operation on the entire central carpet (about 5800 m^2). In this configuration the detector is able to measure the particle density at the core position where it ranges from tens to many thousands of particles per m^2. Thus ARGO-YBJ provides a highly detailed image of the charge component at the core of air showers. In this paper we describe the analog readout of RPCs in ARGO-YBJ and discuss both the performance of the system and the physical impact on the EAS measurements.

The ARGO-YBJ experiment is a full-coverage air shower detector located at the Yangbajing Cosmic Ray Observatory (Tibet, People's Republic of China, 4300 m a.s.l.). The high altitude, combined with the full-coverage technique, allows the detection of extensive air showers in a wide energy range and offer the possibility of measuring the cosmic ray proton plus helium spectrum down to the TeV region, where direct balloon/space-borne measurements are available. The detector has been in stable data taking in its full configuration from November 2007 to February 2013. In this paper the measurement of the cosmic ray proton plus helium energy spectrum is presented in the region 3-300 TeV by analyzing the full collected data sample. The resulting spectral index is $\gamma = -2.64 \pm 0.01$. These results demonstrate the possibility of performing an accurate measurement of the spectrum of light elements with a ground based air shower detector.

The search for Gamma Ray Burst (GRB) emission in the energy range 1-100 GeV in coincidence with the satellite detection has been carried out using the Astrophysical Radiation with Ground-based Observatory at YangBaJing (ARGO-YBJ) experiment. The high altitude location (4300 m a.s.l.), the large active surface ($\sim$ 6700 m$^2$ of Resistive Plate Chambers), the wide field of view ($\sim 2~$sr, limited only by the atmospheric absorption) and the high duty cycle ($>$ 86 %) make the ARGO-YBJ experiment particularly suitable to detect short and unexpected events like GRBs. With the scaler mode technique, i.e., counting all the particles hitting the detector with no measurement of the primary energy and arrival direction, the minimum threshold of $\sim$ 1 GeV can be reached, overlapping the direct measurements carried out by satellites. During the experiment lifetime, from December 17, 2004 to February 7, 2013, a total of 206 GRBs occurring within the ARGO-YBJ field of view (zenith angle $\theta$ $\le$ 45$^{\circ}$) have been analyzed. This is the largest sample of GRBs investigated with a ground-based detector. Two lightcurve models have been assumed and since in both cases no significant excess has been found, the corresponding fluence upper limits in the 1-100 GeV energy region have been derived, with values as low as 10$^{-5}~$erg cm$^{-2}$. The analysis of a subset of 24 GRBs with known redshift has been used to constrain the fluence extrapolation to the GeV region together with possible cutoffs under different assumptions on the spectrum.

The measurement of cosmic ray energy spectra, in particular for individual species, is an essential approach in finding their origin. Locating the "knees" of the spectra is an important part of the approach and has yet to be achieved. Here we report a measurement of the mixed Hydrogen and Helium spectrum using the combination of the ARGO-YBJ experiment and of a prototype Cherenkov telescope for the LHAASO experiment. A knee feature at 640+/-87 TeV, with a clear steepening of the spectrum, is observed. This gives fundamental inputs to galactic cosmic ray acceleration models.

We simulate evolution of cometary H II regions based on several champagne flow models and bow shock models, and calculate the profiles of the [Ne II] fine-structure line at $12.81\mu m$, the $H30\alpha$ recombination line and the [Ne III] fine-structure line at $15.55\mu m$ for these models at different inclinations of $0^o, 30^o \textrm{and} 60^o$. We find that the profiles in the bow shock models are generally different from those in the champagne flow models, but the profiles in the bow shock with lower stellar velocity ($\leq5km s^{-1}$) are similar to those in the champagne flow models. In champagne flow models, both the velocity of peak flux and the flux weighted central velocities of all three lines are pointing outward from molecular clouds. In bow shock models, the directions of these velocities rely on the speed of stars. They have the similar motion in high stellar speed case but opposite directions in low stellar speed case. We notice that the line profiles from the slit along the symmetrical axis of the projected 2D image of these models are useful for distinguishing bow shock models and champagne flow models. It is also confirmed by the calculation that the flux weighted central velocity and the line luminosity of the [Ne III] line can be estimated from the [Ne II] line and the $H30\alpha$ line.

We studied sulfur chemistry of massive star-forming regions through single-dish submillimeter spectroscopy. OCS, O13CS, 13CS, H2S and SO transitions were observed toward a sample of massive star-forming regions with embedded UCH II or CH II regions. These sources could be divided into H II-hot core and H II-only sources based on their CH3CN emission. Our results show that the OCS line of thirteen sources is optically thick, with optical depth ranging from 5 to 16. Column densities of these molecules were computed under LTE conditions. CS column densities were also derived using its optically thin isotopologue 13CS. H2S is likely to be the most abundant gas-phase sulfuretted molecules in hot assive cores. Both the column density and abundance of sulfur-bearing molecules decrease significantly from H II-hot core to H II-only sources. Ages derived from hot core models appear to be consistent with star-formation theories, suggesting that abundance ratios of [CS]/[SO], [SO]/[OCS] and [OCS]/[CS] could be used as chemical clocks in massive star-forming regions.

The Lyman-alpha (Lya) emission has played an important role in detecting high-redshift galaxies, including recently distant ones at redshift z > 7. It may also contain important information on the origin of these galaxies. Here, we investigate the formation of a typical L* galaxy and its observational signatures at the earliest stage, by combining a cosmological hydrodynamic simulation with three-dimensional radiative transfer calculations using the newly improved ART^2 code. Our cosmological simulation uses the Aquila initial condition which zooms in onto a Milky Way-like halo with high resolutions, and our radiative transfer couples multi-wavelength continuum, Lya line, and ionization of hydrogen. We find that the modeled galaxy starts to form at redshift z ~ 24 through efficient accretion of cold gas, which produces a strong Lya line with a luminosity of L(Lya) ~ 10^42 erg/s as early as z ~ 14. The Lya emission appears to trace the cold, dense gas. The lines exhibit asymmetric, single-peak profiles, and are shifted to the blue wing, a characteristic feature of gas inflow. Moreover, the contribution to the total Lya luminosity by excitation cooling increases with redshift, and it becomes dominant at z >~ 6. We predict that L* galaxies such as the modeled one may be detected at z <~ 8 by JWST and ALMA with a reasonable integration time. Beyond redshift 12, however, only Lya line may be observable by spectroscopic surveys. Our results suggest that Lya line is one of the most powerful tools to detect the first generation of galaxies, and to decipher their formation mechanism.

We study the convergence properties of smoothed particle hydrodynamics (SPH) using numerical tests and simple analytic considerations. Our analysis shows that formal numerical convergence is possible in SPH only in the joint limit $N \rightarrow \infty$, $h \rightarrow 0$, and $N_{nb} \rightarrow \infty$, where $N$ is the total number of particles, $h$ is the smoothing length, and $N_{nb}$ is the number of neighbor particles within the smoothing volume used to compute smoothed estimates. Previous work has generally assumed that the conditions $N \rightarrow \infty$ and $h \rightarrow 0$ are sufficient to achieve convergence, while holding $N_{nb}$ fixed. We demonstrate that if $N_{nb}$ is held fixed as the resolution is increased, there will be a residual source of error that does not vanish as $N \rightarrow \infty$ and $h \rightarrow 0$. Formal numerical convergence in SPH is possible only if $N_{nb}$ is increased systematically as the resolution is improved. Using analytic arguments, we derive an optimal compromise scaling for $N_{nb}$ by requiring that this source of error balance that present in the smoothing procedure. For typical choices of the smoothing kernel, we find $N_{nb} \propto N^{1/2}$. This means that if SPH is to be used as a numerically convergent method, the required computational cost does not scale with particle number as $O(N)$, but rather as $O(N^{1+\delta})$, where $\delta \approx 1/2$, with a weak dependence on the form of the smoothing kernel.

Observations have shown that supermassive black holes in nearby elliptical galaxies correlate tightly with the stellar velocity dispersion (the $\MBH - \sigma$ relation) and the stellar mass (the $\MBH - \Mhost$ relation) of their host spheroids. However, the origin of these correlations remains ambiguous. In a previous paper by Zhu et al., we proposed a model which links the M-$\sigma$ relation to the the dynamical state of the system and the $\MBH - \Mhost$ relation to the self-regulation of galaxy growth. To test this model, we compile a sample of observed galaxies with different properties and examine the dependence of the above correlations on these parameters. We find that galaxies that satisfy the the $\MBH - \sigma$ correlation appear to have reached virial equilibrium, as indicated by the ratio between kinetic energy and gravitational potential, 2K/U $\sim$ 1. Furthermore, the ratio of black hole accretion rate to star formation rate remains nearly constant, BHAR /SFR $\sim$ $10^{-3}$, in active galaxies over a wide range of mass in the redshift range z=0 - 3. These results confirm our theoretical model that the observed correlations have different origins: the $\MBH - \sigma$ relation may result from galaxy relaxation, while the $\MBH$ - $\Mhost$ relation may be due to self-regulated black hole accretion and star formation in galaxies.

The extended TeV gamma-ray source ARGO J2031+4157 (or MGRO J2031+41) is positionally consistent with the Cygnus Cocoon discovered by $Fermi$-LAT at GeV energies in the Cygnus superbubble. Reanalyzing the ARGO-YBJ data collected from November 2007 to January 2013, the angular extension and energy spectrum of ARGO J2031+4157 are evaluated. After subtracting the contribution of the overlapping TeV sources, the ARGO-YBJ excess map is fitted with a two-dimensional Gaussian function in a square region of $10^{\circ}\times 10^{\circ}$, finding a source extension $\sigma_{ext}$= 1$^{\circ}$.8$\pm$0$^{\circ}$.5. The observed differential energy spectrum is $dN/dE =(2.5\pm0.4) \times 10^{-11}(E/1 TeV)^{-2.6\pm0.3}$ photons cm$^{-2}$ s$^{-1}$ TeV$^{-1}$, in the energy range 0.2-10 TeV. The angular extension is consistent with that of the Cygnus Cocoon as measured by $Fermi$-LAT, and the spectrum also shows a good connection with the one measured in the 1-100 GeV energy range. These features suggest to identify ARGO J2031+4157 as the counterpart of the Cygnus Cocoon at TeV energies. The Cygnus Cocoon, located in the star-forming region of Cygnus X, is interpreted as a cocoon of freshly accelerated cosmic rays related to the Cygnus superbubble. The spectral similarity with Supernova Remnants indicates that the particle acceleration inside a superbubble is similar to that in a SNR. The spectral measurements from 1 GeV to 10 TeV allows for the first time to determine the possible spectrum slope of the underlying particle distribution. A hadronic model is adopted to explain the spectral energy distribution.

Since asymptotic giant branch (AGB) stars are bright and extended infrared objects, most Galactic AGB stars saturate the Wide-field Infrared Survey Explorer (WISE) detectors and therefore the WISE magnitudes that are restored by applying point-spread-function fitting need to be verified. Statistical properties of circumstellar envelopes around AGB stars are discussed on the basis of a WISE AGB catalog verified in this way. We cross-matched an AGB star sample with the WISE All-Sky Source Catalog and the Two Mircon All Sky Survey catalog. Infrared Space Observatory (ISO) spectra of a subsample of WISE AGB stars were also exploited. The dust radiation transfer code DUSTY was used to help predict the magnitudes in the W1 and W2 bands, the two WISE bands most affected by saturation, for calibration purpose, and to provide physical parameters of the AGB sample stars for analysis. DUSTY is verified against the ISO spectra to be a good tool to reproduce the spectral energy distributions of these AGB stars. Systematic magnitude-dependent offsets have been identified in WISE W1 and W2 magnitudes of the saturated AGB stars, and empirical calibration formulas are obtained for them on the basis of 1877 (W1) and 1558 (W2) AGB stars that are successfully fit with DUSTY. According to the calibration formulae, the corrections for W1 at 5 mag and W2 at 4 mag are $-0.383$ and 0.217 mag, respectively. In total, we calibrated the W1/W2 magnitudes of 2390/2021 AGB stars. The model parameters from the DUSTY and the calibrated WISE W1 and W2 magnitudes are used to discuss the behavior of the WISE color-color diagrams of AGB stars. The model parameters also reveal that O-rich AGB stars with opaque circumstellar envelopes are much rarer than opaque C-rich AGB stars toward the anti-Galactic center direction, which we attribute to the metallicity gradient of our Galaxy.

A large number of high-redshift galaxies have been discovered via their narrow-band Lya line or broad-band continuum colors in recent years. The nature of the escaping process of photons from these early galaxies is crucial to understanding galaxy evolution and the cosmic reionization. Here, we investigate the escape of Lya, non-ionizing UV-continuum (l = 1300 - 1600 angstrom in rest frame), and ionizing photons (l < 912 angstrom) from galaxies by combining a cosmological hydrodynamic simulation with three-dimensional multi-wavelength radiative transfer calculations. The galaxies are simulated in a box of 5^3 h^-3 Mpc^3 with high resolutions using the Aquila initial condition which reproduces a Milky Way-like galaxy at redshift z=0. We find that the escape fraction (fesc) of these different photons shows a complex dependence on redshift and galaxy properties: fesc(Lya) and fesc(UV) appear to evolve with redshift, and they show similar, weak correlations with galaxy properties such as mass, star formation, metallicity, and dust content, while fesc(Ion) remains roughly constant at ~ 0.2 from z ~ 0 - 10, and it does not show clear dependence on galaxy properties. fesc(Lya) correlates more strongly with fesc(UV) than with fesc(Ion). In addition, we find a relation between the emergent Lya luminosity and the ionizing photon emissivity of Lyman Alpha Emitters (LAEs). By combining this relation with the observed luminosity functions of LAEs at different redshift, we estimate the contribution from LAEs to the reionization of intergalactic medium (IGM). Our result suggests that ionizing photons from LAEs alone are not sufficient to ionize IGM at z > 6, but they can maintain the ionization of IGM at z ~ 0 - 5.

The energy spectrum of cosmic Hydrogen and Helium nuclei has been measured, below the so-called "knee", by using a hybrid experiment with a wide field-of-view Cherenkov telescope and the Resistive Plate Chamber (RPC) array of the ARGO-YBJ experiment at 4300 m above sea level. The Hydrogen and Helium nuclei have been well separated from other cosmic ray components by using a multi-parameter technique. A highly uniform energy resolution of about 25% is achieved throughout the whole energy range (100 TeV - 700 TeV). The observed energy spectrum is compatible with a single power law with index gamma=-2.63+/-0.06.

The ARGO-YBJ detector is an extensive air shower array that has been used to monitor the northern $\gamma$-ray sky at energies above 0.3 TeV from 2007 November to 2013 January. In this paper, we present the results of a sky survey in the declination band from $-10^{\circ}$ to $70^{\circ}$, using data recorded over the past five years. With an integrated sensitivity ranging from 0.24 to $\sim$1 Crab units depending on the declination, six sources have been detected with a statistical significance greater than 5 standard deviations. Several excesses are also reported as potential $\gamma$-ray emitters. The features of each source are presented and discussed. Additionally, $95\%$ confidence level upper limits of the flux from the investigated sky region are shown. Specific upper limits for 663 GeV $\gamma$-ray AGNs inside the ARGO-YBJ field of view are reported. The effect of the absorption of $\gamma$-rays due to the interaction with extragalactic background light is estimated.

Recent observations have discovered a population of extended Lya sources, dubbed Lyman-alpha blobs (LABs), at high redshift z ~ 3 - 6.6. These LABs typically have a luminosity of L ~ 10^42-10^44 erg/s, and a size of tens of kiloparsecs, with some giant ones reaching up to D ~ 100 kpc. However, the origin of these LABs is not well understood. In this paper, we investigate a merger model for the formation of LABs by studying Lya emission from interacting galaxies at high redshifts by means of a combination of hydrodynamics simulations with three dimensional radiative transfer calculations. Our galaxy simulations focus on a set of binary major mergers of galaxies with a mass range of 3-7 *10^12 Msun in the redshift range of z ~ 3 -7, and we use the newly improved ART^2 code to perform the radiative transfer calculations which couple multi-wavelength continuum, ionization of hydrogen, and Lya line emission. We find that intense star formation and enhanced cooling induced by gravitational interaction produce strong Lya emission from these merging galaxies. The Lya emission appears to be extended due to the extended distribution of sources and gas. During the close encounter of galaxy progenitors when the star formation rate peaks at ~ 10^3 Msun/yr, our model produces Lya blobs with luminosity of L ~ 10^42-10^44 erg/s, and size of D ~ 10-20 kpc at z>6 and D ~ 20-50 kpc at z ~ 3, in broad agreement with observations in the same redshift range. Our results suggest that merging galaxies may produce some typical LABs as observed, but the giant ones may be produced by mergers more massive than those in our model, or a combination of mergers and cold accretion from filaments on a large scale.

We report the observation of a very high energy \gamma-ray source, whose position is coincident with HESS J1841-055. This source has been observed for 4.5 years by the ARGO-YBJ experiment from November 2007 to July 2012. Its emission is detected with a statistical significance of 5.3 standard deviations. Parameterizing the source shape with a two-dimensional Gaussian function we estimate an extension \sigma=(0.40(+0.32,-0.22}) degree, consistent with the HESS measurement. The observed energy spectrum is dN/dE =(9.0-+1.6) x 10^{-13}(E/5 TeV)^{-2.32-+0.23} photons cm^{-2} s^{-1} TeV^{-1}, in the energy range 0.9-50 TeV. The integral \gamma-ray flux above 1 TeV is 1.3-+0.4 Crab units, which is 3.2-+1.0 times the flux derived by HESS. The differences in the flux determination between HESS and ARGO-YBJ, and possible counterparts at other wavelengths are discussed.

We present the first multi-wavelength follow-up observations of two candidate gravitational-wave (GW) transient events recorded by LIGO and Virgo in their 2009-2010 science run. The events were selected with low latency by the network of GW detectors and their candidate sky locations were observed by the Swift observatory. Image transient detection was used to analyze the collected electromagnetic data, which were found to be consistent with background. Off-line analysis of the GW data alone has also established that the selected GW events show no evidence of an astrophysical origin; one of them is consistent with background and the other one was a test, part of a "blind injection challenge". With this work we demonstrate the feasibility of rapid follow-ups of GW transients and establish the sensitivity improvement joint electromagnetic and GW observations could bring. This is a first step toward an electromagnetic follow-up program in the regime of routine detections with the advanced GW instruments expected within this decade. In that regime multi-wavelength observations will play a significant role in completing the astrophysical identification of GW sources. We present the methods and results from this first combined analysis and discuss its implications in terms of sensitivity for the present and future instruments.

According to the second law of thermodynamics, the arrow of time points to an ever increasing entropy of the Universe. However, exactly how the entropy evolves with time and what drives the growth remain largely unknown. Here, for the first time, we quantify the evolving entropy of cosmic structures using a large-scale cosmological hydrodynamical simulation. Our simulation starts from initial conditions predicted by the leading LambdaCDM cosmology, self-consistently evolves the dynamics of both dark and baryonic matter, star formation, black hole growth and feedback processes, from the cosmic dawn to the present day. Tracing the entropy contributions of these distinct components in the simulation, we find a strong link between entropy growth and structure formation. The entropy is dominated by that of the black holes in all epochs, and its evolution follows the same path as that of galaxies: it increases rapidly from a low-entropy state at high redshift until z~2, then transits to a slower growth. Our results suggest that cosmic entropy may co-evolve with cosmic structure, and that its growth may be driven mainly by the formation of black holes in galaxies. We predict that the entropy will continue to increase in the near future, but likely at a constant rate.

It is well established that supermassive black holes in nearby elliptical galaxies correlate tightly with the kinematic property ($\mbhsigma$ correlation) and stellar mass ($\mbhhost$ correlation) of their host spheroids. However, it is not clear what the relations would be at the low-mass end, and how they evolve. Here, we investigate these relations in low-mass systems ($\MBH \sim \rm{10^{6}- 10^{8}}\, \Msun$) using the Aquila Simulation, a high-resolution cosmological hydrodynamic simulation which follows the formation and evolution of stars and black holes in a Milky Way-size galaxy and its substructures. We find a number of interesting results on the origin and evolution of the scaling relations in these systems: (1) there is a strong redshift evolution in the $\mbhsigma$ relation, but a much weaker one in the $\mbhhost$ relation; (2) there is a close link between the $\mbhsigma$ relation and the dynamical state of the system -- the galaxies that fall on the observed correlation appear to have reached virial equilibrium. (3) the star formation and black hole growth are self-regulated in galaxies -- the ratio between black hole accretion rate and star formation rate remains nearly constant in a wide redshift span $z = 0-6$. These findings suggest that the observed correlations have different origins: the $\mbhsigma$ relation may be the result of virial equilibrium, while the $\mbhhost$ relation may the result of self-regulated star formation and black hole growth in galaxies.

As one of the brightest active blazars in both X-ray and very high energy $\gamma$-ray bands, Mrk 501 is very useful for physics associated with jets from AGNs. The ARGO-YBJ experiment is monitoring it for $\gamma$-rays above 0.3 TeV since November 2007. Starting from October 2011 the largest flare since 2005 is observed, which lasts to about April 2012. In this paper, a detailed analysis is reported. During the brightest $\gamma$-rays flaring episodes from October 17 to November 22, 2011, an excess of the event rate over 6 $\sigma$ is detected by ARGO-YBJ in the direction of Mrk 501, corresponding to an increase of the $\gamma$-ray flux above 1 TeV by a factor of 6.6$\pm$2.2 from its steady emission. In particular, the $\gamma$-ray flux above 8 TeV is detected with a significance better than 4 $\sigma$. Based on time-dependent synchrotron self-Compton (SSC) processes, the broad-band energy spectrum is interpreted as the emission from an electron energy distribution parameterized with a single power-law function with an exponential cutoff at its high energy end. The average spectral energy distribution for the steady emission is well described by this simple one-zone SSC model. However, the detection of $\gamma$-rays above 8 TeV during the flare challenges this model due to the hardness of the spectra. Correlations between X-rays and $\gamma$-rays are also investigated.

Narrow-band Lya line and broad-band continuum have played important roles in the discovery of high-redshift galaxies in recent years. Hence, it is crucial to study the radiative transfer of both Lya and continuum photons in the context of galaxy formation and evolution in order to understand the nature of distant galaxies. Here, we present a three-dimensional Monte Carlo radiative transfer code, All-wavelength Radiative Transfer with Adaptive Refinement Tree (ART^2), which couples Lya line and multi-wavelength continuum, for the study of panchromatic properties of galaxies and interstellar medium. This code is based on the original version of Li et al., and features three essential modules: continuum emission from X-ray to radio, Lya emission from both recombination and collisional excitation, and ionization of neutral hydrogen. The coupling of these three modules, together with an adaptive refinement grid, enables a self-consistent and accurate calculation of the Lya properties. As an example, we apply ART^2 to a cosmological simulation that includes both star formation and black hole growth, and study in detail a sample of massive galaxies at redshifts z=3.1 - 10.2. We find that these galaxies are Lya emitters (LAEs), whose Lya emission traces the dense gas region, and that their Lya lines show a shape characteristic of gas inflow. Furthermore, the Lya properties, including photon escape fraction, emergent luminosity, and equivalent width, change with time and environment. Our results suggest that LAEs evolve with redshift, and that early LAEs such as the most distant one detected at z ~ 8.6 may be dwarf galaxies with a high star formation rate fueled by infall of cold gas, and a low Lya escape fraction.

The Lya emission has been observed from galaxies over a redshift span z ~ 0 - 8.6. However, the evolution of high-redshift Lya emitters (LAEs), and the link between these populations and local galaxies, remain poorly understood. Here, we investigate the Lya properties of progenitors of a local L* galaxy by combining cosmological hydrodynamic simulations with three-dimensional radiative transfer calculations using the new ART^2 code. We find that the main progenitor (the most massive one) of a Milky Way-like galaxy has a number of Lya properties close to those of observed LAEs at z ~ 2 - 6, but most of the fainter ones appear to fall below the detection limits of current surveys. The Lya photon escape fraction depends sensitively on a number of physical properties of the galaxy, such as mass, star formation rate, and metallicity, as well as galaxy morphology and orientation. Moreover, we find that high-redshift LAEs show blue-shifted Lya line profiles characteristic of gas inflow, and that the Lya emission by excitation cooling increases with redshift, and becomes dominant at z > 6. Our results suggest that some observed LAEs at z ~ 2-6 with luminosity of L_Lya ~ 10^{42-43} ergs/s may be similar to the main progenitor of the Milky Way at high redshift, and that they may evolve into present-day L* galaxies.

We report the observation of TeV gamma-rays from the Cygnus region using the ARGO-YBJ data collected from 2007 November to 2011 August. Several TeV sources are located in this region including the two bright extended MGRO J2019+37 and MGRO J2031+41. According to the Milagro data set, at 20 TeV MGRO J2019+37 is the most significant source apart from the Crab Nebula. No signal from MGRO J2019+37 is detected by the ARGO-YBJ experiment, and the derived flux upper limits at 90% confidence level for all the events above 600 GeV with medium energy of 3 TeV are lower than the Milagro flux, implying that the source might be variable and hard to be identified as a pulsar wind nebula. The only statistically significant (6.4 standard deviations) gamma-ray signal is found from MGRO J2031+41, with a flux consistent with the measurement by Milagro.

ARGO-YBJ is an air shower detector array with a fully covered layer of resistive plate chambers. It is operated with a high duty cycle and a large field of view. It continuously monitors the northern sky at energies above 0.3 TeV. In this paper, we report a long-term monitoring of Mrk 421 over the period from 2007 November to 2010 February. This source was observed by the satellite-borne experiments Rossi X-ray Timing Explorer and Swift in the X-ray band. Mrk 421 was especially active in the first half of 2008. Many flares are observed in both X-ray and gamma-ray bands simultaneously. The gamma-ray flux observed by ARGO-YBJ has a clear correlation with the X-ray flux. No lag between the X-ray and gamma-ray photons longer than 1 day is found. The evolution of the spectral energy distribution is investigated by measuring spectral indices at four different flux levels. Hardening of the spectra is observed in both X-ray and gamma-ray bands. The gamma-ray flux increases quadratically with the simultaneously measured X-ray flux. All these observational results strongly favor the synchrotron self-Compton process as the underlying radiative mechanism.

The sun blocks cosmic ray particles from outside the solar system, forming a detectable shadow in the sky map of cosmic rays detected by the ARGO-YBJ experiment in Tibet. Because the cosmic ray particles are positive charged, the magnetic field between the sun and the earth deflects them from straight trajectories and results in a shift of the shadow from the true location of the sun. Here we show that the shift measures the intensity of the field which is transported by the solar wind from the sun to the earth.

We present a multi-wavelength photometric and spectroscopic study of a newly discovered candidate cluster [DBS2003] 45. Our H, Ks photometry confirms that [DBS2003] 45 is a cluster. An average visual extinction Av 7.1+/-0.5 is needed to fit the cluster sequence with a model isochrone. Low resolution spectroscopy indicates that half a dozen early B and at least one late O type giant stars are present in the cluster. We estimate the age of the cluster to be between 5 and 8 Myr based on spectroscopic analysis. Assuming an age of 6 Myr, we fit the observed mass function with a power law, N(M) M^(-Gamma), and find an index Gamma 1.27+/-0.15, which is consistent with the Salpeter value. We estimate the total cluster mass is around 1000 solar masses by integrating the derived mass function between 0.5 and 45 solar masses. Both mid-infrared and radio wavelength observations show that a bubble filled with ionized gas is associated with the cluster. The total ionizing photon flux estimated from radio continuum measurements is consistent with the number of hot stars we detected. Infrared bright point sources along the rim of the bubble suggest that there is triggered star formation at the periphery of the HII region.

We present high spatial and spectral resolution observations of sixteen Galactic compact and ultracompact H II regions in the [Ne II] 12.8 microns fine structure line. The small thermal width of the neon line and the high dynamic range of the maps provide an unprecedented view of the kinematics of compact and ultracompact H II regions. These observations solidify an emerging picture of the structure of ultracompact H II regions suggested in our earlier studies of G29.96-0.02 and Mon R2 IRS1; systematic surface flows, rather than turbulence or bulk expansion, dominate the gas motions in the H II regions. The observations show that almost all of the sources have significant (5-20 km/s) velocity gradients and that most of the sources are limb-brightened. In many cases, the velocity pattern implies tangential flow along a dense shell of ionized gas. None of the observed sources clearly fits into the categories of filled expanding spheres, expanding shells, filled blister flows, or cometary H II regions formed by rapidly moving stars. Instead, the kinematics and morphologies of most of the sources lead to a picture of H II regions confined to the edges of cavities created by stellar wind ram pressure and flowing along the cavity surfaces. In sources where the radio continuum and [Ne II] morphologies agree, the majority of the ionic emission is blue-shifted relative to nearby molecular gas. This is consistent with sources lying on the near side of their natal clouds being less affected by extinction and with gas motions being predominantly outward, as is expected for pressure-driven flows.

We present results from K band slit scan observations of a ~20''x20'' region of the Galactic center (GC) in two separate epochs more than five years apart. The high resolution (R>=14,000) observations allow the most accurate radial velocity and acceleration measurements of the stars in the central parsec of the Galaxy. Detected stars can be divided into three groups based on the CO absorption band heads at ~2.2935 microns and the He I lines at ~2.0581 microns and ~2.112, 2.113 microns: cool, narrow-line hot and broad-line hot. The radial velocities of the cool, late-type stars have approximately a symmetrical distribution with its center at ~-7.8(+/-10.3) km/s and a standard deviation ~113.7(+/-10.3) km/s. Although our statistics are dominated by the brightest stars, we estimate a central black hole mass of 3.9(+/-1.1) million solar masses, consistent with current estimates from complete orbits of individual stars. Our surface density profile and the velocity dispersion of the late type stars support the existence of a low density region at the Galactic center suggested by earlier observations. Many hot, early-type stars show radial velocity changes higher than maximum values allowed by pure circular orbital motions around a central massive object, suggesting that the motions of these stars greatly deviate from circular orbital motions around the Galactic center. The correlation between the radial velocities of the early type He I stars and their declination offsets from Sagittarius A* suggests that a systematic rotation is present for the early-type population. No figure rotation around the Galactic center for the late type stars is supported by the new observations.

We present arcsecond-scale mid-ir photometry (in the 10.5 micron N band and at 24.8 microns), and low resolution spectra in the N band (R~100) of a candidate high mass protostellar object (HMPO) in IRAS 18151-1208 and of two HMPO candidates in IRAS 20343+4129, IRS 1 and IRS 3. In addition we present high resolution mid-ir spectra (R~80000) of the two HMPO candidates in IRAS 20343+4129. These data are fitted with simple models to estimate the masses of gas and dust associated with the mid-ir emitting clumps, the column densities of overlying absorbing dust and gas, the luminosities of the HMPO candidates, and the likely spectral type of the HMPO candidate for which [Ne II] 12.8 micron emission was detected (IRAS 20343+4129 IRS 3). We suggest that IRAS 18151-1208 is a pre-ultracompact HII region HMPO, IRAS 20343+4129 IRS 1 is an embedded young stellar object with the luminosity of a B3 star, and IRAS 20343+4129 IRS 3 is a B2 ZAMS star that has formed an ultracompact HII region and disrupted its natal envelope.

We have mapped [Ne II] (12.8um) and [S IV] (10.5um) emission from W51 IRS 2 with TEXES on Gemini North, and we compare these data to VLA free-free observations and VLT near-infrared images. With 0.5" spatial and 4 km/s spectral resolution we are able to separate the ionized gas into several components: an extended H II region on the front surface of the molecular cloud, several embedded compact H II regions, and a streamer of high velocity gas. We interpret the high velocity streamer as a precessing or fan-like jet, which has emerged from the molecular cloud into an OB star cluster where it is being ionized.

A new type of cosmic-ray experiment is under construction in the Pyh\"asalmi mine in the underground laboratory of the University of Oulu, Finland. It aims to study the composition of cosmic rays at and above the knee region. The experiment, called EMMA, will cover approximately 150 square-metres of detector area. The array is capable of measuring the multiplicity and the lateral distribution of underground muons, and the arrival direction of the air shower. The full-size detector is expected to run by the end of 2007.

High spectral and spatial resolution, mid-infrared fine structure line observations toward two ultracompact HII (UCHII) regions (G29.96 -0.02 and Mon R2) allow us to study the structure and kinematics of cometary UCHII regions. In our earlier study of Mon R2, we showed that highly organized mass motions accounted for most of the velocity structure in that UCHII region. In this work, we show that the kinematics in both Mon R2 and G29.96 are consistent with motion along an approximately paraboloidal shell. We model the velocity structure seen in our mapping data and test the stellar wind bow shock model for such paraboloidal like flows. The observations and the simulation indicate that the ram pressures of the stellar wind and ambient interstellar medium cause the accumulated mass in the bow shock to flow along the surface of the shock. A relaxation code reproduces the mass flow's velocity structure as derived by the analytical solution. It further predicts that the pressure gradient along the flow can accelerate ionized gas to a speed higher than that of the moving star. In the original bow shock model, the star speed relative to the ambient medium was considered to be the exit speed of ionized gas in the shell.

Ultracompact HII regions are signposts of massive star formation and their properties provide diagnostics for the characteristics of very young O stars embedded in molecular clouds. While radio observations have given us a good picture of the morphology of these regions, they have not provided clear information about the kinematics. Using high spectral resolution observations of the 12.8 micron [NeII] line, it has been possible for the first time to trace the internal kinematics of several ultracompact HII regions. We find that the motions in the cometary ultracompact HII regions MonR2 and G29.96-0.02 are highly organized. The velocity patterns are consistent with parabolic ionized flows along a neutral boundary layer.

We present the first results of a study of the kinematics and morphology of ultracompact HII regions using a new observational technique. We used very high spatial (1.5") and spectral (3.4 km/s)resolution observations of the [NeII] line at 12.8 microns to study the ionized gas in Monoceros R2. The [NeII] emission shows an HII region with highest emission measure in a ~24" shell. Line widths are as narrow as 8 km/s at some positions. In places where the lines are complex and broader, the additional width is most likely due to overlap of narrower features along the line of sight. The global kinematics suggest that the 24" shell is expanding at ~10 km/s. However, the spectral profiles toward the southeast side of the shell are not consistent with a simple expansion picture.

The temporal and spectral properties of terrestrial gamma-ray flashes (TGFs) are studied. The delay of low energy photons relative to high energy ones in the gamma-ray variations of the TGFs with high signal to noise ratio has been revealed by an approach of correlation analysis in the time domain on different time scales. The temporal structures of the TGFs in high energy band are usually narrower than that in low energy band. The spectral hardness has a general trend of decreasing with time in a flash. The observed temporal and spectral characteristics give constraints and valuable hints on the flash production mechanism.

We discuss the design and performance of TEXES, the Texas Echelon Cross Echelle Spectrograph. TEXES is a mid-infrared (5-25 um) spectrograph with several operating modes: high-resolution, cross-dispersed with a resolving power of R=100,000, 0.5% spectral coverage, and a ~1.5 by 8" slit; medium-resolution long-slit with R~15,000, 0.5% coverage, and a ~1.5 by 45" slit; low-resolution long-slit with (delta lambda)~0.004 um, 0.25 um coverage, and a ~1.5 by 45" slit; and source acquisition imaging with 0.33" pixels and a 25 by 25" field of view on a 3m telescope. TEXES has been used at the McDonald Observatory 2.7m and the IRTF 3m telescopes, and has proven to be both sensitive and versatile.

We broadly discuss mid-infrared spectroscopy and detail our new high spectral resolution instrument, the Texas Echelon-cross-Echelle Spectrograph (TEXES).