Showing 1–28 of 28 results for author: Zhai, X

Accretion of the relativistic Vlasov gas onto a Kerr black hole

Authors: Ping Li, Yong-Qiang Liu, Xiang-Hua Zhai

Abstract: We study the accretion of relativistic Vlasov gas onto a Kerr black hole, regarding the particles as distributed throughout all the space, other than just in the equatorial plane. We solve the relativistic Liouville equation in the full $3+1$ dimensional framework of Kerr geometry. For the flow that is stationary and axial symmetric, we prove that the distribution function is independent of the co… ▽ More We study the accretion of relativistic Vlasov gas onto a Kerr black hole, regarding the particles as distributed throughout all the space, other than just in the equatorial plane. We solve the relativistic Liouville equation in the full $3+1$ dimensional framework of Kerr geometry. For the flow that is stationary and axial symmetric, we prove that the distribution function is independent of the conjugate coordinates. For an explicit distribution that can approximate to Maxwell-Jüttner distribution, we further calculate the particle current density, the stress energy momentum tensor and the unit accretion rates of mass, energy and angular momentum. The analytic results at large distance are shown to be consistent with the limits of the numerical ones computed at finite distance. Especially, we show that the unit mass accretion rate agrees with the Schwarzschild result in the case of low temperature limit. Furthermore, we find from the numerical results that the three unit accretion rates vary with the angle in Kerr metric and the accretion of Vlasov gas would slow down the Kerr black hole. The closer to the equator, the faster it slows down the black hole. △ Less

Submitted 18 November, 2023; originally announced November 2023.

Comments: 28 pages, 11 figures, accepted for publication in Physical Review D

Latest data constraint of some parameterized dark energy models

Authors: Jing Yang, Xin-Yan Fan, Chao-Jun Feng, Xiang-Hua Zhai

Abstract: Using various latest cosmological datasets including Type-Ia supernovae, cosmic microwave background radiation, baryon acoustic oscillations, and estimations of the Hubble parameter, we test some dark energy models with parameterized equations of state and try to distinguish or select observation-preferred models. We obtain the best fitting results of the six models and calculate their values of t… ▽ More Using various latest cosmological datasets including Type-Ia supernovae, cosmic microwave background radiation, baryon acoustic oscillations, and estimations of the Hubble parameter, we test some dark energy models with parameterized equations of state and try to distinguish or select observation-preferred models. We obtain the best fitting results of the six models and calculate their values of the Akaike Information Criteria and Bayes Information Criterion. And we can distinguish these dark energy models from each other by using these two information criterions. However, the $Λ$CDM model remains the best fit model. Furthermore, we perform geometric diagnostics including statefinder and Om diagnostics to understand the geometric behaviour of the dark energy models. We find that the six DE models can be distinguished from each other and from $Λ$CDM, Chaplygin gas, quintessence models after the statefinder and Om diagnostics were performed. Finally, we consider the growth factor of the dark energy models with comparison to $Λ$CDM model. Still, we find the models can be distinguished from each other and from $Λ$CDM model through the growth factor approximation. △ Less

Submitted 28 November, 2022; originally announced November 2022.

Comments: 23 pages, 6 figures, to be published in CPL

Flux Variations of Cosmic Ray Air Showers Detected by LHAASO-KM2A During a Thunderstorm on 10 June 2021

Authors: LHAASO Collaboration, F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Zhe Cao, Zhen Cao, J. Chang, J. F. Chang, E. S. Chen, Liang Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, S. H. Chen, S. Z. Chen, T. L. Chen, X. J. Chen , et al. (248 additional authors not shown)

Abstract: The Large High Altitude Air Shower Observatory (LHAASO) has three sub-arrays, KM2A, WCDA and WFCTA. The flux variations of cosmic ray air showers were studied by analyzing the KM2A data during the thunderstorm on 10 June 2021. The number of shower events that meet the trigger conditions increases significantly in atmospheric electric fields, with maximum fractional increase of 20%. The variations… ▽ More The Large High Altitude Air Shower Observatory (LHAASO) has three sub-arrays, KM2A, WCDA and WFCTA. The flux variations of cosmic ray air showers were studied by analyzing the KM2A data during the thunderstorm on 10 June 2021. The number of shower events that meet the trigger conditions increases significantly in atmospheric electric fields, with maximum fractional increase of 20%. The variations of trigger rates (increases or decreases) are found to be strongly dependent on the primary zenith angle. The flux of secondary particles increases significantly, following a similar trend with that of the shower events. To better understand the observed behavior, Monte Carlo simulations are performed with CORSIKA and G4KM2A (a code based on GEANT4). We find that the experimental data (in saturated negative fields) are in good agreement with simulations, assuming the presence of a uniform upward electric field of 700 V/cm with a thickness of 1500 m in the atmosphere above the observation level. Due to the acceleration/deceleration and deflection by the atmospheric electric field, the number of secondary particles with energy above the detector threshold is modified, resulting in the changes in shower detection rate. △ Less

Submitted 6 December, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: 18 pages, 11 figures

Journal ref: Chinese Phys. C 47 015001 (2023)

Peta-electron volt gamma-ray emission from the Crab Nebula

Authors: The LHAASO Collaboration, Zhen Cao, F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, H. Cai, J. T. Cai, Zhe Cao, J. Chang, J. F. Chang, B. M. Chen, E. S. Chen, J. Chen, Liang Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen , et al. (250 additional authors not shown)

