Papers for Wednesday, May 19 2021

The TNG50 cosmological simulation produces X-ray emitting bubbles, shells, and cavities in the circumgalactic gas above and below the stellar disks of Milky Way- and Andromeda-like galaxies with morphological features reminiscent of the eROSITA and Fermi bubbles in the Galaxy. Two-thirds of the 198 MW/M31 analogues simulated within TNG50 and inspected at z=0 show one or more large-scale, coherent features of over-pressurized gas that impinge into the gaseous halo. Some of the galaxies include a succession of bubbles or shells of increasing size, ranging from a few to many tens of kpc in height. These are prominent in gas pressure, X-ray emission and gas temperature, and often exhibit sharp boundaries indicative of shocks with typical Mach numbers of 2-4. The gas in the bubbles outflows with maximum (95th pctl) radial velocities of 100-1500 km/s. Across our sample, the bubbles expand with speeds as high as 1000-2000 km/s (about 1-2 kpc/Myr), but with a great diversity and with larger bubbles expanding at slower speeds. The bubble gas is typically at 10^6.4-7.2 K temperatures and is enriched to metallicities of 0.5-2 solar. In TNG50, the bubbles are a manifestation of episodic kinetic energy injections from the supermassive black holes at the galaxy centers that accrete at low Eddington ratios. According to TNG50, X-ray, and possibly gamma-ray, bubbles similar to those observed in the Milky Way should be a frequent feature of disk-like galaxies prior to, or on the verge of, being quenched. They should be within the grasp of eROSITA with a few ks observations of the local Universe.

E+A galaxies are believed to be a short phase connecting major merger ULIRGs with red and dead elliptical galaxies. Their optical spectrum suggests a massive starburst that was quenched abruptly, and their bulge-dominated morphologies with tidal tails suggest that they are merger remnants. AGN-driven winds are believed to be the process responsible for the sudden quenching of star formation and for the expulsion and/or destruction of the remaining molecular gas. So far, little was known observationally about AGN-driven winds in this short-lived phase. In this paper we present the first sample of E+A galaxies with AGN and indications of ionized outflows. Using IRAS-FIR observations, we study the star formation in these systems and find that, contrary to common belief, many E+A galaxies are not fully-quenched, with some showing powerful starbursts that are completely obscured at optical wavelengths. Using SDSS spectroscopy, we study the stationary and outflowing ionized gas. We also detect neutral gas outflows in 40\% of the sources. In these objects, the neutral outflow phase is 10--100 times more massive than the ionized phase. The mass outflow rate and kinetic power of the ionized outflows in E+A galaxies ($\dot{M}\sim 1\, \mathrm{M_{\odot}/yr}$, $\dot{E}\sim 10^{41}\, \mathrm{erg/sec}$) are larger than those derived for typical active galaxies. For the neutral outflow ($\dot{M}\sim 10\, \mathrm{M_{\odot}/yr}$, $\dot{E}\sim 10^{42}\, \mathrm{erg/sec}$) , they are smaller than those observed in (U)LIRGs with and without AGN. Our partial correlation analysis suggests that both AGN and SF contribute to the observed winds.

We show that relic vector fields can significantly impact a spectrum of primordial gravitational waves in the post-inflationary era. We consider a triplet of U(1) fields in a homogeneous, isotropic configuration. The interaction between the gravitational waves and the vector fields, from the end of reheating to the present day, yields novel spectral features. The amplitude, tilt, shape, and net chirality of the gravitational wave spectrum are shown to depend on the abundance of the electric- and magnetic-like vector fields. Our results show that even a modest abundance can have strong implications for efforts to detect the imprint of gravitational waves on the cosmic microwave background polarization. We find that a vector field comprising less than two percent of the energy density during the radiation dominated era can have a greater than order unity effect on the predicted inflationary gravitational wave spectrum.

We investigate the possibility of detecting artificial lights from Proxima b's dark side by computing light curves from the planet and its host star. The two different scenarios we consider are artificial illumination with the same spectrum as commonly used LEDs on Earth, and a narrower spectrum which leads to the same proportion of light as the total artificial illumination on Earth. We find that the James Webb Space Telescope (JWST) will be able to detect LED type artificial lights making up 5% of stellar power with 85% confidence, assuming photon-limited precision. In order for JWST to detect the current level of artificial illumination on Earth, the spectral band must be 10^3 times narrower. Our predictions require optimal performance from the NIRSpec instrument, and even if not possible with JWST, future observatories like LUVOIR might be able to detect this artificial illumination.

The field of galaxy evolution will make a great leap forward in the next decade as a consequence of the huge effort by the scientific community in multi-object spectroscopic facilities. To maximise the impact of such incoming data, the analysis methods must also step up, extracting reliable information from the available spectra. In this paper, we aim to investigate the limits and the reliability of different spectral synthesis methods in the estimation of the mean stellar age and metallicity. The main question this work aims to address is which signal-to-noise ratios (S/N) are needed to reliably determine the mean stellar age and metallicity from a galaxy spectrum and how this depends on the tool used to model the spectra. To address this question we built a set of realistic simulated spectra containing stellar and nebular emission, reproducing the evolution of a galaxy in two limiting cases: a constant star formation rate and an exponentially declining star formation. We degraded the synthetic spectra built from these two star formation histories (SFHs) to different S/N and analysed them with three widely used spectral synthesis codes, namely FADO, STECKMAP, and STARLIGHT. For S/N < 5 all three tools show a large diversity in the results. The FADO and STARLIGHT tools find median differences in the light-weighted mean stellar age of ~0.1 dex, while STECKMAP shows a higher value of ~0.2 dex. Detailed investigations of the best-fit spectrum for galaxies with overestimated mass-weighted quantities point towards the inability of purely stellar models to fit the observed spectra around the Balmer jump. Our results imply that when a galaxy enters a phase of high specific star formation rate the neglect of the nebular continuum emission in the fitting process has a strong impact on the estimation of its SFH when purely stellar fitting codes are used, even in presence of high S/N spectra.

Pulsar timing array (PTA) experiments are becoming increasingly sensitive to gravitational waves (GWs) in the nanohertz frequency range, where the main astrophysical sources are the supermassive black hole binaries (SMBHBs) which are expected to reside in the centers of galaxies. Some of these individual SMBHBs may power active galactic nuclei (AGNs), and thus their binary parameters could be obtained electromagnetically, which makes it possible to apply electromagnetic (EM) information to aid the search for a GW signal in PTA data. In this work, we investigate the effects of such an EM-informed search on binary detection and parameter estimation by performing mock data analyses on simulated PTA datasets. We find that by applying EM priors, the Bayes factor of the signal increases by a factor of a few to an order of magnitude, and thus an EM-informed targeted search is able to find hints of a signal when an uninformed search fails to find any. By combining EM and GW data, one can achieve an overall improvement in parameter estimation, regardless of the source's sky location or GW frequency. We discuss the implications for the multimessenger studies of SMBHBs with PTAs.

Cosmic ray (CR) hydrodynamics is a (re-)emerging field of high interest due to the importance of CRs for the dynamical evolution of the interstellar, the circumgalactic, and the intracluster medium. In these environments, CRs with GeV energies can influence large-scale dynamics by regulating star formation, driving galactic winds or by altering the pressure balance of galactic halos. Recent efforts have moved the focus of the community from a one-moment description of CR transport towards a two-moment model as this allows for a more accurate description of the microphysics impacting the CR population. Like all hydrodynamical theories, these two-moment methods require a closure relation for a consistent and closed set of evolution equations. The goal of this paper is to quantify the impact of different closure relations on the resulting solutions. To this end, we review the common P1 and M1 closure relations, derive a new four-moment H1 description for CR transport and describe how to incorporate CR scattering by Alfv\'en waves into these three hydrodynamical models. While there are significant differences in the transport properties of radiation in the P1 and M1 approximations in comparison to more accurate radiative transfer simulations using the discrete ordinates approximation, we only find small differences between the three hydrodynamical CR transport models in the free streaming limit when we neglect CR scattering. Most importantly, for realistic applications in the interstellar, circumgalactic or intracluster medium where CR scattering is frequent, these differences vanish and all presented hydrodynamical models produce the same results.

