Papers for Tuesday, Mar 21 2023

We probe the outer dynamical mass profile of 16 galaxy clusters from the HIghest X-ray FLUx Galaxy Cluster Sample (HIFLUGCS) using galaxy kinematics. Numerically solving the spherical Jeans equation, we parameterize the dynamical mass profile and the galaxy velocity anisotropy profile using two general functions to ensure that our results are not biased towards any specific model. The mass-velocity anisotropy degeneracy is broken by using two "virial shape parameters" that depend on the fourth moment of velocity distribution. The resulting velocity anisotropy estimates consistently show a nearly isotropic distribution in the inner regions, with an increasing radial anisotropy towards large radii. We compare our derived dynamical masses with those calculated from X-ray gas data assuming hydrostatic equilibrium, finding that massive and rich relaxed clusters generally present consistent mass measurements, while unrelaxed or low-richness clusters have systematically larger total mass than hydrostatic mass by an average of 50\%. This might help alleviate current tensions in the measurement of $\sigma_8$, but it also leads to cluster baryon fractions below the cosmic value. Finally, our approach probes accelerations as low as $10^{-11}$ m s$^{-2}$, comparable to the outskirts of individual late-type galaxies. We confirm that galaxy clusters deviate from the radial acceleration relation defined by galaxies.

Launched on December 9, 2021, the Imaging X-ray Polarimetry Explorer (IXPE) is the first imaging polarimeter ever flown, providing sensitivity in the 2--8 keV range, and during the 2-year initial phase of the mission will sample tens of X-ray sources among different source classes. While most of the measurements will be statistics-limited, for some of the brightest objects observed and long integration times, the systematic uncertainties in the detector response (primarily the effective area, the modulation factor and the absolute energy scale) will be important. In this contribution, we describe a framework to propagate on high-level observables (e.g.: spectro-polarimetric fit parameters) the systematic uncertainties connected with the response of the detector, that we estimate from relevant ground calibrations and from observations of celestial point sources.

We present the discovery of the most distant, dynamically relaxed cool core cluster, SPT-CL J2215-3537 (SPT2215) and its central brightest cluster galaxy (BCG) at z=1.16. Using new X-ray observations, we demonstrate that SPT2215 harbors a strong cool core, with a central cooling time of 200 Myr (at 10 kpc) and a maximal intracluster medium cooling rate of 1900+/-400 Msun/yr. This prodigious cooling may be responsible for fueling extended, star-forming filaments observed in Hubble Space Telescope imaging. Based on new spectrophotometric data, we detect bright [O II] emission in the BCG, implying an unobscured star formation rate (SFR) of 320^{+230}_{-140} Msun/yr. The detection of a weak radio source (2.0+/-0.8 mJy at 0.8 GHz) suggests ongoing feedback from an active galactic nucleus (AGN), though the implied jet power is less than half the cooling luminosity of the hot gas, consistent with cooling overpowering heating. The extreme cooling and SFR of SPT2215 is rare among known cool core clusters, and it is even more remarkable that we observe these at such high redshift, when most clusters are still dynamically disturbed. The high mass of this cluster, coupled with the fact that it is dynamically relaxed with a highly-isolated BCG, suggests that it is an exceptionally rare system that must have formed very rapidly in the early Universe. Combined with the high SFR, SPT2215 may be a high-z analog of the Phoenix cluster, potentially providing insight into the limits of AGN feedback and star formation in the most massive galaxies.

Massive binaries are vital sources of various transient processes, including gravitational-wave mergers. However, large uncertainties in the evolution of massive stars, both physical and numerical, present a major challenge to the understanding of their binary evolution. In this paper, we upgrade our interpolation-based stellar evolution code METISSE to include the effects of mass changes, such as binary mass transfer or wind-driven mass loss, not already included within the input stellar tracks. METISSE's implementation of mass loss (applied to tracks without mass loss) shows excellent agreement with the SSE fitting formulae and with detailed MESA tracks, except in cases where the mass transfer is too rapid for the star to maintain equilibrium. We use this updated version of METISSE within the binary population synthesis code BSE to demonstrate the impact of varying stellar evolution parameters, particularly core overshooting, on the evolution of a massive (25M$_\odot$ and 15M$_\odot$) binary system with an orbital period of 1800 days. Depending on the input tracks, we find that the binary system can form a binary black hole or a black hole-neutron star system, with primary(secondary) remnant masses ranging between 4.47(1.36)M$_\odot$ and 12.30(10.89)M$_\odot$, and orbital periods ranging from 6 days to the binary becoming unbound. Extending this analysis to a population of isolated binaries uniformly distributed in mass and orbital period, we show that the input stellar models play an important role in determining which regions of the binary parameter space can produce compact binary mergers, paving the way for current and future gravitational-wave progenitor predictions.

Certain debris disks have non-axisymmetric shapes in scattered light which are unexplained. The appearance of a disk depends on how its constituent Keplerian ellipses are arranged. The more the ellipses align apsidally, the more non-axisymmetric the disk. Apsidal alignment is automatic for fragments released from a catastrophic collision between solid bodies. We synthesize scattered light images, and thermal emission images, of such giant impact debris. Depending on the viewing geometry, and if and how the initial apsidal alignment is perturbed, the remains of a giant impact can appear in scattered light as a one-sided or two-sided "fork", a lopsided "needle", or a set of "double wings". Double wings are difficult to reproduce in other scenarios involving gravitational forcing or gas drag, which do not align orbits as well. We compare our images with observations and offer a scorecard assessing whether the scattered light asymmetries in HD 15115, HD 32297, HD 61005, HD 111520, HD 106906, beta Pic, and AU Mic are best explained by giant impacts, gravitational perturbations, or sculpting by the interstellar medium.

The arrangement of dwarf galaxies in a thin plane surrounding the Milky Way has been thought to contradict the prevailing cosmological model of cold dark matter in the Universe. New work suggests that this arrangement may just be a temporary alignment, bringing our galaxy back into agreement with theoretical expectations once the radial distribution of satellites is taken into account.

We report the confirmation of HIP 67506 C, a new stellar companion to HIP 67506 A. We previously reported a candidate signal at 2$\lambda$/D (240~mas) in L$^{\prime}$ in MagAO/Clio imaging using the binary differential imaging technique. Several additional indirect signals showed that the candidate signal merited follow-up: significant astrometric acceleration in Gaia DR3, Hipparcos-Gaia proper motion anomaly, and overluminosity compared to single main sequence stars. We confirmed the companion, HIP 67506 C, at 0.1" with MagAO-X in April, 2022. We characterized HIP 67506 C MagAO-X photometry and astrometry, and estimated spectral type K7-M2; we also re-evaluated HIP 67506 A in light of the close companion. Additionally we show that a previously identified 9" companion, HIP 67506 B, is a much further distant unassociated background star. We also discuss the utility of indirect signposts in identifying small inner working angle candidate companions.

We present a new high-precision, JWST-based, strong lensing model for the galaxy cluster Abell 2744 at $z=0.3072$. By combining the deep, high-resolution JWST imaging from the GLASS-JWST and UNCOVER programs and a Director's Discretionary Time program, with newly obtained VLT/MUSE data, we identify 32 multiple images from 11 background sources lensed by two external sub-clusters at distances of ~160" from the main cluster. The new MUSE observations enable the first spectroscopic confirmation of a multiple image system in the external clumps. Moreover, the re-analysis of the spectro-photometric archival and JWST data yields 27 additional multiple images in the main cluster. The new lens model is constrained by 149 multiple images ($\sim66\%$ more than in our previous Bergamini et al. 2023 model) covering an extended redshift range between 1.03 and 9.76. The subhalo mass component of the cluster includes 177 member galaxies down to $m_{\rm F160W}=21$, 163 of which are spectroscopically confirmed. Internal velocity dispersions are measured for 85 members. The new lens model is characterized by a remarkably low scatter between predicted and observed positions of the multiple images (0.43"). This precision is unprecedented given the large multiple image sample, the complexity of the cluster mass distribution, and the large modeled area. The improved accuracy and resolution of the cluster total mass distribution provides a robust magnification map over a $\sim\!45$ arcmin$^2$ area, which is critical for inferring the intrinsic physical properties of the highly magnified, high-$z$ sources. The lens model and the new MUSE redshift catalog are released with this publication.

Context: Recently, Doroshenko and collaborators reported a very low-mass compact star, a Central Compact Object named XMMU J173203.3-344518 inside the supernova remnant HESS J1731-347. Its tiny mass is at odds with all calculations of minimum masses of neutron stars generated by iron cores, therefore (and even if not compellingly) it has been suggested to be a {\it strange star}. Besides the mass, a radius and surface temperature have been extracted from data, and the whole body of information should ultimately reveal if this object is truly consistent with an exotic composition. Aims: To understand the status of the compact object XMMU J173203.3-344518 in HESS J1731-347 within the existing models of strange stars, including its prompt formation. Methods: The information obtained on the mass, radius and surface temperature are compared to theoretical calculations performed within usual theoretical models using General Relativity as the assumed theory of gravitation and a handful of cooling scenarios. A qualitative discussion showing the consistency of the strange-matter driven supernova scenario with a low-mass compact star is provided. Results: We found that the object HESS J1731-347 fits within the same quark star models recently employed to explain the masses and radii of the NICER objects PSR J040+6620 and PSR J0030+0451, in which both quantities were simultaneously determined. It is also remarkable to find that a simple cooling scenario devised $30\text{ yr}$ ago with superconducting quarks provides an overall good explanation of the surface temperature. Conclusions: We conclude that XMMU J173203.3-344518 in the remnant HESS J1731-347 fits into a ``strange star'' scenario that is also consistent with heavier compact stars, which can also belong to the same class and constitute an homogeneous type of self-bound objects produced in Nature.