Abstract: The Crab pulsar and the surrounding nebula powered by the pulsar's rotational energy through the formation and termination of a relativistic electron-positron wind is a bright source of gamma-rays carrying crucial information about this complex conglomerate. We report the detection of $γ$-rays with a spectrum showing gradual steepening over three energy decades, from $5\times 10^{-4}$ to $1.1$ pet… ▽ More The Crab pulsar and the surrounding nebula powered by the pulsar's rotational energy through the formation and termination of a relativistic electron-positron wind is a bright source of gamma-rays carrying crucial information about this complex conglomerate. We report the detection of $γ$-rays with a spectrum showing gradual steepening over three energy decades, from $5\times 10^{-4}$ to $1.1$ petaelectronvolt (PeV). The ultra-high-energy photons exhibit the presence of a PeV electron accelerator (a pevatron) with an acceleration rate exceeding 15% of the absolute theoretical limit. Assuming that unpulsed $γ$-rays are produced at the termination of the pulsar's wind, we constrain the pevatron's size, between $0.025$ and $0.1$ pc, and the magnetic field $\approx 110 μ$G. The production rate of PeV electrons, $2.5 \times 10^{36}$ erg $\rm s^{-1}$, constitutes 0.5% of the pulsar's spin-down luminosity, although we do not exclude a non-negligible contribution of PeV protons to the production of the highest energy $γ$-rays. △ Less

Submitted 11 November, 2021; originally announced November 2021.

Comments: 43 pages, 13 figures, 2 tables; Published in Science

Journal ref: Science, 2021, Vol 373, Issue 6553, pp. 425-430

Realistic neutron star models in $f(T)$ gravity

Authors: Rui-Hui Lin, Xiao-Ning Chen, Xiang-Hua Zhai

Abstract: We investigate the nonrotating neutron stars in $f(T)$ gravity with $f(T)=T+α{T}^2$, where $T$ is the torsion scalar in the teleparallel formalism of gravity. In particular, we utilize the SLy and BSk family of equations of state for perfect fluid to describe the neutron stellar matter and search for the effects of the $f(T)$ modification on the models of neutron stars. For positive $α$, the modif… ▽ More We investigate the nonrotating neutron stars in $f(T)$ gravity with $f(T)=T+α{T}^2$, where $T$ is the torsion scalar in the teleparallel formalism of gravity. In particular, we utilize the SLy and BSk family of equations of state for perfect fluid to describe the neutron stellar matter and search for the effects of the $f(T)$ modification on the models of neutron stars. For positive $α$, the modification results in a smaller stellar mass in comparison to general relativity, while the neutron stars will contain larger amount of matter for negative $α$. Moreover, there seems to be an upper limit for the central density of the neutron stars with $α>0$, beyond which the effective $f(T)$ fluid would have a steplike phase transition in density and pressure profiles, collapsing the numerical system. We obtain the mass-radius relations of the realistic models of neutron stars and subject them to the joint constraints from the observed massive pulsars PSR J0030+0451, PSR J0740+6620, and PSR J2215+5135, and gravitational wave events GW170817 and GW190814. For the neutron star model in $f(T)$ gravity to be able to accommodate all the mentioned data, the model parameter $α$ needs to be smaller than $-4.295$, $-6.476$, $-4.4$, and $-2.12$ (in the unit of ${G}^2M_\odot^2/c^4$) for SLy, BSk19, BSk20, and BSk21 equations of state, respectively. If one considers the unknown compact object in the event GW190814 not to be a neutron star and hence excludes this dataset, the constraints can be loosened to $α<-0.594$, $-3.5$, $0.4$ and $1.9$ (in the unit of ${G}^2M_\odot^2/c^4$), respectively. △ Less

Submitted 28 March, 2022; v1 submitted 1 September, 2021; originally announced September 2021.

Comments: version accepted for publication in EPJC

Superradiant instability of the Kerr-like black hole in Einstein-bumblebee gravity

Authors: Rui Jiang, Rui-Hui Lin, Xiang-Hua Zhai

Abstract: An exact Kerr-like solution has been obtained recently in Einstein-bumblebee gravity model where Lorentz symmetry is spontaneously broken. In this paper, we investigate the superradiant instability of the Kerr-like black hole under the perturbation of a massive scalar field. We find the Lorentz breaking parameter $L$ does not affect the superradiance regime or the regime of the bound states. Howev… ▽ More An exact Kerr-like solution has been obtained recently in Einstein-bumblebee gravity model where Lorentz symmetry is spontaneously broken. In this paper, we investigate the superradiant instability of the Kerr-like black hole under the perturbation of a massive scalar field. We find the Lorentz breaking parameter $L$ does not affect the superradiance regime or the regime of the bound states. However, since $L$ appears in the metric and its effect cannot be erased by redefining the rotation parameter $\tilde{a}=\sqrt{1+L}a$, it indeed affects the bound state spectrum and the superradiance. We calculate the bound state spectrum via the continued-fraction method and show the influence of $L$ on the maximum binding energy and the damping rate. The superradiant instability could occur since the superradiance condition and the bound state condition could be both satisfied. Compared with Kerr black hole, the nature of the superradiant instability of this black hole depends non-monotonously not only on the rotation parameter of the black hole $\tilde{a}$ and the product of the black hole mass $M$ and the field mass $μ$, but also on the Lorentz breaking parameter $L$. Through the Monte Carlo method, we find that for $l=m=1$ state the most unstable mode occurs at $L=-0.79637$, $\tilde{a}/M=0.99884$ and $Mμ=0.43920$, with the maximum growth rate of the field $ω_{I}M=1.676\times10^{-6}$, which is about 10 times of that in Kerr black hole. △ Less

Submitted 5 November, 2021; v1 submitted 10 August, 2021; originally announced August 2021.

Comments: 18 pages, 3 figures, 1 table, version accepted for publication in Physical Review D

Calibration of the Air Shower Energy Scale of the Water and Air Cherenkov Techniques in the LHAASO experiment

Authors: F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, H. Cai, J. T. Cai, Z. Cao Z. Cao, J. Chang, J. F. Chang, X. C. Chang, B. M. Chen, J. Chen, L. Chen, L. Chen, L. Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (233 additional authors not shown)

