Papers for Wednesday, Jun 08 2022

We measure homogeneous distances to M31 and 38 associated stellar systems ($-$16.8$\le M_V \le$ $-$6.0), using time-series observations of RR Lyrae stars taken as part of the Hubble Space Telescope Treasury Survey of M31 Satellites. From $>700$ orbits of new/archival ACS imaging, we identify $>4700$ RR Lyrae stars and determine their periods and mean magnitudes to a typical precision of 0.01 days and 0.04 mag. Based on Period-Wesenheit-Metallicity relationships consistent with the Gaia eDR3 distance scale, we uniformly measure heliocentric and M31-centric distances to a typical precision of $\sim20$ kpc (3%) and $\sim10$ kpc (8%), respectively. We revise the 3D structure of the M31 galactic ecosystem and: (i) confirm a highly anisotropic spatial distribution such that $\sim80$% of M31's satellites reside on the near side of M31; this feature is not easily explained by observational effects; (ii) affirm the thin (rms $7-23$ kpc) planar "arc" of satellites that comprises roughly half (15) of the galaxies within 300 kpc from M31; (iii) reassess physical proximity of notable associations such as the NGC 147/185 pair and M33/AND XXII; and (iv) illustrate challenges in tip-of-the-red-giant branch distances for galaxies with $M_V > -9.5$, which can be biased by up to 35%. We emphasize the importance of RR Lyrae for accurate distances to faint galaxies that should be discovered by upcoming facilities (e.g., Rubin Observatory). We provide updated luminosities and sizes for our sample. Our distances will serve as the basis for future investigation of the star formation and orbital histories of the entire known M31 satellite system.

The gravitational waves emitted by massive black hole binaries in the LISA band can be lensed. Wave-optics effects in the lensed signal are crucial when the Schwarzschild radius of the lens is smaller than the wavelength of the radiation. These frequency-dependent effects can enable us to infer the lens parameters, possibly with a single detection alone. In this work, we assess the observability of wave-optics effects with LISA by performing an information-matrix analysis using analytical solutions for both point-mass and singular isothermal sphere lenses. We use gravitational-waveform models that include the merger, ringdown, higher harmonics, and aligned spins to study how waveform models and source parameters affect the measurement errors in the lens parameters. We find that previous work underestimated the observability of wave-optics effects and that LISA can detect lensed signals with higher impact parameters and lower lens masses.

About a hundred tidal disruption events (TDEs) have been observed and they exhibit a wide range of emission properties both at peak and over their lifetimes. Some TDEs peak predominantly at X-ray energies while others radiate chiefly at UV and optical wavelengths. While the peak luminosities across TDEs show distinct properties, the evolutionary behavior can also vary between TDEs with similar peak emission properties. At late time, some optical TDEs rebrighten in X-rays, while others maintain strong UV/optical emission components. In this Letter, we conduct three-dimensional general relativistic radiation magnetohydrodynamics simulations of TDE accretion disks at varying accretion rates ranging from a few to a few tens of the Eddington accretion rate. We make use of Monte Carlo radiative transfer simulations to calculate the reprocessed spectra at various inclinations and at different evolutionary stages. We confirm the unified model proposed by Dai et al. (2018), which predicts that the observed emission largely depends on the viewing angle of the observer with respect to the disk orientation (X-ray strong when viewed face-on and UV/optically strong when viewed disk-on). What is more, we find that disks with higher accretion rates have elevated wind and disk densities, which increases the reprocessing of the high-energy radiation and thus generally augments the optical-to-X-ray flux ratio along a particular viewing angle. This implies that at later times, as the accretion level declines, we expect the funnel to effectively open up and allow more X-rays to leak out along intermediate viewing angles. Such dynamical model for TDEs can provide a natural explanation for the diversity in the emission properties observed in TDEs at peak and along their temporal evolution.

The aim of this work is to improve models for the gamma-ray discrete or small-scale structure related to H2 interstellar gas. Reliably identifying this contribution is important to disentangle gamma-ray point sources from interstellar gas, and to better characterize extended gamma-ray signals. Notably, the Fermi-LAT Galactic center (GC) excess, whose origin remains unclear, might be smooth or point-like. If the data contain a point-like contribution that is not adequately modeled, a smooth GC excess might be erroneously deemed to be point-like. We improve models for the H2-related gamma-ray discrete emission for a $50^\circ \times 1^\circ$ region along the Galactic plane via H2 proxies better suited to trace these features. We find that these are likely to contribute significantly to the gamma-ray Fermi-LAT data in this region, and the brightest ones are likely associated with detected Fermi-LAT sources, a compelling validation of this methodology. We discuss prospects to extend this methodology to other regions of the sky and implications for the characterization of the GC excess.

Type Ia supernovae (SNe Ia) play a crucial role as standardizable candles in measurements of the Hubble constant and dark energy. Increasing evidence points towards multiple possible explosion channels as the origin of normal SNe Ia, with possible systematic effects on the determination of cosmological parameters. We present, for the first time, a comprehensive comparison of publicly-available SN Ia model nucleosynthetic data with observations of late-time light curve observations of SN Ia events. These models span a wide range of white dwarf (WD) progenitor masses, metallicities, explosion channels, and numerical methodologies. We focus on the influence of $^{57}$Ni and its isobaric decay product $^{57}$Co in powering the late-time ($t > 1000$ d) light curves of SNe Ia. $^{57}$Ni and $^{57}$Co are neutron-rich relative to the more abundant radioisotope $^{56}$Ni, and are consequently a sensitive probe of neutronization at the higher densities of near-Chandrashekhar (near-$M_{\rm Ch}$) progenitor WDs. We demonstrate that observations of one SN Ia event, SN 2015F is only consistent with a sub-$M_{\rm Ch}$ WD progenitor. Observations of four other events (SN 2011fe, SN 2012cg, SN 2014J, SN2013aa) are consistent with both near-$M_{\rm Ch}$ and sub-$M_{\rm Ch}$ progenitors. Continued observations of late-time light curves of nearby SNe Ia will provide crucial information on the nature of the SN Ia progenitors.

The delay time distribution of neutron star mergers provides critical insights into binary evolution processes and the merger rate evolution of compact object binaries. However, current observational constraints on this delay time distribution rely on the small sample of Galactic double neutron stars (with uncertain selection effects), a single multimessenger gravitational-wave event, and indirect evidence of neutron star mergers based on $r$-process enrichment. We use a sample of 67 host galaxies of short gamma-ray bursts to place novel constraints on the delay time distribution, and leverage this result to infer the merger rate evolution of compact object binaries containing neutron stars. We recover a power-law slope of $\alpha = -1.82^{+0.40}_{-0.41}$ and a minimum delay time of $t_\mathrm{min} = 171^{+74}_{-91}~\mathrm{Myr}$ (median and 90% credible interval), with the maximum delay time constrained to $t_\mathrm{max} > 7.73~\mathrm{Gyr}$ at 99% credibility. We find these constraints to be broadly consistent with theoretical expectations, although our recovered power-law slope is substantially steeper than the conventional value of $t_\mathrm{d}^{-1}$, and our minimum delay time is larger than the typically assumed value of $10~\mathrm{Myr}$. Pairing this cosmological probe of the fate of compact object binary systems with the Galactic population of double neutron stars will be crucial for understanding the unique selection effects governing both of these populations. In addition to probing a significantly larger redshift regime of neutron star mergers than possible with current gravitational-wave detectors, complementing our results with future multimessenger gravitational-wave events will also help determine if short gamma-ray bursts ubiquitously result from compact object binary mergers.

We present Hubble Space Telescope (HST) imaging of the site of SN 2015bh in the nearby spiral galaxy NGC 2770 taken between 2017 and 2019, nearly four years after the peak of the explosion. In 2017-2018, the transient fades steadily in optical filters before declining more slowly to $F814W = -7.2$ mag in 2019, nearly 4 mag below the level of its eruptive luminous blue variable (LBV) progenitor observed with HST in 2008-2009. The source fades at a constant color of $F555W - F814W = 0.3$ mag until 2018, similar to SN 2009ip and consistent with a spectrum dominated by continued interaction of the ejecta with circumstellar material (CSM). A deep optical spectrum obtained in 2021 lacks signatures of ongoing interaction ($L_{\mathrm{H}\alpha} \lesssim 10^{39}$ erg s$^{-1}$ for broadened emission $\lesssim 2000$ km s$^{-1}$), but but indicates the presence of a nearby H II region ($\lesssim 300$ pc). The color evolution of the fading source makes it unlikely that emission from a scattered light echo or binary OB companion of the progenitor contributes significantly to the flattening of the late-time light curve. The remaining emission in 2019 may plausibly be attributed to an unresolved ($\lesssim 3$ pc), young stellar cluster. Importantly, the color evolution of SN 2015bh also rules out scenarios in which the surviving progenitor evolves back to a hot, quiescent, optically faint state or is obscured by nascent dust. The simplest explanation is that the massive progenitor did not survive. SN 2015bh, therefore, likely represents a remarkable example of the terminal explosion of a massive star preceded by decades of end-stage eruptive variability.

