Research in Astronomy and Astrophysics

In order to evaluate how much Total Solar Irradiance (TSI) has influenced Northern Hemisphere surface air temperature trends, it is important to have reliable estimates of both quantities. Sixteen different estimates of the changes in TSI since at least the 19^th century were compiled from the literature. Half of these estimates are "low variability" and half are "high variability". Meanwhile, five largely-independent methods for estimating Northern Hemisphere temperature trends were evaluated using: 1) only rural weather stations; 2) all available stations whether urban or rural (the standard approach); 3) only sea surface temperatures; 4) tree-ring widths as temperature proxies; 5) glacier length records as temperature proxies. The standard estimates which use urban as well as rural stations were somewhat anomalous as they implied a much greater warming in recent decades than the other estimates, suggesting that urbanization bias might still be a problem in current global temperature datasets – despite the conclusions of some earlier studies. Nonetheless, all five estimates confirm that it is currently warmer than the late 19^th century, i.e., there has been some "global warming" since the 19^th century. For each of the five estimates of Northern Hemisphere temperatures, the contribution from direct solar forcing for all sixteen estimates of TSI was evaluated using simple linear least-squares fitting. The role of human activity on recent warming was then calculated by fitting the residuals to the UN IPCC's recommended "anthropogenic forcings" time series. For all five Northern Hemisphere temperature series, different TSI estimates suggest everything from no role for the Sun in recent decades (implying that recent global warming is mostly human-caused) to most of the recent global warming being due to changes in solar activity (that is, that recent global warming is mostly natural). It appears that previous studies (including the most recent IPCC reports) which had prematurely concluded the former, had done so because they failed to adequately consider all the relevant estimates of TSI and/or to satisfactorily address the uncertainties still associated with Northern Hemisphere temperature trend estimates. Therefore, several recommendations on how the scientific community can more satisfactorily resolve these issues are provided.

Since 2007, the Intergovernmental Panel on Climate Change (IPCC) has heavily relied on the comparison between global climate model hindcasts and global surface temperature (ST) estimates for concluding that post-1950s global warming is mostly human-caused. In Connolly et al., we cautioned that this approach to the detection and attribution of climate change was highly dependent on the choice of Total Solar Irradiance (TSI) and ST data sets. We compiled 16 TSI and five ST data sets and found by altering the choice of TSI or ST, one could (prematurely) conclude anything from the warming being "mostly human-caused" to "mostly natural." Richardson and Benestad suggested our analysis was "erroneous" and "flawed" because we did not use a multilinear regression. They argued that applying a multilinear regression to one of the five ST series re-affirmed the IPCC's attribution statement. They also objected that many of the published TSI data sets were out-of-date. However, here we show that when applying multilinear regression analysis to an expanded and updated data set of 27 TSI series, the original conclusions of Connolly et al. are confirmed for all five ST data sets. Therefore, it is still unclear whether the observed warming is mostly human-caused, mostly natural or some combination of both.

Observing and timing a group of millisecond pulsars with high rotational stability enables the direct detection of gravitational waves (GWs). The GW signals can be identified from the spatial correlations encoded in the times-of-arrival of widely spaced pulsar-pairs. The Chinese Pulsar Timing Array (CPTA) is a collaboration aiming at the direct GW detection with observations carried out using Chinese radio telescopes. This short article serves as a "table of contents" for a forthcoming series of papers related to the CPTA Data Release 1 (CPTA DR1) which uses observations from the Five-hundred-meter Aperture Spherical radio Telescope. Here, after summarizing the time span and accuracy of CPTA DR1, we report the key results of our statistical inference finding a correlated signal with amplitude $\mathrm{log}{A}_{{\rm{c}}}=-14.4{\,}_{-2.8}^{+1.0}$ for spectral index in the range of α ∈ [ − 1.8, 1.5] assuming a GW background (GWB) induced quadrupolar correlation. The search for the Hellings–Downs (HD) correlation curve is also presented, where some evidence for the HD correlation has been found that a 4.6σ statistical significance is achieved using the discrete frequency method around the frequency of 14 nHz. We expect that the future International Pulsar Timing Array data analysis and the next CPTA data release will be more sensitive to the nHz GWB, which could verify the current results.

