The Large High Altitude Air Shower Observatory (LHAASO) project is a new generation multi-component instrument, to be built at 4410 meters of altitude in the Sichuan province of China, with the aim to study with unprecedented sensitivity the spec trum, the composition and the anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$ eV, as well as to act simultaneously as a wide aperture (one stereoradiant), continuously-operated gamma ray telescope in the energy range between 10$^{11}$ and $10^{15}$ eV. The experiment will be able of continuously surveying the TeV sky for steady and transient sources from 100 GeV to 1 PeV, t hus opening for the first time the 100-1000 TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov/fluorescence components) that will be measured in LHAASO will allow to investigate origin, acceleration and propagation of the radiation through a measurement of energy spec trum, elemental composition and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO in cosmic rays physics and gamma astronomy would play a key-role in the comprehensive general program to explore the High Energy Universe. LHAASO will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. In this document we introduce the concept of LHAASO and the main science goals, providing an overview of the project.

We present a catalogue of galaxy photometric redshifts for the Sloan Digital Sky Survey (SDSS) Data Release 12. We use various supervised learning algorithms to calculate redshifts using photometric attributes on a spectroscopic training set. Two training sets are analysed in this paper. The first training set consists of 995,498 galaxies with redshifts up to $z \approx 0.8$. On the first training set, we achieve a cost function of 0.00501 and a root mean squared error value of 0.0707 using the XGBoost algorithm. We achieved an outlier rate of 2.1\% and 86.81\%, 95.83\%, 97.90\% of our data points lie within one, two, and three standard deviation of the mean respectively. The second training set consists of 163,140 galaxies with redshifts up to $z\approx0.2$ and is merged with the Galaxy Zoo 2 full catalog. We also experimented on convolutional neural networks to predict five morphological features (Smooth, Features/Disk, Star, Edge-on, Spiral). We achieve a root mean squared error of 0.117 when validated against an unseen dataset with over 200 epochs. Morphological features from the Galaxy Zoo, trained with photometric features are found to consistently improve the accuracy of photometric redshifts.

The physics behind the origin and composition of the Cosmic Microwave Background (CMB) is a well-established topic in the field of Cosmology. Literature on CMB anisotropies reveal consistency with Gaussianity, but these were conducted on full multi-frequency temperature maps. In this thesis, we utilise clustering algorithms to specifically conduct statistical analyses on the distribution of hotspots in the CMB. We describe a series of data processing and clustering methodologies conducted, with results that conclusively show that the counts-in-cells distribution of hotspots in the CMB does not follow a Poisson distribution. Rather, the distribution exhibits a much closer fit to both the Negative Binomial Distribution (NBD) and the Gravitational Quasi-Equilibrium Distribution (GQED). From this result, we conclude that structure likely existed in the early universe, from the period of the recombination Epoch, possibly opening new insights in the field of galaxy formation.

We take advantage of the red clump stars to build the relation of the optical extinction (AV) and distance in each direction of supernova remnants (SNRs) with known extinction in the fourth Galactic quadrant. The distances of 9 SNRs are well determined by this method. Their uncertainties range from 10% to 30%, which is significantly improved for 8 SNRs, G279.0+1.1, G284.3-1.8, G296.1-0.5, G299.2-2.9, G308.4-1.4, G309.2-0.6, G309.8- 2.6, G332.4-0.4. In addition, SNR G284.3-1.8 with the new distance of 5.5 kpc is not likely associated with the PSR J1016-5857 at 3 kpc.