We present a broadband plasma spray anti-reflection (AR) coating technology for millimeter-wave astrophysics experiments with large-format, cryogenic optics. By plasma spraying alumina- and silica-based powders, we have produced coatings of tunable index of refraction and thickness, low loss, and coefficient of thermal expansion matched to alumina substrates. We demonstrate two-layer AR coatings on alumina with reflection below 5% over 82% and 69% fractional bandwidths for 90/150 and 220/280 GHz passband designs, respectively, and band-averaged absorption loss reduced to $\sim$1\% at 100 K for both AR coatings. We describe the design, tolerance, fabrication process, and optical measurements of these AR coatings.

We use $z=1$ mock galaxy catalogues produced with the semi-analytic code GALACTICUS to study the dependence of the non-Gaussian bias parameter $b_\phi$ on the mass assembly history of the host halos. We generate large sets of merger trees and measure the non-Gaussian assembly bias $\Delta b_\phi$ for galaxies selected by color magnitude and emission line luminosities. For galaxies selected by $g-r$ color, we find a large assembly bias consistent with the analysis of Barreira et al. (2020) based on hydro-dynamical simulations of galaxy formation. This effect arises from the fact that a larger value of the normalization amplitude $\sigma_8$ implies a faster mass assembly (at fixed halo mass) and, therefore, older and redder galaxies. On the contrary, for galaxies selected by their H$\alpha$ luminosity, we do not detect a significant assembly bias, at least at $z=1$ and in the halo mass range $3\times10^{10} < M < 10^{12}\ M_\odot$ considered here. This is presumably due to the fact that emission line strengths are mainly sensitive to the instantaneous star formation rate, which appears to depend weakly on $\sigma_8$ at $z=1$. This indicates that the non-Gaussian assembly bias should be less of a concern for future emission line galaxy surveys. We also investigate, for the first time, the sensitivity of the non-Gaussian assembly bias to a change in the parameters of the galaxy formation model that control the AGN and stellar feedback as well as the star formation rate. When these parameters change within a factor of two from their fiducial value, they induce variations up to order unity in the measured $\Delta b_\phi$, but the overall trends with color or luminosity remain the same. However, since these results may be sensitive to the choice of galaxy formation model, it will be prudent to extend this analysis to other semi-analytic models in addition to halo mass and redshift.

Broadband photometric reverberation mapping (PRM) have been investigated for AGNs in recent years, but mostly on accretion disk continuum RM. Due to the small fraction of broad emission lines in the broadband, PRM for emission lines is very challenging. Here we present an ICCF-Cut method for broadband PRM to obtain the H$\alpha$ broad line lag and apply it to four Seyfert 1 galaxies, MCG+08-11-011, NGC 2617, 3C 120 and NGC 5548. All of them have high quality broadband lightcurves with daily/sub-daily cadence, which enables us to extract H$\alpha$ lightcurves from the line band by subtracting the contributions from the continuum and host galaxy. Their extracted H$\alpha$ lightcurves are compared with the lagged continuum band lightcurves, as well as the lagged H$\beta$ lightcurves obtained by spectroscopic RM (SRM) at the same epochs. The consistency of these lightcurves and the comparison with the SRM H$\beta$ lags provide supports to the H$\alpha$ lags of these AGNs, in a range from 9 to 19 days, obtained by the ICCF-Cut, JAVELIN and $\chi^2$ methods. The simulations to evaluate the reliability of H$\alpha$ lags and the comparisons between SRM H$\beta$ and PRM H$\alpha$ lags indicate that the consistency of the ICCF-Cut, JAVELIN and $\chi^2$ results can ensure the reliability of the derived H$\alpha$ lags. These methods may be used to estimate the broad line region sizes and black hole masses of a large sample of AGNs in the large multi-epoch high cadence photometric surveys such as LSST in the future.

Studying open clusters (OCs) is essential for a comprehensive understanding of the structure and evolution of the Milky Way. Many previous studies have systematically searched for OCs near the solar system within 1.2 kpc or 20 degrees of galactic latitude. However, few studies searched for OCs at higher galactic latitudes and deeper distances. In this study, based on a hybrid unsupervised clustering algorithm (Friends-of-Friends and pyUPMASK) and a binary classification algorithm (Random Forest), we extended the search region (i.e., galactic latitude |b|>=20 degrees) and performed a fine-grained blind search of Galactic clusters in Gaia DR3. After cross-matching, the newly discovered cluster candidates are fitted using isochrone fitting to estimate the main physical parameters (age and metallicity) of these clusters. These cluster candidates were then checked using manual visual inspection. Their statistical properties were compared with previously exposed cluster catalogs as well. In the end, we found 1,179 new clusters with considerable confidence within 5kpc.

Computational Fluid Dynamics (CFD)-based global solar coronal simulations are slowly making their way into the space weather modeling toolchains to replace the semi-empirical methods such as the Wang-Sheeley-Arge (WSA) model. However, since they are based on CFD, if the assumptions in them are too strong, these codes might experience issues with convergence and unphysical solutions. Particularly the magnetograms corresponding to solar maxima can pose problems as they contain active regions with strong magnetic fields, resulting in large gradients. Combined with the approximate way in which the inner boundary is often treated, this can lead to non-physical features or even a complete divergence of the simulation in these cases. Here, we show some of the possible approaches to handle this inner boundary in our global coronal model COolfluid COrona uNstrUcTured (COCONUT) in a way that improves both convergence and accuracy. Since we know that prescribing the photospheric magnetic field for a region that represents the lower corona is not entirely physical, first, we look at the ways in which we can adjust the input magnetograms to remove the highest magnitudes and gradients. Secondly, since in the default setup we also assume a constant density, here we experiment with changing these values locally and globally to see the effect on the results. We conclude, through comparison with observations and convergence analysis, that modifying the density locally in active regions is the best way to improve the performance both in terms of convergence and physical accuracy from the tested approaches.

We propose a new evolutionary process of protoplanetary disks "co-evolution of dust grains and protoplanetary disks", revealed by dust-gas two-fluid non-ideal magnetohydrodynamics simulations considering the growth of dust and associated changes in magnetic resistivity. We found that the dust growth significantly affects disk evolution by changing the coupling between the gas and magnetic field. Moreover, once the dust grains sufficiently grow, the physical quantities (e.g., density and magnetic field) of the disk are well described by nontrivial power laws, regardless of the details of the dust model. In this disk structure, the radial profile of density is steeper and the disk mass is smaller than those of the model ignoring dust growth and they are more consistent with the disk observations. We analytically derive these power laws from the basic equations of non-ideal magnetohydrodynamics. The analytical power laws are determined only by observable physical quantities, e.g., central stellar mass and mass accretion rate, and do not include difficult-to-determine parameters e.g., viscous parameter $\alpha$. Therefore, they are applicable to various stages of disk evolution. We believe that the disk structure provides a new basis for future studies on star and planet formation.

Shock wave heating is a leading candidate for the mechanisms of chondrule formation. This mechanism forms chondrules when the shock velocity is in a certain range. If the shock velocity is lower than this range, dust particles smaller than chondrule precursors melt, while chondrule precursors do not. We focus on the low-velocity shock waves as the igneous rim accretion events. Using a semi-analytical treatment of the shock-wave heating model, we found that the accretion of molten dust particles occurs when they are supercooling. The accreted igneous rims have two layers, which are the layers of the accreted supercooled droplets and crystallized dust particles. We suggest that chondrules experience multiple rim-forming shock events.

With the aim of improving our knowledge about their nature, we conduct a comparative study on a sample of long-rising Type II supernovae (SNe) resembling SN 1987A. To do so, we deduce various scaling relations from different analytic models of H-rich SNe, discussing their robustness and feasibility. Then we use the best relations in terms of accuracy to infer the SN progenitor's physical properties at the explosion for the selected sample of SN 1987A-like objects, deriving energies of $\sim 0.5$-$15$ foe, radii of $\sim 0.2$-$100 \times 10^{12}$ cm, and ejected masses of $\sim 15$-$55$\msun. Although the sample may be too small to draw any final conclusion, these results suggest that (a) SN 1987A-like objects have parameters at explosion covering a wide range of values; (b) the main parameter determining their distribution is the explosion energy; (c) a high-mass ($\gtrsim 30$\,\Msun), high-energy ($\gtrsim 10$\,foe) tail of events, linked to extended progenitors with radii at explosion $\sim 10^{13}$-$10^{14}$\,cm, challenge standard theories of neutrino-driven core-collapse and stellar evolution. We also find a correlation between the amount of $^{56}$Ni in the ejecta of the SN 1987A-like objects and the spectrophotometric features of the SN at maximum, that may represent a tool for estimating the amount of $^{56}$Ni in the SN ejecta whitout having information on the tail luminosity.

