Application of Regularization Methods in the Sky Map Reconstruction of the Tianlai Cylinder Pathfinder Array

The Tianlai cylinder pathfinder is a radio interferometer array to test 21 cm intensity mapping techniques in the post-reionization era.It works in passive drift scan mode to survey the sky visible in the northern hemisphere.To deal with the large instantaneous field of view and the spherical sky,we decompose the drift scan data into m-modes,which are linearly related to the sky intensity.The sky map is reconstructed by solving the linear interferometer equations.Due to incomplete uv coverage of the interferometer baselines,this inverse problem is usually ill-posed,and regularization method is needed for its solution.In this paper,we use simulation to investigate two frequently used regularization methods,the Truncated Singular Value Decomposition(TSVD),and the Tikhonov regularization techniques.Choosing the regularization parameter is very important for its application.We employ the generalized cross validation method and the L-curve method to determine the optimal value.We compare the resulting maps obtained with the different regularization methods,and for the different parameters derived using the different criteria.While both methods can yield good maps for a range of regularization parameters,in the Tikhonov method the suppression of noisy modes are more gradually applied,produce more smooth maps which avoids some visual artefacts in the maps generated with the TSVD method.

The FRB-searching Pipeline of the Tianlai Cylinder Pathfinder Array

This paper presents the design,calibration,and survey strategy of the Fast Radio Burst(FRB)digital backend and its real-time data processing pipeline employed in the Tianlai Cylinder Pathfinder Array.The array,consisting of three parallel cylindrical reflectors and equipped with 96 dual-polarization feeds,is a radio interferometer array designed for conducting drift scans of the northern celestial semi-sphere.The FRB digital backend enables the formation of 96 digital beams,effectively covering an area of approximately 40 square degrees with the 3 dB beam.Our pipeline demonstrates the capability to conduct an automatic search of FRBs,detecting at quasi-realtime and classifying FRB candidates automatically.The current FRB searching pipeline has an overall recall rate of88%.During the commissioning phase,we successfully detected signals emitted by four well-known pulsars:PSR B0329+54,B2021+51,B0823+26,and B2020+28.We report the first discovery of an FRB by our array,designated as FRB 20220414A.We also investigate the optimal arrangement for the digitally formed beams to achieve maximum detection rate by numerical simulation.

The electrical design of a membrane antenna for a lunar-based low-frequency radio telescope

Detecting primordial fluctuations from the cosmic dark ages requires extremely large low-frequency radio telescope arrays deployed on the far side of the Moon.The antenna of such an array must be lightweight,easily storable and transportable,deployable on a large scale,durable,and capable of good electrical performance.A membrane antenna is an excellent candidate to meet these criteria.We study the design of a low-frequency membrane antenna for a lunar-based low-frequency(<30 MHz)radio telescope constructed from polyimide film widely used in aerospace applications,owing to its excellent dielectric properties and high stability as a substrate material.We first design and optimize an antenna in free space through dipole deformation and coupling principles,then simulate an antenna on the lunar surface with a simple lunar soil model,yielding an efficiency greater than 90%in the range of 12-19 MHz and greater than 10%in the range of 5-35 MHz.The antenna inherits the omni-directional radiation pattern of a simple dipole antenna in the 5-30 MHz frequency band,giving a large field of view and allowing detection of the 21 cm global signal when used alone.A demonstration prototype is constructed,and its measured electrical property is found to be consistent with simulated results using|S11|measurements.This membrane antenna can potentially fulfill the requirements of a lunar low-frequency array,establishing a solid technical foundation for future large-scale arrays for exploring the cosmic dark ages.

Ultrasonic positioning system for the calibration of central detector

A thorough detector response calibration using radioactive sources is necessary for the Jiangmen Underground Neutrino Observatory. Herein, we discuss the design of a source positioning system based on ultrasonic technology, aiming for a 3-cm precision over the entire35-m diameter detector sphere. A prototype system is constructed and demonstrated for the experiment.

QCD Axial Anomaly Enhances theηη'Decay of the Hybrid Candidateη1(1855)

