Papers for Monday, Mar 01 2021

Accurate simulations of flows in stellar interiors are crucial to improving our understanding of stellar structure and evolution. Because the typically slow flows are but tiny perturbations on top of a close balance between gravity and pressure gradient, such simulations place heavy demands on numerical hydrodynamics schemes. We demonstrate how discretization errors on grids of reasonable size can lead to spurious flows orders of magnitude faster than the physical flow. Well-balanced numerical schemes can deal with this problem. Three such schemes are applied in the implicit, finite-volume code SLH in combination with a low-Mach-number numerical flux function. We compare how the schemes perform in four numerical experiments addressing some of the challenges imposed by typical problems in stellar hydrodynamics. We find that the $\alpha$-$\beta$ and Deviation well-balancing methods can accurately maintain hydrostatic solutions provided that gravitational potential energy is included in the total energy balance. They accurately conserve minuscule entropy fluctuations advected in an isentropic stratification, which enables the methods to reproduce the expected scaling of convective flow speed with the heating rate. The Deviation method also substantially increases accuracy of maintaining stationary orbital motions in a Keplerian disk on long time scales. The Cargo-LeRoux method fares substantially worse in our tests, although its simplicity may still offer some merits in certain situations. Overall, we find the well-balanced treatment of gravity in combination with low Mach number flux functions essential to reproducing correct physical solutions to challenging stellar slow-flow problems on affordable collocated grids.

We study the radial acceleration relation (RAR) between the total ($a_{\rm tot}$) and baryonic ($a_{\rm bary}$) centripetal acceleration profiles of central galaxies in the cold dark matter (CDM) paradigm. We analytically show that the RAR is intimately connected with the physics of the quasi-adiabatic relaxation of dark matter in the presence of baryons in deep potential wells. This cleanly demonstrates how a near-universal mean RAR and its scatter emerges in the low-acceleration regime ($10^{-12}\,{\rm m\,s}^{-2}\lesssim a_{\rm bary}\lesssim10^{-10}\,{\rm m\,s}^{-2}$) from an interplay between baryonic feedback processes and the distribution of CDM in dark halos. Our framework allows us to go further and study both higher and lower accelerations in detail, using analytical approximations and a realistic mock catalog of $\sim342,000$ low-redshift central galaxies with $M_r\leq-19$. We show that, while the RAR in the baryon-dominated, high-acceleration regime ($a_{\rm bary}\gtrsim10^{-10}\,{\rm m\,s}^{-2}$) is very sensitive to details of the relaxation physics, a simple `baryonification' prescription matching the relaxation results of hydrodynamical CDM simulations is remarkably successful in reproducing the observed RAR without any tuning. And in the (currently unobserved) ultra-low-acceleration regime ($a_{\rm bary}\lesssim 10^{-12}\,{\rm m\,s}^{-2}$), the RAR is sensitive to the abundance of diffuse gas in the halo outskirts, with our default model predicting a distinctive break from a simple power-law-like relation for HI-deficient, diffuse gas-rich centrals. Our mocks also show that the RAR provides more robust, testable predictions of the $\Lambda$CDM paradigm at galactic scales, with implications for alternative gravity theories, than the baryonic Tully-Fisher relation.

We present new WFC3/UVIS observations of UGC 4483, the closest example of a metal-poor blue compact dwarf galaxy, with a metallicity of $Z \simeq 1/15\ Z_{\odot}$ and located at a distance of $D \simeq 3.4$ Mpc. The extremely high quality of our new data allows us to clearly resolve the multiple stellar evolutionary phases populating the color-magnitude diagram (CMD), to reach more than 4 mag deeper than the tip of the red giant branch, and to detect for the first time core He-burning stars with masses $\lesssim 2$ M$_{\odot}$, populating the red clump and possibly the horizontal branch (HB) of the galaxy. By applying the synthetic CMD method to our observations, we determine an average star formation rate over the whole Hubble time of at least $(7.01\pm0.44) \times 10^{-4}$ $\mathrm{M_{\odot}/yr}$, corresponding to a total astrated stellar mass of $(9.60\pm0.61)\times 10^6$ $\mathrm{M_{\odot}}$, 87% of which went into stars at epochs earlier than 1 Gyr ago. With our star formation history recovery method we find the best fit with a distance modulus of DM = $27.45\pm0.10$, slightly lower than previous estimates. Finally, we find strong evidence of an old ($\gtrsim 10$ Gyr) stellar population in UGC 4483 thanks to the detection of an HB phase and the identification of six candidate RR Lyrae variable stars.

J-PAS will soon start imaging 8000 deg2 of the northern sky with its unique set of 56 filters (R $\sim$ 60). Before, we observed 1 deg2 on the AEGIS field with an interim camera with all the J-PAS filters. With this data (miniJPAS), we aim at proving the scientific potential of J-PAS to identify and characterize the galaxy populations with the goal of performing galaxy evolution studies across cosmic time. Several SED-fitting codes are used to constrain the stellar population properties of a complete flux-limited sample (rSDSS <= 22.5 AB) of miniJPAS galaxies that extends up to z = 1. We find consistent results on the galaxy properties derived from the different codes, independently of the galaxy spectral-type or redshift. For galaxies with SNR>=10, we estimate that the J-PAS photometric system allows to derive stellar population properties with a precision that is equivalent to that obtained with spectroscopic surveys of similar SNR. By using the dust-corrected (u-r) colour-mass diagram, a powerful proxy to characterize galaxy populations, we find that the fraction of red and blue galaxies evolves with cosmic time, with red galaxies being $\sim$ 38% and $\sim$ 18% of the whole population at z = 0.1 and z = 0.5, respectively. At all redshifts, the more massive galaxies belong to the red sequence and these galaxies are typically older and more metal rich than their counterparts in the blue cloud. Our results confirm that with J-PAS data we will be able to analyze large samples of galaxies up to z $\sim$ 1, with galaxy stellar masses above of log(M$_*$/M$_{\odot}$) $\sim$ 8.9, 9.5, and 9.9 at z = 0.3, 0.5, and 0.7, respectively. The SFH of a complete sub-sample of galaxies selected at z $\sim$ 0.1 with log(M$_*$/M$_{\odot}$) > 8.3 constrain the cosmic evolution of the star formation rate density up to z $\sim$ 3 in good agreement with results from cosmological surveys.

Galaxy clusters identified from the Sunyaev Zel'dovich (SZ) effect are a key ingredient in multi-wavelength cluster-based cosmology. We present a comparison between two methods of cluster identification: the standard Matched Filter (MF) method in SZ cluster finding and a method using Convolutional Neural Networks (CNN). We further implement and show results for a `combined' identifier. We apply the methods to simulated millimeter maps for several observing frequencies for an SPT-3G-like survey. There are some key differences between the methods. The MF method requires image pre-processing to remove point sources and a model for the noise, while the CNN method requires very little pre-processing of images. Additionally, the CNN requires tuning of hyperparameters in the model and takes as input, cutout images of the sky. Specifically, we use the CNN to classify whether or not an 8 arcmin $\times$ 8 arcmin cutout of the sky contains a cluster. We compare differences in purity and completeness. The MF signal-to-noise ratio depends on both mass and redshift. Our CNN, trained for a given mass threshold, captures a different set of clusters than the MF, some of which have SNR below the MF detection threshold. However, the CNN tends to mis-classify cutouts whose clusters are located near the edge of the cutout, which can be mitigated with staggered cutouts. We leverage the complementarity of the two methods, combining the scores from each method for identification. The purity and completeness of the MF alone are both 0.61, assuming a standard detection threshold. The purity and completeness of the CNN alone are 0.59 and 0.61. The combined classification method yields 0.60 and 0.77, a significant increase for completeness with a modest decrease in purity. We advocate for combined methods that increase the confidence of many lower signal-to-noise clusters.