In this work, we make use of a supervised machine learning algorithm based on Logistic Regression (LR) to select TeV blazar candidates from the 4FGL-DR2 / 4LAC-DR2, 3FHL, 3HSP, and 2BIGB catalogs. LR constructs a hyperplane based on a selection of optimal parameters, named features, and hyper-parameters whose values control the learning process and determine the values of features that a learning algorithm ends up learning, to discriminate TeV blazars from non-TeV blazars. In addition, it gives the probability (or logistic) that a source may be considered as a TeV blazar candidate. Non-TeV blazars with logistics greater than 80% are considered high-confidence TeV candidates. Using this technique, we identify 40 high-confidence TeV candidates from the 4FGL-DR2 / 4LAC-DR2 blazars and we build the feature hyper-plane to distinguish TeV and non-TeV blazars. We also calculate the hyper-planes for the 3FHL, 3HSP, and 2BIGB. Finally, we construct the broadband spectral energy distributions (SED) for the 40 candidates, testing for their detectability with various instruments. We find that 7 of them are likely to be detected by existing or upcoming IACT observatories, while 1 could be observed with EAS particle detector arrays.

GRB221009A is the brightest gamma-ray burst ever detected. To probe the very-high-energy (VHE, $>$\!100 GeV) emission, the High Energy Stereoscopic System (H.E.S.S.) began observations 53 hours after the triggering event, when the brightness of the moonlight no longer precluded observations. We derive differential and integral upper limits using H.E.S.S. data from the third, fourth, and ninth nights after the initial GRB detection, after applying atmospheric corrections. The combined observations yield an integral energy flux upper limit of $\Phi_\mathrm{UL}^{95\%} = 9.7 \times 10^{-12}~\mathrm{erg\,cm^{-2}\,s^{-1}}$ above $E_\mathrm{thr} = 650$ GeV. The constraints derived from the H.E.S.S. observations complement the available multiwavelength data. The radio to X-ray data are consistent with synchrotron emission from a single electron population, with the peak in the SED occurring above the X-ray band. Compared to the VHE-bright GRB190829A, the upper limits for GRB221009A imply a smaller gamma-ray to X-ray flux ratio in the afterglow. Even in the absence of a detection, the H.E.S.S. upper limits thus contribute to the multiwavelength picture of GRB221009A, effectively ruling out an IC dominated scenario.

We present a high-density relativistic reflection analysis of 21 spectra of six black hole X-ray binaries in the hard state with data from \textit{NuSTAR} and \textit{Swift}. We find that 76\% of the observations in our sample require a disk density higher than the 10$^{15}$~cm$^{-3}$ assumed in the previous reflection analysis. Compared with the measurements from active galactic nuclei, stellar mass black holes have higher disk densities. Our fits indicate that the inner disk radius is close to the innermost stable circular orbit in the hard state. The coronal temperatures are significantly lower than the prediction of a purely thermal plasma, which can be explained with a hybrid plasma model. If the disk density is fixed at 10$^{15}$~cm$^{-3}$, the disk ionization parameter would be overestimated while the inner disk radius is unaffected.

Results from Juno's microwave radiometer indicate non-uniform mixing of ammonia vapor in Jupiter's atmosphere down to tens of bars, far beneath the cloud level. Helioseismic observations suggest solar convection may require narrow, concentrated downdrafts called entropy rain to accommodate the Sun's luminosity. Both observations suggest some mechanism of non-local convective transport. We seek to predict the depth that a concentrated density anomaly can reach before efficiently mixing with its environment in bottomless atmospheres. We modify classic self-similar analytical models of entraining thermals to account for the compressibility of an abyssal atmosphere. We compare these models to the output of high resolution three dimensional fluid dynamical simulations to more accurately model the chaotic influence of turbulence. We find that localized density anomalies propagate down to ~3-8 times their initial size without substantially mixing with their environment. Our analytic model accurately predicts the initial flow, but the self-similarity assumption breaks down after the flow becomes unstable at a characteristic penetration depth. In the context of Jupiter, our findings suggest that precipitation concentrated into localized downdrafts of size ~20km can coherently penetrate to on the order of a hundred kilometers (tens of bars) beneath its initial vaporization level without mixing with its environment. This finding is consistent with expected convective storm length-scales, and Juno MWR measurements of ammonia depletion. Compositional gradients of volatiles beneath their cloud levels may be common on stormy giant planets. In the context of the Sun, we find that turbulent downdrafts in abyssal atmospheres cannot maintain their coherence through the Sun's convective layer, a potential challenge for the entropy rain hypothesis.

Lorentz invariance is one of the fundamental tenets of Special Relativity, and has been extensively tested with laboratory and astrophysical observations. However, many quantum gravity models and theories beyond the Standard Model of Particle Physics predict a violation of Lorentz invariance at energies close to Planck scale. This article reviews observational and experimental tests of Lorentz invariance violation (LIV) with photons, neutrinos and gravitational waves. Most astrophysical tests of LIV using photons are based on searching for a correlation of the spectral lag data with redshift and energy. These have been primarily carried out using compact objects such as pulsars, Active Galactic Nuclei (AGN), and Gamma-ray bursts (GRB). There have also been some claims for LIV from some of these spectral lag observations with GRBs, which however are in conflict with the most stringent limits obtained from other LIV searches. Searches have also been carried out using polarization measurements from GRBs and AGNs. For neutrinos, tests have been made using both astrophysical observations at MeV energies (from SN 1987A) as well as in the TeV-PeV energy range based on IceCube observations, atmospheric neutrinos, and long-baseline neutrino oscillation experiments. Cosmological tests of LIV entail looking for a constancy of the speed of light as a function of redshift using multiple observational probes, as well as looking for birefringence in Cosmic Microwave background observations. This article will review all of these aforementioned observational tests of LIV, including results which are in conflict with each other.

The {\it Linear Point} (LP), defined as the midpoint between the BAO peak and the associated left dip of the two-point correlation function (2PCF), $\xi(s)$, is proposed as a new standard ruler which is insensitive to nonlinear effects. In this paper, we use a Bayesian sampler to measure the LP and estimate the corresponding statistical uncertainty, and then perform cosmological parameter constraints with LP measurements. Using the Patchy mock catalogues, we find that the measured LPs are consistent with theoretical predictions at 0.6 per cent level. We find constraints with midpoints identified from the rescaled 2PCF ($s^2 \xi$) more robust than those from the traditional LP based on $\xi$, as the BAO peak is not always prominent when scanning the cosmological parameter space, with the cost of 2--4 per cent increase of statistical uncertainty. This problem can also be solved by an additional dataset that provides strong parameter constraints. Measuring LP from the reconstructed data slightly increases the systematic error but significantly reduces the statistical error, resulting in more accurate measurements. The 1$\,\sigma$ confidence interval of distance scale constraints from LP measurements are 20--30 per cent larger than those of the corresponding BAO measurements. For the reconstructed SDSS DR12 data, the constraints on $H_0$ and $\Omega_{\rm m}$ in a flat-$\Lambda$CDM framework with the LP are generally consistent with those from BAO. When combined with Planck cosmic microwave background data, we obtain $H_0=68.02_{-0.37}^{+0.36}$ ${\rm km}\,{\rm s}^{-1}\,{\rm Mpc}^{-1}$ and $\Omega_{\rm m}=0.3055_{-0.0048}^{+0.0049}$ with the LP.

Third-generation (3G) gravitational wave detectors, in particular Einstein Telescope (ET) and Cosmic Explorer (CE), will explore unprecedented cosmic volumes in search for compact binary mergers, providing us with tens of thousands of detections per year. In this study, we simulate and employ binary black holes detected by 3G interferometers as dark sirens, to extract and infer cosmological parameters by cross-matching gravitational wave data with electromagnetic information retrieved from a simulated galaxy catalog. Considering a standard $\Lambda$CDM model, we apply a suitable Bayesian framework to obtain joint posterior distributions for the Hubble constant $H_0$ and the matter energy density parameter $\Omega_m$ in different scenarios. Assuming a galaxy catalog complete up to $z=1$ and dark sirens detected with a network signal-to-noise ratio greater than 300, we show that a network made of ET and two CEs can constrain $H_0$ ($\Omega_m$) to a promising $0.7\%$ ($9.0\%$) at $90\%$ confidence interval within one year of continuous observations. Additionally, we find that most of the information on $H_0$ is contained in local, single-host dark sirens, and that dark sirens at $z>1$ do not substantially improve these estimates. Our results imply that a sub-percent measure of $H_0$ can confidently be attained by a network of 3G detectors, highlighting the need for characterising all systematic effects to a higher accuracy.

According to our understanding of stellar evolution, early-type stars have radiative envelopes and convective cores due to a steep temperature gradient produced by the CNO cycle. Some of these stars (mainly, the subclasses Ap and Bp) have strong magnetic fields, enough to be directly observed using the Zeeman effect. Here, we present 3D magnetohydrodynamic simulations of an $2 ~M_{\odot}$ A-type star using the star-in-a-box model. Our goal is to explore if the modeled star is able to maintain a magnetic field as strong as the observed ones, via a dynamo driven by its convective core, or via maintaining a stable fossil field configuration coming from its early evolutionary stages, using different rotation rates. We created two models, a partially radiative and a fully radiative one, which are determined by the value of the heat conductivity. Our model is able to explore both scenarios, including convection-driven dynamos.

We consider inflation driven by an axion-like particle coupled to an SU(2) gauge sector via a Chern-Simons term. Known as chromo-natural inflation, this scenario is in tension with CMB observations. In order to remedy this fact and preserve both the symmetries and the intriguing gravitational wave phenomenology exhibited by the model, we explore the non-minimal coupling of the axion-inflaton to the Einstein tensor. We identify regions of parameter space corresponding to a viable cosmology at CMB scales. We also highlight the possibility of a large and chiral gravitational wave signal at small scales. This is of particular interest for gravitational wave interferometers.

