Papers for Wednesday, Apr 27 2022

We explore how diffuse stellar light and globular clusters (GCs) can be used to trace the matter distribution of their host halo using an observational methodology. For this, we use 117 simulated dark matter (DM) haloes from the $(34.4~\rm cMpc)^3$ periodic volume of the E-MOSAICS project. For each halo, we compare the stellar surface brightness and GC projected number density maps to the surface densities of DM and total mass. We find that the dominant structures identified in the stellar light and in the GCs correspond closely with those from the DM and total mass. Our method is unaffected by the presence of satellites and its precision improves with fainter GC samples. We recover tight relations between the profiles of stellar surface brightness and GC number density to those of the DM, suggesting that the profile of DM can be accurately recovered from the stars and GCs ($\sigma\leq0.5~$dex). We quantify the projected morphology of DM, stars and GCs, and find that the stars and GCs are more flattened than the DM. Additionally, the semi-major axes of the distribution of stars and GCs are typically misaligned by $\sim 10~$degrees from that of DM. We demonstrate that deep imaging of diffuse stellar light and GCs can place constraints on the shape, profile and orientation of their host halo. These results extend down to haloes with central galaxies $M_{\star}\geq10^{10}~M_{\odot}$, and the analysis will be applicable to future data from the Euclid, Roman and the Rubin observatories.

Starting from first principles, we study radiative transfer by new feebly-interacting bosons (FIBs) such as axions, axion-like particles (ALPs), dark photons, and others. Our key simplification is to include only boson emission or absorption (including decay), but not scattering between different modes of the radiation field. Based on a given distribution of temperature and FIB absorption rate in a star, we derive explicit volume-integral expressions for the boson luminosity, reaching from the free-streaming to the strong-trapping limit. The latter is seen explicitly to correspond to quasi-thermal emission from a "FIB sphere" according to the Stefan-Boltzmann law. Our results supersede expressions and approximations found in the recent literature on FIB emission from a supernova core and, for radiatively unstable FIBs, provide explicit expressions for the nonlocal ("ballistic") transfer of energy recently discussed in horizontal-branch stars.

The possible existence of primordial black holes in the stellar mass window has received considerable attention because their mergers may contribute to current and future gravitational-wave detections. Primordial black hole mergers, together with mergers of black holes originating from Population~III stars, are expected to dominate at high redshifts ($z\gtrsim 10$). However the primordial black hole merger rate density is expected to rise monotonically with redshift, while Population~III mergers can only occur after the birth of the first stars. Next-generation gravitational-wave detectors such as Cosmic Explorer~(CE) and Einstein Telescope~(ET) can access this distinctive feature in the merger rates as functions of redshift, allowing for a direct measurement of the abundance of the two populations, and hence for robust constraints on the abundance of primordial black holes. We simulate four-months worth of data observed by a CE-ET detector network and perform hierarchical Bayesian analysis to recover the merger rate densities. We find that if the Universe has no primordial black holes with masses of $\mathcal{O}(10M_{\odot})$, the projected upper limit on their abundance $f_{\rm PBH}$ as a fraction of dark matter energy density may be as low as $f_{\rm PBH}\sim \mathcal{O}({10^{-5}})$, about two orders of magnitude lower than current upper limits in this mass range. If instead $f_{\rm PBH}\gtrsim 10^{-4}$, future gravitational wave observations would exclude $f_{\rm PBH}=0$ at the 95\% credible interval.

Transit surveys indicate that there is a deficit of Neptune-sized planets on close-in orbits. If this ``Neptune desert' is entirely cleared out by atmospheric mass loss, then planets at its upper edge should only be marginally stable against photoevaporation, exhibiting strong outflow signatures in tracers like the metastable helium triplet. We test this hypothesis by carrying out a 12-night photometric survey of the metastable helium feature with Palomar/WIRC, targeting seven gas-giant planets orbiting K-type host stars. Eight nights of data are analyzed here for the first time along with reanalyses of four previously-published datasets. We strongly detect helium absorption signals for WASP-69b, HAT-P-18b, and HAT-P-26b; tentatively detect signals for WASP-52b and NGTS-5b; and do not detect signals for WASP-177b and WASP-80b. We interpret these measured excess absorption signals using grids of Parker wind models to derive mass-loss rates, which are in good agreement with predictions from the hydrodynamical outflow code ATES for all planets except WASP-52b and WASP-80b, where our data suggest that the outflows are much smaller than predicted. Excluding these two planets, the outflows for the rest of the sample are consistent with a mean energy-limited outflow efficiency of $\varepsilon = 0.41^{+0.16}_{-0.13}$. Even when we make the relatively conservative assumption that gas-giant planets experience energy-limited outflows at this efficiency for their entire lives, photoevaporation would still be too inefficient to carve the upper boundary of the Neptune desert. We conclude that this feature of the exoplanet population is a pristine tracer of giant planet formation and migration mechanisms.

We aim to identify and characterise binary systems containing red supergiant (RSG) stars in the Small Magellanic Cloud (SMC) using a newly available ultraviolet (UV) point source catalogue obtained using the Ultraviolet Imaging Telescope (UVIT) on board AstroSat. We select a sample of 560 SMC RSGs based on photometric and spectroscopic observations at optical wavelengths and cross-match this with the far-UV point source catalogue using the UVIT F172M filter, finding 88 matches down to m$_{F172M}$=20.3 ABmag, which we interpret as hot companions to the RSGs. Stellar parameters (luminosities, effective temperatures and masses) for both components in all 88 binary systems are determined and we find mass distributions in the ranges 6.1 to 22.3 Solar masses for RSGs and 3.7 to 15.6 Solar masses for their companions. The most massive RSG binary system in the SMC has a combined mass of 32 $\pm$4 M$_\odot$, with a mass ratio (q) of 0.92. By simulating observing biases, we find an intrinsic multipliciy fraction of 18.8 $\pm$ 1.5% for mass ratios in the range 0.3 < q < 1.0 and orbital periods approximately in the range 3 < log P [days] < 8. By comparing our results with those of a similar mass on the main-sequence, we determine the fraction of single stars to be ~20% and argue that the orbital period distribution declines rapidly beyond log P ~ 3.5. We study the mass-ratio distribution of RSG binary systems and find that a uniform distribution best describes the data below 14 M$_\odot$. Above 15 M$_\odot$, we find a lack of high mass-ratio systems.

We present the first model-agnostic analysis of the complete set of Sloan Digital Sky Survey III (BOSS) and -IV (eBOSS) catalogues of luminous red galaxy and quasar clustering in the redshift range $0.2\leq z \leq 2.2$ (10 billion years of cosmic evolution), which consistently includes the baryon acoustic oscillations (BAO), redshift space distortions (RSD) and the shape of the transfer function signatures, from pre- and post-reconstructed catalogues in Fourier space. This approach complements the standard analyses techniques which only focus on the BAO and RSD signatures, and the full-modeling approaches which assume a specific underlying cosmology model to perform the analysis. These model-independent results can then easily be interpreted in the context of the cosmological model of choice. In particular, when combined with $z>2.1$ Ly-$\alpha$ BAO measurements, the clustering BAO, RSD and {\it Shape} parameters can be interpreted within a flat-$\Lambda$CDM model yielding $h=0.6816\pm0.0067$, $\Omega_{\rm m}=0.3001\pm0.0057$ and $10^{9}\times A_s= 2.43\pm0.20$ (or $\sigma_8=0.858\pm0.036$) with a Big Bang Nucleosynthesis prior on the baryon density. Without any external dataset, the BOSS and eBOSS data alone imply $\Omega_{\rm m}=0.2971\pm 0.0061$ and $10^{9}\times A_s=2.39^{+0.24}_{-0.43}$ (or $\sigma_8=0.857\pm0.040$). For models beyond $\Lambda$CDM, eBOSS data alone (in combination with Planck) constrain the sum of neutrino mass to be $\Sigma m_\nu< 0.40$ eV with a BBN prior ($\Sigma m_\nu <0.082$ eV) at 95\% CL, the curvature energy density to $\Omega_\mathrm{k} = -0.022_{-0.038}^{+0.032}$ ($\Omega_\mathrm{k} = 0.0015\pm 0.0016$) and the dark energy equation of state parameter to $w=-0.998_{-0.073}^{+0.085}$ ($w=-1.093_{-0.044}^{+0.048}$) at 68\% CL without a BBN prior.