The Large sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) general survey is a spectroscopic survey that will eventually cover approximately half of the celestial sphere and collect 10 million spectra of stars, galaxies and QSOs. Objects in both the pilot survey and the first year regular survey are included in the LAMOST DR1. The pilot survey started in October 2011 and ended in June 2012, and the data have been released to the public as the LAMOST Pilot Data Release in August 2012. The regular survey started in September 2012, and completed its first year of operation in June 2013. The LAMOST DR1 includes a total of 1202 plates containing 2 955 336 spectra, of which 1 790 879 spectra have observed signal-to-noise ratio (SNR) ≥ 10. All data with SNR ≥ 2 are formally released as LAMOST DR1 under the LAMOST data policy. This data release contains a total of 2 204 696 spectra, of which 1 944 329 are stellar spectra, 12 082 are galaxy spectra and 5017 are quasars. The DR1 not only includes spectra, but also three stellar catalogs with measured parameters: late A,FGK-type stars with high quality spectra (1061 918 entries), A-type stars (100 073 entries), and M-type stars (121 522 entries). This paper introduces the survey design, the observational and instrumental limitations, data reduction and analysis, and some caveats. A description of the FITS structure of spectral files and parameter catalogs is also provided.

In this work, the γ-ray photon flux, photon spectral index (α_ph), variability index ($\mathrm{log}\,V.I.$), and the synchrotron peak frequency ($\mathrm{log}\,{\nu }_{{\rm{p}}}$) are compiled for 851 common blazars from the 3FGL and 4FGL catalogs and Fan et al. to investigate variability properties for Fermi blazars. Our calculations and analyses reach following results: (1) the averaged luminosity, spectral index, and variability index of FSRQs are higher than those of BL Lacs for the whole sample. (2) It is found that the spectral index variation is closely anti-correlated with the luminosity variation implying that the spectrum becomes harder when the source becomes brighter in the γ-ray band. (3) Positive correlations are found between the photon spectral index and both γ-ray luminosity and variability index ($\mathrm{log}\,V.I.$) for the whole sample, but anti-correlations are found in the two correlations for FSRQs. For BL Lac subclass, there is a marginal anti-correlation between the photon spectral index and both γ-ray luminosity, and a positive correlation between the photon spectral index and the variability index ($\mathrm{log}\,V.I.$). We think those two positive correlations found for the whole sample are apparent. (4) We adopted the SVM machine learning method to separate BL Lacs and FSRQs in the ${\alpha }_{\mathrm{ph}}\,\,\mathrm{versus}\,\,\mathrm{log}\,{L}_{\gamma }$ and ${\alpha }_{\mathrm{ph}}\,\,\mathrm{versus}\,\,\mathrm{log}\,V.I.$ plots and proposed that a BCU is an FSRQ candidate if it satisfies $\mathrm{log}\,V.I.\,\gt \,-10.119\alpha +24.855$, or ${\alpha }_{\mathrm{ph}}\gt -0.048\,\mathrm{log}\,{L}_{\gamma }+4.498$, otherwise, it is a BL Lac candidate. Our classification results are quite consistent with those by Kang et al (2019).

The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST, also called the Guo Shou Jing Telescope) is a special reflecting Schmidt telescope. LAMOST's special design allows both a large aperture (effective aperture of 3.6 m–4.9 m) and a wide field of view (FOV) (5°). It has an innovative active reflecting Schmidt configuration which continuously changes the mirror's surface that adjusts during the observation process and combines thin deformable mirror active optics with segmented active optics. Its primary mirror (6.67 m × 6.05 m) and active Schmidt mirror (5.74m × 4.40m) are both segmented, and composed of 37 and 24 hexagonal sub-mirrors respectively. By using a parallel controllable fiber positioning technique, the focal surface of 1.75 m in diameter can accommodate 4000 optical fibers. Also, LAMOST has 16 spectrographs with 32 CCD cameras. LAMOST will be the telescope with the highest rate of spectral acquisition. As a national large scientific project, the LAMOST project was formally proposed in 1996, and approved by the Chinese government in 1997. The construction started in 2001, was completed in 2008 and passed the official acceptance in June 2009. The LAMOST pilot survey was started in October 2011 and the spectroscopic survey will launch in September 2012. Up to now, LAMOST has released more than 480000 spectra of objects. LAMOST will make an important contribution to the study of the large-scale structure of the Universe, structure and evolution of the Galaxy, and cross-identification of multi-waveband properties in celestial objects.