Context: Coronal mass ejections (CMEs) are eruptions of plasma from the Sun that travel through interplanetary space and may encounter Earth. CMEs often enclose a magnetic flux rope (MFR), the orientation of which largely determines the CME's geoeffectiveness. Current operational CME models do not model MFRs, but a number of research ones do, including the Open Solar Physics Rapid Ensemble Information (OSPREI) model. Aims: We report the sensitivity of OSPREI to a range of user-selected photospheric and coronal conditions. Methods: We model four separate CMEs observed in situ by Parker Solar Probe (PSP). We vary the input photospheric conditions using four input magnetograms (HMI Synchronic, HMI Synoptic, GONG Synoptic Zero-Point Corrected, and GONG ADAPT). To vary the coronal field reconstruction, we employ the Potential-Field Source-Surface (PFSS) model and we vary its source-surface height in the range 1.5--3.0 R$_{\odot}$ with 0.1 R$_{\odot}$ increments. Results: We find that both the input magnetogram and PFSS source surface often affect the evolution of the CME as it propagates through the Sun's corona into interplanetary space, and therefore the accuracy of the MFR prediction compared to in-situ data at PSP. There is no obvious best combination of input magnetogram and PFSS source surface height. Conclusions: The OSPREI model is moderately sensitive to the input photospheric and coronal conditions. Based on where the source region of the CME is located on the Sun, there may be best practices when selecting an input magnetogram to use.

The gamma-ray burst GRB 221009A, known as the ``brightest-of-all-time" (BOAT), is the closest energetic burst detected so far, with an energy of $E_{\gamma,\rm iso} \sim 10^{55}$ ergs. This study aims to assess its compatibility with known GRB energy and luminosity distributions. Our analysis indicates that the energy/luminosity function of GRBs is consistent across various redshift intervals, and that the inclusion of GRB 221009A does not significantly impact the function at low redshifts. Additionally, our evaluation of the best-fitting result of the entire GRB sample suggests that the expected number of GRBs with energy greater than $10^{55}$ ergs at a low redshift is 0.2, so that the emergence of GRB 221009A is consistent with expected energy/luminosity functions within $\sim 2\sigma$ Poisson fluctuation error, still adhering to the principles of small number statistics. Furthermore, we find that GRB 221009A and other energetic bursts, defined as $E_{\gamma,\rm iso} \gtrsim10^{54} {\rm ergs}$, exhibit no significant differences in terms of distributions of $T_{90}$, minimum timescale, Amati relation, $E_{\rm \gamma,iso}$-$E_{\rm X,iso}$ relation, $L_{\gamma,\rm iso}-\Gamma_0$ relation, $E_{\gamma,\rm iso}-\Gamma_0$ relation, $L_{\gamma,\rm iso}-E_{\rm p,i}-\Gamma_0$ relation, and host galaxy properties, compared to normal long GRBs. This suggests that energetic GRBs (including GRB 221009A) and other long GRBs likely have similar progenitor systems and undergo similar energy dissipation and radiation processes. The generation of energetic GRBs may be due to more extreme central engine properties or, more likely, a rarer viewing configuration of a quasi-universal structured jet.

Winds can be launched in tidal disruption event (TDE). It has been proposed that the winds can interact with the cloud surrounding the black hole, produce bow shocks, accelerate electrons, and produce radio emission. We restudy the wind-cloud interaction model. We employ the properties of winds found by the radiation hydrodynamic simulations of super-Eddington circularized accretion flow in TDEs. We can calculate the peak radio emission frequency, the luminosity at the peak frequency, and their time-evolution based on the TDEs wind-cloud interaction model. We find that the model predicted peak radio emission frequency, the luminosity at peak frequency, and their time evolution can be well consistent with those in TDEs AT2019dsg and ASASSN-14li. This indicates that in these two radio TDEs, the wind-cloud interaction mechanism may be responsible for the radio emission.

We develop a self-consistent and accurate halo model by partitioning matter according to the depletion radii of haloes. Unlike conventional models that define haloes with the virial radius while relying on a separate exclusion radius or ad-hoc fixes to account for halo exclusion, our model distributes mass across all scales self-consistently. Using a cosmological simulation, we show that our halo definition leads to very simple and intuitive model components, with the one-halo term given by the Einasto profile with no truncation needed, and the halo-halo correlation function following a universal power-law form down to the halo boundary. The universal halo-halo correlation also allows us to easily model the distribution of unresolved haloes as well as diffuse matter. Convolving the halo profile with the halo-halo correlation function, we obtain a complete description of the halo-matter correlation across all scales, which self-consistently accounts for halo exclusion on the transition scale. Mass conservation is explicitly maintained in our model, and the scale dependence of the classical halo bias is easily reproduced. Our model can successfully reconstruct the halo-matter correlation function with percent level accuracy for halo virial masses in the range of $10^{11.5}h^{-1}{\rm M}_{\odot}<M_{\rm vir}<10^{15.35}h^{-1}{\rm M}_{\odot}$ at $z=0$, and covers the radial range of $0.01h^{-1}{\rm Mpc}<r<20h^{-1}{\rm Mpc}$. We also show that our model profile can accurately predict the characteristic depletion radius at the minimum bias and the splash-back radius at the steepest density slope locations.

We investigate the build-up of the halo profile out to large scale in a cosmological simulation, focusing on the roles played by the recently proposed depletion radii. We explicitly show that halo growth is accompanied by the depletion of the environment, with the inner depletion radius demarcating the two. This evolution process is also observed via the formation of a trough in the bias profile, with the two depletion radii identifying key scales in the evolution. The ratio between the inner depletion radius and the virial radius is approximately a constant factor of 2 across redshifts and halo masses. The ratio between their enclosed densities is also close to a constant of 0.18. These simple scaling relations reflect the largely universal scaled mass profile on these scales, which only evolves weakly with redshift. The overall picture of the boundary evolution can be broadly divided into three stages according to the maturity of the depletion process, with cluster halos lagging behind low mass ones in the evolution. We also show that the traditional slow and fast accretion dichotomy of halo growth can be identified as accelerated and decelerated depletion phases respectively.

As well as the galaxy number density and peculiar velocity, the galaxy intrinsic alignment can be used to test the cosmic isotropy. We study distinctive impacts of the isotropy breaking on the configuration-space two-point correlation functions (2PCFs) composed of the spin-2 galaxy ellipticity field. For this purpose, we build a methodology to efficiently compute general types of the isotropy-violating 2PCFs by generalizing the polypolar spherical harmonic decomposition approach to the spin-weighted version. As a demonstration, we analyze the 2PCFs when the matter power spectrum has a well-known $g_*$-type isotropy-breaking term (induced by, e.g., dark vector fields). We then confirm that some anisotropic distortions indeed appear in the 2PCFs and their shapes rely on a preferred direction causing the isotropy violation, $\hat{d}$. Such a feature can be a distinctive indicator for testing the cosmic isotropy. Comparing the isotropy-violating 2PCFs computed with and without the plane parallel (PP) approximation, we find that, depending on $\hat{d}$, the PP approximation is no longer valid when an opening angle between the directions towards target galaxies is ${\cal O}(1^\circ)$ for the density-ellipticity and velocity-ellipticity cross correlations and around $10^\circ$ for the ellipticity auto correlation. This suggests that an accurate test for the cosmic isotropy requires the formulation of the 2PCF without relying on the PP approximation.

The spectroheliograph is a spectroscopic instrument designed to produce monochromatic images of the photosphere (the visible layer) and the chromosphere of the Sun. It was invented at the same time (1892), but independently, by Hale in the USA and Deslandres in France and was dedicated to long-term surveys of the solar cycles. For that purpose, systematic observations of the CaII K and H$\alpha$ lines started in Meudon observatory in 1908 and continue today, so that a huge collection of more than 100000 spectroheliograms, spanning 115 years of solar activity, was recorded. We present in this paper the optical characteristics and the capabilities of the successive versions of the instrument, from 1908 to now.

Context. A solar eruption on 17 April 2021 produced a widespread Solar Energetic Particle (SEP) event that was observed by five longitudinally well-separated observers in the inner heliosphere at heliocentric distances of 0.42 to 1 au: BepiColombo, Parker Solar Probe, Solar Orbiter, STEREO A, and near-Earth spacecraft. The event produced relativistic electrons and protons. It was associated with a long-lasting solar hard X-ray flare and a medium fast Coronal Mass Ejection (CME) with a speed of 880 km/s driving a shock, an EUV wave as well as long-lasting radio burst activity showing four distinct type III burst. Methods. A multi-spacecraft analysis of remote-sensing and in-situ observations is applied to attribute the SEP observations at the different locations to the various potential source regions at the Sun. An ENLIL simulation is used to characterize the interplanetary state and its role for the energetic particle transport. The magnetic connection between each spacecraft and the Sun is determined. Based on a reconstruction of the coronal shock front we determine the times when the shock establishes magnetic connections with the different observers. Radio observations are used to characterize the directivity of the four main injection episodes, which are then employed in a 2D SEP transport simulation. Results. Timing analysis of the inferred SEP solar injection suggests different source processes being important for the electron and the proton event. Comparison among the characteristics and timing of the potential particle sources, such as the CME-driven shock or the flare, suggests a stronger shock contribution for the proton event and a more likely flare-related source of the electron event. Conclusions. We find that in this event an important ingredient for the wide SEP spread was the wide longitudinal range of about 110 degrees covered by distinct SEP injections.