Forthcoming observational facilities will make the exploration of the early universe routine, likely probing large populations of galaxies at very low metallicities. It will therefore be important to have diagnostics that can solidly identify and distinguish different classes of objects in such low metallicity regimes. We use new photoionisation models to develop diagnostic diagrams involving various nebular lines. We show that combinations of these diagrams allow the identification and discrimination of the following classes of objects in the early universe: PopIII and Direct Collapse Black Holes (DCBH) in pristine environments, PopIII and DCBH embedded in slightly enriched ISM (Z~10^{-5}-10^{-4}), (metal poor) PopII and AGN in enriched ISM. Diagnostics involving rest-frame optical lines (that will be accessible by JWST) have a better discriminatory power, but also rest-frame UV diagnostics can provide very useful information. Interestingly, we find that metal lines such as [OIII] 5007A and CIV 1549A can remain relatively strong (about a factor of 0.1-1 relative H-beta and HeII 1640A, respectively), even in extremely metal poor environments (Z~10^{-5}-10^{-4}), which could be embedding PopIII galaxies and DCBH.

Fast neutrino flavor conversion can occur in core-collapse supernovae or compact binary merger remnants when non-forward collisions are also at play, and neutrinos are not fully decoupled from matter. This work aims to shed light on the conditions under which fast flavor conversion is enhanced or suppressed by collisions. By relying on a neutrino toy model with three angular bins in the absence of spatial inhomogeneities, we consider two angular configurations: one with the distributions of $\nu_e$ and $\bar\nu_e$ being almost isotropic (as expected before neutrino decoupling) and the other one with the distributions of $\nu_e$ and $\bar\nu_e$ being strongly forward peaked (as expected in the free-streaming regime). By including angle-independent, direction-changing collisions, we find that collisions are responsible for an overall enhancement (damping) of flavor conversion in the former (latter) angular configuration. These opposite outcomes are due to the non-trivial interplay between collisions, flavor conversion, and the initial angular distributions of the electron flavors.

We explore how observations relate to the physical properties of the emitting galaxies by post-processing a pair of merging $z\sim2$ galaxies from the cosmological, hydrodynamical simulation NewHorizon using LCARS (Light from Cloudy Added to RAMSES) to encode the physical properties of the simulated galaxy into H$\alpha$ emission line. By carrying out mock observations and analysis on these data cubes we ascertain which physical properties of the galaxy will be recoverable with the HARMONI spectrograph on the European Extremely Large Telescope (ELT). We are able to estimate the galaxy's star formation rate and dynamical mass to a reasonable degree of accuracy, with values within a factor of $1.81$ and $1.38$ of the true value. The kinematic structure of the galaxy is also recovered in mock observations. Furthermore, we are able to recover radial profiles of the velocity dispersion and are therefore able to calculate how the dynamical ratio varies as a function of distance from the galaxy centre. Finally, we show that when calculated on galaxy scales the dynamical ratio does not always provide a reliable measure of a galaxy's stability against gravity or act as an indicator of a minor merger.

We present a comparison of the interstellar medium traced by [CII] (ALMA), and ionised halo gas traced by Lya (MUSE), in and around QSO host galaxies at z~6. To date, 18 QSOs at this redshift have been studied with both MUSE and high-resolution ALMA imaging; of these, 8 objects display a Lya halo. Using datacubes matched in velocity resolution, we compare and contrast the spatial and kinematic information of the Lya halos and the host galaxies' [CII] (and dust-continuum) emission. We find that the Lya halos extend typically 3-30 times beyond the interstellar medium of the host galaxies. The majority of the Lya halos do not show ordered motion in their velocity fields, whereas most of the [CII] velocity fields do. In those cases where a velocity gradient can be measured in Lya, the kinematics do not align with those derived from the [CII] emission. This implies that the Lya emission is not tracing the outskirts of a large rotating disk that is a simple extension of the central galaxy seen in [CII] emission. It rather suggests that the kinematics of the halo gas are decoupled from those of the central galaxy. Given the scattering nature of Lya, these results need to be confirmed with JWST IFU observations that can constrain the halo kinematics further using the non-resonant Ha line.

Post-starburst E+A galaxies are believed to be systems in a rapid transition between major merger starbursts and quiescent ellipticals. Their optical spectrum is dominated by A-type stars, suggesting a significant starburst that was quenched recently. While optical observations of post-starburst galaxies suggest little ongoing star formation, they have been shown to host significant molecular gas reservoirs. This led to the suggestion that gas consumption or expulsion are not required to end the starburst, and that star formation is suppressed by turbulent heating of the molecular gas. We present NOEMA continuum and CO(1-0) observations of 15 post-starburst galaxies, and collect CO measurements in post-starburst galaxies from the literature. Using archival far-infrared observations, we show that the majority of these systems host obscured star formation, with some showing far-infrared emission that is comparable to those of luminous and ultraluminous infrared galaxies. Once far-infrared star formation rates are used, these systems show similar SFR-$M_{\mathrm{H_2}}$ and Kennicutt-Schmidt relations to those observed in star-forming and starburst galaxies. In particular, there is no need to hypothesize star formation quenching by processes other than the consumption of molecular gas by star formation. The combination of optical, far-infrared, and CO observations indicates that some regions within these galaxies have been recently quenched, while others are still forming stars at a high rate. We find little connection between the post-burst age of the optically-thin quenched regions and the star formation rate in the obscured regions. All this calls into question the traditional classification of E+A galaxies.

We study the evolution of dust in a cosmological volume using a hydrodynamical simulation in which the dust production is coupled with the MUPPI (MUlti Phase Particle Integrator) sub-resolution model of star formation and feedback. As for the latter, we keep as reference the model setup calibrated previously to match the general properties of Milky Way like galaxies in zoom-in simulations. However, we suggest that an increase of the star formation efficiency with the local dust to gas ratio would better reproduce the observed evolution of the cosmic star formation density. Moreover, the paucity of quenched galaxies at low redshift demands a stronger role of AGN feedback. We tune the parameters ruling direct dust production from evolved stars and accretion in the inter stellar medium to get scaling relations involving dust, stellar mass and metallicity in good agreement with observations. In low mass galaxies the accretion process is inefficient. As a consequence, they remain poorer in silicate and small grains than higher mass ones. We reproduce reasonably well the few available data on the radial distribution of dust outside the galactic region, supporting the assumption that the dust and gas dynamics are well coupled at galactic scales.

Multi-planet systems are valuable arenas for investigating exoplanet architectures and comparing planetary siblings. TOI-1246 is one such system, with a moderately bright K dwarf ($\rm{V=11.6,~K=9.9}$) and four transiting sub-Neptunes identified by TESS with orbital periods of $4.31~\rm{d},~5.90~\rm{d},~18.66~\rm{d}$, and $~37.92~\rm{d}$. We collected 130 radial velocity observations with Keck/HIRES and TNG/HARPS-N to measure planet masses. We refit the 14 sectors of TESS photometry to refine planet radii ($\rm{2.97 \pm 0.06~R_\oplus},\rm{2.47 \pm 0.08~R_\oplus}, \rm{3.46 \pm 0.09~R_\oplus}$, $\rm{3.72 \pm 0.16~R_\oplus}$), and confirm the four planets. We find that TOI-1246 e is substantially more massive than the three inner planets ($\rm{8.1 \pm 1.1 M_\oplus}$, $\rm{8.8 \pm 1.2 M_\oplus}$, $\rm{5.3 \pm 1.7 M_\oplus}$, $\rm{14.8 \pm 2.3 M_\oplus}$). The two outer planets, TOI-1246 d and TOI-1246 e, lie near to the 2:1 resonance ($\rm{P_{e}/P_{d}=2.03}$) and exhibit transit timing variations. TOI-1246 is one of the brightest four-planet systems, making it amenable for continued observations. It is one of only six systems with measured masses and radii for all four transiting planets. The planet densities range from $\rm{0.70 \pm 0.24}$ to $3.21 \pm 0.44 \rm{g/cm^3}$, implying a range of bulk and atmospheric compositions. We also report a fifth planet candidate found in the RV data with a minimum mass of 25.6 $\pm$ 3.6 $\rm{M_\oplus}$. This planet candidate is exterior to TOI-1246 e with a candidate period of 93.8 d, and we discuss the implications if it is confirmed to be planetary in nature.

