The interstellar medium (ISM) exhibits complex, multi-scale structures that are challenging to study due to their projection into two-dimensional (2D) column density maps. We present the Volume Density Mapper, a novel algorithm based on constrained diffusion to reconstruct three-dimensional (3D) density distributions of molecular clouds from 2D observations. This method decomposes the column density into multi-scale components, reconstructing a 3D density field that preserves key physical properties such as mean density, maximum density, and standard deviation along the line of sight. Validated against numerical simulations (FLASH and ENZO), the algorithm achieves high accuracy, with mean density estimates within 0.1 dex and dispersions of 0.2 to 0.3 dex across varied cloud structures. The reconstructed 3D density fields enable the derivation of critical parameters, including volume density, cloud thickness, and density probability distribution functions, offering insights into star formation and ISM evolution. The versatility of the method is demonstrated by applying diverse systems from galaxies (NGC 628) to protostellar disks. The code is available at this https URL.

We report the discovery of a dense molecular ring-like structure in a dense (10$^5$ cm$^{-3}$), cold (pc-scale CO depletion at a factor of 5), and young (10$^4$ year) star-forming region G34.74-0.12, revealed by C$^{18}$O (2-1), HNC (1-0), and N$_2$H$^+$ (1-0) observations with the Atacama Large Millimeter/submillimeter Array (ALMA). The ring-like structure is redshifted with respect to the clump, spanning from $V_{\rm sys,lsr} + 0.9$ to $V_{\rm sys,lsr} + 2.9$ km s$^{-1}$, with a total mass of 109 $M_{\odot}$. It is spatially coincident with 1.3 mm and 3.0 mm dust continuum emission from cores, and several protostellar outflows. However, no free-free emission or H\textsc{ii} region is detected in association with this structure. With a slow expansion speed indicated by the position-velocity diagram, this ring structure differs from rings previously identified in more evolved star-forming regions. Possible explanations for the ring-like structure include a relic wind-blown bubble produced by a deeply embedded young stellar object, a hollow cavity formed by cloud-cloud interactions, a gas ring resulting from a temperature gradient, or a line-of-sight superposition of multiple outflows or dense clouds. This discovery offers a rare observational glimpse into the earliest dynamical processes involved in massive star formation.

We present results on a multi-wavelength analysis of SDSS J025214.67-002813.7, a system which has been previously classified as a binary AGN candidate based on periodic signals detected in the optical light curves. We use available radio-X-ray observations of the system to investigate the true accretion nature. Analyzing new observations from XMM-Newton and NuSTAR, we characterize the X-ray emission and search for evidence of circumbinary accretion. Although the 0.5-10 keV spectrum shows evidence of an additional soft emission component, possibly due to extended emission from hot nuclear gas, we find the spectral shape consistent with a single AGN. Compiling a full multi-wavelength SED, we also search for signs of circumbinary accretion, such as a "notch" in the continuum due to the presence of minidisks. We find that the radio-optical emission agrees with the SED of a standard, radio-quiet, AGN, however there is a large deficit in emission blueward of ~1400 A. Although this deficit in emission can plausibly be attributed to a binary AGN system, we find that the SED of SDSS J0252-0028 is better explained by emission from a reddened, single AGN. However, future studies on the expected hard X-ray emission associated with binary AGN (especially in the unequal-mass regime), will allow for more rigorous analyses of the binary AGN hypothesis.

Gravitational waves (GWs) accompanied by electromagnetic (EM) counterparts provide a novel methodology to measure the Hubble constant ($H_0$), known as bright sirens. However, the rarity of such multi-messenger events limits the precision of the $H_0$ constraint. Recently, the newly-discovered quasi-periodic eruptions (QPEs) show intriguing evidence of a stellar-mass companion captured by a supermassive black hole (SMBH) in an extreme mass-ratio inspiral (EMRI), which is the most promising sources of the space-based GW detectors, such as Laser Interferometer Space Antenna (LISA). Here we study the secular orbital evolution of QPE systems under different theoretical frameworks, and assess their GWs detectability by LISA. We find that LISA can detect the EMRI signals of two or three known QPE systems with a four-year observation. One EMRI event can measure the $H_0$ with an uncertainty of $3\%-8\%$ at a 68.3 percent confidence level, while a few events will reduce it below $2\%$. The EMRI events surpass current bright siren limitations and offer an independent pathway to resolve the Hubble constant tension.