We have carried out the largest and most unbiased search for hypercompact (HC) HII regions. Our method combines four interferometric radio continuum surveys (THOR, CORNISH, MAGPIS and White2005) with far-infrared and sub-mm Galactic Plane surveys to identify embedded HII regions with positive spectral indices. 120 positive spectrum HII regions have been identified from a total sample of 534 positive spectral index radio sources. None of these HII regions, including the known HCHII regions recovered in our search, fulfills the canonical definition of an HCHII region at 5 GHz. We suggest that the current canonical definition of HCHII regions is not accurate and should be revised to include a hierarchical structure of ionized gas that results in an extended morphology at 5 GHz. Correlating our search with known ultracompact (UC) HII region surveys, we find that roughly half of detected UCHII regions have positive spectral indices, instead of more commonly assumed flat and optically thin spectra. This implies a mix of optically thin and thick emission and has important implications for previous analyses which have so far assumed optically thin emission for these objects. Positive spectrum HII regions are statistically more luminous and possess higher Lyman continuum fluxes than HII regions with flat or negative indices. Positive spectrum HII regions are thus more likely to be associated with more luminous and massive stars. No differences are found in clump mass, linear diameter or luminosity-to-mass ratio between positive spectrum and non-positive spectrum HII regions.

We carry out a project to independently measure the distances of supernova remnants (SNRs) in the first quadrant of the Galaxy. In this project, red clump (RC) stars are used as standard candles and extinction probes to build the optical extinction (A$_V$) - distance(D) relation in each direction of extinction-known SNRs. 15 SNRs' distances are well determined. Among them, the distances of G65.8-0.5, G66.0-0.0 and G67.6+0.9 are given for the first time. We also obtain 32 upper/lower limits of distances, and the distances to G5.7-0.1, G15.1-1.6, G28.8+1.5 and G78.2+2.1 are constrained. Most of the distances measured by the RC method are consistent with previous results. The RC method provides an independent access to the distances of SNRs.

We have undertaken the largest survey for outflows within the Galactic Plane using simultaneously observed 13CO and C18O data. 325 out of a total of 919 ATLASGAL clumps have data suitable to identify outflows, and 225 (69+-3%) of them show high velocity outflows. The clumps with detected outflows show significantly higher clump masses (M_{clump}), bolometric luminosities (L_{bol}), luminosity-to-mass ratios (L_{bol}/M_{clump}) and peak H_2 column densities (N_{H_2}) compared to those without outflows. Outflow activity has been detected within the youngest quiescent clump (i.e.,70um weak) in this sample and we find that the outflow detection rate increases with M_{clump},L_{bol},L_{bol}/M_{clump} and N_{H_2},approaching 90% in some cases(uchii regions=93+-3%;masers=86+-4%;hchii regions=100%). This high detection rate suggests that outflows are ubiquitous phenomena of massive star formation. The mean outflow mass entrainment rate implies a mean accretion rate of ~10^{-4}M_\odot\,yr^{-1}, in full agreement with the accretion rate predicted by theoretical models of massive star formation. Outflow properties are tightly correlated with M_{clump},L_{bol} and L_{bol}/M_{clump},and show the strongest relation with the bolometric clump luminosity. This suggests that outflows might be driven by the most massive and luminous source within the clump. The correlations are similar for both low-mass and high-mass outflows over 7 orders of magnitude, indicating that they may share a similar outflow mechanism. Outflow energy is comparable to the turbulent energy within the clump, however, we find no evidence that outflows increase the level of clump turbulence as the clumps evolve. This implies that the origin of turbulence within clumps is fixed before the onset of star formation.

We construct \HI~absorption spectra for 18 planetary nebulae (PNe) and their background sources using the data from the International Galactic Plane Survey. We estimate the kinematic distances of these PNe, among which 15 objects' kinematic distances are obtained for the first time. The distance uncertainties of 13 PNe range from 10% to 50%, which is a significant improvement with uncertainties of a factor two or three smaller than most of previous distance measurements. We confirm that PN G030.2-00.1 is not a PN because of its large distance found here.