Winds from massive stars have recently been deemed promising sites for investigating relativistic particle acceleration. Particularly, the resulting bow shock from the interaction of the winds of runaway stars with interstellar matter has been observed at multiple wavelengths. Here we investigate the O4If star, BD+433654, the bow shock of which is, so far, the only one proven to radiate both thermally and non-thermally at radio frequencies. We also consider NGC7635 as a bow shock candidate and examine its apex for indications of thermal and non-thermal radio emission. We observed both sources with the VLA at 4-8 GHz and 8-12 GHz, and with the Effelsberg telescope at 4-8 GHz. We analysed data from both telescopes individually and combined, obtained their spectral index maps and calculated their Spectral Energy Distributions. We present the first high-resolution maps of radio emission from NGC7635. We find that both emit non-thermal emission in the radio regime, with the clearest evidence for NGC7635. Our results are less conclusive for BD+433654, as the emission from its bow shock becomes weaker and fainter at higher radio frequencies. Our results extend the previous radio results for the BD+433654 bow shock to higher frequencies. Modelling of our data for both sources shows that accelerated electrons at the wind termination shock are a plausible source for the non-thermal radio emission, but energetics arguments suggest that any non-thermal X-ray and $\gamma$-ray emission could be significantly below existing upper limits. Enhanced synchrotron emission from compressed Galactic cosmic rays in the radiative bow shock could also explain the radio emission from the BD+433654 bow shock but not NGC7635. Non-detection of point-like radio emission from BD+433654 puts an upper limit on the mass-loss rate of the star that is lower than values quoted in the literature. [abridged]

The realization of fundamental relations between supermassive black holes and their host galaxies would have profound implications in astrophysics. To add further context to studies of their co-evolution, an investigation is carried out to gain insight as to whether quasars and their hosts at earlier epochs follow the local relation between black hole (BH) mass and stellar velocity dispersion. We use 584 SDSS quasars at 0.2 < z < 0.8 with black hole measurements, and properties of their hosts from the Hyper Suprime-Cam Subaru Strategic Program. An inference of the stellar velocity dispersion is achieved for each based on the stellar mass and size of the host galaxy by using the galaxy mass fundamental plane for inactive galaxies at similar redshifts. In agreement with past studies, quasars occupy an elevated position from the local M_BH-sigma relation, considered as a flattening, while maintaining ratios of M_BH/M* consistent with local values. Based on a forward-modeling of the sample, we demonstrate that an evolving intrinsic M_BH-sigma relation can match the observations. However, we hypothesize that these changes may be a reflection of a non-evolving intrinsic relationship between M_BH and M*. Reassuringly, there are signs of migration onto the local M_BH-sigma for galaxies that are either massive, quiescent or compact. Thus, the majority of the bulges of quasar hosts at high redshift are in a development stage and likely to align with their black holes onto the mass scaling relation at later times.

We interpret the 1.3mm VLBI observations made by the Event Horizon Telescope of the black hole in M87. It is proposed that, instead of being a torus of accreting gas, the ring is a rotating, magnetically-dominated ergomagnetosphere that can transmit electromagnetic angular momentum and energy outward to the disc through a combination of large scale magnetic torque and small scale instabilities. It is further proposed that energy can be extracted by magnetic flux threading the ergosphere through the efficient emission of long wavelength electromagnetic disturbances onto negative energy orbits, when the invariant $B^2-E^2$ becomes negative. In this way, the spinning black hole and its ergosphere not only power the jets but also the ejection disc so as to drive away most of the gas supplied near the Bondi radius. This outflow takes the form of a MHD wind, extending over many decades of radius, with a unidirectional magnetic field, that is collimated by the infalling gas across a magnetopause. This wind, in turn, collimates the relativistic jets and the emission observed from the jet sheath may be associated with a return current. A model for the global flow of mass, angular momentum, energy and current, on scales from the horizon to the Bondi radius, is presented and discussed.

Astrophysical lensing has typically been studied in two regimes: diffractive optics and refractive optics. Diffractive optics is characterized by a perturbative expansion of the Kirchhoff-Fresnel diffraction integral, while refractive optics is characterized by the stationary phase approximation. Previously, it has been assumed that the Fresnel scale, $R_F$ , is the relevant physical scale that separates these two regimes. With the recent introduction of Picard-Lefschetz theory to the field of lensing, it has become possible to generalize the refractive description of discrete images to all wave parameters, and, in particular, exactly evaluate the diffraction integral at all frequencies. In this work, we assess the regimes of validity of refractive and diffractive approximations for a simple one-dimensional lens model through comparison with this exact evaluation. We find that, contrary to previous assumptions, the true separation scale between these regimes is given by $R_F / \sqrt{\kappa}$, where $\kappa$ is the convergence of the lens. Thus, when the lens is strong, refractive optics can hold for arbitrarily small scales. We also argue that intensity variations in diffractive optics are generically small, which has implications for the study of strong diffractive scintillation (DISS).

We present a catalog of solar energetic particle (SEP) events covering solar cycles 22, 23 and 24. We correlate and integrate three existing catalogs based on Geostationary Operational Environmental Satellite (GOES) integral proton flux data. We visually verified and labeled each event in the catalog to provide a homogenized data set. We have identified a total of 342 SEP events of which 246 cross the space weather prediction center (SWPC) threshold of a significant proton event. The metadata consists of physical parameters and observables concerning the possible source solar eruptions, namely flares and coronal mass ejections for each event. The sliced time series data of each event, along with intensity profiles of proton fluxes in several energy bands, have been made publicly available. This data set enables researchers in machine learning (ML) and statistical analysis to understand the SEPs and the source eruption characteristics useful for space weather prediction.

In this paper, wavelet analysis is used to study spectral-timing properties of MAXI J1535-571 observed by Insight-HXMT. The low-frequency quasi-periodic oscillations (QPOs) are detected in nine observations. Based on wavelet analysis, the time intervals with QPO and non-QPO are isolated separately, and the corresponding spectra with QPO and non-QPO are analyzed. We find that the spectra with QPO (hereafter QPO spectra) are softer than those without QPO (hereafter non-QPO spectra) in the hard intermediate state (HIMS). While in the soft intermediate state (SIMS), the QPO spectra are slightly harder. The disk temperature of QPO regime is slightly lower during HIMS, but becomes higher during SIMS. The cutoff energies of QPO spectra and non-QPO spectra do not show significant differences. The flux ratio of the disk to total flux is higher for the time intervals with non-QPO than that of QPO regime. We propose that these differences in the spectral properties between QPO and non-QPO regimes could be explained in the scenario of Lense-Thirring precession, and the reversal of the QPO/non-QPO behavior between HIMS and SIMS may be associated with appearance/disappearance of a type-B QPO which might origin from the precession of the jet.

The process of migration into resonance capture has been well studied for planetary systems where the gravitational potential is generated exclusively by the star and planets. However, massive protoplanetary disks add a significant perturbation to these models. In this paper we consider two limiting cases of disk-induced precession on migrating planets and find that small amounts of precession significantly affect the equilibrium reached by migrating planets. We investigate these effects with a combination of semi-analytic models of the resonance and numerical integrations. We also consider the case of the disk's dispersal, which can excite significant libration amplitude and can cause ejection from resonance for large enough precession rates. Both of these effects have implications for interpreting the known exoplanet population and may prove to be important considerations as the population of well-characterized exoplanet systems continues to grow.

The Fermi Bubbles are giant, gamma-ray emitting lobes emanating from the nucleus of the Milky Way discovered in 1-100 GeV data collected by the Large Area Telescope on board the Fermi Gamma-Ray Space Telescope. Previous work has revealed substructure within the Fermi Bubbles that has been interpreted as a signature of collimated outflows from the Galaxy's super-massive black hole. Here we show that much of the gamma-ray emission associated to the brightest region of substructure - the so-called cocoon - is actually due to the Sagittarius dwarf spheroidal (Sgr dSph) galaxy. This large Milky Way satellite is viewed through the Fermi Bubbles from the position of the Solar System. As a tidally and ram-pressure stripped remnant, the Sgr dSph has no on-going star formation, but we demonstrate that its gamma-ray signal is naturally explained by inverse Compton scattering of cosmic microwave background photons by high-energy electron-positron pairs injected by the dwarf's millisecond pulsar (MSP) population, combined with these objects' magnetospheric emission. This finding suggests that MSPs likely produce significant gamma-ray emission amongst old stellar populations, potentially confounding indirect dark matter searches in regions such as the Galactic Centre, the Andromeda galaxy, and other massive Milky Way dwarf spheroidals.

PY Per has been known as a Z Cam star. Using VSOLJ, VSNET, AAVSO, ASAS-SN and ZTF observations, I found that this object experienced faint states and classified it to be a Z Cam + VY Scl star. Furthermore, the object changed the outburst behavior since 2020 and it has been showing regularly recurring long, bright outbursts with cycles of 110-160 d and short, faint outbursts between them. The 2021 December--2022 January outburst particularly resembled a superoutburst with a gradually fading plateau phase having a duration more than 25 d and two rebrightenings in the fading part. This object has a long orbital period of 0.15468(5) d based on the TESS data and the previous radial-velocity study. The existence of a superoutburst in a system with a long orbital period would be unusual, but not unprecedented. Since the phenomenon may still be ongoing, I draw attention to this object when it emerges from in the morning sky after the solar conjunction.