Abstract: The Wide Field-of-View Cherenkov Telescope Array (WFCTA) and the Water Cherenkov Detector Arrays (WCDA) of LHAASO are designed to work in combination for measuring the energy spectra of various cosmic ray species over a very wide energy range from a few TeV to 10 PeV. The energy calibration of WCDA can be achieved with a proven technique of measuring the westward shift of the Moon shadow of galact… ▽ More The Wide Field-of-View Cherenkov Telescope Array (WFCTA) and the Water Cherenkov Detector Arrays (WCDA) of LHAASO are designed to work in combination for measuring the energy spectra of various cosmic ray species over a very wide energy range from a few TeV to 10 PeV. The energy calibration of WCDA can be achieved with a proven technique of measuring the westward shift of the Moon shadow of galactic cosmic rays due to the geomagnetic field. This deflection angle $Δ$ is inversely proportional to the energy of the cosmic rays. The precise measurements of the shifts by WCDA allows us to calibrate its energy scale for energies as high as 35 TeV. The energy scale measured by WCDA can be used to cross calibrate the energy reconstructed by WFCTA, which spans the whole energy range up to 10 PeV. In this work, we will demonstrate the feasibility of the method using the data collected from April 2019 to January 2020 by the WFCTA array and WCDA-1 detector, the first of the three water Cherenkov ponds, already commissioned at LHAASO site. △ Less

Submitted 13 April, 2021; v1 submitted 11 April, 2021; originally announced April 2021.

Journal ref: Phys. Rev. D 104, 062007 (2021)

Construction and On-site Performance of the LHAASO WFCTA Camera

Authors: F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, H. Cai, J. T. Cai, Z. Cao, Z. Cao, J. Chang, J. F. Chang, X. C. Chang, B. M. Chen, J. Chen, L. Chen, L. Chen, L. Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen , et al. (234 additional authors not shown)

Abstract: The focal plane camera is the core component of the Wide Field-of-view Cherenkov/fluorescence Telescope Array (WFCTA) of the Large High-Altitude Air Shower Observatory (LHAASO). Because of the capability of working under moonlight without aging, silicon photomultipliers (SiPM) have been proven to be not only an alternative but also an improvement to conventional photomultiplier tubes (PMT) in this… ▽ More The focal plane camera is the core component of the Wide Field-of-view Cherenkov/fluorescence Telescope Array (WFCTA) of the Large High-Altitude Air Shower Observatory (LHAASO). Because of the capability of working under moonlight without aging, silicon photomultipliers (SiPM) have been proven to be not only an alternative but also an improvement to conventional photomultiplier tubes (PMT) in this application. Eighteen SiPM-based cameras with square light funnels have been built for WFCTA. The telescopes have collected more than 100 million cosmic ray events and preliminary results indicate that these cameras are capable of working under moonlight. The characteristics of the light funnels and SiPMs pose challenges (e.g. dynamic range, dark count rate, assembly techniques). In this paper, we present the design features, manufacturing techniques and performances of these cameras. Finally, the test facilities, the test methods and results of SiPMs in the cameras are reported here. △ Less

Submitted 4 July, 2021; v1 submitted 29 December, 2020; originally announced December 2020.

Comments: 45 pages, 21 figures, article

Journal ref: Eur. Phys. J. C 81, 657 (2021)

The observation of the Crab Nebula with LHAASO-KM2A for the performance study

Authors: F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, H. Cai, J. T. Cai, Z. Cao, Z. Cao, J. Chang, J. F. Chang, X. C. Chang, B. M. Chen, J. Chen, L. Chen, L. Chen, L. Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen , et al. (234 additional authors not shown)

Abstract: As a sub-array of the Large High Altitude Air Shower Observatory (LHAASO), KM2A is mainly designed to cover a large fraction of the northern sky to hunt for gamma-ray sources at energies above 10 TeV. Even though the detector construction is still underway, a half of the KM2A array has been operating stably since the end of 2019. In this paper, we present the pipeline of KM2A data analysis and the… ▽ More As a sub-array of the Large High Altitude Air Shower Observatory (LHAASO), KM2A is mainly designed to cover a large fraction of the northern sky to hunt for gamma-ray sources at energies above 10 TeV. Even though the detector construction is still underway, a half of the KM2A array has been operating stably since the end of 2019. In this paper, we present the pipeline of KM2A data analysis and the first observation on the Crab Nebula, a standard candle in very high energy gamma-ray astronomy. We detect gamma-ray signals from the Crab Nebula in both energy ranges of 10$-$100 TeV and $>$100 TeV with high significance, by analyzing the KM2A data of 136 live days between December 2019 and May 2020. With the observations, we test the detector performance including angular resolution, pointing accuracy and cosmic ray background rejection power. The energy spectrum of the Crab Nebula in the energy range 10-250 TeV fits well with a single power-law function dN/dE =(1.13$\pm$0.05$_{stat}$$\pm$0.08$_{sys}$)$\times$10$^{-14}$$\cdot$(E/20TeV)$^{-3.09\pm0.06_{stat}\pm0.02_{sys}}$ cm$^{-2}$ s$^{-1}$ TeV$^{-1}$. It is consistent with previous measurements by other experiments. This opens a new window of gamma-ray astronomy above 0.1 PeV through which ultrahigh-energy gamma-ray new phenomena, such as cosmic PeVatrons, might be discovered. △ Less

Submitted 13 October, 2020; originally announced October 2020.

Comments: 13 pages, 15 figures,submitted to CPC

Journal ref: Chinese Phys. C 45 025002 (2021)

Artificial Neural Network Spectral Light Curve Template for Type Ia Supernovae and its Cosmological Constraints

Authors: Qiao-Bin Cheng, Chao-Jun Feng, Xiang-Hua Zhai, Xin-Zhou Li

