Papers for Thursday, May 20 2021

We investigate the kinematic properties of a large (N=998) sample of COSMOS spectroscopic galaxy members distributed among 79 groups. We identify the Brightest Group Galaxies (BGGs) and cross-match our data with the VLA-COSMOS Deep survey at 1.4 GHz, classifying our parent sample into radio/non-radio BGGs and radio/non-radio satellites. The radio luminosity distribution spans from $L_R\sim2\times10^{21}$ W Hz$^{-1}$ to $L_R\sim3\times$10$^{25}$ W Hz$^{-1}$. A phase-space analysis, performed by comparing the velocity ratio (line-of-sight velocity divided by the group velocity dispersion) with the galaxy-group centre offset, reveals that BGGs (radio and non-radio) are mostly ($\sim$80\%) ancient infallers. Furthermore, the strongest ($L_R>10^{23}$ W Hz$^{-1}$) radio galaxies are always found within 0.2$R_{\rm vir}$ from the group centre. Comparing our samples with HORIZON-AGN, we find that the velocities and offsets of simulated galaxies are more similar to radio BGGs than to non-radio BGGs, albeit statistical tests still highlight significant differences between simulated and real objects. We find that radio BGGs are more likely to be hosted in high-mass groups. Finally, we observe correlations between the powers of BGG radio galaxies and the X-ray temperatures, $T_{\rm x}$, and X-ray luminosities, $L_{\rm x}$, of the host groups. This supports the existence of a link between the intragroup medium and the central radio source. The occurrence of powerful radio galaxies at group centres can be explained by Chaotic Cold Accretion, as the AGN can feed from both the galactic and intragroup condensation, leading to the observed positive $L_{\rm R}-T_{\rm x}$ correlation.

Most stars in the Galaxy, including the Sun, were born in high-mass star-forming regions. It is hence important to study the chemical processes in these regions to better understand the chemical heritage of both the Solar System and most stellar systems in the Galaxy. The molecular ion HCNH+ is thought to be a crucial species in ion-neutral astrochemical reactions, but so far it has been detected only in a handful of star-forming regions, and hence its chemistry is poorly known. We have observed with the IRAM-30m Telescope 26 high-mass star-forming cores in different evolutionary stages in the J=3-2 rotational transition of HCNH+. We report the detection of HCNH+ in 16 out of 26 targets. This represents the largest sample of sources detected in this molecular ion so far. The fractional abundances of HCNH+, [HCNH+], w.r.t. H2, are in the range 0.9 - 14 X $10^{-11}$, and the highest values are found towards cold starless cores. The abundance ratios [HCNH+]/[HCN] and [HCNH+]/[HCO+] are both < 0.01 for all objects except for four starless cores, for which they are well above this threshold. These sources have the lowest gas temperature in the sample. We run two chemical models, a "cold" one and a "warm" one, which attempt to match as much as possible the average physical properties of the cold(er) starless cores and of the warm(er) targets. The reactions occurring in the latter case are investigated in this work for the first time. Our predictions indicate that in the warm model HCNH+ is mainly produced by reactions with HCN and HCO+, while in the cold one the main progenitor species of HCNH+ are HCN+ and HNC+. The results indicate that the chemistry of HCNH+ is different in cold/early and warm/evolved cores, and the abundance ratios [HCNH+]/[HCN] and [HCNH+]/[HCO+] is a useful astrochemical tool to discriminate between different evolutionary phases in the process of star formation.

The Primordial Inflation Polarization Explorer (PIPER) is a stratospheric balloon payload to measure polarization of the cosmic microwave background. Twin telescopes mounted within an open-aperture bucket dewar couple the sky to bolometric detector arrays. We reduce detector loading and photon noise by cooling the entire optical chain to 1.7 K or colder. A set of fountain-effect pumps sprays superfluid liquid helium onto each optical surface, producing helium flows of 50--100 cm^3 / s at heights up to 200 cm above the liquid level. We describe the fountain-effect pumps and the cryogenic performance of the PIPER payload during two flights in 2017 and 2019.

Motivated by recent theoretical work on tidal disruption events and other peculiar transients, we present moving-mesh radiation-hydrodynamic simulations of radiative luminosity emitted by a central source being reprocessed by a wind-like outflow. We couple the moving-mesh hydrodynamic code JET with our newly-developed radiation module based on mixed-frame grey flux-limited diffusion with implicit timestep update. This allows us to study the self-consistent multi-dimensional radiation-hydrodynamic evolution over more than ten orders of magnitude in both space and time in a single run. We simulate an optically-thick spherical wind with constant or evolving mass-loss rate, which is irradiated by a central isotropic or angularly-dependent radiation source. Our spherically-symmetric simulations confirm previous analytic results by identifying different stages of radiation reprocessing: radiation trapped in the wind, diffusing out through the wind, and reaching constant maximum attenuation. We find that confining the central radiation source in a cone with moderate opening angles decreases significantly the early flux along sightlines oriented away from the direction of radiation injection but that the reprocessed radiation becomes isotropic roughly after one lateral diffusion time through the ejecta. We discuss further applications and guidelines for the use of our novel radiation-hydrodynamics tool in the context of transient modelling.

The first generation (Pop-III) stars can ionize 1-10% of the universe by $z=15$, when the metal-enriched (Pop-II) stars may contribute negligibly to the ionization. This low ionization tail might leave detectable imprints on the large-scale CMB E-mode polarization. However, we show that physical models for reionization are unlikely to be sufficiently extended to detect any parameter beyond the total optical depth through reionization. This result is driven in part by the total optical depth inferred by Planck 2018, indicating a reionization midpoint around $z=8$, which in combination with the requirement that reionization completes by $z\approx 5.5$ limits the amplitude of an extended tail. To demonstrate this, we perform semi-analytic calculations of reionization including Pop-III star formation in minihalos with Lyman-Werner feedback. We find that standard Pop-III models need to produce very extended reionization at $z>15$ to be distinguishable at 2-$\sigma$ from Pop-II-only models, assuming a cosmic variance-limited measurement of the low-$\ell$ EE power spectrum. However, we show that unless appealing to extreme Pop-III scenarios, structure formation makes it quite challenging to produce high enough Thomson scattering optical depth from $z>15$, $\tau(z>15)$, and still be consistent with other observational constraints on reionization. We also find that the Planck likelihood is mostly sensitive to the total $\tau$ and has minor information content on $\tau(z>15)$. We comment on the tension with the Planck 2018 $\tau(z>15)$ constraints whose results suggest more information in $\tau(z>15)$ than we find.

Our understanding of the formation and evolution of binary black holes (BBHs) is significantly impacted by the recent discoveries made by the LIGO/Virgo collaboration. Of utmost importance is the detection of the most massive BBH system, GW190521. Here we investigate what it takes for field massive stellar binaries to account for the formation of such massive BBHs. Whether the high mass end of the BH mass function is populated by remnants of massive stars that either formed at extremely low metallicities and avoid the pair-instability mass gap or increase their birth mass beyond the pair-instability mass gap through the accretion of gas from the surrounding medium. We show that assuming that massive stars at very low metallicities can form massive BHs by avoiding pair-instability supernova, coupled with a correspondingly high formation efficiency for BBHs, can explain the observed BH mass function. To this end, one requires a relation between the initial and final mass of the progenitor stars at low metallicities that is shallower than what is expected from wind mass loss alone. On the other hand, assuming pair-instability operates at all metallicities, one can account for the observed BH mass function if at least about 10% of the BHs born at very low metallicities double their mass before they merge because of accretion of ambient gas. Such BBHs will have to spend about a Gyr within a parsec length-scale of their parent atomic cooling halos or a shorter timescale if they reside in the inner sub-parsecs of their host dark matter halos. Future stellar evolution calculations of massive stars at very low metallicity and hydrodynamical simulations of gas accretion onto BBHs born in atomic cooling halos can shed light on this debate.