With the rapid development of deep learning technology, more and more researchers apply it to gravitational wave (GW) data analysis. Previous studies focused on a single deep learning model. In this paper we design an ensemble algorithm combining a set of convolutional neural networks (CNN) for GW signal recognition. The whole ensemble model consists of two sub-ensemble models. Each sub-ensemble model is also an ensemble model of deep learning. The two sub-ensemble models treat data of Hanford and Livinston detectors respectively. Proper voting scheme is adopted to combine the two sub-ensemble models to form the whole ensemble model. We apply this ensemble model to all reported GW events in the first observation and second observation runs (O1/O2) by LIGO-VIRGO Scientific Collaboration. We find that the ensemble algorithm can clearly identify all binary black hole merger events except GW170818. We also apply the ensemble model to one month (August 2017) data of O2. There is no false trigger happens although only O1 data are used for training. Our test results indicate that the ensemble learning algorithms can be used in real-time GW data analysis.

We compare infrared surface brightness fluctuation (IR SBF) distances measured in galaxies that have hosted type Ia supernovae (SNIa) to distances estimated from SNIa light curve fits. We show that the properties of SNIa found in IR SBF hosts are very different from those exploding in Cepheid calibrators, therefore, this is a direct test of systematic uncertainties on estimation of the Hubble constant (Ho) using supernovae. The IR SBF results from Jensen et al. (2021; arXiv:2105.08299) provide a large and uniformly measured sample of IR SBF distances which we directly compare with distances to 25 SNIa host galaxies. We divide the Hubble flow SNIa into sub-samples that best match the divergent supernova properties seen in the IR SBF hosts and Cepheid hosts. We further divide the SNIa into a sample with light curve widths and host masses that are congruent to those found in the SBF-calibrated hosts. We refit the light curve stretch and color correlations with luminosity, and use these revised parameters to calibrate the Hubble flow supernovae with IR SBF calibrators. Relative to the Hubble flow, the average calibrator distance moduli vary by 0.03mag depending on the SNIa subsamples examined and this adds a 1.8% systematic uncertainty to our Hubble constant estimate. Based on the IR SBF calibrators, Ho=74.6$\pm$0.9(stat)$\pm$ 2.7(syst) km/s/Mpc, which is consistent with the Hubble constant derived from supernovae calibrated from Cepheid variables. We conclude that IR SBF provides reliable calibration of SNIa with a precision comparable to Cepheid calibrators, and with a significant saving in telescope time.

We present the results of a search for core-collapse supernova neutrinos, using long-term KamLAND data from 2002 March 9th to 2020 April 25th. We focus on the electron anti-neutrinos emitted from supernovae in the energy range of 1.8--111 MeV. Supernovae will make a neutrino event cluster with the length of $\sim$10 s in the KamLAND data. We find no neutrino clusters and give the upper limit on the supernova rate as to be 0.15 yr$^{-1}$ with a 90% confidence level. The detectable range, which corresponds to a >95% detection probability, is 40--59 kpc and 65--81 kpc for core-collapse supernovae and failed core-collapse supernovae, respectively. This paper proposes to convert the supernova rate obtained by the neutrino observation to the galactic star formation rate. Assuming a modified Salpter-type initial mass function, the upper limit on the galactic star formation rate is <17.5--22.7 (8.1--10.5) $M_{\odot} \mathrm{yr}^{-1}$ with a 90% (68.3%) confidence level.

Here we greatly improve Artificial Intelligence (AI)-generated solar farside magnetograms using data sets of Solar Terrestrial Relations Observatory (STEREO) and Solar Dynamics Observatory (SDO). We modify our previous deep learning model and configuration of input data sets to generate more realistic magnetograms than before. First, our model, which is called Pix2PixCC, uses updated objective functions which include correlation coefficients (CCs) between the real and generated data. Second, we construct input data sets of our model: solar farside STEREO extreme ultraviolet (EUV) observations together with nearest frontside SDO data pairs of EUV observations and magnetograms. We expect that the frontside data pairs provide the historic information of magnetic field polarity distributions. We demonstrate that magnetic field distributions generated by our model are more consistent with the real ones than before in view of several metrics. The averaged pixel-to-pixel CC for full disk, active regions, and quiet regions between real and AI-generated magnetograms with 8 by 8 binning are 0.88, 0.91, and 0.70, respectively. Total unsigned magnetic flux and net magnetic flux of the AI-generated magnetograms are consistent with those of real ones for test data sets. It is interesting to note that our farside magnetograms produce consistent polar field strengths and magnetic field polarities with those of nearby frontside ones for solar cycle 24 and 25. Now we can monitor the temporal evolution of active regions using solar farside magnetograms by the model together with the frontside ones. Our AI-generated Solar Farside Magnetograms (AISFMs) are now publicly available at Korean Data Center (KDC) for SDO.

Evidence for wave-like corrugations are well established in the Milky Way and in nearby disc galaxies. These were originally detected as a displacement of the interstellar medium about the midplane, either in terms of vertical distance or vertical velocity. Over the past decade, similar patterns have emerged in the Milky Way's stellar disc. We investigate how these vertical waves are triggered by a passing satellite. For the first time, we use high-resolution N-body/hydrodynamical simulations to study how the corrugations set up and evolve jointly in the stellar and gaseous discs. We find that the gas corrugations follow the stellar corrugations, i.e. they are initially in phase although, after a few rotation periods (500-700 Myr), the distinct waves separate and thereafter evolve in different ways. The spatial and kinematic amplitudes (and thus the energy) of the corrugations dampen with time, with the gaseous corrugation settling at a faster rate (~800 Myr vs. ~1 Gyr). In contrast, the vertical energy of individual disc stars is fairly constant throughout the galaxy's evolution. This difference arises because corrugations are an emergent phenomenon supported by the collective, ordered motions of co-spatial ensembles of stars. We show that the damping of the stellar corrugations can be understood as a consequence of incomplete phase mixing, while the damping of the gaseous corrugations is a natural consequence of the dissipative nature of the gas. We suggest that the degree of correlation between the stellar and gaseous waves may help to age-date the phenomenon.

We present the results of the high-resolution (R$\sim$60,000) spectroscopic analysis of the star LAMOSTJ045019.27+394758.7 (hereafter J045) from the list of carbon stars of LAMOST DR2. From our analysis, we find that J045 does not exhibit the spectral characteristics of carbon stars. It is found to be a metal-poor ( [Fe/H] = $-$1.05) giant that shows very unusual elemental abundances, particularly for N, Na, V, and Zn. J045 shows ${\alpha}$-elements (Mg, Si, Ca) with near-solar values ($<$[$\alpha$/Fe]$>$ = 0.09) in contrast to Galactic stars that show [$\alpha$/Fe] in the range 0.2 to 0.3 dex. In J045, Sc and Ti are under abundant with [X/Fe] $\le$ $-$0.25. Vanadium gives [V/Fe] = 0.51 and zinc is under-abundant with [Zn/Fe] = $-$0.62. The object exhibits near-solar abundances for Sr, Y, Ba, Pr, and Sm. The La is marginally enhanced, and Ce and Nd are marginally under-abundant in J045. With [Ba/Eu] = $-$0.38, the object falls into the category of neutron-capture rich r-I stars. The estimated abundances of various elements show that the observed abundance pattern is not compatible with the abundances characteristic of Galactic metal-poor stars but matches quite closely with the abundance pattern of Sculptor Dwarf galaxy stars of similar metallicity. Based on the above observational evidences, we suggest that the object is a possible Sculptor Dwarf Galaxy escapee.

We present a comprehensive metallicity analysis of the Sagittarius dwarf spheroidal galaxy (Sgr dSph) using $Pristine\,CaHK$ photometry. We base our member selection on $Gaia$ EDR3 astrometry applying a magnitude limit at $G_{0} = 17.3$, and our population study on the metallicity-sensitive photometry from the $Pristine$ Inner Galaxy Survey (PIGS). Working with photometric metallicities instead of spectroscopic metallicities allows us to cover an unprecedented large area ($\sim 100$ square degrees) of the dwarf galaxy, and to study the spatial distribution of its members as function of metallicity with little selection effects. Our study compares the spatial distributions of a metal-poor population of 9719 stars with [Fe/H] $< -1.3$ and a metal rich one of 30115 stars with [Fe/H] $> -1.0$. The photometric Sgr sample also allows us to assemble the largest sample of 1150 very metal-poor Sgr candidates ([Fe/H] $< -2.0$). By investigating and fitting the spatial properties of the metal-rich and metal-poor population, we find a negative metallicity gradient which extends up to 12 degrees from the Sgr center (or $\sim 5.5$ kpc at the distance of Sgr), the limit of our footprint. We conclude that the relative number of metal-poor stars increases in the outer areas of the galaxy, while the central region is dominated by metal-rich stars. These finding suggest an outside-in formation process and are an indication of the extended formation history of Sgr, which has been affected by the tidal interaction between Sgr and the Milky Way.