The East Asian Very Long Baseline Interferometry (VLBI) Network (EAVN) is a rapidly evolving international VLBI array that is currently promoted under joint efforts among China, Japan and Korea. EAVN aims at forming a joint VLBI Network by combining a large number of radio telescopes distributed over East Asian regions. After the combination of the Korean VLBI Network (KVN) and the VLBI Exploration of Radio Astrometry (VERA) into KaVA, further expansion with the joint array in East Asia is actively promoted. Here we report the first imaging results (at 22 and 43 GHz) of bright radio sources obtained with KaVA connected to Tianma 65-m and Nanshan 26-m Radio Telescopes in China. To test the EAVN imaging performance for different sources, we observed four active galactic nuclei (AGN) having different brightness and morphology. As a result, we confirmed that the Tianma 65-m Radio Telescope (TMRT) significantly enhances the overall array sensitivity, a factor of 4 improvement in baseline sensitivity and 2 in image dynamic range compared to the case of KaVA only. The addition of the Nanshan 26-m Radio Telescope (NSRT) further doubled the east-west angular resolution. With the resulting high-dynamic-range, high-resolution images with EAVN (KaVA+TMRT+NSRT), various fine-scale structures in our targets, such as the counter-jet in M87, a kink-like morphology of the 3C 273 jet and the weak emission in other sources are successfully detected. This demonstrates the powerful capability of EAVN to study AGN jets and to achieve other science goals in general. Ongoing expansion of EAVN will further enhance the angular resolution, detection sensitivity and frequency coverage of the network.

Discovery of pulsars is one of the main goals for large radio telescopes. The Five-hundred-meter Aperture Spherical radio Telescope (FAST), that incorporates an L-band 19-beam receiver with a system temperature of about 20 K, is the most sensitive radio telescope utilized for discovering pulsars. We designed the snapshot observation mode for a FAST key science project, the Galactic Plane Pulsar Snapshot (GPPS) survey, in which every four nearby pointings can observe a cover of a sky patch of 0.1575 square degrees through beam-switching of the L-band 19-beam receiver. The integration time for each pointing is 300 seconds so that the GPPS observations for a cover can be made in 21 minutes. The goal of the GPPS survey is to discover pulsars within the Galactic latitude of ± 10° from the Galactic plane, and the highest priority is given to the inner Galaxy within ± 5°. Up to now, the GPPS survey has discovered 201 pulsars, including currently the faintest pulsars which cannot be detected by other telescopes, pulsars with extremely high dispersion measures (DMs) which challenge the currently widely used models for the Galactic electron density distribution, pulsars coincident with supernova remnants, 40 millisecond pulsars, 16 binary pulsars, some nulling and mode-changing pulsars and rotating radio transients (RRATs). The follow-up observations for confirmation of new pulsars have polarization-signals recorded for polarization profiles of the pulsars. Re-detection of previously known pulsars in the survey data also leads to significant improvements in parameters for 64 pulsars. The GPPS survey discoveries are published and will be updated at http://zmtt.bao.ac.cn/GPPS/.

The Milky Way is a spiral galaxy with the Schechter characteristic luminosity L_*, thus an important anchor point of the Hubble sequence of all spiral galaxies. Yet the true appearance of the Milky Way has remained elusive for centuries. We review the current best understanding of the structure and kinematics of our home galaxy, and present an updated scientifically accurate visualization of the Milky Way structure with almost all components of the spiral arms, along with the COBE image in the solar perspective. The Milky Way contains a strong bar, four major spiral arms, and an additional arm segment (the Local arm) that may be longer than previously thought. The Galactic boxy bulge that we observe is mostly the peanut-shaped central bar viewed nearly end-on with a bar angle of ∼ 25° – 30° from the Sun-Galactic center line. The bar transitions smoothly from a central peanut-shaped structure to an extended thin part that ends around R ∼ 5 kpc. The Galactic bulge/bar contains ∼ 30% – 40% of the total stellar mass in the Galaxy. Dynamical modelling of both the stellar and gas kinematics yields a bar pattern rotation speed of ∼ 35 – 40 km s⁻¹kpc ⁻¹, corresponding to a bar rotation period of ∼ 160 – 180 Myr. From a galaxy formation point of view, our Milky Way is probably a pure-disk galaxy with little room for a significant merger-made, "classical" spheroidal bulge, and we give a number of reasons why this is the case.