An empirical model is presented that links, for the first time, the demographics of AGN to their ensemble X-ray variability properties. Observations on the incidence of AGN in galaxies are combined with (i) models of the Power Spectrum Density (PSD) of the flux variations of AGN and (ii) parameterisations of the black-hole mass vs stellar-mass scaling relation, to predict the mean excess variance of active black-hole populations in cosmological volumes. We show that the comparison of the model with observational measurements of the ensemble excess variance as a function of X-ray luminosity provides a handle on both the PSD models and the black-hole mass vs stellar mass relation. We find strong evidence against a PSD model that is described by a broken power-law and a constant overall normalisation. Instead our analysis indicates that the amplitude of the PSD depends on the physical properties of the accretion events, such as the Eddington ratio and/or the black hole mass. We also find that current observational measurements of the ensemble excess variance are consistent with the black-hole mass vs stellar mass relation of local spheroids based on dynamically determined black-hole masses. We also discuss future prospects of the proposed approach to jointly constrain the PSD of AGN and the black-hole mass vs stellar mass relation as a function of redshift.

A recent paper by Siraj & Loeb (2021) entitled "Breakup of a long-period comet as the origin of the dinosaur extinction" attempts to revive the perennial debate about what type of body hit the Earth 66 million years ago, triggering the end-Cretaceous extinction. Here we critique the paper and assess the evidence it presents. To consider a comet more likely than an asteroid requires extreme assumptions about how comets fragment, conflation of carbonaceous chondrites with specific types of carbonaceous chondrites, and a blind eye to the evidence of the iridium layer.

The virtual meeting was a success. Several people told us that this was "the best virtual meeting they had seen so far", which, a year into the pandemic and without a commercial provider in the back, is a great success. The biggest point of criticism was the timing: We had programming from UTC 17:00-22:00 (evening and night in central Europe, afternoon on the US East Coast, during the day in South America and on the US West coast, but in the middle of the night for Asia and Australia). There is no good solution, but at least some variation in session time might go a long way to make it easier for all to attend at least some sessions. Feedback also indicates that the schedule was too compressed. Poster sessions and social contacts with the tool Gathertown worked out really well for all that used it. Our way of combining several services (Zoom for plenary and break-out rooms, Zenodo for uploading and viewing posters and proceedings, Google forms for registration and abstract submission, gathertown) allowed for a very low-cost meeting with little overhead (total cost: 600 $ for gathertown, zoom was provided through an institutional subscription, just 4 people on the LOC).

Small scale transients occur in the Solar corona at much higher frequencies than flares and play a significant role in coronal dynamics. Here we study three well-identified transients discovered by Hi-C and also detected by the EUV channels of Atmospheric Imaging Assembly (AIA) on board Solar Dynamics Observatory (SDO). We use 0-D enthalpy-based hydrodynamical simulations and produce synthetic light curves to compare with AIA observations. We have modeled these transients as loops of ~ 1.0~Mm length depositing energies ~ 10^23 ergs in ~ 50 seconds. The simulated synthetic light curves show reasonable agreement with the observed light curves. During the initial phase, conduction flux from the corona dominates over the radiation, like impulsive flaring events. Our results further show that the time-integrated net enthalpy flux is positive, hence into the corona. The fact that we can model the observed light curves of these transients reasonably well by using the same physics as those for nanoflares, microflares, and large flares, suggests that these transients may have a common origin.

We present 42 rapidly evolving (time spent above half-maximum brightness $t_{1/2}<12$d) extragalactic transients from Phase I of the Zwicky Transient Facility (ZTF), of which 22 have spectroscopic classifications. This is one of the largest systematically selected samples of day-timescale transients, and the first with spectroscopic classifications. Most can be classified as core-collapse supernovae (SNe), and we identify several predominant subtypes: (1) subluminous Type IIb or Type Ib SNe; (2) luminous Type Ibn or hybrid IIn/Ibn SNe; and (3) radio-loud, short-duration luminous events similar to AT2018cow. We conclude that rates quoted in the literature for rapidly evolving extragalactic transients are dominated by the subluminous events (mostly Type IIb SNe). From our spectroscopic classifications and radio, X-ray, and millimeter-band upper limits, we are motivated to consider the AT2018cow-like objects a distinct class, and use ZTF's systematic classification experiments to calculate that their rate does not exceed 0.1% of the local core-collapse SN rate, in agreement with previous work. By contrast, most other events are simply the extreme of a continuum of established SN types extending to ordinary timescales. The light curves of our objects are very similar to those of unclassified events in the literature, illustrating how spectroscopically classified samples of low-redshift objects in shallow surveys like ZTF can be used to photometrically classify larger numbers of events at higher redshift.

We present a detailed analysis for a subset of the high resolution (~35 mas, or 5 au) ALMA observations from the Disk Substructures at High Angular Resolution Project (DSHARP) to search for faint 1.3 mm continuum emission associated with dusty circumplanetary material located within the narrow annuli of depleted emission (gaps) in circumstellar disks. This search used the Jennings et al. (2020) $\tt{frank}$ modeling methodology to mitigate contamination from the local disk emission, and then deployed a suite of injection-recovery experiments to statistically characterize point-like circumplanetary disks in residual images. While there are a few putative candidates in this sample, they have only marginal local signal-to-noise ratios and would require deeper measurements to confirm. Associating a 50% recovery fraction with an upper limit, we find these data are sensitive to circumplanetary disks with flux densities $\gtrsim 50-70$ $\mu$Jy in most cases. There are a few examples where those limits are inflated ($\gtrsim 110$ $\mu$Jy) due to lingering non-axisymmetric structures in their host circumstellar disks, most notably for a newly identified faint spiral in the HD 143006 disk. For standard assumptions, this analysis suggests that these data should be sensitive to circumplanetary disks with dust masses $\gtrsim 0.001-0.2$ M$_\oplus$. While those bounds are comparable to some theoretical expectations for young giant planets, we discuss how plausible system properties (e.g., relatively low host planet masses or the efficient radial drift of solids) could require much deeper observations to achieve robust detections.

