Papers for Wednesday, Mar 31 2021

Variable active galactic nuclei showing periodic light curves have been proposed as massive black hole binary (MBHB) candidates. In such scenarios the periodicity can be due to relativistic Doppler-boosting of the emitted light. This hypothesis can be tested through the timing of scattered polarized light. Following the results of polarization studies in type I nuclei and of dynamical studies of MBHBs with circumbinary discs, we assume a coplanar equatorial scattering ring, whose elements contribute differently to the total polarized flux, due to different scattering angles, levels of Doppler boost, and line-of-sight time delays. We find that in the presence of a MBHB, both the degree of polarization and the polarization angle have periodic modulations. The minimum of the polarization degree approximately coincides with the peak of the light curve, regardless of the scattering ring size. The polarization angle oscillates around the semi-minor axis of the projected MBHB orbital ellipse, with a frequency equal either to the binary's orbital frequency (for large scattering screen radii), or twice this value (for smaller scattering structures). These distinctive features can be used to probe the nature of periodic MBHB candidates and to compile catalogs of the most promising sub-pc MBHBs. The identification of such polarization features in gravitational-wave detected MBHBs would enormously increase the amount of physical information about the sources, allowing the measurement of the individual masses of the binary components, and the orientation of the line of nodes on the sky, even for monochromatic gravitational wave signals.

We present full-sky maps of the Integrated Sachs-Wolfe effect (ISW) for the MICE Grand Challenge lightcone simulation up to redshift 1.4. The maps are constructed in the linear regime using spherical Bessel transforms. We compare and contrast this procedure against analytical approximations found in the literature. By computing the ISW in the linear regime we remove the substantial computing and storage resources required to calculate the non-linear Rees-Sciama effect. Since the linear ISW is $~10^{2}\times$ larger in $\Delta T/T$ this has a neglegible impact on the maps produced and only becomes relevant on scales which are dominated by cosmic microwave background (CMB) anisotropies. The MICE simulation products have been extensively used for studies involving current and future galaxy surveys. The availability of these maps will allow MICE to be used for future galaxy and CMB cross-correlation studies, ISW reconstruction studies and ISW void-stacking studies probed by galaxy surveys such as DES, DESI, Euclid and Rubin LSST. The pipeline developed in this study is provided as a public Python package pyGenISW. This could be used in future studies for constructing the ISW from existing and future simulation suites probing vast sets of cosmological parameters and models.

We study radiation pressure due to Lyman alpha line photons, obtaining and exploring analytical expressions for the force-multiplier, $M_F(N_H, Z) = F_\alpha/(L_\alpha/c)$, as a function of gas column density, $N_H$, and metallicity, $Z$, for both dust-free and dusty media, employing a WKB approach for the latter case. Solutions for frequency offset emission to emulate non-static media moving with a bulk velocity $v$, have also been obtained. We find that, in static media, Ly$\alpha$ pressure dominates over both photoionization and dust-mediated UV radiation pressure in a very wide parameter range ($16 < \log N_H < 23$; $-4 < \log[Z/Z_\odot] < 0$). For example, it overwhelms the other two forces by 10 (300) times in standard (low-$Z$) star-forming clouds. Thus, in agreement with previous studies, we conclude that Ly$\alpha$ pressure plays a dominant role in the initial acceleration of the gas around luminous sources, and must be implemented in galaxy formation, evolution and outflow models and simulations.

We present DES14X2fna, a high-luminosity, fast-declining type IIb supernova (SN IIb) at redshift $z=0.0453$, detected by the Dark Energy Survey (DES). DES14X2fna is an unusual member of its class, with a light curve showing a broad, luminous peak reaching $M_r\simeq-19.3$ mag 20 days after explosion. This object does not show a linear decline tail in the light curve until $\simeq$60 days after explosion, after which it declines very rapidly (4.38$\pm$0.10 mag 100 d$^{-1}$ in $r$-band). By fitting semi-analytic models to the photometry of DES14X2fna, we find that its light curve cannot be explained by a standard $^{56}$Ni decay model as this is unable to fit the peak and fast tail decline observed. Inclusion of either interaction with surrounding circumstellar material or a rapidly-rotating neutron star (magnetar) significantly increases the quality of the model fit. We also investigate the possibility for an object similar to DES14X2fna to act as a contaminant in photometric samples of SNe Ia for cosmology, finding that a similar simulated object is misclassified by a recurrent neural network (RNN)-based photometric classifier as a SN Ia in $\sim$1.1-2.4 per cent of cases in DES, depending on the probability threshold used for a positive classification.

The dynamical ellipticity of a planet expresses the departure of its mass distribution from spherical symmetry. It enters as a parameter in the description of a planet's precession and nutation, as well as other rotational normal modes. In the case of the Earth, uncertainties in this quantity's history produce an uncertainty in the solutions for the past evolution of the Earth-Moon system. Constraining this history has been a target of interdisciplinary efforts as it represents an astro-geodetic parameter whose variation is shaped by geophysical processes, and whose imprints can be found in the geological signal. We revisit the classical problem of its variation during ice ages, where glacial cycles exerted a varying surface loading that had altered the shape of the geoid. In the framework of glacial isostatic adjustment, and with the help of a recent paleoclimatic proxy of ice volume, we present the evolution of the dynamical ellipticity over the Cenozoic ice ages. We map out the problem in full generality identifying major sensitivities to surface loading and internal variations in parameter space. This constrained evolution is aimed to be used in future astronomical computations of the orbital and insolation quantities of the Earth.

