We investigate the differences in the stellar population properties, the structure, and the environment between massive compact star-forming galaxies (cSFGs) with or without active galactic nucleus (AGN) at $2<z<3$ in the five 3D-HST/CANDELS fields. In a sample of 221 massive cSFGs, we constitute the most complete AGN census so far, identifying 66 AGNs by the X-ray detection, the mid-infrared color criterion, and/or the SED fitting, while the rest (155) are non-AGNs. Further dividing these cSFGs into two redshift bins, i.e., $2<z<2.5$ and $2.5 \leq z<3$, we find that in each redshift bin the cSFGs with AGNs have similar distributions of the stellar mass, the specific star formation rate, and the ratio of $L_{\rm IR}$ to $L_{\rm UV}$ to those without AGNs. After having performed a two-dimensional surface brightness modeling for those cSFGs with X-ray-detected AGNs (37) to correct for the influence of the central point-like X-ray AGN on measuring the structural parameters of its host galaxy, we find that in each redshift bin the cSFGs with AGNs have comparable distributions of all concerned structural parameters, i.e., the Sersic index, the 20\%-light radius, the Gini coefficient, and the concentration index, to those without AGNs. With a gradual consumption of available gas and dust, the structure of cSFGs, indicated by the above structural parameters, seem to be slightly more concentrated with decreasing redshift. At $2<z<3$, the similar environment between cSFGs with and without AGNs suggests that their AGN activities are potentially triggered by internal secular processes, such as gravitational instabilities or/and dynamical friction.

Galaxy properties are known to correlate most tightly with the galaxy effective stellar velocity dispersion $\sigma_{\rm e}$. Here we look for {\em additional} trends at fixed $\sigma_{\rm e}$ using 1339 galaxies ($M_\ast \gtrsim 6\times10^9$ M$_\odot$) with different morphologies in the MaNGA (DR14) sample with integral-field spectroscopy data. We focus on the gradients ($\gamma_{\rm rms} \equiv \sigma(R_{\rm e}/4)/\sigma_{\rm e}$) of the stellar root-mean-square velocity ($V_{\rm rms} \equiv \sqrt{V^2 + \sigma^2}$), which we show traces the total mass density gradient $\gamma_{\rm tot}$ derived from dynamical models and, more weakly, the bulge fraction. We confirm that $\gamma_{\rm rms}$ increases with $\sigma_{\rm e}$, age and metallicity. We additionally find that these correlations still exist at fixed $\sigma_{\rm e}$, where galaxies with larger $\gamma_{\rm rms}$ are found to be older and more metal-rich. It means that mass density gradients contain information of the stellar population which is not fully accounted for by $\sigma_{\rm e}$. This result puts an extra constraint on our understanding of galaxy quenching. We compare our results with galaxies in the IllustrisTNG hydrodynamical simulations and find that, at fixed $\sigma_{\rm e}$, similar trends exist with age, the bulge fraction, and the total mass density slope but, unlike observations, no correlation with metallicity can be detected in the simulations.

It is believed that massive galaxies have quenched their star formation because of active galactic nucleus feedback. However, recent studies have shown that some massive galaxies are still forming stars. We analyze the morphology of star formation regions for galaxies of stellar mass larger than 10$^{11.3}$ M$_{\odot}$ at around redshift $z_r=0.5$ using $u-z$ color images. We find that about $20\%$ of the massive galaxies are star-forming (SF) galaxies, and most of them ($\sim 85\%$) have asymmetric structures induced by recent mergers. Moreover, for these asymmetric galaxies, we find that the asymmetry of the SF regions becomes larger for bluer galaxies. Using the Illustris simulation, we can qualitatively reproduce the observed relation between asymmetry parameter and color. Furthermore, using the merger trees in the simulation, we find a correlation between the color of the main branch galaxies at $z_r=0.5$ and the sum of the Star Formation Rates (SFRs) of the recently accreted galaxies, which implies that star formation of the accreted galaxies has contributed to the observed star formation of the massive (host) galaxies (ex situ star formation). Furthermore, we find two blue and symmetric galaxies, candidates for massive blue disks, in our observed sample, which indicates that about $\sim 10\%$ of massive SF galaxies are forming stars in the normal mode of disk star formation (in situ star formation). With the simulation, we find that the disk galaxies at $z_r \approx 0.5$ should have experienced few major mergers during the last 4.3 Gyrs.