PSR J1023+0038 is a rapidly-spinning neutron star with a low-mass-binary companion that switches between a radio pulsar and low-luminosity disk state. In 2013, it transitioned to its current disk state accompanied by brightening at all observed wavelengths. In this state, PSR J1023+0038 now shows optical and X-ray pulsations and abrupt X-ray luminosity switches between discrete 'low' and 'high' modes. Continuum radio emission, denoting an outflow, is also present and brightens during the X-ray low modes. Here, we present a simultaneous optical, ultraviolet (UV) and X-ray campaign comprising Kepler ($400-800$ nm), Hubble Space Telescope ($180-280$ nm), XMM-Newton ($0.3-10$ keV) and NuSTAR ($3 - 79$ keV). We demonstrate that low and high luminosity modes in the UV band are strictly simultaneous with the X-ray modes and change the UV brightness by a factor of $\sim25$\% on top of a much brighter persistent UV component. We find strong evidence for UV pulsations (pulse fraction of $0.82\pm0.19$\%) in the high-mode, with a similar waveform as the X-ray pulsations making it the first known UV millisecond pulsar. Lastly, we find that the optical mode changes occur synchronously with the UV/X-ray mode changes, but optical modes are inverted compared to the higher frequencies. There appear to be two broad-band emission components: one from radio to near-infrared/optical that is brighter when the second component from optical to hard X-rays is dimmer (and vice-versa). We suggest that these components trace switches between accretion into the neutron star magnetosphere (high-energy high-mode) versus ejection of material (low-energy high-mode). Lastly, we propose that optical/UV/X-ray pulsations can arise from a shocked accretion flow channeled by the neutron star's magnetic field.

Gas phase Elemental abundances in Molecular CloudS (GEMS) is an IRAM 30m Large Program designed to estimate the S, C, N, and O depletions and gas ionization degree, X(e-), in a set of star-forming filaments of Taurus, Perseus and Orion. Our immediate goal is to build up a complete database of molecular abundances that can serve as an observational basis for estimating X(e-) and the C, O, N, and S depletions through chemical modeling. We observed and derived the abundances of 14 species (13CO, C18O, HCO+, H13CO+, HC18O+, HCN, H13CN, HNC, HCS+, CS, SO, 34SO, H2S, and OCS) in 244 positions, covering the AV 3 to 100 mag, n(H2) a few 10$^{3}$ to 10$^6$ cm$^{-3}$, and Tk 10 to 30 K ranges in these clouds, avoiding protostars, HII regions, and outflows. A statistical analysis is carried out to identify general trends between different species and with physical parameters. Relations between molecules reveal strong linear correlations which define three different families: (1) 13CO and C18O; (2) H13CO+, HC18O+, H13CN, and HNC; and (3) the S-bearing molecules. The abundances of the CO isotopologs increase with the gas kinetic temperature until TK 15 K. For higher temperatures, the abundance remains constant with a scatter of a factor of 3. The abundances of H13CO+, HC18O+, H13CN, and HNC are well correlated with each other, and all of them decrease with molecular hydrogen density, following the law n(H2)$^{-0.8\pm0.2}$. The abundances of S-bearing species also decrease with n(H2) at a rate of (S-bearing/H)gas n(H2)$^{-0.6\pm0.1}$. The abundances of molecules belonging to groups 2 and 3 do not present any clear trend with gas temperature. At scales of molecular clouds, the C18O abundance is the quantity that better correlates with the cloud mass. We discuss the utility of the 13CO/C18O, HCO+/H13CO+, and H13CO+/H13CN abundance ratios as chemical diagnostics of star formation in external galaxies.

The soft MeV gamma-ray sky, from a few hundred keV up to several MeV, is one of the least explored regions of the electromagnetic spectrum. The most promising technology to access this energy range is a telescope that uses Compton scattering to detect the gamma rays. Going from the measured data to all-sky images ready for scientific interpretation, however, requires a well-understood detector setup and a multi-step data-analysis pipeline. We have developed these capabilities for the Compton Spectrometer and Imager (COSI). Starting with a deep understanding of the many intricacies of the Compton measurement process and the Compton data space, we developed the tools to perform simulations that match well with instrument calibrations and to reconstruct the gamma-ray path in the detector. Together with our work to create an adequate model of the measured background while in flight, we are able to perform spectral and polarization analysis, and create images of the gamma-ray sky. This will enable future telescopes to achieve a deeper understanding of the astrophysical processes that shape the gamma-ray sky from the sites of star formation (26-Al map), to the history of core-collapse supernovae (e.g. 60-Fe map) and the distributions of positron annihilation (511-keV map) in our Galaxy.

Long gamma-ray bursts (GRB), explosions of very massive stars, provide crucial information on stellar and galaxy evolution, even at redshifts z ~ 8 - 9.5, when the Universe was only 500-600 million years old. Recently, during observations of a galaxy at a redshift of z ~ 11 (400 million years after the Big Bang), a bright signal, named GN-z11-flash, shorter than 245 s was detected and interpreted as an ultraviolet flash associated with a GRB in this galaxy, or a shock-breakout in a Population III supernova. Its resulting luminosity would be consistent with that of other GRBs, but a discussion based on probability arguments started on whether this is instead a signal from a man-made satellite or a Solar System object. Here we show a conclusive association of GN-z11-flash with Breeze-M upper stage of a Russian Proton rocket on a highly elliptical orbit. This rules out GN-z11-flash as the most distant GRB ever detected. It also implies that monitoring of a larger sample of very high redshift galaxies is needed to detect such distant GRBs. This also highlights the importance of a complete database of Earth satellites and debris, which can allow proper interpretation of astronomical observations.

We perform the first suite of fully general relativistic magnetohydrodynamic simulations of spinning massive black hole binary mergers. We consider binary black holes with spins of different magnitudes aligned to the orbital angular momentum, which are immersed in a hot, magnetized gas cloud. We investigate the effect of the spin and degree of magnetization (defined through the fluid parameter $\beta^{-1}\equiv p_{\mathrm{mag}}/p_{\mathrm{fluid}}$) on the properties of the accretion flow. We find that magnetized accretion flows are characterized by more turbulent dynamics, as the magnetic field lines are twisted and compressed during the late inspiral. Post-merger, the polar regions around the spin axis of the remnant Kerr black hole are magnetically dominated, and the magnetic field strength is increased by a factor $\sim$10$^2$ (independently from the initial value of $\beta^{-1}$). The magnetized gas in the equatorial plane acquires higher angular momentum, and settles in a thin circular structure around the black hole. We find that mass accretion rates of magnetized configurations are generally smaller than in the unmagnetized cases by up to a factor $\sim$3. Black hole spins have also a suppressing effect on the accretion rate, as large as $\sim$48\%. As a potential driver for electromagnetic emission we follow the evolution of the Poynting luminosity, which increases after merger up to a factor $\sim2$ with increasing spin, regardless of the initial level of magnetization of the fluid. Our results stress the importance of taking into account both spins and magnetic fields when studying accretion processes onto merging massive black holes.