A subset of type Ic supernovae (SNe Ic), broad-lined SNe Ic (SNe Ic-bl), show unusually high kinetic energies ($\sim 10^{52}$ erg) which cannot be explained by the energy supplied by neutrinos alone. Many SNe Ic-bl have been observed in coincidence with long gamma-ray bursts (GRBs) which suggests a connection between SNe and GRBs. A small fraction of core-collapse supernovae (CCSNe) form a rapidly-rotating and strongly-magnetized protoneutron star (PNS), a proto-magnetar. Jets from such magnetars can provide the high kinetic energies observed in SNe Ic-bl and also provide the connection to GRBs. In this work we use the jetted outflow produced in a 3D CCSN simulation from a consistently formed proto-magnetar as the central engine for full-star explosion simulations. We extract a range of central engine parameters and find that the extracted engine energy is in the range of $6.231 \times 10^{51}-1.725 \times 10^{52}$ erg, the engine time-scale in the range of $0.479-1.159$ s and the engine half-opening angle in the range of $\sim 9-19^{\circ}$. Using these as central engines, we perform 2D special-relativistic (SR) hydrodynamic (HD) and radiation transfer simulations to calculate the corresponding light curves and spectra. We find that these central engine parameters successfully produce SNe Ic-bl which demonstrates that jets from proto-magnetars can be viable engines for SNe Ic-bl. We also find that only the central engines with smaller opening angles ($\sim 10^{\circ}$) form a GRB implying that GRB formation is likely associated with narrower jet outflows and Ic-bl's without GRBs may be associated with wider outflows.

In this article, the temperature-density relation of the intergalactic medium was studied in the region $1.6 \leq z < 2.0$ divided into two bins. For this purpose, the Ly-$\alpha$ forest decomposition into individual absorption profiles was used for the study of 35 publicly available quasar spectra obtained by the Ultraviolet and Visual Echelle Spectrograph (UVES) on the Very Large Telescope (ESO) and by the High Resolution Echelle Spectrometer (HIRES) on the Keck Telescope. For the determination of the thermal state sensitive cut-off position in the $b - N_{HI}$ distribution, the iterative fitting procedure was adopted. The measurements were calibrated using mock Ly-$\alpha$ forest data generated by 23 hydrodynamical simulations with different thermal histories. The value of the temperature at mean density corresponds to the decreasing trend predicted by various models at the lower redshifts. In the case of power law index, determined values are close to 1.6, which is expected after all reionization events in various models assuming the balance of photoheating with adiabatic cooling.

We present a comprehensive multi-frequency study of the HBL 1ES 1959+650 using data from various facilities during the period 2016-2017, including X-ray data from {\it AstroSat} and {\it Swift} during the historically high X-ray flux state of the source observed until February 2021. The unprecedented quality of X-ray data from high cadence monitoring with the {\it AstroSat} during 2016-2017 enables us to establish a detailed description of X-ray flares in 1ES 1959+650. The synchrotron peak shifts significantly between different flux states, in a manner consistent with a geometric (changing Doppler factor) interpretation. A time-dependent leptonic diffusive-shock-acceleration and radiation transfer model is used to reproduce the spectral energy distributions (SEDs) and X-ray light curves, to provide insight into the particle acceleration during the major activity periods observed in 2016 and 2017. The extensive data of {\it Swift}-XRT from December 2015 to February 2021 (Exp. = 411.3 ks) reveals a positive correlation between flux and peak position.

Despite the fact that different particle species can diffuse with respect to each other in neutron star (NS) cores, the effect of particle diffusion on various phenomena associated with NS oscillations is usually ignored. Here we demonstrate that the diffusion can be extremely powerful dissipative mechanism in superconducting NSs. In particular, it can be much more efficient than the shear and bulk viscosities. This result has important implications for the damping times of NS oscillations, development and saturation of dynamical instabilities in NSs, and for the excitation and coupling of oscillation modes during the late inspiral of binary NSs.

The flux of ultrahigh energy cosmic rays reaching the Earth is affected by the interactions with the cosmic radiation backgrounds as well as with the magnetic fields that are present along their trajectories. We combine the SimProp cosmic ray propagation code with a routine that allows to account for the average effects of a turbulent magnetic field on the direction of propagation of the particles. We compute in this way the modification of the spectrum which is due to the magnetic horizon effect, both for primary nuclei as well as for the secondary nuclei resulting from the photo-disintegration of the primary ones. We also provide analytic parameterizations of the attenuation effects, as a function of the magnetic field parameters and of the density of cosmic ray sources, which make it possible to obtain the expected spectra in the presence of the magnetic fields from the spectra that would be obtained in the absence of magnetic fields. The discrete nature of the distribution of sources with finite density also affects the spectrum of cosmic rays at the highest energies where the flux is suppressed due to the interactions with the radiation backgrounds, and parameterizations of these effects are obtained.

We use a sample of $78,340$ star-forming galaxies at $z\simeq 0.04-0.1$ from the SDSS DR8 survey to calculate the average nebular dust attenuation curve and its variation with the physical properties of galaxies. Using the first four low-order Balmer emission lines (H$\alpha$, H$\beta$, H$\gamma$, H$\delta$) detected in the composite spectrum of all galaxies in the sample, we derive a nebular attenuation curve in the range of $0.41\mu$m to $0.66\mu$m that has a similar shape and normalization to that of the Galactic extinction curve (Milky Way curve), the SMC curve and the nebular curve derived in a previous study for a high redshift sample. We divide the galaxies into bins of stellar mass, gas-phase metallicity, and specific star-formation rate, and derive the nebular attenuation curve in each of these bins. This analysis indicates that there is very little variation in the shape of the nebular dust attenuation curve with the properties used to bin the galaxies. Although this suggests that the combined effect of the dust composition and geometry along the average nebular sightlines of those bins results in the nebular curves with identical shapes in the optical region, but it is not possible to draw any conclusion about the similarities in the dust properties (geometry and composition) alongside those sightlines.

Magnetic arcades in the solar atmosphere, or coronal loops, are common structures known to host magnetohydrodynamic (MHD) waves and oscillations. Of particular interest are the observed properties of transverse loop oscillations, such as their frequency and mode of oscillation, which have received significant attention in recent years because of their seismological capability. Previous studies have relied on standard data analysis techniques, such as a fast Fourier transform (FFT) and wavelet transform (WT), to correctly extract periodicities and identify the MHD modes. However, how these methods can lead to artefacts requires investigation. We assess whether these two common spectral analysis techniques in coronal seismology can successfully identify high-frequency waves from an oscillating coronal loop. We examine extreme ultraviolet images of a coronal loop observed by the Atmospheric Imaging Assembly in the 171 \AA waveband on board the Solar Dynamics Observatory. We perform a spectral analysis of the loop waveform and compare our observation with a basic simulation. The spectral FFT and WT power of the observed loop waveform is found to reveal a significant signal with frequency 2.67 mHz superposed onto the dominant mode of oscillation of the loop (1.33 mHz), that is, the second harmonic of the loop. The simulated data show that the second harmonic is completely artificial even though both of these methods identify this mode as a real signal. This artificial harmonic, and several higher modes, are shown to arise owing to the periodic but non-uniform brightness of the loop. We further illustrate that the reconstruction of the 2.67 mHz component in the presence of noise yields a false perception of oscillatory behaviour that does not otherwise exist. We suggest that additional techniques such as a forward model of a 3D coronal arcade are necessary to verify such high-frequency waves.