The extragalactic background light (EBL) contains all the radiation emitted by nuclear and accretion processes in stars and compact objects since the epoch of recombination. Measuring the EBL density directly is challenging, especially in the near- to far-infrared waveband, mainly due to the zodiacal light foreground. Instead, gamma-ray astronomy offers the possibility to indirectly set limits on the EBL by studying the effects of gamma-ray absorption in the very high energy (VHE:$>$100 GeV) spectra of distant blazars. The High Altitude Water Cherenkov gamma ray observatory (HAWC) is one of the few instruments sensitive to gamma rays with energies above 10 TeV. This offers the opportunity to probe the EBL in the near/mid IR region: $\lambda$ = 1 $\mu$m - 100 $\mu$m. In this study, we fit physically motivated emission models to \textit{Fermi Large Area Telescope (LAT)} GeV data to extrapolate the intrinsic TeV spectra of blazars. We then simulate a large number of absorbed spectra for different randomly generated EBL model shapes and calculate Bayesian credible bands in the EBL intensity space by comparing and testing the agreement between the absorbed spectra and HAWC extragalactic observations of two blazars. The resulting bands are in agreement with current EBL lower and upper limits, showing a downward trend towards higher wavelength values $\lambda>10 \mu$m also observed in previous measurements.

We explore the interplay between supernovae and the ionizing radiation of their progenitors in star forming regions. The relative contributions of these stellar feedback processes are not well understood, particularly on scales greater than a single star forming cloud. We focus predominantly on how they affect the interstellar medium. We re-simulate a 500 pc^2 region from previous work that included photoionization and add supernovae. Over the course of 10 Myr more than 500 supernovae occur in the region. The supernovae remnants cool very quickly in the absence of earlier photoionization, but form much larger and more spherical hot bubbles when photoionization is present. Overall, the photoionization has a significantly greater effect on gas morphology and the sites of star formation. However, the two processes are comparable when looking at their effect on velocity dispersion. When combined, the two feedback processes increase the velocity dispersions by more than the sum of their parts, particularly on scales above 5 pc.

Although the size of the cosmic void is much larger than the size of local scales and much smaller than the size of cosmic scale, it will be shown that they can be considered a very good representative for studying both local and global scales. For this goal, we will first consider the cosmic voids in the cosmic web as interconnected spherical bubbles. We will then show by heuristic calculations that for each cosmic void we obtain different mass densities and cosmological constants that are the same order of magnitude as the entire universe. It will also be shown that the slight difference between the surface tension of the bubbles may be the source of the $H_{\rm 0}$ tension between local and global measurements. As a necessary consequence of this study is that by examining a single cosmic void more seriously, both theoretically and observationally, interesting possible propositions can be made to solve important challenges in physical cosmology at local and global scales.

The Milky Way stellar halo contains relics of ancient mergers that tell the story of our Galaxy's formation. Some of them are identified due to their similarity in energy, actions and chemistry, referred to as the "chemodynamical space", and are often attributed to distinct merger events. It is also known that our Galaxy went through a significant merger event that shaped the local stellar halo during its first Gyr. Previous studies using $N$-body only and cosmological hydrodynamical simulations have shown that such single massive merger can produce several "signatures" in the chemodynamical space, which can potentially be misinterpreted as distinct merger events. Motivated by these, in this work we use a subset of the GASTRO, library which consists of several SPH+$N$-body models of single accretion event in a Milky Way-like galaxy. Here, we study models with orbital properties similar to the main merger event of our Galaxy and explore the implications to known stellar halo substructures. We find that: $i.$ supernova feedback efficiency influences the satellite's structure and orbital evolution, resulting in distinct chemodynamical features for models with the same initial conditions, $ii.$ very retrograde high energy stars are the most metal-poor of the accreted dwarf galaxy and could be misinterpreted as a distinct merger $iii.$ the most bound stars are more metal-rich in our models, the opposite of what is observed in the Milky Way, suggesting a secondary massive merger, and finally $iv.$ our models can reconcile other known substructures to an unique progenitor.

Aims: We investigate the atmospheric response to coronal heating driven by random flows with different characteristic time scales and amplitudes. Methods: We conducted a series of 3D MHD simulations of random driving imposed on a gravitationally stratified model of the solar atmosphere. In order to understand differences between alternating current (AC) and direct current (DC) heating, we considered the effects of changing the characteristic time scales of the imposed velocities. We also investigated the effects of the magnitude of the velocity driving. Results: Complex foot point motions lead to a proliferation of current sheets and energy dissipation throughout the corona. For a given amplitude, DC driving typically leads to a greater rate of energy injection when compared to AC driving. This leads to the formation of larger currents, increased heating rates and higher temperatures in DC simulations. There is no difference in the spatial distribution of energy dissipation across simulations, however, energy release events in AC cases tend to be more frequent and last for less time. Higher velocity driving is associated with larger currents, higher temperatures and the corona occupying a larger fraction of the simulation volume. In all cases, most of heating is associated with small energy release events, which occur much more frequently than large events. Conclusions: When combined with observational results showing a greater abundance of power in low frequency modes, these findings suggest that energy release in the corona is likely to be driven by long time scale motions. In the corona, AC and DC driving will occur concurrently and their effects remain difficult to isolate. The distribution of field line temperatures and the asymmetry of temperature profiles may reveal the frequency and longevity of energy release events and therefore the relative importance of AC and DC heating.

Active region evolution plays an important role in the generation and variability of magnetic fields on the surface of lower main-sequence stars. However, determining the lifetime of active region growth and decay as well as their evolution is a complex task. We aim to test whether the lifetime for active region evolution shows any dependency on the stellar parameters. We identify a sample of stars with well-defined ages via their kinematics. We made use of high-resolution spectra to compute rotational velocities, activity levels, and emission excesses. We use these data to revisit the activity-rotation-age relationship. The time-series of the main optical activity indicators were analysed together with the available photometry by using Gaussian processes to model the stellar activity of these stars. Autocorrelation functions of the available photometry were also analysed. We use the derived lifetimes for active region evolution to search for correlations with the stellar age, the spectral type, and the level of activity. We also use the pooled variance technique to characterise the activity behaviour of our targets. Our analysis confirms the decline of activity and rotation as the star ages. We also confirm that the rotation rate decays with age more slowly for cooler stars and that, for a given age, cooler stars show higher levels of activity. We show that F- and G-type young stars also depart from the inactive stars in the flux-flux relationship. The gaussian process analysis of the different activity indicators does not seem to provide any useful information on active region's lifetime and evolution. On the other hand, active region's lifetimes derived from the light-curve analysis might correlate with the stellar age and temperature. Although we caution the small number statistics, our results suggest that active regions seem to live longer on younger, cooler, and more active stars.

This contribution reviews recent advances in the possible identification of blazars as potential sources of at least some of the very-high-energy neutrinos detected by the IceCube neutrino detector at the South Pole. The basic physical requirements for neutrino production and physics constraints that may be drawn from neutrino - blazar associations are reviewed. Several individual cases of possible associations will be discussed in more detail. It is emphasized that due to $\gamma\gamma$ opacity constraints in efficiently neutrino-producing blazars, an association between X-ray -- soft $\gamma$-ray activity and very-high-energy neutrino production is more naturally expected than a connection between neutrino and high-energy/very-high-energy $\gamma$-ray activity.

In this paper we discuss and confront recent results on metallicity variations in the local interstellar medium, obtained from observations of HII regions and neutral clouds of the Galactic thin disk, and compare them with recent high-quality metallicity determinations of other tracers of the chemical composition of the interstellar medium as B-type stars, classical Cepheids and young clusters. We find that the metallicity variations obtained for these last kinds of objects are consistent with each other and with that obtained for HII regions but significantly smaller than those obtained for neutral clouds. We also discuss the presence of a large population of low-metallicity clouds as the possible origin for large metallicity variations in the local Galactic thin disk. We find that such hypothesis does not seem compatible with: (a) what is predicted by theoretical studies of gas mixing in galactic disks, and (b) the models and observations on the metallicity of high-velocity clouds and its evolution as they mix with the surrounding medium in their fall onto the Galactic plane. We conclude that that most of the evidence favors that the chemical composition of the interstellar medium in the solar neighborhood is highly homogeneous.