We have started a systematic survey of molecular clumps with infall motions to study the very early phase of star formation. Our first step is to utilize the data products by MWISP to make an unbiased survey for blue asymmetric line profiles of CO isotopical molecules. Within a total area of $\sim$ 2400 square degrees nearby the Galactic plane, we have found 3533 candidates showing blue-profiles, in which 3329 are selected from the $^{12}$CO&$^{13}$CO pair and 204 are from the $^{13}$CO&C$^{18}$O pair. Exploration of the parametric spaces suggests our samples are in the cold phase with relatively high column densities ready for star formation. Analysis of the spatial distribution of our samples suggests that they exist virtually in all major components of the Galaxy. The vertical distribution suggest that the sources are located mainly in the thick disk of $\sim$ 85 parsec, but still a small part are located far beyond Galactic midplane. Our follow-up observation indicates that these candidates are a good sample to start a search for infall motions, and to study the condition of very early phase of star formation.

The origin of the slow solar wind is still an open issue. One possibility that has been suggested is that upflows at the edge of an active region can contribute to the slow solar wind. We aim to explain how the plasma upflows are generated, which mechanisms are responsible for them, and what the upflow region topology looks like. We investigated an upflow region using imaging data with the unprecedented temporal (3s) and spatial (2 pixels = 236km) resolution that were obtained on 30 March 2022 with the 174{\AA} of the Extreme-Ultraviolet Imager (EUI)/High Resolution Imager (HRI) on board Solar Orbiter. During this time, the EUI and Earth-orbiting satellites (Solar Dynamics Observatory, Hinode, and the Interface Region Imaging Spectrograph, IRIS) were located in quadrature (92 degrees), which provides a stereoscopic view with high resolution. We used the Hinode/EIS (Fe XII) spectroscopic data to find coronal upflow regions in the active region. The IRIS slit-jaw imager provides a high-resolution view of the transition region and chromosphere. For the first time, we have data that provide a quadrature view of a coronal upflow region with high spatial resolution. We found extended loops rooted in a coronal upflow region. Plasma upflows at the footpoints of extended loops determined spectroscopically through the Doppler shift are similar to the apparent upward motions seen through imaging in quadrature. The dynamics of small-scale structures in the upflow region can be used to identify two mechanisms of the plasma upflow: Mechanism I is reconnection of the hot coronal loops with open magnetic field lines in the solar corona, and mechanism II is reconnection of the small chromospheric loops with open magnetic field lines in the chromosphere or transition region. We identified the locations in which mechanisms I and II work.

We investigate the mass-metallicity relationship of star forming galaxies by analysing the absorption line spectra of $\sim$200,000 galaxies in the Sloan Digital Sky Survey. The galaxy spectra are stacked in bins of stellar mass and a population synthesis technique is applied yielding metallicities, ages and star formation history of the young and old stellar population together with interstellar reddening and extinction. We adopt different lengths of the initial starbursts and different initial mass functions for the calculation of model spectra of the single stellar populations contributing to the total integrated spectrum. We also allow for deviations of the ratio of extinction to reddening RV from 3.1 and determine the value from the spectral fit. We find that burst length and RV have a significant influence on the determination of metallicities whereas the effect of the initial mass function is small. RV values are larger than 3.1. The metallicities of the young stellar population agree with extragalactic spectroscopic studies of individual massive supergiant stars and are significantly higher than those of the older stellar population. This confirms galaxy evolution models where metallicity depends on the ratio of gas to stellar mass and where this ratio decreases with time. Star formation history is found to depend on galaxy stellar mass. Massive galaxies are dominated by stars formed at early times.

Nebular emission lines are powerful diagnostics for the physical processes at play in galaxy formation and evolution. Moreover, emission-line galaxies (ELGs) are one of the main targets of current and forthcoming spectroscopic cosmological surveys. We investigate the contributions to the line luminosity functions (LFs) of different galaxy populations in the local Universe, providing a benchmark for future surveys of earlier cosmic epochs. The large statistics of the observations from the SDSS DR7 Main galaxy sample and the MPA-JHU spectral catalogue enables us to precisely measure the H$\alpha$, H$\beta$, [OII], [OIII], [NII], and [SII] emission-line LFs over ~2.4 Gyrs in the low-z Universe, 0.02<z<0.22. We present a generalised 1/Vmax LF estimator capable of simultaneously correcting for spectroscopic, r-band magnitude, and emission-line incompleteness. We study the contribution to the LF of different types of ELGs classified using two methods: (i) the specific star formation rate (sSFR), and (ii) the line ratios on the Baldwin-Phillips-Terlevich (BPT) and WHAN (i.e. H$\alpha$ equivalent width versus the [NII]/H$\alpha$ ratio) diagrams. The ELGs in our sample are mostly star forming, with 83.6 per cent having sSFR>10$^{-11}$/yr. When classifying ELGs using the BPT+WHAN diagrams, we find that 63 per cent are star forming, 1.5 are passively evolving, and 2.9 have nuclear activity (Seyfert). The rest are LINERs and composite ELGs. We find that a Saunders function is the most appropriate to describe all the emission-line LFs. They are dominated by star-forming regions, except for the bright end of the [OIII] and [NII] LFs (i.e. L[NII]>10$^{42}$ erg/s, L[OIII]>10$^{43}$ erg/s), where the contribution of Seyfert galaxies is not negligible. Besides the star-forming population, composite galaxies and LINERs are the ones that contribute the most to the ELG numbers at L<10$^{41}$ erg/s.

The Quasi-thermal noise (QTN) technique is a reliable tool to yield accurate measurements of the electron parameters in the solar wind. We apply this method on Parker Solar Probe (PSP) observations to derive the total electron temperature ($T_e$) from the linear fit of the high-frequency part of the QTN spectra acquired by the RFS/FIELDS instrument, and present a combination of 12-day period of observations around each perihelion from Encounter One (E01) to Ten (E10) (with E08 not included) with the heliocentric distance varying from about 13 to 60 solar radii ($R_\odot{}$). We find that the total electron temperature decreases with the distance as $\sim$$R^{-0.66}$, which is much slower than adiabatic. The extrapolated $T_e$ based on PSP observations is consistent with the exospheric solar wind model prediction at $\sim$10 $R_\odot{}$, Helios observations at $\sim$0.3 AU and Wind observations at 1 AU. Also, $T_e$, extrapolated back to 10 $R_\odot{}$, is almost the same as the strahl electron temperature $T_s$ (measured by SPAN-E) which is considered to be closely related to or even almost equal to the coronal electron temperature. Furthermore, the radial $T_e$ profiles in the slower solar wind (or flux tube with larger mass flux) are steeper than those in the faster solar wind (or flux tube with smaller mass flux). More pronounced anticorrelated $V_p$-$T_e$ is observed when the solar wind is slower and closer to the Sun.

We present integral field, far-infrared (FIR) spectroscopy of Mrk 54, a local Lyman Continuum Emitter (LCE), obtained with FIFI-LS on the Stratospheric Observatory for Infrared Astronomy. This is only the second time, after Haro 11, that [C II] 158 $\mu$m and [O III] 88 $\mu$m spectroscopy of the known LCEs have been obtained. We find that Mrk 54 has a strong [C II] emission that accounts for $\sim1$% of the total FIR luminosity, whereas it has only moderate [O III] emission, resulting in the low [O III]/[C II] luminosity ratio of $0.22\pm0.06$. In order to investigate whether [O III]/[C II] is a useful tracer of $f_{\rm esc}$ (LyC escape fraction), we examine the correlations of [O III]/[C II] and (i) the optical line ratio of $\rm O_{32} \equiv$ [O III] 5007 \AA/[O II] 3727 \AA, (ii) specific star formation rate, (iii) [O III] 88 $\mu$m/[O I] 63 $\mu$m ratio, (iv) gas phase metallicity, and (v) dust temperature based on a combined sample of Mrk 54 and the literature data from the Herschel Dwarf Galaxy Survey and the LITTLE THINGS Survey. We find that galaxies with high [O III]/[C II] luminosity ratios could be the result of high ionization (traced by $\rm O_{32}$), bursty star formation, high ionized-to-neutral gas volume filling factors (traced by [O III] 88 $\mu$m/[O I] 63 $\mu$m), and low gas-phase metallicities, which is in agreement with theoretical predictions. We present an empirical relation between the [O III]/[C II] ratio and $f_{\rm esc}$ based on the combination of the [O III]/[C II] and $\rm O_{32}$ correlation, and the known relation between $\rm O_{32}$ and $f_{\rm esc}$. The relation implies that high-redshift galaxies with high [O III]/[C II] ratios revealed by ALMA may have $f_{\rm esc}\gtrsim0.1$, significantly contributing to the cosmic reionization.

GRB~230307A is one of the birghtest gamma-ray bursts detected so far. Its prompt emission has been analyzed by the high-energy detection data comprehensive analysis tool ({\tt HEtools}) developed by us. With the excellent observation of GRB~230307A by Fermi-GBM in keV-MeV energy range, we can reveal the details about the evolution of this relation. As found in the high-time-resolution spectral analysis, the early low-energy spectral indices ($\alpha$) of this burst exceed the limit of synchrotron radiation ($\alpha=-2/3$), and gradually decreases with the energy flux ($F$). A tight $E_{\rm p}\propto F^{0.44}$ correlation anyhow holds within the whole duration of the burst, where $E_{\rm p}$ is the spectral peak energy. Such evolution pattern of $\alpha$ and $E_{\rm p}$ with intensity is called ``double tracking". In addition, for the relation between the $\alpha$ and $F$, we find a log Bayes factor $\sim$ 160 in favor of a smoothly broken power-law function over a linear function in log-linear space. We call this particular $\alpha-F$ relation as broken ``$\alpha$-intensity", and interpret it as the evolution of the ratio of thermal and non-thermal components, which is also the evolution of the photosphere.