The Five-hundred-meter Aperture Spherical radio Telescope (FAST) has passed national acceptance and finished one pilot cycle of 'Shared-Risk' observations. It will start formal operation soon. In this context, this paper describes testing results of key fundamental parameters for FAST, aiming to provide basic support for observation and data reduction of FAST for scientific researchers. The 19-beam receiver covering 1.05–1.45 GHz was utilized for most of these observations. The fluctuation in electronic gain of the system is better than 1% over 3.5 hours, enabling enough stability for observations. Pointing accuracy, aperture efficiency and system temperature are three key parameters for FAST. The measured standard deviation of pointing accuracy is 7.9'', which satisfies the initial design of FAST. When zenith angle is less than 26.4°, the aperture efficiency and system temperature around 1.4 GHz are ∼0.63 and less than 24 K for central beam, respectively. The sensitivity and stability of the 19-beam backend are confirmed to satisfy expectation by spectral Hi observations toward NGC 672 and polarization observations toward 3C 286. The performance allows FAST to take sensitive observations for various scientific goals, from studies of pulsars to galaxy evolution.

Attitude is one of the crucial parameters for space objects and plays a vital role in collision prediction and debris removal. Analyzing light curves to determine attitude is the most commonly used method. In photometric observations, outliers may exist in the obtained light curves due to various reasons. Therefore, preprocessing is required to remove these outliers to obtain high quality light curves. Through statistical analysis, the reasons leading to outliers can be categorized into two main types: first, the brightness of the object significantly increases due to the passage of a star nearby, referred to as "stellar contamination," and second, the brightness markedly decreases due to cloudy cover, referred to as "cloudy contamination." The traditional approach of manually inspecting images for contamination is time-consuming and labor-intensive. However, we propose the utilization of machine learning methods as a substitute. Convolutional Neural Networks and SVMs are employed to identify cases of stellar contamination and cloudy contamination, achieving F1 scores of 1.00 and 0.98 on a test set, respectively. We also explore other machine learning methods such as ResNet-18 and Light Gradient Boosting Machine, then conduct comparative analyses of the results.

Determining asteroid properties provides valuable physical insights but inverting them from photometric lightcurves remains computationally intensive. This paper presents a new approach that combines a simplified Cellinoid shape model with the Parallel Differential Evolution (PDE) algorithm to accelerate inversion. The PDE algorithm is more efficient than the Differential Evolution algorithm, achieving an extraordinary speedup of 37.983 with 64 workers on multicore CPUs. The PDE algorithm accurately derives period and pole values from simulated data. The analysis of real asteroid lightcurves validates the method's reliability: in comparison with results published elsewhere, the PDE algorithm accurately recovers the rotational periods and, given adequate viewing geometries, closely matches the pole orientations. The PDE approach converges to solutions within 20,000 iterations and under one hour, demonstrating its potential for large-scale data analysis. This work provides a promising new tool for unveiling asteroid physical properties by overcoming key computational bottlenecks.

Infrared dust bubbles play an important role in the study of star formation and the evolution of the interstellar medium. In this work, we study the infrared dust bubble N75 and the infrared dark cloud G38.93 mainly using the tracers C¹⁸O, HCO⁺, HNC and N₂H⁺ observed by the 30 m IRAM telescope. We also study the targets using data from large-scale surveys: GLIMPSE, MIPSGAL, GRS, NRAO VLA Sky Survey and Bolocam Galactic Plane Survey. We found that the C¹⁸O emission is morphologically similar to the Spitzer IRAC 8.0 μm emission. The 1.1 mm cold dust emission of G38.93 shows an elongated structure from southwest to northeast. The ionized gas from G38.93 is surrounded by polycyclic aromatic hydrocarbon emission, which may be excited by radiation from G38.93. We found that the identified young stellar objects tend to cluster around G38.93 and are mostly in class II, with several class I cases distributed around N75, but no class II examples. We also found evidence of expanding feedback, which could have triggered star formation.