We study the hybrid mesons with the exotic quantum number I^(G)J^(PG)=0+1-+and investigate their decays into theηη’,a1(1260)π,f1(1285)η,f1(1420)η,K*(892)■,K1(1270)■,and K1(1400)■channels.It is found that the QCD axial anomaly enhances the decay width of theηη’channel although this mode is strongly suppressed by the small p-wave phase space.Our results support the interpretation of theη1(1855)recently observed by BESⅢas the ssg hybrid meson of I^(G)J^(PC)=0^(+)1^(-+).The QCD axial anomaly ensures theηη’decay mode to be a characteristic signal of the hybrid nature of theη1(1855).

A Simulation of Calibration and Map-making Errors of the Tianlai Cylinder Pathfinder Array

The Tianlai cylinder array is a pathfinder for developing and testing 21 cm intensity mapping techniques.In this paper,we use numerical simulation to assess how its measurement is affected by thermal noise and the errors in calibration and map-making process,and the error in the sky map reconstructed from a drift scan survey.Here we consider only the single frequency,unpolarized case.The beam is modeled by fitting to the electromagnetic simulation of the antenna,and the variations of the complex gains of the array elements are modeled by Gaussian processes.Mock visibility data are generated and run through our data processing pipeline.We find that the accuracy of the current calibration is limited primarily by the absolute calibration,where the error comes mainly from the approximation of a single dominating point source.We then studied the m-mode map-making with the help of Moore-Penrose inverse.We find that discarding modes with singular values smaller than a threshold could generate visible artifacts in the map.The impacts of the residue variation of the complex gain and thermal noise are also investigated.The thermal noise in the map varies with latitude,being minimum at the latitude passing through the zenith of the telescope.The angular power spectrum of the reconstructed map show that the current Tianlai cylinder pathfinder,which has a shorter maximum baseline length in the North-South direction,can measure modes up to l■2πb_(NS)/λ~200 very well,but would lose a significant fraction of higher angular modes when noise is present.These results help us to identify the main limiting factors in our current array configuration and data analysis procedure,and suggest that the performance can be improved by reconfiguration of the array feed positions.

Numerical simulation of radio signal from extended air showers

The burst of radio emission by an extensive air shower provides a promising alternative for detecting ultra-high energy cosmic rays. We have developed an inde- pendent numerical program to simulate these radio signals. Our code is based on a microscopic treatment, with both the geosynchrotron radiation and charge included. Here we give the first presentation of our basic program and its results. When the time-domain signals for different polarizations are computed, we find that the pulses take on a bipolar pattern and the spectrum is suppressed towards the lower frequencies. We investigate how showers at different heights in the atmosphere contribute to the total signal, and examine the signal strength and distribution at sites with different elevations. We also study the signal from showers with different inclination angles and azimuth directions. In all these cases we find the charge excess effect is important.

Classifying Globular Clusters and Applying them to Estimate the mass of the Milky Way

We combine the kinematics of 159 globular clusters(GCs) provided by the Gaia Early Data Release 3(EDR3) with other observational data to classify the GCs,and to estimate the mass of the Milky Way(MW).We use the agemetallicity relation,integrals of motion,action space and the GC orbits to identify the GCs as either formed in situ(Bulge and Disk) or ex situ(via accretion).We find that 45.3% have formed in situ,while 38.4% may be related to known merger events:Gaia-Sausage-Enceladus,the Sagittarius dwarf galaxy,the Helmi streams,the Sequoia galaxy and the Kraken galaxy.We also further identify three new sub-structures associated with the Gaia-Sausage-Enceladus.The remaining 16.3% of GCs are unrelated to the known mergers and thought to be from small accretion events.We select 46 GCs which have radii 8.0 <r<37.3 kpc and obtain the anisotropy parameter β=0.315_(-0.049)^(+0.055),which is lower than the recent result using the sample of GCs in Gaia Data Release 2,but still in agreement with it by considering the error bar.By using the same sample,we obtain the MW mass inside the outermost GC as M(<37.3 kpc)=0.423_(-0.02)^(+0.02)×10^(12)M_(⊙),and the corresponding M_(200)=1.11_(-0.18)^(+0.25)×10^(12)M_(⊙).The estimated mass is consistent with the results in many recent studies.We also find that the estimated β and mass depend on the selected sample of GCs.However,it is difficult to determine whether a GC fully traces the potential of the MW.

The Electromagnetic Characteristics of the Tianlai Cylindrical Pathfinder Array