Pulsar nulling is a phenomenon of sudden cessation of pulse emission for a number of periods. The nulling fraction was often used to characterize the phenomenon. We propose a new method to analyse pulsar nulling phenomenon, by involving two key parameters, the nulling degree, $\chi$, which is defined as the angle in a rectangular coordinates for the numbers of nulling periods and bursting periods, and the nulling scale, $ N $, which is defined as the effective length of the consecutive nulling periods and bursting periods. The nulling degree $\chi$ can be calculated by $\tan \chi = N_{\rm nulling} / N_{\rm bursting} $ and the mean is related to the nulling fraction, while the nulling scale, $ N $, is also a newly defined fundamental parameter which indicates how often the nulling occurs. We determined the distributions of $\chi$ and $ N $ for 10 pulsars by using the data in literature. We found that the nulling degree $\chi$ indicates the relative length of nulling to that of bursting, and the nulling scale $ N $ is found to be related to the derivative of rotation frequency and hence the loss rate of rotational energy of pulsars. Their deviations reflect the randomness of the nulling process.

We identify low redshift clusters and groups in the Sloan Digital Sky Survey (SDSS) and estimate their kinetic and correlation potential energies. We compare the distribution of these energies to the predictions by Yang and Saslaw (2012) and in the process estimate a measure of an average 3-dimensional velocity and spatial anisotropy of a sample of clusters. We find that the inferred velocity anisotropy is correlated with the inferred spatial anisotropy. We also find that the general shape of the energy distribution agrees with theory over a wide range of scales from small groups to superclusters once the uncertainties and fluctuations in the estimated energies are included.

We develop a general theory for estimating the probability that a galaxy cluster of a given shape exists. The theory is based on the observed result that the distribution of galaxies is very close to quasi-equilibrium, in both its linear and nonlinear regimes. This places constraints on the spatial configuration of a cluster of galaxies in quasi-equilibrium. In particular, we show that that a cluster of galaxies may be described as a collection of nearly virialized subclusters of approximately the same mass. Clusters that contain more than 10 subclusters are very likely to be completely virialized. Using our theory, we develop a method for comparing probabilities of different spatial configurations of subclusters. As an illustrative example, we show that a cluster of galaxies arranged in a line is more likely to occur than a cluster of galaxies arranged in a ring.

We determine the galaxy counts-in-cells distribution from the Sloan Digital Sky Survey (SDSS) for 3D spherical cells in redshift space as well as for 2D projected cells. We find that cosmic variance in the SDSS causes the counts-in-cells distributions in different quadrants to differ from each other by up to 20%. We also find that within this cosmic variance, the overall galaxy counts-in-cells distribution agrees with both the gravitational quasi-equilibrium distribution and the negative binomial distribution. We also find that brighter galaxies are more strongly clustered than if they were randomly selected from a larger complete sample that includes galaxies of all luminosities. The results suggest that bright galaxies could be in dark matter haloes separated by less than ~10 Mpc/h.

We examine the evolution of the spatial counts-in-cells distribution of galaxies and show that the form of the galaxy distribution function does not change significantly as galaxies merge and evolve. In particular, bound merging pairs follow a similar distribution to that of individual galaxies. From the adiabatic expansion of the universe we show how clustering, expansion and galaxy mergers affect the clustering parameter b. We also predict the evolution of b with respect to redshift.

Proper interpretation and understanding of astronomical data requires good knowledge of the data acquisition process. The increase in remote observing, queue observing, and the availability of large archived data products risk insulating astronomers from the telescope, potentially reducing their familiarity with the observational techniques crucial in understanding the data. Learning fundamental observing techniques can be done in at least three ways: 1) College and university operated observing facilities, 2) Student involvement in national facilities through competitive proposals, 3) Programs at national facilities to increase upper-level undergraduate and graduate student exposure to telescopes. We encourage both national organizations and universities to include programs and funding aimed at supporting hands-on experience with telescopes through the three methods mentioned.

The cosmological many-body problem is effectively an infinite system of gravitationally interacting masses in an expanding universe. Despite the interactions' long-range nature, an analytical theory of statistical mechanics describes the spatial and velocity distribution functions which arise in the quasi-equilibrium conditions that apply to many cosmologies. Consequences of this theory agree well with the observed distribution of galaxies. Further consequences such as thermodynamics provide insights into the physical properties of this system, including its robustness to mergers, and its transition from a grand canonical ensemble to a collection of microcanonical ensembles with negative specific heat.