Aims: We provide new insights into the gamma-ray emission from HESS J1912+101, a TeV supernova remnant candidate probably associated with the radio pulsar PSR J1913+1011. Methods: We obtained new observations at 1.5 GHz using the VLA in the D configuration, with the purpose of detecting the radio shell of the putative remnant. In addition, we observed a single pointing at 6.0 GHz toward PSR J1913+1011 to look for a radio pulsar wind nebula. We also studied the properties of the surrounding interstellar medium using data of the 13CO, HI, and infrared emissions, obtained from public surveys. Results: We did not find evidence of a radio shell down to the sensitivity of the new image at 1.5 GHz. We detect faint diffuse emission around PSR J1913+1011 at 6.0 GHz, which could represent a radio pulsar wind nebula powered by the pulsar. We find dense ambient gas at 60 km/s, which shows a good spatial correspondence with the TeV emission only in the western and eastern directions. There is also dense gas near the center of HESS J1912+101, where the TeV emission is weak. Using infrared data, we identify an active star-forming region in the western part of the shell. Conclusions: Based on the poor spatial match between the ambient gas and the TeV emission (which shows a good correlation in the western and eastern directions and an anticorrelation in the other directions), we conclude that the hadronic mechanism alone does not give a satisfactory explanation of the gamma rays from HESS J1912+101. Additional contributions may come from leptonic processes in the shell of the supernova remnant, together with contributions from PSR J1913$+$1011 and its pulsar wind nebula and/or from the star-forming region. A confident determination of the distance to the putative remnant is necessary to determine whether these sources are associated or just appear superimposed in the line of sight.

Braginskii MHD provides a more accurate description of many plasma environments than classical MHD since it actively treats the stress tensor using a closure derived from physical principles. Stress tensor effects nonetheless remain relatively unexplored for solar MHD phenomena, especially in nonlinear regimes. This paper analytically examines nonlinear damping and longitudinal flows of propagating shear Alfv\'en waves. Most previous studies of MHD waves in Braginskii MHD considered the strict linear limit of vanishing wave perturbations. We show that those former linear results only apply to Alfv\'en wave amplitudes in the corona that are so small as to be of little interest, typically a wave energy less than $10^{-11}$ times the energy of the background magnetic field. For observed wave amplitudes, the Braginskii viscous dissipation of coronal Alfv\'en waves is nonlinear and a factor around $10^9$ stronger than predicted by the linear theory. Furthermore, the dominant damping occurs through the parallel viscosity coefficient $\eta_0$, rather than the perpendicular viscosity coefficient $\eta_2$ in the linearized solution. This paper develops the nonlinear theory, showing that the wave energy density decays with an envelope $(1+z/L_d)^{-1}$. The damping length $L_d$ exhibits an optimal damping solution, beyond which greater viscosity leads to lower dissipation as the viscous forces self-organise the longitudinal flow to suppress damping. Although the nonlinear damping greatly exceeds the linear damping, it remains negligible for many coronal applications.

[Abridged] We use high-resolution $HST$ imaging and $e$-MERLIN 1.5-GHz observations of galaxy cores from the LeMMINGs survey to investigate the relation between optical structural properties and nuclear radio emission for a large sample of galaxies. We perform accurate, multi-component decompositions of new surface brightness profiles extracted from $HST$ images for 163 LeMMINGs galaxies and fit up to six galaxy components (e.g., bulges, discs, AGN, bars, rings, spiral arms, and nuclear star clusters) simultaneously with S\'ersic and/or core-S\'ersic models. By adding such decomposition data for 10 LeMMINGs galaxies from our past work, the final sample of 173 nearby galaxies (102 Ss, 42 S0s, 23 Es plus 6 Irr) with bulge stellar mass (typically) M_*, bulge ~ 10^6-10^12.5 M_sun, encompasses all optical spectral classes (LINER, Seyfert, ALG and H II). We show that the bulge mass can be significantly overestimated in many galaxies when components such as bars, rings and spirals are not included in the fits. We additionally implement a Monte Carlo method to determine errors on bulge, disc and other fitted structural parameters. Moving (in the opposite direction) across the Hubble sequence, i.e., from the irregular to elliptical galaxies, we confirm that bulges become larger, more prominent and round. Such bulge dominance is associated with a brighter radio core luminosity. We also find that the radio detection fraction increases with bulge mass. At M_*,bulge > 10^11 M_sun, the radio detection fraction is 77%, declining to 24% for M_bulge < 10^10 M_sun. Furthermore, we observe core-S\'ersic bulges tend to be systematically round and to possess high radio core luminosities and boxy-distorted or pure elliptical isophotes.

The formation of star clusters involves the growth of smaller, gas-rich subclusters through accretion of gas from the giant molecular cloud within which the subclusters are embedded. The two main accretion mechanisms responsible for this are accretion of gas from dense filaments, and from the ambient background of the cloud. We perform simulations of both of these accretion processes onto gas-rich star clusters using coupled smoothed particle hydrodynamics to model the gas, and N-body dynamics to model the stars. We find that, for both accretion processes, the accreting star cluster loses some of its original mass while gaining mass from either the ambient background or the dense filament. The amount of mass lost from both these processes is small compared to the total mass of the cluster. However, in the case of accretion from a background medium, the net effect can be a decrease in the total mass of the cluster if it is travelling fast enough through the ambient medium ($> 4$kms$^{-1}$). We find that the amount of mass lost from the cluster through filamentary accretion is independent of the density, width, or number of filaments funneling gas into the cluster and is always such that the mass of the cluster is constantly increasing with time. We compare our results to idealized prescriptions used to model star cluster formation in larger scale GMC simulations and find that such prescriptions act as an upper limit when describing the mass of the star cluster they represent.

Herbig Ae/Be stars are intermediate-mass pre-main sequence stars undergoing accretion through their circumstellar disk. The optical and infrared (IR) spectra of HAeBe stars show HI emission lines belonging to Balmer, Paschen and Brackett series. We use the archival X-Shooter spectra available for 109 HAeBe stars from Vioque et al. (2018) and analyse the various HI lines present in them. We segregated the stars into different classes based on the presence of higher-order lines in different HI series. We discuss the dependence of the appearance of higher-order lines on the stellar parameters. We find that most massive and younger stars show all the higher-order lines in emission. The stars showing only lower-order lines have Teff < 12000 K and an age range of 5-10 Myr. We perform a Case B line ratio analysis for a sub-sample of stars showing most of the HI lines in emission. We note that all but four stars belonging to the sub-sample show lower HI line ratios than theoretical values, owing to the emitting medium being optically thick. The HI line flux ratios do not depend on the star's spectral type. Further, from the line ratios of lower-order lines and Paschen higher-order lines, we note that line ratios of most HAeBe stars match with electron density value in the range 10^9 - 10^11 cm^-3. The electron temperature, however, could not be ascertained with confidence using the line ratios studied in this work.

Three-nucleon forces are crucial for the accurate description of nuclear systems, including dense matter probed in neutron stars. We explore nuclear Hamiltonians that reproduce two-nucleon scattering data and properties of light nuclei, but differ in the three-nucleon interactions among neutrons. While no significantly improved constraints can be obtained from current astrophysical data, we show that observations of neutron star mergers by next-generation detectors like the proposed Einstein Telescope could provide strong evidence to distinguish between these Hamiltonians.

The WN3/O3 Wolf-Rayet (WR) stars were discovered as part of our survey for WRs in the Magellanic Clouds. The WN3/O3s show the emission lines of a high-excitation WN star and the absorption lines of a hot O-type star, but our prior work has shown that the absorption spectrum is intrinsic to the WR star. Their place in the evolution of massive stars remains unclear. Here we investigate the possibility that they are the products of binary evolution. Although these are not WN3+O3~V binaries, they could still harbor unseen companions. To address this possibility, we have conducted a multi-year radial velocity study of six of the nine known WN3/O3s. Our study finds no evidence of statistically significant radial velocity variations, and allows us to set stringent upper limits on the mass of any hypothetical companion star: for probable orbital inclinations, any companion with a period less than 100 days must have a mass less than 2Mo. For periods less than 10 days, any companion would have to have a mass less than 1Mo. We argue that scenarios where any such companion is a compact object are unlikely. The absorption lines indicate a normal projected rotational velocity, making it unlikely that these stars evolved with the aid of a companion star that has since merged. The modest rotation also suggests that these stars are not the result of homogenous evolution. Thus it is likely that these stars are a normal but short-lived stage in the evolution of massive stars.

The inner 300-500 pc of the Milky Way has some of the most extreme gas conditions in our Galaxy. Physical properties of the Central Molecular Zone (CMZ), including temperature, density, thermal pressure, and turbulent pressure, are key factors for characterizing gas energetics, kinematics, and evolution. The molecular gas in this region is more than an order of magnitude hotter than gas in the Galactic disk, but the mechanism responsible for heating the gas remains uncertain. We characterize the temperature for 16 regions, extending out to a projected radius of $\sim$450 pc. We observe \am\, J,K=(1,1)-(6,6) inversion transitions from SWAG (Survey of Water and Ammonia in the Galactic Center) using the Australia Telescope Compact Array (ATCA), and ammonia lines (J,K) = (8,8)-(14,14) using the 100\,m Green Bank Telescope. Using these two samples we create full Boltzmann plots for every source and fit two rotational temperature components to the data. For the cool component we detect rotational temperatures ranging from 20-80\,K, and for the hot component we detect temperature ranging from 210-580\,K. With this sample of 16 regions, we identify some of the most extreme molecular gas temperatures detected in the Galactic center thus far. We do not find a correlation between gas temperature and Galactocentric radius, and we confirm that these high temperatures are not exclusively associated with actively star-forming clouds. We also investigate temperature and line widths and find (1) no correlation between temperature and line width and (2) the lines are non-thermally broadened indicating that non-thermal motions are dominant over thermal.