A great challenge for 21 cm intensity mapping experiments is the strong foreground radiation which is orders of magnitude brighter than the 21 cm signal.Removal of the foreground takes advantage of the fact that its frequency spectrum is smooth while the redshifted 21 cm signal spectrum is stochastic.However,a complication is the nonsmoothness of the instrument response.This paper describes the electromagnetic simulation of the Tianlai cylinder array,a pathfinder for 21 cm intensity mapping experiments.Due to the vast scales involved,a direct simulation requires a large amount of computing resources.We have made the simulation practical by using a combination of methods:first simulate a single feed,then an array of feed units,finally with the feed array and a cylindrical reflector together,obtain the response for a single cylinder.We studied its radiation pattern,bandpass response and the effects of mutual coupling between feed units,and compared the results with observation.Many features seen in the measurement result are reproduced well in the simulation,especially the oscillatory features which are associated with the standing waves on the reflector.The mutual coupling between feed units is quantified with Sparameters,which decrease as the distance between the two feeds increases.Based on the simulated S-parameters,we estimate the correlated noise which has been seen in the visibility data,and the results show very good agreement with the data in both magnitude and frequency structures.These results provide useful insights on the problem of 21 cm signal extraction for real instruments.

Semi-numerical simulation of reionization with semi-analytical modeling of galaxy formation

In a semi-numerical model of reionization, the evolution of ionization fraction is approximately simulated by the criterion of ionizing photon to baryon ratio. We incorporate a semi-analytical model of galaxy formation based on the Millennium II N-body simulation into the semi-numerical modeling of reionization. The semianalytical model is used to predict the production of ionizing photons, then we use the semi-numerical method to model the reionization process. Such an approach allows more detailed modeling of the reionization, and also connects observations of galaxies at low and high redshifts to the reionization history. The galaxy formation model we use was designed to match the low-z observations, and it also fits the high redshift luminosity function reasonably well, but its prediction about star formation falls below the observed value, and we find that it also underpredicts the stellar ionizing photon production rate, hence the reionization cannot be completed at z ~ 6. We also consider simple modifications of the model with more top heavy initial mass functions, which can allow the reionization to occur at earlier epochs. The incorporation of the semi-analytical model may also affect the topology of the HI regions during the epoch of reionization, and the neutral regions produced by our simulations with the semi-analytical model, which appeared less poriferous than the simple halo-based models.

Evaluation of growth characteristics of Aspergillus parasiticus inoculated in different culture media by shortwave infrared(SWIR) hyperspectral imaging

The growth characteristics of Aspergillus parasitic us incubated on two culture media were ex-amined using shortwave infrared(SWIR,1000-2500 nm)hyperspectral imaging(HSI)in this work.HSI images of the A.parasiticus colonies growing on rose bengal medium(RBM)and maize agar medium(MAM)were recorded daily for 6 days.The growth phases of A.parasiticus were indicated through the pixel number and average spectra of colonies.On score plot of the first principal component(PC1)and PC2,four growth zones with varying mycelium densities were identified.Eight characteristic wavelengths(1095,1145,1195,1279,1442,1655,1834 and 1929 nm)were selected from PC1 loading,average spectra of each colony as well as each growth zone.F urthermore,support vector machine(S VM)classifier based on the eight wavelengths was built,and the classification accuracies for the four zones(from outer to inner zones)on the colonies on RBM were 99.77%,9935%,99.75%and 99.60%and 99.77%,9939%,99.31%and 98.22%for colonies on MAM.In addition,a new score plot of PC2 and PC3 was used to differ-entiate the colonies incubated on RBM and MAM for 6 days.Then characteristic wavelengths of 1067,1195,1279,1369,1459,1694,1834 and 1929 nm were selected from the loading of PC2 and PCg.Based on them,a new SVM model was developed to diferentiate colonies on RBM and MAM with accuracy of 100.00%and 9999%,respectively.In conclusion,SWIR hyperspectral image is a powerful tool for evaluation of growth characteristics of A.parasiticus incubated in diferent culture media.

Perfect DD^(*) Molecular Prediction Matching the T_(cc) Observation at LHCb

