sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Song, C.

Normalized to: Song, C.

20 article(s) in total. 139 co-authors, from 1 to 11 common article(s). Median position in authors list is 3,0.

[1]  oai:arXiv.org:2004.00043  [pdf] - 2077153
Neutrinos and gravitational waves from magnetized neutrino-dominated accretion discs with magnetic coupling
Comments: 10 pages, 7 figures, accepted for publication in MNRAS
Submitted: 2020-03-31
Gamma-ray bursts (GRBs) might be powered by a black hole (BH) hyperaccretion systems via the Blandford-Znajek (BZ) mechanism or neutrino annihilation from neutrino-dominated accretion flows (NDAFs). Magnetic coupling (MC) between the inner disc and BH can transfer angular momentum and energy from the fast-rotating BH to the disc. The neutrino luminosity and neutrino annihilation luminosity are both efficiently enhanced by the MC process. In this paper, we study the structure, luminosity, MeV neutrinos, and gravitational waves (GWs) of magnetized NDAFs (MNDAFs) under the assumption that both the BZ and MC mechanisms are present. The results indict that the BZ mechanism will compete with the neutrino annihilation luminosity to trigger jets under the different partitions of the two magnetic mechanisms. The typical neutrino luminosity and annihilation luminosity of MNDAFs are definitely higher than those of NDAFs. The typical peak energy of neutrino spectra of MNDAFs is higher than that of NDAFs, but similar to those of core-collapse supernovae. Moreover, if the MC process is dominant, then the GWs originating from the anisotropic neutrino emission will be stronger particularly for discs with high accretion rates.
[2]  oai:arXiv.org:2003.00121  [pdf] - 2081013
Large-scale multiconfiguration Dirac-Hartree-Fock calculations for astrophysics: n=4 levels in P-like ions from Mn~XI to Ni~XIV
Comments:
Submitted: 2020-02-28
Using the multiconfiguration Dirac-Hartree-Fock and the relativistic configuration interaction methods, a consistent set of transition energies and radiative transition data for the lowest 546 (623, 701, 745) states of the $3p^4 3d$, $3s 3p^2 3d^2$, $3s 3p^3 4p$, $3s 3p^4$, $3s^2 3d^3$, $3s^2 3p^2 3d$, $3s^2 3p^2 4d$, $3s^2 3p^2 4s$, $3p^3 3d^2$, $3p^5$, $3s 3p 3d^3$, $3s 3p^3 3d$, $3s 3p^3 4s$, $3s^2 3p 3d^2$, %$3s^2 3p^2 4f$, $3s^2 3p^2 4p$, $3s^2 3p^3$ configurations in Mn~XI (Fe~XII, Co~XIII, Ni~XIV) is provided. The comparison between calculated excitation energies for the $n=4$ states and available experimental values for Fe XII indicate that the calculations are highly accurate, with uncertainties of only a few hundred cm$^{-1}$. Lines from these states are prominent in the soft X-rays. With the present calculations, several recent new identifications are confirmed. Other identifications involving $3p^2 4d$ levels in Fe~XII that were found questionable are discussed and a few new assignments are recommended. As some $n=4$ states of the other ions also show large discrepancies between experimental and calculated energies, we reassess their identification. The present study provides highly accurate atomic data for the $n=4$ states of P-like ions of astrophysical interest, for which experimental data are scarce.
[3]  oai:arXiv.org:1905.04850  [pdf] - 1929780
Black hole hyperaccretion in collapsars. I. MeV neutrinos
Comments: 9 pages, 4 figures, accepted for publication in ApJ
Submitted: 2019-05-12
As the plausible central engine of gamma-ray bursts (GRBs), a black hole (BH) hyperaccretion disk should be in a state of neutrino-dominated accretion flow (NDAF) if the accretion rate is larger than the ignition rate of an NDAF. A rotating stellar-mass BH surrounded by a hyperaccretion disk might be born in the center of a massive core collapsar. In the initial hundreds of seconds of the accretion process, the mass supply rate of the massive progenitor is generally higher than the ignition accretion rate, but the jets are generally choked in the envelope. Thus, neutrinos can be emitted from the center of a core collapsar. In this paper, we study the effects of the masses and metallicities of progenitor stars on the time-integrated spectra of electron neutrinos from NDAFs. The peak energies of the calculated spectra are approximately 10-20 MeV. The mass of a collapsar has little influence on the neutrino spectrum, and a low metallicity is beneficial to the production of low-energy ($\lesssim$ 1 MeV) neutrinos. We also investigate the differences in the electron neutrino spectra between NDAFs and proto-neutron stars. Combining with the electromagnetic counterparts and multi-messenger astronomy, one may verify the possible remnants of the core collapse of massive stars with future neutrino detectors.
[4]  oai:arXiv.org:1902.03464  [pdf] - 1829737
A possible feedback mechanism of outflows from a black hole hyperaccretion disk in the center of jet-driven iPTF14hls
Comments: 15 pages, 2 figures, accepted for publication in Journal of High Energy Astrophysics
Submitted: 2019-02-09
iPTF14hls is an unusually bright, long-lived II-P supernova (SN), whose light curve has at least five peaks. We propose that the outflows from the black hole hyperaccretion systems in the center of the collapsars should continuously inject into the envelope. For a jet-driven core-collapsar model, the outflow feedback results in prolonging the accretion timescale and fluctuating accretion rates in our analytic solutions. Thus, the long period of luminous, varying SN iPTF14hls might originate from the choked jets, which are regulated by the feedback of the strong disk outflows in a massive core-collapsar. One can expect that jet-driven iPTF14hls may last no more than approximately 3,000 days, and the luminosity may quickly decrease in the later stages. Moreover, the double-peak light curves in some SNe might be explained by the outflow feedback mechanism.
[5]  oai:arXiv.org:1812.01708  [pdf] - 1838277
Black hole hyperaccretion inflow-outflow model. II. Long-duration gamma-ray bursts and supernova $\rm ^{56}Ni$ bumps
Comments: 9 pages, 1 table, 2 figures, accepted for publication in ApJ
Submitted: 2018-12-04
Long-duration gamma-ray bursts (LGRBs) associated with supernovae (SNe) are possibly born out of the death of a massive star. After the star collapses, a stellar-mass black hole (BH) is formed, surrounded by a hyperaccretion disk with outflows. Blandford-Znajek jets can be launched and then break out from the envelope to power LGRBs. The jet luminosity depends on the net inflow accretion rate at the inner radius of the disk. Furthermore, $\rm ^{56}Ni$ synthesis should occur in the strong outflows from the accretion disk. The decay of $\rm ^{56}Ni$ is considered to be the possible origin of SN bumps in the subsequent optical afterglows of LGRBs. If $\rm ^{56}Ni$ originates entirely from the outflows, there is competition between the luminosities of LGRBs and those of the corresponding $\rm ^{56}Ni$ bumps because of the material distribution between the disk inflows and outflows. In this paper, we investigated these two luminosities based on 15 cases of LGRB-SN in the framework of the BH hyperaccretion inflow-outflow model. Then, one can constrain the characteristics of the progenitor stars of these LGRBs. The results indicate that these LGRBs may originate from the low-metallicity ($Z\lesssim 10^{-2}Z_{\odot}$, where $Z$ and $Z_{\odot}$ are the metallicities of the stars and the Sun, respectively) stars or some massive solar-metallicity stars. For ultra-LGRBs (ULGRBs), such as GRB 111209A, most of the massive low-metallicity stars with $Z \lesssim 10^{-2}Z_{\odot}$ could be progenitors only if very strong outflows are launched from the disks. When the contributions of nucleosynthesis in the disk outflows are considered, there is no shortage of $\rm ^{56}Ni$ mass for luminous SNe associated with ULGRBs.
[6]  oai:arXiv.org:1710.00142  [pdf] - 1666011
Outflows from black hole hyperaccretion systems: short and long-short gamma-ray bursts and "quasi-supernovae"
Comments: 11 pages, 2 figures, 1 table, accepted for publication in MNRAS
Submitted: 2017-09-29, last modified: 2018-04-16
The detections of some long gamma-ray bursts (LGRBs) relevant to mergers of neutron star (NS)-NS or black hole (BH)-NS, as well as some short gamma-ray bursts (SGRBs) probably produced by collapsars, muddle the boundary of two categories of gamma-ray bursts (GRBs). In both cases, a plausible candidate of central engine is a BH surrounded by a hyperaccretion disc with strong outflows, launching relativistic jets driven by Blandford-Znajek mechanism. In the framework of compact binary mergers, we test the applicability of the BH hyperaccretion inflow-outflow model on powering observed GRBs. We find that, for a low outflow ratio, $\sim 50\%$, postmerger hyperaccretion processes could power not only all SGRBs but also most of LGRBs. Some LGRBs might do originate from merger events in the BH hyperaccretion scenario, at least on the energy requirement. Moreover, kilonovae might be produced by neutron-rich outflows, and their luminosities and timescales significantly depend on the outflow strengths. GRBs and their associated kilonovae are competitive with each other on the disc mass and total energy budgets. The stronger the outflow, the more similar the characteristics of kilonovae to supernovae (SNe). This kind of `nova' might be called `quasi-SN'.
[7]  oai:arXiv.org:1710.04355  [pdf] - 1641301
X-ray light curve in GRB 170714A: evidence for quark star?
Comments: 6 pages, 1 figure, 1 table, accepted for publication in ApJ
Submitted: 2017-10-11, last modified: 2018-01-29
Two plateaus and one following bump in the X-ray light curve of GRB 170714A have been detected by the \textit{Swift}/X-Ray Telescope, which could be very meaningful for the central engine of gamma-ray bursts (GRBs), implying that the origin of this burst might be different from that of other ultra-long GRBs. We propose that merging two neutron stars into a hyper-massive quark star (QS) and then collapsing into a black hole (BH), with a delay time around $10^4$~s, could be responsible for those X-ray components. The hyper-massive QS is initially in a fluid state, being turbulent and differentially rotating, but would be solidified and release its latent heat injected into the GRB fireball (lasting about $10^3$~s during the liquid-solid phase transition). Magnetic field as high as $\sim 10^{15}$~G could be created by dynamo action of the newborn liquid QS, and a magnetar-like central engine (after solidification) supplies significant energy for the second plateau. More energy could be released during a fall-back accretion after the post-merger QS collapses to a BH, and the X-ray bump forms. This post-merger QS model might be tested by future observations, with either advanced gravitational wave detectors (e.g., advanced LIGO and VIRGO) or X-ray/optical telescopes.
[8]  oai:arXiv.org:1801.09717  [pdf] - 1634370
Compact binary merger and kilonova: outflows from remnant disc
Comments: 8 pages, 6 figures, accepted for publication in MNRAS
Submitted: 2018-01-29
Outflows launched from a remnant disc of compact binary merger may have essential contribution to the kilonova emission. Numerical calculations are conducted in this work to study the structure of accretion flows and outflows. By the incorporation of limited-energy advection in the hyper-accretion discs, outflows occur naturally from accretion flows due to imbalance between the viscous heating and the sum of the advective and radiative cooling. Following this spirit, we revisit the properties of the merger outflow ejecta. Our results show that around $10^{-3} \sim 10^{-1} M_\odot$ of the disc mass can be launched as powerful outflows. The amount of unbound mass varies with the disc mass and the viscosity. The outflow-contributed peak luminosity is around $10^{40} \sim 10^{41}$ erg/s. Such a scenario can account for the observed kilonovae associated with short gamma-ray bursts, including the recent event AT2017gfo (GW170817).
[9]  oai:arXiv.org:1710.00141  [pdf] - 1614833
Black hole hyperaccretion inflow-outflow model. I. long and ultra-long gamma-ray bursts
Comments: 7 pages, 5 figures, accepted for publication in ApJ
Submitted: 2017-09-29, last modified: 2017-11-27
Long-duration gamma-ray bursts (LGRBs) and ultra-LGRBs (ULGRBs) originate from collapsars, in the center of which a newborn rotating stellar-mass black hole (BH) surrounded by a massive accretion disk may form. In the scenario of BH hyperaccretion inflow-outflow model and Blandford-Znajek (BZ) mechanism to trigger gamma-ray bursts (GRBs), the real accretion rate to power a BZ jet is far lower than the mass supply rate from the progenitor star. The characteristics of the progenitor stars can be constrained by GRB luminosity observations, and the results exceed usual expectations. LGRBs lasting from several seconds to tens of seconds in the rest frame may originate from solar-metallicity ($Z \sim 1~ Z_\odot$, where $Z$ and $Z_\odot$ are the metallicities of progenitor stars and the Sun), massive ($M \gtrsim 34 ~M_\odot$, where $M$ and $M_\odot$ are the masses of progenitor stars and the Sun) stars or some zero-metallicity ($Z \sim 0$) stars. A fraction of low-metallicity ($Z \lesssim 10^{-2}~Z_\odot$) stars, including Population III stars, can produce ULGRBs such as GRB 111209A. The fraction of LGRBs lasting less than tens of seconds in the rest frame is more than 40$\%$, which cannot conform to the fraction of the demanded type of progenitor star. It possibly implies that the activity timescale of central engine may be much longer than the observed timescale of prompt emission phase, as indicated by X-ray late-time activities. Alternatively, LGRBs and ULGRBs may be powered by a millisecond magnetar central engine.
[10]  oai:arXiv.org:1709.09810  [pdf] - 1593653
Comparison of gravitational waves from central engines of gamma-ray bursts: neutrino-dominated accretion flows, Blandford-Znajek mechanisms, and millisecond magnetars
Comments: 8 pages, 3 figures, accepted for publication in ApJ
Submitted: 2017-09-28, last modified: 2017-10-10
Neutrino-dominated accretion flow (NDAF) around a rotating stellar-mass black hole (BH) is one of the plausible candidates for the central engines of gamma-ray bursts (GRBs). Because the time-variant and anisotropic emission of neutrinos from NDAFs leads to GRB variability, NDAFs can be regarded as the sources of the strong gravitational waves (GWs). We calculate the dependences of the GW strains on both the BH spin and the accretion rate. We demonstrate that for typical GRBs with either single pulse or multiple pulses, the GWs from NDAFs might be detected at a distance of $\sim 100$ kpc/$\sim 1$ Mpc by advanced LIGO/Einstein Telescope with a typical frequency of $\sim 10-100$ Hz. Besides NDAFs, the other two competitive candidates for GRB central engine are Blandford-Znajek (BZ) mechanism and millisecond magnetars. We explore the GW signals from these two as well, and compare the corresponding results with NDAFs'. We find that for a certain GRB, the possible detected distance from NDAFs is about two orders of magnitude higher than that from BZ mechanism, but at least two orders of magnitude lower than that from magnetars. The typical GW frequency for BZ mechanism is the same with that of NDAFs, $\sim 10-100$ Hz, while the typical frequency for magnetars is $\sim 2000$ Hz. Therefore, the GWs released by the central engines of adjacent GRBs might be used to determine whether there is an NDAF, a BZ jet or a magnetar in GRB center.
[11]  oai:arXiv.org:1602.06691  [pdf] - 1370424
Testing black hole neutrino-dominated accretion discs for long-duration gamma-ray bursts
Comments: 6 pages, 1 figure, accepted for publication in MNRAS
Submitted: 2016-02-22
Long-duration gamma-ray bursts (LGRBs) are generally considered to originate from the massive collapsars. It is believed that the central engine of gamma-ray bursts (GRBs) is a neutrino-dominated accretion flow (NDAF) around a rotating stellar-mass black hole (BH). The neutrino annihilation above the NDAF is a feasible mechanism to power GRB. In this work, we analyse the distributions of the isotropic gamma-ray radiated energy and jet kinetic energy of 48 LGRBs. According to the NDAF and fireball models, we estimate the mean accreted masses of LGRBs in our sample to investigate whether the NDAFs can power LGRBs with the reasonable BH parameters and conversion efficiency of neutrino annihilation. The results indicate that most of the values of the accreted masses are less than $5~M_\odot$ for the extreme Kerr BHs and high conversion efficiency. It suggests that the NDAFs may be suitable for most of LGRBs except for some extremely high energy sources.
[12]  oai:arXiv.org:1512.08028  [pdf] - 1332641
Evolutions of stellar-mass black hole hyperaccretion systems in the center of gamma-ray bursts
Comments: 20 pages, 6 figures
Submitted: 2015-12-25
A neutrino-dominated accretion disk around a stellar-mass black hole (BH) can power a gamma-ray burst (GRB) via annihilation of neutrinos launched from the disk. For the BH hyperaccretion system, high accretion rate should trigger the violent evolution of the BH's characteristics, which further leads to the evolution of the neutrino annihilation luminosity. In this paper, we consider the evolution of the accretion system to analyze the mean time-dependent neutrino annihilation luminosity with the different mean accretion rates and initial BH parameters. By time-integrating the luminosity, the total neutrino annihilation energy with the reasonable initial disk mass can satisfy the most of short-duration GRBs and about half of long-duration GRBs. Moreover, the extreme Kerr BH should exist in the cental engines of some high-luminosity GRBs. GRBs with higher energy have to request the alternative magnetohydrodynamics processes in the centers, such as the Blandford-Znajek jet from the accretion system or the millisecond magnetar.
[13]  oai:arXiv.org:1511.01243  [pdf] - 1392739
Structures of GMC W 37
Comments: accepted for publication in RAA; any comment is welcome
Submitted: 2015-11-04
We carried out observations toward the giant molecular cloud W 37 with the $J = 1 - 0$ transitions of $^{12}$CO, $^{13}$CO, and C$^{18}$O using the 13.7 m single-dish telescope at the Delingha station of Purple Mountain Observatory. Based on the three CO lines, we calculated the column densities, cloud masses for the molecular clouds with radial velocities at around $+20 \mathrm{km s}^{-1}$. The gas mass of W 37, calculated from $^{13}$CO emission, is $1.7\times10^5 M_\odot$, above the criteria of giant molecular cloud. The dense ridge of W 37 is a dense filament, which is supercritical in linear mass ratio. Dense clumps found by C$^{18}$O emission are aligned along the dense ridge with a regular interval about 2.8 pc, similar to the clump separation caused by large-scale `sausage instability'. We confirm the identification of the giant molecular filament (GMF) G 18.0-16.8 by \cite{2014A&A...568A..73R} and find a new giant filament, G16.5-15.8, located in the west 0.8 degree of G 18.0-16.8. Both GMFs are not gravitationally bound, as indicated by their low linear mass ratio ($\sim80 M_\odot \mathrm{pc}^{-1}$). We compared the gas temperature map with the dust temperature map from \emph{Herschel} images, and find similar structures. The spatial distributions of class I objects and the dense clumps is reminiscent of triggered star formation occurring in the northwestern part of W 37, which is close to NGC 6611.
[14]  oai:arXiv.org:0910.4184  [pdf] - 157172
Indications of Proton-Dominated Cosmic Ray Composition above 1.6 EeV
Comments: 11 pages, 4 figures, accepted for publication in Physical Review Letters
Submitted: 2009-10-21, last modified: 2010-04-03
We report studies of ultra-high energy cosmic ray composition via analysis of depth of airshower maximum (Xmax), for airshower events collected by the High Resolution Fly's Eye (HiRes) observatory. The HiRes data are consistent with a constant elongation rate d<Xmax>/d(log(E)) of 47.9 +- 6.0 (stat.) +- 3.2 (syst.)g/cm^2/decade for energies between 1.6 EeV and 63 EeV, and are consistent with a predominantly protonic composition of cosmic rays when interpreted via the QGSJET01 and QGSJET-II high-energy hadronic interaction models. These measurements constrain models in which the galactic-to-extragalactic transition is the cause of the energy spectrum "ankle' at 4 EeV.
[15]  oai:arXiv.org:1002.1444  [pdf] - 1025033
Analysis of large-scale anisotropy of ultra-high energy cosmic rays in HiRes data
Comments: 9 pages, 5 Postscript figures
Submitted: 2010-02-07
Stereo data collected by the HiRes experiment over a six year period are examined for large-scale anisotropy related to the inhomogeneous distribution of matter in the nearby Universe. We consider the generic case of small cosmic-ray deflections and a large number of sources tracing the matter distribution. In this matter tracer model the expected cosmic ray flux depends essentially on a single free parameter, the typical deflection angle theta. We find that the HiRes data with threshold energies of 40 EeV and 57 EeV are incompatible with the matter tracer model at a 95% confidence level unless theta is larger than 10 degrees and are compatible with an isotropic flux. The data set above 10 EeV is compatible with both the matter tracer model and an isotropic flux.
[16]  oai:arXiv.org:0804.0382  [pdf] - 11438
Search for Correlations between HiRes Stereo Events and Active Galactic Nuclei
Comments: 13 pages, 1 table, 5 figures
Submitted: 2008-04-02, last modified: 2008-08-15
We have searched for correlations between the pointing directions of ultrahigh energy cosmic rays observed by the High Resolution Fly's Eye experiment and Active Galactic Nuclei (AGN) visible from its northern hemisphere location. No correlations, other than random correlations, have been found. We report our results using search parameters prescribed by the Pierre Auger collaboration. Using these parameters, the Auger collaboration concludes that a positive correlation exists for sources visible to their southern hemisphere location. We also describe results using two methods for determining the chance probability of correlations: one in which a hypothesis is formed from scanning one half of the data and tested on the second half, and another which involves a scan over the entire data set. The most significant correlation found occurred with a chance probability of 24%.
[17]  oai:arXiv.org:0803.0554  [pdf] - 10663
An upper limit on the electron-neutrino flux from the HiRes detector
Comments: 13 pages, 3 figures. submitted to Astrophysical Journal
Submitted: 2008-03-04, last modified: 2008-06-03
Air-fluorescence detectors such as the High Resolution Fly's Eye (HiRes) detector are very sensitive to upward-going, Earth-skimming ultrahigh energy electron-neutrino-induced showers. This is due to the relatively large interaction cross sections of these high-energy neutrinos and to the Landau-Pomeranchuk-Migdal (LPM) effect. The LPM effect causes a significant decrease in the cross sections for bremsstrahlung and pair production, allowing charged-current electron-neutrino-induced showers occurring deep in the Earth's crust to be detectable as they exit the Earth into the atmosphere. A search for upward-going neutrino-induced showers in the HiRes-II monocular dataset has yielded a null result. From an LPM calculation of the energy spectrum of charged particles as a function of primary energy and depth for electron-induced showers in rock, we calculate the shape of the resulting profile of these showers in air. We describe a full detector Monte Carlo simulation to determine the detector response to upward-going electron-neutrino-induced cascades and present an upper limit on the flux of electron-neutrinos.
[18]  oai:arXiv.org:astro-ph/0010652  [pdf] - 1943537
Measurement of the Cosmic Ray Energy Spectrum and Composition from 10^{17} to 10^{18.3} eV Using a Hybrid Fluorescence Technique
Comments: 39 pages, 15 figures, in revtex4 epsf style, submited to APJ
Submitted: 2000-10-31
We study the spectrum and average mass composition of cosmic rays with primary energies between 10^{17} eV and 10^{18} eV using a hybrid detector consisting of the High Resolution Fly's Eye (HiRes) prototype and the MIA muon array. Measurements have been made of the change in the depth of shower maximum as a function of energy. A complete Monte Carlo simulation of the detector response and comparisons with shower simulations leads to the conclusion that the cosmic ray intensity is changing f rom a heavier to a lighter composition in this energy range. The spectrum is consistent with earlier Fly's Eye measurements and supports the previously found steepening near 4 \times 10^{17} eV .
[19]  oai:arXiv.org:astro-ph/9910195  [pdf] - 108749
Energy Estimation of UHE Cosmic Rays using the Atmospheric Fluorescence Technique
Comments: 16 pages, 4 figures, 1999 ICRC OG_4.5.09, accepted by Astroparticle Physics
Submitted: 1999-10-11, last modified: 1999-12-22
We use the CORSIKA air shower simulation program to review the method for assigning energies to ultra-high energy (UHE) cosmic rays viewed with the air fluorescence technique. This technique uses the atmosphere as a calorimeter, and we determine the corrections that must be made to the calorimetric energy to yield the primary cosmic ray energy.
[20]  oai:arXiv.org:astro-ph/9911144  [pdf] - 109272
A Multi-Component Measurement of the Cosmic Ray Composition Between 10^{17} eV and 10^{18} eV
Comments: 14 pages including 3 figures in revtex epsf style, submited to PRL
Submitted: 1999-11-09
The average mass composition of cosmic rays with primary energies between $10^{17}$eV and $10^{18}$eV has been studied using a hybrid detector consisting of the High Resolution Fly's Eye (HiRes) prototype and the MIA muon array. Measurements have been made of the change in the depth of shower maximum, $X_{max}$, and in the change in the muon density at a fixed core location, $\rho_\mu(600m)$, as a function of energy. The composition has also been evaluated in terms of the combination of $X_{max}$ and $\rho_\mu(600m)$. The results show that the composition is changing from a heavy to lighter mix as the energy increases.