The follow-up of external science alerts received from Gamma-Ray Bursts (GRB) and Gravitational Waves (GW) detectors is one of the AGILE Team's current major activities. The AGILE team developed an automated real-time analysis pipeline to analyse AGILE Gamma-Ray Imaging Detector (GRID) data to detect possible counterparts in the energy range 0.1-10 GeV. This work presents a new approach for detecting GRBs using a Convolutional Neural Network (CNN) to classify the AGILE-GRID intensity maps improving the GRBs detection capability over the Li&Ma method, currently used by the AGILE team. The CNN is trained with large simulated datasets of intensity maps. The AGILE complex observing pattern due to the so-called 'spinning mode' is studied to prepare datasets to test and evaluate the CNN. A GRB emission model is defined from the Second Fermi-LAT GRB catalogue and convoluted with the AGILE observing pattern. Different p-value distributions are calculated evaluating with the CNN millions of background-only maps simulated varying the background level. The CNN is then used on real data to analyse the AGILE-GRID data archive, searching for GRB detections using the trigger time and position taken from the Swift-BAT, Fermi-GBM, and Fermi-LAT GRB catalogues. From these catalogues, the CNN detects 21 GRBs with a significance $\geq 3 \sigma$, while the Li&Ma method detects only two GRBs. The results shown in this work demonstrate that the CNN is more effective in detecting GRBs than the Li&Ma method in this context and can be implemented into the AGILE-GRID real-time analysis pipeline.

The Crab Pulsar's radio emission is unusual, consisting predominantly of giant pulses, with durations of about a micro-second but structure down to the nano-second level, and brightness temperatures of up to $10^{37}\,$K. It is unclear how giant pulses are produced, but they likely originate near the pulsar's light cylinder, where corotating plasma approaches the speed of light. Here, we report observations in the 400-800 MHz frequency band, where the pulses are broadened by scattering in the surrounding Crab nebula. We find that some pulse frequency spectra show strong bands, which vary during the scattering tail, in one case showing a smooth upward drift. While the banding may simply reflect interference between nano-second scale pulse components, the variation is surprising, as in the scattering tail the only difference is that the source is observed via slightly longer paths, bent by about an arcsecond in the nebula. The corresponding small change in viewing angle could nevertheless reproduce the observed drift by a change in Doppler shift, if the plasma that emitted the giant pulses moved highly relativistically, with a Lorentz factor $\gamma\sim10^4$. If so, this would support models that appeal to highly relativistic plasma to transform ambient magnetic structures to coherent GHz radio emission, be it for giant pulses or for potentially related sources, such as fast radio bursts.

Magnetic fields provide an important probe of the thermal, material, and structural history of planetary and sub-planetary bodies. Core dynamos are a potential source of magnetic field amplification in differentiated bodies, but evidence of magnetization in undifferentiated bodies requires a different mechanism. Here we study stellar wind-induced magnetization (WIM) of an initially unmagnetized body using analytic theory and numerical simulations, employing the resistive MHD AstroBEAR adaptive mesh refinement (AMR) multiphysics code. We obtain a broadly applicable scaling relation for the peak magnetization achieved once a wind advects, piles-up, and drapes a body with magnetic field, reaching a quasi-steady state. We find that the dayside magnetic field for a sufficiently conductive body saturates when it balances the sum of incoming solar wind ram, magnetic, and thermal pressures. Stronger amplification results from pileup by denser and faster winds. Careful quantification of numerical diffusivity is required for accurately interpreting the peak magnetic field strength from simulations and corroborating with theory. As specifically applied to the Solar System, we find that early solar wind-induced magnetization is a viable explanation for observed paleointensities in some undifferentiated bodies. This magnetism mechanism may also be applicable for other Solar System bodies, including metal-rich bodies to be visited in future space missions such as the asteroid (16) Psyche.

We provide new ionization correction factors (ICFs) for carbon, nitrogen, neon, sulfur, chlorine, and argon in giant H II regions. The ICFs were computed using the most representative photoionization models from a large initial grid. The models were selected using an observational sample of 985 giant H II regions (GHR) in spiral galaxies and blue compact galaxies (BCG). The observational sample was also used to assign a weight to each model describing how well it agrees with observations in the [O III]/Hbeta versus [N II]/Halpha diagram. In addition to the ICFs we provide, for the first time, analytical expressions for their formal uncertainties. We use our ICFs to compute the abundances of nitrogen, neon, sulfur, and argon in our samples. Our abundances are robust within the adopted framework, but may require revision in the case of important changes in atomic data or in the spectral energy distribution of the ionizing radiation in H II regions. Considering the abundance patterns we obtained for the BCG sample (abundances for the GHR sample are less reliable) we find that oxygen is depleted into dust grains at a rate increasing with metallicity and reaching 0.12 dex at solar abundances. The discussion of possible depletion of sulfur and argon requires considering recent Type Ia Supernova yields, which are still uncertain.

We present a detailed investigation of millimeter-wave line emitters ALMA J010748.3-173028 (ALMA-J0107a) and ALMA J010747.0-173010 (ALMA-J0107b), which were serendipitously uncovered in the background of the nearby galaxy VV114 with spectral scan observations at $\lambda$ = 2 - 3 mm. Via Atacama Large Millimeter/submillimeter Array (ALMA) detection of CO(4-3), CO(3-2), and [CI](1-0) lines for both sources, their spectroscopic redshifts are unambiguously determined to be $z= 2.4666\pm0.0002$ and $z=2.3100\pm0.0002$, respectively. We obtain the apparent molecular gas masses $M_{\rm gas}$ of these two line emitters from [CI] line fluxes as $(11.2 \pm 3.1) \times 10^{10} M_\odot$ and $(4.2 \pm 1.2) \times 10^{10} M_\odot$, respectively. The observed CO(4-3) velocity field of ALMA-J0107a exhibits a clear velocity gradient across the CO disk, and we find that ALMA-J0107a is characterized by an inclined rotating disk with a significant turbulence, that is, a deprojected maximum rotation velocity to velocity dispersion ratio $v_{\rm max}/\sigma_{v}$ of $1.3 \pm 0.3$. We find that the dynamical mass of ALMA-J0107a within the CO-emitting disk computed from the derived kinetic parameters, $(1.1 \pm 0.2) \times 10^{10}\ M_\odot$, is an order of magnitude smaller than the molecular gas mass derived from dust continuum emission, $(3.2\pm1.6)\times10^{11}\ M_{\odot}$. We suggest this source is magnified by a gravitational lens with a magnification of $\mu \gtrsim10$, which is consistent with the measured offset from the empirical correlation between CO-line luminosity and width.

Phosphine is now well established as a biosignature, which has risen to prominence with its recent tentative detection on Venus. To follow up this discovery and related future exoplanet biosignature detections, it is important to spectroscopically detect the presence of phosphorus-bearing atmospheric molecules that could be involved in the chemical networks producing, destroying or reacting with phosphine. We start by enumerating phosphorus-bearing molecules (P-molecules) that could potentially be detected spectroscopically in planetary atmospheres and collecting all available spectral data. Gaseous P-molecules are rare, with speciation information scarce. Very few molecules have high accuracy spectral data from experiment or theory; instead, the best available data is from the RASCALL approach and obtained using functional group theory. Here, we present a high-throughput approach utilising established computational quantum chemistry methods (CQC) to produce a database of approximate infrared spectra for 958 P-molecules. These data are of interest for astronomy and astrochemistry (importantly identifying potential ambiguities in molecular assignments), improving RASCALL's underlying data, big data spectral analysis and future machine learning applications. However, this data will probably not be sufficiently accurate for secure experimental detections of specific molecules within complex gaseous mixtures in laboratory or astronomy settings.