We present the results of a deep K_s-band (2.1 um) imaging survey of the Spitzer/HETDEX Exploratory Large Area (SHELA) field using the NEWFIRM near-infrared (NIR) camera on the KPNO Mayall 4-m telescope. This NEWFIRM HETDEX Survey (NHS) reaches a 5-sigma depth of 22.4 AB mag (2"-diameter apertures corrected to total), is ~50% and 90% complete at K~22.65 and K~22.15, respectively, and covers 22 deg^2 of the 24 deg^2 SHELA Spitzer/ IRAC footprint (within ``Stripe 82' of the Sloan Digital Sky Survey). We present a K_s-band-selected catalog which includes deep ugriz imaging from the Dark Energy Camera and 3.6 and 4.5 um imaging from Spitzer/IRAC, with forced-photometry of 1.7 million sources across 17.5 deg^2. The large area and moderate depth of this catalog enables the study of the most massive galaxies at high redshift, and minimizes uncertainties associated with counting statistics and cosmic variance. As a demonstration, we derive stellar masses (M*) and star-formation rates (SFRs) for candidate galaxies at 3 < z < 5, and select a conservative sample of nine candidate massive (M* > 10^11 Msol) quiescent galaxies, which have measured SFRs significantly below the main-sequence at this redshift. Five are ultra-massive with M* > 10^12, though uncertainties in IRAC blending, gravitational lensing, or AGN emission could result in true masses which are lower. Simulations predict that these galaxies should be extremely rare, thus we discuss what physical processes in models could be altered to allow the formation of such massive quiescent galaxies at such early times.

Context: Precise chemical abundances coupled with reliable ages are key ingredients to understand the chemical history of our Galaxy. Open Clusters (OCs) are useful for this purpose because they provide ages with good precision. Aims: The aim of this work is to investigate the relations of different chemical abundance ratios vs age traced by red clump (RC) stars in OCs. Methods: We analyze a large sample of 209 reliable members in 47 OCs with available high-resolution spectroscopy. We applied a differential line-by-line analysis to provide a comprehensive chemical study of 25 chemical species. This sample is among the largest samples of OCs homogeneously characterized in terms of atmospheric parameters, detailed chemistry, and ages. Results: In our metallicity range (-0.2<[M/H]<+0.2) we find that while most Fe-peak and alpha elements have flat dependence with age, the s-process elements show decreasing trends with increasing age with a remarkable knee at 1 Gyr. For Ba, Ce, Y, Mo and Zr we find a plateau at young ages (< 1 Gyr). We investigate the relations of all possible combinations among the computed chemical species with age. We find 19 combinations with significant slopes, including [Y/Mg] and [Y/Al]. The ratio [Ba/alpha] is the one with the most significant correlations found. Conclusions: We find that the [Y/Mg] relation found in the literature using Solar twins is compatible with the one found here in the Solar neighbourhood. The age-abundance relations show larger scatter for clusters at large distances (d>1 kpc) than for the Solar neighbourhood, particularly in the outer disk. We conclude that these relations need to be understood also in terms of the complexity of the chemical space introduced by the Galactic dynamics, on top of pure nucleosynthetic arguments, especially out of the local bubble.

A promising channel for producing binary black hole mergers is the Lidov-Kozai orbital resonance in hierarchical triple systems. While this mechanism has been studied in isolation, the distribution of such mergers in time and across star-forming environments is not well characterized. In this work, we explore Lidov-Kozai-induced black hole mergers in open clusters, combining semi-analytic and Monte Carlo methods to calculate merger rates and delay times for eight different population models. We predict a merger rate density of $\sim$1--10\,Gpc$^{-3}$\,yr$^{-1}$ for the Lidov-Kozai channel in the local universe, and all models yield delay-time distributions in which a significant fraction of binary black hole mergers (e.g., $\sim$20\%--50\% in our baseline model) occur during the open cluster phase. Our findings suggest that a substantial fraction of mergers from hierarchical triples occur within star-forming regions in spiral galaxies.

We present optical spectroscopy for 18 halo white dwarfs identified using photometry from the Canada-France Imaging Survey and Pan-STARRS1 DR1 3$\pi$ survey combined with astrometry from Gaia DR2. The sample contains 13 DA, 1 DZ, 2 DC, and two potentially exotic types of white dwarf. We fit both the spectrum and the spectral energy distribution in order to obtain the temperature and surface gravity, which we then convert into a mass, and then an age, using stellar isochrones and the initial-to-final mass relation. We find a large spread in ages that is not consistent with expected formation scenarios for the Galactic halo. We find a mean age of 9.03$^{+2.13}_{-2.03}$ Gyr and a dispersion of 4.21$^{+2.33}_{-1.58}$ Gyr for the inner halo using a maximum likelihood method. This result suggests an extended star formation history within the local halo population.

Based on the assumption that the dark energy possessing bulk viscosity is homogenously and isotropically permeated in the universe, we propose three new viscous dark energy (VDE) models to characterize the accelerating universe. By constraining these three models with the latest cosmological observations, we find that they just deviate very slightly from the standard cosmological model and can alleviate effectively the current $H_0$ tension between the local observation by the Hubble Space Telescope and the global measurement by the Planck Satellite. Interestingly, we conclude that a spatially flat universe in our VDE model with cosmic curvature is still supported by current data, and the scale invariant primordial power spectrum is strongly excluded at least at the $5.5\sigma$ confidence level in three VDE models as the Planck result. We also give the $95\%$ upper limits of the typical bulk viscosity parameter $\eta$ in three VDE scenarios.

We perform the first comprehensive constraints on two Finslerian models, i.e., the simplest Finslerian $\Lambda$CDM (F$\Lambda$) and non-flat F$\Lambda$ models using the Type Ia supernovae, baryonic acoustic oscillations, cosmic microwave background and lensing observations. Using the most stringent constraints we can provide, we find that the constrained typical parameters of both Finslerian models are all consistent with zero at the $2\sigma$ confidence level (C.L.). This means that both Finslerian models just deviate very slightly from the $\Lambda$CDM one, which is also verified by using two geometrical diagnostics to distinguish three different models from each other. We also find that a spatially flat universe is still preferred in the framework of Finsler geometry, that our measured values of the spectral index $n_s$ of primordial power spectrum in both Fisnlerian models exclude the scale invariance at more than $8\sigma$ C.L., and that the current $H_0$ tension can be relieved from $3.4\sigma$ to $2.9\sigma$ and $2.8\sigma$ in the F$\Lambda$ and non-flat F$\Lambda$ models, respectively.

We report on the Swift/XRT Deep Galactic Plane Survey discovery and multi-wavelength follow-up observations of a new intermediate polar Cataclysmic Variable, Swift J183920.1-045350. A 449.7 s spin period is found in Xmm-Newton and NuSTAR data, accompanied by a 459.9 s optical period that is most likely the synodic, or beat period, produced from a 5.6 h orbital period. The orbital period is seen with moderate significance in independent long-baseline optical photometry observations with ZTF and SAAO. We find that the source X-ray pulsed fraction decreases with increasing energy. The X-ray spectra are consistent with the presence of an Fe emission line complex with both local and interstellar absorption. In the optical spectra, strong H$\alpha{}$, H I, He I and He II emission lines are observed, all common features in magnetic CVs. The source properties are thus typical of known intermediate polars, with the exception of its estimated distance of 2.26$^{+1.93}_{-0.83}$ kpc, which is larger than typical, extending the reach of the CV population in our Galaxy.

Recent studies show that the inner Galactic regions host genuine bulge globular clusters, but also halo intruders, complex remnants of primordial building blocks, and objects likely accreted during major merging events. In this study we focus on the properties of M 28, a very old and massive cluster currently located in the Galactic bulge. We analysed wide-field infrared photometry collected by the VVV survey, VVV proper motions, and intermediate-resolution spectra in the calcium triplet range for 113 targets in the cluster area. Our results in general confirm previous estimates of the cluster properties available in the literature. We find no evidence of differences in metallicity between cluster stars, setting an upper limit of Delta[Fe/H]<0.08 dex to any internal inhomogeneity. We confirm that M 28 is one of the oldest objects in the Galactic bulge (13-14 Gyr). From this result and the literature data, we find evidence of a weak age-metallicity relation among bulge globular clusters that suggests formation and chemical enrichment. In addition, wide-field density maps show that M 28 is tidally stressed and that it is losing mass into the general bulge field. Our study indicates that M 28 is a genuine bulge globular cluster, but its very old age and its mass loss suggest that this cluster could be the remnant of a larger structure, possibly a primeval bulge building block.

We revisit constraints on annihilating dark matter based on the global 21cm signature observed by EDGES. For this purpose, we used the numerical data of the latest N-body simulation performed by state-of-art standard in order to estimate the boost factor at high redshifts ($z$ = 10 - 100), which enhances the annihilation of dark matter in course of structure formation. By taking into account to what fraction injected energy from dark matter annihilation contributes to ionization, excitation and heating of intergalactic medium during dark ages, we estimated how large the global 21cm absorption can be. By assuming the thermal freezeout scenario, we find that $m_{\rm DM} < 15$ GeV and $m_{\rm DM} < 3$ GeV have been excluded at 95$\%$ C.L. for the annihilation modes into $b\bar{b}$ and $e^+ e^-$, respectively.

Relativistic magnetic reconnection is a potential particle acceleration mechanism for high-frequency BL Lacs (HBLs). The {\it Imaging X-ray Polarimetry Explorer} ({\it IXPE}) scheduled to launch in 2021 has the capability to probe the magnetic field evolution in HBLs, examining the magnetic reconnection scenario for the HBL flares. In this paper, we make the first attempt to self-consistently predict HBL X-ray polarization signatures arising from relativistic magnetic reconnection via combined particle-in-cell (PIC) and polarized radiation transfer simulations. We find that although the intrinsic optical and X-ray polarization degrees are similar on average, the X-ray polarization is much more variable in both polarization degree and angle (PD and PA). Given the sensitivity of the {\it IXPE}, it may obtain one to a few polarization data points for one flaring event of nearby bright HBLs Mrk~421 and 501. However, it may not fully resolve the highly variable X-ray polarization. Due to the temporal depolarization, where the integration of photons with variable polarization states over a finite period of time can lower the detected PD, the measured X-ray PD can be considerably lower than the optical counterpart or even undetectable. The lower X-ray PD than the optical thus can be a characteristic signature of relativistic magnetic reconnection. For very bright flares where the X-ray polarization is well resolved, relativistic magnetic reconnection predicts smooth X-ray PA swings, which originate from large plasmoid mergers in the reconnection region.

Spatially resolved spectroscopy from SDSS-IV MaNGA survey has revealed a class of quiescent, relatively common early-type galaxies, termed "red geysers", that possibly host large scale active galactic nuclei driven winds. Given their potential importance in maintaining low level of star formation at late times, additional evidence confirming that winds are responsible for the red geyser phenomenon is critical. In this work, we present follow-up observations with the Echellette Spectrograph and Imager (ESI) at the Keck telescope of two red geysers (z$<$0.1) using multiple long slit positions to sample different regions of each galaxy. Our ESI data with a spectral resolution (R) $\sim$ 8000 improves upon MaNGA's resolution by a factor of four, allowing us to resolve the ionized gas velocity profiles along the putative wind cone with an instrumental resolution of $\rm \sigma = 16~km~s^{-1}$. The line profiles of H$\alpha$ and [NII]$\rm \lambda 6584$ show asymmetric shapes that depend systematically on location $-$ extended blue wings on the red-shifted side of the galaxy and red wings on the opposite side. We construct a simple wind model and show that our results are consistent with geometric projections through an outflowing conical wind oriented at an angle towards the line of sight. An alternative hypothesis that assigns the asymmetric pattern to "beam-smearing" of a rotating, ionized gas disk does a poor job matching the line asymmetry profiles. While our study features just two sources, it lends further support to the notion that red geysers are the result of galaxy-scale winds.

Convolution models are powerful tools in many fields of spectral and image analysis owing to their wide applicability, and X-ray astrophysical spectral analysis is no exception. We found that relativistically broadened Fe K${\alpha}$ line profiles obtained through many convolution models both within and without Xspec show deviations from the profiles produced by their non-convolution counterparts. These discrepancies depend on the energy grid considered and on the shape of both the kernel and the underlying spectrum, but can reach as high as 10% of the flux in certain energy bins. We believe that this effect should be taken into consideration, considering how often these models are used to study spectral features of lower relative intensity, and advise great discretion in using them.

Recent observations of the protoplanetary disc surrounding AB Aurigae have revealed the possible presence of two giant planets in the process of forming. The young measured age of $1-4$Myr for this system allows us to place strict time constraints on the formation histories of the observed planets. Hence we may be able to make a crucial distinction between formation through core accretion (CA) or the gravitational instability (GI), as CA formation timescales are typically Myrs whilst formation through GI will occur within the first $\approx10^4-10^5$yrs of disc evolution. We focus our analysis on the $4-13$M$_{\rm Jup}$ planet observed at $R\approx30$AU. We find CA formation timescales for such a massive planet typically exceed the system's age. The planet's high mass and wide orbit may instead be indicative of formation through GI. We use smoothed particle hydrodynamic simulations to determine the system's critical disc mass for fragmentation, finding $M_{\rm d,crit}=0.3$M$_{\odot}$. Viscous evolution models of the disc's mass history indicate that it was likely massive enough to exceed $M_{\rm d,crit}$ in the recent past, thus it is possible that a young AB Aurigae disc may have fragmented to form multiple giant gaseous protoplanets. Calculations of the Jeans mass in an AB Aurigae-like disc find that fragments may initially form with masses $1.6-13.3$M$_{\rm Jup}$, consistent with the planets which have been observed. We therefore propose that the inferred planets in the disc surrounding AB Aurigae may be evidence of planet formation through GI.

The estimation of the main parameters of star clusters is significant in astrophysical studies. The most important aspect of using the Gaia DR2 survey lies in the positions, parallax, and proper motions of cluster stars with homogeneous photometry that make the membership probability determine with high accuracy. In this respect, depending on Gaia DR2 database, an analysis of the open star cluster Melotte 72 is taking place here. It is located at a distance of 2345+/-108 pc with an age of 1.0+/-0.5 Gyr. In studying the radial density profile, the radius is found to be 5.0+/-0.15 arcmin. The reddening, the luminosity and mass functions, the total mass of the cluster, and the galactic geometrical distances (X_Sun, Y_Sun, Z_Sun), and the distance from the galactic center (R_g ) have been estimated as well. Our study has shown a dynamical relaxation behavior of Melotte 72.

We use the Hyper Suprime-Cam Subaru Strategic Program S19A shape catalog to construct weak lensing shear-selected cluster samples. From aperture mass maps covering $\sim 510$~deg$^2$ created using a truncated Gaussian filter, we construct a catalog of 187 shear-selected clusters that correspond to mass map peaks with the signal-to-noise ratio larger than 4.7. Most of the shear-selected clusters have counterparts in optically-selected clusters, from which we estimate the purity of the catalog to be higher than 95\%. The sample can be expanded to 418 shear-selected clusters with the same signal-to-noise ratio cut by optimizing the shape of the filter function and by combining weak lensing mass maps created with several different background galaxy selections. We argue that dilution and obscuration effects of cluster member galaxies can be mitigated by using background source galaxy samples and adopting the filter function with its inner boundary larger than about $2'$. The large samples of shear-selected clusters that are selected without relying on any baryonic tracer are useful for detailed studies of cluster astrophysics and cosmology.

We propose magnetically arrested disks (MADs) in quiescent black-hole (BH) binaries as the origin of the multiwavelength emission, and argue that this class of sources can dominate the cosmic-ray spectrum around the knee. X-ray luminosities of Galactic BH binaries in the quiescent state are far below the Eddington luminosity, and thus, radiatively inefficient accretion flows (RIAFs) are formed in the inner region. Strong thermal and turbulent pressures in RIAFs produce outflows, which can create large-scale poloidal magnetic fields. These fields are carried to the vicinity of the BH by the rapid inflow motion, forming a MAD. Inside the MAD, non-thermal protons and electrons are naturally accelerated by magnetic reconnections or stochastic acceleration by turbulence. Both thermal and non-thermal electrons emit broadband photons via synchrotron emission, which are broadly consistent with the optical and X-ray data of the quiescent BH X-ray binaries. Moreover, protons are accelerated up to PeV energies and diffusively escape from these MADs, which can account for the cosmic-ray intensity around the knee energy.

We derive the non-thermal velocities (NTVs) in the transition region of an active region using the \ion{Si}{4}~1393.78~{\AA} line observed by the Interface Region Imaging Spectrograph (IRIS) and compare them with the line-of-sight photospheric magnetic fields obtained by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). The active region consists of two strong field regions with opposite polarity, separated by a weak field corridor, that widened as the active region evolved. The means of the NTV distributions in strong-field regions (weak field corridors) range between $\sim$18{--}20 (16{--}18)~km~s$^{-1}$, albeit the NTV maps show much larger range. In addition, we identify a narrow lane in the middle of the corridor with significantly reduced NTV. The NTVs do not show a strong center-to-limb variation, albeit somewhat larger values near the disk center. The NTVs are well correlated with redshifts as well as line intensities. The results obtained here and those presented in our companion paper on Doppler shifts suggest two populations of plasma in the active region emitting in \ion{Si}{4}. The first population exists in the strong field regions and extends partway into the weak field corridor between them. We attribute this plasma to spicules heated to $\sim$0.1 MK (often called type II spicules). They have a range of inclinations relative to vertical. The second population exists in the center of the corridor, is relatively faint, and has smaller velocities, likely horizontal. These results provide further insights into the heating of the transition region.