The association between long gamma-ray bursts (LGRBs) and core-collapse supernovae (SNe) has been well established since the discovery of SN 1998bw, which was linked to the low-luminosity LGRB 980425. However, long-term monitoring of several well-localized, low-redshift LGRBs has yielded compelling evidence for the absence of accompanying supernovae. Notably, two long bursts, GRB 211211A and GRB 230307A, show signatures consistent with kilonova emission from compact binary mergers, indicating that at least some long events may originate from progenitors other than core-collapse supernovae. In this study, we conduct a comparative analysis of two samples of LGRBs, i.e., LGRBs with and without supernova associations, to investigate the differences that may reveal intrinsic distinctions in their progenitors. A detailed examination of their prompt emission properties, host galaxy environments, and event rates is performed. While the two samples exhibit considerable overlap in most observed properties, a significant discrepancy in their event rate is revealed. LGRBs without supernova association have an event rate that aligns well with the star formation rate, whereas that of SN-associated LGRBs differs significantly. It indicates that LGRBs without a supernova association may constitute a distinct subclass with intrinsically different progenitors.

Gamma-ray bursts (GRBs) are the most explosive phenomena and can be used to study the expansion of Universe. In this paper, we compile a long GRB sample for the $E_{\mathrm{iso}}$-$E_{\mathrm{p}}$ correlation from Swift and Fermi observations. The sample contains 221 long GRBs with redshifts from 0.03 to 8.20. From the analysis of data in different redshift intervals, we find no statistically significant evidence for the redshift evolution of this correlation. Then we calibrate the correlation in six sub-samples and use the calibrated one to constrain cosmological parameters. Employing a piece-wise approach, we study the redshift evolution of dark energy equation of state (EOS), and find that the EOS tends to be oscillating at low redshift, but consistent with $-1$ at high redshift. It hints a dynamical dark energy at $2\sigma$ confidence level at low redshift.

The disc instability mechanism (DIM) is widely accepted to account for the transient behaviour of dwarf novae (DNe), which experience short outbursts separated by long quiescence. The duty cycle (the ratio between the outburst duration and the recurrence time) determines the amount of accreted mass by the white dwarf (WDs) during outbursts, thus playing an important role in the long-term binary evolution. Employing the code of Modules for Experiments in Stellar Astrophysics, we systemically investigate the influence of the duty cycles on the evolution of DNe and the mass growth of accreting carbon-oxygen (CO) WDs. Our calculations show that, while the DIM can considerably influence the accretion process, efficient WD-mass growth requires a particular range of the duty cycle. For WDs with the initial masses of 0.6, 0.7 and 1.1 $M_\odot$, these duty cycles are 0.006$\,\leq$$d$$\,\leq$0.007, $d$\,=\,0.005 and $d$\,=\,0.003, and the accumulated mass of the WDs can reach 0.1, 0.13 and 0.21 $M_\odot$, respectively. In all of our simulations, no CO WDs can grow their masses to the explosion mass of Type Ia supernovae of about $1.38~M_\odot$. Because of a much short timescale of the outburst state, the final donor-star masses and orbital periods are insensitive to the duty cycles. Therefore, we propose that the DIM in DNe could alleviate the WD mass problem to some extent.