Several attempts have been made to observe whether solar flares excite acoustic modes since Wolff (1972) suggested this possibility. Moreover, the rapid progress of asteroseismology and the study of stellar flares makes the study of these phenomena in the Sun important to inform our study of the influence of the more energetic stellar flares on asteroseismic acoustic modes. We look for the impact of flares on the amplitude of solar acoustic modes and other effects that are also affecting the mode amplitude. Solar acoustic mode amplitudes are known to be sensitive to magnetic fields. As flares usually occur in the presence of strong magnetic fields and most likely are the by-product of magnetic reconnection, we show how the magnetic field in and around the flaring region affects the mode amplitude. The mode amplitudes were obtained using ring-diagram analysis, which was first applied to a single event, the largest flare in the space age (the `Halloween Flare', SOL2003-10-28T11:00), using MDI data. Then, using HMI data, the analysis was applied to the regions corresponding to the flares observed during the high activity phase of cycle 24 and that fall into two groups. These two groups consist of small (10-60 erg cm$^{-2}$ s$^{-1}$) and large ($>$1200 erg cm$^{-2}$ s$^{-1}$) peak-flux flares, based on the Heliophysics Event Knowledgebase (HEK). After applying several corrections in order to take into account several sources of bias, we did not find any amplification in the inferred mode amplitude due to flaring activity, within a 10% uncertainty.

Understanding the transport of cosmic rays is challenging our models of propagation in the Galaxy. A good characterization of the secondary cosmic rays (B, Be, Li and sub-iron species) is crucial to constrain these models and exploit the precision of modern CR experiments. In this work, a Markov chain Monte Carlo analysis has been implemented to fit the experimental flux ratios between B, Be and Li and their flux ratios to the primary cosmic-ray nuclei C and O. We have fitted the data using two different parametrizations for the spallation cross sections. The uncertainties in the evaluation of the spectra of these secondary cosmic rays, due to spallation cross sections, have been taken into account introducing a scale factor as a nuisance parameter in the fits, assuming that this uncertainty is mostly due to the normalization of the cross sections parametrizations. We have also tested two different kind of diffusion coefficients, which differ in the origin of the high energy hardening ($\sim 200$ GeV/n) of cosmic rays. Additionally, two different approaches are used to scale the cross sections, one based on a combined analysis of all the species ("combined" analysis) and the other reproducing the high energy spectra of the secondary-to-secondary flux ratios of Be/B, Li/B, Li/Be ("corrected" analysis). This allows us to make a better comparison between the propagation parameters inferred from both cross sections parametrizations. This novel analysis has been successfully implemented using the numerical code DRAGON2 dedicated to cosmic-ray propagation to reproduce the cosmic-ray nuclei data up to $Z=14$ from the AMS-02 experiment. We report the main results, comparing the different cross sections parametrizations and discussing the impact of these uncertainties.

In our earlier study of this series (Park et al. 2020, Paper I), we examined the hemispheric sign preference (HSP) of magnetic helicity flux $dH/dt$ across photospheric surfaces of 4802 samples of 1105 unique active regions (ARs) observed during solar cycle 24. Here, we investigate any association of the HSP, expressed as a degree of compliance, with flaring activity, analyzing the same set of $dH/dt$ estimates as used in Paper I. The AR samples under investigation are assigned to heliographic regions (HRs) defined in the Carrington longitude-latitude plane with a grid spacing of 45$^\circ$ in longitude and 15$^\circ$ in latitude. For AR samples in each of the defined HRs, we calculate the degree of HSP compliance and the average soft X-ray flare index. The strongest flaring activity is found to be in one distinctive HR with an extremely low HSP compliance of 41% as compared to the mean and standard deviation of 62% and 7%, respectively, over all HRs. This sole HR shows an anti-HSP (i.e., less than 50%) and includes the highly flare-productive AR NOAA 12673, however this AR is not uniquely responsible for the HR's low HSP. We also find that all HRs with the highest flaring activity are located in the southern hemisphere, and they tend to have lower degrees of HSP compliance. These findings point to the presence of localized regions of the convection zone with enhanced turbulence, imparting a greater magnetic complexity and a higher flaring rate to some rising magnetic flux tubes.

Coronal mass ejections (CMEs) and coronal jets are two types of common solar eruptive phenomena, which often independently happen at different spatial scales. In this work, we present a stereoscopic observation of a large-scale CME flux rope arising from an unwinding blowout jet in a multipolar complex magnetic system. Based on a multi-band observational analysis, we find that this whole event starts with a small filament whose eruption occurs at a coronal geyser site after a series of homologous jets. Aided by magnetic field extrapolations, it reveals that the coronal geyser site forms above an elongate opposite-polarity interface, where the emergence-driven photospheric flux cancellation and repetitive reconnection are responsible for those preceding recurrent jets and also contribute to the ultimate filament destabilization. By interacting with overlying fields, the erupting filament breaks one of its legs and results in an unwinding blowout jet. Our estimation suggests that around 1.4$-$2.0 turns of twist release in its jet spire. This prominent twist transport in jet spire rapidly creates a newborn larger-scale flux rope from the jet base to a remote site. Soon after its formation, this large-scale flux rope erupts towards the outer coronae causing an Earth-directed CME. In its source region, two sets of distinct post-flare loops form in succession, indicating this eruption involves two-stage of flare magnetic reconnection. This work not only reveals a real magnetic coupling process between different eruptive activities but provides a new hint for understanding the creation of large-scale CME flux ropes during the solar eruption.

Galaxy internal structure growth has long been accused of inhibiting star formation in disc galaxies. We investigate the potential physical connection between the growth of dispersion-supported stellar structures (e.g. classical bulges) and the position of galaxies on the star-forming main sequence at $z\sim0$. Combining the might of the SAMI and MaNGA galaxy surveys, we measure the $\lambda_{Re}$ spin parameter for 3781 galaxies over $9.5 < \log M_{\star} [\rm{M}_{\odot}] < 12$. At all stellar masses, galaxies at the locus of the main sequence possess $\lambda_{Re}$ values indicative of intrinsically flattened discs. However, above $\log M_{\star}[\rm{M}_{\odot}]\sim10.5$ where the main sequence starts bending, we find tantalising evidence for an increase in the number of galaxies with dispersion-supported structures, perhaps suggesting a connection between bulges and the bending of the main sequence. Moving above the main sequence, we see no evidence of any change in the typical spin parameter in galaxies once gravitationally-interacting systems are excluded from the sample. Similarly, up to 1 dex below the main sequence, $\lambda_{Re}$ remains roughly constant and only at very high stellar masses ($\log M_{\star}[\rm{M}_{\odot}]>11$), do we see a rapid decrease in $\lambda_{Re}$ once galaxies decline in star formation activity. If this trend is confirmed, it would be indicative of different quenching mechanisms acting on high- and low-mass galaxies. The results suggest that while a population of galaxies possessing some dispersion-supported structure is already present on the star-forming main sequence, further growth would be required after the galaxy has quenched to match the kinematic properties observed in passive galaxies at $z\sim0$.

We present Tails, an open-source deep-learning framework for the identification and localization of comets in the image data of the Zwicky Transient Facility (ZTF), a robotic optical time-domain survey currently in operation at the Palomar Observatory in California, USA. Tails employs a custom EfficientDet-based architecture and is capable of finding comets in single images in near real time, rather than requiring multiple epochs as with traditional methods. The system achieves state-of-the-art performance with 99% recall, 0.01% false positive rate, and 1-2 pixel root mean square error in the predicted position. We report the initial results of the Tails efficiency evaluation in a production setting on the data of the ZTF Twilight survey, including the first AI-assisted discovery of a comet (C/2020 T2) and the recovery of a comet (P/2016 J3 = P/2021 A3).

Context. The diffusion of volatile species on amorphous solid water ice affects the chemistry on dust grains in the interstellar medium as well as the trapping of gases enriching planetary atmospheres or present in cometary material. Aims. The aim of the work is to provide diffusion coefficients of CH$_4$ on amorphous solid water (ASW), and to understand how they are affected by the ASW structure. Methods. Ice mixtures of H$_2$O and CH$_4$ were grown in different conditions and the sublimation of CH$_4$ was monitored via infrared spectroscopy or via the mass loss of a cryogenic quartz crystal microbalance. Diffusion coefficients were obtained from the experimental data assuming the systems obey Fick's law of diffusion. Monte Carlo simulations modeled the different amorphous solid water ice structures investigated and were used to reproduce and interpret the experimental results. Results. Diffusion coefficients of methane on amorphous solid water have been measured to be between 10$^{-12}$ and 10$^{-13}$ cm$^2$ s$^{-1}$ for temperatures ranging between 42 K and 60 K. We showed that diffusion can differ by one order of magnitude depending on the morphology of amorphous solid water. The porosity within water ice, and the network created by pore coalescence, enhance the diffusion of species within the pores.The diffusion rates derived experimentally cannot be used in our Monte Carlo simulations to reproduce the measurements. Conclusions. We conclude that Fick's law can be used to describe diffusion at the macroscopic scale, while Monte Carlo simulations describe the microscopic scale where trapping of species in the ices (and their movement) is considered.

The Extreme-ultraviolet Variability Experiment (EVE) onboard the Solar Dynamics Observatory (SDO) detects the solar EUV spectra with high temporal cadence and spectral resolution. The wavelength shifts of emission lines provide key information of the dynamics of the Sun. However, some of EVE spectral observations are influenced by the non-uniformly distributed irradiance on the Sun, which may prevent us from correctly understanding the physical processes happened in the solar corona. Here, based on the only published on-orbit calibration data of EVE He II 30.38 nm line on 27 Jan 2011 (Chamberlin, 2016), we develop a method to correct the He II 30.38 nm line by using AIA 304 imaging data. This correction method is then applied to EVE He II 30.38 nm data from 29 Oct 2010 to 3 Mar 2011 to study the Doppler oscillations of the solar He II 30.38 nm line, in which we show that the half-month periodic Doppler oscillation is caused by non-uniformly distributed irradiance mainly due to the presence of active regions. Other EVE coronal lines also present similar Doppler oscillations, suggesting that an appropriate correction must be implemented before interpret the oscillation phenomena appearing in these lines.

King 2, one of the oldest clusters in the Milky Way, with an age of $\sim$ 6 Gyr and distance of $\sim 5700$ pc, has been observed with UVIT payload on the \textit{ASTROSAT}. With membership information derived from {\it Gaia} EDR3, the cluster is found to have 39 blue straggler stars (BSSs). We created multi-wavelength spectra-energy distributions (SED) of all the BSSs. Out of 10 UV detected BSSs, 6 bright ones fitted with double component SEDs and were found to have hotter companions with properties similar to extreme horizontal branch (EHB)/subdwarf B (sdB) stars, with a range in luminosity and temperature, suggesting a diversity among the hot companions. We suggest that at least 15\% of BSSs in this cluster are formed via mass-transfer pathway. When we compared their properties to EHBs and hotter companions to BSS in open and globular clusters, we suggest that EHB/sdBs like companions can form in binaries of open clusters as young as 6 Gyr.

Using high-resolution Hubble Space Telescope imaging data, we perform a visual morphological classification of $\sim 36,000$ galaxies at $z < 1$ in the DEVILS/COSMOS region. As the main goal of this study, we derive the stellar mass function (SMF) and stellar mass density (SMD) sub-divided by morphological types. We find that visual morphological classification using optical imaging is increasingly difficult at $z > 1$ as the fraction of irregular galaxies and merger systems (when observed at rest-frame UV/blue wavelengths) dramatically increases. We determine that roughly two-thirds of the total stellar mass of the Universe today was in place by $z \sim 1$. Double-component galaxies dominate the SMD at all epochs and increase in their contribution to the stellar mass budget to the present day. Elliptical galaxies are the second most dominant morphological type and increase their SMD by $\sim 2.5$ times, while by contrast, the pure-disk population significantly decreases by $\sim 85\%$. According to the evolution of both high- and low-mass ends of the SMF, we find that mergers and in-situ evolution in disks are both present at $z < 1$, and conclude that double-component galaxies are predominantly being built by the in-situ evolution in disks (apparent as the growth of the low-mass end with time), while mergers are likely responsible for the growth of ellipticals (apparent as the increase of intermediate/high-mass end).

Results are reported from the first full-scale search for transient signals from exotic fields of astrophysical origin using data from a newly constructed Earth-scale detector: the Global Network of Optical Magnetometers for Exotic physics searches (GNOME). Data collected by the GNOME consist of correlated measurements from optical atomic magnetometers located in laboratories all over the world. GNOME data are searched for patterns of signals propagating through the network consistent with exotic fields composed of ultralight bosons such as axion-like particles (ALPs). Analysis of data from a continuous month-long operation of the GNOME finds no statistically significant signals consistent with those expected due to encounters with topological defects (axion domain walls), placing new experimental constraints on such dark matter scenarios.

Coronal holes are large-scale structures in the solar atmosphere that feature a reduced temperature and density in comparison to the surrounding quiet Sun and are usually associated with open magnetic fields. We perform a differential emission measure analysis on the 707 non-polar coronal holes collected in the Collection of Analysis Tools for Coronal Holes (CATCH) catalog to derive and statistically analyze their plasma properties (i.e. temperature, electron density, and emission measure). We use intensity filtergrams of the six coronal EUV filters from the \textit{Atmospheric Imaging Assembly} onboard of the \textit{Solar Dynamics Observatory}, which cover a temperature range from $ \approx 10^{5.5}$ to $10^{7.5}$\,K. Correcting the data for stray and scattered light, we find that all coronal holes have very similar plasma properties with an average temperature of $0.94 \pm 0.18$ MK, a mean electron density of $(2.4 \pm 0.7) \times 10^{8}$\,cm$^{-3}$, and a mean emission measure of $(2.8 \pm 1.6) \times 10^{26}$\,cm$^{-5}$. The temperature distribution within the coronal hole was found to be largely uniform, whereas the electron density shows a $40\,\%$ linear decrease from the boundary towards the inside of the coronal hole. At distances greater than \SI{20}{\arcsecond} ($\approx 15$\,Mm) from the nearest coronal hole boundary, the density also becomes statistically uniform. The coronal hole temperature may show a weak solar cycle dependency, but no statistically significant correlation of plasma properties to solar cycle variations could be determined throughout the observed time period between 2010 and 2019.

We present the photometric and spectroscopic analysis of four W UMa binaries J015829.5+260333 (hereinafter as J0158), J030505.1+293443 (hereinafter as J0305), J102211.7+310022 (hereinafter as J1022) and KW Psc. The VRcIc band photometric observations are carried out with the 1.3-m Devasthal Fast Optical Telescope. For low resolution spectroscopy, we used 2-m Himalayan Chandra Telescope as well as the archival data from 4-m LAMOST survey. The systems J0158 and J0305 show a period increase rate of 5.26(+/-1.72)x10^-7 days/yr and 1.78(+/-1.52)x10^-6 days/yr, respectively. The period of J1022 is found to be decreasing with a rate of 4.22(+/-1.67)x10^-6 days/yr. The period analysis of KW Psc displays no change in its period. PHOEBE package is used for the light curve modeling and basic parameters are evaluated with the help of GAIA parallax. The asymmetry of light curves is explained with the assumption of cool spots at specific positions on one of the components of the system. On the basis of temperatures, mass ratios, fill-out factors and periods, the system J1022 is identified as W-subtype systems while the others show some mixed properties. To probe the chromospheric activities in these W UMa binaries, their spectra are compared with the known inactive stars spectra. The comparison shows emission in H\alpha, H\beta and CaII. To understand the evolutionary status of these systems, the components are drawn in mass-radius and mass-luminosity planes with other well characterized binary systems. The secondary components of all the systems are away from ZAMS which indicates that secondary is more evolved than the primary component.

Interstellar dust (ISD) penetrates into the heliosphere due to the relative motion of the Sun and the local interstellar medium (LISM). Inside the heliosphere and at the boundaries, where solar wind interacts with the LISM, distribution of ISD is modified due to the action of the electromagnetic forces, the solar gravitation and the radiation pressure. These forces make the distribution of the ISD particles in the heliosphere inhomogeneous. In previous work we demonstrated the existence of singularities in the ISD density distribution at 0.03 - 10 a.u. north and south with respect to the heliospheric current sheet. In this paper we show that dispersion in the ISD velocity distribution strongly affects the singularities. Even small values of dispersion have the drastic impact on the density distribution and smooth the high density layers discovered previously.

Context. Radio relics in galaxy clusters are known to be good laboratories for verification of the applicability of the diffusive shock acceleration (DSA) model in its canonical version. The need for such verification stems from the inconsistencies in the shock properties resulting from radio observations compared to X-ray observations. Aims. In this article we aim to explore how the presence of a second shock in the vicinity of a relic modifies the spectrum of accelerated electrons and decipher which of the involved parameters can have a significant impact on their shape. Methods. We analytically studied DSA of cosmic rays in two stationary shocks aiming to investigate the change of the distribution function. The latter eventually leads to spectrum slope deviations visible in different observations and simulations that do not appear to be explained by the case wherein cosmic rays interact with a single shock wave. Results. We obtain a complex distribution function $f(x,p)$ depending on many parameters (distance between two shocks, compression ratios, spatial diffusion coefficients, injection value, etc.). This function reveals modifications that occur because of the coupled acceleration in both shocks. Apparently, deviations in the particle spectrum from the pure power law depend on a few parameters such as $Q_{1}/Q_{2} $, $\kappa_{1}/\kappa_{2} $, $r_{1}/r_{2} $, and $L$. Although we do not verify this idea by taking a particular cluster as an example, we demonstrate a potential cause of spectral disturbances in radio relics. In general terms, our findings appear to correlate with results from the literature when the distance between the shocks is of the order of the width of a radio relic and $\kappa_{1}/\kappa_{2} \propto 3$.

We investigate the sources of parity asymmetry in the CMB temperature maps using a pixel domain approach. We demonstrate that this anomaly is mainly associated with the presence of two pairs of high asymmetry regions. The first pair of peaks with Galactic coordinates $(l, b) = (212^\circ, -21^\circ)$ and $(32^\circ, 21^\circ)$ is associated with the Northern Galactic Spur and the direction of the dipole modulation of the power spectrum of the CMB anisotropy. The other pair ($(l, b)=(332^\circ, -8^\circ)$ and $(152^\circ, 8^\circ)$) is located within the Galactic plane (the Galactic Cold Spot and its antipodal partner). Similar asymmetric peaks, but with smaller amplitudes, belong to the WMAP/Planck Cold Spot and its partner in the Northern Galactic Spur. These local anomalies increase the odd-multipole power to a level consistent with Gaussian simulations. In contrast, the deficit of symmetric peaks is accompanied by a deficit in the even-multipole power and is the source of the parity asymmetry of the CMB temperature maps at the level of about 3 sigma. We also evaluate the influence of the quadrupole, which is another source of the even-multipole deficit. If the low quadrupole is an intrinsic feature of the theoretical model, it will reduce the significance of the parity asymmetry to around the 2 sigma level. We also investigate the relationship between the asymmetry of the power spectrum and the level of the parity asymmetry in the framework of a model with dipole modulation of a statistically uniform Gaussian signal. We show that these two anomalies are innately linked to each other.

Supernova remnants (SNRs) are observable for about 6-15x10^4 years before they fade into the Galactic interstellar medium. With a Galactic supernova rate of approximately two per century, we can expect to have of the order of 1200 SNRs in our Galaxy. However, only about 300 of them are known to date, with the majority having been discovered in Galactic plane radio surveys. Given that these SNRs represent the brightest tail of the distribution and are mostly located close to the plane, they are not representative of the complete sample. Here we report findings from the search for new SNRs in the eROSITA all-sky survey data which led to the detection of one of the largest SNRs discovered at wavelengths other than the radio: G249.5+24.5. This source is located at a relatively high Galactic latitude, where SNRs are not usually expected to be found. The remnant, 'Hoinga', has a diameter of about 4.4 degrees and shows a circular shaped morphology with diffuse X-ray emission filling almost the entire remnant. Spectral analysis of the remnant emission reveals that an APEC spectrum from collisionally ionised diffuse gas and a plane-parallel shock plasma model with non-equilibrium ionisation are both able to provide an adequate description of the data, suggesting a gas temperature of the order of kT = 0.1 keV and an absorbing column density of N_H=3.6 x 10^20 cm^-2. Subsequent searches for a radio counterpart of the Hoinga remnant identified its radio emission in archival data from the Continuum HI Parkes All-Sky Survey (CHIPASS) and the 408-MHz `Haslam' all-sky survey. The radio spectral index alpha=-0.69 +- 0.08 obtained from these data definitely confirms the SNR nature of Hoinga. From its size and X-ray and radio spectral properties we conclude that Hoinga is a middle-aged Vela-like SNR located at a distance of about twice that of the Vela SNR, i.e. at ~500 pc.

We introduce a family of equations of state (EoS) for hybrid neutron star (NS) matter that is obtained by a two-zone parabolic interpolation between a soft hadronic EoS at low densities and a set of stiff quark matter EoS at high densities within a finite region of chemical potentials $\mu_H < \mu < \mu_Q$. Fixing the hadronic EoS as the APR one and chosing the color-superconductiong, nonlocal NJL model with two free parameters for the quark phase, we perform Bayesian analyses with this two-parameter family of hybrid EoS. Using three different sets of observational constraints that include the mass of PSR J0740+6620, the tidal deformability for GW170817 and the mass-radius relation for PSR J0030+0451 from NICER as obligatory (set 1), while set 2 uses the possible upper limit on the maximum mass from GW170817 as additional constraint and set 3 instead the possibility that the lighter object in the asymmetric binary merger GW190814 is a neutron star. We confirm that in any case the quark matter phase has to be color superconducting with the dimensionless diquark coupling approximately fulfilling the Fierz relation $\eta_D=0.75$ and the most probable solutions exhibiting a proportionality between $\eta_D$ and $\eta_V$, the coupling of the repulsive vector interaction that is required for a sufficiently large maximum mass. We anticipate the outcome of the NICER radius measurement on PSR J0740+6220 as a fictitious constraint and find evidence for claiming that GW190814 was a binary black hole merger if the radius will be 11 km or less.

The spatial-dependent propagation (SDP) model with a nearby source works well to reproduce the co-evolving features of both cosmic ray (CR) nuclei spectra and anisotropy. However, it is well known that the Sun is offset from the Galactic plane. This will lead to a dominating anisotropy in perpendicular direction, which is discrepant with observations. Thus it is necessary to reboot further investigation into the effect of the solar offset. In this work, for the first time the combined studies of the solar offset, nuclei spectra and anisotropy are performed based on the SDP model. As a result, to reproduce CR spectra and anisotropy, the thickness of the halo $\rm (\xi z_h)$ increases linearly with the displacement of the Sun. When the offset is $\rm \sim8~pc$ as estimated from the matter-borne methods, $\rm \xi z_h$ is about 0.9 kpc, which is a much thicker halo than usually. Moreover, the PeV anisotropy could estimate the value of diffusion coefficient, thus breaking the degeneracy of diffusion coefficient and halo thickness. Therefore it is a good messenger to constrain the halo thickness. On the other hand, the anisotropy in PeV energy region, as a new probe, might also shed new light to constrain the solar offset. We hope that the anisotropy at the energies of $\rm \sim TeV$ to $\rm PeV$ can be finely measured by LHAASO experiment, leading to a better understanding about the thick halo.

We present the discovery and follow-up observations of two CCSNe that occurred in the luminous infrared galaxy (LIRG), NGC3256. The first, SN2018ec, was discovered using the ESO HAWK-I/GRAAL adaptive optics seeing enhancer, and was classified as a Type Ic with a host galaxy extinction of $A_V=2.1^{+0.3}_{-0.1}$ mag. The second, AT2018cux, was discovered during the course of follow-up observations of SN2018ec, and is consistent with a sub-luminous Type IIP classification with an $A_V=2.1 \pm 0.4$ mag of host extinction. A third CCSN, PSNJ10275082-4354034 in NGC3256, has previously been reported in 2014, and we recovered the source in late time archival HST imaging. Based on template light-curve fitting, we favour a Type IIn classification for it with modest host galaxy extinction of $A_V=0.3^{+0.4}_{-0.3}$ mag. We also extend our study with follow-up data of the recent Type IIb SN2019lqo and Type Ib SN2020fkb that occurred in the LIRG system Arp299 with host extinctions of $A_V=2.1^{+0.1}_{-0.3}$ and $A_V=0.4^{+0.1}_{-0.2}$ mag, respectively. Motivated by the above, we inspected, for the first time, a sample of 29 CCSNe located within a projected distance of 2.5 kpc from the host galaxy nuclei in a sample of 16 LIRGs. We find that, if star formation within these galaxies is modelled assuming a global starburst episode and normal IMF, there is evidence of a correlation between the starburst age and the CCSN subtype. We infer that the two subgroups of 14 H-poor (Type IIb/Ib/Ic/Ibn) and 15 H-rich (Type II/IIn) CCSNe have different underlying progenitor age distributions, with the H-poor progenitors being younger at 3$\sigma$ significance. However, we do note that the available sample sizes of CCSNe and host LIRGs are so far small, and the statistical comparisons between subgroups do not take into account possible systematic or model errors related to the estimated starburst ages. (abridged)

We study the 2018 Martian Global DustStorm (GDS 2018) over the Southern Polar Region using images obtained by the Visual Monitoring Camera (VMC) on board Mars Express during June and July 2018. Dust penetrated into the polar cap region but never covered the cap completely, and its spatial distribution was nonhomogeneous and rapidly changing. However, we detected long but narrow aerosol curved arcs with a length of 2,000-3,000 km traversing part of the cap and crossing the terminator into the night side. Tracking discrete dust clouds allowed measurements of their motions that were towards the terminator with velocities up to 100 m/s. The images of the dust projected into the Martian limb show maximum altitudes of around 70 km but with large spatial and temporal variations. We discuss these results in the context of the predictions of a numerical model for dust storm scenario.

We report a large difference in neutral hydrogen (H I) and metal column densities between the two sight lines probing opposite sides of the lensing galaxy at $z_\mathrm{lens}$ = 0.83 toward the doubly lensed quasar SBS 0909+532. Using archival HST-STIS and Keck HIRES spectra of the lensed quasar images, we measure log $N_\mathrm{H\;I}$ = 18.77 $\pm$ 0.12 cm$^{-2}$ toward the brighter image ($A$) at an impact parameter of $r_A$ = 3.15 kpc and log $N_\mathrm{H\;I}$ = 20.38 $\pm$ 0.20 cm$^{-2}$ toward the fainter image ($B$) at an impact parameter of $r_B$ = 5.74 kpc. This difference by a factor of $\sim$41 is the highest difference between sight lines for a lens galaxy in which H I has been measured, suggesting patchiness and/or anisotropy on these scales. We estimate an average Fe abundance gradient between the sight lines to be $\geq$ +0.35 dex kpc$^{-1}$. The $N_\mathrm{Fe\;II}$/$N_\mathrm{Mg\;II}$ ratios for the individual components detected in the Keck HIRES spectra have supersolar values for all components in sight line $A$ and for 11 out of 18 components in sight line $B$, suggesting that Type Ia supernovae may have contributed to the chemical enrichment of the galaxy's environment. Additionally, these observations provide complementary information to detections of cold gas in early-type galaxies and the tension between these and some models of cloud survival.

Light bridges (LBs) are bright lanes that divide an umbra into multiple parts in some sunspots. Persistent oscillatory bright fronts at a temperature of $\sim$$10^5$ K are commonly observed above LBs in the 1400/1330 \AA~passbands of the Interface Region Imaging Spectrograph (IRIS). Based on IRIS observations, we report small-scale bright blobs from the oscillating bright front above a light bridge. Some of these blobs reveal a clear acceleration, whereas the others do not. The average speed of these blobs projected onto the plane of sky is $71.7\pm14.7$ km s$^{-1}$, with an initial acceleration of $1.9\pm1.3$ km s$^{-2}$. These blobs normally reach a projected distance of 3--7 Mm from their origin sites. From the transition region images we find an average projected area of $0.57\pm0.37$ Mm$^{2}$ for the blobs. The blobs were also detected in multi-passbands of the Solar Dynamics Observatory, but not in the H$\alpha$ images. These blobs are likely to be plasma ejections, and we investigate their kinematics and energetics. Through emission measure analyses, the typical temperature and electron density of these blobs are found to be around $10^{5.47}$ K and $10^{9.7}$ cm$^{-3}$, respectively. The estimated kinetic and thermal energies are on the order of $10^{22.8}$ erg and $10^{23.3}$ erg, respectively. These small-scale blobs appear to show three different types of formation process. They are possibly triggered by induced reconnection or release of enhanced magnetic tension due to interaction of adjacent shocks, local magnetic reconnection between emerging magnetic bipoles on the light bridge and surrounding unipolar umbral fields, and plasma acceleration or instability caused by upward shocks, respectively.

Some luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) host extremely compact and dusty nuclei. The intense infrared radiation arising from warm dust in these sources is prone to excite vibrational levels of molecules such as HCN. This results in emission from the rotational transitions of vibrationally excited HCN (HCN-vib), with the brightest emission found in compact obscured nuclei (CONs). We aim to establish how common CONs are in the local Universe, and whether their prevalence depends on the luminosity or other properties of the host galaxy. We have conducted an Atacama Large Millimeter/submillimeter Array (ALMA) survey of the rotational J=3-2 transition of HCN-vib in a sample of 46 far-infrared luminous galaxies. Compact obscured nuclei are identified in 38 percent of ULIRGs, 21 percent of LIRGs, and 0 percent of lower luminosity galaxies. We find no dependence on the inclination of the host galaxy, but strong evidence of lower IRAS 25 to 60 {\mu}m flux density ratios (f25/f60) in CONs compared to the rest of the sample. Furthermore, we find that CONs have stronger silicate features (s9.7{\mu}m) but similar PAH equivalent widths (EQW6.2{\mu}m) compared to other galaxies. In the local Universe, CONs are primarily found in (U)LIRGs. High resolution continuum observations of the individual nuclei are required to determine if the CON phenomenon is related to the inclinations of the nuclear disks. The lower f25/f60 ratios in CONs as well as the results for the mid-infrared diagnostics investigated are consistent with large dust columns shifting the nuclear radiation to longer wavelengths, making the mid- and far-infrared "photospheres" significantly cooler than the interior regions. To assess the importance of CONs in the context of galaxy evolution, it is necessary to extend this study to higher redshifts where (U)LIRGs are more common.