The formation of gas giant planets must occur during the first few Myr of a star's lifetime, when the protoplanetary disc still contains sufficient gas to be accreted onto the planetary core. The majority of protoplanetary discs are exposed to strong ultraviolet irradiation from nearby massive stars, which drives winds and depletes the mass budget for planet formation. It remains unclear to what degree external photoevaporation affects the formation of massive planets. In this work, we present a simple one dimensional model for the growth and migration of a massive planet under the influence of external FUV fields. We find that even moderate FUV fluxes $F_\mathrm{FUV}\gtrsim 100 \, G_0$ have a strong influence on planet mass and migration. By decreasing the local surface density and shutting off accretion onto the planet, external irradiation suppresses planet masses and halts migration early. The distribution of typical stellar birth environments can therefore produce an anti-correlation between semi-major axis and planet mass, which may explain the apparent decrease in planet occurrence rates at orbital periods $P_\mathrm{orb}\gtrsim 10^3$ days. Even moderate fluxes $F_\mathrm{FUV}$ strongly suppress giant planet formation and inward migration for any initial semi-major axis if the stellar host mass $M_*\lesssim 0.5\, M_\odot$, consistent with findings that massive planet occurrence is much lower around such stars. The outcomes of our prescription for external disc depletion show significant differences to the current approximation adopted in state-of-the-art population synthesis models, motivating future careful treatment of this important process.

A significant point-like component from the small scale (or discrete) structure in the H2 interstellar gas might be present in the Fermi-LAT data, but modeling this emission relies on observations of rare gas tracers only available in limited regions of the sky. Identifying this contribution is important to discriminate gamma-ray point sources from interstellar gas, and to better characterize extended gamma-ray sources. We design and train convolutional neural networks to predict this emission where observations of these rare tracers do not exist and discuss the impact of this component on the analysis of the Fermi-LAT data. In particular, we evaluate prospects to exploit this methodology in the characterization of the Fermi-LAT Galactic center excess through accurate modeling of point-like structures in the data to help distinguish between a point-like or smooth nature for the excess. We show that deep learning may be effectively employed to model the gamma-ray emission traced by these rare H2 proxies within statistical significance in data-rich regions, supporting prospects to employ these methods in yet unobserved regions.

Elucidating dark matter density profiles in the Galactic dwarf satellites is essential to understanding not only the quintessence of dark matter but also the evolution of the satellites themselves. In this work, we present the current constraints on dark matter densities in the ultra-faint dwarf (UFD) and diffuse galaxies in the Milky Way. Applying our constructed non-spherical mass models to the currently available kinematic data of the 25 UFDs and 2 diffuse satellites, we find that whereas most of the galaxies have huge uncertainties on the inferred dark matter density profiles, Eridanus~II, Segue~I, and Willman~1 favor cusped central profiles even considering effects of a prior bias. We compare our results with the simulated subhalos on the plane between the dark matter density at 150~pc and the pericenter distance. We find that the most observed satellites and the simulated subhalos are similarly distributed on this plane, except for Antlia~2, Crater~2, and Tucana~3 which deviate considerably from the simulated ones. Despite considering tidal effects, the subhalos detected by commonly-used subhalo finders can be difficult to explain such a huge deviation. We also estimate the dynamical mass-to-light ratios of the satellites and confirm the ratio is linked to stellar mass and metallicity. Using these scaling relations and the empirical tidal evolution models, we infer that the scaling laws can be sustained, and furthermore, a metallicity scatter at the faint-end of the mass-metallicity relation can be explained by tidal stripping of cusped halos.

We present optical-line gas metallicity diagnostics established by the combination of local SDSS galaxies and the largest compilation of extremely metal-poor galaxies (EMPGs) including new EMPGs identified by the Subaru EMPRESS survey. A total of 103 EMPGs are included that cover a large parameter space of magnitude (Mi=-19 to -7) and H-beta equivalent width (10-600 Ang), i.e., wide ranges of stellar mass and star-formation rate. Using reliable metallicity measurements of the direct method for these galaxies, we derive the relationships between strong optical-line ratios and gas-phase metallicity over the range of 12+log(O/H)=~6.9-8.9 corresponding to 0.02-2 solar metallicity Zsun. We confirm that R23-index, ([OIII]+[OII])/H-beta, is the most accurate metallicity indicator with the metallicity uncertainty of 0.14 dex over the range among various popular metallicity indicators. The other metallicity indicators show large scatters in the metal-poor range (<0.1 Zsun). It is explained by our CLOUDY photoionization modeling that, unlike R23-index, the other metallicity indicators do not use a sum of singly and doubly ionized lines and cannot trace both low and high ionization gas. We find that the accuracy of the metallicity indicators is significantly improved, if one uses H-beta equivalent width measurements that tightly correlate with ionization states. In this work, we also present the relation of physical properties with UV-continuum slope beta and ionization production rate xi_ion derived with GALEX data for the EMPGs, and provide local anchors of galaxy properties together with the optical-line metallicity indicators that are available in the form of ASCII table and useful for forthcoming JWST spectroscopic studies.

The NASA TESS mission has discovered many transiting planets orbiting bright nearby stars, and high-resolution imaging studies have revealed that a number of these exoplanet hosts reside in binary or multiple star systems. In such systems, transit observations alone cannot determine which star in the binary system actually hosts the orbiting planet. The knowledge of which star the planet orbits is necessary for determining accurate physical properties for the planet, especially its true radius and mean bulk density. We derived the mean stellar densities for the components of 23 binary systems using the light curve transit shape and the binary flux ratio from speckle imaging, then tested the consistency with stellar models to determine which component is the host star. We found that 70% of the TESS transiting planets in our sample orbit the primary star.

In the context of the S-PLUS Fornax Project (S+FP), we analyze the galaxy population in the direction of the Fornax cluster ($D\approx 20$~Mpc). We have 23 fields of size $1.4^{\circ}\times 1.4^{\circ}$, covering the projected positions of 999 Fornax galaxies reported in the literature. 244 of those galaxies are detected with confident photometry in our fields which were observed simultaneously in 12 photometric bands. Besides those of Fornax galaxies, we obtained confident structural and photometric parameters for $\approx 3\times10^5$ additional galaxies detected in our fields. In this work we present preliminary results on the characterization of the galaxy population of the Fornax cluster with respect to the background galaxy population. Among other goals, we expect that such a characterization provides photometric criteria to identify new candidate members of the cluster.