Abstract: The spectral energy distribution (SED) sequence for type Ia supernovae (SN Ia) is modeled by an artificial neural network. The SN Ia luminosity is characterized as a function of phase, wavelength, a color parameter and a decline rate parameter. After training and testing the neural network, the SED sequence could give both the spectrum with wavelength range from 3000Å~to 8000Å~ and the light curve… ▽ More The spectral energy distribution (SED) sequence for type Ia supernovae (SN Ia) is modeled by an artificial neural network. The SN Ia luminosity is characterized as a function of phase, wavelength, a color parameter and a decline rate parameter. After training and testing the neural network, the SED sequence could give both the spectrum with wavelength range from 3000Å~to 8000Å~ and the light curve with phase from 20 days before to 50 days after the maximum luminosity for the supernovae with different colors and decline rates. Therefore, we call this the Artificial Neural Network Spectral Light Curve Template (ANNSLCT) model. We retrain the Joint Light-curve Analysis (JLA) supernova sample by using the ANNSLCT model and obtain the parameters for each supernova to make a constraint on the cosmological $Λ$CDM model. We find that the best fitting values of these parameters are almost the same as those from the JLA sample trained with the Spectral Adaptive Lightcurve Template 2 (SALT2) model. So we believe that the ANNSLCT model could be used to analyze a large number of SN Ia multi-color light curves measured in the current and future observational projects. △ Less

Submitted 9 April, 2020; originally announced April 2020.

Comments: 8 pages two columns, 17 figures

Dynamics of Viscous Phantom Universe

Authors: Ji-Yao Wang, Chao-Jun Feng, Xiang-Hua Zhai, Xin-Zhou Li

Abstract: The phantom dark energy remarkably boosts our prehension of the accelerating Universe. Various models are widely discussed in the phantom Universe without bulk viscosity. From the hydrodynamics' point of view, it is natural to introduce the nonperfect fluid in the study of the Universe, as an ideal fluid just an approximation to the real world after all and using the generalized equation of state… ▽ More The phantom dark energy remarkably boosts our prehension of the accelerating Universe. Various models are widely discussed in the phantom Universe without bulk viscosity. From the hydrodynamics' point of view, it is natural to introduce the nonperfect fluid in the study of the Universe, as an ideal fluid just an approximation to the real world after all and using the generalized equation of state (EoS) with bulk viscosity, the early inflationary universe and the accelerated expansion of the late-time universe are described by many authors. In this paper, in order to investigate how the viscosity will influence the evolution of the Universe, we study a class of phantom dark energy models with bulk viscosity by the method of dynamical analysis technique. We show that there are different cosmic late-time behaviors and the stability also brings some constraints on the models. We also plot the evolutionary trajectories of this model in the statefinder parameter-planes to see the different behaviors of the models from the statefinder viewpoint. △ Less

Submitted 9 June, 2020; v1 submitted 12 March, 2020; originally announced March 2020.

Comments: 14 pages, 7 figures

Solar System Tests of a New Class of $f(z)$ Theory

Authors: Ji-Yao Wang, Chao-Jun Feng, Xiang-Hua Zhai, Xin-Zhou Li

Abstract: Recently, a new kind of $f(z)$ theory is proposed to provide a different perspective for the development of reliable alternative models of gravity in which the $f(R)$ Lagrangian terms are reformulated as polynomial parameterizations $f(z)$. In the previous study, the parameters in the $f(z)$ models have been constrained by using cosmological data. In this paper, these models will be tested by the… ▽ More Recently, a new kind of $f(z)$ theory is proposed to provide a different perspective for the development of reliable alternative models of gravity in which the $f(R)$ Lagrangian terms are reformulated as polynomial parameterizations $f(z)$. In the previous study, the parameters in the $f(z)$ models have been constrained by using cosmological data. In this paper, these models will be tested by the observations in the solar system. After solving the Ricci scalar as a function of the redshift, one could obtain $f(R)$ that could be used to calculate the standard Parameterized-Post-Newtonian (PPN) parameters. First, we fit the parametric models with the latest cosmological observational data. Then the tests are performed by solar system observations. And last we combine the constraints of solar system and cosmology together and reconstruct the $f(R)$ actions of the $f(z)$ parametric models. △ Less

Submitted 7 May, 2020; v1 submitted 29 May, 2019; originally announced May 2019.

Comments: 15 pages

Journal ref: International Journal of Modern Physics D. Vol. 29, No. 08, 2050060 (2020)

Generalized Birkhoff theorem and its applications in mimetic gravity

Authors: Xin-zhou Li, Xiang-hua Zhai, Ping Li

Abstract: There is undetermined potential function $V(φ)$ in the action of mimetic gravity which should be resolved through physical means. In general relativity(GR), the static spherically symmetric(SSS) solution to the Einstein equation is a benchmark and its deformation also plays a crucial role in mimetic gravity. The equation of motion is provided with high nonlinearity, but we can reduce primal nonlin… ▽ More There is undetermined potential function $V(φ)$ in the action of mimetic gravity which should be resolved through physical means. In general relativity(GR), the static spherically symmetric(SSS) solution to the Einstein equation is a benchmark and its deformation also plays a crucial role in mimetic gravity. The equation of motion is provided with high nonlinearity, but we can reduce primal nonlinearity to a frequent Riccati form in the SSS case of mimetic gravity. In other words, we obtain an expression of solution to the functional differential equation of motion with any potential function. Remarkably, we proved rigorously that there is a zero point of first order for the metric function $β(r)$ if another metric function $α(r)$ possesses a pole of first order within mimetic gravity. The zero point theorem may be regarded as the generalization of Birkhoff theorem $αβ=1$ in GR. As a corollary, we show that there is a modified black hole solution for any given $V(φ)$, which can pass the test of solar system. As another corollary, the zero point theorem provides a dynamical mechanism for the maximum size of galaxies. Especially, there are two analytic solutions which provide good fits to the rotation curves of galaxies without the demand for particle dark matter. △ Less

Submitted 22 July, 2018; originally announced July 2018.

Comments: 11 pages