The magnetic beta Cep pulsator xi^1 CMa has the longest rotational period of any known magnetic B-type star. It is also the only magnetic B-type star with magnetospheric emission that is known to be modulated by both rotation and pulsation. We report here the first unambiguous detection of a negative longitudinal magnetic field in xi^1 CMa (<Bz>=-87 +/- 2 G in 2019 and <Bz>=-207 +/- 3 G in 2020), as well as the results of ongoing monitoring of the star's Halpha variability. We examine evidence for deviation from a purely dipolar topology. We also report a new HST UV spectrum of xi^1 CMa obtained near magnetic null that is consistent with an equatorial view of the magnetosphere, as evidenced by its similarity to the UV spectrum of beta Cep obtained near maximum emission. The new UV spectrum of xi^1 CMa provides additional evidence for the extremely long rotation period of this star via comparison to archival data.

We extend the SDSS Stripe 82 Standard Stars Catalog with post-2007 SDSS imaging data. This improved version lists averaged SDSS ugriz photometry for nearly a million stars brighter than r~22 mag. With 2-3x more measurements per star, random errors are 1.4-1.7x smaller than in the original catalog, and about 3x smaller than for individual SDSS runs. Random errors in the new catalog are ~< 0.01 mag for stars brighter than 20.0, 21.0, 21.0, 20.5, and 19.0 mag in u, g, r, i, and z-bands, respectively. We achieve this error threshold by using the Gaia Early Data Release 3 (EDR3) Gmag photometry to derive gray photometric zeropoint corrections, as functions of R.A. and Declination, for the SDSS catalog, and use the Gaia BP-RP colour to derive corrections in the ugiz bands, relative to the r-band. The quality of the recalibrated photometry, tested against Pan-STARRS1, DES, CFIS and GALEX surveys, indicates spatial variations of photometric zeropoints <=0.01 mag (RMS), with typical values of 3-7 millimag in the R.A., and 1-2 millimag in the Declination directions, except for <~6 millimag scatter in the u-band. We also report a few minor photometric problems with other surveys considered here, including a magnitude-dependent ~0.01 mag bias between 16 < G_Gaia < 20 in the Gaia EDR3. Our new, publicly available catalog offers robust calibration of ugriz photometry below 1% level, and will be helpful during the commissioning of the Vera C. Rubin Observatory Legacy Survey of Space and Time.

We present a new method to assess the properties of transiting planet candidates by multicolor photometry. By analyzing multicolor transit/eclipse light curves and apparent magnitudes of the target in parallel, this method attempts to identify the nature of the system and provide a quantitative constraint on the properties of unresolved companion(s). We demonstrate our method by observing the six systems hosting candidate transiting planets, identified by the K2 mission (EPIC 206036749, EPIC 206500801, EPIC 210513446, EPIC 211800191, EPIC 220621087, and EPIC 220696233). Applying our analysis code to the six targets, we find that EPIC 206036749, EPIC 210513446, and EPIC 211800191 are likely to be triple-star systems including eclipsing binaries, and EPIC 220696233 is likely a planetary system, albeit further observations are required to confirm the nature. Additionally, we confirm that the systematic errors in the derived system parameters arising from adopting specific isochrone models and observing instruments (passbands) are relatively small. While this approach alone is not powerful enough to validate or refute planet candidates, the technique allows us to constrain the properties of resolved/unresolved companions, and prioritize the planet candidates for further follow-up observations (e.g., radial-velocity measurements).

The unknown $\gamma$-ray excess in the northwest region of Arp 220 was revisited by analyzing $\sim$11.8 years of the \textit{Fermi} Large Area Telescope (\textit{Fermi}-LAT) data in this study. We found that its photon flux was approximately three times higher than that of the previous study in the 0.2-100 GeV band, and the corresponding significance level$\sim8.15\sigma$ was approximately four times higher than before. The light curves of 15 and 45 time bins from the whole time all showed two active periods, and the variability of the second period was more significant than that of the first period. The spectral indices from the two active periods were not statistically different and were close to the range of $\gamma$-ray flat-spectrum radio quasars observed by \textit{Fermi}-LAT. Because the position of CRATES J153246+234400 was consistent with the best-fit position of our analysis, we suggest that CRATES J153246+234400 is more likely a $\gamma$-ray counterpart for the variational region. For Arp 220, there was no significant variability in the $\gamma$-ray emission.

We present analytic and numerical models of a cluster wind flow resulting from the interaction of stellar winds of massive stars, with a super massive black hole (SMBH). We consider the motion of the stars as well as the gravitational force of the SMBH. In the numerical simulations we consider two cases: the first one with the stars is in circular orbits, and the second one with the stars in eccentric orbits around the SMBH. We found that after the system reaches an equilibrium, the circular and elliptical cases are very similar. We found a very good agreement between the analytical and numerical results, not only from our numerical simulations but also from other high resolution numerical calculations. The analytical models are very interesting, since the properties of such complex systems involving strong winds and a massive compact object, can be rapidly inferred without the need of a numerical calculation.

Aims. We present the first detailed analysis of the astrophysical parameters of the poorly studied Sco-Cen member HD 152384 and its circumstellar environment. Methods. We analyze newly obtained optical-near-IR XSHOOTER spectra, as well as archival TESS data, of HD 152384. In addition, we use literature photometric data to construct a detailed spectral energy distribution (SED) of the star. Results. The photospheric absorption lines in the spectrum of HD 152384 are characteristic of a A0 V star, for which we derive a stellar mass of 2.1 +/- 0.1 M_sun and a stellar age > 4.5 Myr. Superimposed on the photospheric absorption, the optical spectrum also displays double-peaked emission lines of Ca II, Fe I, Mg I and Si I, typical of circumstellar disks. Notably, all Hydrogen and Helium lines appear strictly in absorption. A toy model shows that the observed emission line profiles can be reproduced by emission from a compact (radius < 0.3 au) disk seen at an inclination of ~24 degrees. Further evidence for the presence of circumstellar material comes from the detection of a moderate infrared excess in the SED, similar to those found in extreme debris disk systems. Conclusions. We conclude that HD 152384 is surrounded by a tenuous circumstellar disk which, although rich in refractory elements, is highly depleted of volatile elements. To the best of our knowledge such a disk is unique within the group of young stars. However, it is reminiscent of the disks seen in some white dwarfs, which have been attributed to the disruption of rocky planets. We suggest that the disk around HD 152384 may have a similar origin and may be due to collisions in a newly formed planetary system.

The precision measurements of the cosmic microwave background power spectrum put a strong constraint on the dark matter annihilation cross section since the electromagnetic energy injection by the dark matter annihilation affects the ionization history of the universe. In this paper, we update our previous simulation code for calculating the ionization history with the effect of dark matter annihilation by including Helium interactions and improving the precision of calculations. We give an updated constraint on the annihilation cross section and mass of dark matter using the modified CosmoMC code with the Planck 2018 datasets.