The detection of gravitational waves from a neutron star merger, GW170817, marked the dawn of a new era in time-domain astronomy. Monitoring of the radio emission produced by the merger, including high-resolution radio imaging, enabled measurements of merger properties including the energetics and inclination angle. In this work we compare the capabilities of current and future gravitational wave facilities to the sensitivity of radio facilities to quantify the prospects for detecting the radio afterglows of gravitational wave events. We consider three observing strategies to identify future mergers -- widefield follow-up, targeting galaxies within the merger localisation and deep monitoring of known counterparts. We find that while planned radio facilities like the Square Kilometre Array will be capable of detecting mergers at gigaparsec distances, no facilities are sufficiently sensitive to detect mergers at the range of proposed third-generation gravitational wave detectors that would operate starting in the 2030s.

To explain X-ray spectra of active galactic nuclei (AGN), non-thermal activity in AGN coronae such as pair cascade models has been extensively discussed in the past literature. Although X-ray and gamma-ray observations in the 1990s disfavored such pair cascade models, recent millimeter-wave observations of nearby Seyferts establish the existence of weak non-thermal coronal activity. Besides, the IceCube collaboration reported NGC 1068, a nearby Seyfert, as the hottest spot in their 10-yr survey. These pieces of evidence are enough to investigate the non-thermal perspective of AGN coronae in depth again. This article summarizes our current observational understandings of AGN coronae and describes how AGN coronae generate high-energy particles. We also provide ways to test the AGN corona model with radio, X-ray, MeV gamma-ray, and high-energy neutrino observations.

Surveys of galaxy distances and radial peculiar velocities can be used to reconstruct the large scale structure. Other than systematic errors in the zero-point calibration of the galaxy distances the main source of uncertainties of such data are errors on the distance moduli, assumed here to be Gaussian and thus turn into lognormal errors on distances and velocities. Naively treated, it leads to spurious nearby outflow and strong infall at larger distances. The lognormal bias is corrected here and tested against mock data extracted from a $\Lambda$CDM simulation, designed to statistically follow the grouped Cosmicflows-3 (CF3) data. Considering a subsample of data points, all of which have the same true distances or same redshifts, the lognormal bias arises because the means of the distributions of observed distances and velocities are skewed off the means of the true distances and velocities. Yet, the medians are invariant under the lognormal transformation. That invariance allows the Gaussianization of the distances and velocities and the removal of the lognormal bias. This Bias Gaussianization correction (BGc) algorithm is tested against mock CF3 catalogs. The test consists of a comparison of the BGC estimated with the simulated distances and velocities and of an examination of the Wiener filter reconstruction from the BGc data. Indeed, the BGc eliminates the lognormal bias. The estimation of Hubble's ($H_{0}$) constant is also tested. The residual of the BGc estimated $H_{0}$ from the simulated values is $0.6 \pm 0.7 {\rm kms}^{-1}{\rm Mpc}^{-1}$ and is dominated by the cosmic variance. The BGc correction of the actual CF3 data yields $H_{0} = 75.8 \pm 1.1 {\rm kms}^{-1}{\rm Mpc}^{-1}$ .

Ellerman bombs (EBs) and Ultraviolet (UV) bursts are common brightening phenomena which are usually generated in the low solar atmosphere of emerging flux regions. In this paper, we have investigated the emergence of an initial un-twisted magnetic flux rope based on three-dimensional (3D) magneto-hydrodynamic (MHD) simulations. The EB-like and UV burst-like activities successively appear in the U-shaped part of the undulating magnetic fields triggered by Parker Instability. The EB-like activity starts to appear earlier and lasts for about 80 seconds. Six minutes later, a much hotter UV burst-like event starts to appear and lasts for about 60 seconds. Along the direction vertical to the solar surface, both the EB and UV burst start in the low chromosphere, but the UV burst extends to a higher altitude in the up chromosphere. The regions with apparent temperature increase in the EB and UV burst are both located inside the small twisted flux ropes generated in magnetic reconnection processes, which are consistent with the previous 2D simulations that most hot regions are usually located inside the magnetic islands. However, the twisted flux rope corresponding to the EB is only strongly heated after it floats up to an altitude much higher than the reconnection site during that period. Our analyses show that the EB is heated by the shocks driven by the strong horizontal flows at two sides of the U-shaped magnetic fields. The twisted flux rope corresponding to the UV burst is heated by the driven magnetic reconnection process.

Under certain conditions, stellar radial velocities can be determined from astrometry, without any use of spectroscopy. This enables to identify phenomena other than the Doppler effect, that are displacing spectral lines. We aim to use the change of stellar proper motions over time (perspective acceleration) to determine radial velocities from accurate astrometric data, now available from the Gaia and Hipparcos missions. Positions and proper motions at the epoch of Hipparcos are compared with values propagated back from the epoch of the Gaia Early Data Release 3. This propagation depends on the radial velocity, which obtains its value from an optimal fit, assuming space motion to be uniform relative to the solar system barycentre. For 926 nearby stars we obtain astrometric radial velocities with formal uncertainties better than 100 km/s; for 55 stars the uncertainty is below 10 km/s, and for seven it is below 1 km/s. Most stars that are not components of double or multiple systems show good agreement with available spectroscopic radial velocities. Astrometry offers geometric methods to determine stellar radial velocity, irrespective of complexities in stellar spectra. This enables to segregate wavelength displacements caused by the radial motion of the stellar centre-of-mass from those induced by other effects, such as gravitational redshifts in white dwarfs.

We reobserved in the $R_C$ and $i'_{Sloan}$ bands, during the years 2020-2021, seven Mira variables in Cassiopeia, for which historical $i'_{Sloan}$ light curves were available from Asiago Observatory plates taken in the years 1967-84. The aim was to check if any of them had undergone a substantial change in the period or in the light curve shape. Very recent public data form ZTF-DR5 were also used to expand our time base window. A marked color change was detected for all the stars along their variability cycle. The star V890 Cas showed a significant period decrease of 12\% from 483 to 428 days, one of the largest known to date. All the stars, save AV Cas, showed a smaller variation amplitude in the recent CCD data, possibly due to a photometric accuracy higher than that of the photographic plates.

The Copernican principle (CP), i.e. the assumption that we are not privileged observers of the universe, is a fundamental tenet of the standard cosmological model. A violation of this postulate implies the possibility that the apparent cosmic acceleration could be explained without the need of a cosmological constant, dark energy or covariant modifications of gravity. In this letter we present a new test of the CP relating the distance and the expansion rate, derived via Noether's theorem, which is complementary to other tests found in the literature. We also simulate fiducial data based on upcoming stage IV galaxy surveys and use them to reconstruct the Hubble rate $H(z)$ and the angular diameter distance $d_A(z)$ in order to reconstruct our null test and forecast how well it can constrain deviations from the cosmological constant model. We find that our new test can easily rule out several realistic scenarios based on the Lemaitre-Tolman-Bondi void model at confidence of $\gtrsim 3\sigma$ at mid to high redshifts ($z>0.5$).