In 2012,we investigated the possible molecular states composed of two charmed mesons[Phys.Rev.D 88(2013)114008;2012 arXiv:1211.5007[hep-ph]].The D^(*)D𝐷system with the quantum numbers of I(J^(P))=0(1^(+))was found to be a good candidate of the loosely bound molecular state.This state is very close to theD*D threshold with a binding energy around 0.47 MeV.This prediction was confirmed by the new LHCb observation of𝑇T^(+)_(cc)[see Franz Muheim’s talk at the European Physical Society conference on high energy physics 2021].

Reflections and standing waves on the Tianlai cylinder array

In 21 cm intensity mapping,the spectral smoothness of the foreground is exploited to separate it from the much weaker 21 cm signal.However,the non-smooth frequency response of the instrument complicates this process.Reflections and standing waves generate modulations on the frequency response.Here we report the analysis of the standing waves in the bandpass of the signal channels of the Tianlai Cylinder Array.By Fourier transforming the bandpass into the delay time domain,we find various standing waves generated on the telescope.A standing wave with time delay at~142 ns is most clearly identified which is produced in the 15-meter feed cable.We also identify a strong peak at a shorter delay ofτ<50 ns,which may be a mix of the standing wave between the reflector and feed,and the standing wave on the 4 m intermediate frequency(IF)cable.We also show that a smoother frequency response could be partially recovered by removing the reflection-inducted modulations.However,the standing wave on the antenna is direction-dependent,which poses a more difficult challenge for high precision calibration.

Reionizing Islands with Inhomogeneous Recombinations

Observations are beginning to constrain the history of the epoch of reionization(EoR).Modeling the reionization process is indispensable to interpret the observations,to infer the properties of ionizing sources,and to probe the various astrophysical processes from the observational data.Here we present an improved version of the seminumerical simulation islandFAST,by incorporating inhomogeneous recombinations and a corresponding inhomogeneous ionizing background,and simulate the reionization process of neutral islands during the late EoR.We find that the islands are more fragmented in models with inhomogeneous recombinations than the case with a homogeneous recombination number.In order to investigate the effects of basic assumptions in the reionization modeling,we compare the results from islandFAST with those from 21cmFAST for the same assumptions on the ionizing photon sources and sinks,to find how the morphology of the ionization field and the reionization history depend on the different treatments of these two models.Such systematic bias should be noted when interpreting the upcoming observations.

TeV cosmic-ray proton and helium spectra in the myriad model Ⅱ

Recent observations show that the spectra of cosmic ray nuclei start to harden above ～ 102 GeV, which contradicts the conventional steady-state cosmic ray model. We had suggested that this anomaly is due to the propagation effect of cosmic rays released from local young cosmic ray sources; the total flux of cosmic rays should be computed with the Myriad Model, where a contribution from sources in the local catalog is added to the background. However, although the hardening could be elegantly explained in this model, the model parameters obtained from the fit are skewed toward a region with fast diffusion and a low supernova rate in the Galaxy, in disagreement with other observations. We further explore this model in order to set up a concordant picture. Two possible improvements related to cosmic ray sources are considered. First, instead of the usual axisymmetric disk model, we examine a spiral model for the source distribution. Second, for nearby and young sources which are necessary to explain the hardening, we allow for an energy-dependent escape. We find that a major improvement comes from incorporating an energy-dependent escape time for local sources, and with both modifications not only are the cosmic ray proton and helium anomalies resolved, but also the parameters attain values in a reasonable range compatible with other analyses.

A Fast Radio Burst Backend for the Tianlai Dish Pathfinder Array

The Tianlai Dish Pathfinder Array is a radio interferometer array consisting of 16 six-meter dish antennas.The original digital backend integration time is at the seconds level,designed for H I intensity mapping experiment.A new digital backend with millisecond response is added to enable it to search for fast radio burst during its observations.The design and calibration of this backend,and the real time search pipeline for it are described in this paper.It is capable of forming 16 digital beams for each linear polarization,covering an area of 19.6 square degrees.The search pipeline is capable of searching for,recording and classifying FRBs automatically in real time.In commissioning,we succeeded in capturing the signal pulses from the pulsars PSR B0329+54 and B2021+51.

Cosmological constraint on Brans-Dicke Model