We have conducted a comprehensive investigation into the bright single pulse emission from PSR B1133+16 using the Giant Metrewave Radio Telescope. High time resolution data (61 μs) were obtained at a center frequency of 322 MHz with a bandwidth of 32 MHz over a continuous observation period of 7.45 hr. A total of 1082 bright pulses were sporadically detected with peak flux densities ranging from 10 to 23 times stronger than the average pulse profile. However, no giant pulse-like emission with a relative pulse energy larger than 10 and extremely short duration was detected, indicating that these bright pulses cannot be categorized as giant pulse emission. The majority of these bright pulses are concentrated in pulse phases at both the leading and trailing windows of the average pulse profile, with an occurrence ratio of approximately 2.74. The pulse energy distribution for all individual pulses can be described by a combination of two Gaussian components and a cutoff power-law with an index of α = − 3.2. An updated nulling fraction of 15.35% ± 0.45% was determined from the energy distribution. The emission of individual pulses follows a log-normal distribution in peak flux density ratio. It is imperative that regular phase drifting in bright pulse sequence is identified in both the leading and trailing components for the first time. Possible physical mechanisms are discussed in detail to provide insights into these observations.

The initial condition of high-mass star formation is a complex area of study because of the high densities (${n}_{{{\rm{H}}}_{2}}\gt {10}^{6}$ cm⁻³) and low temperatures (T_dust < 18 K) involved. Under such conditions, many molecules become depleted from the gas phase by freezing out onto dust grains. However, the N-bearing and deuterated species could remain gaseous under these extreme conditions, suggesting that they may serve as ideal tracers. In this paper, using the Plateau de Bure Interferometer and Very Large Array observations at 1.3 mm, 3.5 mm, and 1.3 cm, we investigate the possible habitats for NH₃, NH₂D, H¹³CN, HC¹⁵N, SO, and C¹⁸O in eight massive precluster and protocluster clumps G18.17, G18.21, G23.97N, G23.98, G23.44, G23.97S, G25.38, and G25.71. We found that the NH₃ cores are in good agreement with the 3.5 mm peak emission, but the NH₃ is much more extended than the 3.5 mm emission structure. The SO distributions agree well with the 3.5 mm peaks for the evolved star formation stage, but we did not detect any SO emission in the four earliest star formation sources. C¹⁸O is a poor tracer in conditions of the cold (≲18 K) and dense (≳10⁴ cm⁻³) cores, e.g., the prestellar cores. We also found that the NH₂D cores are mainly located in the temperature range of 13.0–20.0 K, and the NH₂D lines may be strongly depleted above 20 K.

Density functional theory (DFT) is the most versatile electronic structure method used in quantum chemical calculations, and is increasingly applied in astrochemical research. This mini-review provides an overview of the applications of DFT calculations in understanding the chemistry that occurs in star-forming regions. We survey investigations into the formation of biologically relevant compounds such as nucleobases in the interstellar medium, and also cover the formation of both achiral and chiral amino acids, as well as biologically relevant molecules such as sugars, and nitrogen-containing polycyclic aromatic hydrocarbons. Additionally, DFT calculations are used to estimate the potential barriers for chemical reactions in astronomical environments. We conclude by noting several areas that require more research, such as the formation pathways of chiral amino acids, complex sugars, and other biologically important molecules, and the role of environmental factors in the formation of interstellar biomolecules.