In the last 20 years, over 600 impact flashes have been documented on the lunar surface. This wealth of data presents a unique opportunity to study the meteoroid flux of the Earth-Moon environment, and in recent years the physical properties of the impactors. However, other than through serendipitous events, there has not been yet a systematic search and discovery of the craters associated to these events. Such a meteoroid-crater link would allow us to get insight into the crater formation via these live observations of collisions. Here we present the PyNAPLE (Python NAC Automated Pair Lunar Evaluator) software pipeline for locating newly formed craters using the location and epoch of an observed impact flash. We present the first results from PyNAPLE, having been implemented on the 2017-09-27 impact flash. A rudimentary analysis on the impact flash and linked impact crater is also performed, finding that the crater's ejecta pattern indicates an impact angle between 10-30\degree, and although the rim-to-rim diameter of the crater is not resolvable in current LRO NAC images, using crater scaling laws we predict this diameter to be 24.1-55.3~m, and using ejecta scaling predict a diameter of 27.3-37.7~m. We discuss how PyNAPLE will enable large scale analyses of sub-kilometer scale cratering rates and refinement of both scaling laws, and the luminous efficiency.

Compton scattering is a key process shaping spectra formation and accretion flow dynamics in accreting strongly magnetized neutron stars. A strong magnetic field affects the scattering cross section and makes it dependent on photon energy, momentum, and polarization state. Using Monte Carlo simulations, we investigate statistical features of Compton scattering of polarized X-ray radiation in a strong magnetic field. {Our analysis is focused on photon gas behaviour well inside the scattering region.} We take into account the resonant scattering at the fundamental cyclotron frequency, thermal distribution of electrons at the ground Landau level, and bulk velocity of the electron gas. We show that (i) the photons scattered around the cyclotron energy by the electron gas at rest tend to acquire the final energy close to the cyclotron one with a very small dispersion measure; (ii) the redistribution of photons within the Doppler core of cyclotron resonance differs significantly from the complete redistribution; (iii) the efficiency of momentum transfer from photons to the electron gas is affected by the temperature of electron gas both for photons at cyclotron energy and below it; (iv) the momentum transfer from photons to the electron gas of non-zero bulk velocity is more efficient in the case of magnetic scattering.

Protostellar jets and outflows are pointers of star-formation and serve as important sources of momentum and energy transfer to the interstellar medium. Radio emission from ionized jets have been detected towards a number of protostellar objects. In few cases, negative spectral indices and polarized emission have also been observed suggesting the presence of synchrotron emission from relativistic electrons. In this work, we develop a numerical model that incorporates both thermal free-free and non-thermal synchrotron emission mechanisms in the jet geometry. The flux densities include contribution from an inner thermal jet, and a combination of emission from thermal and non-thermal distributions along the edges and extremities, where the jet interacts with the interstellar medium. We also include the effect of varying ionization fraction laterally across the jet. An investigation of radio emission and spectra along the jet shows the dependence of the emission process and optical depth along the line of sight. We explore the effect of various parameters on the turnover frequencies and the radio spectral indices (between 10 MHz and 300 GHz) associated with them.

We study the general-relativistic dynamics of matter being accreted onto and ejected by a magnetised and nonrotating neutron star. The dynamics is followed in the framework of fully general relativistic magnetohydrodynamics (GRMHD) within the ideal-MHD limit and in two spatial dimensions. More specifically, making use of the numerical code BHAC, we follow the evolution of a geometrically thick matter torus driven into accretion by the development of a magnetorotational instability. By making use of a number of simulations in which we vary the strength of the stellar dipolar magnetic field, we can determine self-consistently the location of the magnetospheric (or Alfv\'en) radius $r_{\rm msph}$ and study how it depends on the magnetic moment $\mu$ and on the accretion rate. Overall, we recover the analytic Newtonian scaling relation, i.e. $r_{\rm msph} \propto B^{4/7}$, but also find that the dependence on the accretion rate is very weak. Furthermore, we find that the material torque correlates linearly with the mass-accretion rate, although both of them exhibit rapid fluctuations. Interestingly, the total torque fluctuates drastically in strong magnetic field simulations and these unsteady torques observed in the simulations could be associated with the spin fluctuations observed in X-ray pulsars.

(Abridged) Context: In the previous paper in this series, we identified that a pentagonal arrangement of five telescopes, using a kernel-nulling beam combiner, shows notable advantages for some important performance metrics for a space-based mid-infrared nulling interferometer over several other considered configurations for the detection of Earth-like exoplanets around solar-type stars. Aims: We aim to produce a practical implementation of a kernel-nulling beam combiner for such a configuration, as well as a discussion of systematic and stochastic errors associated with the instrument. Methods: We develop the mathematical framework around a beam-combiner based on a nulling combiner first suggested by Guyon et al. (2013), and then use it along with the simulator developed in the previous paper to identify the effects of systematic uncertainties. Results: We find that errors in the beam combiner optics, systematic phase errors and the RMS fringe tracking errors result in instrument limited performance at $\sim$4-7 $\mu$m, and zodiacal limited at $\gtrsim$10 $\mu$m. Assuming a beam splitter reflectance error of $|\Delta R| = 5\%$ and phase shift error of $\Delta\phi = 3$ degrees, we find that the fringe tracking RMS should be kept to less than 3 nm in order to be photon limited, and the systematic piston error be less than 0.5 nm to be appropriately sensitive to planets with a contrast of 1$\times 10^{-7}$ over a 4-19 $\mu$m bandpass. We also identify that the beam combiner design, with the inclusion of a well positioned shutter, provides an ability to produce robust kernel observables even if one or two collecting telescopes were to fail. The resulting four telescope combiner, when put into an X-array formation, results in a transmission map with an relative SNR equivalent to 80% of the fully functioning X-array combiner.

We report on the detection of transit timing variations (TTV) of WASP-161b by using the combination of TESS data and archival data. The midpoint of the transits in TESS data are offset by $\sim$67 minutes in Jan. 2019, and $\sim$203 minutes in Jan. 2021, based on the ephemeris published in previous work. We are able to reproduce the transit timings from the archival light curve (SSO-Europa; Jan. 2018) and find SSO-Europa timing is consistent with the published ephemeris under a constant period assumption. Conversely, we find that the SSO-Europa transit midpoint indicates a 6.62-minute variation at 4.40 $\sigma$ compared to the prediction obtained from TESS timings, and a constant orbit period assumption. The TTVs could be modeled with a quadratic function, yielding a constant period change. The period derivative $\dot{P}$ is -1.16$\times$10$^{-7}\pm$2.25$\times$10$^{-8}$ days per day (or $-3.65$ s/year), derived from SSO-Europa and TESS timing. Different scenarios, including a decaying period and apsidal precession can potentially explain these TTVs but they both introduce certain inconsistencies. We have obtained CHEOPS observations for two transits in Jan. 2022 to distinguish between different TTV scenarios. We expect the timing to vary by 5 minutes, compared to the timing predicted from SSO-Europa and TESS with a constant period assumption.

We conduct 24.4~fps optical observations of repeating Fast Radio Burst (FRB) 20190520B using Tomo-e Gozen, a high-speed CMOS camera mounted on the Kiso 105-cm Schmidt telescope, simultaneously with radio observations carried out using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). We succeeded in the simultaneous optical observations of 11 radio bursts that FAST detected. However, no corresponding optical emission was found. The optical fluence limits as deep as 0.068 Jy ms are obtained for the individual bursts (0.029 Jy ms on the stacked data) corrected for the dust extinction in the Milky Way. The fluence limit is deeper than those obtained in the previous simultaneous observations for an optical emission with a duration $\gtrsim 0.1$ ms. Although the current limits on radio--optical spectral energy distribution (SED) of FRBs are not constraining, we show that SED models based on observed SEDs of radio variable objects such as optically detected pulsars, and a part of parameter spaces of theoretical models in which FRB optical emission is produced by inverse-Compton scattering in a pulsar magnetosphere or a strike of a magnetar blastwave into a hot wind bubble, can be ruled out once a similar fluence limit as in our observation is obtained for a bright FRB with a radio fluence $\gtrsim 5$ Jy ms.

We review dark energy models which can present non-negligible fluctuations on scales smaller than Hubble radius. Both linear and nonlinear evolutions of dark energy fluctuations are discussed. The linear evolution has a well-established framework, based on linear perturbation theory in General Relativity, and is well studied and implemented in numerical codes. We highlight the main results from linear theory to explain how dark energy perturbations become important on the scales of interest for structure formation. Next, we review some attempts to understand the impact of clustering dark energy models in the nonlinear regime, usually based on generalizations of the Spherical Collapse Model. We critically discuss the proposed generalizations of the Spherical Collapse Model that can treat clustering dark energy models and their shortcomings. Proposed implementations of clustering dark energy models in halo mass functions are reviewed. We also discuss some recent numerical simulations capable of treating dark energy fluctuations. Finally, we make an overview of the observational predictions based on these models.