Interferometric measurements with arrays of radio antennas are a powerful and widely used technique in astronomy. Recently, this technique has been revisited for the reconstruction of extensive air showers [1]. This radio-interferometric technique exploits the coherence in the radio emission emitted by billions of secondary shower particles to reconstruct the shower parameters, in particular the shower axis and depth of the shower maximum $X_\mathrm{max}$. The accuracy previously demonstrated on simulations with an idealized detector is very promising. In this article we evaluate the potential of interferometric $X_\mathrm{max}$ measurements using (simulated) inclined air showers with sparse antenna arrays under realistic conditions. To determine prerequisites for the application of the radio-interferometric technique with various antenna arrays, the influence of inaccuracies in the time synchronisation between antennas and its inter-dependency with the antenna density is investigated in detail. We find a strong correlation between the antenna multiplicity (per event) and the maximum acceptable time jitter, i.e., inaccuracy in the time synchronisation. For data recorded with a time synchronisation accurate to within 1 ns in the commonly used frequency band of 30 to 80 MHz, an antenna multiplicity of $> 50$ is needed to achieve an $X_\mathrm{max}$ resolution of $\sigma_{X_\mathrm{max}} \lesssim 20$ g cm$^{-2}$. For data recorded with 2 ns accuracy, already $\gtrsim 200$ antennas are needed to achieve this $X_\mathrm{max}$ resolution. Furthermore, we find no advantage reconstructing $X_\mathrm{max}$ from data simulated at higher observation frequencies, i.e., up to several hundred MHz. Finally, we provide a generalisation of our results from very inclined air showers to vertical geometries.