The 11-year solar cycle is the dominant pattern of solar activity reflecting the oscillatory dynamo mechanism in the Sun. Solar cycles were directly observed since 1700, while indirect proxies suggest their existence over a much longer period of time but generally without resolving individual cycles and their continuity. Here we reconstruct individual cycles for the last millennium using recent 14C data and state-of-the-art models. Starting with the 14C production rate determined from the so far most precise measurements of radiocarbon content in tree rings, solar activity is reconstructed in three physics-based steps: (1) Correction of the 14C production rate for the changing geomagnetic field; (2) Computation of the open solar magnetic flux; and (3) Conversion into sunspot numbers outside of grand minima. Solar activity is reconstructed for the period 971-1900 (85 individual cycles). This more than doubles the number of solar cycles known from direct solar observations. We found that lengths and strengths of well-defined cycles outside grand minima are consistent with those obtained from the direct sunspot observations after 1750. The validity of the Waldmeier rule is confirmed at a highly significant level. Solar activity is found to be in a deep grand minimum when the activity is mostly below the sunspot formation threshold, during about 250 years. Therefore, although considerable cyclic variability in 14C is seen even during grand minima, individual solar cycles can hardly be reliably resolved therein. Three potential solar particle events, ca. 994, 1052 and 1279 AD, are shown. A new about 1000-year long solar activity reconstruction, in the form of annual (pseudo) sunspot numbers with full assessment of uncertainties, is presented based on new high-precision 14C measurements and state-of-the-art models, more than doubling the number of individually resolved solar cycles.

We perform a Bayesian analysis of the maximum mass $M_{\rm TOV}$ of neutron stars with a quark core, incorporating the observational data from tidal deformability of the GW170817 binary neutron star merger as detected by LIGO/Virgo and the mass and radius of PSR J0030+0451 as detected by \nicer. The analysis is performed under the assumption that the hadron-quark phase transition is of first order, where the low-density hadronic matter described in a unified manner by the soft QMF or the stiff DD2 equation of state (EOS) transforms into a high-density phase of quark matter modeled by the generic "Constant-sound-speed" (CSS) parameterization. The mass distribution measured for the $2.14 \,{\rm M}_{\odot}$ pulsar, MSP J0740+6620, is used as the lower limit on $M_{\rm TOV}$. We find the most probable values of the hybrid star maximum mass are $M_{\rm TOV}=2.36^{+0.49}_{-0.26}\,{\rm M}_{\odot}$ ($2.39^{+0.47}_{-0.28}\,{\rm M}_{\odot}$) for QMF (DD2), with an absolute upper bound around $2.85\,{\rm M}_{\odot}$, to the $90\%$ posterior credible level. Such results appear robust with respect to the uncertainties in the hadronic EOS. We also discuss astrophysical implications of this result, especially on the post-merger product of GW170817, short gamma-ray bursts, and other likely binary neutron star mergers.

We analyse emission spectra of WASP-12b from a partial phase curve observed over three epochs with the Hubble Space Telescope, covering eclipse, quadrature, and transit, respectively. As the 1.1-day period phase curve was only partially covered over three epochs, traditional methods to extract the planet flux and instrument systematic errors cannot recover the thermal emission away from the secondary eclipse. To analyse this partial phase curve, we introduce a new method, which corrects for the wavelength-independent component of the systematic errors. Our new method removes the achromatic instrument and stellar variability, and uses the measured stellar spectrum in eclipse to then retrieve a relative planetary spectrum in wavelength at each phase. We are able to extract the emission spectrum of an exoplanet at quadrature outside of a phase curve for the first time; we recover the quadrature spectrum of WASP-12b up to an additive constant. The dayside emission spectrum is extracted in a similar manner, and in both cases we are able to estimate the brightness temperature, albeit at a greatly reduced precision. We estimate the brightness temperature from the dayside (Tday=3186+-677 K) and from the quadrature spectrum (Tquad=2124+-417 K) and combine them to constrain the energy budget of the planet. We compare our extracted relative spectra to global circulation models of this planet, which are generally found to be a good match. However, we do see tentative evidence of a steeper spectral slope in the measured dayside spectrum compared to our models. We find that we cannot match this increased slope by increasing optical opacities in our models. We also find that this spectral slope is unlikely to be explained by a non-equilibrium water abundance, as water advected from the nightside is quickly dissociated on the dayside.

In the recent years, sub/mm observations of protoplanetary disks have discovered an incredible diversity of substructures in the dust emission. An important result was the finding that dust grains of mm size are embedded in very thin dusty disks. This implies that the dust mass fraction in the midplane becomes comparable to the gas, increasing the importance of the interaction between the two components there. We address this problem by means of numerical 2.5D simulations in order to study the gas and dust interaction in fully global stratified disks. To this purpose, we employ the recently developed dust grain module in the PLUTO code. Our model focuses on a typical T Tauri disk model, simulating a short patch of the disk at 10 au which includes grains of constant Stokes number of $St=0.01$ and $St=0.1$, corresponding to grains with sizes of 0.9 cm and 0.9 mm, respectively, for the given disk model. By injecting a constant pebble flux at the outer domain, the system reaches a quasi steady state of turbulence and dust concentrations driven by the streaming instability. For our given setup and using resolutions up to 2500 cells per scale height we resolve the streaming instability, leading to local dust clumping and concentrations. Our results show dust density values of around 10-100 times the gas density with a steady state pebble flux between $3.5 \times 10^{-4}$ and $2.5 \times 10^{-3} M_{\rm Earth}/\mathit{year}$ for the models with $\mathit{St}=0.01$ and $\mathit{St}=0.1$. The grain size and pebble flux for model $\mathit{St}=0.01$ compares well with dust evolution models of the first million years of disk evolution. For those grains the scatter opacity dominates the extinction coefficient at mm wavelengths. These types of global dust and gas simulations are a promising tool for studies of the gas and dust evolution at pressure bumps in protoplanetary disks.

The peculiar morphology of Head-Tail (HT) radio galaxies indicates strong interactions between the radio jets and their intra-cluster medium. We systematically search for HT radio galaxies from LOFAR Two-metre Sky Survey first data release (LoTSS DR1) at 144 MHz frequency. We present here a catalogue of fifty new HT radio sources, among them, five are Narrow-Angle Tailed sources (NATs) and forty-five are Wide Angle Tailed sources (WATs). NATs are characterized by tails bent in a narrow V like shape with less than a ninety-degree opening angle. For WAT radio galaxies, the opening angle between jets is more than ninety degrees which exhibit wide C like morphologies. We found that thirty-one out of fifty HT sources are associated with known galaxy clusters. The various physical properties and statistical studies of these HT sources are also presented in this paper.