We present an in-depth analysis and results of eleven XMM-Newton datasets, spanning 2000 to 2016, of the anomalous X-ray Pulsar CXOU J010043.1$-$721134 which has been classified as a magnetar. We find a spin-period of 8.0275(1) s as of December 2016 and calculate the period derivative to be $(1.76\pm 0.02) \times 10^{-11}$ s s$^{-1}$, which translate to a dipolar magnetic field strength of $3.8\times 10^{14}$ G and characteristic age of $\sim 7200$ yr for the magnetar. It has a double-peaked pulse profile, with one broad and one narrow peak, in both soft ($0.3-1.3$ keV) and hard ($1.3-8$ keV) energy bands. The pulse fractions in the two energy bands are found to be consistent with constant values. These results are in agreement with previously published results for this source. Although two-component models produce acceptable fits to its energy spectra, single component models are much simpler and are able to explain the similarity of the pulse profiles in the low and high energy bands. We attempt fitting with four different single-component models and find that the best fit to the spectra is obtained by fitting a thermal Comptonization model with the photon index $(\Gamma)$ between $2.0-2.7$ and the electron temperature $(kT_e)$ between $0.5-0.9$ keV, for a seed blackbody photon distribution of 0.2 keV. Finally, we conclude by discussing our results briefly.

On page 10 of the 2018 National Academies Exoplanet Science Strategy document (NASEM 2018), 'Expect the unexpected' is described as a general principle of the exoplanet field. But for the next 150 pages, this principle is apparently forgotten, as strategy decisions are repeatedly put forward based on our expectations. This paper explores what exactly it might mean to 'expect the unexpected', and how this could possibly be achieved by the space science community. An analogy with financial investment strategies is considered, where a balanced portfolio of low/medium/high-risk investments is recommended. Whilst this kind of strategy would certainly be advisable in many scientific contexts (past and present), in certain contexts, especially exploratory science, a significant disanalogy needs to be factored in: financial investors cannot choose low-risk high-reward investments, but sometimes scientists can. The existence of low-risk high-impact projects in cutting-edge space science significantly reduces the warrant for investing in high-risk projects, at least in the short term. However, high-risk proposals need to be fairly judged alongside medium- and low-risk proposals, factoring in both the degree of possible reward and the expected cost of the project. Attitudes towards high-risk high-impact projects within NASA since 2009 are critically analysed.

Here we consider the perspective to detect sub-pc super-massive binary black-hole (SMBBH) systems using long-term photometric and spectroscopic monitoring campaigns of active galactic nuclei. This work explores the nature of long-term spectral variability caused by the dynamical effects of SMBBH systems. We describe in great detail a model of SMBBH system which considers that both black holes have their accretion disc and additional line emitting region(s). We simulate the H$\beta$ spectral band (continuum+broad H$\beta$ line) for different mass ratios of components and different total masses of the SMBBH systems ($10^6-10^8\mathrm{M\odot}$). We analyze the set of continuum and broad line light curves for several full orbits of SMBBHs with different parameters, to test the possibility to extract the periodicity of the system. We consider different levels of the signal-to-noise ratio, which is added to the simulated spectra. Our analysis showed that the continuum and broad line profiles emitted from an SMBBH system are strongly dependent, not only on the mass ratio of the components but also on the total mass of the system. We found that the mean broad line profile and its rms could indicate the presence of an SMBBH. However, some effects caused by the dynamics of a binary system could be hidden due to a low signal-to-noise ratio. Finally, we can conclude that the long-term AGN monitoring campaigns could be beneficial for the detection of SMBBH candidates.

The measurement of the energy spectrum of cosmic ray helium nuclei from 70 GeV to 80 TeV using 4.5 years of data recorded by the DArk Matter Particle Explorer (DAMPE) is reported in this work. A hardening of the spectrum is observed at an energy of about 1.3 TeV, similar to previous observations. In addition, a spectral softening at about 34 TeV is revealed for the first time with large statistics and well controlled systematic uncertainties, with an overall significance of $4.3\sigma$. The DAMPE spectral measurements of both cosmic protons and helium nuclei suggest a particle charge dependent softening energy, although with current uncertainties a dependence on the number of nucleons cannot be ruled out.

Population synthesis studies of binary black-hole mergers often lack robust black-hole spin estimates as they cannot accurately follow tidal spin-up during the late black-hole-Wolf-Rayet evolutionary phase. We provide an analytical approximation of the dimensionless second-born black-hole spin given the binary orbital period and Wolf-Rayet stellar mass at helium depletion or carbon depletion. These approximations are obtained from fitting a sample of around $10^5$ detailed MESA simulations that follow the evolution and spin up of close black-hole--Wolf-Rayet systems with metallicities in the range $[10^{-4},1.5Z_\odot]$. Following the potential spin up of the Wolf-Rayet progenitor, the second-born black-hole spin is calculated using up-to-date core collapse prescriptions that account for any potential disk formation in the collapsing Wolf-Rayet star. The fits for second-born black hole spin provided in this work can be readily applied to any astrophysical modeling that relies on rapid population synthesis, and will be useful for the interpretation of gravitational-wave sources using such models.

PSR J0540$-$6919 is the second-most energetic radio pulsar known and resides in the Large Magellanic Cloud. Like the Crab pulsar it is observed to emit giant radio pulses (GPs). We used the newly-commissioned PTUSE instrument on the MeerKAT radio telescope to search for GPs across three observations. In a total integration time of 5.7 hrs we detected 865 pulses above our 7$\sigma$ threshold. With full polarisation information for a subset of the data, we estimated the Faraday rotation measure, $\rm{RM}=-245.8 \pm 1.0$ rad m$^{-2}$ toward the pulsar. The brightest of these pulses is $\sim$ 60% linearly polarised but the pulse-to-pulse variability in the polarisation fraction is significant. We find that the cumulative GP flux distribution follows a power law distribution with index $-2.75 \pm 0.02$. Although the detected GPs make up only $\sim$ 10% of the mean flux, their average pulse shape is indistinguishable from the integrated pulse profile, and we postulate that there is no underlying emission. The pulses are scattered at L-band frequencies with the brightest pulse exhibiting a scattering time-scale of $\tau = 0.92 \pm 0.02$ ms at 1.2 GHz. We find several of the giants display very narrow-band "flux knots" similar to those seen in many Fast Radio Bursts, which we assert cannot be due to scintillation or plasma lensing. The GP time-of-arrival distribution is found to be Poissonian on all but the shortest time-scales where we find four GPs in six rotations, which if GPs are statistically independent is expected to occur in only 1 of 7000 observations equivalent to our data.

State of the art radial velocity (RV) exoplanet searches are limited by the effects of stellar magnetic activity. Magnetically active spots, plage, and network regions each have different impacts on the observed spectral lines, and therefore on the apparent stellar RV. Differentiating the relative coverage, or filling factors, of these active regions is thus necessary to differentiate between activity-driven RV signatures and Doppler shifts due to planetary orbits. In this work, we develop a technique to estimate feature-specific magnetic filling factors on stellar targets using only spectroscopic and photometric observations. We demonstrate linear and neural network implementations of our technique using observations from the solar telescope at HARPS-N, the HK Project at the Mt. Wilson Observatory, and the Total Irradiance Monitor onboard SORCE. We then compare the results of each technique to direct observations by the Solar Dynamics Observatory (SDO). Both implementations yield filling factor estimates that are highly correlated with the observed values. Modeling the solar RVs using these filling factors reproduces the expected contributions of the suppression of convective blueshift and rotational imbalance due to brightness inhomogeneities. Both implementations of this technique reduce the overall activity-driven RMS RVs from 1.64 m/s to 1.02 m/s, corresponding to a 1.28 m/s reduction in the RMS variation. The technique provides an additional 0.41 m/s reduction in the RMS variation compared to traditional activity indicators.