Polarization of electromagnetic waves carries a large amount of information about their astrophysical emitters and the media they passed through, and hence is crucial in various aspects of astronomy. Here we demonstrate an important but long-overlooked depolarization mechanism in astrophysics: when the polarization vector of light travels along a non-planar curve, it experiences an additional rotation, in particular for radio waves. The process leads to depolarization, which we call `geometric' depolarization (GDP). We give a concise theoretical analysis of the GDP effect on the transport of radio waves in a randomly inhomogeneous plasma under the geometrical optics approximation. In the case of isotropic scattering in the coronal plasma, we show that the GDP of the angle-of-arrival of the linearly polarized radio waves propagating through the turbulent plasma cannot be ignored. The GDP effect of linearly polarized radio waves can be generalized to astrophysical phenomena, such as fast radio bursts and stellar radio bursts, etc. Our findings may have a profound impact on the analysis of astrophysical depolarization phenomena.

In this paper, we use two model-independent methods to standardize long gamma-ray bursts (GRBs) using the $E_{\rm iso}-E_{\rm p}$ correlation, where $E_{\rm iso}$ is the isotropic-equivalent gamma-ray energy and $E_{\rm p}$ is the spectral peak energy. We update 42 long GRBs and try to make constraint on cosmological parameters. The full sample contains 151 long GRBs with redshifts from 0.0331 to 8.2. The first method is the simultaneous fitting method. The extrinsic scatter $\sigma_{\rm ext}$ is taken into account and assigned to the parameter $E_{\rm iso}$. The best-fitting values are $a=49.15\pm0.26$, $b=1.42\pm0.11$, $\sigma_{\rm ext}=0.34\pm0.03$ and $\Omega_m=0.79$ in the flat $\Lambda$CDM model. The constraint on $\Omega_m$ is 0.55<\Omega_m<1 at the 1$\sigma$ confidence level. If reduced $\chi^2$ method is used, the best-fit results are $a=48.96\pm0.18$, $b=1.52\pm0.08$ and $\Omega_m=0.50\pm0.12$. The second method is using type Ia supernovae (SNe Ia) to calibrate the $E_{\rm iso}-E_{\rm p}$ correlation. We calibrate 90 high-redshift GRBs in the redshift range from 1.44 to 8.1. The cosmological constraints from these 90 GRBs are $\Omega_m=0.23^{+0.06}_{-0.04}$ for flat $\Lambda$CDM, and $\Omega_m=0.18\pm0.11$ and $\Omega_{\Lambda}=0.46\pm0.51$ for non-flat $\Lambda$CDM. For the combination of GRB and SNe Ia sample, we obtain $\Omega_m=0.271\pm0.019$ and $h=0.701\pm0.002$ for the flat $\Lambda$CDM, and for the non-flat $\Lambda$CDM, the results are $\Omega_m=0.225\pm0.044$, $\Omega_{\Lambda}=0.640\pm0.082$ and $h=0.698\pm0.004$. These results from calibrated GRBs are consistent with that of SNe Ia. Meanwhile, the combined data can improve cosmological constraints significantly, comparing to SNe Ia alone. Our results show that the $E_{\rm iso}-E_{\rm p}$ correlation is promising to probe the high-redshift universe.

The Kepler mission has discovered thousands of exoplanets around various stars with different spectral types (M, K, G, and F) and thus different masses and effective temperatures. Previous studies have shown that the planet occurrence rate, in terms of average number of planets per star, drops with increasing stellar effective temperature (Teff). In this paper, with the final Kepler Data Release (DR25) catalog, we revisit the relation between stellar effective temperature (as well as mass) and planet occurrence, but in terms of the fraction of stars with planets and the number of planets per planetary system (i.e. planet multiplicity). We find that both the fraction of stars with planets and planet multiplicity decrease with increasing stellar temperature and mass. Specifically, about 75% late-type stars (Teff<5000 K) have Kepler-like planets with an average planet multiplicity of ~2.8; while for early-type stars (Teff>6500 K), this fraction and the average multiplicity fall down to ~35% and ~1.8, respectively. The decreasing trend in the fraction of stars with planets is very significant with ${\Delta}$AIC> 30, though the trend in planet multiplicity is somewhat tentative with ${\Delta}$AIC~5. Our results also allow us to derive the dispersion of planetary orbital inclinations in relationship with stellar effective temperature. Interestingly, it is found to be similar to the well-known trend between obliquity and stellar temperature, indicating that the two trends might have a common origin.

The definition of the Galactic coordinate system was announced by the IAU Sub-Commission 33b on behalf of the IAU in 1958. For more than 50 years the definition of the Galactic coordinate system has remained unchanged from this IAU1958 version. On the basis of deep and all-sky catalogs, the position of the Galactic plane can be revised and updated definitions of the Galactic coordinate systems can be proposed. We re-determine the position of the Galactic plane based on modern large catalogs, such as the Two Micron All-Sky Survey (2MASS) and the SPECFIND v2.0. This paper also aims to propose a possible definition of the optimal Galactic coordinate system by adopting the ICRS position of the Sgr A* at the Galactic center. The near-infrared 2MASS point source catalog and the SPECFIND v2.0 catalog of radio continuum spectra are used to determine the mean position of the Galactic plane on the celestial sphere. By fitting the data to an ideal Galactic equator, the parameters defining the Galactic coordinate system are obtained. We find that the obliquity of the Galactic equator on the ICRS principal plane is about $0.4^\circ$ (2MASS) and $0.6^\circ$ (SPECFIND v2.0) larger than the J2000.0 value, which is widely used in coordinate transformations between the equatorial $(\alpha, \delta)$ and the Galactic $(\ell, b)$. Depending on the adopted parameters, data, and methods, the largest difference between the resulting Galactic coordinate systems is several arcminutes. We derive revised transformation matrices and parameters describing the orientation of the Galactic coordinate systems in the ICRS at the 1 milli-arcsecond level to match the precision of modern observations. For practical applications, we propose that a revised definition of the Galactic coordinate system should be required in the near future.

As a new kind of radio transient sources detected at $\sim 1.4$ GHz, fast radio bursts are specially characterized by their short durations and high intensities. Although only ten events are detected so far, fast radio bursts may actually frequently happen at a rate of $\sim 10^{3}$ --- $10^4~\rm{sky}^{-1}~\rm{day}^{-1}$. We suggest that fast radio bursts can be produced by the collisions between neutron stars and asteroids. This model can naturally explain the millisecond duration of fast radio bursts. The energetics and event rate can also be safely accounted for. Fast radio bursts thus may be one side of the multifaces of the neutron star-small body collision events, which are previously expected to lead to X-ray/gamma-ray bursts or glitch/anti-glitches.

The current discrepancy between the Hubble constant $H_0$ derived from the local distance ladder and from the cosmic microwave background is one of the most crucial issues in cosmology, as it possibly indicates unknown systematics or new physics. Here we present a novel non-parametric method to estimate Hubble constant as a function of redshift. We establish independent estimates of the evolution of Hubble constant by diagonalizing the covariance matrix. From type Ia supernovae and the observed Hubble parameter data, a decreasing trend of Hubble constant with a significance of 5.6$\sigma$ confidence level is found. At low redshift, its value is dramatically consistent with that measured from the local distance ladder, and it drops to the value measured from the cosmic microwave background at high redshift. Our results can relieve the Hubble tension, and prefer the late-time solutions of it, especially the new physics.

Kinematics of solar eruptive filaments is one of the important diagnostic parameters for predicting whether solar eruptions would induce geomagnetic storms. Particularly, some geomagnetic storms might be induced by solar filament eruptions originating from unexpected surface source regions because of non-radial ejection. The non-radial ejection of filaments has received widespread attention but remains inconclusive. We select two eruptive filaments, both of which are supported by flux ropes, as indicated by the hot channel structures seen in the 94 Å images and the hook-shaped brightenings where the filament material falls back. We measure the three-dimensional ejection trajectory of the eruptive filaments by integrating the simultaneous observations from SDO and STEREO. Furthermore, we calculate the distribution of the poloidal field along the ejection path and compare it to the ejection acceleration. It is revealed that the reinforcement of the poloidal magnetic field may lead to the suppression of the acceleration, with the acceleration resuming its increase only when the poloidal field diminishes to a certain level. Additionally, we compute the spatial distribution of the poloidal field in various directions and find that the poloidal magnetic field above the filaments is asymmetric. For both investigated events, the filaments appear to eject towards the side where the poloidal magnetic field is weaker, indicating that the eruptive filaments tend to propagate along the side with weaker strapping force. This may provide a new explanation for the inclined ejection of filaments.

Elemental abundances hold important information about the star formation history in the Galactic Center. The thermal X-ray spectra of certain stars can provide a robust probe of elemental abundances, mainly through the presence of K-shell emission lines. In this work, based on deep archival {\it Chandra} observations, we obtain X-ray measurements of five heavy elements (Si, S, Ar, Ca and Fe) for three sources in the Arches cluster, one source in the Quintuplet cluster, as well as a field source known as Edd 1, which are all probable WR stars exhibiting a high quality X-ray spectrum. A two-temperature, non-equilibrium ionization plasma model is employed for the spectral fit, taking into account light element compositions characteristic of WR star winds, which is substantially depleted in hydrogen but enriched in nitrogen and/or carbon. It is found that the Arches and Quintuplet WR stars share similar abundances of Si, S, and Ar, while exhibiting distinct Ca and Fe abundances, which may be understood as due to dust depletion of the latter two elements in Quintuplet. The observed near-solar or sub-solar metallicity of the WR star winds can be naturally understood as the result of nucleosynthesis and internal mixing of the parent star, which have a supersolar initial metallicity as expected for the Galactic center in general. Implications of our findings on the origin of the young star clusters and isolated massive stars in the Galactic center, as well as the elemental composition of the accretion flow onto Sgr A*, are addressed.

The radius valley, a dip in the radius distribution of exoplanets at ~1.9 Earth radii separates compact rocky Super-Earths and Sub-Neptunes with lower density. Various hypotheses have been put forward to explain the radius valley. Characterizing the radius valley morphology and its correlation to stellar properties will provide crucial observation constraints on its origin mechanism and deepen the understanding of planet formation and evolution. In this paper, the third part of the Planets Across the Space and Time (PAST) series, using the LAMOST-Gaia-Kepler catalog, we perform a systematical investigation into how the radius valley morphology varies in the Galactic context, i.e., thin/thick galactic disks, stellar age and metallicity abundance ([Fe/H] and [alpha/Fe]). We find that (1) The valley becomes more prominent with the increase of both age and [Fe/H]. (2) The number ratio of super-Earths to sub-Neptunes monotonically increases with age but decreases with [Fe/H] and [alpha/Fe]. (3) The average radius of planets above the valley (2.1-6 Earth radii) decreases with age but increases with [Fe/H]. (4) In contrast, the average radius of planets below the valley (R < 1.7 Earth radii) is broadly independent on age and metallicity. Our results demonstrate that the valley morphology as well as the whole planetary radius distribution evolves on a long timescale of giga-years, and metallicities (not only Fe but also other metal elements, e.g., Mg, Si, Ca, Ti) play important roles in planet formation and in the long term planetary evolution.

High-mass stars, born in massive dense cores (MDCs), profoundly impact the cosmic ecosystem through feedback processes and metal enrichment, yet little is known about how MDCs assemble and transfer mass across scales to form high-mass young stellar objects (HMYSOs). Using multi-scale (40-2500 au) observations of an MDC hosting an HMYSO, we identify a coherent dynamical structure analogous to barred spiral galaxies: three 20,000 au spiral arms feed a 7,500 au central bar, which channels gas to a 2,000 au pseudodisk. Further accretion proceeds through the inner structures, including a Keplerian disk and an inner disk (100 au), which are thought to be driving a collimated bipolar outflow. This is the first time that these multi-scale structures (spiral arms, bar, streamers, envelope, disk, and outflow) have been simultaneously observed as a physically coherent structure within an MDC. Our discovery suggests that well-organized hierarchical structures play a crucial role during the gas accretion and angular momentum build-up of a massive disk.

As part of our comprehensive, ongoing characterisation of the low-mass end of the main sequence in the Solar neighbourhood, we used the OSIRIS instrument at the 10.4 m Gran Telescopio Canarias to acquire low- and mid-resolution (R$\approx$300 and R$\approx$2500) optical spectroscopy of 53 late-M and L ultracool dwarfs. Most of these objects are known but poorly investigated and lacking complete kinematics. We measured spectral indices, determined spectral types (six of which are new) and inferred effective temperature and surface gravity from BT-Settl synthetic spectra fits for all objects. We were able to measure radial velocities via line centre fitting and cross correlation for 46 objects, 29 of which lacked previous radial velocity measurements. Using these radial velocities in combination with the latest Gaia DR3 data, we also calculated Galactocentric space velocities. From their kinematics, we identified two candidates outside of the thin disc and four in young stellar kinematic groups. Two further ultracool dwarfs are apparently young field objects: 2MASSW J1246467$+$402715 (L4$\beta$), which has a potential, weak lithium absorption line, and G 196$-$3B (L3$\beta$), which was already known as young due to its well-studied primary companion.