Context. The TOPoS project has the goal to find and analyse Turn-Off (TO) stars of extremely low metallicity. To select the targets for spectroscopic follow-up at high spectral resolution, we have relied on low-resolution spectra from the Sloan Digital Sky Survey. Aims. In this paper we use the metallicity estimates we have obtained from our analysis of the SDSS spectra to construct the metallicity distribution function (MDF) of the Milky Way, with special emphasis on its metal-weak tail. The goal is to provide the underlying distribution out of which the TOPoS sample was extracted. Methods. We make use of SDSS photometry, Gaia photometry and distance estimates derived from the Gaia parallaxes to derive a metallicity estimate for a large sample of over 24 million TO stars. This sample is used to derive the metallicity bias of the sample for which SDSS spectra are available. Results. We determined that the spectroscopic sample is strongly biased in favour of metal-poor stars, as intended. A comparison with the unbiased photometric sample allows to correct for the selection bias. We select a sub-sample of stars with reliable parallaxes for which we combine the SDSS radial velocities with Gaia proper motions and parallaxes to compute actions and orbital parameters in the Galactic potential. This allows us to characterize the stars dynamically, and in particular to select a sub-sample that belongs to the Gaia-Sausage-Enceladus (GSE) accretion event. We are thus able to provide also the MDF of GSE. Conclusions. The metal-weak tail derived in our study is very similar to that derived in the H3 survey and in the Hamburg/ESO Survey. This allows us to average the three MDFs and provide an error bar for each metallicity bin. Inasmuch the GSE structure is representative of the progenitor galaxy that collided with the Milky Way, that galaxy appears to be strongly deficient in metal-poor stars compared to the Milky Way, suggesting that the metal-weak tail of the latter has been largely formed by accretion of low mass galaxies rather than massive galaxies, such as the GSE progenitor.

Observations carried out toward starless and pre-stellar cores have revealed that complex organic molecules are prevalent in these objects, but it is unclear what chemical processes are involved in their formation. Recently, it has been shown that complex organics are preferentially produced at an intermediate-density shell within the L1544 pre-stellar core at radial distances of ~4000 au with respect to the core center. However, the spatial distribution of complex organics has only been inferred toward this core and it remains unknown whether these species present a similar behaviour in other cores. We report high-sensitivity observations carried out toward two positions in the L1498 pre-stellar core, the dust peak and a position located at a distance of ~11000 au from the center of the core where the emission of CH$_3$OH peaks. Similarly to L1544, our observations reveal that small O-bearing molecules and N-bearing species are enhanced by factors ~4-14 toward the outer shell of L1498. However, unlike L1544, large O-bearing organics such as CH3CHO, CH3OCH3 or CH3OCHO are not detected within our sensitivity limits. For N-bearing organics, these species are more abundant toward the outer shell of the L1498 pre-stellar core than toward the one in L1544. We propose that the differences observed between O-bearing and N-bearing species in L1498 and L1544 are due to the different physical structure of these cores, which in turn is a consequence of their evolutionary stage, with L1498 being younger than L1544.

Nowadays, we know thousands of exoplanets, some of them potentially habitable. Next technological facilities (JWST, for example) have exoplanet atmosphere analysis capabilities, but they also have limits in terms of how many targets can be studied. Therefore, there is a need to rank and prioritize these exoplanets with the aim of searching for biomarkers. Some criteria involved, such as the habitability potential of a dry-rock planet versus a water-rich planet, or a potentially-locked planet versus a tidally-locked planet, are often vague and the use of the fuzzy set theory is advisable. We have applied a combination of Multi-Criteria Decision-Making methodologies with fuzzy logic, the Fuzzy Reference Ideal Method (FRIM), to this problem. We have analyzed the habitability potential of 1798 exoplanets from TEPCat database based on a set of criteria (composition, atmosphere, energy, tidal locking, type of planet and liquid water), in terms of their similarity to the only ideal alternative, The Earth. Our results, when compared with the probability index SEPHI, indicate that Kepler-442b, Kepler-062e/f, and LHS_1140b are the best exoplanets for searching for biomarkers, regardless its technical difficulty. If we take into account current technical feasibility, the best candidate is TRAPPIST-1e.

MOA-2006-BLG-074 was selected as one of the most promising planetary candidates in a retrospective analysis of the MOA collaboration: its asymmetric high-magnification peak can be perfectly explained by a source passing across a central caustic deformed by a small planet. However, after a detailed analysis of the residuals, we have realized that a single lens and a source orbiting with a faint companion provides a more satisfactory explanation for all the observed deviations from a Paczynski curve and the only physically acceptable interpretation. Indeed the orbital motion of the source is constrained enough to allow a very good characterization of the binary source from the microlensing light curve. The case of MOA-2006-BLG-074 suggests that the so-called xallarap effect must be taken seriously in any attempts to obtain accurate planetary demographics from microlensing surveys.

X-ray polarimetry promises us an unprecedented look at the structure of magnetic fields and on the processes at the base of acceleration of particles up to ultrarelativistic energies in relativistic jets. Crucial pieces of information are expected from observations of blazars (that are characterized by the presence of a jet pointing close to the Earth), in particular of the subclass defined by a synchrotron emission extending to the X-ray band (so-called high synchrotron peak blazars, HSP). In this review, I give an account of some of the models and numerical simulations developed to predict the polarimetric properties of HSP at high energy, contrasting the predictions of scenarios assuming particle acceleration at shock fronts with those that are based on magnetic reconnection, and I discuss the prospects for the observations of the upcoming Imaging X-ray Polarimetry Explorer (IXPE) satellite.

We present SOFIA/FIFI-LS observations of three Class 0 and one Class I outflows (Cep E, HH 1, HH 212, and L1551 IRS5) in the far-infrared [O I]63mum and [O I]145mum transitions. Spectroscopic [O I]63mum maps enabled us to infer the spatial extent of warm, low-excitation atomic gas within these protostellar outflows. If proper shock conditions prevail, the instantaneous mass-ejection rate is directly connected to the [O I]63mum luminosity. In order to unravel evolutionary trends, we analysed a set of 14 Class 0/I outflow sources that were spatially resolved in the [O I]63mum emission. We compared these data with a sample of 72 Class 0/I/II outflow sources that have been observed with Herschel (WISH, DIGIT, WILL, GASPS surveys) without spatially resolving the [O I]63mum line.

The X-ray transient MAXI J1631-479 went into outburst on 2018 December 21 and remained active for about seven months. Owing to various constraints it was monitored by NICER only during the decay phase of the outburst for about four months. The NICER observations were primarily in the soft state with a brief excursion to the hard intermediate state. While the soft state spectrum was dominated by thermal disk emission, the hard intermediate state spectrum had maximum contribution from the thermal Comptonization. Almost all intermediate-state power spectra had a Type-C low frequency quasi-periodic oscillation (within 4 - 10 Hz), often accompanied by a harmonic component. The frequency of these oscillations increased and the fractional rms decreased with inner-disk temperature suggesting a geometric origin. One observation in the middle of the outburst during the hard intermediate state had two non-harmonically related peaks. While one of them was definitely a Type-C QPO, the identification of the other one is uncertain. The rms spectra during the intermediate state had a hard shape from above 1 keV. Below 1 keV the shape could not be constrained in most cases, while only a few observations showed a rise in amplitude.

(abridged) Context: The ratio of kinematic viscosity to thermal diffusivity, the Prandtl number, is much smaller than unity in stellar convection zones. Aims: To study the statistics of convective flows and energy transport as functions of the Prandtl number. Methods: Three-dimensional numerical simulations convection in Cartesian geometry are used. The convection zone (CZ) is embedded between two stably stratified layers. Statistics and transport properties of up- and downflows are studied separately. Results: The rms velocity increases with decreasing Prandtl number. At the same time the filling factor of downflows decreases and leads to stronger downflows at lower Prandtl numbers, and to a strong dependence of overshooting on the Prandtl number. Velocity power spectra do not show marked changes as a function of Prandtl number. At the highest Reynolds numbers the velocity power spectra are compatible with the Bolgiano-Obukhov $k^{-11/5}$ scaling. The horizontally averaged convected energy flux ($\overline{F}_{\rm conv}$) is independent of the Prandtl number within the CZ. However, the upflows (downflows) are the dominant contribution to the convected flux at low (high) Prandtl number. These results are similar to those from Rayleigh-Ben\'ard convection in the low Prandtl number regime where convection is vigorously turbulent but inefficient at transporting energy. Conclusions: The current results indicate a strong dependence of convective overshooting and energy flux on the Prandtl number. Numerical simulations of astrophysical convection often use Prandtl number of unity. The current results suggest that this can lead to misleading results and that the astrophysically relevant low Prandtl number regime is qualitatively different from the parameters regimes explored in typical simulations.