Diagnostics for generalized power-law torsion-matter coupling $f(T)$ model

Authors: Xiang-Hua Zhai, Qiang Wen, Rui-Hui Lin, Xin-Zhou Li

Abstract: The currently accelerated expansion of our Universe is unarguably one of the most intriguing problems in today's physics research. Two realistic non-minimal torsion-matter coupling $f(T)$ models have been established and studied in our previous papers [Phys. Rev. D92, 104038(2015) and Eur. Phys. J. C77, 504(2017)] aiming to explain this "dark energy" problem. In this paper, we study the generalize… ▽ More The currently accelerated expansion of our Universe is unarguably one of the most intriguing problems in today's physics research. Two realistic non-minimal torsion-matter coupling $f(T)$ models have been established and studied in our previous papers [Phys. Rev. D92, 104038(2015) and Eur. Phys. J. C77, 504(2017)] aiming to explain this "dark energy" problem. In this paper, we study the generalized power-law torsion-matter coupling $f(T)$ model. Dynamical system analysis shows that the three expansion phases of the Universe, i.e. the radiation dominated era, the matter dominated era and the dark energy dominated era, can all be reproduced in this generalized model. By using the statefinder and $Om$ diagnostics, we find that the different cases of the model can be distinguished from each other and from other dark energy models such as the two models in our previous papers, $Λ$CDM, quintessence and Chaplygin gas. Furthermore, the analyses also show that all kinds of generalized power-law torsion-matter coupling model are able to cross the $w=-1$ divide from below to above, thus the decrease of the energy density resulting from the crossing of $w$ will make the catastrophic fate of the Universe avoided and a de Sitter expansion fate in the future will be approached. △ Less

Submitted 12 March, 2018; originally announced March 2018.

Comments: 9 pages,12 figures

Artificial Neural Network for Constructing Type Ia Supernovae Spectrum Evolution Model

Authors: Qiao-Bin Cheng, Chao-Jun Feng, Xiang-Hua Zhai, Xin-Zhou Li

Abstract: We construct and train an artificial neural network called the back-propagation neural network to describe the evolution of the type Ia supernova spectrum by using the data from the CfA Supernova Program. This network method has many attractive features, and one of them is that the constructed model is differentiable. Benefitting from this, we calculate the absorption velocity and its variation. T… ▽ More We construct and train an artificial neural network called the back-propagation neural network to describe the evolution of the type Ia supernova spectrum by using the data from the CfA Supernova Program. This network method has many attractive features, and one of them is that the constructed model is differentiable. Benefitting from this, we calculate the absorption velocity and its variation. The model we constructed can well describe not only the spectrum of SNe Ia with wavelength range from $3500Å$ to $8000Å$, but also the light-curve evolution with phase time from $-15$ to $50$ with different colors. Moreover, the number of parameters needed during the training process is much less than the usual methods. △ Less

Submitted 14 January, 2018; v1 submitted 5 January, 2018; originally announced January 2018.

Comments: v2 refs added, two columns, 7pages, 12 figures

Journal ref: Phys. Rev. D 97, 123530 (2018)

Action functional of the Cardassian universe

Authors: Xiang-hua Zhai, Rui-hui Lin, Chao-jun Feng, Xin-zhou Li

Abstract: It is known that the Cardassian universe is successful in describing the accelerated expansion of the universe, but its dynamical equations are hard to get from the action principle. In this paper, we establish the connection between the Cardassian universe and $f(T, \mathcal{T})$ gravity, where $T$ is the torsion scalar and $\mathcal{T}$ is the trace of the matter energy-momentum tensor. For dust… ▽ More It is known that the Cardassian universe is successful in describing the accelerated expansion of the universe, but its dynamical equations are hard to get from the action principle. In this paper, we establish the connection between the Cardassian universe and $f(T, \mathcal{T})$ gravity, where $T$ is the torsion scalar and $\mathcal{T}$ is the trace of the matter energy-momentum tensor. For dust matter, we find that the modified Friedmann equations from $f(T, \mathcal{T})$ gravity can correspond to those of Cardassian models, and thus, a possible origin of Cardassian universe is given. We obtain the original Cardassian model, the modified polytropic Cardassian model, and the exponential Cardassian model from the Lagrangians of $f(T,\mathcal{T})$ theory. Furthermore, by adding an additional term to the corresponding Lagrangians, we give three generalized Cardassian models from $f(T,\mathcal{T})$ theory. Using the observation data of type Ia supernovae, cosmic microwave background radiation, and baryon acoustic oscillations, we get the fitting results of the cosmological parameters and give constraints of model parameters for all of these models. △ Less

Submitted 26 May, 2017; originally announced May 2017.

Comments: 7 pages, 3 figures

Journal ref: Physical Review D 95, 104030 (2017)

Spherical collapse and virialization in $f(T)$ gravities

Authors: Rui-hui Lin, Xiang-hua Zhai, Xin-zhou Li

Abstract: Using the classical top-hat profile, we study the non-linear growth of spherically symmetric density perturbation and structure formation in $f(T)$ gravities. In particular, three concrete models, which have been tested against the observation of large-scale evolution and linear perturbation of the universe in the cosmological scenario, are investigated in this framework, covering both minimal and… ▽ More Using the classical top-hat profile, we study the non-linear growth of spherically symmetric density perturbation and structure formation in $f(T)$ gravities. In particular, three concrete models, which have been tested against the observation of large-scale evolution and linear perturbation of the universe in the cosmological scenario, are investigated in this framework, covering both minimal and nonminimal coupling cases of $f(T)$ gravities. Moreover, we consider the virialization of the overdense region in the models after they detach from the background expanding universe and turn around to collapse. We find that there are constraints in the magnitude and occurring epoch of the initial perturbation. The existence of these constraints indicates that a perturbation that is too weak or occurs too late will not be able to stop the expanding of the overdense region. The illustration of the evolution of the perturbation shows that in $f(T)$ gravities, the initial perturbation within the constraints can eventually lead to clustering and form structure. The evolution also shows that nonminimal coupling models collapse slower than the minimal coupling one. △ Less