The wind collision region (WCR) in a colliding wind binary (CWB) is a particularly violent place, as such, it is surprising that it is also a region where significant quantities of interstellar dust can form. In extreme cases, approximately 30% of the total mass loss rate of a system can be converted into dust. These regions are poorly understood, as observation and simulation of these systems are difficult. In our previous paper we simulated dust growth in CWB systems using an advected scalar model and found our model to be suitable for qualitative study. For this paper we simulated the periodic dust forming CWB (WCd) system WR140 with our dust model, to determine how dust growth changes over the systems periastron passage. We found that dust production increases significantly at periastron passage, which is consistent with IR emission of the surrounding dusty shell. We also find that the dust production rate of the system decreases rapidly as the stars recede from each other, though the rate of decrease is significantly lower than the rate of increase during periastron passage. This was found to be due to strong cooling and its associated thermal instabilities, resulting in cool, high-density pockets of gas in the WCR where dust forms. The WCR also shows a degree of hysteresis, resulting in a radiative post-shock flow even when the stars are separated enough for the region to behave adiabatically.

Applying the suggestion that vigorous core convection during core helium flash on the tip of the red giant branch (RGB) of low mass stars excites waves that carry energy to the envelope and inflate it for few years, we find that in binary systems this process leads to a sharp peak of extreme horizontal branch (EHB; sdB and sdO) stars with masses of ~0.47Mo and ease the formation of elliptical planetary nebulae (PNe) by such low mass stars. Using the open-source mesa-binary we follow the evolution of a number of eccentric binary systems with an initial primary stellar mass of 1.6Mo. The energy that the waves carry to the envelope leads to envelope expansion at the tip of the RGB. The inflated RGB star engulfs many secondary stars to start a common envelope evolution (CEE) that otherwise would not occur. If the secondary star manages to remove most of the RGB envelope the primary evolves to become an EHB star with a mass of ~0.47Mo. In cases where the secondary star does not have time to spiral-in to close orbits, it ends at a large orbit and leaves a massive enough envelope for the primary star to later evolve along the asymptotic giant branch and to engulf the secondary star, therefore forming a non-spherical planetary nebula.

We present an assessment of the greenhouse gases emissions of the Institute for Research in Astrophysics and Planetology (IRAP), located in Toulouse (France). It was performed following the established "Bilan Carbone" methodology, over a large scope compared to similar previous studies, including in particular the contribution from the purchase of goods and services as well as IRAP's use of external research infrastructures, such as ground-based observatories and space-borne facilities. The carbon footprint of the institute for the reference year 2019 is 7400 +/- 900 tCO2e. If we exclude the contribution from external research infrastructures to focus on a restricted perimeter over which the institute has some operational control, IRAP's emissions in 2019 amounted to 3300 +/- 400 tCO2e. Over the restricted perimeter, the contribution from purchasing goods and services is dominant, about 40% of the total, slightly exceeding the contribution from professional travel including hotel stays, which accounts for 38%. Local infrastructures make a smaller contribution to IRAP's carbon footprint, about 25% over the restricted perimeter. We note that this repartition may be specific to IRAP, since the energy used to produce the electricity and heating has a relatively low carbon footprint. Over the full perimeter, the large share from the use of ground-based observatories and space-borne facilities and the fact that the majority of IRAP purchases are related to instrument development indicate that research infrastructures represent the most significant challenge for reducing the carbon footprint of research at our institute. With ~260 staff members employed, our results imply that performing research in astronomy and astrophysics at IRAP according to the standards of 2019 produces average GHG emissions of 28 tCO2e/yr per person involved in that activity (Abridged).

We analyze the distribution of known multi-planet systems ($N \geq 3$) in the plane of mean-motion ratios, and compare it with the resonance web generated by two-planet mean-motion resonances (2P-MMR) and pure 3-planet commensurabilities (pure 3P-MMR). We find intriguing evidence of a statistically significant correlation between the observed distribution of compact low-mass systems and the resonance structure, indicating a possible causal relation. While resonance chains such as Kepler-60, Kepler-80 and TRAPPIST-1 are strong contributors, most of the correlation appears to be caused by systems not identified as resonance chains. Finally, we discuss their possible origin through planetary migration during the last stages of the primordial disc and/or an eccentricity damping process.

The high-resolution transmission molecular spectroscopic database (HITRAN) has recently introduced line-by-line broadening parameters (and their temperature dependence) appropriate for the dominant constituents of planetary atmospheres. The latest version of the database, HITRAN2020, has recently been released and constitutes a large and expansive update. In this work, line shape codes suitable for calculating microwave spectra have been implemented within the HITRAN Application Programming Interface (HAPI). These new additions allow for spectroscopic calculations of microwave absorbing species pertinent to current and future studies of the atmospheres of Jupiter and Venus, and more generally for the atmospheres of gas giants and rocky planets. The inversion spectrum of the NH$_3$ molecule broadened by H$_2$, He and H$_2$O dominates the microwave region of Jupiter. Whereas for Venus, accurate spectroscopic data of SO$_2$ broadened by CO$_2$ is necessary in order to determine its significance, if any, on the reported detection of PH$_3$ in the Venusian upper cloud deck. Comparisons have been made to available microwave laboratory opacities and the following results illustrate that HITRAN data can be used in conjunction with HAPI to reproduce the existing experimental measurements and provide reliable calculation of planetary opacities.

We obtained extensive narrowband photoelectric photometry of Comet 21P/Giacobini-Zinner with observations spanning 33 years. The original data from 1985 (Schleicher et al. 1987) were re-reduced and are presented along with data from three additional apparitions including 2018/19. The original conclusion regarding Giacobini-Zinner's chemical composition remains unchanged, with it having a 4-6x depletion in the carbon-chain molecules C2 and C3, and in NH, as compared with both OH and CN. The comet continues to exhibit a large asymmetry in production rates as a function of time and heliocentric distance, with production reaching a peak 3-5 weeks prior to perihelion. All species, including dust, follow the same general production rate curve each apparition, and the carbon-bearing species are always very similar to one another. However, OH and NH each differ in detail from the carbon-bearing species, implying somewhat varied composition between source regions. Longer term, there are only small secular changes among the apparitions before and near perihelion, but larger changes are evident as the comet recedes from the Sun, suggestive of a progressive precession of the rotation axis.

Large-scale volcanism has played a critical role in the long-term habitability of Earth. Contrary to widely held belief, volcanism, rather than impactors, has had the greatest influence on and bears most of the responsibility for large-scale mass extinction events throughout Earth's history. We examine the timing of large igneous provinces (LIPs) throughout Earth's history to estimate the likelihood of nearly simultaneous events that could drive a planet into an extreme moist or runaway greenhouse, leading to the end of volatile cycling and causing the heat death of formerly temperate terrestrial worlds. In one approach, we make a conservative estimate of the rate at which sets of near-simultaneous LIPs (pairs, triplets, and quartets) occur in a random history statistically the same as Earth's. We find that LIPs closer in time than 0.1-1 million yr are likely; significantly, this is less than the time over which terrestrial LIP environmental effects are known to persist. In another approach, we assess the cumulative effects with simulated time series consisting of randomly occurring LIP events with realistic time profiles. Both approaches support the conjecture that environmental impacts of LIPs, while narrowly avoiding grave effects on the climate history of Earth, could have been responsible for the heat death of our sister world Venus.

Experimental refinements and technical innovations in the field of extensive air shower telescopes have enabled measurements of Galactic cosmic-ray interactions in the sub-PeV (100 TeV to 1 PeV) range, providing new avenues for the search for new physics and dark matter. For the first time, we exploit sub-PeV (from 10 TeV to 1 PeV) observations of Galactic diffuse gamma rays by Tibet AS$\gamma$ and HAWC to search for an axion-like-particle (ALP) induced gamma-ray signal directly linked to the origin of the IceCube extragalactic high-energy neutrino flux. Indeed, the production of high-energy neutrinos in extragalatic sources implies the concomitant production of gamma rays at comparable energies. Within the magnetic field of the neutrino emitting sources, gamma rays may efficiently convert into ALPs and escape their host galaxy unattenuated, propagate through intergalactic space, and reconvert into gamma rays in the magnetic field of the Milky Way. Such a scenario creates an all-sky diffuse high-energy gamma-ray signal in the sub-PeV range. Accounting for the guaranteed Galactic astrophysical gamma-ray contributions from cosmic-ray interactions with gas and radiation and from sub-threshold sources, we set competitive upper limits on the photon-ALP coupling constant $g_{a\gamma\gamma}$. We find $g_{a\gamma\gamma} < 2.1\times10^{-11}$ GeV$^{-1}$ for ALP masses $m_a \leq 2\times10^{-7}$ eV at a 95$\%$ confidence level, progressively closing the mass gap towards ADMX limits.