Discovery of pulsars is one of the main goals for large radio telescopes. The Five-hundred-meter Aperture Spherical radio Telescope (FAST), that incorporates an L-band 19-beam receiver with a system temperature of about 20~K, is the most sensitive radio telescope utilized for discovering pulsars. We designed the {\it snapshot} observation mode for a FAST key science project, the Galactic Plane Pulsar Snapshot (GPPS) survey, in which every four nearby pointings can observe {\it a cover} of a sky patch of 0.1575 square degrees through beam-switching of the L-band 19-beam receiver. The integration time for each pointing is 300 seconds so that the GPPS observations for a cover can be made in 21 minutes. The goal of the GPPS survey is to discover pulsars within the Galactic latitude of $\pm10^{\circ}$ from the Galactic plane, and the highest priority is given to the inner Galaxy within $\pm5^{\circ}$. Up to now, the GPPS survey has discovered 201 pulsars, including currently the faintest pulsars which cannot be detected by other telescopes, pulsars with extremely high dispersion measures (DMs) which challenge the currently widely used models for the Galactic electron density distribution, pulsars coincident with supernova remnants, 40 millisecond pulsars, 16 binary pulsars, some nulling and mode-changing pulsars and rotating radio transients (RRATs). The follow-up observations for confirmation of new pulsars have polarization-signals recorded for polarization profiles of the pulsars. Re-detection of previously known pulsars in the survey data also leads to significant improvements in parameters for 64 pulsars. The GPPS survey discoveries are published and will be updated at this http URL .

The dark matter (DM) can consist of the primordial black holes (PBHs) in addition to the conventional weakly interacting massive particles (WIMPs). The Poisson fluctuations of the PBH number density produce the isocurvature perturbations which can dominate the matter power spectrum at small scales and enhance the early structure formation. We study how the WIMP annihilation from those early formed structures can affect the CMB (in particular the E-mode polarization anisotropies and $y$-type spectral distortions) and global 21cm signals. Our studies would be of particular interest for the light (sub-GeV) WIMP scenarios which have been less explored compared with the mixed DM scenarios consisting of PBHs and heavy ($\gtrsim 1$ GeV) WIMPs. For instance, for the self-annihilating DM mass $m_{\chi}=1$ MeV and the thermally averaged annihilation cross section $\langle \sigma v \rangle \sim 10^{-30} \rm cm^3/s$, the latest Planck CMB data requires the PBH fraction with respect to the whole DM to be at most ${\cal O}(10^{-3})$ for the sub-solar mass PBHs and an even tighter bound (by a factor $\sim 5$) can be obtained from the global 21-cm measurements.

The Italian AGILE space mission, with its Gamma-Ray Imaging Detector (GRID) instrument sensitive in the 30 MeV-50 GeV gamma-ray energy band, has been operating since 2007. Agilepy is an open-source Python package to analyse AGILE/GRID data. The package is built on top of the command-line version of the AGILE Science Tools, developed by the AGILE Team, publicly available and released by ASI/SSDC. The primary purpose of the package is to provide an easy to use high-level interface to analyse AGILE/GRID data by simplifying the configuration of the tasks and ensuring straightforward access to the data. The current features are the generation and display of sky maps and light curves, the access to gamma-ray sources catalogues, the analysis to perform spectral model and position fitting, the wavelet analysis. Agilepy also includes an interface tool providing the time evolution of the AGILE off-axis viewing angle for a chosen sky region. The Flare Advocate team also uses the tool to analyse the data during the daily monitoring of the gamma-ray sky. Agilepy (and its dependencies) can be easily installed using Anaconda.

How simple organic matter appeared on Earth and the processes by which it transformed into more evolved organic compounds, which ultimately led to the emergence of life, is still an open topic. Different scenarios have been proposed, the main one assumes that simple organic compounds were synthesized, either in the gas phase or on the surfaces of dust grains, during the process of star formation, and were incorporated into larger bodies in the protoplanetary disk. Transformation of these simple organic compounds in more complex forms is still a matter of debate. Recent discoveries point out to catalytic properties of dust grains present in the early stellar envelope, which can nowadays be found in the form of chondrites. The huge infall of chondritic meteorites during the early periods of Earth suggests that the same reactions could have taken place in certain environments of the Earth surface, with conditions more favorable for organic synthesis. This work attempts the synthesis of simple organic molecules, such as hydrocarbons and alcohols, via Fischer-Tropsch Type reactions supported by different chondritic materials under early-Earth conditions, to investigate if organic synthesis can likely occur in this environment and which are the differences in selectivity when using different types of chondrites. Fischer-Tropsch-type reactions are investigated from mixtures of CO and H2 at 1 atm of pressure on the surfaces of different chondritic samples. The different products obtained are analyzed in situ by gas chromatography. Different Fischer-Tropsch reaction products are obtained in quantitative amounts. The formation of alkanes and alkenes being the main processes. Formation of alcohols also takes place in a smaller amount. Other secondary products were obtained in a qualitative way.

Space debris is a major problem in space exploration. International bodies continuously monitor a large database of orbiting objects and emit warnings in the form of conjunction data messages. An important question for satellite operators is to estimate when fresh information will arrive so that they can react timely but sparingly with satellite maneuvers. We propose a statistical learning model of the message arrival process, allowing us to answer two important questions: (1) Will there be any new message in the next specified time interval? (2) When exactly and with what uncertainty will the next message arrive? The average prediction error for question (2) of our Bayesian Poisson process model is smaller than the baseline in more than 3 hours in a test set of 50k close encounter events.

Observations have confirmed the existence of multiple-planet systems containing a hot Jupiter and smaller planetary companions. Examples include WASP-47, Kepler-730, and TOI-1130. We examine the plausibility of forming such systems in situ using $N$-body simulations that include a realistic treatment of collisions, an evolving protoplanetary disc and eccentricity/inclination damping of planetary embryos. Initial conditions are constructed using two different models for the core of the giant planet: a 'seed-model' and an 'equal-mass-model'. The former has a more massive protoplanet placed among multiple small embryos in a compact configuration. The latter consists only of equal-mass embryos. Simulations of the seed-model lead to the formation of systems containing a hot Jupiter and super-Earths. The evolution consistently follows four distinct phases: early giant impacts; runaway gas accretion onto the seed protoplanet; disc damping-dominated evolution of the embryos orbiting exterior to the giant; a late chaotic phase after dispersal of the gas disc. Approximately 1% of the equal-mass simulations form a giant and follow the same four-phase evolution. Synthetic transit observations of the equal-mass simulations provide an occurrence rate of 0.26% for systems containing a hot Jupiter and an inner super-Earth, similar to the 0.2% occurrence rate from actual transit surveys, but simulated hot Jupiters are rarely detected as single transiting planets, in disagreement with observations. A subset of our simulations form two close-in giants, similar to the WASP-148 system. The scenario explored here provides a viable pathway for forming systems with unusual architectures, but does not apply to the majority of hot Jupiters.

Using data from the literature, we made a list of individual estimates of the solar Galactocentric distance, which were performed after 2017 by different methods. These values have not yet been used to calculate the best value of mean $R_0$. For the sample containing 21 estimates, based on the standard approach, we found the weighted mean ${\overline R_0}=8.14$ kpc with the dispersion $\sigma=0.16$ kpc, and using the median statistics, we obtained the estimate $R_0=8.15\pm0.11$ kpc. For practical use, the value $R_0=8.1\pm0.1$ kpc can be recommended.