The Cadmium Zinc Telluride (CZT) Imager on board AstroSat is a hard X-ray imaging spectrometer operating in the energy range of 20 $-$ 100 keV. It also acts as an open hard X-ray monitor above 100 keV capable of detecting transient events like the Gamma-ray Bursts (GRBs). Additionally, the instrument has the sensitivity to measure hard X-ray polarization in the energy range of 100 $-$ 400 keV for bright on-axis sources like Crab and Cygnus X-1 and bright GRBs. As hard X-ray instruments like CZTI are sensitive to cosmic rays in addition to X-rays, it is required to identify and remove particle induced or other noise events and select events for scientific analysis of the data. The present CZTI data analysis pipeline includes algorithms for such event selection, but they have certain limitations. They were primarily designed for the analysis of data from persistent X-ray sources where the source flux is much less than the background and thus are not best suited for sources like GRBs. Here, we re-examine the characteristics of noise events in CZTI and present a generalized event selection method that caters to the analysis of data for all types of sources. The efficacy of the new method is reviewed by examining the Poissonian behavior of the selected events and the signal to noise ratio for GRBs.

The Cadmium Zinc Telluride Imager (CZTI) onboard AstroSat is designed for hard X-ray imaging and spectroscopy in the energy range of 20 - 100 keV. The CZT detectors are of 5 mm thickness and hence have good efficiency for Compton interactions beyond 100 keV. The polarisation analysis using CZTI relies on such Compton events and have been verified experimentally. The same Compton events can also be used to extend the spectroscopy up to 380 keV. Further, it has been observed that about 20% pixels of the CZTI detector plane have low gain, and they are excluded from the primary spectroscopy. If these pixels are included, then the spectroscopic capability of CZTI can be extended up to 500 keV and further up to 700 keV with a better gain calibration in the future. Here we explore the possibility of using the Compton events as well as the low gain pixels to extend the spectroscopic energy range of CZTI for ON-axis bright X-ray sources. We demonstrate this technique using Crab observations and explore its sensitivity.