Dense aperture arrays provide key benefits in modern astrophysical research. They are flexible, employing cheap receivers, while relying on the ever more sophisticated compute back-end to deal with the complexities of signal processing required for their optimal use. Their advantage is that they offer very large fields of view and are readily scalable to any size, all other things being equal. Since they represent "software telescopes", the science cases these arrays can be applied to are quite broad. Here, we describe the calibration and performance of the AARTFAAC-12 instrument, which is composed of the twelve centrally located stations of the LOFAR array. We go into the details of the data acquisition and pre-processing, we describe the newly developed calibration pipeline as well as the noise parameters of the resulting images. We also present the derived radio source counts at 41.7 MHz and 61 MHz.

The study of Head Tail (HT) radio galaxies track the information of associated galaxy clusters. With the help of the VLA FIRST survey at 1.4 GHz, we detected 607 new HT radio sources, among them, 398 are Wide Angle Tail (WAT) and 216 are Narrow-Angle Tail (NAT) sources. NAT sources generally have `V' shaped structure with an opening angle less than ninety degrees and for WAT sources opening angle between the jets is more than ninety degrees. We found that almost 80 per cent of our sources are associated with a known galaxy cluster. We mentioned various useful physical properties of these HT sources. Taking advantage of a large sample of newly discovered HT sources, various statistical studies have been done. The luminosity range of sources presented in the current paper is $10^{39}$ $\leq$ $L_{1.4GHz}$ $\leq$ $10^{43}$ erg sec$^{-1}$. We identified optical counterparts for 193 WAT and 104 NAT sources. The sources are found up to redshift 2.08.

Models of chemical enrichment and inhomogeneity in high-redshift galaxies are challenging to constrain observationally. In this work, we discuss a novel approach to probe chemical inhomogeneities within long Gamma-Ray Burst (GRB) host galaxies, by comparing the absorption metallicity, Z_abs, from the GRB afterglow (which probes the environment along the line of sight) with the emission-line metallicity, Z_emiss, measured via slit spectroscopy. Using the IllustrisTNG simulation, the theoretical relationship between these metallicity metrics is explored for a range of GRB formation models, varying the GRB progenitor metallicity threshold. For galaxies with fixed Z_emiss, the median value of Z_abs depends strongly on the GRB progenitor threshold metallicity, with Z_abs significantly lower than Z_emiss for high metallicity hosts. Conversely, at fixed Z_abs, the median value of Z_emiss depends primarily on the metallicity distribution of galaxies in IllustrisTNG and their chemical inhomogeneities, offering a GRB-model-independent way to constrain these processes observationally. Currently, only one host galaxy has data for both absorption and emission metallicities (GRB121014A). We re-analyse the emission spectrum and compare the inferred metallicity Z_emiss to a recent Bayesian determination of Z_abs, finding $\log(Z_{\rm emiss}/Z_{\odot}) = \log(Z_{\rm abs}/Z_{\odot}) +0.35^{+ 0.14}_{- 0.25}$, within ~2 standard deviations of predictions from the IllustrisTNG simulation. Future observations with the James Webb Space Telescope will be able to measure Z_emiss for 4 other GRB hosts with known Z_abs values, using ~2 hour observations. While small, the sample will provide preliminary constraints on the Z_abs-Z_emiss relation to test chemical enrichment schemes in cosmological simulations.

Binary neutron star (NS) mergers have been expected to synthesize r-process elements and emit radioactively powered radiation, called kilonova. Although r-process nucleosynthesis was confirmed by the observations of GW170817/AT2017gfo, no trace of individual elements has been identified except for strontium. In this paper, we perform systematic calculations of line strength for bound-bound transitions and radiative transfer simulations in NS merger ejecta toward element identification in kilonova spectra. We find that Sr II triplet lines appear in the spectrum of a lanthanide-poor model, which is consistent with the absorption feature observed in GW170817/AT2017gfo. The synthetic spectrum also shows the strong Ca II triplet lines. This is natural because Ca and Sr are co-produced in the material with relatively high electron fraction and their ions have similar atomic structures with only one s-electron in the outermost shell. The line strength, however, highly depends on the abundance distribution and temperature in the ejecta. For our lanthanide-rich model, the spectra show the features of doubly ionized heavy elements, such as Ce, Tb and Th. Our results suggest that the line forming region of GW170817/AT2017gfo was lanthanide-poor. We show that the Sr II and Ca II lines can be used as a probe of physical conditions in NS merger ejecta. Absence of the Ca II line features in GW170817/AT2017gfo implies that the Ca/Sr ratio is < 0.002 in mass fraction, which is consistent with nucleosynthesis for electron fraction >= 0.40 and entropy per nucleon (in units of Boltzmann constant) >= 25.

Photon is the fundamental quantum of electromagnetic fields, whose mass, $m_{\gamma}$, should be strictly zero in Maxwell's theory. But not all theories adopt this hypothesis. If the rest mass of the photon is not zero, there will be an additional time delay between photons of different frequencies after they travel through a fixed distance. By analyzing the time delay, we can measure or constrain the photon mass. Fast radio bursts (FRBs) -- transient radio bursts characterized by millisecond duration and cosmological propagation -- are excellent astrophysical laboratories to constrain $m_{\gamma}$. In this work we use a catalog of 129 FRBs in a Bayesian framework to constrain $m_{\gamma}$. As a result, we obtain a new bound on the photon mass, $m_{\gamma} \leq 3.1\times 10^{-51}\rm\,kg\simeq 1.7 \times 10^{-15}\,eV/c^2$ ($m_{\gamma} \leq 3.9\times 10^{-51}\rm\,kg \simeq 2.2 \times 10^{-15}\,eV/c^2$) at the $68\%$ $(95\%$) confidence level. The result represents the best limit purely from kinematic analysis of light propagation. The bound on the photon mass will be tighter in the near future with increment in the number of FRBs, more accurate measurement of the redshift for FRBs, and refinement in the knowledge about the origin of dispersion measures (DMs).

In the currently accepted paradigm, dark matter is hypothesized as an explanation of the flat rotation curves of galaxies under the assumption of virialized orbits. The use of millisecond pulsar timing as a probe of Galactic dark matter content is explored as a means of relaxing this assumption. A method of inference of the Galactic potential using the frequency derivative $\dot{\nu}$ is produced, and an estimate for a virialized Galactic rotation curve is given through direct observation of acceleration. The data set used includes 210 pulsars with known $\dot{\nu}$ and astrometric properties, a subset of which also have measured $\ddot{\nu}$. In principle, this enables the exploration of kinematic effects, but in practice, $\ddot{\nu}$ values are found to be too imprecise at present to adequately constrain radial velocities of pulsars. Additionally, surface magnetic field strengths are inferred from $\dot{\nu}$ and the magnetic spin-down contribution to $\ddot{\nu}$ is estimated. For several pulsars the radial velocity is known, and the kinematic contribution to $\ddot{\nu}$ is estimated accordingly. The binary orbital periods of PSR J1713+0747 and other binary pulsars are also used to constrain Galactic mass density models.

The multi-color passband CCD light curves of ASAS J124343+1531.7 and LINEAR 2323566 were first obtained by the 0.84-m Ritchey-Chr\'{e}tien telescope with follow up observations by the WIYN 0.90m Cassegrain telescope. The data from the $Transiting\quad Exoplanet\quad Survey\quad Satellite\quad (TESS)$ of ASAS J124343+1531.7 was also applied for subsequent analysis. By analyzing the data through the W-D program, their mass ratios and fill-out factors were determined as 3.758, 1.438 and 31.8$\%$, 14.9$\%$, respectively. ASAS J124343+1531.7 is a W-subtype median contact binary, while LINEAR 2323566 is a W-subtype shallow contact binary, and the asymmetric light curves prove that they both have the O'Connell effect, which is generally explained by magnetic activity. The equivalent widths (EWs) of H$_\alpha$ lines were calculated, which show they certainly have magnetic activity. Moreover, LINEAR 2323566 has a stronger magnetic activity. The analysis of orbital period changes shows that ASAS J124343+1531.7 has a trend of secular period increase, which is generally explained by the mass transfer from the less massive to the more massive star. According to the estimated absolute parameters, their evolutionary states are discussed. The two components of ASAS J124343+1531.7 are both main sequence stars. While for LINEAR 2323566, the more massive star is a main sequence star, the less massive star has evolved out of main sequence and is over-luminous and over-sized.

Prompt extra power-law (PL) spectral component is discovered in some bright gamma-ray bursts (GRBs), which usually dominates the spectral energy distribution below tens of keV or above $\sim$10 MeV. However, its origin is still unclear. In this letter, we present a systematic analysis for 13 \textit{Fermi} short GRBs as of August 2020, with significant keV--MeV and GeV detections at the prompt emission phase. We find that the extra PL component is a ubiquitous spectral feature for short GRBs, showing up in all 13 analyzed GRBs. The PL indices are mostly harder than $-$2.0, which may be well reproduced by considering the electromagnetic cascade induced by ultra-relativistic protons or electrons accelerated in the prompt emission phase. The average flux of these extra PL components positively correlates with that of the main spectral components, which implies they may share the same physical origin.

Aims.Searching for GeV $\gamma$-ray emission from supernova remnant (SNR) G272.2-3.2, and analyzing the features of its GeV $\gamma$-ray emission. Methods. We analyzed features of the GeV $\gamma$-ray emission from the region of SNR G272.2-3.2 by using Fermitools with 12.4 years of observations from the Fermi Large Area Telescope (Fermi-LAT). These features include the $\gamma$-ray spatial distribution, spectral energy distribution (SED), and light curve (LC). Results. A significant $\gamma$-ray new source with approximately 5$\sigma$ significance level is found from the region of SNR G272.2-3.2; its $\gamma$-ray spatial distribution does not exist extended feature; it has a soft spectrum with a spectral index of 2.56$\pm$0.01; no significant variability of its LC is found; its spatial positions in the X-ray and GeV bands well overlap; we suggest that the new $\gamma$-ray source is likely to be a counterpart of SNR G272.2-3.2.

Very young asteroid families may record processes that accompanied their formation in the most pristine way. This makes analysis of this special class particularly interesting. We studied the very young Adelaide family in the inner part of the main belt. This cluster is extremely close to the previously known Datura family in the space of proper orbital elements and their ages overlap. As a result, we investigated the possibility of a causal relationship between the two families. We identified Adelaide family members in the up-to-date catalogue of asteroids. By computing their proper orbital elements we inferred the family structure. Backward orbital integration of selected members allowed us to determine the age of the family. The largest fragment (525) Adelaide, an S-type asteroid about $10$ km in size, is accompanied by 50 sub-kilometre fragments. This family is a typical example of a cratering event. The very tiny extent in the semi-major axis minimises chances that some significant mean motion resonances influence the dynamics of its members, though we recognise that part of the Adelaide family is affected by weak, three-body resonances. Weak chaos is also produced by distant encounters with Mars. Simultaneous convergence of longitude of node for the orbits of six selected members to that of (525) Adelaide constrains the Adelaide family age to $536\pm 12$ kyr (formal solution). While suspiciously overlapping with the age of the Datura family, we find it unlikely that the formation events of the two families are causally linked. In all likelihood, the similarity of their ages is just a coincidence.

If gamma-ray bursts are at cosmological distances, they must be gravitationally lensed occasionally. The detection of lensed images with millisecond-to-second time delays provides evidence for intermediate-mass black holes, a population that has been difficult to observe. Several studies have searched for these delays in gamma-ray burst light curves, which would indicate an intervening gravitational lens. Among the $\sim 10^4$ gamma-ray bursts observed, there have been a handful of claimed lensing detections, but none have been statistically robust. Here we present a Bayesian analysis identifying gravitational lensing in the light curve of GRB950830. The inferred lens mass depends on the unknown lens redshift $z_l$, and is given by $(1+z_l)M_l = 5.5^{+1.7}_{-0.9}\times 10^4 $ M$_\odot$ (90% credibility), which we interpret as evidence for an intermediate-mass black hole. The most probable configuration, with a lens redshift $z_l\sim 1$ and a gamma-ray burst redshift $z_s\sim 2$, yields a present day number density of $n_\text{imbh}\approx 2.3^{+4.9}_{-1.6}\times10^{3} \text{ Mpc}^{-3}$ (90% credibility) with a dimensionless energy density $\Omega_\text{imbh} \approx 4.6^{+9.8}_{-3.3}\times10^{-4}$. The false alarm probability for this detection is $\sim0.6\%$ with trial factors. While it is possible that GRB950830 was lensed by a globular cluster, it is unlikely since we infer a cosmic density inconsistent with predictions for globular clusters $\Omega_\text{gc} \approx 8 \times 10^{-6}$ at 99.8% credibility. If a significant intermediate-mass black hole population exists, it could provide the seeds for the growth of supermassive black holes in the early Universe.

In recent years, survey facilities have started to provide remarkable new opportunities for transient astronomy. This is made possible by combining unprecedented fields of view with high sensitivity in uninvestigated wavelength regimes. The result of this combination is an extremely high volume of data to be analysed. In this work, we have opened up the way to real-time automated analysis of the image data-stream. We present a fully automated GPU-based machine-learning backed pipeline for analysis of radio images at multiple frequencies. It consists of four consecutive steps: quality control, source detection, association, flux measurement and physical parameter inference. At the end of the pipeline, an alert of a significant detection can be sent out and data will be saved for further investigation. In our current implementation, the entire pipeline analyses images from AARTFAAC (Amsterdam Astron Radio Transients Facility And Analysis Centre) with sizes upwards of $16 \times 1024 \times 1024$ pixels. AARTFAAC is a transient detector based on the Low-Frequency Array, an interferometer based in the Netherlands with stations across Europe. First results show that dispersed signals were found on which follow-up analysis can be performed. The identification and response to transients is the key science goal for AARTFAAC, and the current work brings us one step closer to meeting it.

GRB160203A is a high redshift long gamma-ray burst presenting a collection of unusual features in the afterglow light curve. We study its optical and X-ray data. We find this event to occur within a constant density medium during the first part of the afterglow. However, after 13 ks we spot some flaring activities in the optical and X-ray light curves. We explain these flares by fluctuation of densities of the surrounding medium. Other scenarios, such as energy injection from a magnetar or variation of microphysical parameters are not supported by the data. We tentatively link these fluctuations to an unusual host galaxy, with gas density similar to the Milky Way and a dense cocoon of matter around a stellar progenitor similar to a Wolf-Rayet star. A termination shock scenario is found to be less likely.

Observational astronomy of tidal disruption events (TDEs) began with the detection of X-ray flares from quiescent galaxies during the ROSAT all-sky survey of 1990-1991. The flares complied with theoretical expectations, having high peak luminosities ($L_{\rm x}$ up to $\ge4\times 10^{44}$ erg/s), a thermal spectrum with $kT\sim$few$\times10^5$ K, and a decline on timescales of months to years, consistent with a diminishing return of stellar debris to a black hole of mass $10^{6-8}$ solar masses. These measurements gave solid proof that the nuclei of quiescent galaxies are habitually populated by a super-massive black hole. Beginning in 2000, XMM-Newton, Chandra and Swift have discovered further TDEs which have been monitored closely at multiple wavelengths. A general picture has emerged of, initially near-Eddington accretion, powering outflows of highly-ionised material, giving way to a calmer sub-Eddington phase, where the flux decays monotonically, and finally a low accretion rate phase with a harder X-ray spectrum indicative of the formation of a disk corona. There are exceptions to this rule though which at the moment are not well understood. A few bright X-ray TDEs have been discovered in optical surveys but in general X-ray TDEs show little excess emission in the optical band, at least at times coincident with the X-ray flare. X-ray TDEs are powerful new probes of accretion physics down to the last stable orbit, revealing the conditions necessary for launching jets and winds. Finally we see that evidence is mounting for nuclear and non-nuclear intermediate mass black holes based on TDE flares which are relatively hot and/or fast.

Energetic particles, such as stellar cosmic rays, produced at a heightened rate by active stars (like the young Sun) may have been important for the origin of life on Earth and other exoplanets. Here we compare, as a function of stellar rotation rate ($\Omega$), contributions from two distinct populations of energetic particles: stellar cosmic rays accelerated by impulsive flare events and Galactic cosmic rays. We use a 1.5D stellar wind model combined with a spatially 1D cosmic ray transport model. We formulate the evolution of the stellar cosmic ray spectrum as a function of stellar rotation. The maximum stellar cosmic ray energy increases with increasing rotation i.e., towards more active/younger stars. We find that stellar cosmic rays dominate over Galactic cosmic rays in the habitable zone at the pion threshold energy for all stellar ages considered ($t_*=0.6-2.9\,$Gyr). However, even at the youngest age, $t_*=0.6\,$Gyr, we estimate that $\gtrsim\,80$MeV stellar cosmic ray fluxes may still be transient in time. At $\sim1\,$Gyr when life is thought to have emerged on Earth, we demonstrate that stellar cosmic rays dominate over Galactic cosmic rays up to $\sim$4$\,$GeV energies during flare events. Our results for $t_*=0.6\,$Gyr ($\Omega = 4\Omega_\odot$) indicate that $\lesssim$GeV stellar cosmic rays are advected from the star to 1$\,$au and are impacted by adiabatic losses in this region. The properties of the inner solar wind, currently being investigated by the Parker Solar Probe and Solar Orbiter, are thus important for accurate calculations of stellar cosmic rays around young Sun-like stars.

The Kelvin-Helmholtz instability (KHI) is a nonlinear shear-driven instability that develops at the interface between shear flows in plasmas. KHI has been inferred in various astrophysical plasmas and has been observed in situ at the magnetospheric boundaries of solar-system planets and through remote sensing at the boundaries of coronal mass ejections. While it was hypothesized to play an important role in the mixing of plasmas and in triggering solar wind fluctuations, its direct and unambiguous observation in the solar wind was still lacking. We report in-situ observations of ongoing KHI in the solar wind using Solar Orbiter during its cruise phase. The KHI is found in a shear layer in the slow solar wind in the close vicinity of the Heliospheric Current Sheet, with properties satisfying linear theory for its development. An analysis is performed to derive the local configuration of the KHI. A 2-D MHD simulation is also set up with empirical values to test the stability of the shear layer. In addition, magnetic spectra of the KHI event are analyzed. We find that the observed conditions satisfy the KHI onset criterion from the linear theory analysis, and its development is further confirmed by the simulation. The current sheet geometry analyses are found to be consistent with KHI development. Additionally, we report observations of an ion jet consistent with magnetic reconnection at a compressed current sheet within the KHI interval. The KHI is found to excite magnetic and velocity fluctuations with power-law scalings that approximately follow $k^{-5/3}$ and $k^{-2.8}$ in the inertial and dissipation ranges, respectively. These observations provide robust evidence of KHI development in the solar wind. This sheds new light on the process of shear-driven turbulence as mediated by the KHI with implications for the driving of solar wind fluctuations.

We present results from an analysis of $\sim$ 29,000 RR Lyrae stars located in the Large Magellanic Cloud (LMC). For these objects, near-infrared time-series photometry from the VISTA survey of the Magellanic Clouds system (VMC) and optical data from the OGLE (Optical Gravitational Lensing Experiment) IV survey and the Gaia Data Release 2 catalogue of confirmed RR Lyrae stars were exploited. Using VMC and OGLE IV magnitudes we derived period-luminosity (PL), period-luminosity-metallicity (PLZ), period-Wesenheit (PW) and period-Wesenheit-metallicity (PWZ) relations in all available bands. More that ~7,000 RR Lyrae were discarded from the analysis because they appear to be overluminous with respect to the PL relations. The $PL_{K_{\mathrm{s}}}$ relation was used to derive individual distance to $\sim 22,000$ RR Lyrae stars, and study the three-dimensional structure of the LMC. The distribution of the LMC RR Lyrae stars is ellipsoidal with the three axis $S_1$=6.5 kpc, $S_2$=4.6 kpc and $S_3$=3.7 kpc, inclination i=$22\pm4^{\circ }$ relative to the plane of the sky and position angle of the line of nodes $\theta=167\pm7^{\circ }$ (measured from north to east). The north-eastern part of the ellipsoid is closer to us and no particular associated substructures are detected as well as any metallicity gradient.

Isolated massive elliptical galaxies, or that are present at the center of cool-core clusters, are believed to be powered by hot gas accretion directly from their surrounding hot X-ray emitting gaseous medium. This leads to a giant Bondi-type spherical/quasi-spherical accretion flow onto their host SMBHs, with the accretion flow region extending well beyond the Bondi radius. In this work, we present a detailed study of Bondi-type spherical flow in the context of these massive ellipticals by incorporating the effect of entire gravitational potential of the host galaxy in the presence of cosmological constant $\Lambda$, considering a five-component galactic system (SMBH + stellar + dark matter + hot gas + $\Lambda$). The current work is an extension of Ghosh \& Banik (2015), who studied only the cosmological aspect of the problem. The galactic contribution to the potential renders the (adiabatic) spherical flow to become {\it multi-transonic} in nature, with the flow topology and flow structure significantly deviating from that of classical Bondi solution. More notably, corresponding to moderate to higher values of galactic mass-to-light ratios, we obtain Rankine-Hugoniot shocks in spherical wind flows. Galactic potential enhances the Bondi accretion rate. Our study reveals that there is a strict lower limit of ambient temperature below which no Bondi accretion can be triggered; which is as high as $\sim 9 \times 10^6 \, K$ for flows from hot ISM-phase, indicating that the hot phase tightly regulates the fueling of host nucleus. Our findings may have wider implications, particularly in the context of outflow/jet dynamics, and radio-AGN feedback, associated with these massive galaxies in the contemporary Universe.

We present a method for the in-flight relative flux self-calibration of a spectro-photometer instrument, general enough to be applied to any upcoming galaxy survey on satellite. The instrument response function, that accounts for a smooth continuous variation due to telescope optics, on top of a discontinuous effect due to the segmentation of the detector, is inferred with a $\chi^2$ statistics. The method provides unbiased inference of the sources count rates and of the reconstructed relative response function, in the limit of high count rates. We simulate a simplified sequence of observations with realistic distributions of sources and count rates, with the purpose of quantifying the relative importance of the number of sources and exposures for the correct reconstruction of the instrument response. We present a validation of the method, with the definition of figures of merit to quantify the expected performance, in plausible scenarios.

On October 1, 2019, the IceCube Collaboration detected a muon track neutrino with high probability of being of astrophysical origin, IC191001A. After a few hours, the tidal disruption event (TDE) AT2019dsg, observed by the Zwicky Transient Facility (ZTF), was indicated as the most likely counterpart of the IceCube track. More recently, the follow-up campaign of the IceCube alerts by ZTF suggested a second TDE, AT2019fdr, as a promising counterpart of another IceCube muon track candidate, IC200530A, detected on May 30, 2020. These are the second and third associations between astrophysical sources and high-energy neutrinos after the compelling identification of the blazar TXS 0506+056. Here, the search for ANTARES neutrinos from the directions of AT2019dsg and AT2019fdr using a time-integrated approach is presented. As no significant evidence for space clustering is found in the ANTARES data, upper limits on the one-flavour neutrino flux and fluence are set.

Modern wide-field high-cadence surveys have revealed the significant diversity of optical transient phenomena in their luminosity and timescale distributions, which led to the discovery of some mysterious fast optical transients (FOTs). These FOTs can usually rise and decline remarkably in a timescale of a few days to weeks, which are obviously much rapider than ordinary supernovae. SN 2019bkc/ATLAS19dqr is one of the fastest detected FOTs so far and, meanwhile, it was found to be un-associated with a host galaxy. These discoveries provide a good chance to explore the possible origins of FOTs. So, we model the light curves of SN 2019bkc in details. It is found that SN 2019bkc can be well explained by the thermal emission of an explosion ejecta that is powered by a long-lasting central engine. The engine could be a spinning-down millisecond magnetar or a fallback accretion onto a compact object. Combining the engine property, the mass of the ejecta, and the hostlessness of SN 2019bkc, we suggest that this FOT is likely to originate from a merger of a white dwarf and a neutron star.

Artificial neural networks (ANN) have been successfully used in the last years to identify patterns in astronomical images. The use of ANN in the field of asteroid dynamics has been, however, so far somewhat limited. In this work we used for the first time ANN for the purpose of automatically identifying the behaviour of asteroid orbits affected by the M1:2 mean-motion resonance with Mars. Our model was able to perform well above 85% levels for identifying images of asteroid resonant arguments in term of standard metrics like accuracy, precision and recall, allowing to identify the orbital type of all numbered asteroids in the region. Using supervised machine learning methods, optimized through the use of genetic algorithms, we also predicted the orbital status of all multi-opposition asteroids in the area. We confirm that the M1:2 resonance mainly affects the orbits of the Massalia, Nysa, and Vesta asteroid families.

Solar activity in all its varied manifestations is driven by the magnetic field. Particularly important for many purposes are two global quantities, the Sun's total and open magnetic flux, which can be computed from sunspot number records using models. Such sunspot-driven models, however, do not take into account the presence of magnetic flux during grand minima, such as the Maunder minimum. Here we present a major update of a widely used simple model, which now takes into account the observation that the distribution of all magnetic features on the Sun follows a single power law. The exponent of the power law changes over the solar cycle. This allows for the emergence of small-scale magnetic flux even when no sunspots are present for multiple decades and leads to non-zero total and open magnetic flux also in the deepest grand minima, such as the Maunder minimum, thus overcoming a major shortcoming of the earlier models. The results of the updated model compare well with the available observations and reconstructions of the solar total and open magnetic flux. This opens up the possibility of improved reconstructions of sunspot number from time series of cosmogenic isotope production rate.

An improved high-resolution ground-to-thermosphere version of the Institut Pierre-Simon Laplace (IPSL) Venus General Circulation Model (VGCM), including non-orographic gravity waves (GW) parameterization and fine-tuned non-LTE parameters, is presented here. We focus on the validation of the model built from a collection of data mostly from Venus Express experiments and coordinated ground-based telescope campaigns, in the upper mesosphere/lower thermosphere of Venus. These simulations result in an overall better agreement with temperature observations above 90 km, compared with previous versions of the VGCM. Density of CO2 and light species (CO and O) are also comparable with observations in terms of trend and order of magnitude. Systematic biases in the temperature structure are found at about 80-100 km and above 130 km at the terminator, possibly due to assumptions on the CO2 mixing ratio made for stellar/solar occultation retrievals and uncertainties on the collisional rate coefficients used in the non-LTE parameterization, respectively. Diurnal and latitudinal distribution of dynamical tracers are also analyzed. Overall, our simulations indicate that a weak westward retrograde wind is present up to about 120 km, producing the CO bulge displacement toward 2h-3h in the morning, instead of piling up at the anti-solar point, as for an idealised SS-AS circulation. This retrograde imbalance is suggested to be produced by perturbations of about 5 days Kelvin wave impacting the mesosphere up to 110 km (see the companion paper Navarro et al. 2021), combined with GW westward acceleration mostly above 110 km. On the whole, these model developments point to the importance of the inclusion of the lower atmosphere, higher resolution and finely tuned parameterizations in GCM of the Venusian upper atmosphere, in order to shed light on existing observations.

(Abridged) Short-period gas giant exoplanets are susceptible to intense atmospheric escape due to their large scale heights and strong high-energy irradiation. This process is thought to occur ubiquitously, but to date we have only detected direct evidence of atmospheric escape in hot Jupiters and warm Neptunes. The paucity of cases for intermediate, Saturn-sized exoplanets at varying levels of irradiation precludes a detailed understanding of the underlying physics in atmospheric escape of hot gas giants. Our objectives here are to assess the high-energy environment of the warm ($T_\mathrm{eq} = 970$ K) Saturn WASP-29 b and search for signatures of atmospheric escape. We used far-ultraviolet (FUV) observations from the Hubble Space Telescope to analyze the flux time series of H I, C II, Si III, Si IV, and N V during the transit of WASP-29 b. At 3$\sigma$ confidence, we rule out any in-transit absorption of H Ilarger than 92% in the Lyman-$\alpha$ blue wing and 19% in the red wing. We found an in-transit flux decrease of $39\%^{+12\%}_{-11\%}$ in the ground-state C II emission line at 133.45 nm. But due to moderate variability in the other C, N and Si lines, it is difficult to attribute a planetary or stellar origin for the ground-state C II signal. We place 3$\sigma$ absorption upper limits of 40%, 49% and 24% for Si III, Si IV, and for excited-state C II at 133.57 nm, respectively. Low activity levels and the faint X-ray luminosity suggest that WASP-29 is an old, inactive star. An energy-limited approximation combined with the reconstructed EUV spectrum of the host suggests that the planet is losing its atmosphere at a rate of $4 \times 10^9$ g s$^{-1}$. The non-detection at Lyman-$\alpha$ could be partly explained by a low fraction of escaping neutral hydrogen, or by the state of fast radiative blow-out we infer from the reconstructed stellar Lyman-$\alpha$ line.

The cool white dwarf WD 1856+534 was found to be transited by a Jupiter-sized object with a mass at or below 14 M$_{\rm{Jup}}$. We used the GTC telescope to obtain and analyse photometry and low resolution spectroscopy of six transits of WD 1856+534 b, with the intention to derive the slope of the transmission spectrum, towards an eventual detection of Rayleigh scattering of the particles in its atmosphere. Such a slope, assuming a cloud-free atmosphere dominated by Rayleigh scattering, could be translated into an estimation of the mass of WD 1856+534 b. However, the resultant transmission spectrum is essentially flat, and therefore permits only the determination of lower mass limits of 2.4 M$_{\rm{Jup}}$ at the 2-$\sigma$ level, or 1.6 M$_{\rm{Jup}}$ at 3-$\sigma$. These limits have implications for some of the proposed formation scenarios for the object. We elaborate on the potential effects of clouds and hazes in our estimations, based on previous studies of Jupiter and Titan. In addition, we detected an H$\alpha$ absorption feature in the combined spectrum of the host white dwarf, that leads to the assignation of a DA classification and allows derivation of an independent set of atmospheric parameters. Furthermore, the epochs of five transits were measured with sub-second precision, which demonstrates that additional objects more massive than $\approx$5 M$_{\rm{Jup}}$ and with periods longer than $O(100)$ days could be detected through the light travel time effect

Deuterium represents the only bound isotope in the universe with atomic mass number $A=2$. Motivated by the possibility of other universes, where the strong force could be stronger, this paper considers the effects of bound diprotons and dineutrons on stars. We find that the existence of additional stable nuclei with $A=2$ has relatively modest effects on the universe. Previous work indicates that Big Bang Nucleosynthesis (BBN) produces more deuterium, but does not lead to catastrophic heavy element production. This paper revisits BBN considerations and confirms that the universe is left with an ample supply of hydrogen and other light nuclei for typical cosmological parameters. Using the $MESA$ numerical package, we carry out stellar evolution calculations for universes with stable diprotons, with nuclear cross sections enhanced by large factors $X$. This work focuses on $X=10^{15}-10^{18}$, but explores the wider range $X$ = $10^{-3}-10^{18}$. For a given stellar mass, the presence of stable diprotons leads to somewhat brighter stars, with the radii and photospheric temperatures roughly comparable to thoese of red giants. The central temperature decreases from the characteristic value of $T_c\approx1.5\times10^7$ K for hydrogen burning down to the value of $T_c\approx10^6$ K characteristic of deuterium burning. The stellar lifetimes are smaller for a given mass, but with the extended possible mass range, the smallest stars live for trillions of years, far longer than the current cosmic age. Finally, the enhanced cross sections allow for small, partially degenerate objects with mass $M_\ast=1-10M_J$ to produce significant steady-state luminosity and thereby function as stars.

How has the solar wind evolved to reach what it is today? In this review, I discuss the long-term evolution of the solar wind, including the evolution of observed properties that are intimately linked to the solar wind: rotation, magnetism and activity. Given that we cannot access data from the solar wind 4 billion years ago, this review relies on stellar data, in an effort to better place the Sun and the solar wind in a stellar context. I overview some clever detection methods of winds of solar-like stars, and derive from these an observed evolutionary sequence of solar wind mass-loss rates. I then link these observational properties (including, rotation, magnetism and activity) with stellar wind models. I conclude this review then by discussing implications of the evolution of the solar wind on the evolving Earth and other solar system planets. I argue that studying exoplanetary systems could open up new avenues for progress to be made in our understanding of the evolution of the solar wind.

The InfraRed Imaging Spectrograph (IRIS) is a first-light instrument for the Thirty Meter Telescope (TMT) that will be used to sample the corrected adaptive optics field by the Narrow-Field Infrared Adaptive Optics System (NFIRAOS) with a near-infrared (0.8 - 2.4 $\mu$m) imaging camera and integral field spectrograph. To better understand IRIS science specifications we use the IRIS data simulator to characterize relative photometric precision and accuracy across the IRIS imaging camera 34"x34" field of view. Because the Point Spread Function (PSF) varies due to the effects of anisoplanatism, we use the Anisoplanatic and Instrumental Reconstruction of Off-axis PSFs for AO (AIROPA) software package to conduct photometric measurements on simulated frames using PSF-fitting as the PSF varies in single-source, binary, and crowded field use cases. We report photometric performance of the imaging camera as a function of the instrumental noise properties including dark current and read noise. Using the same methods, we conduct comparisons of photometric performance with reconstructed PSFs, in order to test the veracity of the current PSF-Reconstruction algorithms for IRIS/TMT.

Abstract The microquasar SS 433 exhibits in H alpha intermittent flares, Doppler shifted to both the red and the blue. The mean remembers the orbital phase of the compact object. I show that the flares are not intermittent sightings of an accretion disk; rather, plasma must be expelled through the L2 point, thus remembering the phase of the orbit as it invades the space beyond the system. That space has been mapped with GRAVITY observations of a similar flare, revealing a strong rotation component.

The Gaia optical reference frame is intrinsically undefined with respect to global orientation and spin, so it needs to be anchored in the radio-based International Celestial Reference Frame (ICRF) to provide a referenced and quasi-inertial celestial coordinate system. The link between the two fundamental frames is realized through two samples of distant extragalactic sources, mostly AGNs and quasars, but only the smaller sample of radio-loud ICRF sources with optical counterparts is available to determine the mutual orientation. The robustness of this link can be mathematically formulated in the framework of functional principal component analysis using a set of vector spherical harmonics to represent the differences in celestial positions of the common objects. The weakest eigenvectors are computed, which describe the greatest deficiency of the link. The deficient or poorly determined terms are specific vector fields on the sphere which carry the largest errors of absolute astrometry using Gaia in reference to the ICRF. This analysis provides guidelines to the future development of the ICRF maximizing the accuracy of the link over the entire celestial sphere. A measure of robustness of a least-squares solution, which can be applied to any linear model fitting problem, is introduced to help discriminate between reference frame tie models of different degrees.

The cosmological constant if considered as a fundamental constant, provides an information treatment for gravitation problems, both cosmological and of black holes. The efficiency of that approach is shown via gedanken experiments for the information behavior of the horizons for Schwarzschild-de Sitter and Kerr-de Sitter metrics. A notion of entropy regarding any observer and in all possible non-extreme black hole solutions is suggested, linked also to Bekenstein bound. The suggested information approach forbids the existence of naked singularities.

Conventional indirect dark matter (DM) searches look for an excess in the electromagnetic emission from the sky that cannot be attributed to known astrophysical sources. Here, we argue that the photon polarisation is an important feature to understand new physics interactions and can be exploited to improve our sensitivity to DM. In particular, circular polarisation can be generated from Beyond the Standard Model interactions if they violate parity and there is an asymmetry in the number of particles which participate in the interaction. In this work, we consider a simplified model for fermionic (Majorana) DM and study the circularly polarised gamma rays below 10 GeV from the scattering of cosmic ray electrons on DM. We calculate the differential flux of positive and negative polarised photons from the Galactic Center and show that the degree of circular polarization can reach up to 90%. Finally, once collider and DM constraints have been taken into account, we estimate the required sensitivity from future experiments to detect this signal finding that, although a distinctive peak will be present in the photon flux spectrum, a near future observation is unlikely. However, different sources or models not considered in this work could provide higher intensity fluxes, leading to a possible detection by e-ASTROGAM. In the event of a discovery, we argue that the polarisation fraction is a valuable characterisation feature of the new sector.

Knowledge about the shape of the mass spectrum of compact stars helps in breaking the degeneracy between the mass and redshift of the gravitational wave (GW) sources and can be used to infer cosmological parameters in the absence of redshift measurements obtained from electromagnetic observations. In this paper we study the achievable accuracy and the limits of this approach. We perform cosmological inference from GW data assuming parametric compact binary population models. We consider two representative models for the mass spectrum, namely a power-law model between two hard cut-offs at a minimum and maximum mass, and a similar model combined with a Gaussian peak. Both models exhibit characteristic scales that allow an indirect estimate of the source redshift. In the case of the LIGO-Virgo detector network with current and future sensitivities we perform the analysis of an extensive set of simulated data using a hierarchical Bayesian scheme that jointly fit the source population and cosmological parameter. We also re-analyse the LIGO-Virgo O2 data. Those analyses all evidence the tight interplay between source population and cosmological parameters and the influence of initial assumptions formulated on the ones or the others. We find that: $(i)$ the upper mass cut-off and the position of the Gaussian peak display the largest correlation with the cosmological parameters; $(ii)$ incorrect population models may bias the Hubble constant estimate by 40%, or incorrect value for $\Omega_{m,0}$ may lead a significant bias on $H_0$; $(iii)$ the estimates enter the large sample regime with asymptotic normality and $1/\sqrt{N}$ error decay from $N\sim200$ GW events. Overall, our results suggest that the inference on source population and cosmological parameters should be addressed \textit{jointly} and not separately as it is in most studies so far.

In the light of the recently observed XENON1T electronic recoil (ER) data, we investigate the possibility of constraining the parameter space of a generic fermionic warm dark matter (WDM), decaying into a standard model (SM) neutrino and a photon. The photon as a decay product, when produced inside the XENON1T chamber, interacts with an electron of a xenon (Xe) atom, leading to a contribution in the observed ER data. We add this dark matter (DM) induced signal over four distinct backgrounds (taking one at a time) and perform a single parameter $\chi^2$ fit against the XENON1T data to obtain the best-fit values of the DM decay width and the associated 95$\%$ confidence level (C.L.) bands for DM mass varied in the range $2 - 18$ keV. By comparing the constraints, obtained by fitting the XENON1T data, with the upper limits, arising from various existing astrophysical and cosmological observations, we find that a fair amount of the DM parameter space is allowed at $95\%$ C.L., for each of the background models considered.

The aim of the present work is to investigate the effects of strong magnetic fields on the hadron-quark phase transition point at zero temperature. To describe the hadronic phase, a relativistic mean field (RMF) model is used and to describe the quark phase a density dependent quark mass model (DDQM) is employed. As compared with the results obtained with non-magnetised matter, we observe a shift of the transition point towards higher pressures and, generally also towards higher chemical potentials. An investigation of the phase transitions that could sustain hybrid stars is also performed.

The traditional approach to perturbations of nonrotating black holes in General Relativity uses the reformulation of the equations of motion into a radial second-order Schr\"odinger-like equation, whose asymptotic solutions are elementary. Imposing specific boundary conditions at spatial infinity and near the horizon defines, in particular, the quasi-normal modes of black holes. For more complicated equations of motion, as encountered for instance in modified gravity models with different background solutions and/or additional degrees of freedom, such a convenient Schr\"odinger-like reformulation might be unavailable, even in a generalised matricial form. In order to tackle such cases, we present a new approach that analyses directly the first-order differential system in its original form and extracts the asymptotic behaviour of perturbations. As a pedagogical illustration, we apply this treatment to the perturbations of Schwarzschild black holes and then show that the standard quasi-normal modes can be obtained numerically by solving this first-order system with a spectral method. This new approach paves the way for a generic treatment of the asymptotic behaviour of black hole perturbations and the identification of quasi-normal modes in theories of modified gravity.

We study the linear perturbations about nonrotating black holes in the context of degenerate higher-order scalar-tensor (DHOST) theories, using a systematic approach that extracts the asymptotic behaviour of perturbations (at spatial infinity and near the horizon) directly from the first-order radial differential system governing these perturbations. For axial (odd-parity) modes, this provides an alternative to the traditional approach based on a second-order Schr\"odinger-like equation with an effective potential, which we also discuss for completeness. By contrast, for polar (even-parity) modes, which contain an additional degree of freedom in DHOST theories, a similar generalised second-order Schr\"odinger-like matricial system does not seem available in general, which leaves only the option of a direct treatment of the four-dimensional first-order differential system. We illustrate our study with two specific types of black hole solutions: "stealth" Schwarzschild black holes, with a non trivial scalar hair, as well as a class of non-stealth black holes whose metric is distinct from Schwarzschild. The knowledge of the asymptotic behaviours of the pertubations enables us to compute numerically quasi-normal modes, as we show explicitly for the non-stealth solutions. Finally, the asymptotic form of the modes also signals some pathologies in the stealth and non-stealth solutions considered here.

Third-generation gravitational wave detectors, such as Einstein Telescope and Cosmic Explorer, will detect a bunch of gravitational wave (GW) signals originated from the coalescence of neutron star and black hole binary systems out to the higher redshifts, $z\sim 5-10$. There is a potential concern that some of the GW signals detected at a high statistical significance eventually overlap with each other, and the parameter estimation of such an overlapping system can differ from the one expected from the single event. Also, there is certainly the overlapping systems in which one of overlapping events has low signal-to-noise ratio unable to clearly detect. Those system will be potentially misidentified with a single GW event, and the estimated parameters of binary GWs can be biased. In this paper, we estimate the occurrence rate of those overlapping events, and study their statistical impacts on the parameter estimation based on the Fisher matrix analysis. Our finding is that the overlapping signals produce neither large statistical errors nor serious systematic biases on the parameters of binary systems, unless the coalescence time and the redshifted chirp masses of the two overlapping GWs are very close to each other. The occurrence rate of such a closely overlapping event is shown to be very small with the third-generation detectors.

In the modern search for life elsewhere in the Universe, we are broadly looking for the following: the planets similar to Earth - physical indicators of habitability, and the manifestation of life - the biological signatures. A biosignature is a measured parameter that has a high probability of being caused by the living organisms, either atmospheric gas species or some surface features. Therefore, the focus of a search is on a product or phenomena produced by the living systems, mostly by microorganisms as these are the most abundant on our planet like, say, methane. However, we may need to distinguish the terms `biosignature' and `bioindicator'. A biosignature is what living organisms produce - a bioproduct, while a bioindicator may be anything necessary for life as we know it, such as water or a rocky planet. Oxygen in this case is a double biomarker; first, it is a byproduct of oxygenic photosynthesis and, second, it is a signature of a complex life, because complex highly organized life requires high levels of oxygen. It is possible that there are other such bioindicators. For example, in the atmospheric compositions of terrestrial planets in our Solar System (including Titan), argon is one of the major constituents, moreover it was recently acknowledged to be a `biologically' active gas, exhibiting organprotective and neuroprotective properties, especially under hypoxic conditions. Here we propose that argon in the atmosphere of a rocky planet is a bioindicator of a highly organized life, provided that the planet is already deemed potentially habitable: with water, atmosphere, and of a certain age allowing for the complex life to evolve. We also delineate its possible detection methods.

In this paper, we show the equivalency between the existence of fast neutrino flavor instability and that of neutrino flavor lepton number (NFLN) crossings. The veracity of this proposition has been uncertain and sometimes controversial despite it being essential in the flavor evolutions of dense neutrinos. This study clarifies the conditions under which fast instability occurs and contributes to the elucidation of collective neutrino oscillations.

In this work, we investigate gravitational baryogenesis in the framework of $f(P)$ gravity to understand the applicability of this class of modified gravity in addressing the baryon asymmetry of the Universe. For the analysis, we set $f(P) = \alpha P$ where $\alpha$ is the model parameter. We found that in $f(P)$ gravity, the CP-violating interaction acquires a modification through the addition of the nontopological cubic term $P$ in addition to the Ricci scalar $R$ and the mathematical expression of the baryon-to-entropy ratio depends not only on the time derivative of $R$ but also the time derivative of $P$. Additionally, we also investigate the consequences of a more complete and generalized CP-violating interaction proportional to $f(P)$ instead of $P$ in addressing the baryon asymmetry of the Universe. For this type of interaction, we report that the baryon-to-entropy ratio is proportional to $\dot{R}$, $\dot{P}$ and $f^{'}(P)$. We report that for both of these cases, rational values of $\alpha$ and $\chi$ generate acceptable baryon-to-entropy ratios compatible with observations.

Dark energy is the constituent with an enormous abundance of the present universe, responsible for the universe's accelerated expansion. Therefore, it is plausible that dark energy may interact within any compact astrophysical objects. The author in Ref. [Phys. Rev. D 83, 127501 (2011)], constructs an exact star solution consisting of an ordinary matter and phantom field from a constant density star (CDS) known as Schwarzschild interior solution. The star denotes a dark energy star (DES). The author claims that the phantom field represents dark energy within the star. So far, the role of the phantom field as dark energy in DES is not systematically studied yet. Related to this issue, we analyze the energy condition of DES. We expect that DES shall violate the strong energy condition (SEC) for a particular condition. We discover that SEC is fully violated only when the compactness reaches the Buchdahl limit. Furthermore, we also investigate the causal conditions and stabilities due to the convective motion and gravitational cracking. We also find that those conditions are violated. These results indicate that DES is not physically stable. However, we may consider DES as an ultra-compact object of which we can calculate the gravitational wave echo time and echo frequency and compare them to those of CDS. We find that the contribution of the phantom field delays the gravitational wave echoes. The effective potential of the perturbed DES is also studied. The potential also enjoys a potential well like CDS but with a deeper well. We also investigate the possibility that DES could form a gravastar when $ C=1 $. It is found that gravastar produced from DES possesses no singularity with a dS-like phase as the interior. These results could open more opportunities for the observational study of dark energy in the near future, mostly from the compact astrophysical objects.

Vacuum-UV (VUV) photodesorption from water-rich ice mantles coating interstellar grains is known to play an important role in the gas-to-ice ratio in star- and planet-forming regions. Quantitative photodesorption yields from water ice are crucial for astrochemical models. We aim to provide the first quantitative photon-energy dependent photodesorption yields from water ice in the VUV. This information is important to understand the photodesorption mechanisms and to account for the variation of the yields under interstellar irradiation conditions. Experiments have been performed on the DESIRS beamline at the SOLEIL synchrotron, delivering tunable VUV light, using the SPICES (Surface Processes and ICES) set-up. Compact amorphous solid water ice (H$_2$O and D$_2$O) has been irradiated from 7 to 13.5 eV. Quantitative yields have been obtained by detection in the gas phase with mass-spectrometry for sample temperatures ranging from 15 K to 100 K. Photodesorption spectra of H$_2$O (D$_2$O), OH (OD), H$_2$ (D$_2$) and O$_2$ peak around 9-10 eV and decrease at higher energies. Average photodesorption yields of intact water at 15 K are 5 $\times$ 10$^{-4}$ molecule/photon for H$_2$O and 5 $\times$ 10$^{-5}$ molecule/photon for D$_2$O over the 7-13.5 eV range. The strong isotopic effect can be explained by a differential chemical recombination between OH (OD) and H (D) photofragments originating from lower kinetic energy available for the OH photofragments upon direct water photodissociation and/or possibly by an electronic relaxation process. It is expected to contribute to water fractionation during the building-up of the ice grain mantles in molecular clouds and to favor OH-poor chemical environment in comet-formation regions of protoplanetary disks. The yields of all the detected species except OH (OD) are enhanced above (70 $\pm$10) K, suggesting an ice restructuration at this temperature.

We examine the propagation and flavor oscillations of neutrinos under the influence of gravitational waves (GWs) with an arbitrary polarization. We rederive the effective Hamiltonian for the system of three neutrino flavors using the perturbative approach. Then, using this result, we consider the evolution of neutrino flavors in stochastic GWs with the general energy density spectrum.The equation for the density matrix for neutrino flavors is obtained and solved analytically. As an application, we study the flavor content of a neutrino beam emitted in a supernova type-II explosion. We obtain the analytical expressions for the contributions of GWs to the neutrino fluxes and for the damping decrement, which describes the attenuation of the fluxes to their asymptotic values. We are able to evaluate qualitatively the effect of various sources of stochastic GWs on the evolution of neutrino fluxes. We prove that the major contribution is by GWs emitted by merging supermassive black holes. The implication of the obtained results for the measurement of astrophysical neutrinos with neutrino telescopes is discussed.

We explore equilibrium solutions of spherically symmetric boson stars in the Palatini formulation of $f(\mathcal{R})$ gravity. We account for the modifications introduced in the gravitational sector by using a recently established correspondence between modified gravity with scalar matter and general relativity with modified scalar matter. We focus on the quadratic theory $f(\mathcal{R})=R+\xi R^2$ and compare its solutions with those found in general relativity, exploring both positive and negative values of the coupling parameter $\xi$. As matter source, a complex, massive scalar field with and without self-interaction terms is considered. Our results show that the existence curves of boson stars in Palatini $f(\mathcal{R})$ gravity are fairly similar to those found in general relativity. Major differences are observed for negative values of the coupling parameter which results in a repulsive gravitational component for high enough scalar field density distributions. Adding self-interactions makes the degeneracy between $f(\mathcal{R})$ and general relativity even more pronounced, leaving very little room for observational discrimination between the two theories.