Studies of cluster mass and velocity anisotropy profiles are useful tests of dark matter models, and of the assembly history of clusters of galaxies. These studies might be affected by unknown systematics caused by projection effects. We aim at testing observational methods for the determination of mass and velocity anisotropy profiles of clusters of galaxies. Particularly, we focus on the MAMPOSSt technique (Mamon et al. 2013). We use results from two semi-analytic models of galaxy formation coupled with high-resolution N-body cosmological simulations, the catalog of De Lucia & Blaizot (2007) and the FIRE catalog based on the new GAlaxy Evolution and Assembly model. We test the reliability of the Jeans equation in recovering the true mass profile when full projected phase-space information is available. We examine the reliability of the MAMPOSSt method in estimating the true mass and velocity anisotropy profiles of the simulated halos when only projected phase-space information is available, as in observations. The spherical Jeans equation provides a reliable tool for the determination of cluster mass profiles, also for subsamples of tracers separated by galaxy color. Results are equally good for prolate and oblate clusters. Using only projected phase-space information, MAMPOSSt provides estimates of the mass profile with a standard deviation of 35-69 %, and a negative bias of 7-17 %, nearly independent of radius, and that we attribute to the presence of interlopers in the projected samples. The bias changes sign, that is, the mass is over-estimated, for prolate clusters with their major axis aligned along the line-of-sight. MAMPOSSt measures the velocity anisotropy profiles accurately in the inner cluster regions, with a slight overestimate in the outer regions, both for the whole sample of observationally-identified cluster members and separately for red and blue galaxies.