Cadmium Zinc Telluride Imager (CZTI) onboard AstroSat has been a prolific Gamma-Ray Burst (GRB) monitor. While the 2-pixel Compton scattered events (100 - 300 keV) are used to extract sensitive spectroscopic information, the inclusion of the low-gain pixels (around 20% of the detector plane) after careful calibration extends the energy range of Compton energy spectra to 600 keV. The new feature also allows single-pixel spectroscopy of the GRBs to the sub-MeV range which is otherwise limited to 150 keV. We also introduced a new noise rejection algorithm in the analysis ('Compton noise'). These new additions not only enhances the spectroscopic sensitivity of CZTI, but the sub-MeV spectroscopy will also allow proper characterization of the GRBs not detected by Fermi. This article describes the methodology of single, Compton event and veto spectroscopy in 100 - 600 keV for the GRBs detected in the first year of operation. CZTI in last five years has detected around 20 bright GRBs. The new methodologies, when applied on the spectral analysis for this large sample of GRBs, has the potential to improve the results significantly and help in better understanding the prompt emission mechanism.

We aim at identifying very low-mass isolated planetary-mass member candidates in the nearest OB association to the Sun, Upper Scorpius (145 pc; 5-10 Myr), to constrain the form and shape of the luminosity function and mass spectrum in this regime. We conducted a deep multi-band ($Y$=21.2, $J$=20.5, $Z$=22.0 mag) photometric survey of six square degrees in the central region of Upper Scorpius. We extend the current sequence of astrometric and spectroscopic members by about two magnitudes in $Y$ and one magnitude in $J$, reaching potentially T-type free-floating members in the association with predicted masses below 5 Jupiter masses, well into the planetary-mass regime. We extracted a sample of 57 candidates in this area and present infrared spectroscopy confirming two of them as young L-type members with characteristic spectral features of 10 Myr-old brown dwarfs. Among the 57 candidates, we highlight 10 new candidates fainter than the coolest members previously confirmed spectroscopically. We do not see any obvious sign of decrease in the mass spectrum of the association, suggesting that star processes can form substellar objects with masses down to 4-5 Jupiter masses.