TIGRE (Telescopio Internacional de Guanajuato Rob\'otico Espectrosc\'opico) has been operating in fully robotic mode in the Observatory La Luz (Guanajuato, Mexico) since the end of 2013. With its sole instrument, HEROS, an \'echelle spectrograph with a spectral resolution R~20000, TIGRE has collected more than 48000 spectra of 1151 different sources with a total exposure time of more than 11000 hours in these eight years. Here we briefly describe the system and the upgrades performed during the last years. We present the statistics of the weather conditions at the La Luz Observatory, emphasizing the characteristics that affect the astronomical observations. We evaluate the performance and efficiency of TIGRE, both optical and operational, and describe the improvements of the system implemented to optimize the telescope's performance and meet the requirements of the astronomer in terms of timing constraints for the observations and the quality of the spectra. We describe the actions taken to slow down the optical efficiency loss due to the aging of the optical surfaces as well as the upgrades of the scheduler and the observing procedures to minimize the time lost due to interrupted observations or observations that do not reach the required quality. Finally, we highlight a few of the main scientific results obtained with TIGRE data.

Observations of low order CO transitions represent the most direct way to study galaxies' cold molecular gas, the fuel of star formation. Here we present the first detection of CO(2-1) in a galaxy lying on the main-sequence of star-forming galaxies at z>6. Our target, G09-83808 at z=6.03, has a short depletion time-scale of T_dep~50Myr and a relatively low gas fraction of M_gas/M_star=0.30 that contrasts with those measured for lower redshift main-sequence galaxies. We conclude that this galaxy is undergoing a starburst episode with a high star formation efficiency that might be the result of gas compression within its compact rotating disk. Its starburst-like nature is further supported by its high star formation rate surface density, thus favoring the use of the Kennicutt-Schmidt relation as a more precise diagnostic diagram for starbursts. Without further significant gas accretion, this galaxy would become a compact, massive quiescent galaxy at z~5.5. In addition, we find that the calibration for estimating ISM masses from dust continuum emission satisfactorily reproduces the gas mass derived from the CO(2-1) transition (within a factor of ~2). This is in line with previous studies claiming a small redshift evolution in the gas-to-dust ratio of massive, metal-rich galaxies. In the absence of gravitational amplification, this detection would have required of order ~1000h of observing time. The detection of cold molecular gas in unlensed star-forming galaxies at high redshifts is thus prohibitive with current facilities and requires a ten-fold improvement in sensitivity, such as that envisaged for the ngVLA.

Galaxy mergers have been observed to trigger nuclear activity by feeding gas to the central supermassive black hole. One such class of objects are Ultra Luminous InfraRed Galaxies (ULIRGs), which are mostly late stage major mergers of gas-rich galaxies. Recently, large-scale ($\sim$100 kpc) radio continuum emission has been detected in a select number of ULIRGs, all of which also harbour powerful Active Galactic Nuclei (AGN). This hints at the presence of large-scale radio emission being evidence for nuclear activity. Exploring the origin of this radio emission and its link to nuclear activity requires high sensitivity multi-frequency data. We present such an analysis of the ULIRG Mrk 273. Using the International LOFAR telescope (ILT), we detected spectacular large-scale arcs in this system. This detection includes, for the first time, a giant $\sim$190 kpc arc in the north. We propose these arcs are fuelled by a low power radio AGN triggered by the merger. We also identified a bright $\sim$45 kpc radio ridge, which is likely related to the ionised gas nebula in that region. We combined this with high sensitivity data from APERture Tile In Focus (Apertif) and archival data from the Very Large Array (VLA) to explore the spectral properties. The ILT simultaneously allowed us to probe the nucleus at a resolution of $\sim$0.3 arcsec, where we detected three components, and, for the first time, diffuse emission around these components. Combining this with archival high frequency VLA images of the nucleus allowed us to detect absorption in one component, and a steep spectrum radio AGN in another. We then extrapolate from this case study to the importance of investigating the presence of radio emission in more ULIRGs and what it can tell us about the link between mergers and the presence of radio activity.

Multi-messenger observations of transient astrophysical sources have the potential to characterize the highest energy accelerators and the most extreme environments in the Universe. Detection of neutrinos, in particular tau neutrinos generated by neutrino oscillations in transit from their sources to Earth, is possible for neutrino energies above 10 PeV using optical Cherenkov detectors imaging upward-moving extensive air showers (EAS). These EAS are produced from Earth interacting tau neutrinos leading to tau leptons that subsequently decay in the atmosphere. We compare neutrino detection sensitivities for generic short- and long-burst transient neutrino sources and sensitivities to the diffuse neutrino flux for the second generation Extreme Universe Space Observatory on a Super Pressure Balloon (EUSO-SPB2) balloon-borne mission and the proposed space-based Probe of Extreme Multi-Messenger Astrophysics (POEMMA) mission.

We consider scalar-vector-tensor (SVT) theories with second-order equations of motion and tensor propagation speed equivalent to the speed of light. Under the sub-horizon and the quasi-static approximations we find analytical formulae for an effective dark energy fluid, i.e., sound speed, anisotropic stress as well as energy density and pressure. We took advantage of our general, analytical fluid description and showed that it is possible to design SVT cosmological models which are degenerate with $\Lambda$CDM at the background level while having gravity strength $G_{\rm eff}<G_{\rm N}$ at late-times as well as non-vanishing dark energy perturbations. We implemented SVT designer models in the widely used Boltzmann solver CLASS thus making it possible to test SVT models against astrophysical observations. Our effective fluid approach to SVT models reveals non trivial behaviour in the sound speed and the anisotropic stress well worth an investigation in light of current discrepancies in cosmological parameters such as $H_0$ and $\sigma_8$.

In the last 15 years different ground-based spectroscopic surveys have been started (and completed) with the general aim of delivering stellar parameters and elemental abundances for large samples of Galactic stars, complementing Gaia astrometry. Among those surveys, the Gaia-ESO Public Spectroscopic Survey (GES), the only one performed on a 8m class telescope, was designed to target 100,000 stars using FLAMES on the ESO VLT (both Giraffe and UVES spectrographs), covering all the Milky Way populations, with a special focus on open star clusters. This article provides an overview of the survey implementation (observations, data quality, analysis and its success, data products, and releases), of the open cluster survey, of the science results and potential, and of the survey legacy. A companion article (Gilmore et al.) reviews the overall survey motivation, strategy, Giraffe pipeline data reduction, organisation, and workflow. The GES has determined homogeneous good-quality radial velocities and stellar parameters for a large fraction of its more than 110,000 unique target stars. Elemental abundances were derived for up to 31 elements for targets observed with UVES. Lithium abundances are delivered for about 1/3 of the sample. The analysis and homogenisation strategies have proven to be successful; several science topics have been addressed by the Gaia-ESO consortium and the community, with many highlight results achieved. The final catalogue has been released through the ESO archive at the end of May 2022, including the complete set of advanced data products. In addition to these results, the Gaia-ESO Survey will leave a very important legacy, for several aspects and for many years to come.

We report the discovery of three new "near-pristine" Lyman Limit Systems (LLSs), with metallicities ~1/1000 solar, at redshifts 2.6, 3.8 and 4.0, with a targeted survey at the Keck Observatory. High resolution echelle spectra of eight candidates yielded precise column densities of hydrogen and weak, but clearly detected, metal lines in seven LLSs; we previously reported the one remaining, apparently metal-free LLS, to have metallicity <1/10000 solar. Robust photoionisation modelling provides metallicities [Si/H] = -3.05 to -2.94, with 0.26 dex uncertainties (95% confidence) for three LLSs, and [Si/H] >~ -2.5 for the remaining four. Previous simulations suggest that near-pristine LLSs could be the remnants of PopIII supernovae, so comparing their detailed metal abundances with nucleosynthetic yields from supernovae models is an important goal. Unfortunately, at most two different metals were detected in each new system, despite their neutral hydrogen column densities (10^{19.2-19.4} cm^{-2}) being two orders of magnitude larger than the two previous, serendipitously discovered near-pristine LLSs. Nevertheless, the success of this first targeted survey for near-pristine systems demonstrates the prospect that a much larger, future survey could identify clear observational signatures of PopIII stars. With a well-understood selection function, such a survey would also yield the number density of near-pristine absorbers which, via comparison to future simulations, could reveal the origin(s) of these rare systems.

Spatially extended halos of Lyman Alpha (Ly$\alpha$) emission are now ubiquitously found around high-redshift star-forming galaxies. But our understanding of the nature and powering mechanisms of these halos is still hampered by the complex radiative transfer effects of the Ly$\alpha$ line and limited angular resolution. In this paper, we present resolved Multi Unit Spectroscopic Explorer (MUSE) observations of SGAS J122651.3+215220, a strongly-lensed pair of L* galaxies at z=2.92 embedded in a Ly$\alpha$ halo of $L_{\mathrm{Ly}\alpha}=(6.2\pm1.3)\times 10^{42}\,\mathrm{erg}\,\mathrm{s}^{-1}$. Globally, the system shows a line profile that is markedly asymmetric and redshifted, but its width and peak shift vary significantly across the halo. By fitting the spatially binned Ly$\alpha$ spectra with a collection of radiative transfer galactic wind models, we infer a mean outflow expansion velocity of $\approx$ 211 km/s, with higher values preferentially found on both sides of the system's major axis. The velocity of the outflow is validated with the blueshift of low ionization metal absorption lines in the spectra of the central galaxies. We also identify a faint ($M_{1500}\approx -16.7$) companion detected in both Ly$\alpha$ and continuum, whose properties are in agreement with a predicted population of satellite galaxies that contribute to the extended Ly$\alpha$ emission. Finally, we briefly discuss the impact of the interaction between the central galaxies on the properties of the halo and the possibility of in situ fluorescent Ly$\alpha$ production.

Enormous Ly$\alpha$ nebulae (ELANe) around quasars have provided unique insights into the formation of massive galaxies and their associations with super-massive black holes since their discovery. However, their detection remains highly limited. This paper introduces a systematic search for extended Ly$\alpha$ emission around 8683 quasars at $z=$2.34-3.00 using a simple but very effective broad-band $gri$ selection based on the Third Public Data Release of the Hyper Suprime-Cam Subaru Strategic Program. Although the broad-band selection detects only bright Ly$\alpha$ emission ($\gtrsim1\times10^{-17}$ erg~s$^{-1}$cm$^{-2}$arcsec$^{-2}$) compared with narrow-band imaging and integral field spectroscopy, we can apply this method to far more sources than such common approaches. We first generated continuum $g$-band images without contributions from Ly$\alpha$ emission for host and satellite galaxies using $r$- and $i$-bands. Then, we established Ly$\alpha$ maps by subtracting them from observed $g$-band images with Ly$\alpha$ emissions. Consequently, we discovered extended Ly$\alpha$ emission (with masked area $>40$ arcsec$^2$) for 7 and 32 out of 366 and 8317 quasars in the Deep and Ultra-deep (35 deg$^2$) and Wide (890 deg$^2$) layers, parts of which may be potential candidates of ELANe. However, none of them seem to be equivalent to the largest ELANe ever found. We detected higher fractions of quasars with large nebulae around more luminous or radio-loud quasars, supporting previous results. Future applications to the forthcoming big data from the Vera C. Rubin Observatory will help us detect more promising candidates. The source catalogue and obtained Ly$\alpha$ properties for all the quasar targets are accessible as online material.

When the first CubeSats were launched nearly two decades ago, few people believed that the miniature satellites would likely prove to be a useful scientific tool. Skeptics abounded. However, the last decade has seen the highly successful implementation of space missions that make creative and innovative use of fast-advancing CubeSat and small satellite technology to carry out important science experiments and missions. Several projects now have used CubeSats to obtain first-of-their-kind observations and findings that have formed the basis for high-profile engineering and science publications, thereby establishing without doubt the scientific value and broad utility of CubeSats. In this paper, we describe recent achievements and lessons learned from a representative selection of successful CubeSat missions with a space weather focus. We conclude that these missions were successful in part because their limited resources promoted not only mission focus but also appropriate risk-taking for comparatively high science return. Quantitative analysis of refereed publications from these CubeSat missions and several larger missions reveals that mission outcome metrics compare favorably when publication number is normalized by mission cost or if expressed as a weighted net scientific impact of all mission publications.

In this paper, we extend the work of Freundlich et al. 2020 who showed how to obtain a Dekel-Zhao density profile with mass dependent shape parameters in the case of galaxies. In the case of Freundlich et al. 2020, the baryonic dependence was obtained using the NIHAO set of simulations. In our case, we used simulations based on a model of ours. Following Freundlich et al. 2020, we obtained the dependence from baryon physics of the two shape parameters, obtaining in this way a mass dependent Dekel-Zhao profile describing the dark matter profiles from galaxies to clusters of galaxies. The extension to the Dekel-Zhao mass dependent profile to clusters of galaxies is the main result of the paper. In the paper, we show how the Dekel-Zhao mass dependent profile gives a good description of the density profiles of galaxies, already shown by Freundlich et al. 2020, but also to a set of clusters of galaxies.

We consider the optical properties of the solar gravitational lens (SGL) treating the Sun as a massive compact body. Using our previously developed wave-optical treatment of the SGL, we convolve it with a thin-lens representing an optical telescope, and estimate the power spectral density and associated photon flux at individual pixel locations on the image sensor at the focal plane of the telescope. We also consider the solar corona, which is the dominant noise source when imaging faint objects with the SGL. We evaluate the signal-to-noise ratio at individual pixels as a function of wavelength. To block out the solar light, we contrast the use of a conventional internal coronagraph with a Lyot-stop to an external occulter (i.e., starshade). An external occulter, not being a subject to the diffraction limit of the observing telescope, makes it possible to use small telescopes (e.g., $\sim 40$~cm) for spatially and spectrally resolved imaging with the SGL in a broad range of wavelengths from optical to mid-infrared (IR) and without the substantial loss of optical throughput that is characteristic to internal devices. Mid-IR observations are especially interesting as planets are self-luminous at these wavelengths, producing a strong signal, while there is significantly less noise from the solar corona. This part of the spectrum contains numerous features of interest for exobiology and biosignature detection. We develop tools that may be used to estimate instrument requirements and devise optimal observing strategies to use the SGL for high-resolution, spectrally resolved imaging, ultimately improving our ability to confirm and study the presence of life on a distant world.

Context. An essential facet of turbulence is the space-time intermittency of the cascade of energy that leads to coherent structures of high dissipation. Aims. In this work, we attempt to investigate systematically the physical nature of the intense dissipation regions in decaying isothermal magnetohydrodynamical (MHD) turbulence. Methods. We probe the turbulent dissipation with grid based simulations of compressible isothermal decaying MHD turbulence. We take unprecedented care at resolving and controlling dissipation: we design methods to locally recover the dissipation due to the numerical scheme. We locally investigate the geometry of the gradients of the fluid state variables. We develop a method to assess the physical nature of the largest gradients in simulations and to estimate their travelling velocity. Finally we investigate their statistics. Results. We find that intense dissipation regions mainly correspond to sheets: locally, density, velocity and magnetic fields vary primarily across one direction. We identify these highly dissipative regions as fast/slow shocks or Alfv{\'e}n discontinuities (Parker sheets or rotational discontinuities). On these structures, we find the main deviation from 1D planar steady-state is mass loss in the plane of the structure. We investigate the effect of initial conditions which yield different imprints at early time on the relative distributions between these four categories. However, these differences fade out after about one turnover time, when they become dominated by weakly compressible Alfv{\'e}n discontinuities. We show that the magnetic Prandtl number has little influence on the statistics of these discontinuities, but it controls the Ohmic vs viscous heating rates within them. Finally, we find the entrance characteristics of the structures (such as entrance velocity and magnetic pressure) are strongly correlated. Conclusions. These new methods allow to consider developed compressible turbulence as a statistical collection of intense dissipation structures. This can be used to post-process 3D turbulence with detailed 1D models apt for comparison with observations. It could also reveal useful as a framework to formulate new dynamical properties of turbulence.

Radial velocities and direct imaging observations of exoplanets suggest that the frequency of giant planets may decrease for intermediate-mass stars (2.5-8Msun) and it is unclear what is the key mechanism that may hinder their formation. From a theoretical point of view, planet formation around intermediate-mass stars may take longer times, which together with fast migration and efficient photoevaporation may prevent planet formation in these environments. In this letter, we investigate the temporal evolution of the radial drift for dust particles in disks when stellar evolution is taken into account. We demonstrate that the particles drift velocity around intermediate-mass stars sharply increases after 1-2 Myr thus as to potentially become a difficult barrier to overcome in the first steps of planet formation. This high radial drift could explain the lack of disk detections around intermediate-mass stars older than 3-4 Myr, as opposed to low-mass stars (<2.5sun) where the drift may not be the most impacting factor for the disk evolution. Future high-resolution images of these disks can help us explain why planets around intermediate-mass stars may be rare, and whether the role of efficient dust radial drift is to hinder or not planet formation around intermediate-mass stars.

We analyse the clustering of matter on large scales in an extension of the concordance model that allows for spatial curvature. We develop a consistent approach to curvature and wide-angle effects on the galaxy 2-point correlation function in redshift space. In particular we derive the Alcock-Paczynski distortion of $f\sigma_{8}$, which differs significantly from empirical models in the literature. A key innovation is the use of the `Clustering Ratio', which probes clustering in a different way to redshift-space distortions, so that their combination delivers more powerful cosmological constraints. We use this combination to constrain cosmological parameters, without CMB information. In a curved Universe, we find that $\Omega_{{\rm m}, 0}=0.26\pm 0.04$ (68\% CL). When the clustering probes are combined with low-redshift background probes -- BAO and SNIa -- we obtain the first CMB-independent constraint on curvature: $\Omega_{K,0}=0.0041\,_{-0.0504}^{+0.0500}$. We find no Bayesian evidence that the flat concordance model can be rejected. In addition we show that the sound horizon at decoupling is $r_{\rm d} = 144.57 \pm 2.34 \; {\rm Mpc}$, in agreement with its measurement from CMB anisotropies. As a consequence, the late-time Universe is compatible with flat $\Lambda$CDM and a standard sound horizon, leading to a small value of $H_{0}$, {\em without} assuming any CMB information. Clustering Ratio measurements produce the only low-redshift clustering data set that is not in disagreement with the CMB, and combining the two data sets we obtain $\Omega_{K,0}= -0.023 \pm 0.010$.

The third Gaia data release provides low-resolution spectra for around 200 million sources. It is expected that a sizeable fraction of them contain a white dwarf (WD), either isolated, or in a binary system with a main-sequence (MS) companion, i.e. a white dwarf-main sequence (WDMS) binary. Taking advantage of a consolidated Random Forest algorithm used in the classification of WDs, we extend it to study the feasibility of classifying Gaia WDMS binary spectra. The Random Forest algorithm is first trained with a set of synthetic spectra generated by combining individual WD and MS spectra for the full range of effective temperatures and surface gravities. Moreover, with the aid of a detailed population synthesis code, we simulate the Gaia spectra for the above mentioned populations. For evaluating the performance of the models, a set of metrics are applied to our classifications. Our results show that for resolving powers above ~300 the accuracy of the classification depends exclusively on the SNR of the spectra, while below that value the SNR should be increased as the resolving power is reduced to maintain a certain accuracy. The algorithm is then applied to the already classified SDSS WDMS catalogue, revealing that the automated classification exhibits an accuracy comparable (or even higher) to previous classification methods. Finally, we simulate the {\it Gaia} spectra, showing that our algorithm is able to correctly classify nearly 80% the synthetic WDMS spectra. Our algorithm represents a useful tool in the analysis and classification of real Gaia WDMS spectra. Even for those spectra dominated by the flux of the MS stars, the algorithm reaches a high degree of accuracy (60%).

Switchbacks are sudden, large radial deflections of the solar wind magnetic field, widely revealed in interplanetary space by the Parker Solar Probe. The switchbacks' formation mechanism and sources are still unresolved, although candidate mechanisms include Alfv\'enic turbulence, shear-driven Kelvin-Helmholtz instabilities, interchange reconnection, and geometrical effects related to the Parker spiral. This Letter presents observations from the Metis coronagraph onboard Solar Orbiter of a single large propagating S-shaped vortex, interpreted as first evidence of a switchback in the solar corona. It originated above an active region with the related loop system bounded by open-field regions to the East and West. Observations, modeling, and theory provide strong arguments in favor of the interchange reconnection origin of switchbacks. Metis measurements suggest that the initiation of the switchback may also be an indicator of the origin of slow solar wind.

In this study, we report on a detailed single pulse analysis of the radio emission from a rotating radio transient (RRAT) J1918$-$0449 which is the first RRAT discovered with the five hundred meter aperture spherical radio telescope (FAST). The sensitive observations were carried out on 30 April 2021 using the FAST with a central frequency of 1250 MHz and a short time resolution of 49.152 $\mu$s, which forms a reliable basis to probe single pulse emission properties in detail. The source was successively observed for around 2 hours. A total of 83 dispersed bursts with significance above 6$\sigma$ are detected over 1.8 hours. The source's DM and rotational period are determined to be 116.1$\pm$0.4 \pcm \ and 2479.21$\pm$0.03 ms, respectively. The share of registered pulses from the total number of observed period is 3.12\%. No underlying emission is detected in the averaged off pulse profile. For bursts with fluence larger than 10 Jy ms, the pulse energy follows a power-law distribution with an index of $-3.1\pm0.4$, suggesting the existence of bright pulse emission. We find that the distribution of time between subsequent pulses is consistent with a stationary Poisson process and find no evidence of clustering over the 1.8 h observations, giving a mean burst rate of one burst every 66 s. Close inspection of the detected bright pulses reveals that 21 pulses exhibit well-defined quasi-periodicities. The subpulse drifting is present in non-successive rotations with periodicity of $2.51\pm0.06$ periods. Finally, possible physical mechanisms are discussed.

Quasi-periodic eruptions (QPEs) are found in the center of five galaxies, where a tidal disruption event (TDE)-like event is also reported in GSN 069 that happened a couple of years before the QPEs. We explain the connection of these phenomena based on a model of a highly eccentric white-dwarf (WD)-$10^{4-6}M_{\odot}$ massive black hole (MBH) binary formed by the Hill mechanism. In this system, the tidal induced internal oscillation of WD can heat the WD envelope thereby induces the tidal nova and inflates the WD envelope, which can be captured by the MBH and form a TDE. The tidal stripping of the surviving WD in the eccentric orbit can produce the QPEs. We also apply this model to the other four QPE sources. Based on the estimated fallback rate, we find that the remaining time after the QPE-observed time for these QPEs is only around 1-2 years based on our simple model estimation, and then the WD will be fully disrupted. We also show that the accretion rate can be much higher than Eddington accretion rate in final stage of these QPE sources. The peak frequency of spectral energy distribution of disk stays in soft X-ray band ($\sim 0.1-1$ keV), which is consistent with observational results.

There is an unceasing incoming flux of extraterrestrial materials reaching the Earth's atmosphere. Some of these objects produce luminous columns when they ablate during the hypersonic encounter with air molecules. A few fireballs occur each year bright enough to be detected from space. The source of these events is still a matter of debate, but it is generally accepted that they are of sporadic origin. We studied the NASA-JPL Center for NEOs Studies (CNEOS) fireball database to infer the dynamic origin of large bolides produced by meter-sized projectiles that impacted our planet. These likely meteorite-dropping events were recorded by the US Government satellite sensors. We estimated the false-positive rate and analyzed the time evolution of multiple orbit dissimilarity criteria concerning potential associations with near-Earth objects and meteoroid streams. We found that at least 16% of the large bolides could be associated with meteoroid streams, about 4% are likely associated with near-Earth asteroids, and 4% may be linked to near-Earth comets. This implies that a significant fraction of meter-sized impactors producing large bolides may have an asteroidal or cometary origin. In addition, we found at least three bolides having hyperbolic orbits with high tensile strength values. Meter-sized meteoroids of interstellar origin could be more common than previously thought, representing about 1% of the flux of large bolides. The inferred bulk physical properties suggest that the interstellar medium could bias these projectiles towards high strength rocks with the ability to survive prolonged exposure to the harsh interstellar space conditions.

We searched for diffuse interstellar bands (DIBs) in the 0.91$<\lambda<$1.33 $\mu$m region by analyzing the near-infrared (NIR) high-resolution ($R=20,000$ and 28,000) spectra of 31 reddened early-type stars ($0.04<E(B-V)<4.58$) and an unreddened reference star. The spectra were collected using the WINERED spectrograph, which was mounted on the 1.3-m Araki telescope at Koyama Astronomical Observatory, Japan, in 2012 to 2016, and on the 3.58-m New Technology Telescope at La Silla Observatory, Chile, in 2017 to 2018. We detected 54 DIBs, 34 of which are newly detected by this study, and eight DIB candidates. Using this updated list, the DIB distributions in a wide wavelength range from optical to NIR are investigated. The full width at half maximum values of NIR DIBs are found to be narrower than those of optical DIBs on average, which suggests that DIBs at longer wavelengths tend to be caused by larger molecules. Assuming that the larger carriers are responsible for the DIBs at longer wavelengths and have the larger oscillator strengths, we found that the total column densities of DIB carriers tend to decrease with increasing DIB wavelength. The candidate molecules and ions for NIR DIBs are also discussed.

We present a model for one cycle of a classical nova outburst based on a self-consistent wind mass loss accelerated by the gradient of radiation pressure, i.e., the so-called optically thick winds. Evolution models are calculated by a Henyey code for a 1.0 $M_\odot$ white dwarf (WD) with a mass accretion rate of $5 \times 10^{-9}~M_\odot$ yr$^{-1}$. The outermost part of hydrogen-rich envelope is connected to a steadily moving envelope when optically thick winds occur. We confirm that no internal shock waves occur at the thermonuclear runaway. The wind mass loss rate reaches a peak of $1.4 \times 10^{-4}~M_\odot$ yr$^{-1}$ at the epoch of the maximum photospheric expansion, where the photospheric temperature decreases to $\log T_{\rm ph}$ (K)=3.90. Almost all of the accreted mass is lost in the wind. The nuclear energy generated in hydrogen burning is lost in a form of photon emission (64 %), gravitational energy (lifting-up the wind matter against the gravity, 35 %), and kinetic energy of the wind (0.23 %). A classical nova should be very bright in a far-UV (100 - 300 \AA) band, during a day just after the onset of thermonuclear runaway ($\sim$25 d before the optical maximum). In the decay phase of the nova outburst, the envelope structure is very close to that of a steady state solution.

We propose that primordial black hole (PBH) binary systems can work as standard timers in tracking the evolution of the Universe. Through gravitational waves from monochromatic PBH binaries, the probability distribution on major axis and eccentricity from the same redshift is obtained. By studying the dynamical evolution of PBH binaries from the initial probability distribution to observed redshifted ones, the redshift-time calibration can be extracted, which can constrain cosmological models. A general formalism of the standard timer is further concluded based on the evolution of statistical distribution in dynamical systems.

The third Antarctic Survey Telescope array instrument {at} Dome A in Antarctica, the AST3-3 telescope, has been in commissioning from March 2021. We deployed {AST3-3} at the Yaoan astronomical station in Yunnan Province for an automatic time-domain survey and follow-up observations with an optimised observation and protection system. The telescope system of AST3-3 is similar to that of AST3-1 and AST3-2, except that it is equipped with a {14K~$ \times$~10K} QHY411 CMOS camera. AST3-3 has a field of view of $1.65^\circ \times 1.23^\circ$ and {is currently using} the $g$ band filter. During commissioning at Yaoan, AST3-3 aims to conduct an extragalactic transient survey, coupled with prompt follow-ups of opportunity targets. {In this paper, we present the architecture of the AST3-3 automatic observation system.} We demonstrate the data processing of observations by representatives SN 2022eyw and GRB 210420B.

A recent analysis of the 100 pc white dwarf sample in the SDSS footprint demonstrated for the first time the existence of a well defined ultracool -- or IR-faint -- white dwarf sequence in the Hertzsprung-Russell diagram. Here we take advantage of this discovery to enlarge the IR-faint white dwarf sample threefold. We expand our selection to the entire Pan-STARRS survey footprint as well as the Montreal White Dwarf Database 100 pc sample, and identify 37 candidates with strong flux deficits in the optical. We present follow-up Gemini optical spectroscopy of 30 of these systems, and confirm all of them as IR-faint white dwarfs. We identify an additional set of 33 objects as candidates based on their colors and magnitudes. We present a detailed model atmosphere analysis of all 70 newly identified IR-faint white dwarfs together with 35 previously known objects reported in the literature. We discuss the physics of model atmospheres and show that the key physical ingredient missing in our previous generation of model atmospheres was the high-density correction to the He-minus free-free absorption coefficient. With new model atmospheres calculated for the purpose of this analysis, we now obtain significantly higher effective temperatures and larger stellar masses for these IR-faint white dwarfs than the Teff and M values reported in previous analyses, thus solving a two decade old problem. In particular, we identify in our sample a group of ultramassive white dwarfs in the Debye cooling phase with stellar parameters never measured before.

CONTEXT: V1309 Sco is the only certain noncompact stellar merger, due to its indisputable preoutburst light curve matching that of a contact binary of almost equal mass stars. Therefore, anything that can be deduced from the existing observations serves as benchmark constraints for models. AIMS: We present some observational evidences to guide future hydrodynamical simulations and common envelope studies. METHODS: Using archive spectra taken at high and mid spectral resolution during the V1309 Sco outburst and late decline, together with the inferential methods we developed to study nova ejecta through panchromatic high resolution spectroscopic follow ups, we constrain the physical state, structure, dynamics and geometry of the transient originated in the stellar merger. RESULTS: We found that the emitted spectra arise from two distinct contributions: matter expelled during the 2008 outburst and circumbinary gas produced during historic mass loss episodes. These two components likely have orthogonal geometry with the 2008 mass loss displaying a dust-laden bipolar ejecta produced by a time limited rapidly accelerating wind and the circumbinary gas having a donut-like shape. A central source powers them both, having produced a fluorescent light pulse, but we cannot precisely determine the time it started or its spectral energy distribution. We can, however, place its upper energy cutoff at about 54 eV and the bulk of its emission at $<$20 eV. We also know that the central source turned off within months from the outburst and before the ejecta turned optically thin.

Radioactive decay of unstable atomic nuclei leads to liberation of nuclear binding energy in the forms of gamma-ray photons and secondary particles (electrons, positrons); their energy then energises surrounding matter. Unstable nuclei are formed in nuclear reactions, which can occur either in hot and dense extremes of stellar interiors or explosions, or from cosmic-ray collisions. In high-energy astronomy, direct observations of characteristic gamma-ray lines from the decay of radioactive isotopes are important tools to study the process of cosmic nucleosynthesis and its sources, as well as tracing the flows of ejecta from such sources of nucleosynthesis. These observations provide a valuable complement to indirect observations of radioactive energy deposits, such as the measurement of supernova light in the optical. Here we present basics of radioactive decay in astrophysical context, and how gamma-ray lines reveal details about stellar interiors, about explosions on stellar surfaces or of entire stars, and about the interstellar-medium processes that direct the flow and cooling of nucleosynthesis ashes once having left their sources. We address radioisotopes such as $^{56}$Ni, $^{44}$Ti, $^{26}$Al, $^{60}$Fe, $^{22}$Na, $^{7}$Be, and also how characteristic gamma-ray emission from the annihilation of positrons is connected to these.

We present the results of a spectral-line survey of the W51e1/e2 star-forming region at 68-88 GHz. 79 molecules and their isotopologues were detected, from simple diatomic or triatomic molecules, such as SO, SiO, and CCH, to complex organic compounds, such as CH$_3$OCH$_3$ or CH$_3$COCH$_3$. A number of lines that are absent from the Lovas list of molecular lines observed in space were detected, and most of these were identified. A significant number of the detected molecules are typical for hot cores. These include the neutral molecules HCOOCH$_3$, CH$_3$CH$_2$OH, CH$_3$COCH$_3$ etc, which are currently believed to exist in the gas phase only in hot cores and shock-heated gas. In addition, vibrationally excited C$_4$H and HC$_3$N lines with upper-level energies of several hundred Kelvins were found. Such lines can arise only in hot gas with temperatures on the order of 100 K or higher. Apart from neutral molecules, various molecular ions were also detected. Some of these (HC$^{18}$O$^+$, H$^{13}$CO$^+$, and HCS$^+$) usually exist in molecular clouds with high visual extinctions. Potential formation pathways of complex organic molecules (COMs) and of hydrocarbons, along with nitriles, are considered. These formation routes are first discussed in the context of laboratory experiments elucidating the synthesis of organic molecules in interstellar ices in cold molecular clouds, followed by sublimation into the gas phase in the hot core stage. Thereafter, we discuss the predominant formation of hydrocarbons and their nitriles in the gas phase through bimolecular neutral-neutral reactions.

The interpretation of molecular-line data in terms of hydro dynamical simulations of planet-disk interactions fosters new hopes for the indirect detection of protoplanets. In a model-independent approach, embedded protoplanets should be found at the roots of abrupt Doppler flips in velocity centroid maps. However, the largest velocity perturbation known for an unwarped disk, in the disk of HD100546, leads to a conspicuous Doppler flip that coincides with a thick dust ring, in contradiction with an interpretation in terms of a >~ 1Mjup body. Here we present new ALMA observations of the 12CO(2-1) kinematics in HD\,100546, with a factor of two finer angular resolutions. We find that the disk rotation curve is consistent with a central mass 2.1 < M* /Msun < 2.3, and that the blue-shifted side of the Doppler flip is due to vertical motions, reminiscent of the disk wind proposed previously from blue-shifted SO lines. We tentatively propose a qualitative interpretation in terms of a surface disturbance to the Keplerian flow, i.e. a disk eruption, driven by an embedded outflow launched by a ~10Mearth body. Another interpretation involves a disk-mass-loading hot-spot at the convergence of an envelope accretion streamer.

In this work we present a novel compute framework for reconstructing Faraday depth signals from noisy and incomplete spectro-polarimetric radio datasets. This framework is based on a compressed-sensing approach that addresses a number of outstanding issues in Faraday depth reconstruction in a systematic and scaleable manner. We apply this framework to early-release data from the MeerKAT MIGHTEE polarisation survey.

BL Lac object PKS 0735+178 showed some complex multi-wavelength variation phenomena in the previous studies, especially for its color behavior. Bluer-when-brighter, redder-when-brighter and achromatic behavior were all found to be possible long-term trends of PKS 0735+178. In this work, we collected the long-term multi-wavelength data of PKS 0735+178, and also performed a multi-color optical monitoring on intraday timescale. The intraday variability was detected on one night. On the long timescale, a possible 22-day time lag was found between the $R$ and $\gamma$-ray bands. The results of the cross-correlation analysis exhibited strong correlations between various optical bands on both intraday and long timescale. However, only a mild correlation was found between the long-term $\gamma$-ray and $R$-band light curves, which could be interpreted by the different emission mechanisms of the $\gamma$-ray and optical emissions. PKS 0735+178 showed a significant harder-when-brighter in the $\gamma$-ray band, which is consistent with the observed optical bluer-when-brighter trend on both long-term and intraday timescales. We found that the HWB and BWB trends will be enhanced during the active states, especially for the historical low state. Such phenomenon indicates a special activity-dependent color behavior of PKS 0735+178, and it could be well interpreted by the jet emission model.

We present an assessment of the accuracy of common operations performed in $21$-cm spectral line stacking experiments. To this end, we generate mock interferometric data surveying the 21-cm emission at frequency $1310<\nu<1420$ MHz ($0.005<z<0.084$) and covering an area $\sim 6$ deg$^2$ of the sky, mimicking the observational characteristics of real MeerKAT observations. We find that the primary beam correction accounts for just few per cent ($\sim8\%$ at 0 primary beam power, $\sim 3\%$ at 0.6 primary beam power) deviations from the true $M_{\rm HI}$ signal, and that weighting schemes based on noise properties provide unbiased results. On the contrary, weighting schemes based on distance can account for significant systematic mass differences when applied to a flux-limited sample ($\Delta M_{\rm HI}\sim 40-50\%$ in the studied case). We find no significant difference in the final $\braket{M_{\rm HI}}$ obtained when spectroscopic redshift uncertainties are accounted for in the stacking procedure ($ \Delta z\sim 0.00035$, i.e. $\Delta v \sim 100\,{\rm km\, s}^{-1}$). We also present a novel technique to increase the effective size of the galaxy sample by exploiting the geometric symmetries of galaxy cubelets, potentially enhancing the SNR by a factor $\sim\sqrt{2}$ when analyzing the final stacked spectrum (a factor 4 in a cubelet). This procedure is found to be robustly unbiased, while efficiently increasing the SNR, as expected. We argue that an appropriate framework employing detailed and realistic simulations is required to exploit upcoming datasets from SKA pathfinders in an accurate and reliable manner.

We fit the multi-band lightcurves of 40 fast blue optical transients (FBOTs) with the magnetar engine model. The mass of the FBOT ejecta, the initial spin period and polar magnetic field of the FBOT magnetars are respectively constrained to $M_{\rm{ej}}=0.18^{+0.52}_{-0.13}\,M_\odot$, $P_{\rm{i}}=9.4^{+8.1}_{-3.9}\,{\rm{ms}}$, and $B_{\rm p}=7^{+16}_{-5}\times10^{14}\,{\rm{G}}$. The wide distribution of the value of $B_{\rm p}$ spreads the parameter ranges of the magnetars from superluminous supernovae (SLSNe) to broad-line Type Ic supernovae (SNe Ic-BL; some are observed to be associated with long-duration gamma-ray bursts), which are also suggested to be driven by magnetars. Combining FBOTs with the other transients, we find a strong universal anti-correlation as $P_{\rm{i}}\propto{M_{\rm{ej}}^{-0.45}}$, indicating them could share a common origin. To be specific, it is suspected that all of these transients originate from collapse of extreme-stripped stars in close binary systems, but with different progenitor masses. As a result, FBOTs distinct themselves by their small ejecta masses with an upper limit of ${\sim}1\,M_\odot$, which leads to an observational separation in the rise time of the lightcurves $\sim12\,{\rm d}$. In addition, the FBOTs together with SLSNe can be separated from SNe Ic-BL by an empirical line in the $M_{\rm peak}-t_{\rm rise}$ plane corresponding to an energy requirement of a mass of $^{56}$Ni of $\sim0.3M_{\rm ej}$, where $M_{\rm peak}$ is the peak absolute magnitude of the transients and $t_{\rm rise}$ is the rise time.

The star formation efficiency (SFE) has been shown to vary across different environments, particularly within galactic starbursts and deep within the bulges of galaxies. Various quenching mechanisms may be responsible, ranging from galactic dynamics to feedback from active galactic nuclei (AGN). Here, we use spatially-resolved observations of warm ionised gas emission lines from the imaging Fourier transform spectrograph SITELLE at the Canada-France-Hawaii Telescope (CFHT) and cold molecular gas (CO(2-1)) from the Atacama Large Millimeter/sub-millimeter Array (ALMA) to study the SFE in the bulge of the AGN-host galaxy NGC 3169. After distinguishing star-forming regions from AGN-ionised regions using emission-line ratio diagnostics, we measure spatially-resolved molecular gas depletion times (\tau_dep = 1/SFE) with a spatial resolution of \approx 100 pc within a galactocentric radius of 1.8 kpc. We identify a star-forming ring located at radii 1.25 \pm 0.6 kpc with an average \tau_dep of 0.3 Gyr. At radii < 0.9 kpc, however, the molecular gas surface densities and depletion times increase with decreasing radius, the latter reaching approximately 2.3 Gyr at a radius \approx 500 pc. Based on analyses of the gas kinematics and comparisons with simulations, we identify AGN feedback, bulge morphology and dynamics as the possible causes of the radial profile of SFE observed in the central region of NGC 3169.

We show that the Southern Crab (aka Hen2-104) presents an auspicious opportunity to study the form and speed of the invisible winds that excavate and shock the lobes of various types of bipolar nebulae associated with close and highly evolved binary stars. A deep three-color image overlay of Hen2-104 reveals that the ionization state of its lobe edges, or "claws", increases steadily from singly to doubly ionized values with increasing wall latitude. This "reverse" ionization pattern is unique among planetary nebulae (and similar objects) and incompatible with UV photoionization from a central source. We show that the most self-consistent explanation for the ionization pattern is shock ionization by a fast (~600 km s^-1) "tapered" stellar wind in which the speed and momentum flux of the wind increase with equatorial latitude. We present a hydrodynamic simulation that places the latitude-dependent form, the knotty walls, and the reverse ionization of the outer lobes of Hen2-104 into a unified context.

We report results of the first detailed spectral and temporal studies of the recently discovered Be/X-ray binary eRASSU J050810.4$-$660653 in LMC based on the data from the SRG/ART-XC, NuSTAR and Swift/XRT instruments obtained in December 2021 - May 2022 in a wide energy range of 0.5-79 keV. Pulsations with the period of $40.5781 \pm 0.0004$ s were found in the source light curve with the pulsed fraction monotonically increasing with the energy. An estimate of the orbital period of $\sim38$ days was obtained based on the long-term monitoring of the system. The source spectrum can be well approximated with a power-law model modified by an exponential cutoff at high energies. The pulse phase-resolved spectroscopy shows a strong variation of spectral parameters depending on the phase of a neutron star rotation. We have not found any features connected with the cyclotron absorption line both in the phase-averaged and phase-resolved spectra of eRASSU J050810.4$-$660653. However, the neutron star magnetic field was estimated around several $10^{13}$ G using different indirect methods. Discovered variations of the hardness ratio over the pulse phase is discussed in terms of physical and geometrical properties of the emitting region.

Large-scale magnetic fields reveal themselves through diffuse synchrotron sources observed in galaxy clusters such as radio halos. Total intensity filaments of these sources have been observed in polarization as well, but only in three radio halos out of about one hundred currently known. In this paper we analyze new polarimetric Very Large Array data of the diffuse emission in the galaxy cluster Abell 523 in the frequency range 1-2 GHz. We find for the first time evidence of polarized emission on scales of ~ 2.5 Mpc. Total intensity emission is observed only in the central part of the source, likely due to observational limitations. To look for total intensity emission beyond the central region, we combine these data with single-dish observations from the Sardinia Radio Telescope and we compare them with multi-frequency total intensity observations obtained with different instruments, including the LOw Frequency ARray and the Murchison Widefield Array. By analysing the rotation measure properties of the system and utilizing numerical simulations, we infer that this polarized emission is associated with filaments of the radio halo located in the outskirts of the system, in the peripheral region closest to the observer.

We introduce a new prescription for the evolution of globular clusters (GCs) during the initial embedded gas phase into a Monte Carlo method. With a simplified version of the Monte Carlo MOCCA code embedded in the AMUSE framework, we study the survival of GCs after the removal of primordial gas. We first test our code and show that our results for the evolution of mass and Lagrangian radii are in good agreement with those obtained with N-body simulations. The Monte Carlo code enables a more rapid exploration of the evolution of systems with a larger number of stars than N-body simulations. We have carried out a new survey of simulations to explore the evolution of globular clusters with up to $N = 500000$ stars for a range of different star formation efficiencies and half-mass radii. Our study shows the range of initial conditions leading to the clusters' dissolution and those for which the clusters can survive this early evolutionary phase.

Two new diagnostics to break the dark sector degeneracy based on the intrinsic shape alignments of group/cluster size dark matter halos are presented. Using the snapshot data from a series of the DUSTGRAIN-pathfinder $N$-body simulations for the Planck $\Lambda$CDM cosmology and three $f(R)$ gravity models with massive neutrinos ($\nu$), we first determine the probability density functions of the alignment angles between the shape orientations of massive halos and the minor principal axes of the local tidal fields. The numerically obtained results turn out to agree very well with the analytic formula derived under the assumption that the anisotropic merging along the cosmic web induces the halo shape alignments. The four cosmologies, which several standard diagnostics failed to discriminate, are found to yield significantly different best-fit values of the single parameter that characterizes the analytic formula. We also numerically determine the spatial cross-correlations between the shape orientations of neighbor group/cluster halos, and find them to be in good agreements with a fitting formula characterized by two parameters, whose best-fit values are found to substantially differ among the four models. It is concluded that our new diagnostics can disentangle the effect of $f(R)$ gravity from that of massive neutrinos on the occurrence of anisotropic merging along the cosmic web, which is believed to originate the intrinsic shape alignments of massive halos.

We investigate the spatial regularity in the distribution of the young stellar population along spiral arms of three late type spiral galaxies: NGC 895, NGC 5474, and NGC 6946. This study is based on an analysis of photometric properties of spiral arms using {\it GALEX} ultraviolet, optical UBVRI, Halpha, and 8mum {IRAC} infrared surface photometry data. Using the Fourier analysis approach, we found features of spatial regularity or quasi-regularity in the distribution of the young stellar population or (and) regular chains of star formation regions in all arms of NGC 895, NGC 5474, and NGC 6946 with characteristic scales of spacing from 350 to 500 pc in different arms, and (or) scales which are multiples of them. These characteristic scales are close to the those found earlier in NGC 628, NGC 6217, and M100.

The transit technique is responsible for the majority of exoplanet discoveries to date. Characterizing these planets involves careful modeling of their transit profiles. A common technique involves expressing the transit duration using a density-like parameter, $\tilde{\rho}$, often called the "circular density." Most notably, the Kepler project -- the largest analysis of transit lightcurves to date -- adopted a linear prior on $\tilde{\rho}$. Here, we show that such a prior biases measurements of impact parameter, $b$, due to the non-linear relationship between $\tilde{\rho}$ and transit duration. This bias slightly favors low values ($b \lesssim 0.3$) and strongly disfavors high values ($b \gtrsim 0.7$) unless transit signal-to-noise ratio is sufficient to provide an independent constraint on $b$, a criterion that is not satisfied for the majority of Kepler planets. Planet-to-star radius ratio, $r$, is also biased due to $r{-}b$ covariance. Consequently, the median Kepler DR25 target suffers a $1.6\%$ systematic underestimate of $r$. We present a techniques for correcting these biases and for avoiding them in the first place.

V392 Persei is a known dwarf nova (DN) that underwent a classical nova eruption in 2018. Here we report ground-based optical, Swift UV and X-ray, and Fermi-LAT \gamma-ray observations following the eruption for almost three years. V392 Per is one of the fastest evolving novae yet observed, with a $t_2$ decline time of 2 days. Early spectra present evidence for multiple and interacting mass ejections, with the associated shocks driving both the \gamma-ray and early optical luminosity. V392 Per entered Sun-constraint within days of eruption. Upon exit, the nova had evolved to the nebular phase, and we saw the tail of the super-soft X-ray phase. Subsequent optical emission captured the fading ejecta alongside a persistent narrow line emission spectrum from the accretion disk. Ongoing hard X-ray emission is characteristic of a standing accretion shock in an intermediate polar. Analysis of the optical data reveals an orbital period of 3.230 \pm 0.003 days, but we see no evidence for a white dwarf (WD) spin period. The optical and X-ray data suggest a high mass WD, the pre-nova spectral energy distribution (SED) indicates an evolved donor, and the post-nova SED points to a high mass accretion rate. Following eruption, the system has remained in a nova-like high mass transfer state, rather than returning to the pre-nova DN low mass transfer configuration. We suggest that this high state is driven by irradiation of the donor by the nova eruption. In many ways, V392 Per shows similarity to the well-studied nova and DN GK Persei.

Photodisintegration is a main energy loss process of ultrahigh-energy cosmic-ray (UHECR) nuclei in intergalactic space at the highest energies. Therefore, it is crucial to understand photodisintegration's systematic uncertainty to simulate the propagation of UHECR nuclei. In this work, we calculated the cross sections using the random phase approximation (RPA) of density functional theory (DFT), a microscopic nuclear model. We calculated the $E1$ strength of 29 nuclei using three different density functionals. We obtained the cross sections of photonuclear reactions, including photodisintegration, with the $E1$ strength. Then, we implemented the cross sections in a cosmic ray propagation code CRPropa. We found that the difference between the RPA calculations and TALYS in CRPropa in the energy spectrum can be more than the statistical uncertainty of UHECR energy spectrum assuming some astrophysical parameters. We also found that the difference of some astrophysical parameters obtained by a combined fit of UHECR energy spectrum and composition data can be more than the uncertainty of the data between the RPA calculations and TALYS assuming a phenomenological model of UHECR sources.

The models currently used for the analysis of the reflection spectra of black holes usually assume a disk with constant ionization and electron density. However, there is some debate on the impact of these assumptions on the estimate of the properties of the sources, in particular when the fits suggest very steep emissivity profiles in the inner part of the accretion disk. In this work, we re-analyze a selected set of high-quality NuSTAR and Suzaku data of Galactic black holes and we fit the reflection component with three different models: RELXILL_NK, in which the ionization parameter and the electron density are constant, RELXILLION_NK, where the electron density is still constant but the ionization profile is described by a power law, and RELXILLDGRAD_NK, where the electron density profile is described by a power law and the ionization profile is calculated self-consistently from the electron density and the emissivity. While RELXILLION_NK may provide slightly better fits, we do not find any substantial difference in the estimate of the properties of the sources among the three models. Our conclusion is that models with constant electron density and ionization parameter are probably sufficient, in most cases, to fit the currently available X-ray data of accreting black holes.

A toy-model is studied, which considers two flat directions meeting at an enhanced symmetry point such that they realise the usual hybrid inflation mechanism. The kinetic term of the waterfall field features a pole at its Planckian vacuum expectation value (VEV). Consequently, after the phase transition which terminates hybrid inflation, the waterfall field never rolls to its VEV. Instead, it drives a period of "kination", where the stiff barotropic parameter of the Universe $w\approx 1/2$ results in a mild spike in the spectrum of primordial gravitational waves, which will be observable by the forthcoming LISA mission.

Gravitational production of massive particles due to cosmic expansion can be significant during the inflationary and reheating period of the Universe. If the particle also has non-gravitational interactions that do not significantly affect its production, numerous observational probes open up, including cosmological probes. In this work, we focus on the gravitational production of light vector bosons that couple feebly to the Standard Model (SM) particles. Due to the very feeble coupling, the light vector bosons never reach thermal equilibrium, and if the Hubble scale at the end of inflation is above $10^8$ GeV, the gravitational production can overwhelm the thermal production via the freeze-in mechanism by many orders of magnitude. As a result, much stronger constraints from the Big Bang Nucleosynthesis (BBN) can be placed on the lifetime and mass of the vector bosons compared to the scenario where only thermal production is considered. As an example, we study the sub-GeV scale dark photons, which couple to the SM only through kinetic mixing, and derive constraints on the mass and kinetic mixing parameter of the dark photon from the photodisintegration effects on the light element abundances relevant at the end of the BBN when the cosmic age was around $10^4$ s.

By applying a relativistic mean-field description of neutron star matter with density dependent couplings, we analyse the properties of two different matter compositions: nucleonic matter with delta baryons and nucleonic matter with hyperons and delta baryons. The delta-meson couplings are allowed to vary within a wide range of values obtained by experimental data, while the hyperon-meson couplings are fitted to hypernuclear properties. Neutron star properties with no deconfinement phase transition are studied. It is verified that many models are excluded because the effective nucleon mass becomes zero before the maximum mass configuration is attained. Hyperon-free with delta-dominated composition compact stars are possible, the deltic stars. It is found that with a convenient choice of parameters the existence of deltic stars with 80% of delta baryons at the center of the star is possible. However, the presence of hyperons lowers the delta baryon fraction to values below 20% at the center and below 30% at 2-3 saturation densities. It is discussed that in the presence of delta baryons, the hyperon softening is not so drastic because deltas couple more strongly to the $\omega$-meson, and the stiffness of the equation of state is determined by the $\omega$-dominance at high densities. The speed of sound reflects very well this behavior. The compactness of the pulsar RX J0720.4-3125 imposes $x_{\sigma\Delta}>x_{\omega\Delta}>1$ and favors $x_{\rho\Delta}>1$.

We propose a novel experiment, the Canfranc Axion Detection Experiment (CADEx), to probe dark matter axions with masses in the range 330-460 $\mu$eV, within the W-band (80-110 GHz), an unexplored parameter space in the well-motivated dark matter window of Quantum ChromoDynamics (QCD) axions. The experimental design consists of a microwave resonant cavity haloscope in a high static magnetic field coupled to a highly sensitive detecting system based on Kinetic Inductance Detectors via optimized quasi-optics (horns and mirrors). The experiment is in preparation and will be installed in the dilution refrigerator of the Canfranc Underground Laboratory. Sensitivity forecasts for axion detection with CADEx, together with the potential of the experiment to search for dark photons, are presented.

In this short note, we briefly comment on the analytical bounds that must be imposed on the parameter space of the Rezzolla-Zhidenko (RZ) metric-parametrization approach introduced in Ref. [1]. We hope this will clarify some of the confusion recently emerged on this issue [2].

In this work, we introduce a class of extended Minimal Theories of Massive Gravity (eMTMG), without requiring a priori that the theory should admit the same homogeneous and isotropic cosmological solutions as the de Rham-Gabadadze-Tolley massive gravity. The theory is constructed as to have only two degrees of freedom in the gravity sector. In order to perform this step we first introduce a precursor theory endowed with a general graviton mass term, to which, at the level of the Hamiltonian, we add two extra constraints as to remove the unwanted degrees of freedom, which otherwise would typically lead to ghosts and/or instabilities. On analyzing the number of independent constraints and the properties of tensor mode perturbations, we see that the gravitational waves are the only propagating gravitational degrees of freedom which do acquire a non-trivial mass, as expected. In order to understand how the effective gravitational force works for this theory we then investigate cosmological scalar perturbations in the presence of a pressureless fluid. We then restrict the whole class of models by imposing the following conditions at all times: 1) it is possible to define an effective gravitational constant, $G_{{\rm eff}}$; 2) the value $G_{\text{eff}}/G_{N}$ is always finite but not always equal to unity (as to allow some non-trivial modifications of gravity, besides the massive tensorial modes); and 3) the square of mass of the graviton is always positive. These constraints automatically make also the ISW-effect contributions finite at all times. Finally we focus on a simple subclass of such theories, and show they already can give a rich and interesting phenomenology.

The collision of two equilibrium ground state solutions of the Schr\"odinger-Poisson (SP) system, in orthogonal states, is proposed as a formation mechanism of mixed state solutions of the SP system with spherical and first dipolar components. The collisions are simulated by solving numerically the SP system for two orthogonal states, considering head-on encounters, and using various mass ratios between the initial configurations with different head-on momentum. The results indicate that the less massive of the configurations pinches-off the more massive one, and redistributes its density along the axis of collision. The averaged in time density of the two states resembles the distribution of matter of bi-state equilibrium configurations with monopolar and dipolar contributions.

The energy deposition in stellar explosions due to the pair annihilation of neutrinos $(\nu\overline{\nu}\rightarrow e^+e^-)$ can energize events such as Type II supernovae, merging neutron stars, gamma ray bursts (GRBs) etc. This neutrino heating can be further enhanced by modifying the background geometry over that of the Newtonian spacetime. However, even then the observed energy in GRBs cannot be achieved. In this paper, we explore if the inclusion of the contribution due to an extra $Z^\prime$ gauge boson in the neutrino pair annihilation process can explain the energy required for a GRB. We compute the expression for the energy deposition coming from the $Z^\prime$ mediated process in Newtonian, Schwarzschild, and Hartle-Thorne spacetime backgrounds. We show that the contribution due to the $Z^\prime$ mediated process can enhance the energy deposition rate significantly even in the Newtonian background. The inclusion of the general relativistic corrections and rotation in the background spacetime near neutron stars can further enhance the energy deposition rates. From the observed energy of GRBs, we obtain constraints on the mass and the gauge coupling of the extra $Z^\prime$ for the three background spacetimes mentioned.