We present the discovery of NGTS-19b, a high mass transiting brown dwarf discovered by the Next Generation Transit Survey (NGTS). We investigate the system using follow up photometry from the South African Astronomical Observatory, as well as sector 11 TESS data, in combination with radial velocity measurements from the CORALIE spectrograph to precisely characterise the system. We find that NGTS-19b is a brown dwarf companion to a K-star, with a mass of $69.5 ^{+5.7}_{-5.4}$ M$_{Jup}$ and radius of $1.034 ^{+0.055}_{-0.053}$ R$_{Jup}$. The system has a reasonably long period of 17.84 days, and a high degree of eccentricity of $0.3767 ^{+0.0061}_{-0.0061}$. The mass and radius of the brown dwarf imply an age of $0.46 ^{+0.26}_{-0.15}$ Gyr, however this is inconsistent with the age determined from the host star SED, suggesting that the brown dwarf may be inflated. This is unusual given that its large mass and relatively low levels of irradiation would make it much harder to inflate. NGTS-19b adds to the small, but growing number of brown dwarfs transiting main sequence stars, and is a valuable addition as we begin to populate the so called brown dwarf desert.

Context. There are typically two different approaches to infer the mass formation history (MFH) of a given galaxy from its luminosity in different bands. Non-parametric methods are known for their flexibility and accuracy, while parametric models are more computationally efficient. Aims. In this work we propose an alternative that combines the advantages of both techniques, based on a polynomial expansion around the present time. Methods. In our approach, the MFH is decomposed through an orthonormal basis of N polynomia in lookback time. To test the proposed framework, synthetic observations are generated from models based on common analytical approximations (exponential, delayed-tau and Gaussian star formation histories). A normalized distance is used to measure the quality of the fit, and the input MFH are compared with the polynomial reconstructions both at the present time as well as through cosmic evolution. Results. The observed luminosities are reproduced with an accuracy of around 10 per cent for a constant star formation rate (N=1) and better for higher-order polynomia. Our method provides good results on the reconstruction of the total stellar mass, star formation rate and even its first derivative for smooth star formation histories, but it has difficulties in reproducing variations on short timescales and/or star formation histories peaking at the earliest times of the Universe. Conclusions. The polynomial expansion appears to be a promising alternative to other analytical functions used in parametric methods, combining both speed and flexibility.

X-ray observations of kilo-parsec scale jets indicate that a synchrotron origin of the sustained non-thermal emission is likely. This requires distributed acceleration of electrons up to near PeV energies along the jet. The underlying acceleration mechanism is still unclear. Shear acceleration is a promising candidate, as velocity-shear stratification is a natural consequence of the collimated flow of a jet. We study the details of shear acceleration by solving the steady-state Fokker-Planck-type equation and provide a simple general solution for trans-relativistic jets for a range of magnetohydrodynamic turbulent power-law spectra. In general, the accelerated particle population is a power-law spectrum with an exponential-like cut-off, where the power-law index is determined by the turbulence spectrum and the balance of escape and acceleration of particles. Adopting a simple linearly decreasing velocity profile in the boundary of large-scale jets, we find that the multi-wavelength spectral energy distribution of X-ray jets, such as Centaurus A and 3C 273, can be reproduced with electrons that are accelerated up to $\sim$ PeV. In kpc-scale jets, protons may be accelerated up to $\sim$ EeV, supporting the hypothesis that large-scale jets are strong candidates for ultra-high-energy-cosmic-ray sources within the framework of shear acceleration.

The analysis of individual X-ray sources that appear in a crowded field can easily be compromised by the misallocation of recorded events to their originating sources. Even with a small number of sources, that nonetheless have overlapping point spread functions, the allocation of events to sources is a complex task that is subject to uncertainty. We develop a Bayesian method designed to sift high-energy photon events from multiple sources with overlapping point spread functions, leveraging the differences in their spatial, spectral, and temporal signatures. The method probabilistically assigns each event to a given source. Such a disentanglement allows more detailed spectral or temporal analysis to focus on the individual component in isolation, free of contamination from other sources or the background. We are also able to compute source parameters of interest like their locations, relative brightness, and background contamination, while accounting for the uncertainty in event assignments. Simulation studies that include event arrival time information demonstrate that the temporal component improves event disambiguation beyond using only spatial and spectral information. The proposed methods correctly allocate up to 65% more events than the corresponding algorithms that ignore event arrival time information. We apply our methods to two stellar X-ray binaries, UV Cet and HBC515 A, observed with Chandra. We demonstrate that our methods are capable of removing the contamination due to a strong flare on UV Cet B in its companion approximately 40 times weaker during that event, and that evidence for spectral variability at timescales of a few ks can be determined in HBC515 Aa and HBC515 Ab.

In the multi-messenger era, astronomical projects share information about transients phenomena issuing science alerts to the Scientific Community through different communications networks. This coordination is mandatory to understand the nature of these physical phenomena. For this reason, astrophysical projects rely on real-time analysis software pipelines to identify as soon as possible transients (e.g. GRBs), and to speed up external alerts' reaction time. These pipelines can share and receive the science alerts through the Gamma-ray Coordinates Network. This work presents a framework designed to simplify the development of real-time scientific analysis pipelines. The framework provides the architecture and the required automatisms to develop a real-time analysis pipeline, allowing the researchers to focus more on the scientific aspects. The framework has been successfully used to develop real-time pipelines for the scientific analysis of the AGILE space mission data. It is planned to reuse this framework for the Super-GRAWITA and AFISS projects. A possible future use for the Cherenkov Telescope Array (CTA) project is under evaluation.

We report the first discovery of a transiting circumbinary planet detected from a single sector of TESS data. During Sector 21, the planet TIC 172900988b transited the primary star and then 5 days later it transited the secondary star. The binary is itself eclipsing, with a period of P = 19.7 days and an eccentricity of e = 0.45. Archival data from ASAS-SN, Evryscope, KELT, and SuperWASP reveal a prominent apsidal motion of the binary orbit, caused by the dynamical interactions between the binary and the planet. A comprehensive photodynamical analysis of the TESS, archival and follow-up data yields stellar masses and radii of M1 = 1.2377 +/- 0.0009 MSun and R1 = 1.3827 +/- 0.0025 RSun for the primary and M2 = 1.2021 +/- 0.0008 MSun and R2 = 1.3141 +/- 0.0013 RSun for the secondary. The radius of the planet is R3 = 11.07 +/- 0.47 REarth (1.009 +/- 0.043 RJup). The planet's mass and orbital properties are not uniquely determined - there are six solutions with nearly equal likelihood. Specifically, we find that the planet's mass is in the range of 822 < M3 < 981 MEarth (2.59 < M3 < 3.09 MJup), its orbital period could be 188.8, 190.4, 194.0, 199.0, 200.4, or 204.1 days, and the eccentricity is between 0.02 and 0.09. At a V = 10.141 mag, the system is accessible for high-resolution spectroscopic observations, e.g. Rossiter-McLaughlin effect and transit spectroscopy.

We examine the possibility of "soft cosmology", namely small deviations from the usual framework due to the effective appearance of soft-matter properties in the Universe sectors. One effect of such a case would be the dark energy to exhibit a different equation-of-state parameter at large scales (which determine the universe expansion) and at intermediate scales (which determine the sub-horizon clustering and the large scale structure formation). Concerning soft dark matter, we show that it can effectively arise due to the dark-energy clustering, even if dark energy is not soft. We propose a novel parametrization introducing the "softness parameters" of the dark sectors. As we see, although the background evolution remains unaffected, due to the extreme sensitivity and significant effects on the global properties even a slightly non-trivial softness parameter can improve the clustering behavior and alleviate e.g. the $f\sigma_8$ tension. Lastly, an extension of the cosmological perturbation theory and a detailed statistical mechanical analysis, in order to incorporate complexity and estimate the scale-dependent behavior from first principles, is necessary and would provide a robust argumentation in favour of soft cosmology.