We present measurements of the morphology of star-forming gas in galaxies from the EAGLE simulations, and its alignment relative to stars and dark matter (DM). Imaging of such gas in the radio continuum enables weak lensing experiments that complement traditional optical approaches. Star-forming gas is typically more flattened than its associated stars and DM, particularly for present-day subhaloes of total mass $\sim$$10^{ 12-12.5} \mathrm{M_{ \odot}}$, which preferentially host star-forming galaxies with rotationally-supported stellar discs. Such systems have oblate, spheroidal star-forming gas distributions, but in both less- and more-massive subhaloes the distributions tend to be prolate, and its morphology correlates positively and significantly with that of its host galaxy's stars, both in terms of sphericity and triaxiality. The minor axis of star-forming gas most commonly aligns with the minor axis of its host subhalo's DM, but often aligns more closely with one of the other two principal axes of the DM distribution in prolate subhaloes. Star-forming gas aligns with DM less strongly than is the case for stars, but its morphological minor axis aligns closely with its kinematic axis, affording a route to observational identification of the unsheared morphological axis. The projected ellipticities of star-forming gas in EAGLE are consistent with shapes inferred from high-fidelity radio continuum images, and they exhibit greater shape noise than is the case for images of the stars, owing to the greater characteristic flattening of star-forming gas with respect to stars.

We reconstruct the neutrino mass as a function of redshift, $z$, from current cosmological data using both standard binned priors and linear spline priors with variable knots. Using cosmic microwave background temperature, polarization and lensing data, in combination with distance measurements from baryonic acoustic oscillations and supernovae, we find that the neutrino mass is consistent with $\sum m_\nu(z)=$ const. We obtain a larger bound on the neutrino mass at low redshifts coinciding with the onset of dark energy domination, $\sum m_\nu(z=0)<1.41$ eV (95% CL). This result can be explained either by the well-known degeneracy between $\sum m_\nu$ and $\Omega_\Lambda$ at low redshifts, or by models in which neutrino masses are generated very late in the Universe. We convert our results into cosmological limits for models with post-recombination neutrino decay and find $\sum m_\nu <0.19$ eV (95% CL), which is below the sensitivity of the KATRIN experiment. Thus, a neutrino mass discovery by KATRIN would hint towards models predicting both post-recombination neutrino mass generation and subsequent relic neutrino annihilation.