We calculate the quasinormal modes (QNM) frequencies of a test massless scalar field around static black holes in f(T) gravity. Focusing on quadratic f(T) modifications, which is a good approximation for every realistic f(T) theory, we first extract the spherically symmetric solutions using the perturbative method, imposing two ansatze for the metric functions, which suitably quantify the deviation from the Schwarzschild solution. Moreover, we extract the effective potential and then we calculate the QNM frequency of the obtained solutions. Firstly, we solve the Schrodinger-like equation numerically using the discretization method, and we extract the frequency and the time evolution of the dominant mode applying the function fit method. Secondly, we perform a semi-analytical calculation by applying the third-order WKB approximation method. We show that the results for f(T) gravity are different comparing to General Relativity, and in particular we obtain a different slope and period of the field decay behavior for different values of the model parameter. Hence, under the light of gravitational-wave observations of increasing accuracy from binary systems, the whole analysis could be used as an additional tool in order to test General Relativity and examine whether torsional gravitational modifications are possible.

False vacuum decay, a quantum mechanical first-order phase transition in scalar field theories, is an important phenomenon in early universe cosmology. Recently, a new real-time semi-classical technique based on ensembles of lattice simulations was introduced to describe false vacuum decay. In this context, or any other lattice simulation, the effective potential experienced by long-wavelength modes is not the same as the bare potential. To make quantitative predictions using the real-time semi-classical techniques, it is therefore necessary to understand the redefinition of model parameters and the corresponding deformation of the vacuum state, as well as stochastic contributions that require modeling of unresolved subgrid modes. In this work, we focus on the former corrections and compute the expected modification of the true and false vacuum effective mass, which manifests as a modified dispersion relationship for linear fluctuations about the vacuum. We compare these theoretical predictions to numerical simulations and find excellent agreement. Motivated by this, we use the effective masses to fix the shape of a parameterized effective potential, and explore the modeling uncertainty associated with non-linear corrections. We compute the decay rates in both the Euclidean and real-time formalisms, finding qualitative agreement in the dependence on the UV cutoff. These calculations further demonstrate that a quantitative understanding of the rates requires additional corrections.

Large-field inflation is a major class of inflation models featuring a near- or super-Planckian excursion of the inflaton field. We point out that the large excursion generically introduces significant scale dependence to spectator fields through inflaton couplings, which in turn induces characteristic distortions to the oscillatory shape dependence in the primordial bispectrum mediated by a spectator field. This so-called cosmological collider signal can thus be a useful indicator of large field excursions. We show an explicit example with signals from the "tower states" motivated by the swampland distance conjecture.

We demonstrate that ab-initio calculations in QCD at high densities offer significant and nontrivial information about the equation of state of matter in the cores of neutron stars, going beyond that which is obtainable from current astrophysical observations. We do so by extrapolating the equation of state to neutron-star densities using a Gaussian process and conditioning it sequentially with astrophysical observations and QCD input. Using our recent work, imposing the latter does not require an extrapolation to asymptotically high density. We find the QCD input to be complementary to the astrophysical observations, offering strong additional constraints at the highest densities reached in the cores of neutron stars; with the QCD input, the equation of state is no longer prior dominated at any density. The QCD input reduces the pressure and speed of sound at high densities, and it predicts that binary collisions of equal-mass neutron stars will produce a black hole with greater than $95\%$ ($68\%$) credence for masses $M \geq 1.38 M_\odot$ ($M \geq 1.25 M_\odot$). We provide a Python implementation of the QCD likelihood function so that it can be conveniently used within other inference setups.

Neutrinos might interact among themselves through forces that have so far remained hidden. Throughout the history of the Universe, such secret interactions could lead to scatterings between the neutrinos from supernova explosions and the non-relativistic relic neutrinos left over from the Big Bang. Such scatterings can boost the cosmic neutrino background (C$\nu$B) to energies of ${\cal O}$(MeV), making it, in principle, observable in experiments searching for the diffuse supernova neutrino background. Assuming a model-independent four-Fermi interaction, we determine the upscattered cosmic neutrino flux, and derive constraints on such secret interactions from the latest results from Super-Kamiokande. Furthermore, we also study prospects for detection of the boosted flux in future lead-based coherent elastic neutrino-nucleus scattering experiments.

We study the gravitational production of dark photon dark matter during inflation, when dark photons acquire mass by the Higgs mechanism. In the previous study, it was assumed that the dark photon has a St\"uckelberg mass, or a mass generated by the Higgs mechanism with a sufficiently heavy Higgs boson. In this paper we consider a case in which the Higgs boson is not fully decoupled; the Higgs field changes its vacuum expectation value after inflation. Then, the dark photon mass also changes with time after inflation, and the time evolution of the longitudinal mode is different from the case with a St\"{u}ckelberg mass. Consequently, the spectrum of the dark photon energy density can have two peaks at an intermediate scale and a small scale. We show that the dark photon can explain the dark matter if its current mass is larger than $6 \, \mu {\rm eV} \times (H_I / 10^{14} \, {\rm GeV})^{-4}$ and smaller than $0.8 \, {\rm GeV} \times (H_I / 10^{14} \, {\rm GeV})^{-3/2}$, with $H_I$ being the Hubble parameter during inflation. A higher mass is required if one considers a larger gauge coupling constant. The result for the St\"uckelberg mass can be reproduced in the limit of a small gauge coupling constant. We also comment on the constraints set by various conjectures in quantum gravity theory.

The future space-based gravitational-wave detector LISA will deliver rich and information-dense data by listening to the milliHertz Universe. The measured time series will contain the imprint of tens of thousands of detectable Galactic binaries constantly emitting, tens of supermassive black hole merger events per year, tens of stellar-origin black holes, and possibly thousands of extreme mass-ratio inspirals. On top of that, we expect to detect the presence of stochastic gravitational wave backgrounds and bursts. Finding and characterizing many such sources is a vast and unsolved task. The LISA Data Challenges (LDCs) are an open and collaborative effort to tackle this exciting problem. A new simulated data set, nicknamed Sangria, has just been released with the purpose of tackling mild source confusion with idealized instrumental noise. This presentation will describe the LDC strategy, showcase the available datasets and analysis tools, and discuss future efforts to prepare LISA data analysis.

The search for correlations between secondary cosmic ray detection rates and seismic effects has long been a subject of investigation motivated by the hope of identifying a new precursor type that could feed a global early warning system against earthquakes. Here we show for the first time that the average variation of the cosmic ray detection rates correlates with the global seismic activity to be observed with a time lag of approximately two weeks, and that the significance of the effect varies with a periodicity resembling the undecenal solar cycle, with a shift in phase of around three years, exceeding 6 sigma at local maxima. The precursor characteristics of the observed correlations point to a pioneer perspective of an early warning system against earthquakes.

We compute the equation of state for an ensemble of degenerate fermions by using the curved spacetime of a slowly rotating axially symmetric star. We show that the equation of state computed in such curved spacetime depends on the gravitational time dilation as well as on the dragging of inertial frames, unlike an equation of state computed in a globally flat spacetime. The effect of gravitational time dilation leads to a significant enhancement of the maximum mass limit of a degenerate neutron star. However, such an enhancement due to the frame-dragging effect is extremely small.

We show that a primeval seed magnetic field arises as a natural consequence of spin-degeneracy breaking of fermions caused by the dragging of inertial frames in the curved spacetime of spinning astrophysical bodies. This seed magnetism would arise even due to electrically neutral fermions such as neutrons. As an example, we show that an ideal neutron star spinning at $500$ revolutions per second, having mass $0.83$ M$_{\odot}$ and described by an ensemble of degenerate neutrons would have $0.12$ Gauss seed magnetic field at its center arising through the breaking of spin-degeneracy.

We consider decay of a particle with some energy $E_{0}>0$ inside the ergosphere of a black hole. After the first decay one of particles with the energy $E_{1}<0$ falls towards a black hole while the second one with $% E_{2}>E_{0}\,\ $moves in the outward direction. It bounces back from a reflecting shell and, afterwards, the process repeats. For radial motion of charged particles in the Reissner-Nordst\"{o}m metric, the result depends strongly on a concrete scenario. In particular, an indefinitely large growth of energy inside a shell is possible that gives rise to a black-hole bomb. We also consider a similar multiple process with neutral particles in the background of a rotating axially symmetric stationary black hole. We demonstrate that, if decay occurs in the turning point, a black-hole bomb in this case is impossible at all. For a generic point inside the ergoregion, there is condition for a black-hole bomb to exist. It relates the ratio of masses before and after decay and the velocity of a fragment in the center of mass frame.