The Cherenkov Telescope Array (CTA) is an initiative that is currently building the largest gamma-ray ground Observatory that ever existed. A Science Alert Generation (SAG) system, part of the Array Control and Data Acquisition (ACADA) system of the CTA Observatory, analyses online the telescope data - arriving at an event rate of tens of kHz - to detect transient gamma-ray events. The SAG system also performs an online data quality analysis to assess the instruments' health during the data acquisition: this analysis is crucial to confirm good detections. A Python and a C++ software library to perform the online data quality analysis of CTA data, called rta-dq-lib, has been proposed for CTA. The Python version is dedicated to the rapid prototyping of data quality use cases. The C++ version is optimized for maximum performance. The library allows the user to define, through XML configuration files, the format of the input data and, for each data field, which quality checks must be performed and which types of aggregations and transformations must be applied. It internally translates the XML configuration into a direct acyclic computational graph that encodes the dependencies of the computational tasks to be performed. This model allows the library to easily take advantage of parallelization at the thread level and the overall flexibility allow us to develop generic data quality analysis pipelines that could also be reused in other applications.

It is now well established that the quiet Sun contains in total more magnetic flux than active regions and represents an important reservoir of magnetic energy. But the nature and evolution of these fields remain largely unknown. We investigate the solar-cycle and center-to-limb variations of magnetic-flux structures at small scales in internetwork regions of the quiet Sun. We used Hinode SOT/SP data from the irradiance program between 2008 and 2016. Maps of the magnetic-flux density are derived from the center-of gravity method applied to the FeI 630.15 nm and FeI 630.25 nm lines. To correct the maps from the instrumental smearing, we applied a deconvolution method based on a principal component analysis of the line profiles and on a Richardson-Lucy deconvolution of their coefficients. We then performed a spectral analysis of the spatial fluctuations of the magnetic-flux density in 10'' x 10'' internetwork regions spanning a wide range of latitudes. At low and mid latitudes the power spectra do not vary significantly with the solar cycle. However at solar maximum for one scan in the activity belt showing an enhanced network, a marginal increase in the power of the magnetic fluctuations is observed at granular and larger scales in the internetwork. At high latitudes, we observe variations at granular and larger scales where the power decreases at solar maximum. At all the latitudes the power of the magnetic fluctuations at scales smaller than 0.5''remain constant throughout the solar cycle. Our results favor a small-scale dynamo that operates in the internetwork, but they show that the global dynamo also contributes to the internetwork fields.

The largest anisotropy in the cosmic microwave background (CMB) is the 3 mK kinematic dipole reflecting our motion with respect to the CMB frame and pointed in the direction $(l, b) = (264^\circ, +48^\circ)$ in Galactic coordinates. We introduce the concept of the ring of attraction (RA), which is orthogonal to the axis of the kinematic dipole. These directions overlap with the zone of percolation for the kinematic dipole, where its amplitude almost vanishes. We show that along this ring are oriented the directions of the dipole modulation of the CMB, and positions of the peaks responsible for generation of parity asymmetry. This coincidence is peculiar at around the 3 sigma level. We analyzed the "interaction" of low multipoles of the CMB with RA and showed that for odd modes there is a sequence of peaks in the RA direction. These peaks correlate with each other for different multipoles and result in mutual amplification of the odd $\ell$ signal for the first 30 multipoles. Our method sheds new light on the nature of parity asymmetry. It consists of the deficit of symmetrically located and equal in amplitude peaks in the CMB map in comparison with asymmetric peaks.

In a search for eclipsing white dwarfs using the Zwicky Transient Facility lightcurves, we identified a deep eclipsing white dwarf with a dark, substellar companion. The lack of an infrared excess and an orbital period of 10 hours made this a potential exoplanet candidate. We obtained high-speed photometry and radial velocity measurements to characterize the system. The white dwarf has a mass of $0.50\pm0.02\,\mathrm{M_{\odot}}$ and a temperature of $10900\pm200\,$K. The companion has a mass of $0.059\pm0.004\,\mathrm{M_{\odot}}$ and a small radius of $0.0783\pm0.0013\,\mathrm{R_{\odot}}$. It is one of the smallest transiting brown dwarfs known and likely old, $\gtrsim 8\,$Gyr. The ZTF discovery efficiency of substellar objects transiting white dwarfs is limited by the number of epochs and as ZTF continues to collect data we expect to find more of these systems. This will allow us to measure period and mass distributions and allows us to understand the formation channels of white dwarfs with substellar companions.

In the past few years, new observations of neutron stars and neutron-star mergers have provided a wealth of data that allow one to constrain the equation of state of nuclear matter at densities above nuclear saturation density. However, most observations were based on neutron stars with masses of about 1.4 solar masses, probing densities up to $\sim$ 3-4 times the nuclear saturation density. Even higher densities are probed inside massive neutron stars such as PSR J0740+6620. Very recently, new radio observations provided an update to the mass estimate for PSR J0740+6620 and X-ray observations by the NICER and XMM telescopes constrained its radius. Based on these new measurements, we revisit our previous nuclear-physics multi-messenger astrophysics constraints and derive updated constraints on the equation of state describing the neutron-star interior. By combining astrophysical observations of two radio pulsars, two NICER measurements, the two gravitational-wave detections GW170817 and GW190425, detailed modeling of the kilonova AT2017gfo, as well as the gamma-ray burst GRB170817A, we are able to estimate the radius of a typical 1.4-solar mass neutron star to be $11.94^{+0.76}_{-0.87} \rm{km}$ at 90\% confidence. Our analysis allows us to revisit the upper bound on the maximum mass of neutron stars and disfavours the presence of a strong first-order phase transition from nuclear matter to exotic forms of matter, such as quark matter, inside neutron stars.

Although micro-lensing of macro-lensed quasars and supernovae provides unique opportunities for several kinds of investigations, it can add unwanted and sometimes substantial noise. While micro-lensing flux anomalies may be safely ignored for some observations, they severely limit others. "Worst-case" estimates can inform the decision whether or not to undertake an extensive examination of micro-lensing scenarios. Here, we report "worst-case" micro-lensing uncertainties for point sources lensed by singular isothermal potentials, parameterized by a convergence equal to the shear and by the stellar fraction. The results can be straightforwardly applied to non-isothermal potentials utilizing the mass sheet degeneracy. We use micro-lensing maps to compute fluctuations in image micro-magnifications and estimate the stellar fraction at which the fluctuations are greatest for a given convergence. We find that the worst-case fluctuations happen at a stellar fraction $\kappa_\star=\frac{1}{|\mu_{macro}|}$. For macro-minima, fluctuations in both magnification and demagnification appear to be bounded ($1.5>\Delta m>-1.3$, where $\Delta m$ is magnitude relative to the average macro-magnification). Magnifications for macro-saddles are bounded as well ($\Delta m > -1.7$). In contrast, demagnifications for macro-saddles appear to have unbounded fluctuations as $1/\mu_{macro}\rightarrow0$ and $\kappa_\star\rightarrow0$.

Sub-GeV mass dark matter particles whose collisions with nuclei would not deposit sufficient energy to be detected, could instead be revealed through their interaction with electrons. Analyses of data from direct detection experiments usually require assuming a local dark matter halo velocity distribution. In the halo-independent analysis method, properties of this distribution are instead inferred from direct dark matter detection data, which allows then to compare different data without making any assumption on the uncertain local dark halo. This method has so far been developed for and applied to dark matter scattering off nuclei. Here we demonstrate how this analysis can be applied to scattering off electrons.