Submitted 23 March, 2017; v1 submitted 18 December, 2016; originally announced December 2016.

Comments: 19 pages,5 figures

Journal ref: JCAP03(2017)040

Vaidya solution and its generalization in de Rham-Gabadadze-Tolley massive gravity

Authors: Ping Li, Xin-zhou Li, Xiang-hua Zhai

Abstract: We present a detailed study of the Vaidya solution and its generalization in de Rham-Gabadadze-Tolley (dRGT) theory. Since the diffeomorphism invariance can be restored with the Stückelberg fields $φ^a$ introduced, there is a new invariant $I^{ab}=g^{μν}\partial_μφ^a\partial_νφ^b$ in the massive gravity, which adds to the ones usually encountered in general relativity. There is no conventional Vai… ▽ More We present a detailed study of the Vaidya solution and its generalization in de Rham-Gabadadze-Tolley (dRGT) theory. Since the diffeomorphism invariance can be restored with the Stückelberg fields $φ^a$ introduced, there is a new invariant $I^{ab}=g^{μν}\partial_μφ^a\partial_νφ^b$ in the massive gravity, which adds to the ones usually encountered in general relativity. There is no conventional Vaidya solution if we choose unitary gauge. In this paper, we obtain three types of self-consistent ansatz with some nonunitary gauge, and find accordingly the Vaidya, generalized Vaidya and furry Vaidya solution. As by-products, we obtain a series of furry black hole. The Vaidya solution and its generalization in dRGT massive gravity describe the black holes with a variable horizon. △ Less

Submitted 1 December, 2016; originally announced December 2016.

Comments: 9 pages, to be published in Phys. Rev. D. arXiv admin note: substantial text overlap with arXiv:1603.06039

Journal ref: Phys. Rev. D 94, 124022 (2016)

Three-Dimensional MHD Simulation of Caltech Plasma Jet Experiment: First Results

Authors: Xiang Zhai, Hui Li, Paul M. Bellan, Shengtai Li

Abstract: Magnetic fields are believed to play an essential role in astrophysical jets with observations suggesting the presence of helical magnetic fields. Here, we present three-dimensional (3D) ideal MHD simulationsof the Caltech plasma jet experiment using a magnetic tower scenario as the baseline model. Magnetic fields consist of an initially localized dipole-like poloidal component and a toroidal comp… ▽ More Magnetic fields are believed to play an essential role in astrophysical jets with observations suggesting the presence of helical magnetic fields. Here, we present three-dimensional (3D) ideal MHD simulationsof the Caltech plasma jet experiment using a magnetic tower scenario as the baseline model. Magnetic fields consist of an initially localized dipole-like poloidal component and a toroidal component that is continuously being injected into the domain. This flux injection mimics the poloidal currents driven by the anode-cathode voltage drop in the experiment. The injected toroidal field stretches the poloidal fields to large distances, while forming a collimated jet along with several other key features. Detailed comparisons between 3D MHD simulations and experimental measurements provide a comprehensive description of the interplay among magnetic force, pressure and flow effects. In particular, we delineate both the jet structure and the transition process that converts the injected magnetic energy to other forms. With suitably chosen parameters that are derived from experiments, the jet in the simulation agrees quantitatively with the experimental jet in terms of magnetic/kinetic/inertial energy, total poloidal current, voltage, jet radius, and jet propagation velocity. Specifically, the jet velocity in the simulation is proportional to the poloidal current divided by the square root of the jet density, in agreement with both the experiment and analytical theory. This work provides a new and quantitative method for relating experiments, numerical simulations and astrophysical observation, and demonstrates the possibility of using terrestrial laboratory experiments to study astrophysical jets. △ Less

Submitted 13 July, 2014; originally announced July 2014.

Comments: accepted by ApJ 37 pages, 15 figures, 2 tables

Journal ref: Zhai et al. 2014, Astrophysical Journal, 791, 40

Schwarzschild solution as a result of thermodynamics

Authors: Hongsheng Zhang, Sean A. Hayward, Xiang-Hua Zhai, Xin-Zhou Li

Abstract: We obtain the Schwarzschild solution from thermodynamic considerations using the assumptions of a quasi local mass form (the Misner-Sharp mass) and geometric surface gravity in a spherically symmetric spacetime. The deduction is extended to other cases such as the de Sitter, anti-de Sitter, Reissner-Nordström and higher dimensional spacetimes. This paper demonstrates the simple hypotheses to obtai… ▽ More We obtain the Schwarzschild solution from thermodynamic considerations using the assumptions of a quasi local mass form (the Misner-Sharp mass) and geometric surface gravity in a spherically symmetric spacetime. The deduction is extended to other cases such as the de Sitter, anti-de Sitter, Reissner-Nordström and higher dimensional spacetimes. This paper demonstrates the simple hypotheses to obtain these known spherically symmetric solutions via thermodynamics, where essentially the Misner-Sharp mass is the mass for an adiabatic system. △ Less

Submitted 2 March, 2014; v1 submitted 8 April, 2013; originally announced April 2013.

Comments: two new authors joined, 5pages, no figure, PRD in press

Journal ref: Phys.Rev. D89, 064052 (2014)

Alternative mechanism of avoiding the big rip or little rip for a scalar phantom field

Authors: Ping Xi, Xiang-hua Zhai, Xin-zhou Li