We report first results on a method aimed at simultaneously characterising atmospheric parameters and magnetic properties of M dwarfs from high-resolution nIR spectra recorded with SPIRou in the framework of the SPIRou Legacy Survey. Our analysis relies on fitting synthetic spectra computed from MARCS model atmospheres to selected spectral lines, both sensitive and insensitive to magnetic fields. We introduce a new code, $\texttt{ZeeTurbo}$, obtained by including the Zeeman effect and polarised radiative transfer capabilities to $\texttt{Turbospectrum}$. We compute a grid of synthetic spectra with $\texttt{ZeeTurbo}$ for different magnetic field strengths and develop a process to simultaneously constrain $T_{\rm eff}$, $\log{g}$, [M/H], [$\alpha$/Fe] and the average surface magnetic flux. In this paper, we present our approach and assess its performance using simulations, before applying it to six targets observed in the context of the SPIRou Legacy Survey (SLS), namely AU Mic, EV Lac, AD Leo, CN Leo, PM J18482+0741, and DS Leo. Our method allows us to retrieve atmospheric parameters in good agreement with the literature, and simultaneously yields surface magnetic fluxes in the range 2-4 kG with a typical precision of 0.05 kG, in agreement with literature estimates, and consistent with the saturated dynamo regime in which most of these stars are.

We report on the discovery of the companion star to the millisecond pulsar PSR J1835-3259B in the Galactic globular cluster NGC 6652. Taking advantage of deep photometric archival observations acquired through the Hubble Space Telescope in near-ultraviolet and optical bands, we identified a bright and blue object at a position compatible with that of the radio pulsar. The companion is located along the helium-core white dwarf cooling sequence and the comparison with binary evolution models provides a mass of $0.17 \pm 0.02~M_\odot$, a surface temperature of $11500\pm1900$ K and a very young cooling age of only $200\pm100$ Myr. The mass and the age of the companion are compatible with a progenitor star of about $0.87~M_{\odot}$, which started transferring mass to the primary during its evolution along the sub-giant branch and stopped during the early red giant branch phase. Combining together the pulsar mass function and the companion mass, we found that this system is observed at an almost edge-on orbit and hosts a neutron star with a mass of $1.44 \pm 0.06~M_\odot$, thus suggesting a highly non-conservative mass accretion phase. The young age of the WD companion is consistent with the scenario of a powerful, relatively young MSP indicated by the earlier detection of gamma-rays from this system.

The use of the baryonic acoustic oscillations (BAO) datasets offers a unique opportunity to connect the early universe and the late one. In this proceeding, we discuss recent results that used a marginalised likelihood to remove the $H_0-r_d $ degeneracy and then tested it on different dark energy (DE) models. It was found that this approach which does not rely on calibration on $r_d$ or $H_0$, allows us to obtain results, comparable to the ones calculated with standard likelihoods. Here we emphasize on the major differences that we observed for the two different BAO datasets that we employed -- a transversal one, containing only angular BAO measurements, and a mixed one, containing both angular and radial BAO measurements. We see that the two datasets have different statistical preferences for DE models and also different preference for the curvature of the universe.

Trans-Neptunian objects smaller than a few kilometers are difficult to observe directly. They can be detected when they randomly occult a background star. Close to the ecliptic plane, each star is occulted once every tens of thousands of hours, and occultations typically last for less than a second. We present an algorithm, and companion pipeline, for detection of diffractive occultation events. Our approach includes: cleaning the data; an efficient and optimal matched filtering of the light-curves with a template bank of diffractive occultations; treating the red-noise in the light-curves; injection of simulated events for efficiency estimation; and applying data quality cuts. We discuss human vetting of the candidate events in a blinded way to reduce bias caused by the human-in-the-loop. We present Markov Chain Monte Carlo tools to estimate the parameters of candidate occultations, and test them on simulated events. This pipeline is used by the Weizmann Fast Astronomical Survey Telescope (W-FAST).

Based on the linear mixing approach, we calculate the latent heat for crystallizing fully-ionized $^{12}$C/$^{16}$O and $^{16}$O/$^{20}$Ne mixtures in white dwarf (WD) cores for two different parametrizations of the corrections to the linear-mixing energies and with account of ion quantum effects. We report noticeable composition-dependent deviations of the excess entropy in both directions from the standard value of 0.77 per ion. Within the same framework, we evaluate the excess entropy and released or absorbed heat accompanying the exsolution process in solidified WD layers. The inclusion of this effect is shown to be important for reliable interpretation of WD cooling data. We also analyze the latent heat of crystallizing eutectic $^{12}$C/$^{22}$Ne mixture, where we find a qualitative dependence of both the phase diagram and the latent heat behaviour on ion quantum effects. This may be important for the model with $^{22}$Ne distillation in cooling C/O/$^{22}$Ne WD proposed as a solution for the ultramassive WD multi-Gyr cooling anomaly. Astrophysical implications of our findings for crystallizing WD are discussed.

The radial acceleration relation (RAR) of late-type galaxies relates their dynamical acceleration, $g_\text{obs}$, to that sourced by baryons alone, $g_\text{bar}$, across their rotation curves. Literature fits to the RAR have fixed the galaxy parameters on which the relation depends -- distance, inclination, luminosity and mass-to-light ratios -- to their maximum a priori values with an uncorrelated Gaussian contribution to the uncertainties on $g_\text{bar}$ and $g_\text{obs}$. In reality these are free parameters of the fit, contributing systematic rather than statistical error. Assuming a range of possible functional forms for the relation with or without intrinsic scatter (motivated by Modified Newtonian Dynamics with or without the external field effect), I use Hamiltonian Monte Carlo to perform the full joint inference of RAR and galaxy parameters for the SPARC dataset. This reveals the intrinsic RAR underlying that observed. I find an acceleration scale $a_0=(1.19 \pm 0.04 \, \text{(stat)} \pm 0.09 \, \text{(sys)}) \: \times \: 10^{-10}$ m s$^{-2}$, an intrinsic scatter $\sigma_\text{int}=(0.034 \pm 0.01 \, \text{(stat)} \pm 0.01 \, \text{(sys)})$ dex (assuming the SPARC error model is reliable) and weak evidence for the external field effect. I make summary statistics of all my analyses publicly available for future SPARC studies or applications of a calibrated RAR, for example redshift-independent distance measurement.

Most protoplanetary discs are thought to undergo violent and frequent accretion outbursts, during which the accretion rate and central luminosity are elevated for several decades. This temporarily increases the disc temperature, leading to the sublimation of ice species as snowlines move outwards. In this paper, we investigate how an FUor-type accretion outburst alters the growth and appearance of dust aggregates at different locations in protoplanetary discs. We develop a model based on the Monte Carlo approach to simulate locally the coagulation and fragmentation of icy dust particles and investigate different designs for their structure and response to sublimation. Our main finding is that the evolution of dust grains located between the quiescent and outburst water snowlines is driven by significant changes in composition and porosity. The time required for the dust population to recover from the outburst and return to a coagulation/fragmentation equilibrium depends on the complex interplay of coagulation physics and outburst properties, and can take up to 4500 yr at 5 au. Pebble-sized particles, the building blocks of planetesimals, are either deprecated in water ice or completely destroyed, respectively resulting in drier planetesimals or halting their formation altogether. When accretion outbursts are frequent events, the dust can be far from collisional equilibrium for a significant fraction of time, offering opportunities to track past outbursts in discs at millimetre wavelengths. Our results highlight the importance of including accretion outbursts in models of dust coagulation and planet formation.

We present the first observational infrared luminosity function (IRLF) measurement in the Epoch of Reionization (EoR) based on a UV-selected galaxy sample with ALMA spectroscopic observations. Our analysis is based on the ALMA large program Reionization Era Bright Emission Line Survey (REBELS), which targets 42 galaxies at $\mathrm{z=6.4-7.7}$ with [CII] 158$\micron$ line scans. 16 sources exhibit a dust detection, 15 of which are also spectroscopically confirmed through the [CII] line. The IR luminosities of the sample range from $\log L_{IR}/L_\odot=11.4$ to 12.2. Using the UVLF as a proxy to derive the effective volume for each of our target sources, we derive IRLF estimates, both for detections and for the full sample including IR luminosity upper limits. The resulting IRLFs are well reproduced by a Schechter function with the characteristic luminosity of $\log L_{*}/L_\odot=11.6^{+0.2}_{-0.1}$. Our observational results are in broad agreement with the average of predicted IRLFs from simulations at $z\sim7$. Conversely, our IRLFs lie significantly below lower redshift estimates, suggesting a rapid evolution from $z\sim4$ to $z\sim7$, into the reionization epoch. The inferred obscured contribution to the cosmic star-formation rate density at $z\sim7$ amounts to $\mathrm{log(SFRD/M_{\odot}/yr/Mpc^{3}) = -2.66^{+0.17}_{-0.14} }$ which is at least $\sim$10\% of UV-based estimates. We conclude that the presence of dust is already abundant in the EoR and discuss the possibility of unveiling larger samples of dusty galaxies with future ALMA and JWST observations.