Axion-like particles (ALPs) are pseudo-scalar particles that are candidates for ultralight dark matter. ALPs interact with photons slightly and cause the rotational oscillation of linear polarization. DANCE searches for ALP dark matter by enhancing the rotational oscillation in a bow-tie ring cavity. The signal to noise ratio of DANCE can be improved by long-term observation, and we are planning a year-long observation for the final DANCE. In this document, I will report on the control systems of the ring cavity we developed for the future long-term observation.

The cosmological constant puzzle, traditionally viewed as a naturalness problem, is evidently nullified by the $S$-matrix formulation of quantum gravity/string theory. We point out an implication of this fact for another naturalness puzzle, the Hierarchy Problem between the weak and Planck scales. By eliminating the landscape of de Sitter vacua and eternal inflation, the $S$-matrix formulation exhibits an obvious tension with the explanations based on anthropic selection or cosmological relaxation of the Higgs mass. This sharpens the Hierarchy Problem in a profound way. On one hand, it strengthens the case for explanations based on new physics not far from the weak scale. At the same time, it opens up a question, whether instead the hierarchy is imposed by the $S$-matrix consistency between the Standard Model and gravity.

Galaxy clusters are filled with hot, diffuse X-ray emitting plasma, with a stochastically tangled magnetic field whose energy is close to equipartition with the energy of the turbulent motions \cite{zweibel1997, Vacca}. In the cluster cores, the temperatures remain anomalously high compared to what might be expected considering that the radiative cooling time is short relative to the Hubble time \cite{cowie1977,fabian1994}. While feedback from the central active galactic nuclei (AGN) \cite{fabian2012,birzan2012,churazov2000} is believed to provide most of the heating, there has been a long debate as to whether conduction of heat from the bulk to the core can help the core to reach the observed temperatures \cite{narayan2001,ruszkowski2002,kunz2011}, given the presence of tangled magnetic fields. Interestingly, evidence of very sharp temperature gradients in structures like cold fronts implies a high degree of suppression of thermal conduction \cite{markevitch2007}. To address the problem of thermal conduction in a magnetized and turbulent plasma, we have created a replica of such a system in a laser laboratory experiment. Our data show a reduction of local heat transport by two orders of magnitude or more, leading to strong temperature variations on small spatial scales, as is seen in cluster plasmas \cite{markevitch2003}.

It was recently shown, that in a class of tensor-multi-scalar theories of gravity with a nontrivial target space metric, there exist scalarized neutron star solutions. An important property of these compact objects is that the scalar charge is zero and therefore, the binary pulsar experiments can not impose constraints based on the absence of scalar dipole radiation. Moreover, the structure of the solutions is very complicated. For a fixed central energy density up to three neutron star solutions can exist -- one general relativistic and two scalarized, that is quite different from the scalarization in other alternative theories of gravity. In the present paper we address the stability of these solutions using two independent approaches -- solving the linearized radial perturbation equations and performing nonlinear simulations in spherical symmetry. The results show that the change of stability occurs at the maximum mass models and all solutions before that point are stable. This leads to the interesting consequence that there exists a stable part of the scalarized branch close to the bifurcation point where the mass of the star increases with the decrease of the central energy density.

Since particles obey wave equations, in general one is not free to postulate that particles move on the geodesics associated with test particles. Rather, for this to be the case one has to be able to derive such behavior starting from the equations of motion that the particles obey, and to do this for either massless or massive particles one can employ the eikonal approximation. While for massive particles one does obtain standard geodesic behavior this way, for a conformally coupled massless scalar field the eikonal approximation only leads to geodesic behavior if the Ricci scalar is zero. Similarly, for the propagation of the light waves associated with the conformal invariant Maxwell equations geodesic behavior only holds if the Ricci tensor is zero. While for practical purposes such terms might only be of relevance in regions of high curvature, the point of this paper is only to establish their presence in principle. Thus in principle the standard null-geodesic-based gravitational bending formula and the behavior of light rays in cosmology are in need of modification in regions with high enough curvature. We show how to appropriately modify the geodesic equations in such situations. We show that relativistic eikonalization has an intrinsic light-front structure, and show that eikonalization in a theory with local conformal symmetry leads to trajectories that are only globally conformally symmetric. The modifications to geodesics that we find lead to the propagation of massless particles off the light cone. This is a curved space reflection of the fact that when light travels through a refractive medium in flat spacetime its velocity is modified from its free flat spacetime value. In the presence of gravity spacetime itself acts as a medium, and this medium can then take light waves off the light cone.

We describe a new phenomenon in quantum cosmology: self-organised localisation. When the fundamental parameters of a theory are functions of a scalar field subject to large fluctuations during inflation, quantum phase transitions can act as dynamical attractors. As a result, the theory parameters are probabilistically localised around the critical value and the Universe finds itself at the edge of a phase transition. We illustrate how self-organised localisation could account for the observed near-criticality of the Higgs self-coupling, the naturalness of the Higgs mass, or the smallness of the cosmological constant.

Context: The analysis of the thermal part of velocity distribution functions (VDF) is fundamentally important for understanding the kinetic physics that governs the evolution and dynamics of space plasmas. However, calculating the proton core, beam and alpha-particle parameters for large data sets of VDFs is a time consuming and computationally demanding process that always requires supervision by a human expert. Aims: We developed a machine learning tool that can extract proton core, beam and alpha-particle parameters using images (2-D grid consisting pixel values) of VDFs. Methods: A database of synthetic VDFs is generated, which is used to train a convolutional neural network that infers bulk speed, thermal speed and density for all three particle populations. We generate a separate test data set of synthetic VDFs that we use to compare and quantify the predictive power of the neural network and a fitting algorithm. Results: The neural network achieves significantly smaller root-mean-square errors to infer proton core, beam and alpha-particle parameters than a traditional fitting algorithm. Conclusion: The developed machine learning tool has the potential to revolutionize the processing of particle measurements since it allows the computation of more accurate particle parameters than previously used fitting procedures.

In this work, a bumblebee field is adopted in order to generate cosmological anisotropies. For that purpose, we assume a Bianchi I cosmology, as the background geometry, and a bumblebee field coupled to it. Bumblebee models are examples of a mechanism for the Lorentz symmetry violation by assuming a nonzero vacuum expectation value for the bumblebee field. When coupled to the Bianchi I geometry, which is not in agreement with a cosmological principle, the bumblebee field plays the role of a source of anisotropies and produces a preferred axis. Thus, a fraction of the cosmic anisotropies would come from the Lorentz symmetry violation. In the last part of the article, we try to assume an upper bound on the bumblebee field using the quadrupole and octopole moments of the cosmic microwave background radiation.

We review a testable, the axion quark nugget (AQN) model outside of the standard WIMP paradigm. The model was originally invented to explain the observed similarity between the dark and the visible components, $\Omega_{\rm DM}\approx \Omega_{\rm visible}$ in a natural way as both types of matter are formed during the same QCD transition and proportional to the same dimensional fundamental parameter of the system, $\Lambda_{\rm QCD}$. In this framework the baryogenesis is actually a charge segregation (rather than charge generation) process which is operational due to the $\cal{CP}$-odd axion field,while the global baryon number of the Universe remains zero. The nuggets and anti-nuggets are strongly interacting but macroscopically large objects with approximately nuclear density. We overview several specific recent applications of this framework. First, we discuss the "solar corona mystery" when the so-called nanoflares are identified with the AQN annihilation events in corona. Secondly, we review a proposal that the recently observed by the Telescope Array puzzling events is a result of the annihilation events of the AQNs under thunderstorm. Finally, we overview a broadband strategy which could lead to the discovery the AQN-induced axions representing the heart of the construction.