Abstract: Depending on the choice of its potential, the scalar phantom field $φ$ (the equation of state parameter $w<-1$) leads to various catastrophic fates of the universe including big rip, little rip and other future singularity. For example, big rip results from the evolution of the phantom field with an exponential potential and little rip stems from a quadratic potential in general relativity (GR). B… ▽ More Depending on the choice of its potential, the scalar phantom field $φ$ (the equation of state parameter $w<-1$) leads to various catastrophic fates of the universe including big rip, little rip and other future singularity. For example, big rip results from the evolution of the phantom field with an exponential potential and little rip stems from a quadratic potential in general relativity (GR). By choosing the same potential as in GR, we suggest a new mechanism to avoid these unexpected fates (big and little rip) in the inverse-\textit{R} gravity. As a pedagogical illustration, we give an exact solution where phantom field leads to a power-law evolution of the scale factor in an exponential type potential. We also find the sufficient condition for a universe in which the equation of state parameter crosses $w=-1$ divide. The phantom field with different potentials, including quadratic, cubic, quantic, exponential and logarithmic potentials are studied via numerical calculation in the inverse-\textit{R} gravity with $R^{2}$ correction. The singularity is avoidable under all these potentials. Hence, we conclude that the avoidance of big or little rip is hardly dependent on special potential. △ Less

Submitted 28 November, 2011; originally announced November 2011.

Comments: 9 pages,6 figures

Diffuse X-ray Emissions from Dynamic Planetary Nebulae

Authors: Yu-Qing Lou, Xiang Zhai

Abstract: We present theoretical results of a piecewise isothermal shock wind model devised for predicting the luminosity and surface brightness profile of diffuse X-ray emissions primarily from the inner shocked downstream wind zone of a planetary nebula (PN) surrounded by self-similar shocked dense shell and outer slow AGB wind envelope involving self-gravity and compare/fit our computational model result… ▽ More We present theoretical results of a piecewise isothermal shock wind model devised for predicting the luminosity and surface brightness profile of diffuse X-ray emissions primarily from the inner shocked downstream wind zone of a planetary nebula (PN) surrounded by self-similar shocked dense shell and outer slow AGB wind envelope involving self-gravity and compare/fit our computational model results with available observations of a few grossly spherical X-ray emitting PNe. Matching shocked piecewise isothermal self-similar void (ISSV) solutions with self-gravity of Lou & Zhai (LZ) for the outer zone and a stationary isothermal fast tenuous wind with a reverse shock for the inner zone across an expanding contact discontinuity, we can consistently construct dynamic evolution models of PNe with diffuse X-ray emissions. On the basis of such a chosen dynamic wind interaction model, both X-ray luminosity and radial X-ray brightness profile are determined by three key parameters, namely the so-called X-ray parameter $X$, two radii $R_{rs}$ and $R_c$ of the reverse shock and the contact discontinuity. We find that morphologies of X-ray emissions would appear in the forms of either a central luminous sphere or a bright ring embedded within optically bright shells. In contrast to previous adiabatic models, the X-ray brightness peaks around the reverse shock, instead of the contact discontinuity surface just inside the outer shocked dense shell. Diffuse X-ray emissions of a few observed PNe appear to support this wind-wind dynamic interaction scenario. △ Less

Submitted 14 June, 2010; originally announced June 2010.

Comments: 13 pages, 8 figures, accepted by Mon. Not. R. Astron. Soc. Accepted 2010 June 2. Received 2010 May 28; in original form 2010 May 4

Dynamic Evolution Model of Isothermal Voids and Shocks

Authors: Yu-Qing Lou, Xiang Zhai