Type Ia supernovae (SNe Ia) play a key role in the fields of astrophysics and cosmology. It is widely accepted that SNe Ia arise from thermonuclear explosions of white dwarfs (WDs) in binary systems. However, there is no consensus on the fundamental aspects of the nature of SN Ia progenitors and their actual explosion mechanism. This fundamentally flaws our understanding of these important astrophysical objects. In this review, we outline the diversity of SNe Ia and the proposed progenitor models and explosion mechanisms. We discuss the recent theoretical and observational progress in addressing the SN Ia progenitor and explosion mechanism in terms of the observables at various stages of the explosion, including rates and delay times, pre-explosion companion stars, ejecta–companion interaction, early excess emission, early radio/X-ray emission from circumstellar material interaction, surviving companion stars, late-time spectra and photometry, polarization signals and supernova remnant properties. Despite the efforts from both the theoretical and observational sides, questions of how the WDs reach an explosive state and what progenitor systems are more likely to produce SNe Ia remain open. No single published model is able to consistently explain all observational features and the full diversity of SNe Ia. This may indicate that either a new progenitor paradigm or an improvement in current models is needed if all SNe Ia arise from the same origin. An alternative scenario is that different progenitor channels and explosion mechanisms contribute to SNe Ia. In the next decade, the ongoing campaigns with the James Webb Space Telescope, Gaia and the Zwicky Transient Facility, and upcoming extensive projects with the Vera C. Rubin Observatory's Legacy Survey of Space and Time and the Square Kilometre Array will allow us to conduct not only studies of individual SNe Ia in unprecedented detail but also systematic investigations for different subclasses of SNe Ia. This will advance theory and observations of SNe Ia sufficiently far to gain a deeper understanding of their origin and explosion mechanism.

Magnetic fields play a key role in driving a broad range of dynamic phenomena in the atmospheres of the Sun and other stars. Routine and accurate measurements of the magnetic fields at all the atmospheric layers are of critical importance to understand these magnetic activities, but in the solar and stellar coronae such a measurement is still a challenge due to the weak field strength and the high temperature. Recently, a magnetic-field-induced transition (MIT) of Fe x at 257.26 Å has been proposed for the magnetic field measurements in the solar and stellar coronae. In this review, we present an overview of recent progresses in the application of this method in astrophysics. We start by introducing the theory underlying the MIT method and reviewing the existing atomic data critical for the spectral modeling of Fe x lines. We also discuss the laboratory measurements that verify the potential capability of the MIT technique as a probe for diagnosing the plasma magnetic fields. We then continue by investigating the suitability and accuracy of solar and stellar coronal magnetic field measurements based on the MIT method through forward modeling. Furthermore, we discuss the application of the MIT method to the existing spectroscopic observations obtained by the Extreme-ultraviolet Imaging Spectrometer onboard Hinode. This novel technique provides a possible way for routine measurements of the magnetic fields in the solar and stellar coronae, but still requires further efforts to improve its accuracy. Finally, the challenges and prospects for future research on this topic are discussed.

I review studies of core collapse supernovae (CCSNe) and similar transient events that attribute major roles to jets in powering most CCSNe and in shaping their ejecta. I start with reviewing the jittering jets explosion mechanism that I take to power most CCSN explosions. Neutrino heating does play a role in boosting the jets. I compare the morphologies of some CCSN remnants to planetary nebulae to conclude that jets and instabilities are behind the shaping of their ejecta. I then discuss CCSNe that are descendants of rapidly rotating collapsing cores that result in fixed-axis jets (with small jittering) that shape bipolar ejecta. A large fraction of the bipolar CCSNe are superluminous supernovae (SLSNe). I conclude that modeling of SLSN light curves and bumps in the light curves must include jets, even when considering energetic magnetars and/or ejecta interaction with the circumstellar matter (CSM). I connect the properties of bipolar CCSNe to common envelope jets supernovae (CEJSNe) where an old neutron star or a black hole spirals-in inside the envelope and then inside the core of a red supergiant. I discuss how jets can shape the pre-explosion CSM, as in Supernova 1987A, and can power pre-explosion outbursts (precursors) in binary system progenitors of CCSNe and CEJSNe. Binary interaction also facilitates the launching of post-explosion jets.

The purpose of this paper is to address the question: Using our knowledge of infrared planetary spectroscopy, what can we learn about the atmospheres of exoplanets? In a first part, a simplified classification of exoplanets, assuming thermochemical equilibrium, is presented, based on their masses and their equilibrium temperatures, in order to propose some possible estimations about their atmospheric composition. In the second part, infrared spectra of planets are discussed, in order to see what lessons can be drawn for exoplanetary spectroscopy. In the last part, we consider the solar system as it would appear from a star located in the ecliptic plane. It first appears that the solar system (except in a few specific cases) would not be seen as a multiple system, because, contrary to many exoplanetary systems, the planets are too far from the Sun and the inclinations of their orbits with respect to the ecliptic plane are too high. Primary transit synthetic spectra of solar system planets are used to discuss the relative merits of transmission and direct emission spectroscopy for probing exoplanetary atmospheres.