The continuously growing number of objects orbiting around the Earth is expected to be accompanied by an increasing frequency of objects re-entering the Earth's atmosphere. Many of these re-entries will be uncontrolled, making their prediction challenging and subject to several uncertainties. Traditionally, re-entry predictions are based on the propagation of the object's dynamics using state-of-the-art modelling techniques for the forces acting on the object. However, modelling errors, particularly related to the prediction of atmospheric drag may result in poor prediction accuracies. In this context, we explore the possibility to perform a paradigm shift, from a physics-based approach to a data-driven approach. To this aim, we present the development of a deep learning model for the re-entry prediction of uncontrolled objects in Low Earth Orbit (LEO). The model is based on a modified version of the Sequence-to-Sequence architecture and is trained on the average altitude profile as derived from a set of Two-Line Element (TLE) data of over 400 bodies. The novelty of the work consists in introducing in the deep learning model, alongside the average altitude, three new input features: a drag-like coefficient (B*), the average solar index, and the area-to-mass ratio of the object. The developed model is tested on a set of objects studied in the Inter-Agency Space Debris Coordination Committee (IADC) campaigns. The results show that the best performances are obtained on bodies characterised by the same drag-like coefficient and eccentricity distribution as the training set.

We consider the structure and physical properties of specific classes of neutron, quark, and Bose-Einstein Condensate stars in the conformally invariant Weyl geometric gravity theory. The basic theory is derived from the simplest conformally invariant action, constructed, in Weyl geometry, from the square of the Weyl scalar, the strength of the Weyl vector, and a matter term, respectively. The action is linearized in the Weyl scalar by introducing an auxiliary scalar field. To keep the theory conformally invariant the trace condition is imposed on the matter energy-momentum tensor. The field equations are derived by varying the action with respect to the metric tensor, Weyl vector field and scalar field. By adopting a static spherically symmetric interior geometry, we obtain the field equations, describing the structure and properties of stellar objects in Weyl geometric gravity. The solutions of the field equations are obtained numerically, for different equations of state of the neutron and quark matter. More specifically, constant density stellar models, and models described by the stiff fluid, radiation fluid, quark bag model, and Bose-Einstein Condensate equations of state are explicitly constructed numerically in both general relativity and Weyl geometric gravity, thus allowing an in depth comparison between the predictions of these two gravitational theories. As a general result it turns out that for all the considered equations of state, Weyl geometric gravity stars are more massive than their general relativistic counterparts. As a possible astrophysical application of the obtained results we suggest that the recently observed neutron stars, with masses in the range of 2$M_{\odot}$ and 3$M_{\odot}$, respectively, could be in fact conformally invariant Weyl geometric neutron or quark stars.

Consequences of the consistent exact solution of Einstein-Cartan equation on the time dependence of Hubble parameter is discussed. The torsion leads to a space and time dependent expansion parameter which results into nontrivial windows of Hubble parameter between diverging behaviour. Only one window shows a period of decreasing followed by increasing time dependence. Provided a known cosmological constant and the present values of Hubble and deceleration parameter this changing time can be given in the past as well as the ending time of the windows or universe. From the metric with torsion outside matter it is seen that torsion can feign dark matter.

The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evidence for such signals and set sky averaged 95% upper limits on their amplitude h 95 . The most sensitive frequency is 10nHz with h 95 = 9.1 10-15 . We achieved the best upper limit to date at low and high frequencies of the PTA band thanks to improved effective cadence of observations. In our analysis, we have taken into account the recently discovered common red noise process, which has an impact at low frequencies. We also find that the peculiar noise features present in some pulsars data must be taken into account to reduce the false alarm. We show that using custom noise models is essential in searching for continuous gravitational wave signals and setting the upper limit.

The scaling of magnetic fluctuations provides crucial information for the understanding of solar wind turbulence. However, the observed magnetic fluctuations contain not only turbulence but also magnetic structures, leading to the violation of the time stationarity. This violation would conceal the true scaling and influence the determination of the sampling angle with respect to the local background magnetic field. Here, to investigate the scaling anisotropy, we utilize an easy but effective criterion $\phi<10^\circ$ to ensure the time stationarity of the magnetic field, where $\phi$ is the angle between the two averaged magnetic fields after cutting the interval into two halves. We study the scaling anisotropy using higher-order statistics of structure functions under the condition of stationarity for the near-Sun solar wind turbulence for the first time based on measurements obtained from Parker Solar Probe (PSP) at 0.17 au. We find that the scaling indices $\xi$ of magnetic field show a linear dependence on the order $p$ close to $\xi(p)=p/4$. The multifractal scaling of magnetic-trace structure functions becomes monoscaling close to $\xi(p)=p/3$ with the local magnetic field perpendicular to the sampling direction and close to $\xi(p)=p/4$ with the local magnetic field parallel to the sampling direction when measured with the stationary background magnetic field. The scaling of velocity-trace structure functions has similar but less significant changes. The near-Sun solar wind turbulence displays different scaling anisotropies with the near-Earth solar wind turbulence, suggesting the evolution of the nonlinear interaction process during the solar wind expansion.

We investigate the effect of $\Delta$ baryons on the radial oscillations of neutron and hyperon stars, employing the density-dependent relativistic mean-field model. The spin-$3/2$ baryons are described by the Rarita-Schwinger Lagrangian density. The baryon-meson coupling constants for the spin-3/2 decuplet and the spin-1/2 baryonic octet are calculated using a unified approach relying on the fact that the Yukawa couplings present in the Lagrangian density of the mean-field models must be invariant under the SU(3) and SU(6) group transformations. We calculate the 20 lowest eigenfrequencies and corresponding oscillation functions of $\Delta$-admixtured nuclear (N$\Delta$) and hyperonic matter (NH$\Delta$) by solving the Sturm-Liouville boundary value problem and also verifying its validity. We see that the lowest mode frequencies for N+$\Delta$ and N+H EoSs are higher as compared to the pure nucleonic matter because of the delta and hyperonic admixtures. Furthermore, the separation between consecutive modes increases with the addition of hyperons and $\Delta$s.

We present a comprehensive study on the self-interaction cross-section of puffy dark matter (DM) particles, which have a significant intrinsic size compared to their Compton wavelength. For such puffy DM self-interaction cross-section in the resonant and classical regimes, our study demonstrates the significance of the Yukawa potential and the necessity of partial wave analysis: (i) Due to the finite-size effect of puffy DM particles, the new Yukawa potential of puffy DM is found to enlarge the Born-effective regime for the self-interaction cross-section, compared with the point-like DM; (ii) Our partial wave analysis shows that depending on the value of the ratio between $R_{\chi}$ (radius of a puffy DM particle) and $1/m_{\phi}$ (force range), the three regimes (Born-effective, resonant and classical) for puffy DM self-interaction cross-section can be very different from the point-like DM; (iii) We find that to solve the small-scale anomalies via self-interacting puffy DM, the Born-effective and the resonant regimes exist for dwarf galaxies, while for the cluster and Milky Way galaxy the non-Born regime is necessary.

The short answer is $\textit{probably no}$. Specifically, this paper considers a recent body of work which suggests that general relativity requires neither the support of dark matter halos, nor unconventional baryonic profiles, nor any infrared modification, to be consistent after all with the anomalously rapid orbits observed in many galactic discs. In particular, the gravitoelectric flux is alleged to collapse nonlinearly into regions of enhanced force, in an analogue of the colour-confining chromoelectric flux tube model which has yet to be captured by conventional post-Newtonian methods. However, we show that the scalar gravity model underpinning this proposal is wholly inconsistent with the nonlinear Einstein equations, which themselves appear to prohibit the linear confinement-type potentials which could indicate a disordered gravitational phase. Our findings challenge the fidelity of the previous Euclidean lattice analyses. We confirm by direct calculation using a number of perturbation schemes and gauges that the next-to-leading order gravitoelectric correction to the rotation curve of a reasonable baryonic profile would be imperceptible. The `gravitoelectric flux collapse' programme was also supported by using intragalactic lensing near a specific galactic baryon profile as a field strength heuristic. We recalculate this lensing effect, and conclude that it has been overstated by three orders of magnitude. As a by-product, our analysis suggests fresh approaches to (i) the fluid ball conjecture and (ii) gravitational energy localisation, both to be pursued in future work. In summary, whilst it may be interesting to consider the possibility of confinement-type effects in gravity, we may at least conclude here that confinement-type effects $\textit{cannot play any significant part}$ in explaining flat or rising galactic rotation curves without dark matter halos.

In this paper, we develop a framework to re-examine the weak lensing (including the microlensing) effects of the static spherical symmetric wormhole in terms of the radial equation of state $\eta=\frac{p_r}{\rho}$ (REoS). As for its application, we calculate its magnification, and event rate under this REoS, in which we show that the maximal value of magnification of the Ellis-Bronnikov wormhole is only related to the relative position and intrinsic angle, whose the maximal value is around five. For the event rate, our results indicate that one cannot distinguish the Eillis-Bronnikov wormhole and charged wormhole, but its order is much higher than the vacuum case, in which all these metrics belong to the static spherical symmetric wormhole metric. By calculating the lensing equation of the static spherical symmetric wormhole, we find an explicit formula between the maximal number of images of the wormhole and $\eta$. It shows that this relation is consistent with the classical wormhole, but the case for wormhole with quantum corrections is still mysterious. Our new method may shed new light on distinguishing the wormhole and blackhole via the event rate.