Abstract: We explore self-similar hydrodynamic evolution of central voids embedded in an isothermal gas of spherical symmetry under the self-gravity. More specifically, we study voids expanding at constant radial speeds in an isothermal gas and construct all types of possible void solutions without or with shocks in surrounding envelopes. We examine properties of void boundaries and outer envelopes. Voids… ▽ More We explore self-similar hydrodynamic evolution of central voids embedded in an isothermal gas of spherical symmetry under the self-gravity. More specifically, we study voids expanding at constant radial speeds in an isothermal gas and construct all types of possible void solutions without or with shocks in surrounding envelopes. We examine properties of void boundaries and outer envelopes. Voids without shocks are all bounded by overdense shells and either inflows or outflows in the outer envelope may occur. These solutions, referred to as type $\mathcal{X}$ void solutions, are further divided into subtypes $\mathcal{X}_{\rm I}$ and $\mathcal{X}_{\rm II}$ according to their characteristic behaviours across the sonic critical line (SCL). Void solutions with shocks in envelopes are referred to as type $\mathcal{Z}$ voids and can have both dense and quasi-smooth edges. Asymptotically, outflows, breezes, inflows, accretions and static outer envelopes may all surround such type $\mathcal{Z}$ voids. Both cases of constant and varying temperatures across isothermal shock fronts are analyzed; they are referred to as types $\mathcal{Z}_{\rm I}$ and $\mathcal{Z}_{\rm II}$ void shock solutions. We apply the `phase net matching procedure' to construct various self-similar void solutions. We also present analysis on void generation mechanisms and describe several astrophysical applications. By including self-gravity, gas pressure and shocks, our isothermal self-similar void (ISSV) model is adaptable to various astrophysical systems such as planetary nebulae, hot bubbles and superbubbles in the interstellar medium as well as supernova remnants. △ Less

Submitted 30 May, 2009; originally announced June 2009.

Comments: 24 pages, 13 figuers, accepted by ApSS

Journal ref: Astrophys.Space Sci.323:17-40,2009

Global monopole surrounded by quintessence-like matter

Authors: Xin-zhou Li, Ping Xi, Xiang-hua Zhai

Abstract: We present new static spherically-symmetric solutions of Einstein equations with the quintessence-like matter surrounding a global monopole. These new solutions of the coupling scalar-Einstein equations are more complicated, which depend on the parameter of equation of state $-1 < w_{q} <-{1/3}$. A gravitating global monopole produces a gravitational field of de Sitter kind outside the core in a… ▽ More We present new static spherically-symmetric solutions of Einstein equations with the quintessence-like matter surrounding a global monopole. These new solutions of the coupling scalar-Einstein equations are more complicated, which depend on the parameter of equation of state $-1 < w_{q} <-{1/3}$. A gravitating global monopole produces a gravitational field of de Sitter kind outside the core in addition to a solid angular deficit. In the $w_{q} = -{1/3}$ case, we have proved that the solution cannot exist since the density of quintessence-like tends to zero if $w_{q} \to -{1/3}$. As a new feature, these monopoles have the outer horizon depending on both Goldstone field and quintessence-like. Since current observations constrain $-1.14 < w_{q} < -0.93$, new global monopoles have interesting astrophysical applications. △ Less

Submitted 1 July, 2008; v1 submitted 2 December, 2007; originally announced December 2007.

Comments: 8 pages, 8 figures, added discussion and some references, the form accepted for publication in Physics Letter B

Journal ref: Phys.Lett.B666:125-130,2008

Viscous generalized Chaplygin gas

Authors: Xiang-hua Zhai, You-dong Xu, Xin-zhou Li

Abstract: Viscous GCG(generalized Chaplygin gas) cosmology is discussed, assuming that there is bulk viscosity in the barotropic fluid and GCG. The dynamical analysis indicates that the phase $w_g=-1+\sqrt{3}γκτ_{g}/(γ-\sqrt{3}κτ_γ)$ is a dynamical attractor and the equation of state of GCG approaches it from either $w_g>-1$ or $w_g<-1$ depending on the choice of its initial cosmic density parameter and t… ▽ More Viscous GCG(generalized Chaplygin gas) cosmology is discussed, assuming that there is bulk viscosity in the barotropic fluid and GCG. The dynamical analysis indicates that the phase $w_g=-1+\sqrt{3}γκτ_{g}/(γ-\sqrt{3}κτ_γ)$ is a dynamical attractor and the equation of state of GCG approaches it from either $w_g>-1$ or $w_g<-1$ depending on the choice of its initial cosmic density parameter and the ratio of pressure to critical energy density. Obviously, the equation of state $w_g$ can cross the boundary $w_g=-1$. Also, from the point of view of dynamics, the parameters of viscous GCG should be in the range of $γ>\sqrt{3}κτ_γ/(1-\sqrt{3}κτ_{g})$ and $0<α<1+\sqrt{3} κγτ_{g}/(γ-\sqrt{3}κτ_γ-\sqrt{3}κγ\t au_{g})$. △ Less

Submitted 23 January, 2006; v1 submitted 29 November, 2005; originally announced November 2005.

Comments: 8 pages, 5 eps figures, typos corrected

Journal ref: Int.J.Mod.Phys.D15:1151-1162,2006

Dynamics of Quintessential Inflation

Authors: Xiang-hua Zhai, Yi-bin Zhao

Abstract: In this paper, we study a realistic model of quintessential inflation with radiation and matter. By the analysis of the dynamical system and numerical work about the evolution of the equation of state and cosmic density parameter, we show that this model is a good match for the current astronomical observation. The conclusion we obtain is in favor of the model where the modular part of complex f… ▽ More In this paper, we study a realistic model of quintessential inflation with radiation and matter. By the analysis of the dynamical system and numerical work about the evolution of the equation of state and cosmic density parameter, we show that this model is a good match for the current astronomical observation. The conclusion we obtain is in favor of the model where the modular part of complex field plays the role of the inflaton whereas the argument part is the quintessence field. The numerical calculation shows that a heteroclinic orbit (solution of dynamical system) which interpolates between early-time de Sitter phase (an unstable critical point) and a late-time de Sitter attractor. △ Less

Submitted 29 November, 2005; v1 submitted 16 November, 2005; originally announced November 2005.

Comments: 5 pages, 2 figures

Journal ref: Chin.Phys.15:2465,2006

On the Quintessence with Abelian and Non-abelian Symmetry

Authors: Xin-zhou Li, Jian-gang Hao, Dao-jun Liu, Xiang-hua Zhai

Abstract: We study the perturbations on both "radial" and "angular" components of the quintessence with an internal abelian and non-abelian symmetry. The properties of the perturbation on the "radial" component depend on the specific potential of the model and is similiar for both abelian and non-abelian case. We show that the consine-type potential is very interesting for the O(\textit{N}) quintessence m… ▽ More We study the perturbations on both "radial" and "angular" components of the quintessence with an internal abelian and non-abelian symmetry. The properties of the perturbation on the "radial" component depend on the specific potential of the model and is similiar for both abelian and non-abelian case. We show that the consine-type potential is very interesting for the O(\textit{N}) quintessence model and also give a critical condition of instability for the potential. While the properties of perturbations on "angular" components depend on whether the internal symmetry is abelian or non-abelian, which we have discussed respectively. In the non-abelian case, the fluctuation of the "angular" component will increase rapidly with time while in the abelian case it will not. △ Less

Submitted 18 August, 2003; originally announced August 2003.

Comments: 15 pages, to be published in Int. J. Mod. Phys. A

Journal ref: Int.J.Mod.Phys.A18:5921,2003; Int.J.Mod.Phys.A18:5921-5930,2003

The Tachyon Inflationary Models with Exact Mode Functions

Authors: Xin-zhou Li, Xiang-hua Zhai

Abstract: We show two analytical solutions of the tachyon inflation for which the spectrum of curvature (density) perturbations can be calculated exactly to linear order, ignoring both gravity and the self-interactions of the tachyon field . The main feature of these solutions is that the spectral indices are independent with scale. We show two analytical solutions of the tachyon inflation for which the spectrum of curvature (density) perturbations can be calculated exactly to linear order, ignoring both gravity and the self-interactions of the tachyon field . The main feature of these solutions is that the spectral indices are independent with scale. △ Less

Submitted 9 January, 2003; originally announced January 2003.

Comments: 5 pages, no figure, to appear in Phys. Rev. D

Journal ref: Phys.Rev. D67 (2003) 067501