Dark matter could take the form of dark massive compact halo objects (dMACHOs); i.e., composite objects that are made up of dark-sector elementary particles, that could have a macroscopic mass from the Planck scale to above the solar mass scale, and that also admit a wide range of energy densities and sizes. Concentrating on the gravitational interaction of dMACHOs with visible matter, we map out the mass-radius parameter space that is consistent with gravitational lensing experiments, as well as anisotropies of the cosmic microwave background (CMB) based on the spherical accretion of matter onto a dMACHO in the hydrostatic approximation. For dMACHOs with a uniform-density mass profile and total mass in the range of $\sim 10^{-12} - 10\,M_\odot$, we find that a dMACHO could explain 100% of the dark matter if its radius is above $\approx 3$ times the Einstein radius of the lensing system. For a larger mass above $10\,M_\odot$, a dMACHO with radius above $\sim 1 \times 10^8 \mbox{cm} \times(M/100\,M_\odot)^{9/2}$ is consistent with CMB observables. For a lighter dMACHO with mass below $\sim 10^{-12}\,M_\odot$, there still is not a good experimental probe. Finally, we point out that heavier dMACHOs with masses $\sim 0.1\,M_\odot$ may be observed by X-ray and optical telescopes if they reside at rest in a large molecular cloud, nearby to our solar system, and accrete ordinary matter to emit photons.

