Papers for Friday, Jul 30 2021

We present the results of our search for variable stars using the long-term Las Cumbres Observatory (LCO) monitoring of white dwarf ZTF J0139+5245 with the two 1.0-m telescope nodes located at McDonald Observatory using the Sinistro imaging instrument. In this search, we find 38 variable sources, of which 27 are newly discovered or newly classified (71%) based on comparisons with previously published catalogs, thereby increasing the number of detections in the field-of-view under consideration by a factor of $\approx$ 2.5. We find that the improved photometric precision per-exposure due to longer exposure time for LCO images combined with the greater time-sampling of LCO photometry enables us to increase the total number of detections in this field-of-view. Each LCO image covers a field-of-view of $26' \times 26'$ and observes a region close to the Galactic plane ($b = -9.4^\circ$) abundant in stars with an average stellar density of $\approx 8$ arcmin$^{-2}$. We perform aperture photometry and Fourier analysis on over 2000 stars across 1560 LCO images spanning 537 days to find 28 candidate BY Draconis variables, 3 candidate eclipsing binaries of type EA, and 7 candidate eclipsing binaries of type EW. In assigning preliminary classifications to our detections, we demonstrate the applicability of the Gaia color-magnitude diagram (CMD) as a powerful classification tool for variable star studies.

'Downwards drifting' structures observed in fast radio bursts (FRBs) could naturally arise from a screen occurring after a small initial dispersion measure region. The screen imprints temporally sharp but broadband structure on the pulse that has already been dispersed, and the `structure-maximising' bulk of the dispersion measure is acquired further along the path to observer. If so, scaling of the drift rate of repeating FRBs with frequency suggests that the emission is beamed -- out of our line of sight -- with the circum-burst plasma deflecting the beam towards us. This in turn explains the observed limited, variable bandwidths of bursts despite broadband nature of the underlying emission. We summarise the geometric constraints on this simple model delivering much of the observed FRB spectro-temporal phenomenology, present generic predictions and discuss possible nature of modulation.

The innermost regions of accretion disks around black holes are strongly irradiated by X-rays that are emitted from a highly variable, compact corona, in the immediate vicinity of the black hole. The X-rays that are seen reflected from the disk and the time delays, as variations in the X-ray emission echo or reverberate off the disk provide a view of the environment just outside the event horizon. I Zwicky 1 (I Zw 1), is a nearby narrow line Seyfert 1 galaxy. Previous studies of the reverberation of X-rays from its accretion disk revealed that the corona is composed of two components; an extended, slowly varying component over the surface of the inner accretion disk, and a collimated core, with luminosity fluctuations propagating upwards from its base, which dominates the more rapid variability. Here we report observations of X-ray flares emitted from around the supermassive black hole in I Zw 1. X-ray reflection from the accretion disk is detected through a relativistically broadened iron K line and Compton hump in the X-ray emission spectrum. Analysis of the X-ray flares reveals short flashes of photons consistent with the re-emergence of emission from behind the black hole. The energy shifts of these photons identify their origins from different parts of the disk. These are photons that reverberate off the far side of the disk and bent around the black hole and magnified by the strong gravitational field. Observing photons bent around the black hole confirms a key prediction of General Relativity.

The Submillimetre Common User Bolometer Array 2 (SCUBA-2) is the James Clerk Maxwell Telescope's continuum imager, operating simultaneously at 450 and 850~$\mu$m. SCUBA-2 was commissioned in 2009--2011 and since that time, regular observations of point-like standard sources have been performed whenever the instrument is in use. Expanding the calibrator observation sample by an order of magnitude compared to previous work, in this paper we derive updated opacity relations at each wavelength for a new atmospheric-extinction correction, analyze the Flux-Conversion Factors (FCFs) used to convert instrumental units to physical flux units as a function of date and observation time, present information on the beam profiles for each wavelength, and update secondary-calibrator source fluxes. Between 07:00 and 17:00 UTC, the portion of the night that is most stable to temperature gradients that cause dish deformation, the total-flux uncertainty and the peak-flux uncertainty measured at 450~$\mu$m are found to be 14\% and 17\%, respectively. Measured at 850~$\mu$m, the total-flux and peak-flux uncertainties are 6\%, and 7\%, respectively. The analysis presented in this work is applicable to all SCUBA-2 projects observed since 2011.

The fraction of stars which are in binaries or triples at the time of stellar death and the fraction of these systems which survive the supernova (SN) explosion are crucial constraints for evolution models and predictions for gravitational wave source populations. These fractions are also subject to direct observational determination. Here we search 10 supernova remnants (SNR) containing compact objects with proper motions for unbound binaries or triples using Gaia EDR3 and new statistical methods and tests for false positives. We confirm the one known example of an unbound binary, HD 37424 in G180.0-01.7, and find no other examples. Combining this with our previous searches for bound and unbound binaries, and assuming no bias in favor of finding interacting binaries, we find that 72.0% (52.2%-86.4%, 90% confidence) of SN producing neutron stars are not binaries at the time of explosion, 13.9% (5.4%-27.2%) produce bound binaries and 12.5% (2.8%-31.3%) produce unbound binaries. With a strong bias in favor of finding interacting binaries, the medians shift to 76.0% were not binaries at death, 9.5% leave bound and 13.2% leave unbound binaries. Of explosions that do not leave binaries, <18.9% can be fully unbound triples. These limits are conservatively for M>5Msun stars, although the mass limits for individual systems are significantly stronger. At birth, the progenitor of PSR J0538+2817 was probably a 13-19Msun star, and at the time of explosion it was probably a Roche limited, partially stripped star transferring mass to HD 37424 and then producing a Type IIL or IIb supernova.

[Full abstract in the paper] In recent years, protoplanetary disks with spiral structures have been detected in scattered light, millimeter continuum, and CO gas emission. The mechanisms causing these structures are still under debate. A popular scenario to drive the spiral arms is the one of a planet perturbing the material in the disk. However, if the disk is massive, gravitational instability is usually the favored explanation. Multiwavelength studies could be helpful to distinguish between the two scenarios. So far, only a handful of disks with spiral arms have been observed in both scattered light and millimeter continuum. We aim to perform an in-depth characterization of the protoplanetary disk morphology around WaOph 6 analyzing data obtained at different wavelengths, as well as to investigate the origin of the spiral features in the disk. We present the first near-infrared polarimetric observations of WaOph 6 obtained with SPHERE at the VLT and compare them to archival millimeter continuum ALMA observations. We traced the spiral features in both data sets and estimated the respective pitch angles. We discuss the different scenarios that can give rise to the spiral arms in WaOph 6. We tested the planetary perturber hypothesis by performing hydrodynamical and radiative transfer simulations to compare them with scattered light and millimeter continuum observations.

We address the nature of the giant clumps in high-z galaxies that undergo Violent Disc Instability, attempting to distinguish between long-lived clumps that migrate inward and short-lived clumps that disrupt by feedback. We study the evolution of clumps as they migrate through the disc using an analytic model tested by simulations and confront theory with CANDELS observations. The clump ``bathtub" model, which considers gas and stellar gain and loss, is characterized by four parameters: the accretion efficiency, the star-formation-rate (SFR) efficiency, and the outflow mass-loading factors for gas and stars. The relevant timescales are all comparable to the migration time, two-three orbital times. A clump differs from a galaxy by the internal dependence of the accretion rate on the varying clump mass. The analytic solution, involving exponential growing and decaying modes, reveals a main evolution phase during the migration, where the SFR and gas mass are constant and the stellar mass is rising linearly with time. This makes the inverse of the specific SFR an observable proxy for clump age. Later, the masses and SFR approach an exponential growth with a constant specific SFR, but this phase is hypothetical as the clump disappears in the galaxy center. The model matches simulations with different, moderate feedback, both in isolated and cosmological settings. The observed clumps agree with our theoretical predictions, indicating that the massive clumps are long-lived and migrating. A non-trivial challenge is to model feedback that is non-disruptive in massive clumps but suppresses SFR to match the galactic stellar-to-halo mass ratio.

IRAS~04158+2805 has long been thought to be a very low mass T-Tauri star (VLMS) surrounded by a nearly edge-on, extremely large disc. Recent observations revealed that this source hosts a binary surrounded by an extended circumbinary disc with a central dust cavity. In this paper, we combine ALMA multi-wavelength observations of continuum and $^{12}$CO line emission, with H$\alpha$ imaging and Keck astrometric measures of the binary to develop a coherent dynamical model of this system. The system features an azimuthal asymmetry detected at the western edge of the cavity in Band~7 observations and a wiggling outflow. Dust emission in ALMA Band 4 from the proximity of the individual stars suggests the presence of marginally resolved circumstellar discs. We estimate the binary orbital parameters from the measured arc of the orbit from Keck and ALMA astrometry. We further constrain these estimates using considerations from binary-disc interaction theory. We finally perform three SPH gas + dust simulations based on the theoretical constraints; we post-process the hydrodynamic output using radiative transfer Monte Carlo methods and directly compare the models with observations. Our results suggest that a highly eccentric $e\sim 0.5\textrm{--}0.7$ equal mass binary, with a semi-major axis of $\sim 55$ au, and small/moderate orbital plane vs. circumbinary disc inclination $\theta\lesssim 30^\circ$ provides a good match with observations. A dust mass of $\sim 1.5\times 10^{-4} {\rm M_\odot}$ best reproduces the flux in Band 7 continuum observations. Synthetic CO line emission maps qualitatively capture both the emission from the central region and the non-Keplerian nature of the gas motion in the binary proximity.

The spatial distribution of the HI gas in galaxies holds important clues on the physical processes that shape the structure and dynamics of the interstellar medium (ISM). In this work, we quantify the structure of the HI gas in a sample of 33 nearby galaxies taken from the THINGS Survey using the delta-variance spectrum. The THINGS galaxies display a large diversity in their spectra, however, there are a number of recurrent features. In many galaxies, we observe a bump in the spectrum on scales of a few to several hundred pc. We find the characteristic scales associated with the bump to be correlated with galactic SFR for values of the SFR > 0.5 M$_{sol}$ yr$^{-1}$ and also with the median size of the HI shells detected in those galaxies. On larger scales, we observe the existence of two self-similar regimes. The first one, on intermediate scales is shallow and the power law that describes this regime has an exponent in the range [0.1-1] with a mean value of 0.55 which is compatible with the density field being generated by supersonic turbulence in the cold phase of the HI gas. The second power law is steeper, with a range of exponents between [0.5-1.5] and a mean value of 1.5. These values are associated with subsonic turbulence which is characteristic of the warm phase of the HI gas. The spatial scale at which the transition between the two regimes occurs is found to be $\approx 0.5 R_{25}$ which is similar to the size of the molecular disk in the THINGS galaxies. Overall, our results suggest that on scales < $0.5 R_{25}$, the structure of the ISM is affected by the effects of supernova explosions. On larger scales (> 0.5 $R_{25}$), stellar feedback has no significant impact, and the structure of the ISM is determined by large scale processes that govern the dynamics of the gas in the warm neutral medium such as the flaring of the HI disk and the effects of ram pressure stripping.

Hydrogen emissions of RR Lyrae variables are the imprints of shock waves traveling through their atmospheres. We develop a pattern recognition algorithm, which is then applied to single-epoch spectra of SDSS and LAMOST. These two spectroscopic surveys covered $\sim$ 10,000 photometrically confirmed RR Lyrae stars. We discovered in total 127 RR Lyrae stars with blueshifted Balmer emission feature, including 103 fundamental mode (RRab), 20 first-overtone (RRc), 3 double-mode (RRd), and 1 Blazhko type (temporary classification for RR Lyrae stars with strong Blazhko modulation in Catalina sky survey that cannot be characterized) RR Lyrae variable. This forms the largest database to date of the properties of hydrogen emission in RR Lyrae variables. Based on ZTF DR5, we carried out a detailed light-curve analysis for the Blazhko type RR Lyrae star with hydrogen emission of long-term modulations. We characterize the Blazhko type RR Lyrae star as an RRab and point out a possible Blazhko period. Finally, we set up simulations on mock spectra to test the performance of our algorithm and on the real observational strategy to investigate the occurrence of the "first apparition".

Quantitative morphologies, such as asymmetry and concentration, have long been used as an effective way to assess the distribution of galaxy starlight in large samples. Application of such quantitative indicators to other data products could provide a tool capable of capturing the 2-dimensional distribution of a range of galactic properties, such as stellar mass or star-formation rate maps. In this work, we utilize galaxies from the Illustris and IllustrisTNG simulations to assess the applicability of concentration and asymmetry indicators to the stellar mass distribution in galaxies. Specifically, we test whether the intrinsic values of concentration and asymmetry (measured directly from the simulation stellar mass particle maps) are recovered after the application of measurement uncertainty and a point spread function (PSF). We find that random noise has a non-negligible systematic effect on asymmetry that scales inversely with signal-to-noise, particularly at signal-to-noise less than 100. We evaluate different methods to correct for the noise contribution to asymmetry at very low signal-to-noise, where previous studies have been unable to explore due to systematics. We present algebraic corrections for noise and resolution to recover the intrinsic morphology parameters. Using Illustris as a comparison dataset, we evaluate the robustness of these fits in the presence of a different physics model, and confirm these correction methods can be applied to other datasets. Lastly, we provide estimations for the uncertainty on different correction methods at varying signal-to-noise and resolution regimes.

Giant stars as known exoplanet hosts are relatively rare due to the potential challenges in acquiring precision radial velocities and the small predicted transit depths. However, these giant host stars are also some of the brightest in the sky and so enable high signal-to-noise follow-up measurements. Here we report on new observations of the bright (V ~ 3.3) giant star $\iota$ Draconis ($\iota$ Dra), known to host a planet in a highly eccentric ~511 day period orbit. TESS observations of the star over 137 days reveal asteroseismic signatures, allowing us to constrain the stellar radius, mass, and age to ~2%, ~6%, and ~28%, respectively. We present the results of continued radial velocity monitoring of the star using the Automated Planet Finder over several orbits of the planet. We provide more precise planet parameters of the known planet and, through the combination of our radial velocity measurements with Hipparcos and Gaia astrometry, we discover an additional long-period companion with an orbital period of ~$68^{+60}_{-36}$ years. Mass predictions from our analysis place this sub-stellar companion on the border of the planet and brown dwarf regimes. The bright nature of the star combined with the revised orbital architecture of the system provides an opportunity to study planetary orbital dynamics that evolve as the star moves into the giant phase of its evolution.

Triple-star systems exhibit a phenomenon known as the Triple Evolution Dynamical Instability (TEDI), in which mass loss in evolving triples triggers short-term dynamical instabilities, potentially leading to collisions of stars, exchanges, and ejections. Previous work has shown that the TEDI is an important pathway to head-on stellar collisions in the Galaxy, significantly exceeding the rate of collisions due to random encounters in globular clusters. Here, we revisit the TEDI evolutionary pathway using state-of-the-art population synthesis methods that self-consistently take into account stellar evolution and binary interactions, as well as gravitational dynamics and perturbations from passing stars in the field. We find Galactic TEDI-induced collision rates on the order of 1e-4/yr, consistent with previous studies which were based on more simplified methods. The majority of TEDI-induced collisions involve main sequence stars, potentially producing blue straggler stars. Collisions are also possible involving more evolved stars, potentially producing eccentric post-common-envelope systems, and white dwarfs collisions leading to Type Ia supernovae (although the latter with a negligible contribution to the Galactic rate). In our simulations, the TEDI is not only triggered by adiabatic wind mass loss, but also by Roche lobe overflow in the inner binary: when the donor star becomes less massive than the accretor, the inner binary orbit widens, triggering triple dynamical instability. We find that collision rates are increased by ~17% when fly-bys in the field are taken into account. In addition to collisions, we find that the TEDI produces ~1e-4/yr of unbound stars, although none with escape speeds in excess of 1e3 km/s.

Waiting time distributions allow us to distinguish at least three different types of dynamical systems, such as (i) linear random processes (with no memory); (ii) nonlinear, avalanche-type, nonstationary Poisson processes (with memory during the exponential growth of the avalanche rise time); and (iii) chaotic systems in the state of a nonlinear limit cycle (with memory during the oscillatory phase). We describe the temporal evolution of the flare rate $\lambda(t) \propto t^p$ with a polynomial function, which allows us to distinguish linear ($p \approx 1$) from nonlinear ($p \gapprox 2$) events. The power law slopes $\alpha$ of observed waiting times (with full solar cycle coverage) cover a range of $\alpha=2.1-2.4$, which agrees well with our prediction of $\alpha = 2.0+1/p = 2.3-2.5$. The memory time can also be defined with the time evolution of the logistic equation, for which we find a relationship between the nonlinear growth time $\tau_G = \tau_{rise}/(4p)$ and the nonlinearity index $p$. We find a nonlinear evolution for most events, in particular for the clustering of solar flares ($p=2.2\pm0.1$), partially occulted flare events ($p=1.8\pm0.2$), and the solar dynamo ($p=2.8\pm0.5$). The Sun exhibits memory on time scales of $\lapprox$2 hours to 3 days (for solar flare clustering), 6 to 23 days (for partially occulted flare events), and 1.5 month to 1 year (for the rise time of the solar dynamo).

We study the relationship between the H{\sc i} specific angular momentum (j$_{\rm g}$) and the H{\sc i} mass (M$_{\rm g}$) for a sample of galaxies with well measured H{\sc i} rotation curves. We find that the relation is well described by an unbroken power law \jg $\propto$ \mg$^{\alpha}$ over the entire mass range (10$^{7}$-10$^{10.5}$ M$_{\odot}$), with $\alpha = 0.89 \pm 0.05$ (scatter 0.18 dex). This is in reasonable agreement with models which assume that evolutionary processes maintain H{\sc i} disks in a marginally stable state. The slope we observe is also significantly different from both the $j \propto M^{2/3}$ relation expected for dark matter haloes from tidal torquing models and the observed slope of the specific angular momentum-mass relation for the stellar component of disk galaxies. Our sample includes two H{\sc i}-bearing ultra diffuse galaxies, and we find that their angular momentum follows the same relation as other galaxies. The only discrepant galaxies in our sample are early-type galaxies with large rotating H{\sc i} disks which are found to have significantly higher angular momentum than expected from the power law relation. The H{\sc i} disks of all these early-type galaxies are misaligned or counter-rotating with respect to the stellar disks, consistent with the gas being recently accreted. We speculate that late stage wet mergers, as well as cold flows play a dominant role in determining the kinematics of the baryonic component of galaxies as suggested by recent numerical simulations.

Observations of small-scale brightenings in the low solar atmosphere can provide valuable constraints on possible heating/heat-transport mechanisms. We present a method for the detection and analysis of brightenings and demonstrate its application to IRIS EUV time-series imagery. The method uses band-pass filtering, adaptive thresholding and centroid tracking, and records an event's position, duration, and total/maximum brightness. Area, brightness, and position are also recorded as functions of time throughout their lifetime. Detected brightenings can fragment or merge over time, thus the number of distinct regions constituting a brightening event is recorded over time and the maximum number of regions are recorded as a simple measure of an event's coherence/complexity. A test is made on a synthetic datacube composed of a static background based on IRIS data, Poisson noise and $\approx10^4$ randomly-distributed, moving, small-scale Gaussian brightenings. Maximum brightness, total brightness, area, and duration follow power-law distributions and the results show the range over which the method can extract information. The recorded maximum brightness is a reliable measure for the brightest and most accurately detected events with an error of 6%. Area, duration, and speed are generally underestimated by 15% with an uncertainty of 20-30%. Total brightness is underestimated by 30% with an uncertainty of 30%. Applying this method to real IRIS QS data spanning 19 minutes over a 54.40"$\times$55.23" FOV yields 2997 detections. 1340 of these either remain un-fragmented or fragment to two distinct regions at least once during their lifetime equating to an event density of $3.96\times10^{-4}$arcsec$^{-2}$s$^{-1}$. The method will be used for a future large-scale statistical analysis of several QS data sets from IRIS, other EUV imagers, as well as H-$\alpha$ and visible photospheric imagery.

Upcoming NASA astrophysics missions such as the James Webb Space Telescope will search for signs of life on planets transiting nearby stars. Doing so will require co-adding dozens of transmission spectra to build up sufficient signal to noise while simultaneously accounting for challenging systematic effects such as surface/weather variability, atmospheric refraction, and stellar activity. To determine the magnitude and impacts of both stellar and planet variability on measured transmission spectra, we must assess the feasibility of stacking multiple transmission spectra of exo-Earths around their host stars. Using our own solar system, we can determine if current methodologies are sufficient to detect signs of life in Earth's atmosphere and measure the abundance of habitability indicators, such as H2O and CO2, and biosignature pairs, such as O2 and CH4. We assess the impact on transmission spectra of Earth transiting across the Sun from solar and planetary variability and identify remaining unknowns for understanding exoplanet transmission spectra. We conclude that a satellite observing Earth transits across the Sun from beyond L2 is necessary to address these long-standing concerns about the reliability of co-adding planet spectra at UV, optical, and infrared wavelengths from multiple transits in the face of relatively large astrophysical systematics.

ANTARES is currently the largest undersea neutrino telescope, located in the Mediterranean Sea and taking data since 2007. It consists of a 3D array of photo sensors, instrumenting about 10Mt of seawater to detect Cherenkov light induced by secondary particles from neutrino interactions. The event reconstruction and background discrimination is challenging and machine-learning techniques are explored to improve the performance. In this contribution, two case studies using deep convolutional neural networks are presented. In the first one, this approach is used to improve the direction reconstruction of low-energy single-line events, for which the reconstruction of the azimuth angle of the incoming neutrino is particularly difficult. We observe a promising improvement in resolution over classical reconstruction techniques and expect to at least double our sensitivity in the low-energy range, important for dark matter searches. The second study employs deep learning to reconstruct the visible energy of neutrino interactions of all flavors and for the multi-line setup of the full detector.

The origin of the high-energy gamma-ray emission from the Milky Way center is still unclear and debated because of the impact of systematics afflicting the measurements from current experiments. Several theories and phenomenological models attempt to explain the intricate panorama. The presence of a PeVatron in the Central Molecular Zone or in its vicinity, the contribution of the hard-component of the diffuse gamma-ray emission, and dark matter annihilation scenario are among the most promising mechanisms for describing the observed excess. The development of increasingly precise models able to reproduce the measured gamma-ray emission is the challenge for the scientific community in view of the next generation telescopes. A detailed treatment of phenomenological models for the dubbed Cosmic Rays Sea (CR-sea) characterized by different configurations is scrutinized in comparison with the observed spectrum in the inner Galaxy, using DRAGON and GAMMASKY codes.

The orbital periods of most eclipsing cataclysmic binaries are not undergoing linear secular decreases of order a few parts per billion as expected from simple theory. Instead, they show several parts per million increases and decreases on timescales of years to decades, ascribed to magnetic effects in their donors, triple companions, or both. To directly test the triple companion hypothesis, we carried out a speckle imaging survey of six of the nearest and brightest cataclysmic variables. We found no main sequence companions earlier than spectral types M4V in the separation range 0.02" - 1.2", corresponding to projected linear separations of 2 - 100 AU, and periods of 3 - 1000 years. We conclude that main sequence triple companions to CVs are not very common, but cannot rule out the presence of the faintest M dwarfs or close brown dwarf companions.

Galaxies and groups of galaxies exist in dark-matter halos filled with diffuse gas. The diffuse gas represents up to 80\% of the mass in luminous matter within the halos (1,2), and is difficult to detect because of its low density (particle number densities of $\lesssim10^{-4}$\,cm$^{-3}$) and high temperature (mostly greater than $10^{6}$\,K). The spatial distribution and total mass of this material determines, and is influenced by, the evolution of galaxies and galaxy groups (3-5). Existing observational constraints on these quantities are limited by sensitivity, and the necessity to accurately model the ionization fraction, metal content, and pressure of the gas (6-8). Here we report the detection of diffuse gas associated with nearby galaxies using the dispersion measures (DMs) of extragalactic fast radio bursts (FRBs). FRB DMs provide direct measurements of the total ionized-gas contents along their sightlines. Out of a sample of 474 distant FRBs from the CHIME/FRB Catalog 1 (9), we identify a sample of 24 events that likely intersect the dark-matter halos of galaxies in the local Universe ($<40$\,Mpc). This subset of FRBs has an excess mean DM of $200\pm100$\,pc\,cm$^{-3}$ over those that do not intersect nearby galaxies. The excess is larger than expected for the diffuse gas surrounding isolated galaxies, but may be explained by additional contributions from gas surrounding galaxy groups, including from the Local Group. This result demonstrates the predicted ability of FRBs to be used as sensitive, model-independent measures of the diffuse-gas contents of dark-matter halos (10-13).

Cosmic-rays interacting with nucleons in the solar atmosphere produce a cascade of particles that give rise to a flux of high-energy neutrinos and gamma-rays. Fermi has observed this gamma-ray flux; however, the associated neutrino flux has escaped observation. In this contribution, we put forward two strategies to detect these neutrinos, which, if seen, would push forward our understanding of the solar atmosphere and provide a new testing ground of neutrino properties. First, we will extend the previous analysis, which used high-energy through-going muon events collected in the years of maximum solar activity and yielded only flux upper limits, to include data taken during the solar minimum from 2018 to 2020. Extending the analysis to the solar minimum is important as the gamma-ray data collected during past solar cycles indicates a possible enhancement in the high-energy neutrino flux. Second, we will incorporate sub-TeV events and include contributions from all neutrino flavors. These will improve our analysis sensitivity since the solar atmospheric spectrum is soft and, due to oscillation, contains significant contributions of all neutrino flavors. As we will present in this contribution, these complementary strategies yield a significant improvement in sensitivity, making substantial progress towards observing this flux.

Our knowledge of the birth mass function of neutron stars and black holes is based on observations of binary systems but the binary evolution likely affects the final mass of the compact object. Gravitational microlensing allows us to detect and measure masses of isolated stellar remnants, which are nearly impossible to obtain with other techniques. Here, we analyze a sample of 4360 gravitational microlensing events detected during the third phase of the OGLE survey. We select a subsample of 87 long-timescale low-blending events. We estimate the masses of lensing objects by combining photometric data from OGLE and proper-motion information from OGLE and Gaia EDR3. We find 35 high-probability dark lenses - white dwarfs, neutron stars, and black holes - which we use to constrain the mass function of isolated stellar remnants. In the range 1-100 M_Sun, occupied by neutron stars and black holes, the remnant mass function is continuous and can be approximated as a power-law with a slope of $0.83^{+0.16}_{-0.18}$ with a tentative evidence against a broad gap between neutron stars and black holes. This slope is slightly flatter than the slope of the mass function of black holes detected by gravitational wave detectors LIGO and Virgo, although both values are consistent with each other within the quoted error bars. The measured slope of the remnant mass function agrees with predictions of some population synthesis models of black hole formation.

Gravitational microlensing may detect dark stellar remnants - black holes or neutron stars - even if they are isolated. However, it is challenging to estimate masses of isolated dark stellar remnants using solely photometric data for microlensing events. A recent analysis of OGLE-III long-timescale microlensing events exhibiting the annual parallax effects claimed that a number of bright events were due to "mass-gap" objects (with masses intermediate between those of neutron stars and black holes). Here, we present a detailed description of the updated and corrected method that can be used to estimate masses of dark stellar remnants detected in microlensing events given the light curve data and the proper motion of the source. We use this updated method, in combination with new proper motions from Gaia EDR3, to revise masses of dark remnant candidates previously found in the OGLE-III data. We demonstrate that masses of "mass-gap" and black hole events identified in the previous work are overestimated and, hence, these objects are most likely main-sequence stars, white dwarfs, or neutron stars.

Centaurs are small bodies orbiting in the giant planet region which were scattered inwards from their source populations beyond Neptune. Some members of the population display comet-like activity during their transition through the solar system, the source of which is not well understood. The range of heliocentric distances where the active Centaurs have been observed, and their median lifetime in the region suggest this activity is neither driven by water-ice sublimation, nor entirely by super-volatiles. Here we present an observational and thermo-dynamical study of 13 Centaurs discovered in the Pan-STARRS1 detection database aimed at identifying and characterizing active objects beyond the orbit of Jupiter. We find no evidence of activity associated with any of our targets at the time of their observations with the Gemini North telescope in 2017 and 2018, or in archival data from 2013 to 2019. Upper limits on the possible volatile and dust production rates from our targets are 1-2 orders of magnitude lower than production rates in some known comets, and are in agreement with values measured for other inactive Centaurs. Our numerical integrations show that the orbits of six of our targets evolved interior to r$\sim$15 AU over the past 100,000 years where several possible processes could trigger sublimation and outgassing, but their apparent inactivity indicates their dust production is either below our detection limit or that the objects are dormant. Only one Centaur in our sample -- 2014 PQ$_{70}$ experienced a sudden decrease in semi-major axis and perihelion distance attributed to the onset of activity for some previously known inactive Centaurs, and therefore is a likely candidate for future outburst. This object should be a target of interest for further observational monitoring.

Magnetohydrodynamic (MHD) waves are candidates for heating the solar chromosphere, although it is still unclear which mode of the wave is dominant in heating. We perform two-dimensional radiative MHD simulation to investigate the propagation of MHD waves in the quiet region of the solar chromosphere. We identify the mode of the shock waves by using the relationship between gas pressure and magnetic pressure across the shock front and calculate their corresponding heating rate through the entropy jump to obtain a quantitative understanding of the wave heating process in the chromosphere. Our result shows that the fast magnetic wave is significant in heating the low-beta chromosphere. The low-beta fast magnetic waves are generated from high-beta fast acoustic waves via mode conversion crossing the equipartition layer. Efficient mode conversion is achieved by large attacking angles between the propagation direction of the shock waves and the chromospheric magnetic field.

The discovery of fast radio burst (FRB) 200428 from galactic SGR J1935+2154 makes it possible to measure rotational changes accompanied by FRBs and to test several FRB models which may be simultaneously associated with glitches. Inspired by this idea, we present order of magnitude calculations to the scenarios proposed. FRB models such as global starquakes, crust fractures and collisions between pulsars and asteroids/comets are discussed. For each mechanism, the maximum glitch sizes are constrained by the isotropic energy release during the X-ray burst and/or the SGR J1935+2154-like radio burst rate. Brief calculations show that, the maximum glitch sizes for different mechanisms differ by order(s) of magnitude. If glitches are detected to be coincident with FRBs from galactic magnetars in the future, glitch behaviors (such as glitch size, rise timescale, the recovery coefficient and spin down rate offset) are promising to serve as criterions to distinguish glitch mechanisms and in turn to constrain FRB models.

The rotation powered pulsar loses angular momentum at a rate of the rotation power divided by the angular velocity $\Omega_*$. This means that the length of the lever arm of the angular momentum extracted by the photons, relativistic particles and wind must be on average $c/\Omega_*$, which is known as the light cylinder radius. Therefore, any deposition of the rotation power within the light cylinder causes insufficient loss of angular momentum. In this paper, we investigate two cases of this type of energy release: polar cap acceleration and Ohmic heating in the magnetospheric current inside the star. As for the first case, the outer magnetosphere beyond the light cylinder is found to compensate the insufficient loss of the angular momentum. We argue that the energy flux coming from the sub-rotating magnetic field lines must be larger than the solid-angle average value, and as a result, an enhanced energy flux emanating beyond the light cylinder is observed in different phases in the light curve from those of emission inside the light cylinder. As for the second case, the stellar surface rotates more slowly than the stellar interior. We find that the way the magnetospheric current closes inside the star is linked to how the angular momentum is transferred inside the star. We obtain numerical solutions which show that the magnetospheric current inside the star spreads over the polar cap magnetic flux embedded in the star in such a way that electromotive force is gained efficiently.

Magnetic reconnection in solar flares can efficiently generate non-thermal electron beams. The accelerated electrons can, in turn, cause radio waves through kinetic instabilities as they propagate through the ambient plasma. We aim at investigating the wave emission caused by fast electron beams (FEBs) with characteristic non-thermal electron velocity distribution functions (EVDFs) generated by kinetic magnetic reconnection: bump-on-tail EVDFs along the separatrices and in the diffusion region, and perpendicular crescent-shaped EVDFs close to the diffusion region. For this sake we utilized 2.5D fully kinetic Particle-In-Cell (PIC) code simulations in this study. We found that: (1) the bump-on-tail EVDFs are unstable to cause electrostatic Langmuir waves via bump-on-tail instabilities and then multiple harmonic transverse waves from the diffusion region and the separatrices of reconnection. (2) The perpendicular crescent-shaped EVDFs, on the other hand, can cause multi-harmonic electromagnetic electron cyclotron waves through electron cyclotron maser instabilities in diffusion region of reconnection. Our results are applicable to diagnose the plasma parameters which control reconnection in solar flares by means of radio waves observations.

In this paper we discuss the design, manufacturing and characterization of the feed horn array of the Strip instrument of the Large Scale Polarization Explorer (LSPE) experiment. Strip is a microwave telescope, operating in the Q- and W-band, for the observation of the polarized emissions from the sky in a large fraction (about 37%) of the Northern hemisphere with subdegree angular resolution. The Strip focal plane is populated by forty-nine Q-band and six W-band corrugated horns, each feeding a cryogenically cooled polarimeter for the detection of the Stokes $Q$ and $U$ components of the polarized signal from the sky. The Q-band channel is designed to accurately monitor Galactic polarized synchrotron emission, while the combination of Q- and W-band will allow the study of atmospheric effects at the observation site, the Observatorio del Teide, in Tenerife. In this paper we focus on the development of the Strip corrugated feed horns, including design requirements, engineering and manufacturing, as well as detailed characterization and performance verification.

Time-resolved spectra of six short gamma-ray bursts (sGRBs), measured by the {\em Swift} telescope, are used to estimate the parameters of a plerion-like model of the X-ray afterglow. The unshrouded, optically thin component of the afterglow is modelled as emanating from an expanding bubble of relativistic, shock-accelerated electrons fuelled by a central object. The electrons are injected with a power-law distribution and cool mainly by synchrotron losses. We compute posteriors for model parameters describing the central engine (e.g. spin frequency at birth, magnetic field strength) and shock acceleration (e.g. power-law index, minimum injection energy). It is found that the central engine is compatible with a millisecond magnetar, and the shock physics is compatible with what occurs in Galactic supernova remnants, assuming standard magnetic field models for the magnetar wind. Separately, we allow the magnetic field to vary arbitrarily and infer that it is roughly constant and lower in magnitude than the wind-borne extension of the inferred magnetar field. This may be due to the expansion history of the bubble, or the magnetization of the circumstellar environment of the sGRB progenitor.

Outflows are common in many astrophysical systems. In the Two Component Advective Flow ({\fontfamily{qcr}\selectfont TCAF}) paradigm which is essentially a generalized Bondi flow including rotation, viscosity and cooling effects, the outflow is originated from the hot, puffed up, post-shock region at the inner edge of the accretion disk. We consider this region to be the base of the jet carrying away matter with high velocity. In this paper, we study the spectral properties of black holes using {\fontfamily{qcr}\selectfont TCAF} which includes also a jet ({\fontfamily{qcr}\selectfont JeTCAF}) in the vertical direction of the disk plane. Soft photons from the Keplerian disk are up-scattered by the post-shock region as well as by the base of the jet and are emitted as hard radiation. We also include the bulk motion Comptonization effect by the diverging flow of jet. Our self-consistent accretion-ejection solution shows how the spectrum from the base of the jet varies with accretion rates, geometry of the flow and the collimation factor of the jet. We apply the solution to a jetted candidate GS\,1354-64 to estimate its mass outflow rate and the geometric configuration of the flow during 2015 outburst using {\it NuSTAR} observation. The estimated mass outflow to mass inflow rate is $0.12^{+0.02}_{-0.03}$. From the model fitted accretion rates, shock compression ratio and the energy spectral index, we identify the presence of hard and intermediate spectral states of the outburst. Our model fitted jet collimation factor ($f_{\rm col}$) is found to be $0.47^{+0.09}_{-0.09}$.

We report the results of our pilot millimeter observations of periodic maser sources. Using SMA, we carried out 1.3\,mm observations of G22.357$+$0.066 and G25.411$+$0.105, while ALMA 1.3 mm archival data was used in the case of G9.62+0.19E. Two continuum cores (MM1 and MM2) were detected in G22.357$+$0.066, while 3 cores (MM1 -- MM3) detected in G25.411$+$0.105. Assuming dust-to-gas ratio of 100, we derived the masses of the detected cores. Using the $^{12}$CO (2-1) and $^{13}$CO (2-1) line emission, we observed gas kinematics tracing the presence of bipolar outflows in all three star-forming regions. In the cases of G22.357$+$0.066 and G9.62+0.19E, both with similar periodic maser light curve profiles, the outflowing gas is seen in the north-west south-east direction. This suggest edge-on view of the 2 sources. G25.411$+$0.105, with a contrasting light curve profile, show a spatially co-located blue and red outflow lobes, suggesting a face-on view. Our results suggest that orientation effects may play a role in determining the characteristics of the light curves of periodic methanol masers.

Context: While stars have often been used as laboratories to study dark matter (DM), red giant branch (RGB) stars and all the rich phenomenology they encompass have frequently been overlooked by such endeavors. Aims: We study the capture, evaporation, and annihilation of weakly interacting massive particle (WIMP) DM in low-mass RGB stars ($M=0.8-1.4~\mathrm{M_{\odot}}$). Methods: We used a modified stellar evolution code to study the effects of DM self-annihilation on the structure and evolution of low-mass RGB stars. Results: We find that the number of DM particles that accumulate inside low-mass RGB stars is not only constant during this phase of evolution, but also mostly independent of the stellar mass and to some extent stellar metallicity. Moreover, we find that the energy injected into the stellar core due to DM annihilation can promote the conditions necessary for helium burning and thus trigger an early end of the RGB phase. The premature end of the RGB, which is most pronounced for DM particles with $m_\chi \simeq 100~\mathrm{GeV}$, is thus achieved at a lower helium core mass, which results in a lower luminosity at the tip of the red giant branch (TRGB). Although in the current WIMP paradigm, these effects are only relevant if the number of DM particles inside the star is extremely large, we find that for light WIMPs ($m_\chi \simeq 4~\mathrm{GeV}$), relevant deviations from the standard TRGB luminosity ($\sim 8\%$) can be achieved with conditions that can be realistic in the inner parsec of the Milky Way.

To date, only 18 exoplanets with radial velocity (RV) semi-amplitudes $<2$ m/s have had their masses directly constrained. The biggest obstacle to RV detection of such exoplanets is variability intrinsic to stars themselves, e.g. nuisance signals arising from surface magnetic activity such as rotating spots and plages, which can drown out or even mimic planetary RV signals. We use Kepler-37 - known to host three transiting planets, one of which, Kepler-37d, should be on the cusp of RV detectability with modern spectrographs - as a case study in disentangling planetary and stellar activity signals. We show how two different statistical techniques - one seeking to identify activity signals in stellar spectra, and another to model activity signals in extracted RVs and activity indicators - can enable detection of the hitherto elusive Kepler-37d. Moreover, we show that these two approaches can be complementary, and in combination, facilitate a definitive detection and precise characterisation of Kepler-37d. Its RV semi-amplitude of $1.22\pm0.31$ m/s (mass $5.4\pm1.4$ $M_\oplus$) is formally consistent with TOI-178b's $1.05^{+0.25}_{-0.30}$ m/s, the latter being the smallest detected RV signal of any transiting planet to date, though dynamical simulations suggest Kepler-37d's mass may be on the lower end of our $1\sigma$ credible interval. Its consequent density is consistent with either a water-world or that of a gaseous envelope ($\sim0.4\%$ by mass) surrounding a rocky core. Based on RV modelling and a re-analysis of Kepler-37 TTVs, we also argue that the putative (non-transiting) planet Kepler-37e should probably be stripped of its 'confirmed' status.

The Cosmic Multiperspective Event Tracker (CoMET) R&D project aims to optimize the techniques for the detection of soft-spectrum sources through very-high-energy gamma-ray observations using particle detectors (called ALTO detectors), and atmospheric Cherenkov light collectors (called CLiC detectors). The accurate reconstruction of the energies and maximum depths of gamma-ray events using a surface array only, is an especially challenging problem at low energies, and the focus of the project. In this contribution, we leverage Convolutional Neural Networks (CNNs) using the ALTO detectors only, to try to improve reconstruction performance at lower energies ( < 1 TeV ) as compared to the SEMLA analysis procedure, which is a more traditional method using manually derived features.

The mismatch between the locally measured expansion rate of the universe and the one inferred from observations of the cosmic microwave background (CMB) assuming the canonical $\Lambda$CDM model has become the new cornerstone of modern cosmology, and many new-physics set ups are rising to the challenge. Concomitant with the so-called $H_0$ problem, there is evidence of a growing tension between the CMB-preferred value and the local determination of the weighted amplitude of matter fluctuations $S_8$. It would be appealing and compelling if both the $H_0$ and $S_8$ tensions were resolved at once, but as yet none of the proposed new-physics models have done so to a satisfactory degree. Herein, we adopt a systematic approach to investigate the possible interconnection among the free parameters in several classes of models that typify the main theoretical frameworks tackling the tensions on the universe expansion rate and the clustering of matter. Our calculations are carried out using the publicly available Boltzmann solver CAMB in combination with the sampler CosmoMC. We show that even after combining the leading classes of models sampling modifications of both the early and late time universe a simultaneous solution to the $H_0$ and $S_8$ tensions remains elusive.

Baikal-GVD is a gigaton-scale neutrino observatory under construction in Lake Baikal. It currently produces about 100 GB of data every day. For their automatic processing, the Baikal Analysis and Reconstruction software (BARS) was developed. At the moment, it includes such stages as hit extraction from PMT waveforms, assembling events from raw data, assigning timestamps to events, determining the position of the optical modules using an acoustic positioning system, data quality monitoring, muon track and cascade reconstruction, as well as the alert signal generation. These stages are implemented as C++ programs which are executed sequentially one after another and can be represented as a directed acyclic graph. The most resource-consuming programs run in parallel to speed up processing. A separate Python package based on the luigi package is responsible for program execution control. Additional information such as the program execution status and run metadata are saved into a central database and then displayed on the dashboard. Results can be obtained several hours after the run completion.

We present the results of the kinematic investigations carried out with the use of spatial velocities of red giants and sub-giants containing in the $Gaia$~EDR3 catalogue. The twelve kinematic parameters of the Ogorodnikov--Milne model have been derived for stellar systems with radii 0.5 and 1.0 kpc, located along the direction the Galactic center -- the Sun -- the Galactic anticenter within the range of Galactocentric distances $R$ 0--8--16~kpc. By combining some of the local parameters the information related to the Galaxy as a whole has been received in the distance range 4--12~kpc, in particular the Galactic rotational curve, its slope, velocity gradients. We show that when using this approach, there is an alternative possibility to infer the behaviour of the Galactic rotational curve and its slope without using the Galactocentric distance $R_\odot$. The kinematic parameters derived within the Solar vicinity of 1~kpc radius are in good agreement with those given in literature.

In order to constrain the models describing circumstellar environments, it is necessary to solve the radiative transfer equation in the presence of absorption and scattering, coupled with the equation for radiative equilibrium. However, solving this problem requires much CPU time, which makes the use of automatic minimisation procedures to characterise these environments challenging. The use of approximate methods is then of primary interest. One promising candidate method is the flux-limited diffusion (FLD), which recasts the radiative transfer problem into a non-linear diffusion equation. One important aspect for the accuracy of the method lies in the implementation of appropriate boundary conditions (BCs). We present new BCs for the FLD approximation in circumstellar environments that we apply here to spherically symmetric envelopes. At the inner boundary, the entering flux (coming from the star and from the envelope itself) may be written in the FLD formalism and provides us with an adequate BC. At the free outer boundary, we used the FLD formalism to constrain the ratio of the mean radiation intensity over the emerging flux. In both cases we derived non-linear mixed BCs relating the surface values of the mean specific intensity and its gradient. We implemented these conditions and compared the results with previous benchmarks and the results of a Monte Carlo radiative transfer code. A comparison with results derived from BCs that were previously proposed in other contexts is presented as well. For all the tested cases, the average relative difference with the benchmark results is below 2\% for the temperature profile and below 6\% for the corresponding spectral energy distribution or the emerging flux. We point out that the FLD method together with the new outer BC also allows us to derive an approximation for the emerging flux.

The observations of compact star inspirals from LIGO/Virgo provide a valuable tool to study the highly uncertain equation of state (EOS) of dense matter at the densities in which the compact stars reside. It is not clear whether the merging stars are neutron stars or quark stars containing self-bound quark matter. In this work, we explore the allowed bag-model-like EOSs for the possibility of the latter by assuming the merging stars are strange quark stars (SQSs) from a Bayesian analysis employing the tidal deformability observational data of the GW170817 and GW190425 binary mergers. We consider two extreme states of strange quark matter, either in nonsuperfluid or color-flavor locked (CFL) and find the results in these two cases essentially reconcile. In particular, our results indicate that the sound speed in the SQS matter is approximately a constant close to the conformal limit of $c/\sqrt{3}$. The universal relations between the mass, the tidal defromability and the compactness are provided for the SQSs. The most probable values of the maximum mass are found to be $M_{\rm TOV}=2.10_{-0.12}^{+0.12}~(2.15_{-0.14}^{+0.16})\,M_{\odot}$ for normal (CFL) SQSs at $90\%$ confidence level. The corresponding radius and tidal deformability for a $1.4 \,M_{\odot}$ star are $R_{\rm 1.4}= 11.50_{-0.55}^{+0.52}~({11.42}_{-0.44}^{+0.52})~\rm km$ and $\Lambda_{1.4}= {650}_{-190}^{+230}~({630}_{-150}^{+220})$, respectively. We also investigate the possibility of GW190814's secondary component $m_2$ of mass $2.59_{-0.09}^{+0.08}\,M_{\odot}$ as a SQS, and find that it could be a CFL SQS with the pairing gap $\Delta$ larger than $244~\rm MeV$ and the effective bag parameter $B_{\rm eff}^{1/4}$ in the range of $170$ to $192$ MeV, at $90\%$ confidence level.

In the present work, we study the possibility of assessing the inhomogeneities of the ionic composition along the axis of the magnetic cloud using the method that was recently used by Song et al. (2021). Possible violations of the used assumptions do not allow one to draw reliable conclusions.

We present a GPU-accelerated cosmological simulation code, PhotoNs-GPU, based on algorithm of Particle Mesh Fast Multipole Method (PM-FMM), and focus on the GPU utilization and optimization. A proper interpolated method for truncated gravity is introduced to speed up the special functions in kernels. We verify the GPU code in mixed precision and different levels of interpolated method on GPU. A run with single precision is roughly two times faster that double precision for current practical cosmological simulations. But it could induce a unbiased small noise in power spectrum. Comparing with the CPU version of PhotoNs and Gadget-2, the efficiency of new code is significantly improved. Activated all the optimizations on the memory access, kernel functions and concurrency management, the peak performance of our test runs achieves 48% of the theoretical speed and the average performance approaches to 35% on GPU.

We estimate effective temperature ($T_{\rm eff}$), stellar radius, and luminosity for a sample of 271 M-dwarf stars (M0V-M7V) observed as a part of CARMENES (Calar Alto high-Resolution search for M dwarfs with Exo-earths with Near-infrared and optical Echelle Spectrographs) radial-velocity planet survey. For the first time, using the simultaneously observed high resolution (R$\sim90000$) spectra in the optical (0.52 - 0.96 $\mu$m) and near-infrared (0.96 - 1.71 $\mu$m) bands, we derive empirical calibration relationships to estimate the fundamental parameters of these low-mass stars. We select a sample of nearby and bright M-dwarfs as our calibrators for which the physical parameters are acquired from high-precision interferometric measurements. To identify the most suitable indicators of $T_{\rm eff}$, radius, and luminosity (log $L/L_{\odot}$), we inspect a range of spectral features and assess them for reliable correlations. We perform multivariate linear regression and find that the combination of pseudo equivalent widths and equivalent width ratios of the Ca II at 0.854 $\mu$m and Ca II at 0.866 $\mu$m lines in the optical and the Mg I line at 1.57 $\mu$m in the NIR give the best fitting linear functional relations for the stellar parameters with root mean square errors (RMSE) of 99K, 0.06 $R_{\odot}$ and 0.22 dex respectively. We also explore and compare our results with literature values obtained using other different methods for the same sample of M dwarfs.

To fully constrain the orbits of low mass circumstellar companions, we conduct combined analyses of the radial velocity data as well as the Gaia and Hipparcos astrometric data for eight nearby systems. Our study shows that companion-induced position and proper motion differences between Gaia and Hipparcos are significant enough to constrain orbits of low mass companions to a precision comparable with previous combined analyses of direct imaging and radial velocity data. We find that our method is robust to whether we use Gaia DR2 or Gaia EDR3, as well as whether we use all of the data, or just proper motion differences. In particular, we fully characterize the orbits of HD 190360 b and HD 16160 C for the first time. With a mass of 1.8$\pm$0.2$m_{\rm Jup}$ and an effective temperature of 123-176 K and orbiting around a Sun-like star, HD 190360 b is the smallest Jupiter-like planet with well-constrained mass and orbit, belonging to a small sample of fully characterized Jupiter analogs. It is separated from its primary star by 0.25$''$ and thus may be suitable for direct imaging by the CGI instrument of the Roman Space Telescope.

Next-generation cosmological surveys will observe larger cosmic volumes than ever before, enabling us to access information on the primordial Universe, as well as on relativistic effects. We consider forthcoming 21cm intensity mapping surveys (SKAO) and optical galaxy surveys (DESI and Euclid), combining the information via multi-tracer cross-correlations that suppress cosmic variance on ultra-large scales. In order to fully incorporate wide-angle effects and redshift-bin cross-correlations, together with lensing magnification and other relativistic effects, we use the angular power spectra, $C_\ell(z_i,z_j)$. Applying a Fisher analysis, we forecast the expected precision on $f_{\rm NL}$ and the detectability of lensing and other relativistic effects. We find that the full combination of two pairs of 21cm and galaxy surveys, one pair at low redshift and one at high redshift, could deliver $\sigma(f_{\rm NL})\sim 1.5$, detect the Doppler effect with a signal-to-noise ratio $\sim$8 and measure the lensing convergence contribution at $\sim$2% precision. In a companion paper, we show that the best-fit values of $f_{\rm NL}$ and of standard cosmological parameters are significantly biased if the lensing contribution is neglected.

We present the first measurements of Charon's far-ultraviolet surface reflectance, obtained by the Alice spectrograph on New Horizons. We find no measurable flux shortward of 1650 A, and Charon's geometric albedo is $<0.019$ ($3\sigma$) at 1600 A. From 1650--1725 A Charon's geometric albedo increases to $0.166\pm0.068$, and remains nearly constant until 1850 A. As this spectral shape is characteristic of H$_2$O ice absorption, Charon is the first Kuiper belt object with a H$_2$O ice surface to be detected in the far-ultraviolet. Charon's geometric albedo is $\sim3.7$ times lower than Enceladus' at these wavelengths, but has a very similar spectral shape. We attribute this to similarities in their surface compositions, and the difference in absolute reflectivity to a high concentration or more-absorbing contaminants on Charon's surface. Finally, we find that Charon has different solar phase behavior in the FUV than Enceladus, Mimas, Tethys, and Dione, with a stronger opposition surge than Enceladus and a shallower decline at intermediate solar phase angles than any of these Saturnian satellites.

The upcoming X-ray missions based on the microcalorimeter technology require exquisite precision in spectral simulation codes in order to match the unprecedented spectral resolution. In this work, we improve the fluorescence K$\alpha$ energies for Si II-XI and S II-XIII in the code Cloudy. In particular, we provide here a patch to update the Cloudy fluorescence energy table, originally based on Kaastra & Mewe (1993), with the laboratory energies measured by Hell et al. (2016). The new Cloudy simulations were used to model the Chandra/HETG spectra of the High Mass X-ray Binary Vela X-1 previously presented in Amato et al. (2021), showing a remarkable agreement and a dramatic improvement with respect to the current release version of Cloudy (C17.02).

Dynamical Dark Energy (DDE) is a possible solution to the Hubble tension. This work analyses the combination of Baryon Acoustic Oscillation (BAO) data that include $19$ points from the range $0.11 \le z \le 2.40$ and additional points from the Cosmic Microwave Background (CMB) Distant Priors. We test equation of states with a Linear, CPL and a Logarithm dependence on the redshift. DDE seems to be stronger then the standard $\Lambda$CDM model statistically using different selection criteria. The result calls for new observations and stimulates the investigation of alternative theoretical models beyond the standard model. We also test the same dataset on the $\Omega_k$CDM model but for it, $\Lambda$CDM gives better statistical measures.

The runaway greenhouse represents the ultimate climate catastrophe for rocky, Earth-like worlds: when the incoming stellar flux cannot be balanced by radiation to space, the oceans evaporate and exacerbate heating, turning the planet into a hot wasteland with a steam atmosphere overlying a possibly molten magma surface. The equilibrium state beyond the runaway greenhouse instellation limit depends on the radiative properties of the atmosphere and its temperature structure. Here, we use 1-D radiative-convective models of steam atmospheres to explore the transition from the tropospheric radiation limit to the post-runaway climate state. To facilitate eventual simulations with 3-D global circulation models, a computationally efficient band-grey model is developed, which is capable of reproducing the key features of the more comprehensive calculations. We analyze two factors which determine the equilibrated surface temperature of post-runaway planets. The infrared cooling of the planet is strongly enhanced by the penetration of the dry adiabat into the optically thin upper regions of the atmosphere. In addition, thermal emission of both shortwave and near-IR fluxes from the hot lower atmospheric layers, which can radiate through window regions of the spectrum, is quantified. Astronomical surveys of rocky exoplanets in the runaway greenhouse state may discriminate these features using multi-wavelength observations.

The Einstein equivalence principle in the electromagnetic sector can be violated in modifications of gravity theory generated by a multiplicative coupling of a scalar field to the electromagnetic Lagrangian. In such theories, deviations of the standard result for the cosmic distance duality relation, and a variation of the fine structure constant are expected and are unequivocally intertwined. In this paper, we search for these possible cosmological signatures by using galaxy cluster gas mass fraction measurements and cosmic chronometers. No significant departure from general relativity is found regardless of our assumptions about cosmic curvature or a possible depletion factor evolution in cluster measurements.

The first stage of the construction of the deep underwater neutrino telescope Baikal-GVD is planned to be completed in 2024. The second stage of the detector deployment is planned to be carried out using a data acquisition system based on fibre optic technologies, which will allow for increased data throughput and more flexible trigger conditions. A dedicated test facility has been built and deployed at the Baikal-GVD site to test the new technological solutions. We present the principles of operation and results of tests of the new data acquisition system.

In this letter we introduce the non-linear partial differential equation (PDE) $\partial^2_{\tau} \pi \propto (\vec\nabla \pi)^2$ showing a new type of instability. Such equations appear in the effective field theory (EFT) of dark energy for the $k$-essence model as well as in many other theories based on the EFT formalism. We demonstrate the occurrence of instability in the cosmological context using a relativistic $N$-body code, and we study it mathematically in 3+1 dimensions within spherical symmetry. We show that this term dominates for the low speed of sound limit where some important linear terms are suppressed.

We describe the \texttt{Period detection and Identification Pipeline Suite} (PIPS) -- a new, fast, and user-friendly platform for period detection and analysis of astrophysical time-series data. PIPS is an open-source Python package that provides various pre-implemented methods and a customisable framework for automated, robust period measurements with principled uncertainties and statistical significance calculations. In addition to detailing the general algorithm that underlies PIPS, this paper discusses one of PIPS' central and novel features, the Fourier-likelihood periodogram, and compares its performance to existing methods. The resulting improved performance implies that one can construct deeper, larger, and more reliable sets of derived properties from various observations, including all-sky surveys. We present a comprehensive validation of PIPS against artificially generated data, which demonstrates the reliable performance of our algorithm for a class of periodic variable stars (RR Lyrae stars). We further showcase an application to recently obtained data of variable stars in the globular cluster M15.

We present a Markov chain Monte Carlo pipeline that can be used for robust and unbiased constraints of $f(R)$ gravity using galaxy cluster number counts. This pipeline makes use of a detailed modelling of the halo mass function in $f(R)$ gravity, which is based on the spherical collapse model and calibrated by simulations, and fully accounts for the effects of the fifth force on the dynamical mass, the halo concentration and the observable-mass scaling relations. Using a set of mock cluster catalogues observed through the thermal Sunyaev-Zel'dovich effect, we demonstrate that this pipeline, which constrains the present-day background scalar field $f_{R0}$, performs very well for both $\Lambda$CDM and $f(R)$ fiducial cosmologies. We find that using an incomplete treatment of the scaling relation, which could deviate from the usual power-law behaviour in $f(R)$ gravity, can lead to imprecise and biased constraints. We also find that various degeneracies between the modified gravity, cosmological and scaling relation parameters can significantly affect the constraints, and show how this can be rectified by using tighter priors and better knowledge of the cosmological and scaling relation parameters. Our pipeline can be easily extended to other modified gravity models, to test gravity on large scales using galaxy cluster catalogues from ongoing and upcoming surveys.

The apparent angular size of the shadow of a black hole in an expanding Universe is redshift-dependent. Since cosmological redshifts change with time - known as the redshift drift - all redshift-dependent quantities acquire a time-dependence, and a fortiori so do black hole shadows. We find a mathematical description of the black hole shadow drift and show that the amplitude of this effect is of order $10^{-16}$ per day for M87$^{\star}$. While this effect is small, we argue that its non-detection can be used as a novel probe of the equivalence principle and derive a constraint of $|\dot{G}/G| \leq 10^{-3}-10^{-4}$ per year. The effect of redshift drift on the visibility amplitude and frequency of the universal interferometric signatures of photon rings is also discussed, which we show to be very similar to the shadow drift. This is of particular interest for future experiments involving spectroscopic and interferometric techniques, which could make observations of photon rings and their frequency drifts viable.

We study cosmological dynamics of the energy-momentum squared gravity. By adding the squared of the matter field's energy-momentum tensor ($\zeta\, \textbf{T}^{2}$) to the Einstein Hilbert action, we obtain the Einstein's field equations and study the conservation law. We show that the presence of $\zeta\, \textbf{T}^{2}$ term, breaks the conservation of the energy-momentum tensor of the matter fields. However, an effective energy-momentum tensor in this model is conserved in time. By considering the FRW metric as the background, we find the Friedmann equations and by which we explore the cosmological inflation in $\zeta\,\textbf{T}^{2}$ model. We perform numerical analysis on the perturbation parameters and compare the results with Planck2018 different data sets at $68\%$ and $95\%$ CL, to obtain some constraints on the coupling parameter $\zeta$. We show that \textbf{ for $0< \zeta \leq 2.1\times 10^{-5}$, the $\zeta\, \textbf{T}^{2}$ gravity is an observationally viable model of inflation.

We study axion dark matter production from a misalignment mechanism in scenarios featuring a general nonstandard cosmology. Before the onset of Big Bang nucleosynthesis, the energy density of the universe is dominated by a particle field $\phi$ described by a general equation of state $\omega$. The ensuing enhancement of the Hubble expansion rate decreases the temperature at which axions start to oscillate, opening this way the possibility for axions heavier than in the standard window. This is the case for kination, or in general for scenarios with $\omega > 1/3$. However, if $\omega < 1/3$, as in the case of an early matter domination, the decay of $\phi$ injects additional entropy relative to the case of the standard model, diluting this way the preexisting axion abundance, and rendering lighter axions viable. For a misalignment angle $0.5 < \theta_i < \pi/\sqrt{3}$, the usual axion window becomes expanded to $4 \times 10^{-9}$ eV $\lesssim m_a \lesssim 2 \times 10^{-5}$ eV for the case of an early matter domination, or to $2 \times 10^{-6}$ eV $\lesssim m_a \lesssim 10^{-2}$ eV for the case of kination. Interestingly, the coupling axion-photon in such a wider range can be probed with next generation experiments such as ABRACADABRA, KLASH, ADMX, MADMAX, and ORGAN. Axion dark matter searches may therefore provide a unique tool to probe the history of the universe before Big Bang nucleosynthesis.

Silica, water and hydrogen are known to be the major components of celestial bodies, and have significant influence on the formation and evolution of giant planets, such as Uranus and Neptune. Thus, it is of fundamental importance to investigate their states and possible reactions under the planetary conditions. Here, using advanced crystal structure searches and first-principles calculations in the Si-O-H system, we find that a silica-water compound (SiO2)2(H2O) and a silica-hydrogen compound SiO2H2 can exist under high pressures above 450 and 650 GPa, respectively. Further simulations reveal that, at high pressure and high temperature conditions corresponding to the interiors of Uranus and Neptune, these compounds exhibit superionic behavior, in which protons diffuse freely like liquid while the silicon and oxygen framework is fixed as solid. Therefore, these superionic silica-water and silica-hydrogen compounds could be regarded as important components of the deep mantle or core of giants, which also provides an alternative origin for their anomalous magnetic fields. These unexpected physical and chemical properties of the most common natural materials at high pressure offer key clues to understand some abstruse issues including demixing and erosion of the core in giant planets, and shed light on building reliable models for solar giants and exoplanets.

The understanding of the physical processes that lead to the origin of matter in the early Universe, creating both an excess of matter over anti-matter that survived until the present and a dark matter component, is one of the most fascinating challenges in modern science. The problem cannot be addressed within our current description of fundamental physics and, therefore, it currently provides a very strong evidence of new physics. Solutions can either reside in a modification of the standard model of elementary particle physics or in a modification of the way we describe gravity, based on general relativity, or at the interface of both. We will mainly discuss the first class of solutions. Traditionally, models that separately explain either the matter-antimatter asymmetry of the Universe or dark matter have been proposed. However, in the last years there has also been an accreted interest and intense activity on scenarios able to provide a unified picture of the origin of matter in the early universe. In this review we discuss some of the main ideas emphasising primarily those models that have more chances to be experimentally tested during next years. Moreover, after a general discussion, we will focus on extensions of the standard model that can also address neutrino masses and mixing, since this is currently the only evidence of physics beyond the standard model coming directly from particle physics and it is, therefore, reasonable they might also provide a solution to the problem of the origin of matter in the universe.

After the detection of gravitational waves from compact binary coalescences, the search for transient gravitational-wave signals with less well-defined waveforms for which matched filtering is not well-suited is one of the frontiers for gravitational-wave astronomy. Broadly classified into "short" $ \lesssim 1~$\,s and "long" $ \gtrsim 1~$\,s duration signals, these signals are expected from a variety of astrophysical processes, including non-axisymmetric deformations in magnetars or eccentric binary black hole coalescences. In this work, we present a search for long-duration gravitational-wave transients from Advanced LIGO and Advanced Virgo's third observing run from April 2019 to March 2020. For this search, we use minimal assumptions for the sky location, event time, waveform morphology, and duration of the source. The search covers the range of $2~\text{--}~ 500$~s in duration and a frequency band of $24 - 2048$ Hz. We find no significant triggers within this parameter space; we report sensitivity limits on the signal strength of gravitational waves characterized by the root-sum-square amplitude $h_{\mathrm{rss}}$ as a function of waveform morphology. These $h_{\mathrm{rss}}$ limits improve upon the results from the second observing run by an average factor of 1.8.

We study Quantum Field Theory (QFT) on a background de Sitter spacetime dS$_{d+1}$. Our main tool is the Hilbert space decomposition in irreducible unitarity representations of its isometry group $SO(d+1,1)$. Throughout this work, we focus on the late-time physics of dS$_{d+1}$, in particular on the boundary operators that appear in the late-time expansion of bulk local operators. As a first application of the Hilbert space formalism, we recover the K\"allen-Lehmann spectral decomposition of bulk two-point functions. In the process, we exhibit a relation between poles in the corresponding spectral densities and boundary CFT data. Next, we study the conformal partial wave decomposition of four-point functions of boundary operators. These correlation functions are very similar to the ones of standard conformal field theory, but have different positivity properties that follow from unitarity in de Sitter. We conclude by proposing a non-perturbative conformal bootstrap approach to the study of these late-time four-point functions, and we illustrate our proposal with a concrete example for QFT in dS$_2$.

We study possible spatial domains containing expanding extra dimensions. We show that they are predicted in the framework of $f(R)$ gravity and could appear due to quantum fluctuations during inflation. Their interior is characterized by the multidimensional curvature ultimately tending to zero and a slowly growing size of the extra dimensions.

We include the single graviton loop contribution to the linearized Einstein equation. Explicit results are obtained for one loop corrections to the propagation of gravitational radiation. Although suppressed by a minuscule loop-counting parameter, these corrections are enhanced by the square of the number of inflationary e-foldings. One consequence is that perturbation theory breaks down for a very long epoch of primordial inflation. Another consequence is that the one loop correction to the tensor power spectrum might be observable, in the far future, after the full development of 21cm cosmology.

Atom interferometers (AIs) as gravitational-wave (GW) detector had been proposed a decade ago. Both ground and space-based projects will be in construction and preparation in a near future. In this paper, for the first time, we investigate the potential of the space-borne AIs on detecting GW standard sirens and hence the applications on cosmology. We consider AEDGE as our fiducial AI GW detector and estimate the number of bright sirens that would be obtained within a 5-years data-taking period. We then construct the mock catalogue of bright sirens and predict their ability on constraining such as the Hubble constant, dynamics of dark energy, and modified gravity theory. The preliminary results show that there should be of order $\mathcal{O} (30)$ bright sirens detected within 5 years observation time by AEDGE. The bright sirens alone can measure $H_0$ with precision 2.1\%, which is sufficient to arbitrate the Hubble tension. Combining current most precise electromagnetic experiments, the inclusion of AEDGE bright sirens can improve the measurement of equation of state of dark energy, though marginally. However, by modifying GW propagation on cosmological scales, the deviations from general relativity (modified gravity theory effects ) can be constrained at 5.7\% precision level, which is two times better than by 10-years operation of LIGO, Virgo and KAGRA network.

Astrophysical neutrino fluxes are often modeled as power laws of the energy. This is reasonable in the case of hadronic sources, but it does not capture the behavior in photohadronic sources, where the spectrum depends on the properties of the target photons on which protons collide. This limits the possibility of a unified treatment of different sources. In order to overcome this difficulty, we model the target photons by a blackbody spectrum. This model is sufficiently flexible to reproduce neutrino fluxes from known photohadronic sources; we apply it to study the sensitivity of Dense Neutrino Arrays, Neutrino Telescopes and Neutrino Radio Arrays to photohadronic sources. We also classify the flavor composition of the neutrino spectrum in terms of the parameter space. We discuss the interplay with the experiments, studying the changes in the track-to-shower ratio induced by different flavor compositions, both within and outside the region of the Glashow resonance.

In this paper, we study a non-canonical extension of a supergravity-motivated model acting as a vivid counterexample to the cosmic no-hair conjecture due to its unusual coupling between scalar and electromagnetic fields. In particular, a canonical scalar field is replaced by the string-inspired Dirac-Born-Infeld one in this extension. As a result, exact anisotropic inflationary solutions for this Dirac-Born-Infeld model are figured out under a constant-roll condition. Furthermore, numerical calculations are performed to verify that these anisotropic constant-roll solutions are indeed attractive during their inflationary phase.

We present data on the Baikal water luminescence collected with the Baikal-GVD neutrino telescope. This three-dimensional array of photo-sensors allows the observation of time and spatial variations of the ambient light field. We report on annual increase of luminescence activity in years 2018-2020. We observed a unique event of a highly luminescent layer propagating upwards with a maximum speed of 28 m/day for the first time.