This paper investigates the feasibility of simulating Fuzzy Dark Matter (FDM) with a reduced number of spatial dimensions. Our aim is to set up a realistic, yet numerically inexpensive, toy model in $(1+1)$-dimensional space time, that - under well controlled system conditions - is capable of realizing important aspects of the full-fledged $(3+1)$-FDM phenomenology by means of one-dimensional analogues. Based on the coupled, nonlinear and nonlocal $(3+1)$-Sch\"odinger-Poisson equation under periodic boundary conditions, we derive two distinct one-dimensional models that differ in their transversal matter distribution and consequently in their nonlocal interaction along the single dimension of interest. We show that these discrepancies change the relaxation process of initial states as well as the asymptotic, i.e., thermalized and virialized, equilibrium state. Our investigation includes the dynamical evolution of artificial initial conditions for non-expanding space, as well as cosmological initial conditions in expanding space. The findings of this work are relevant for the interpretation of numerical simulation data modelling nonrelativistic fuzzy cold dark matter in reduced dimensions, in the quest for testing such models and for possible laboratory implementations of them.

Aims. Planetary nebulae (PNe) are a brief phase of stellar evolution and as such some of the rarest objects in our galaxy. Accurate identification of PN central stars (CSPNe) in large surveys such as Gaia is key to exploiting the scientific potential of those surveys for the study of PNe and stellar evolution. Methods. We apply our automated search method to identify source detections corresponding to CSPNe and compact PNe in the recently released Gaia Early Data Release 3 (EDR3). The method is updated to incorporate photometric uncertainties. Results. The new catalogue offers improved completeness and accuracy over the previous one, and will continue to be valid for the upcoming Gaia Data Release 3 (DR3). The identification of EDR3 sources unlocks the improved astrometry and photometry stemming from calibration updates and a longer mission duration, even for sources that have previously been identified.

We propose a new model of Early Dark Energy (EDE) as a solution to the Hubble tension in cosmology, the apparent discrepancy between local measurements of the Hubble constant $H_0\simeq 74$ km s$^{-1}$ Mpc$^{-1}$ and $H_0\simeq 67$ km s$^{-1}$ Mpc$^{-1}$ inferred from the Cosmic Microwave Background (CMB). In Chain EDE, the Universe undergoes a series of first order phase transitions, starting at a high energy vacuum in a potential, and tunneling down through a chain of every lower energy metastable minima. As in all EDE models, the contribution of the vacuum energy to the total energy density of the universe is initially negligible, but reaches $\sim 10\%$ around matter-radiation equality, before cosmological data require it to redshift away quickly -- at least as fast as radiation. We indeed obtain this required behavior with a series of $N$ tunneling events, and show that for $N>20,000$ the phase transitions are rapid enough to allow fast percolation and thereby avoid large scale anisotropies in the CMB. We construct a specific example of Chain EDE featuring a scalar field in a quasiperiodic potential (a tilted cosine), which is ubiquitous in axion physics and, therefore, carries strong theoretical motivation. Interestingly, the energy difference between vacua can be roughly the size of today's Dark Energy (meV scale). Therefore, the end result of Chain EDE could provide a natural explanation of Dark Energy, if the tunneling becomes extremely slow in the final step before the field reaches zero (or negative) energy. We discuss a simple mechanism which can stop the scalar field in the desired minimum. Thus Chain EDE offers the exciting prospect to explain EDE and Dark Energy by the same scalar field.

Geostationary observations offer the unique opportunity to resolve the diurnal cycle of the Earth's Radiation Budget at the top of the atmosphere (TOA), crucial for climate-change studies. However, a drawback of the continuous temporal coverage of the geostationary orbit is the fixed viewing geometry. As a consequence, errors in the angular distribution models (ADMs) used in the radiance-to-flux conversion process can result in systematic errors of the estimated radiative fluxes. In this work, focusing on clear-sky reflected TOA observations, we compare the overlapping views from Meteosat Second Generation satellites at 0{\deg} and 41.5{\deg}E longitude which enable a quantification of viewing-angle-dependent differences. Using data derived from the Spinning Enhanced Visible and InfraRed Imager (SEVIRI), we identify some of the main sources of discrepancies, and show that they can be significantly reduced at the level of one month. This is achieved, separately for each satellite, via a masking procedure followed by an empirical fit at the pixel-level that takes into account all the clear-sky data from that satellite, calculated separately per timeslot of the day, over the month of November 2016. The method is then applied to each month of 2017, and gives a quadratic mean of the albedo root-mean squared difference over the dual-view region which is comparable from month to month, with a 2017 average value of 0.01. Sources of discrepancies include the difficulty to estimate the flux over the sunglint ocean region close to the limbs, the absence of dedicated angular distribution models for the aerosol-over-ocean case in the data processing, and the existence of an observer-dependent diurnal-asymmetry artefact affecting the clear-sky-albedo dependence on the solar zenith angle particularly over land areas.

We consider a situation where a light source orbiting the innermost stable circular orbit (ISCO) of the Kerr black hole is gently falling from the marginally stable orbit due to an infinitesimal perturbation. Assuming that the light source emits photons isotropically, we show that the last radius at which more than 50\% of emitted photons can escape to infinity is approximately halfway between the ISCO radius and the event horizon radius. To evaluate them, we determine emitter orbits from the vicinity of the ISCO, which are uniquely specified for each black hole spin, and identify the conditions for a photon to escape from any point on the equatorial plane of the Kerr spacetime to infinity by specifying regions in the two-dimensional photon impact parameter space completely. We further show that the proper motion of the emitter affects the photon escape probability and blueshifts the energy of emitted photons.

This paper presents an exact solution of the Einstein-Maxwell field equations in a static and spherically symmetric Schwarzschild canonical coordinate system in the presence of charged perfect fluid. We have employed the Vaidya-Tikekar ansatz for the metric potential. Using graphical analysis and tabular information we have shown that our model obeys all the physical requirements and stability conditions required for a realistic stellar model. This theoretical model approximates observations of pulsar PSR B0943+10 to an excellent degree of accuracy.

The destructive interference of the neighbouring field configurations with infinite classical action in the gravitational path integral approach serves as a dynamical mechanism resolving the black hole singularity problem. It also provides an isotropic and homogeneous early universe without the need of inflation. In this work, we elaborate on the finite action in the framework of Horava-Lifshitz gravity -- a ghost-free QFT. Assuming the mixmaster chaotic solutions in the projectable H-L theory, we show that the beginning of the universe is homogeneous and isotropic. Furthermore, we show that the H-L gravity action selects only the regular black-hole spacetimes. We also comment on possibility of traversable wormholes in theories with higher curvature invariants.

The unprecedented image of the M87* supermassive black hole has sparked some controversy over its usefulness as a test of the general relativistic Kerr metric. The criticism is mainly related to the black hole's quasi-circular shadow and advocates that its radius depends not only on the black hole's true spacetime properties but also on the poorly known physics of the illuminating accretion flow. In this paper we take a sober view of the problem and argue that our ability to probe gravity with a black hole shadow is only partially impaired by the matter degrees of freedom and the number of non-Kerr parameters used in the model. As we show here, a much more fatal complication arises from the mass scaling of the dimensional coupling constants that typically appear in non-GR theories of gravity. Existing limits from gravitational wave observations imply that supermassive systems like the M87* black hole would suffer a suppression of all non-GR deviation parameters in their metric, making the spacetime and the produced shadow virtually Kerr. Therefore, a supermassive black hole shadow is likely to probe only those extensions of General Relativity which are endowed with dimensionless coupling constants.