We present a study on structures and physical properties of massive ($M_* >10^{10} M_{\sun} $) compact galaxies at $1.0<z<2.0$ in five 3D--{\it HST}/CANDELS fields. Compared with the extended star-forming galaxies (eSFGs), compact star-forming galaxies (cSFGs) are found to have the lower level of star formation, and mainly distribute in the quiescent region of the {\it UVJ} diagram. The distributions of dust attenuation and S{\'e}rsic index support that the progenitors of cQGs are cSFGs, and cSFGs are at a transitional phase between eSFGs and cQGs. The prevalence of X-ray selected AGNs ($\sim 28\%$) is confirmed in the cSFGs at $1<z<2$ which indicates that the violent gas-rich processes such as merger and disk instability could drive the structure to be more compact, and trigger both star formation and black hole growth in the central regions. Our results support the "two-step" scenario that the cSFGs at $1<z<2$ are the intermediate population after compaction but before a quick quenching. Our analysis of parametric and nonparametric morphologies shows that cQGs (eQGs) are more concentrated and have less substructures than cSFGs (eSFGs), and quenching and compactness should be associated with each other. The cSFGs at $1.5<z<2$ ($1<z<1.5$) prefer to be in higher (lower) density environment, similar as cQGs (eSFGs). It suggests that merger or strong interaction might be the main driving mechanism of compaction at higher redshifts, whereas the disk instability of individual galaxies might play a more important role on the formation of cSFGs at lower redshifts.

We discuss two possible scenarios, namely the curvaton mechanism and the dark matter density modulation, where non-Gaussianity signals of superheavy dark matter produced by gravity can be enhanced and observed. In both scenarios, superheavy dark matter couples to an additional light field as a mediator. In the case of derivative coupling, the resulting non-Gaussianities induced by the light field can be large, which can provide inflationary evidences for these superheavy dark matter scenarios.

We investigate the fundamental plane (FP) evolution of early-type galaxies in the IllustrisTNG-100 simulation (TNG100) from redshift $z=0$ to $z=2$. We find that a tight plane relation already exists as early as $z=2$. Its scatter stays as low as $\sim 0.08$ dex across this redshift range. Both slope parameters $b$ and $c$ (where $R \propto \sigma^b I^c$ with $R$, $\sigma$, and $I$ being the typical size, velocity dispersion, and surface brightness) of the plane evolve mildly since $z=2$, roughly consistent with observations. The FP residual $\rm Res$ ($\equiv\,a\,+\,b\log \sigma\,+\,c\log I\,-\,\log R$, where $a$ is the zero point of the FP) is found to strongly correlate with stellar age, indicating that stellar age can be used as a crucial fourth parameter of the FP. However, we find that $4c+b+2=\delta$, where $\delta \sim 0.8$ for FPs in TNG, rather than zero as is typically inferred from observations. This implies that a tight power-law relation between the dynamical mass-to-light ratio $M_{\rm dyn}/L$ and the dynamical mass $M_{\rm dyn}$ (where $M_{\rm dyn}\equiv 5\sigma^2R/G$, with $G$ being the gravitational constant) is not present in the TNG100 simulation. Recovering such a relation requires proper mixing between dark matter and baryons, as well as star formation occurring with correct efficiencies at the right mass scales. This represents a powerful constraint on the numerical models, which has to be satisfied in future hydrodynamical simulations.

"Dark quark nuggets", a lump of dark quark matter, can be produced in the early universe for a wide range of confining gauge theories and serve as a macroscopic dark matter candidate. The two necessary conditions, a nonzero dark baryon number asymmetry and a first-order phase transition, can be easily satisfied for many asymmetric dark matter models and QCD-like gauge theories with a few massless flavors. For confinement scales from 10 keV to 100 TeV, these dark quark nuggets with a huge dark baryon number have their masses vary from $10^{23}~\mathrm{g}$ to $10^{-7}~\mathrm{g}$ and their radii from $10^{8}~\mathrm{cm}$ to $10^{-15}~\mathrm{cm}$. Such macroscopic dark matter candidates can be searched for by a broad scope of experiments and even new detection strategies. Specifically, we have found that the gravitational microlensing experiments can probe heavier dark quark nuggets or smaller confinement scales around 10 keV; collision of dark quark nuggets can generate detectable and transient electromagnetic radiation signals; the stochastic gravitational wave signals from the first order phase transition can be probed by the pulsar timing array observations and other space-based interferometry experiments; the approximately massless dark mesons can behave as dark radiation to be tested by the next-generation CMB experiments; the free dark baryons, as a subcomponent of dark matter, can have direct detection signals for a sufficiently strong interaction strength with the visible sector.

To investigate the mass dependence of structural transformation and star formation quenching, we construct three galaxy samples using massive ($M_* > 10^{10} M_{\odot}$) red, green, and blue galaxy populations at $0.5<z<2.5$ in five 3D--{\it HST}/CANDELS fields. The structural parameters, including effective radius ($r_{\rm e}$), galaxy compactness ($\Sigma_{1.5}$), and second order moment of 20\% brightest pixels ($M_{20}$) are found to be correlated with stellar mass. S\'{e}rsic index ($n$), concentration ($C$), and Gini coefficient ($G$) seem to be insensitive to stellar mass. The morphological distinction between blue and red galaxies is found at a fixed mass bin, suggesting that quenching processes should be accompanied with transformations of galaxy structure and morphology. Except for $r_e$ and $\Sigma_{1.5}$ at high mass end, structural parameters of green galaxies are intermediate between red and blue galaxies in each stellar mass bin at $z < 2$, indicating green galaxies are at a transitional phase when blue galaxies are being quenched into quiescent statuses. The similar sizes and compactness for the blue and green galaxies at high-mass end implies that these galaxies will not appear to be significantly shrunk until they are completely quenched into red QGs. For the green galaxies at $0.5<z<1.5$, a morphological transformation sequence of bulge buildup can be seen as they are gradually shut down their star formation activities, while a faster morphological transformation is verified for the green galaxies at $1.5<z<2.5$.

Magnetohydrodynamic (MHD) theory and simulations have shown that reconnection is triggered via a fast "ideal" tearing instability in current sheets whose inverse aspect ratio decreases to $a/L\sim S^{-1/3}$, with $S$ is the Lundquist number defined by the half-length $L$ of the current sheet (of thickness $2a$). Ideal tearing, in 2D sheets, triggers a hierarchical collapse via stretching of X-points and recursive instability. At each step, the local Lundquist number decreases, until the subsequent sheet thickness either approaches kinetic scales or the Lundquist number becomes sufficiently small. Here we carry out a series of Hall-MHD simulations to show how the Hall effect modifies recursive reconnection once the ion inertial scale is approached. We show that as the ion inertial length becomes of the order of the inner, singular layer thickness at some step of the recursive collapse, reconnection transits from the plasmoid-dominant regime to an intermediate plasmoid+Hall regime and then to the Hall-dominant regime. The structure around the X-point, the reconnection rate, the dissipation property and the power spectra are also modified significantly by the Hall effect.

The quantum fluctuations of the Higgs field during inflation could be a source of primordial density perturbations through Higgs-dependent inflaton decay. By measuring primordial non-Gaussianities, this so-called Higgs-modulated reheating scenario provides us a unique chance to probe Higgs interactions at extremely high energy scale, which we call the Cosmological Higgs Collider (CHC). We realize CHC in a simple scenario where the inflaton decays into Higgs-portal scalars, taking account of the decay of the Higgs fluctuation amplitude after inflation. We also calculate the CHC signals of Standard Model particles, namely their imprints in the squeezed bispectrum, which can be naturally large. The concept of CHC can be straightforwardly generalized to cosmological isocurvature colliders with other types of isocurvature perturbations.

Based on a large sample of massive ($M_{*}\geq 10^{10} M_{\odot}$) compact galaxies at $1.0 < z < 3.0$ in five 3D-HST/CANDELS fields, we quantify the fractional abundance and comoving number density of massive compact galaxies as a function of redshift. The samples of compact quiescent galaxies (cQGs) and compact star-forming galaxies (cSFGs) are constructed by various selection criteria of compact galaxies in literatures, and the effect of compactness definition on abundance estimate is proved to be remarkable, particularly for the cQGs and cSFGs at high redshifts. Regardless of the compactness criteria adopted, their overall redshift evolutions of fractional abundance and number density are found to be rather similar. Large samples of the cQGs exhibit a sustaining increase in number density from $z \sim 3$ to 2 and a plateau at $1<z<2$. For massive cSFGs, a plateau in the number density at $2<z<3$ can be found, as well as a continuous drop from $z \sim 2$ to 1. The evolutionary trends of the cQG and cSFG abundances support the scenario that the cSFGs at $z \geq 2$ may have been rapidly quenched into quiescent phase via violent dissipational processes such as major merger and disk instabilities. Rarity of the cSFGs at lower redshifts ($z < 1$) can be interpreted by the decrease of gas reservoirs in dark matter halos and the consequent low efficiency of gas-rich dissipation.

Motivated by the recently discovered kinematic "Hubble sequence" shown by the stellar orbit-circularity distribution of 260 CALIFA galaxies, we make use of a comparable galaxy sample at z = 0 with a stellar mass range from 5E9 to 5E11 solar masses, selected from the IllustrisTNG simulation and study their stellar orbit compositions in relation to a number of other fundamental galaxy properties.We find that the TNG100 simulation broadly reproduces the observed fractions of different orbital components and their stellar mass dependencies. In particular, the mean mass dependencies of the luminosity fractions for the kinematically warm and hot orbits are well reproduced within model uncertainties of the observed galaxies. The simulation also largely reproduces the observed peak and trough features at a stellar mass of 1-2E10 solar masses, in the mean distributions of the cold- and hot-orbit fractions, respectively, indicating fewer cooler orbits and more hotter orbits in both more- and less-massive galaxies beyond such a mass range. Several marginal disagreements are seen between the simulation and observations: the average cold-orbit (counter-rotating) fractions of the simulated galaxies below (above) a stellar mass of 6E10 solar masses, are systematically higher than the observational data by < 10% (absolute orbital fraction); the simulation also seems to produce more scatter for the cold-orbit fraction and less so for the non-cold orbits at any given galaxy mass. Possible causes that stem from the adopted heating mechanisms are discussed.

The binary black hole merger events observed by the Advanced LIGO (aLIGO) and VIRGO collaboration can shed light on the origins of black holes. Many studies based on black hole stellar origins have shown a maximum mass for stellar black holes, which can be measured or constrained from the observed black hole mass distribution. In this paper, we point out that the redshift distribution of the observed merger events can provide complementary information for studying the black hole mass distribution, because the detectability correlates the event redshift to the black hole masses. Based on the five observed events and using the Kolmogorov-Smirnov test, we have found that the maximum of the stellar black hole masses are constrained to be below $76\,M_\odot$ at 90\% confidence level, for a negative power-law index of 2.3 for the heavier black hole. With the improved sensitivity of aLIGO, a few dozen merger events may be obtained, for which we estimate that the maximum mass will be constrained to $10\,M_\odot$ accuracy.

We derive the mass weighted total density slopes within the effective (half-light) radius, $\gamma'$, for more than 2000 nearby galaxies from the SDSS-IV MaNGA survey using Jeans-anisotropic-models applied to IFU observations. Our galaxies span a wide range of the stellar mass ($10^9$ $M_{\rm \odot}< M_* < 10^{12}$ M$_{\odot}$) and the velocity dispersion (30 km/s $< \sigma_v <$ 300 km/s). We find that for galaxies with velocity dispersion $\sigma_v>100$ km/s, the density slope has a mean value $\langle \gamma^{\prime} \rangle = 2.24$ and a dispersion $\sigma_{\gamma}=0.22$, almost independent of velocity dispersion. A clear turn over in the $\gamma'-\sigma_v$ relation is present at $\sigma\sim 100$ km/s, below which the density slope decreases rapidly with $\sigma_v$. Our analysis shows that a large fraction of dwarf galaxies (below $M_* = 10^{10}$ M$_{\odot}$) have total density slopes shallower than 1, which implies that they may reside in cold dark matter halos with shallow density slopes. We compare our results with that of galaxies in hydrodynamical simulations of EAGLE, Illustris and IllustrisTNG projects, and find all simulations predict shallower density slopes for massive galaxies with high $\sigma_v$. Finally, we explore the dependence of $\gamma'$ on the positions of galaxies in halos, namely centrals vs. satellites, and find that for the same velocity dispersion, the amplitude of $\gamma'$ is higher for satellite galaxies by about 0.1.

We show how cosmological correlation functions of massless fields can be rewritten in terms of Minkowski correlation functions, by extracting symmetry-breaking operators from the cosmological correlators. This technique simplifies some cosmological calculations. Also, known properties of Minkowski correlation functions can be translated to non-trivial properties of cosmological correlations. To illustrate this idea, inflation to Minkowski and matter bounce to Minkowski relations are presented for the interactions of general single field inflation. And a Minkowski recursion relation is translated to a novel relation for inflation.

We construct a thermal dark matter model with annihilation mediated by a resonance to explain the positron excess observed by PAMELA, Fermi-LAT and AMS-02, while satisfying constraints from cosmic microwave background (CMB) measurements. The challenging requirement is that the resonance has twice the dark matter mass to one part in a million. We achieve this by introducing an $SU(3)_f$ dark flavor symmetry that is spontaneously broken to $SU(2)_f \times U(1)_f$. The resonance is the heaviest state in the dark matter flavor multiplet and the required mass relation is protected by the vacuum structure and supersymmetry from radiative corrections. The pseudo-Nambu Goldstone Bosons (PNGB's) from the dark flavor symmetry breaking can be slightly lighter than one GeV and dominantly decay into two muons just from kinematics, with subsequent decay into positrons. The PNGB's are produced in resonant dark matter semi-annihilation, where two dark matter particles annihilate into an anti-dark matter particle and a PNGB. The dark matter mass in our model is constrained to be below around 1.9 TeV from fitting thermal relic abundance, AMS-02 data and CMB constraints. The superpartners of Standard Model (SM) particles can cascade decay into a light PNGB along with SM particles, yielding a correlated signal of this model at colliders. One of the interesting signatures is a resonance of a SM Higgs boson plus two collimated muons, which has superb discovery potential at LHC Run 2.

A characteristic signature of String Gas Cosmology is primordial power spectra for scalar and tensor modes which are almost scale-invariant but with a red tilt for scalar modes but a blue tilt for tensor modes. This feature, however, can also be realized in the so-called G-inflation model, in which Horndeski operators are introduced which leads to a blue tensor tilt by softly breaking the Null Energy Condition. In this article we search for potential observational differences between these two cosmologies by performing detailed perturbation analyses based on the Effective Field Theory approach. Our results show that, although both two models produce blue tilted tensor perturbations, they behave differently in three aspects. Firstly, String Gas Cosmology predicts a specific consistency relation between the index of the scalar modes ns and that of tensor ones nt, which is hard to be reproduced by G-inflation. Secondly, String Gas Cosmology typically predicts non-Gaussianities which are highly suppressed on observable scales, while G-inflation gives rise to observationally large non-Gaussianities because the kinetic terms in the action become important during inflation. However, after finely tuning the model parameters of G-inflation it is possible to obtain a blue tensor spectrum and negligible non-Gaussianities with a degeneracy between the two models. This degeneracy can be broken by a third observable, namely the scale dependence of the nonlinearity parameter, which vanishes for G-inflation but has a blue tilt in the case of String Gas Cosmology. Therefore, we conclude that String Gas Cosmology is in principle observationally distinguishable from the single field inflationary cosmology, even allowing for modifications such as G-inflation.

Using the data taken from Tibet II High Density (HD) Array (1997 February-1999 September) and Tibet-III array (1999 November-2005 November), our previous northern sky survey for TeV $\gamma-$ray point sources has now been updated by a factor of 2.8 improved statistics. From $0.0^{\circ}$ to $60.0^{\circ}$ in declination (Dec) range, no new TeV $\gamma-$ray point sources with sufficiently high significance were identified while the well-known Crab Nebula and Mrk421 remain to be the brightest TeV $\gamma-$ray sources within the field of view of the Tibet air shower array. Based on the currently available data and at the 90% confidence level (C.L.), the flux upper limits for different power law index assumption are re-derived, which are approximately improved by 1.7 times as compared with our previous reported limits.

The amplitude of the Compton-Getting (CG) anisotropy contains the power-law index of the cosmic-ray energy spectrum. Based on this relation and using the Tibet air-shower array data, we measure the cosmic-ray spectral index to be $-3.03 \pm 0.55_{stat} \pm < 0.62_{syst}$ between 6 TeV and 40 TeV, consistent with $-$2.7 from direct energy spectrum measurements. Potentially, this CG anisotropy analysis can be utilized to confirm the astrophysical origin of the ``knee'' against models for non-standard hadronic interactions in the atmosphere.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We are planning to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV (10 - 1000 TeV) energy region by discriminating them from cosmic-ray hadrons. We report on the possibility of detection of gamma rays in the 100 TeV energy region in our field of view, based on the improved sensitivity of our air shower array deduced from the full Monte Carlo simulation.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We have a plan to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV energy region by discriminating them from cosmic-ray hadrons. We have deduced the attainable sensitivity of the muon detector array using our Monte Carlo simulation. We report here on the detailed procedure of our Monte Carlo simulation.

The magmatic and metamorphic zircons were investigated with Raman spectrum microprobe analysis. We found notable differences between these two kinds of zircons exhibited by the variation trend of Raman peak intensity from core to rim of a crystal. In magmatic zircons, the intensity and the ratio H/W of Raman spectrum peaks gradually decrease from core to rim of a crystal, which is produced by an increase in metamictization degree and suggests an increase in U and Th concentrations from core to rim. In metamorphic zircons, there are two kinds of crystals according to their Raman spectra: the first group of zircons exhibits a variation trend opposite to those of magmatic zircons, tending to increase in the Raman peak intensity and H/W value from core to rim of a crystal, which is produced by a decrease in metamictization degree and indicates a decrease of U and Th concentrations from core to rim of a crystal. The second group of zircons exhibits no change in Raman peak intensity and H/W value through a crystal. The data of infrared and Raman spectra of these crystals show that they are well crystallized and have no lattice destruction induced by metamictization, and are thought to crystallize in high temperature stages of metamorphism. During these stages, the U and Th ions have been removed by metamorphic fluids from the parent rocks of these zircons. The other difference between magmatic and metamorphic zircons is the background level of their Raman spectra, which is high and sloped in magmatic zircons, but low and horizontal in metamorphic zircons. The differences between magmatic and metamorphic zircons can be used to identify the genesis of zircons and understand the origin and evolution history of their parent rocks.