We combine new Cosmic Microwave Background （CMB） data from Planck with Baryon Acoustic Oscillation （BAO） data to constrain the Brans-Dicke （BD） the- ory, in which the gravitational constant G evolves with time. Observations of type Ia supernovae （SNela） provide another important set of cosmological data, as they may be regarded as standard candles after some empirical corrections. However, in theories that include modified gravity like the BD theory, there is some risk and complication when using the SNIa data because their luminosity may depend on G. In this paper, we assume a power law relation between the SNIa luminosity and G, but treat the power index as a free parameter. We then test whether the differ- ence in distances measured with SNIa data and BAO data can be reduced in such a model. We also constrain the BD theory and cosmological parameters by making a global fit with the CMB, BAO and SNIa data set. For the CMB＋BAO＋SNIa data set, we find 0.08 × 10-2 〈 ζ 〈 0.33 ×10-2 at the 68% confidence level （CL） and -0.01 × 10-2 〈 ζ〈 0.43 × 10-2 at the 95% CL, where ζ is related to the BD parameter w by ζ = In（1 ＋ 1/w）.

Primordial non-Gaussianity from LAMOST surveys

The primordial non-Gaussianity （PNG） in the matter density perturbation is a very powerful probe of the physics of the very early Universe. The local PNG can induce a distinct scale-dependent bias on the large scale structure distribution of galaxies and quasars, which could be used for constraining it. We study the detection limits of PNG from the surveys of the LAMOST telescope. The cases of the main galaxy survey, the luminous red galaxy （LRG） survey, and the quasar survey of dif- ferent magnitude limits are considered. We find that the Mainl sample （i.e. the main galaxy survey which is one magnitude deeper than the SDSS main galaxy survey, or r 〈 18.8） could only provide a very weak constraint on PNG. For the Main2 sample （r 〈 19.8） and the LRG survey, the 2or （95.5%） limits on the PNG parameter fNL are ｜ fNL｜〈 145 and ｜ fNL｜ 〈 114 respectively, which are comparable to the current limit from cosmic microwave background （CMB） data. The quasar survey could provide a much more stringent constraint, and we find that the 2a limit for ｜fNL｜ is between 50 and 103, depending on the magnitude limit of the survey. With Planck-like priors on cosmological parameters, the quasar survey with g 〈 21.65 would improve the constraint to I fNLI 〈 43 （2or）. We also discuss the possibility of further tightening the constraint by using the relative bias method proposed by Seljak.

Dark matter direct search sensitivity of the PandaX-4T experiment

The Panda X-4T experiment, a 4-ton scale dark matter direct detection experiment, is being planned at the China Jinping Underground Laboratory. In this paper we present a simulation study of the expected background in this experiment. In a 2.8-ton fiducial mass and the signal region between 1-10 keV electron equivalent energy, the total electron recoil background is found to be 4.9 × 10^(-5) kg^(-1) d^(-1) keV^(-1). The nuclear recoil background in the same region is 2.8 × 10^(-7) kg^(-1) d^(-1) keV^(-1). With an exposure of 5.6 ton-years, the sensitivity of Panda X-4 T could reach a minimum spin-independent dark matter-nucleon cross section of 6 × 10^(-48) cm^2 at a dark matter mass of 40 Ge V/c^2.

Implications of the Z_(cs)(3985) and Z_(cs)(4000) as two different states

Recently,the hidden charm tetraquark states Z_(cs)(3985)and Z_(cs)(4000)with strangeness were observed by the BESⅢ and LHCb collaborations,respectively,which are great breakthroughs for exploring exotic quantum chromodynamics(QCD)structures.The first and foremost question is whether they are the same state.In this work,we explore the implications of the narrower state Z_(cs)(3985)in BESⅢ and the wider one Z_(cs)(4000)in LHCb as two different states.Within a solvable nonrelativistic effective field theory,we include the possible violations of heavy quark spin symmetry and SU(3)flavor symmetry in a comprehensive approach.If Z_(cs)(3985)and Z_(cs)(4000)are two different states,our results show that Z_(cs)(4000)/Z_(cs)(3985)is the pure(|D_(s)^(*)D>+/-|D_(s)D^(*)>)/√2 state,and the SU(3)flavor partner of Zc(3900)is Z_(cs)(4000)rather than the Z_(cs)(3985).Another two important consequences are the existence of a tensor D_(s)^(*)D^(*) resonance with mass about 4126 MeV and width 13 MeV,and the suppression of the decay mode Z_(cs)(3985)J/ψK.The two consequences can be tested in experiments and distinguish the two-state interpretation from the one-state scheme.