To interpret adaptive-optics observations of (216) Kleopatra, we need to describe an evolution of multiple moons, orbiting an extremely irregular body and including their mutual interactions. Such orbits are generally non-Keplerian and orbital elements are not constants. Consequently, we use a modified $N$-body integrator, which was significantly extended to include the multipole expansion of the gravitational field up to the order $\ell = 10$. Its convergence was verified against the `brute-force' algorithm. We computed the coefficients $C_{\ell m},S_{\!\ell m}$ for Kleopatra's shape, assuming a~constant bulk density. For solar-system applications, it was also necessary to implement a variable distance and geometry of observations. Our $\chi^2$ metric then accounts for the absolute astrometry, the relative astrometry (2nd moon with respect to 1st), angular velocities, and also silhouettes, constraining the pole orientation. This allowed us to derive the orbital elements of Kleopatra's two moons. Using both archival astrometric data and new VLT/SPHERE observations (ESO LP 199.C-0074), we were able to identify the true periods of the moons, $P_1 = (1.822359\pm0.004156)\,{\rm d}$, $P_2 = (2.745820\pm0.004820)\,{\rm d}$. They orbit very close to the 3:2 mean-motion resonance, but their osculating eccentricities are too small compared to other perturbations (multipole, mutual), so that regular librations of the critical argument are not present. The resulting mass of Kleopatra, $m_1 = (1.49\pm0.16)\cdot10^{-12}\,M_\odot$ or $2.97\cdot10^{18}\,{\rm kg}$, is significantly lower than previously thought. An implication explained in the accompanying paper (Marchis et al.) is that (216) Kleopatra is a critically rotating body.

We present the third Open Gravitational-wave Catalog (3-OGC) of compact-binary coalescences, based on the analysis of the public LIGO and Virgo data from 2015 through 2019 (O1, O2, O3a). Our updated catalog includes a population of 57 observations, including four binary black hole mergers that had not previously been reported. This consists of 55 binary black hole mergers and the two binary neutron star mergers GW170817 and GW190425. We find no additional significant binary neutron star or neutron star--black hole merger events. The most confident new detection is the binary black hole merger GW190925\_232845 which was observed by the LIGO Hanford and Virgo observatories with $\mathcal{P}_{\textrm{astro}} > 0.99$; its primary and secondary component masses are $20.2^{+3.9}_{-2.5} M_{\odot}$ and $15.6^{+2.1}_{-2.6} M_{\odot}$, respectively. We estimate the parameters of all binary black hole events using an up-to-date waveform model that includes both sub-dominant harmonics and precession effects. To enable deep follow-up as our understanding of the underlying populations evolves, we make available our comprehensive catalog of events, including the sub-threshold population of candidates, and the posterior samples of our source parameter estimates.

Star-forming regions have been proposed as potential Galactic cosmic-ray accelerators for decades. Cosmic-ray acceleration can be probed through observations of gamma-rays produced in inelastic proton-proton collisions, at GeV and TeV energies. In this paper, we analyze more than 11 years of Fermi-LAT data from the direction of Westerlund 2, one of the most massive and best-studied star-forming regions in our Galaxy. In particular, we investigate the characteristics of the bright pulsar PSR J1023-5746 that dominates the gamma-ray emission below a few GeV at the position of Westerlund 2, and the underlying extended source FGES J1023.3-5747. The analysis results in a clear identification of FGES J1023.3-5747 as the GeV counterpart of the TeV source HESS J1023-575, through its morphological and spectral properties. This identification provides new clues about the origin of the HESS J1023-575 gamma-ray emission, favouring a hadronic origin of the emission, powered by Westerlund 2, rather than a leptonic origin related to either the pulsar wind nebula associated with PSR J1023-5746 or the cluster itself. This result indirectly supports the hypothesis that star-forming regions can contribute to the cosmic-ray sea observed in our Galaxy

This Letter reports observations of an event that connects all major classes of solar eruptions: those that erupt fully into the heliosphere versus those that fail and are confined to the Sun, and those that eject new flux into the heliosphere, in the form of a flux rope, versus those that eject only new plasma in the form of a jet. The event originated in a filament channel overlying a circular polarity inversion line (PIL) and occurred on 2013-03-20 during the extended decay phase of the active region designated NOAA 12488/12501. The event was especially well-observed by multiple spacecraft and exhibited the well-studied null-point topology. We analyze all aspects of the eruption using SDO AIA and HMI, STEREO-A EUVI, and SOHO LASCO imagery. One section of the filament undergoes a classic failed eruption with cool plasma subsequently draining onto the section that did not erupt, but a complex structured CME/jet is clearly observed by SOHO LASCO C2 shortly after the failed filament eruption. We describe in detail the slow buildup to eruption, the lack of an obvious trigger, and the immediate reappearance of the filament after the event. The unique mixture of major eruption properties observed during this event places severe constraints on the structure of the filament channel field and, consequently, on the possible eruption mechanism.

We present a method of selecting quasars up to redshift $\approx$ 6 with random forests, a supervised machine learning method, applied to Pan-STARRS1 and WISE data. We find that, thanks to the increasing set of known quasars we can assemble a training set that enables supervised machine learning algorithms to become a competitive alternative to other methods up to this redshift. We present a candidate set for the redshift range 4.8 to 6.3 which includes the region around z = 5.5 where quasars are difficult to select due to photometric similarity to red and brown dwarfs. We demonstrate that under our survey restrictions we can reach a high completeness ($66 \pm 7 \%$ below redshift 5.6 / $83^{+6}_{-9}\%$ above redshift 5.6) while maintaining a high selection efficiency ($78^{+10}_{-8}\%$ / $94^{+5}_{-8}\%$). Our selection efficiency is estimated via a novel method based on the different distributions of quasars and contaminants on the sky. The final catalog of 515 candidates includes 225 known quasars. We predict the candidate catalog to contain an additional $148^{+41}_{-33}$ new quasars below redshift 5.6 and $45^{+5}_{-8}$ above and make the catalog publicly available. Spectroscopic follow-up observations of 37 candidates lead us to discover 20 new high redshift quasars (18 at $4.6\le z\le5.5$, 2 $z\sim5.7$). These observations are consistent with our predictions on efficiency. We argue that random forests can lead to higher completeness because our candidate set contains a number of objects that would be rejected by common color cuts, including one of the newly discovered redshift 5.7 quasars.

Short-lived radioactive nuclei (SLR) with mean lives below 100 Myr provide us with unique insights into current galactic nucleosynthetic events, as well as events that contributed to the material of our Solar System more that 4.6 Gyr ago. Here we present a statistical analysis of the ratios of these radioactive nuclei at the time of early Solar System (ESS) using both analytical derivations and Monte Carlo methods. We aim to understand the interplay between the production frequency and the mean lives of these isotopes, and its impact on their theoretically predicted ratios in the interstellar medium (ISM). We find that when the ratio of two SRLs, instead of the ratios of each single SLR relative to its stable or long-lived isotope, is considered, not only the uncertainties related to the galactic chemical evolution of the stable isotope are completely eliminated, but also the statistical uncertainties are much lower. We identify four ratios, 247Cm/129I, 107Pd/182Hf, 97Tc/98Tc, and 53Mn/97Tc, that have the potential to provide us with new insights into the r-, s-, and p-process nucleosynthesis at the time of the formation of the Sun, and need to be studied using variable stellar yields. Additionally, the latter two ratios need to be better determined in the ESS to allow us to fully exploit them to investigate the galactic sites of the p process.

Primordial gravitational waves (GWs) carry the imprints of the dynamics of the universe during its earliest stages. With a variety of GW detectors being proposed to operate over a wide range of frequencies, there is great expectation that observations of primordial GWs can provide us with an unprecedented window to the physics operating during inflation and reheating. In this work, we closely examine the effects of the regime of reheating on the spectrum of primordial GWs observed today. We consider a scenario wherein the phase of reheating is described by an averaged equation of state (EoS) parameter with an abrupt transition to radiation domination as well as a scenario wherein there is a gradual change in the effective EoS parameter to that of radiation due to the perturbative decay of the inflaton. We show that the perturbative decay of the inflaton leads to oscillations in the spectrum of GWs, which, if observed, can possibly help us decipher finer aspects of the reheating mechanism. We also examine the effects of a secondary phase of reheating arising due to a brief epoch driven possibly by an exotic, non-canonical, scalar field. Interestingly, we find that, for suitable values of the EoS parameter governing the secondary phase of reheating, the GWs can be of the strength as suggested by the recent NANOGrav observations. We conclude with a discussion of the wider implications of our analysis.

On September 2015, a century after Einstein's predictions of their existence, the first gravitational waves (GWs) direct detection was performed by LIGO. On August 17, 2017, the two Advanced LIGO and the Advanced Virgo interferometers detected a GW produced by two merging neutron stars. The subsequent localization of the source in the sky, thanks to the presence of a third detector, led to the detection of the electromagnetic counterpart and follow-up of the event by roughly 70 electromagnetic and neutrino telescopes. After the first two data taking runs (O1 and O2), the LIGO-Virgo network detected 11 GWs from 10 binary black holes and one binary neutron star. On April 1, 2019, Advanced Virgo and Advanced LIGO started their third observing period (O3). After an introduction on GW detection, I will give an overview on the Advanced Virgo detector design, with a description of the technical choices made before O3 and their consequences on the detector sensitivity. Finally, I will describe the planned upgrades for the Advanced Virgo+ project.

In this work we consider a scenario where the dark energy is a dynamical fluid whose energy density can be transferred to the dark matter via a coupling function proportional to the energy density of the dark energy. In particular, we investigate this model's ability to address the $S_8$ tension and find that against data from Planck, BAO and Pantheon the model 1) can significantly reduce the significance of the tension, 2) does so without exacerbating nor introducing any other tension (such as the $H_0$ tension) and 3) without worsening the fit to the considered data sets with respect to the $\Lambda$CDM model. We also test the model against data from weak lensing surveys such as KiDS and DES, and find that the model's ability to address the $S_8$ tension further improves, without a significant impact on any other parameter nor statistical measure.

Galactic cosmic rays are energetic particles important in the context of life. Many works have investigated the propagation of Galactic cosmic rays through the Sun's heliosphere. However, the cosmic ray fluxes in M dwarf systems are still poorly known. Studying the propagation of Galactic cosmic rays through the astrospheres of M dwarfs is important to understand the effect on their orbiting planets. Here, we focus on the planetary system GJ 436. We perform simulations using a combined 1D cosmic ray transport model and 1D Alfv\'en-wave-driven stellar wind model. We use two stellar wind set-ups: one more magnetically-dominated and the other more thermally-dominated. Although our stellar winds have similar magnetic field and velocity profiles, they have mass-loss rates two orders of magnitude different. Because of this, they give rise to two different astrosphere sizes, one ten times larger than the other. The magnetically-dominated wind modulates the Galactic cosmic rays more at distances < 0.2 au than the thermally-dominated wind due to a higher local wind velocity. Between 0.2 and 1 au the fluxes for both cases start to converge. However, for distances > 10 au, spatial diffusion dominates, and the flux of GeV cosmic rays is almost unmodulated. We find, irrespective of the wind regime, that the flux of Galactic cosmic rays in the habitable zone of GJ 436 (0.2 - 0.4 au) is comparable with intensities observed at Earth. On the other hand, around GJ 436 b (0.028 au), both wind regimes predict Galactic cosmic ray fluxes that are approximately $10^4$ times smaller than the values observed at Earth.

Determining accurate velocity measurements from observations of the Sun is of vital importance to solar physicists who are studying the wave dynamics in the solar atmosphere. Weak chromospheric absorption lines, due to dynamic events in the solar atmosphere, often consist of multiple spectral components. Isolating these components allows for the velocity field of the dynamic and quiescent regimes to be studied independently. However, isolating such components is particularly challenging due to the wide variety of spectral shapes present in the same dataset. MCALF provides a novel method and infrastructure to determine Doppler velocities in a large dataset. Each spectrum is fitted with a model adapted to its specific spectral shape.

On 31 August 2019, an interstellar comet was discovered as it passed through the Solar System (2I/Borisov). Based on initial imaging observations, 2I/Borisov appeared to be completely similar to ordinary Solar System comets - an unexpected characteristic after the multiple peculiarities of the only previous known interstellar visitor 1I/'Oumuamua. Spectroscopic investigations of 2I/Borisov identified the familiar cometary emissions from CN, C2, O I, NH2, OH, HCN, and CO, revealing a composition similar to that of carbon monoxide-rich Solar System comets. At temperatures >700 K, comets additionally show metallic vapors produced by the sublimation of metal-rich dust grains. However, due to the high temperature needed, observation of gaseous metals has been limited to bright sunskirting and sungrazing comets and giant star-plunging exocomets. Here we report spectroscopic detection of atomic nickel vapor in the cold coma of 2I/Borisov observed at a heliocentric distance of 2.322 au - equivalent to an equilibrium temperature of 180 K. Nickel in 2I/Borisov seems to originate from a short-lived nickelbearing molecule with a lifetime of $340^{+260}_{-200}$ s at 1 au and is produced at a rate of $0.9 \pm 0.3 \times 10^{22}$ atoms s$^{-1}$, or 0.002% relative to OH and 0.3% relative to CN. The detection of gas-phase nickel in the coma of 2I/Borisov is in line with the concurrent identification of this atom (as well as iron) in the cold comae of Solar System comets.

Neutrino non-standard interactions (NSI) with electrons are known to alter the picture of neutrino decoupling from the cosmic plasma. NSI modify both flavour oscillations through matter effects, and the annihilation and scattering between neutrinos and electrons and positrons in the thermal plasma. In view of the forthcoming cosmological observations, we perform a precision study of the impact of non-universal and flavour-changing NSI on the effective number of neutrinos, $N_{eff}$. We present the variation of $N_{eff}$ arising from the different NSI parameters and discuss the existing degeneracies among them, from cosmology alone and in relation to the current bounds from terrestrial experiments. Even though cosmology is generally less sensitive to NSI than these experiments, we find that future cosmological data would provide competitive and complementary constraints for some of the couplings and their combinations.

The close limit approximation of binary black hole is a powerful method to study gravitational-wave emission from highly non-linear geometries. In this work, we use it as a tool to model black hole spacetimes in theories of gravity with a new fundamental scalar degree of freedom. As an example, we consider Einstein-scalar-Gauss-Bonnet gravity, which admits as solution the Schwarzschild geometry as well as black holes with scalar hair. Accordingly, we find scalar perturbations growing unbounded around binary systems. This "dynamical scalarization" process is easier to trigger (i.e. occurs at lower values of the coupling constant of the theory) than the corresponding process for isolated black holes. Our results and framework highlight the fundamental role of the interaction during the collision of compact objects. They also emphasize the importance of having waveforms for black hole binaries in alternative theories, in order to consistently perform tests beyond General Relativity.

The oscillation of neutrino flavors, due to its interferometry nature, is extremely sensitive to the phase differences developing during the propagation of neutrinos. In this paper we investigate the effect of the Violation of Equivalence Principle (VEP) on the flavor oscillation probabilities of atmospheric and cosmic neutrinos observed at neutrino telescopes such as IceCube. Assuming a general parameterization of VEP, dubbed extended parameter space, we show that the synergy between the collected data of high energy atmospheric and cosmic neutrinos severely constrains the VEP parameters. Also, the projected sensitivity of IceCube-Gen2 to VEP parameters is discussed.

More than 130 pure rotational transitions of $^{46}$TiO, $^{47}$TiO, $^{48}$TiO, $^{49}$TiO, $^{50}$TiO, and $^{48}$Ti$^{18}$O are recorded using a high-resolution mm-wave supersonic jet spectrometer in combination with a laser ablation source. For the first time a mass-independent Dunham-like analysis is performed encompassing rare titanium monoxide isotopologues, and are compared to results from high-accuracy quantum-chemical calculations. The obtained parametrization reveals for titanium monoxide effects due to deviations from the Born-Oppenheimer approximation. Additionally, the dominant titanium properties enable an insight into the electronic structure of TiO by analyzing its hyperfine interactions. Further, based on the mass-independent analysis, the frequency positions of the pure rotational transitions of the short lived rare isotopologue $^{44}$TiO are predicted with high accuracy, i.e., on a sub-MHz uncertainty level. This allows for dedicated radio-astronomical searches of this species in core-collapse environments of supernovae.

Eddington-inspired Born-Infeld gravity is an important modification of Einstein's general relativity, which can give rise to non-singular cosmologies at the classical level, and avoid the end-stage singularity in a gravitational collapse process. In the Newtonian limit, this theory gives rise to a modified Poisson's equation, as a consequence of which stellar observables acquire model dependent corrections, compared to the ones computed in the low energy limit of general relativity. This can in turn be used to establish astrophysical constraints on the theory. Here, we obtain such a constraint using observational data from cataclysmic variable binaries. In particular, we consider the tidal disruption limit of the secondary star by a white dwarf primary. The Roche lobe filling condition of this secondary star is used to compute stellar observables in the modified gravity theory in a numerical scheme. These are then contrasted with the values obtained by using available data on these objects, via a Monte Carlo error progression method. This way, we are able to constrain the theory within $5\sigma$ confidence level.

Ongoing observations in the strong-field regime are in optimal agreement with general relativity, although current errors still leave room for small deviations from Einstein's theory. Here we summarise our recent results on superradiance of scalar and electromagnetic test fields in Kerr-like spacetimes, focusing mainly on the Konoplya--Zhidenko metric. We observe that, while for large deformations with respect to the Kerr case superradiance is suppressed, it can be nonetheless enhanced for small deformations. We also study the superradiant instability caused by massive scalar fields, and we provide a first estimate of the effect of the deformation on the instability timescale.

Although bulk peculiar motions are commonplace in the universe, most theoretical studies either bypass them, or take the viewpoint of the idealised Hubble-flow observers. As a result, the role of these peculiar flows remains largely unaccounted for, despite the fact that relative-motion effects have led to the misinterpretation of the observations in a number of occasions. Here, we examine the implications of large-scale peculiar flows for the interpretation of the deceleration parameter. We compare, in particular, the deceleration parameters measured by the Hubble-flow observers and by their bulk-flow counterparts. In so doing, we use Newtonian theory and general relativity and employ closely analogous theoretical tools, which allows for the direct and transparent comparison of the two studies. We find that the Newtonian relative-motion effects are generally too weak to make a difference between the two measurements. In relativity, however, the deceleration parameters measured in the two frames differ considerably, even at the linear level. This could deceive the unsuspecting observers to a potentially serious misinterpretation of the universe's global kinematic status.

The nuclear fission of very neuron-rich nuclei related to the r-process is essential for the termination of nucleosynthesis flows on the nuclear chart and the final abundances. Nevertheless, most of the available fission data for the r-process calculations are based on theory predictions, including phenomenological treatments. In this study, we calculated a series of nuclear fission distribution for neutron-rich nuclei away from the beta-stability line. As most of these nuclei are experimentally unknown, we are based on theoretical calculations based on the dynamical fission model with the Langevin method. We performed fission distribution calculations for neutron-rich actinoid nuclei, applicable to the r-process nucleosynthesis simulations. In the present paper, we compared the obtained mass and charge distributions with experimental data. We also show the results of the systematic behaviour of mass distribution for neutron-rich U and Fm isotopes.

The standard argument for the Lorentz invariance of the thermodynamic entropy in equilibrium is based on the assumption that it is possible to perform an adiabatic transformation whose only outcome is to accelerate a macroscopic body, keeping its rest mass unchanged. The validity of this assumption constitutes the very foundation of relativistic thermodynamics and needs to be tested in greater detail. We show that, indeed, such a transformation is always possible, at least in principle. The only two assumptions invoked in the proof are that there is at least one inertial reference frame in which the second law of thermodynamics is valid and that the microscopic theory describing the internal dynamics of the body is a field theory, with Lorentz invariant Lagrangian density. The proof makes no reference to the connection between entropy and probabilities and is valid both within classical and quantum physics. To avoid any risk of circular reasoning, we do not postulate that the laws of thermodynamics are the same in every reference frame, but we obtain this fact as a direct consequence of the Lorentz invariance of the entropy.